Infrared observations of Saturn's rings by Cassini CIRS : Phase angle and local time dependence

Since the Saturn orbit insertion (SOI) of the Cassini spacecraft, in July 2004, the Cassini Composite Infrared Spectrometer (CIRS) has obtained a large number of thermal infrared spectra of Saturn's rings. Over the two and a half years of observations to date, ring temperatures were retrieved for a large range of unique geometries, inaccessible from Earth. Understanding their dependencies with phase angle and local time is a clue to understanding the thermal properties and dynamics of Saturn's ring particles.Azimuthal scans of rings, which have been obtained by CIRS at constant radial distance from the planet, have been planned to measure ring temperature variations with local hour angle. Over 47 azimuthal scans for Saturn's main rings (A, B, C and Cassini Division) have been retrieved to date, on both lit and unlit sides, at different phase angles and spacecraft elevations. The first measurements of the transient thermal episode of eclipse cooling in the planetary shadow have also been obtained for all three rings.In this paper, we present an overview of all azimuthal scans obtained by the Cassini/CIRS instrument so far and the dependencies of the temperature and the filling factor with the phase angle and the local hour angle. The ring temperature varies with longitude as the input heating flux coming from Saturn and the Sun changes. The decrease in temperature with the increasing phase angle on both the lit and the unlit sides and for most of the local time also suggests the presence of slowly rotating particles. The crossing of the planet's shadow generates drastic azimuthal variations in temperature, up to 20 K in the C ring. The strong anisotropy of emission observed outside the shadow between low and high phase angles decreases when ring particles cross the shadow, suggesting that particles are almost isothermal in the shadow. This suggests a thermal inertia associated with a rotating rate of particles low enough to have a thermal contrast on their surface.The temperature in the B ring is less sensitive to the phase angle effect on the lit side, suggesting that particles are close enough to form a flat layer at a scale larger than the particle's radius. On the unlit side, particles in the B ring are less sensitive to the lack of solar input than in the C ring or in the A ring. Azimuthal variations of the filling factor in the A ring are also detected with changing ring local time. This effect might be created by the presence of gravitational instabilities (wakes).     

Cassini thermal observations of Saturn's main rings: Implications for particle rotation and vertical mixing

In late 2004 and 2005 the Cassini composite infrared spectrometer (CIRS) obtained spatially resolved thermal infrared radial scans of Saturn's main rings (A, B and C, and Cassini Division) that show ring temperatures decreasing with increasing solar phase angle, α, on both the lit and unlit faces of the ring plane. These temperature differences suggest that Saturn's main rings include a population of ring particles that spin slowly, with a spin period greater than 3.6 h, given their low thermal inertia. The A ring shows the smallest temperature variation with α, and this variation decreases with distance from the planet. This suggests an increasing number of smaller, and/or more rapidly rotating ring particles with more uniform temperatures, resulting perhaps from stirring by the density waves in the outer A ring and/or self-gravity wakes.The temperatures of the A and B rings are correlated with their optical depth, τ, when viewed from the lit face, and anti-correlated when viewed from the unlit face. On the unlit face of the B ring, not only do the lowest temperatures correlate with the largest τ, these temperatures are also the same at both low and high α, suggesting that little sunlight is penetrating these regions.The temperature differential from the lit to the unlit side of the rings is a strong, nearly linear, function of optical depth. This is consistent with the expectation that little sunlight penetrates to the dark side of the densest rings, but also suggests that little vertical mixing of ring particles is taking place in the A and B rings.     

Brightness of Saturn's rings with decreasing solar elevation

Early ground-based and spacecraft observations suggested that the temperature of Saturn's main rings (A, B and C) varied with the solar elevation angle, B^′. Data from the composite infrared spectrometer (CIRS) on board Cassini, which has been in orbit around Saturn for more than five years, confirm this variation and have been used to derive the temperature of the main rings from a wide variety of geometries while B^′ varied from near −24^° to 0^° (Saturn's equinox).Still, an unresolved issue in fully explaining this variation relates to how the ring particles are organized and whether even a simple mono-layer or multi-layer approximation describes this best. We present a set of temperature data of the main rings of Saturn that cover the ∼23^°—range of B^′ angles obtained with CIRS at low (α∼30^°) and high (α≥120^°) phase angles. We focus on particular regions of each ring with a radial extent on their lit and unlit sides. In this broad range of B^′, the data show that the A, B and C rings’ temperatures vary as much as 29-38, 22-34 and 18-23 K, respectively. Interestingly the unlit sides of the rings show important temperature variations with the decrease of B^′ as well. We introduce a simple analytical model based on the well known Froidevaux monolayer approximation and use the ring particles’ albedo as the only free parameter in order to fit and analyze this data and estimate the ring particle's albedo. The model considers that every particle of the ring behaves as a black body and warms up due to the direct energy coming from the Sun as well as the solar energy reflected from the atmosphere of Saturn and on its neighboring particles. Two types of shadowing functions are used. One analytical that is used in the latter model in the case of the three rings and another, numerical, that is applied in the case of the C ring alone. The model lit side albedo values at low phase are 0.59, 0.50 and 0.35-0.38 for the A, B and C rings, respectively.     

Saturn's rings: Infrared brightness variation with solar elevation

The variation in infrared equilibrium brightness temperature of Saturn's A, B, and C rings is modeled as a function of solar elevation B′ with respect to the ring plane. The basic model includes estimates of minimum and maximum interparticle shadowing in a monolayer approximation. Simple laboratory observations of random particle distributions at various illumination angles provide more realistic shadowing functions. Radiation balance calculations yield the physical (kinetic) temperature of particles in equilibrium with radiation from the Sun, Saturn, and neighboring particles. Infrared brightness temperatures as a function of B′ are then computed and compared to the available 20-μm data (Pioneer results are also briefly discussed). The A and B rings are well modeled by an optically thick monolayer, or equivalently, a flat sheet, radiating from one side only. This points to a temperature contrast between the two sides, possibly due to particles with low thermal inertia. Other existing models for the B ring are discussed. The good fit for the monolayer model does not rule out the possibility that the A and B rings are many particles thick. It could well be that a multilayer ring produces an infrared behavior (as a function of tilt angle) similar to that of a monolayer. The C ring brightness increases as B′ decreases. This contrast in behavior can be understood simply in terms of the low C ring optical depth and small amount of interparticle shadowing. High-albedo particles (A?0.5) can fit the C ring infrared data if they radiate mostly from one hemisphere due to slow rotation or low thermal inertia (or both). Alternatively, particles isothermal over their surface (owing to a rapid spin, high inertia, or small size), and significantly darker (A?0.3) than the A and B ring particles, can produce a similar brightness variation with ring inclination. In any case, the C ring particles have significantly hotter physical temperatures than the particles in the A and B rings, whether or not the rings form a monolayer.     

Spinning particles in Saturn's C ring from mid-infrared observations: Pre-Cassini mission results

Saturn's C ring thermal emission has been observed in mid-infrared wavelengths, at three different epochs and solar phase angles, using ground based instruments (CFHT in 1999 and VLT/ESO in 2005) and the Infrared Radiometer Instrument Spectrometer (IRIS) onboard the Voyager 1 spacecraft in 1980. Azimuthal variations of temperature in the C ring's inner region, observed at several phase angles, have been analyzed using our new standard thermal model [Ferrari, C., Leyrat, C., 2006. Astron. Astrophys. 447, 745-760]. This model provides predicted ring temperatures for a monolayer ring composed of spinning icy spherical particles. We confirm the very low thermal inertia (on the order of 10 ) found previously by Ferrari et al. [Ferrari, C., Galdemard, P., Lagage, P.O., Pantin E., Quoirin, C., 2005. Astron. Astrophys. 441, 379-389] that reveals the very porous regolith at the surface of ring particles. We are able to explain both azimuthal variations of temperature and the strong asymmetry of the emission function between low and high phase angles. We show that large particles spinning almost synchronously might be present in the C ring to explain differences of temperature observed between low and high phase angle. Their cross section might represent about 45% of the total cross section. However, their numerical fraction is estimated to only ∼0.1% of all particles. Thermal behavior of other particles can be modeled as isothermal behavior. This work provides an indirect estimation of the particle's rotation rate in Saturn's rings from observations.     

Sputtering of ice grains and icy satellites in Saturn's inner magnetosphere

Icy grains and satellites orbiting in Saturn's magnetosphere are immersed in a plasma that sputters their surfaces. This limits the lifetime of the E-ring grains and ejects neutrals that orbit Saturn until they are ionized and populate its magnetosphere. Here we re-evaluate the sputtering rate of ice in Saturn's inner magnetosphere using the recent Cassini data on the plasma ion density, temperature and composition [Sittler Jr., E.C., et al., 2007a. Ion and neutral sources and sinks within Saturn's inner magnetosphere: Cassini results. Planet. Space Sci. 56, 3-18.] and a recent summary of the relevant sputtering data for ice [Famá, M., Shi, J., Baragiola, R.A., 2008. Sputtering of ice by low-energy ions. Surf. Sci. 602, 156-161.]. Although the energetic (>10 keV) ion component at Saturn is much smaller than was assumed to be the case after Voyager [Jurac, S., Johnson, R.E., Richardson, J.D., Paranicas, C., 2001a. Satellite sputtering in Saturn's magnetosphere. Planet. Space Sci. 49, 319-326; Jurac, S., Johnson, R.E., Richardson, J.D., 2001b. Saturn's E ring and production of the neutral torus. Icarus 149, 384-396.], we show that the sputtering rates are sensitive to the temperature of the thermal plasma and are still robust, so that sputtering likely determines the lifetime of the grains in Saturn's tenuous E-ring.     

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.     

A multilayer model for thermal infrared emission of Saturn’s rings II: Albedo, spins, and vertical mixing of ring particles inferred from Cassini CIRS

Since the Saturn orbit insertion of the Cassini spacecraft in mid-2004, the Cassini composite infrared spectrometer (CIRS) measured temperatures of Saturn’s main rings at various observational geometries. In the present study, we apply our new thermal model (Morishima, R., Salo, H., Ohtsuki, K. [2009]. Icarus 201, 634-654) for fitting to the early phase Cassini data (Spilker, L.J., and 11 colleagues [2006]. Planet. Space Sci. 54, 1167-1176). Our model is based on classical radiative transfer and takes into account the heat transport due to particle motion in the azimuthal and vertical directions. The model assumes a bimodal size distribution consisting of small fast rotators and large slow rotators. We estimated the bolometric Bond albedo, A_V, the fraction of fast rotators in cross section, f_fast, and the thermal inertia, Γ, by the data fitting at every radius from the inner C ring to the outer A ring. The albedo A_V is 0.1-0.4, 0.5-0.7, 0.4, 0.5 for the C ring, the B ring, the Cassini division, and the A ring, respectively. The fraction f_fast depends on the ratio of scale height of fast rotators to that of slow rotators, h_r. When h_r = 1, f_fast is roughly half for the entire rings, except for the A ring, where f_fast increases from 0.5 to 0.9 with increasing saturnocentric radius. When h_r increases from 1 to 3, f_fast decreases by 0.2-0.4 for the B and A rings while no change in f_fast is seen for the optically thin C ring and Cassini division. The large f_fast seen in the outer A ring probably indicates that a large number of small particles detach from large particles in high velocity collisions due to satellite perturbations or self-gravity wakes. The thermal inertia, Γ, is constrained from the efficiency of the vertical heat transport due to particle motion between the lit and unlit faces, and is coupled with the type of vertical motion. We found that in most regions, except for the mid B ring, sinusoidal vertical motion without bouncing is more reasonable than cycloidal motion assuming bouncing at the midplane, because the latter motion gives too large Γ as compared with previous estimations. For the mid B ring, where the optical depth is highest in Saturn’s rings, cycloidal vertical motion is more reasonable than sinusoidal vertical motion which gives too small Γ.     

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].     

Saturn's rings in the thermal infrared

This paper reviews our current knowledge of Saturn's rings’ physical properties as derived from thermal infrared observations. Ring particle composition, surface structure and spin as well as the vertical structure of the main rings can be determined. These properties are the key to understand the origin and evolution of Saturn's rings. Ring composition is mainly constrained by observations in the near-infrared but the signature of some probable contaminants present in water ice may also be found at mid-infrared wavelengths. The absence of the silicate signature limits nowadays their mass fraction to 10^−7±1. Recent measurements on the thermal inertia of the ring particle surface show it is very low, of the order of 5±2 Jm⁻² K⁻¹ s^−1/2. New models and observations of the complete crossing of the planetary shadow are needed to attribute this low value either to compact regoliths covered by cracks due to collisions and thermal stresses or to large fluffy and irregular surfaces. Studies of the energy balance of ring particles show a preference for slowly spinning particles in the main rings. Supplementary observations at different phase angles, showing the temperature contrast between night and day sides of particles, and new models including finite spin and thermal inertia, are needed to constrain the actual spin distribution of ring particles. These results can then be compared to numerical simulations of ring dynamics. Many thermal models have been proposed to reproduce observations of the main rings, including alternative mono- or many-particles-thick layers or vertical heterogeneity, with no definitive answer. Observations on the lit and dark faces of rings as a function of longitude, at many incidence and emission angles, would provide prime information on the vertical thermal gradient due to interparticle shadowing from which constraints on the local vertical structure and dynamics can be produced. Future missions such as Cassini will provide new information to further constrain the ring thermal models.     

Thermal response of Saturn's ring particles during and after eclipse

We present infrared (20 μm) observations of Saturn's rings for a solar elevation angle of 10° and phase angle of 6°. Scans across the rings yield information about the cooling of particles during eclipse and the subsequent heating along their orbits. All three rings exhibit significant cooling during eclipse, as well as a 20-μm brightness asymmetry between east and west ansae, the largest asymmetry occuring in the C ring (the brightest ring). The eclipse cooling is a simple and adequate explanation for 20-μm brightness asymmetries between the ansae of Saturn's rings. The relatively large C ring asymmetry is thought to be primarily due to the short travel time of the particles in that ring from eclipse exit to east ansa. We compare the B ring data to the theoretical models of H.H. Aumann and H.H. Kieffer (1973, Astrophys. J.186, 305–311) in order to set constraints on the average particle size and thermal inertia. The rather rapid heating after exit from eclipse points to low-conductivity-particle surfaces, similar to the water frost surfaces of Galilean satellites. If the surface conductivity is indeed low, one cannot determine an upper limit for the particle size through such infrared observations, since only the uppermost millimeters experience a thermal response during eclipse. However, based on these infrared data alone, it is clear that particles of radius equal to a few millimeters or less cannot occupy a significant fraction of the ring surface area, because-regardless of thermal inertia-their thermal response is much faster than observed.     

Simulations of dense planetary rings: IV. Spinning self-gravitating particles with size distribution

Previous self-gravitating simulations of dense planetary rings are extended to include particle spins. Both identical particles as well as systems with a modest range of particle sizes are examined. For a ring of identical particles, we find that mutual impact velocity is always close to the escape velocity of the particles, even if the total rms velocity dispersion of the system is much larger, due to collective motions associated to wakes induced by near-gravitational instability or by viscous overstability. As a result, the spin velocity (i.e., the product of the particle radius and the spin frequency) maintained by mutual impacts is also of the order of the escape velocity, provided that friction is significant. For the size distribution case, smaller particles have larger impact velocities and thus larger spin velocities, particularly in optically thick rings, since small particles move rather freely between wakes. Nevertheless, the maximum ratio of spin velocities between the smallest and largest particles, as well as the ratio for translational velocities, stays below about 5 regardless of the width of the size distribution. Particle spin state is one of the important factors affecting the temperature difference between the lit and unlit face of Saturn's rings. Our results suggest that, to good accuracy, the spin frequency is inversely proportional to the particle size. Therefore, the mixing ratio of fast rotators to slow rotators on the scale of the thermal relaxation time increases with the width of the particle size distribution. This will offer means to constrain the particle size distribution with the systematic thermal infrared observations carried by the Cassini probe.     

Moonlets and clumps in Saturn's F ring

Cassini UVIS star occultations by the F ring detect 13 events ranging from 27 m to 9 km in width. We interpret these structures as likely temporary aggregations of multiple smaller objects, which result from the balance between fragmentation and accretion processes. One of these features was simultaneously observed by VIMS. There is evidence that this feature is elongated in azimuth. Some features show sharp edges. At least one F ring object is opaque and may be a “moonlet.” This possible moonlet provides evidence for larger objects embedded in Saturn's F ring, which were predicted as the sources of the F ring material by Cuzzi and Burns [Cuzzi, J.N., Burns, J.A., 1988. Icarus 74, 284-324], and as an outcome of tidally modified accretion by Barbara and Esposito [Barbara, J.M., Esposito, L.W., 2002. Icarus 160, 161-171]. We see too few events to confirm the bi-modal distribution which Barbara and Esposito [Barbara, J.M., Esposito, L.W., 2002. Icarus 160, 161-171] predict. These F ring structures and other youthful features detected by Cassini may result from ongoing destruction of small parent bodies in the rings and subsequent aggregation of the fragments. If so, the temporary aggregates are 10 times more abundant than the solid objects. If recycling by re-accretion is significant, the rings could be quite ancient, and likely to persist far into the future.     

Exploring local N-body simulations of Saturn's rings

The dynamical behavior of low and moderately high optical depth regions of Saturn's ring system of discrete, mutually gravitating, and inelastically colliding particles is studied by simplified local N-body simulations in Hill's linearized equations context. The focus is on a statistical analysis of time-evolution of fine-scale structures seen in the simulations and the comparison between theoretical predictions and computer experiments. Prospects for the Cassini spacecraft mission are briefly summarized.     

A procedure to analyze nonlinear density waves in Saturn's rings using several occultation profiles

Cassini radio science experiments have provided multiple occultation optical depth profiles of Saturn's rings that can be used in combination to analyze density waves. This paper establishes an accurate procedure of inversion of the wave profiles to reconstruct the wave kinematic parameters as a function of semi-major axis, in the nonlinear regime. This procedure is established using simulated data in the presence of realistic noise perturbations, to control the reconstruction error. It is then applied to the Mimas 5:3 density wave. There are two important concepts at the basis of this procedure. The first one is that it uses the nonlinear representation of density waves, and the second one is that it relies on a combination of optical depth profiles instead of just one profile. A related method to analyze density waves was devised by Longaretti and Borderies [Longaretti, P.-Y., Borderies, N., 1986. Icarus 67, 211-223] to study the nonlinear density wave associated with the Mimas 5:3 resonance, but the single photopolarimetric profile provided limited constraints. Other studies of density waves analyzing Cassini data [Colwell, J.E., Esposito, L.W., 2007. Bull. Am. Astron. Soc. 39, 461; Tiscareno, M.S., Burns, J.A., Nicholson, P.D., Hedman, M.M., Porco, C.C., 2007. Icarus 189, 14-34] are based on the linear theory and find inconsistent results from profile to profile. Multiple cuts of the rings are helpful in a fundamental way to ensure the accuracy of the procedure by forcing consistency among the various optical depth profiles. By way of illustration we have applied our procedure to the Mimas 5:3 density wave. We were able to recover precisely the kinematic parameters from the radio experiment occultation data in most of the propagation region; a preliminary analysis of the pressure-corrected dispersion allowed us to determine new but still uncertain values for the opacity (K?0.02 cm²/g) and velocity dispersion of (c₀?0.6 cm/s) in the wave region. Our procedure constitutes the first step in our planned analysis of the density waves of Saturn's rings. It is very accurate and efficient in the far-wave region. However, improvements are required within the first wavelength. The ways in which this method can be used to establish diagnostics of ring physics are outlined.     

Radar imaging of Saturn's rings

We present delay-Doppler images of Saturn's rings based on radar observations made at Arecibo Observatory between 1999 and 2003, at a wavelength of 12.6 cm and at ring opening angles of 20.1°?|B|?26.7°. The average radar cross-section of the A ring is ∼77% relative to that of the B ring, while a stringent upper limit of 3% is placed on the cross-section of the C ring and 9% on that of the Cassini Division. These results are consistent with those obtained by Ostro et al. [1982, Icarus 49, 367-381] from radar observations at |B|=21.4°, but provide higher resolution maps of the rings' reflectivity profile. The average cross-section of the A and B rings, normalized by their projected unblocked area, is found to have decreased from 1.25±0.31 to 0.74±0.19 as the rings have opened up, while the circular polarization ratio has increased from 0.64±0.06 to 0.77±0.06. The steep decrease in cross-section is at variance with previous radar measurements [Ostro et al., 1980, Icarus 41, 381-388], and neither this nor the polarization variations are easily understood within the framework of either classical, many-particle-thick or monolayer ring models. One possible explanation involves vertical size segregation in the rings, whereby observations at larger elevation angles which see deeper into the rings preferentially see the larger particles concentrated near the rings' mid-plane. These larger particles may be less reflective and/or rougher and thus more depolarizing than the smaller ones. Images from all four years show a strong m=2 azimuthal asymmetry in the reflectivity of the A ring, with an amplitude of ±20% and minima at longitudes of 67±4° and 247±4° from the sub-Earth point. We attribute the asymmetry to the presence of gravitational wakes in the A ring as invoked by Colombo et al. [1976, Nature 264, 344-345] to explain the similar asymmetry long seen at optical wavelengths. A simple radiative transfer model suggests that the enhancement of the azimuthal asymmetry in the radar images compared with that seen at optical wavelengths is due to the forward-scattering behavior of icy ring particles at decimeter wavelengths. A much weaker azimuthal asymmetry with a similar orientation may be present in the B ring.     

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.     

The opposition effect in the outer Solar system: A comparative study of the phase function morphology

In this paper, we characterize the morphology of the disk-integrated phase functions of satellites and rings around the giant planets of our solar system. We find that the shape of the phase function is accurately represented by a logarithmic model [Bobrov, M.S., 1970. Physical properties of Saturn's rings. In: Dollfus, A. (Ed.), Surfaces and Interiors of Planets and Satellites. Academic, New York, pp. 376-461]. For practical purposes, we also parametrize the phase curves by a linear-exponential model [Kaasalainen, S., Muinonen, K., Piironen, J., 2001. Comparative study on opposition effect of icy solar system objects. Journal of Quantitative Spectroscopy and Radiative Transfer 70, 529-543] and a simple linear-by-parts model [Lumme, K., Irvine, W.M., 1976. Photometry of Saturn's rings. Astronomical Journal 81, 865-893], which provides three morphological parameters: the amplitude A and the half-width at half-maximum (HWHM) of the opposition surge, and the slope S of the linear part of the phase function at larger phase angles.Our analysis demonstrates that all of these morphological parameters are correlated with the single-scattering albedos of the surfaces.By taking more accurately into consideration the finite angular size of the Sun, we find that the Galilean, Saturnian, Uranian and Neptunian satellites have similar HWHMs (?0.5^°), whereas they have a wide range of amplitudes A. The Moon has the largest HWHM (∼2^°). We interpret that as a consequence of the “solar size bias”, via the finite angular size of the Sun which varies dramatically from the Earth to Neptune. By applying a new method that attempts to morphologically deconvolve the phase function to the solar angular size, we find that icy and young surfaces, with active resurfacing, have the smallest values of A and HWHM, whereas dark objects (and perhaps older surfaces) such as the Moon, Nereid and Saturn's C ring have the largest A and HWHM.Comparison between multiple objects also shows that solar system objects belonging to the same planet have comparable opposition surges. This can be interpreted as a “planetary environmental effect” that acts to locally modify the regolith and the surface properties of objects which are in the same environment.     

A multilayer model for thermal infrared emission of Saturn’s rings. III: Thermal inertia inferred from Cassini CIRS

The thermal inertia values of Saturn’s main rings (the A, B, and C rings and the Cassini division) are derived by applying our thermal model to azimuthally scanned spectra taken by the Cassini Composite Infrared Spectrometer (CIRS). Model fits show the thermal inertia of ring particles to be 16, 13, 20, and 11 J m⁻² K⁻¹ s^−1/2 for the A, B, and C rings, and the Cassini division, respectively. However, there are systematic deviations between modeled and observed temperatures in Saturn’s shadow depending on solar phase angle, and these deviations indicate that the apparent thermal inertia increases with solar phase angle. This dependence is likely to be explained if large slowly spinning particles have lower thermal inertia values than those for small fast spinning particles because the thermal emission of slow rotators is relatively stronger than that of fast rotators at low phase and vise versa. Additional parameter fits, which assume that slow and fast rotators have different thermal inertia values, show the derived thermal inertia values of slow (fast) rotators to be 8 (77), 8 (27), 9 (34), 5 (55) J m⁻² K⁻¹ s^−1/2 for the A, B, and C rings, and the Cassini division, respectively. The values for fast rotators are still much smaller than those for solid ice with no porosity. Thus, fast rotators are likely to have surface regolith layers, but these may not be as fluffy as those for slow rotators, probably because the capability of holding regolith particles is limited for fast rotators due to the strong centrifugal force on surfaces of fast rotators. Other additional parameter fits, in which radii of fast rotators are varied, indicate that particles less than ∼1 cm should not occupy more than roughly a half of the cross section for the A, B, and C rings.