The vertical structure of the F ring of Saturn from ring-plane crossings

We present a photometric model of the rings of Saturn which includes the main rings and an F ring, inclined to the main rings, with a Gaussian vertical profile of optical depth. This model reproduces the asymmetry in brightness between the east and west ansae of the rings of Saturn that was observed by the Hubble Space Telescope (HST) within a few hours after the Earth ring-plane crossing (RPX) of 10 August 1995. The model shows that during this observation the inclined F ring unevenly blocked the east and west ansae of the main rings. The brightness asymmetry produced by the model is highly sensitive to the vertical thickness and radial optical depth of the F ring. The F-ring model that best matches the observations has a vertical full width at half maximum of 13 ± 7 km and an equivalent depth of 10 ± 4 km. The model also reproduces the shape of the HST profiles of ring brightness vs. distance from Saturn, both before and after the time of ring-plane crossing. Smaller asymmetries observed before the RPX, when the Earth was on the dark side of the rings, cannot be explained by blocking of the main rings by the F ring or vice versa and are probably instead due to the intrinsic longitudinal variation exhibited by the F ring.     

Far ultraviolet spectral properties of Saturn’s rings from Cassini UVIS

Spectra taken by the Cassini Ultraviolet Imaging Spectrograph (UVIS) of Saturn’s C ring, B ring, Cassini Division, and A ring have been analyzed in order to characterize ring particle surface properties and water ice abundance in the rings. UVIS spectra sense the outer few microns of the ring particles. Spectra of the normalized reflectance (I/F) in all four regions show a characteristic water ice absorption feature near 165 nm. Our analysis shows that the fractional abundance of surface water ice is largest in the outer B ring and decreases by over a factor of 2 across the inner C ring. We calculate the mean path length of UV photons through icy ring particle regolith and the scattering asymmetry parameter using a Hapke reflectance model and a Shkuratov reflectance model to match the location of the water ice absorption edge in the data. Both models give similar retrieved values of the photon mean length, however the retrieved asymmetry (g) values are different. The photon mean path lengths are nearly uniform across the B and A rings. Shortward of 165 nm the rings exhibit a slope that turns up towards shorter wavelengths, while the UV slope of 180/150 nm (reflectance outside the water absorption ratioed to that inside the absorption band) tracks I/F with maxima in the outer B ring and in the central A ring. Retrieved values of the scattering asymmetry parameter show the regolith grains to be highly backscattering in the FUV spectral regime.     

Photometric modeling of Saturn's rings. II. Azimuthal asymmetry in reflected and transmitted light

Dynamical N-body simulations (Salo, 992, Nature 359, 619) suggest the formation of trailing density enhancements in the outer portions of Saturn's rings, due to local gravitational instabilities. These Julian-Toomre type wakes, having a pitch angle of about 20°-25° with respect to the local tangential direction, seem to provide a plausible explanation for the observed quadrupole brightness variation in Saturn's A ring (Salo and Karjalainen, 1999, Bull. Am. Astron. Soc. 31, 1160; French et al., 2000, Bull. Am. Astron. Soc. 32, 806; Porco et al., 2001, Bull. Am. Astron. Soc. 33, 1091). We have carried out systematic photometric modeling of gravitational wake structures seen in dynamical simulations, performed for the parameter values of the A ring, using the Monte Carlo radiative transfer code described in Salo and Karjalainen (2003, Icarus 164, 428). Comparisons to the observed asymmetry in various cases are presented (asymmetry in reflected and transmitted light, ring longitude and opening angle dependence), in all cases confirming the applicability of the wake model. Typically, minimum brightness corresponds to viewing/illumination along the long axis of wakes; however, the sense of modeled asymmetry reverses at small tilt angles in diffuse transmission. Implications of wakes on the occultation optical depth profiles and the A ring overall brightness behavior are also discussed: it is shown that the wake structure needs to be taken into account when the Cassini occultation profiles for the A ring are interpreted in terms of variations in surface density. Also, the presence of wakes offers a plausible explanation for the inverse tilt effect seen in the mid A-ring.     

More Microwave Observations of Saturn: Modeling the Ring with a Monte Carlo Radiative Transfer Code

We present a second epoch of Very Large Array Saturn observations taken in February 1997 spanning wavelengths 1.3-21 cm. These observations complement earlier observations at Saturn's autumnal equinox in November 1995. In this epoch, however, we generally have better signal-to-noise ratios and the ring inclination of the present observations was −5.0°, whereas the previous observations were made with ring inclination +2.7°.Our observations confirm the latitudinal structure on the saturnian disk as seen at 2.0, 3.6, and 6.1 cm. We also see some latitudinal structure at 1.3 cm for the first time. The details of this structure have changed dramatically from those reported by I. de Pater and J. R. Dickel (1991, Icarus94, 474-492) for the 1980s and are consistent with those seen in F. van der Tak et al. (1999, Icarus142, 125-147). The most prominent features are a pair of brightness enhancements just inside the edges of the Equatorial Zone.The rings do not show the east-west asymmetry seen in our previous epoch, perhaps indicative of a viewing angle effect on the scattering properties of the rings. The radial trend in brightness in the ansae is generally consistent with that expected from optical depth variations and increasing distance from the source of scattered light. In particular the increased optical depth towards the center of the C ring is evident. Azimuthal variation in brightness in the C ring shows the forward scattering expected of Mie scattering. By contrast, the A and B rings show little or no azimuthal variation.We present Monte Carlo simulations of the ring brightness under the assumptions of isotropic and Mie scattering. These are the first synthetic maps of Saturn which can be directly compared to the images we obtained. Neither model fits all the data well. However, a hybrid model combining isotropic and Mie scattering does fit well. We interpret the consistency with isotropic scattering in the outer rings as an indication that near-field effects may be important. This in turn implies geometrically thin rings, as predicted by dynamical simulations of these rings.     

Examining the wake structure in Saturn's rings from microwave observations over varying ring opening angles and wavelengths

Over the last 15 to 20 years several high quality, high resolution data have been taken with the very large array (VLA). These data exhibit a wide range of ring opening angles (|B|=0 to 26°) and wavelengths (λ=0.7 to 20 cm). At these wavelengths the primary flux from the rings is scattered saturnian thermal emission, with a small contribution coming from the ring particles' own thermal emission. Much of the data do show signs of asymmetries due to wakes either on the ansae or the portion of the rings which occult the planet. As in previous work, we use our Monte Carlo radiative transfer code including idealized wakes [Dunn, D.E., Molnar, L.A., Fix, J.D., 2002. Icarus 160, 132-160; Dunn, D.E., Molnar, L.A., Niehof, J.T., de Pater, I., Lissauer, J.L., 2004. Icarus 171, 183-198] to model the relative contributions of the scattered and thermal radiation emanating from the rings and compare the results to that seen in the data. Although the models do give satisfactory fits to all of our data, we find that no single model can simulate the data at all different |B| and λ. We find that one model works best for moderate and low |B| and another one at higher |B|. The main difference between these models is the ratio of the wake width to their separation. We similarly find that the 2 cm data require higher density wakes than the longer wavelength data, perhaps caused by a preponderance of somewhat smaller ring material in the wakes. We further find evidence for an increase in the physical temperature of the rings with increasing |B|. Continuous observations are required to determine whether the above results regarding variations in wake parameters with |B| and λ are indeed caused by these parameters, or instead by changes over time.     

HST observations of azimuthal asymmetry in Saturn's rings

From 378 Hubble Space Telescope WFPC2 images obtained between 1996-2004, we have measured the detailed nature of azimuthal brightness variations in Saturn's rings. The extensive geometric coverage, high spatial resolution (), and photometric precision of the UBVRI images have enabled us to determine the dependence of the asymmetry amplitude and longitude of minimum brightness on orbital radius, ring elevation, wavelength, solar phase angle, and solar longitude. We explore a suite of dynamical models of self-gravity wakes for two particle size distributions: a single size and a power law distribution spanning a decade in particle radius. From these N-body simulations, we calculate the resultant wake-driven brightness asymmetry for any given illumination and viewing geometry. The models reproduce many of the observed properties of the asymmetry, including the shape and location of the brightness minimum and the trends with ring elevation and solar longitude. They also account for the “tilt effect” in the A and B rings: the change in mean ring brightness with effective ring opening angle, |Beff|. The predicted asymmetry depends sensitively on dynamical ring particle properties such as the coefficient of restitution and internal mass density, and relatively weakly on photometric parameters such as albedo and scattering phase function. The asymmetry is strongest in the A ring, reaching a maximum amplitude A∼25% near a=128,000 km. Here, the observations are well-matched by an internal particle density near 450 kg m⁻³ and a narrow particle size distribution. The B ring shows significant asymmetry (∼5%) in regions of relatively low optical depth (τ∼0.7). In the middle and outer B ring, where τ?1, the asymmetry is much weaker (∼1%), and in the C ring, A<0.5%. The asymmetry diminishes near opposition and at shorter wavelengths, where the albedo of the ring particles is lower and multiple-scattering effects are diminished. The asymmetry amplitude varies strongly with ring elevation angle, reaching a peak near |Beff|=10° in the A ring and at |Beff|=15-20° in the B ring. These trends provide an estimate of the thickness of the self-gravity wakes responsible for the asymmetry. Local radial variations in the amplitude of the asymmetry within both the A and B rings are probably caused by regional differences in the particle size distribution.     

Self-gravity wakes and radial structure of Saturn's B ring

We analyze stellar occultations by Saturn's rings observed with the Cassini Ultraviolet Imaging Spectrograph and find large variations in the apparent normal optical depth of the B ring with viewing angle. The line-of-sight optical depth is roughly independent of the viewing angle out of the ring plane so that optical depth is independent of the path length of the line-of-sight. This suggests the ring is composed of virtually opaque clumps separated by nearly transparent gaps, with the relative abundance of clumps and gaps controlling the observed optical depth. The observations can be explained with a model of self-gravity wakes like those observed in the A ring. These trailing spiral density enhancements are due to the competing processes of self-gravitational accretion of ring particles and Kepler shear. The B ring wakes are flatter and more closely packed than their neighbors in the A ring, with height-to-width ratios <0.1 for most of the ring. The self-gravity wakes are seen in all regions of the B ring that are not opaque. The observed variation in total B ring optical depth is explained by the amount of relatively empty space between the self-gravity wakes. Wakes are more tightly packed in regions where the apparent normal optical depth is high, and the wakes are more widely spaced in lower optical depth regions. The normal optical depth of the gaps between the wakes is typically less than 0.5 and shows no correlation with position or overall optical depth in the ring. The wake height-to-width ratio varies with the overall optical depth, with flatter, more tightly packed wakes as the overall optical depth increases. The highly flattened profile of the wakes suggests that the self-gravity wakes in Saturn's B ring correspond to a monolayer of the largest particles in the ring. The wakes are canted to the orbital direction in the trailing sense, with a trend of decreasing cant angle with increasing orbital radius in the B ring. We present self-gravity wake properties across the B ring that can be used in radiative transfer modeling of the ring. A high radial resolution (∼10 m) scan of one part of the B ring during a grazing occultation shows a dominant wavelength of 160 m due to structures that have zero cant angle. These structures are seen at the same radial wavelength on both ingress and egress, but the individual peaks and troughs in optical depth do not match between ingress and egress. The structures are therefore not continuous ringlets and may be a manifestation of viscous overstability.     

A multilayer model for thermal infrared emission of Saturn's rings: Basic formulation and implications for Earth-based observations

We present our new model for the thermal infrared emission of Saturn's rings based on a multilayer approximation. In our model, (1) the equation of classical radiative transfer is solved directly for both visible and infrared light, (2) the vertical heterogeneity of spin frequencies of ring particles is taken into account, and (3) the heat transport due to particles motion in the vertical and azimuthal directions is taken into account. We adopt a bimodal size distribution, in which rapidly spinning small particles (whose spin periods are shorter than the thermal relaxation time) with large orbital inclinations have spherically symmetric temperatures, whereas non-spinning large particles (conventionally called slow rotators) with small orbital inclinations are heated up only on their illuminated sides. The most important physical parameters, which control ring temperatures, are the albedo in visible light, the fraction of fast rotators (ffast) in the optical depth, and the thermal inertia. In the present paper, we apply the model to Earth-based observations. Our model can well reproduce the observed temperature for all the main rings (A, B, and C rings), although we cannot determine exact values of the physical parameters due to degeneracy among them. Nevertheless, the range of the estimated albedo is limited to 0-0.52±0.05, 0.55±0.07-0.74±0.03, and 0.51±0.07-0.74±0.06 for the C, B, and A rings, respectively. These lower and upper limits are obtained assuming all ring particles to be either fast and slow rotators, respectively. For the C ring, at least some fraction of slow rotators is necessary (ffast?0.9) in order for the fitted albedo to be positive. For the A and B rings, non-zero fraction of fast rotators (ffast?0.1-0.2) is favorable, since the increase of the brightness temperature with increasing solar elevation angle is enhanced with some fraction of fast rotators.     

Keck infrared observations of Saturn's main rings bracketing Earth's August 1995 ring plane crossing

We present results of near-infrared (2.26 μm) observations of Saturn's main rings taken with the W.M. Keck telescope during August 8-11, 1995, surrounding the time that Earth crossed Saturn's ring plane. These observations provide a unique opportunity to study the evolution of the ring brightness in detail, and by combining our data with Hubble Space Telescope (HST) results (Nicholson et al., 1996, Science 272, 453-616), we extend the 12-hour HST time span to several days around the time of ring plane crossing (RPX). In this paper, we focus on the temporal evolution of the brightness in Saturn's main rings. We examine both edge-on ring profiles and radial profiles obtained by “onion-peeling” the edge-on data. Before RPX, when the dark (unlit) face of the rings was observed, the inner C ring (including the Colombo gap), the Maxwell gap, Cassini Division and F ring region were very bright in transmitted light. After RPX, the main rings brighten rapidly, as expected. The profiles show east-west asymmetries both before and after RPX. Prior to RPX, the evolution in ring brightness of the Keck and HST data match one another quite well. The west side of the rings showed a nonlinear variation in brightness during the last hours before ring plane crossing, suggestive of clumping and longitudinal asymmetries in the F ring. Immediately after RPX, the east side of the rings brightened more rapidly than the west. A quantitative comparison of the Keck and HST data reveals that the rings were redder before RPX than after; we ascribe this difference to the enhanced multiple scattering of photons passing through to the unlit side of the rings.     

Morphology and variability of the Titan ringlet and Huygens ringlet edges

We present a forward modeling approach for determining, in part, the ring particle spatial distribution in the vicinity of sharp ring or ringlet edges. Synthetic edge occultation profiles are computed based on a two-parameter particle spatial distribution model. One parameter, h, characterizes the vertical extent of the ring and the other, δ, characterizes the radial scale over which the ring optical depth transitions from the background ring value to zero. We compare our synthetic occultation profiles to high resolution stellar occultation light curves observed by the Cassini Ultraviolet Imaging Spectrograph (UVIS) High Speed Photometer (HSP) for occultations by the Titan ringlet and Huygens ringlet edges.More than 100 stellar occultations of the Huygens ringlet and Titan ringlet edges were studied, comprising 343 independent occultation cuts of the edges of these two ringlets. In 237 of these profiles the measured light-curve was fit well with our two-parameter edge model. Of the remaining edge occultations, 69 contained structure that could only be fit with extremely large values of the ring-plane vertical thickness (h > 1 km) or by adopting a different model for the radial profile of the ring optical depth. An additional 37 could not be fit by our two-parameter model.Certain occultations at low ring-plane incidence angles as well as occultations nearly tangent to the ring edge allow the direct measurement of the radial scale over which the particle packing varies at the edge of the ringlet. In 24 occultations with these particular viewing geometries, we find a wide variation in the radial scale of the edge. We are able to constrain the vertical extent of the rings at the edge to less than ∼300 m in the 70% of the occultations with appropriate viewing geometry, however tighter constraints could not be placed on h due to the weaker sensitivity of the occultation profile to vertical thickness compared to its sensitivity to δ.Many occultations of a single edge could not be fit to a single value of δ, indicating large temporal or azimuthal variability, although the azimuthal variation in δ with respect to the longitudes of various moons in the system did not show any discernible pattern.     

How much oxygen is too much? Constraining Saturn's ring atmosphere

The discovery of a molecular oxygen atmosphere around Saturn's rings has important implications for the electrodynamics of the ring system. Its existence was inferred from the Cassini in situ detection of molecular oxygen ions above the rings during Saturn Orbit Insertion in 2004. Molecular oxygen is difficult to observe remotely, and theoretical estimates have yielded only a lower limit (Nn?¹⁰¹³ cm⁻²) to the O₂ column density. Comparison with observations has previously concerned matching ion densities at spacecraft altitudes far larger than the scale height of the neutral atmosphere. This is further complicated by charged particle propagation effects in Saturn's offset magnetic field. In this study we adopt a complementary approach, by focusing on bulk atmospheric properties and using additional aspects of the Cassini observations to place an upper limit on the column density. We develop a simple analytic model of the molecular atmosphere and its photo-ionization and dissociation products, with Nn a free parameter. Heating of the neutrals by viscous stirring, cooling by collisions with the rings, and torquing by collisions with pickup ions are all included in the model. We limit the neutral scale height to h?3000 km using the INMS neutral density nondetection over the A ring. A first upper limit to the neutral column is derived by using our model to reassess O₂ production and loss rates. Two further limits are then obtained from Cassini observations: corotation of the observed ions with the planet implies that the height-integrated conductivity of the ring atmosphere is less than that of Saturn's ionosphere; and the nondetection of fluorescent atomic oxygen over the rings constrains the molecular column from which it is produced via photo-dissociation. These latter limits are independent of production and loss rates and are only weakly dependent on temperature. From the three independent methods described, we obtain similar limits: Nn?2×¹⁰¹⁵ cm⁻². The mean free path for collisions between neutrals thus cannot be very much smaller than the scale height.     

Moonlet Collisions and the Effects of Tidally Modified Accretion in Saturn's F Ring

We both test and offer an alternative to a meteoroid bombardment model (M. R. Showalter 1998, Science282, 1099-1102) and suggest that anomalous localized brightenings in the F ring observed by Voyager result from disruptive collisions involving poorly consolidated moonlets, or “rubble piles.” This model can also explain the transient events observed during ring plane crossing. We have developed an evolutionary model that considers both the competing effects of accretion and disruption at the location of the F ring. Our numerical model is a Markov process where probabilities of mass transfer between the states of the system form a “transition matrix.” Successive multiplications of this matrix by the state vector generate expectation values of the distribution after each time step as the system approaches quasi-equilibrium. Competing effects of accretion and disruption in the F ring are found to lead to a bimodal distribution of ring particle sizes. In fact, our simulation predicts the presence of a belt of kilometer-sized moonlets in the F ring. These moonlets may continually disrupt one another and re-accrete on short time scales. We also agree with J. N. Cuzzi and J. A. Burns (1988, Icarus74, 284-324), who suggest that the classical F ring itself may be the consequence of a relatively recent collision between two of the largest of these yet unseen objects. Cassini observations can confirm the existence of the moonlet belt by directly observing these objects or the waves they create in the rings.     

Structure of self-gravity wakes in Saturn's A ring as measured by Cassini CIRS

The CIRS infrared spectrometer onboard the Cassini spacecraft has scanned Saturn's A ring azimuthally from several viewing angles since its orbit insertion in 2004. A quadrupolar asymmetry has been detected in this ring at spacecraft elevations ranging between 16° to 37°. Its fractional amplitude decreases from 22% to 8% from 20° to 37° elevations. The patterns observed in two almost complete azimuthal scans at elevations 20° and 36° strongly favor the self-gravity wakes as the origin of the asymmetry. The elliptical, infinite cylinder model of Hedman et al. [Hedman, M.M., Nicholson, P.D., Salo, H., Wallis, B.D., Buratti, B.J., Baines, K.H., Brown, R.H., Clark, R.N., 2007. Astron. J. 133, 2624-2629] can reproduce the CIRS observations well. Such wakes are found to have an average height-to-spacing ratio H/λ=0.1607±0.0002, a width-over-spacing W/λ=0.3833±0.0008. Gaps between wakes, which are filled with particles, have an optical depth τG=0.1231±0.0005. The wakes mean pitch angle ΦW is 70.70°±0.07°, relative to the radial direction. The comparison of ground-based visible data with CIRS observations constrains the A ring to be a monolayer. For a surface mass density of 40 g cm⁻² [Tiscarino, M.S., Burns, J.A., Nicholson, P.D., Hedman, M.M., Porco, C.C., 2007. Icarus 189, 14-34], the expected spacing of wakes is λ≈60 m. Their height and width would then be H≈10 m and W≈24 m, values that match the maximum size of particles in this ring as determined from ground-based stellar occultations [French, R.G., Nicholson, P.D., 2000. Icarus 145, 502-523].     

Density waves in Cassini UVIS stellar occultations: 1. The Cassini Division

We analyze density waves in the Cassini Division of Saturn's rings revealed by multiple stellar occultations by Saturn's rings observed with the Cassini Ultraviolet Imaging Spectrograph. The dispersion and damping of density waves provide information on the local ring surface mass density and viscosity. Several waves in the Cassini Division are on gradients in the background optical depth, and we find that the dispersion of the wave reflects a change in the underlying surface mass density. We find that over most of the Cassini Division the ring opacity (the ratio of optical depth to surface mass density) is nearly constant and is ∼5 times higher than the opacity in the A ring where most density waves are found. However, the Cassini Division ramp, a 1100-km-wide, nearly featureless region of low optical depth that connects the Cassini Division to the inner edge of the A ring, has an opacity like that of the A ring and significantly less than that in the rest of the Cassini Division. This is consistent with particles in the ramp originating in the A ring and being transported into the Cassini Division through ballistic transport processes. Damping of the waves in the Cassini Division suggests a vertical thickness of 3–6 m. Using a mean opacity of 0.1 cm²/g we find the mass of the Cassini Division, excluding the ramp, is 3.1×¹⁰¹⁶ kg while the mass of the Cassini Division ramp, with an opacity of 0.015 cm²/g, is 1.1×¹⁰¹⁷ kg. Assuming a power-law size distribution for the ring particles, the larger opacity of the main Cassini Division is consistent with the largest ring particles there being ∼5 times smaller than the largest particles in the ramp and A ring.     

Self-gravity wake parameters in Saturn’s A and B rings

An increasing body of evidence shows that, at the sub-km level, Saturn’s main A and B rings are dominated by an ever-changing pattern of elongated, canted structures known as self-gravity wakes. Best known for causing azimuthal variations in the rings’ reflectivity, these structures also have a profound influence on how the transmission of the rings varies with both longitude and opening angle, B (Colwell et al. [2006] Geophys. Res. Lett. 33, 7201; Colwell et al. [2007] Icarus 190, 127-144; Hedman et al. [2007] Astron. J. 133, 2624-2629). We use data from three stellar occultations observed by Cassini’s Visual and Infrared Mapping Spectrometer (VIMS) to measure the transmission of the rings as a function of B, when viewed parallel to the wakes. These data are used to constrain properties of the self-gravity wakes as a function of radius across the A and B rings: specifically the fractional width of the gaps between the wakes, G/λ, and the average normal optical depth in the gaps, τ_G. We find that the overall normal optical depth of the rings, τ_n is primarily controlled by G/λ, which varies between <0.05 and ∼0.70 in the A and B rings. The gaps, however, are not completely empty, being filled by material — possibly cm-sized ring particles — with an average normal optical depth which varies from 0.12 to ∼0.4. In addition to regional variations, local variations in τ_G are seen in the regular structure which dominates the inner B ring, and in the environs of strong density waves in the A ring. The same model applied to the lower optical depth Cassini Division reveals very little evidence of self-gravity wakes, except where τ_n exceeds ∼0.25.     

Saturn's Rings: Azimuthal variations,phase curves,and radial profiles in four colors

Four-color photographic photometry of Saturn for the 1977–1979 apparitions has been analyzed to determine the dependence of ring brightness on wavelength, solar phase angle, ring particle orbital phase angle (azimuthal effect), declination of the Earth relative to the ring plane (tilt angle), and radial distance from Saturn. Azimuthal brightness variations up to ±20% relative to the ansae are clearly apparent for the maximum of ring A, but are not detectable for ring B or the outer portion of ring A. The shape of the intensity (I) versus orbital phase angle (θ) curve varies with ring tilt (B) and probably with wavelength, and shows 180° symmetry. As characterized by its slope near the ansae, this curve suggests that the azimuthal effect increases as B decreases from 26 to ≈11°. The phase curves l(α) for the ansae show very little dependence on ring tilt (26° > B > 6°), on wavelength, or on radial distance from Saturn; possibly the curves are somewhat steeper at the smallest tilt angles and for ring A relative to ring B. The radial profile of both rings becomes flatter with decreasing tilt angle and with decreasing wavelength. The latter effect is a natural result of the classical, many-particle-thick ring model.     

Total particulate mass in Enceladus plumes and mass of Saturn’s E ring inferred from Cassini ISS images

The eclipse mosaic (PIA08329) of the Saturn system, taken on September 15, 2006 when Cassini was in Saturn’s shadow, contains numerous color images of the Enceladus plume and the E ring at phase angles ranging from 173° to 179°. These forward-scattering observations sample the diffraction peak for particle radii in the 1–5 μm range. The phase angle dependence and total brightness are sensitive indicators of the total mass of solid material in the plume. We fit the data with a variety of particle shapes and size distributions, and find that the median radius of the equivalent-volume sphere is 3.1 μm, with an uncertainty of ±0.5 μm. The total mass of particles in the plume is (1.45 ± 0.5) × 10⁵ kg. We have not considered variations with altitude in the particle size and shape distribution, and we leave that for another paper. We find that the brightness of the E ring varies with position in the orbit, not only because of the viewing geometry, e.g., variations in phase angle, but also because of some unknown intrinsic variability. The total mass of solid material in the E ring is (12 ± 5.5) × 10⁸ kg. For the plume, the production rate of particles – the mass per unit time leaving the vents is 51 ± 18 kg s⁻¹. We estimate that 9% of these particles are escaping from Enceladus, implying lifetimes of ∼8 years for the E ring particles. Based on three comparisons with vapor amounts from ultraviolet spectroscopy, the ice/vapor ratio is in the range 0.35–0.70. This high ratio poses a problem for theories in which particles form by condensation from the gas phase, and could indicate that particles are formed as spray from a liquid reservoir.     

Cassini imaging of Saturn's rings: II. A wavelet technique for analysis of density waves and other radial structure in the rings

We describe a powerful signal processing method, the continuous wavelet transform, and use it to analyze radial structure in Cassini ISS images of Saturn's rings. Wavelet analysis locally separates signal components in frequency space, causing many structures to become evident that are difficult to observe with the naked eye. Density waves, generated at resonances with saturnian satellites orbiting outside (or within) the rings, are particularly amenable to such analysis. We identify a number of previously unobserved weak waves, and demonstrate the wavelet transform's ability to isolate multiple waves superimposed on top of one another. We also present two wave-like structures that we are unable to conclusively identify. In a multi-step semi-automated process, we recover four parameters from clearly observed weak spiral density waves: the local ring surface density, the local ring viscosity, the precise resonance location (useful for pointing images, and potentially for refining saturnian astrometry), and the wave amplitude (potentially providing new constraints upon the masses of the perturbing moons). Our derived surface densities have less scatter than previous measurements that were derived from stronger non-linear waves, and suggest a gentle linear increase in surface density from the inner to the mid-A Ring. We show that ring viscosity consistently increases from the Cassini Division outward to the Encke Gap. Meaningful upper limits on ring thickness can be placed on the Cassini Division (3.0 m at r∼118,800 km, 4.5 m at r∼120,700 km) and the inner A Ring (10-15 m for r<127,000 km).     

Moonlet induced wakes in planetary rings: Analytical model including eccentric orbits of moon and ring particles

Saturn’s rings host two known moons, Pan and Daphnis, which are massive enough to clear circumferential gaps in the ring around their orbits. Both moons create wake patterns at the gap edges by gravitational deflection of the ring material (Cuzzi, J.N., Scargle, J.D. [1985]. Astrophys. J. 292, 276-290; Showalter, M.R., Cuzzi, J.N., Marouf, E.A., Esposito, L.W. [1986]. Icarus 66, 297-323). New Cassini observations revealed that these wavy edges deviate from the sinusoidal waveform, which one would expect from a theory that assumes a circular orbit of the perturbing moon and neglects particle interactions. Resonant perturbations of the edges by moons outside the ring system, as well as an eccentric orbit of the embedded moon, may partly explain this behavior (Porco, C.C., and 34 colleagues [2005]. Science 307, 1226-1236; Tiscareno, M.S., Burns, J.A., Hedman, M.M., Spitale, J.N., Porco, C.C., Murray, C.D., and the Cassini Imaging team [2005]. Bull. Am. Astron. Soc. 37, 767; Weiss, J.W., Porco, C.C., Tiscareno, M.S., Burns, J.A., Dones, L. [2005]. Bull. Am. Astron. Soc. 37, 767; Weiss, J.W., Porco, C.C., Tiscareno, M.S. [2009]. Astron. J. 138, 272-286). Here we present an extended non-collisional streamline model which accounts for both effects. We describe the resulting variations of the density structure and the modification of the nonlinearity parameter q. Furthermore, an estimate is given for the applicability of the model. We use the streamwire model introduced by Stewart (Stewart, G.R. [1991]. Icarus 94, 436-450) to plot the perturbed ring density at the gap edges.We apply our model to the Keeler gap edges undulated by Daphnis and to a faint ringlet in the Encke gap close to the orbit of Pan. The modulations of the latter ringlet, induced by the perturbations of Pan (Burns, J.A., Hedman, M.M., Tiscareno, M.S., Nicholson, P.D., Streetman, B.J., Colwell, J.E., Showalter, M.R., Murray, C.D., Cuzzi, J.N., Porco, C.C., and the Cassini ISS team [2005]. Bull. Am. Astron. Soc. 37, 766), can be well described by our analytical model. Our analysis yields a Hill radius of Pan of 17.5 km, which is 9% smaller than the value presented by Porco (Porco, C.C., and 34 colleagues [2005]. Science 307, 1226-1236), but fits well to the radial semi-axis of Pan of 17.4 km. This supports the idea that Pan has filled its Hill sphere with accreted material (Porco, C.C., Thomas, P.C., Weiss, J.W., Richardson, D.C. [2007]. Science 318, 1602-1607). A numerical solution of a streamline is used to estimate the parameters of the Daphnis-Keeler gap system, since the close proximity of the gap edge to the moon induces strong perturbations, not allowing an application of the analytic streamline model. We obtain a Hill radius of 5.1 km for Daphnis, an inner edge variation of 8 km, and an eccentricity for Daphnis of 1.5 × 10⁻⁵. The latter two quantities deviate by a factor of two from values gained by direct observations (Jacobson, R.A., Spitale, J., Porco, C.C., Beurle, K., Cooper, N.J., Evans, M.W., Murray, C.D. [2008]. Astron. J. 135, 261-263; Tiscareno, M.S., Burns, J.A., Hedman, M.M., Spitale, J.N., Porco, C.C., Murray, C.D., and the Cassini Imaging team [2005]. Bull. Am. Astron. Soc. 37, 767), which might be attributed to the neglect of particle interactions and vertical motion in our model.     

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.