A&A 441, 379-389 (2005)
C. Ferrari1 - P. Galdemard1 - P. O. Lagage1 - E. Pantin1 - C. Quoirin2
1 - Laboratoire AIM, Unité Mixte de Recherche No. 7158 CEA-CNRS-Université Paris 7, France
2 - Université Paris 6 - Pierre et Marie Curie, France
Received 21 March 2005 / Accepted 11 May 2005
Thermal inertias of Saturn's B and C ring particles have been derived from infrared observations using the CAMIRAS camera mounted on the Canada-France-Hawaii Telescope. They are respectively and . Such low values might be characteristic of a frosty and porous regolith fractured by cracks or of very porous particle aggregates. Particles have to be slowly spinning to explain the observed ring temperatures. A large azimuthal asymmetry with an amplitude about 1 K is detected on the West ansa of the B ring. It cannot be explained by a model that considers the ring as a slab of low thermal inertia rapidly warming up to the sunlight after its eclipse into the planetary shadow.
Key words: planetary rings - infrared: solar system - conduction - radiation mechanism: thermal
Ring particles cool down when they cross the planetary shadow and heat up when back in the sunlight. During this transient thermal regime, the surface temperature changes at a rate which depends on its thermal inertia. Compact water ice surfaces of high thermal inertia will suffer only small temperature variations. Fluffy, porous, or frosty regoliths will prevent heat propagation into the particle or will hardly be able to store heat. Temperature changes will be large due to their low thermal inertia.
In order to constrain the thermal inertia and the particle size distribution in Saturn's rings from thermal infrared observations, Aumann and Kieffer (1973) proposed a model of heating and cooling of ring particles. Morrison (1974) gave a lower limit of 40 for the thermal inertia in the B ring, assuming that the particles are larger than 2 cm and using Aumann and Kieffer model figures. Froidevaux et al. (1981) have measured a gradient in brightness temperature of for the B ring and for the C ring between entry and exit from shadow. They have also observed how the brightness temperature of the B ring varies along the ring (Fig. 6). Qualitative comparison with Aumann's model has led to rejecting a thermal inertia as high as that of solid water ice (2600 ). A thermal inertia of 13 , measured at the same epoch for the Galilean satellites, has been favored instead, with the condition that ring particles are larger than 1 cm. This is in agreement with Morrison's conclusions on the size limit but supports the idea of a lower value for the thermal inertia thanks to more accurate measurements. The Pioneer infrared radiometer was the first instrument onboard a spacecraft to provide data on ring temperatures in the shadow itself. Froidevaux and Ingersoll (1980) reported temperature drops of about 10-15 K in the A and C rings confirming the idea of a low thermal inertia. Observations of C ring shadow were also performed with the IRIS Voyager infrared spectrometer (Hanel et al. 1981), but no estimate of thermal inertia has been derived.
With such a low thermal inertia, particles should have recovered their equilibrium temperature at the East ansa, a quarter of an orbit after midnight, located at the middle of the planetary shadow. Many observers in the seventies and eighties have not observed any temperature contrast between the East and far West ansae in the A, B, and C rings, with 2 K uncertainty on measurements (Murphy 1972; Nolt et al. 1980; Tokunaga et al. 1980). Nolt et al. (1978) and Froidevaux et al. (1980) did not detect any gradient between ansae except in the C ring, with a values of and , respectively. This difference has been explained simply by the fact that C ring particles had spent less time back in sunlight after their exit from eclipse (Froidevaux et al. 1980). However some observers (Allen et al. 1971; Murphy et al. 1972; Morrison 1974) have mentioned the existence of such a gradient between both ansae, which they relate either to the cooling of particles in the planetary shadow or to a variation in albedo on the particle surface, the leading side being brighter than the trailing side. There is certainly a thermal gradient between exit from shadow and the East morning ansa for the B and C rings, while the temperature asymmetry between East and West ansae remains uncertain for A and B rings. More recently, thermal IR spectroscopy of Saturn's A and B rings has not detected any temperature gradient between East and West ansae (Lynch et al. 2000), confirming that most of the heating of ring particles in these rings happens before they reach the morning East ansa.
Our objective is to derive quantitative estimates of the thermal inertia of Saturn's main rings. We mounted the CAMIRAS instrument built at the SAp/DAPNIA/DSM/CEA Saclay laboratory on the CFH Telescope to obtain images of Saturn's rings in the mid-infrared domain and to measure full longitudinal temperature profiles of the rings. We then built a simple thermal model to derive the thermal inertia from the observed heating period on the morning East ansa. In this paper, we present the observations performed in July 1999 and March 2000 in Sect. 2. In Sect. 3, a model of particles heating and cooling is proposed to constrain the thermal inertia of ring particles. The thermal inertias for both B and C rings are derived and discussed in Sect. 4. Conclusions and future work are in Sect. 5.
|Figure 1: CAMIRAS images of Saturn and its rings at 20.5 . The effective spatial resolution of 1.4 is plotted on the upper right corner. Ring particles motion is clockwise. The East ansa of the ring and the exit of the shadow are on the left side of the planet.|
|Open with DEXTER|
|Figure 2: Radial profiles of Saturn's rings in July 1999 and March 2000. Brightness temperature is given as a function of the distance to Saturn center. Profiles taken on the East and West ansae are overlapping for a better comparison. Error bars represent () average uncertainty over the rings. In March 2000, the East ansa is slightly truncated at the edge of the detector array. Vertical dashed lines represent Saturn radius ( ) and boundaries of main rings.|
|Open with DEXTER|
A, B, and C ring temperatures on the ansae were 82 K, 88.5 K, and 86.5 K in July respectively and 81 K, 89.5 K, and 86 K in March (Fig. 2).
The average uncertainty is about 0.3 K for 1
To be compared with previously cited observations, which are normalized to a Sun-Saturn distance of 9.25 AU,
temperatures have to be corrected by
The temperatures of A and B rings are in very good agreement with expectations of both monolayer and multilayer thermal models for a solar elevation around
(Froidevaux 1981; Kawata 1983).
A and B ring normalized temperatures increase with solar elevation,
which is compatible with the idea that the physical temperature increases as the solar flux penetrates deeper into the ring layers.
The C ring exhibits significantly warmer normalized temperatures (87 K in July 1999 and 88.5 K in March 2000)
compared to the value of 82 K expected by models and previous observations.
But they remain compatible with the Nolt et al. (1978) measurement of 87.7 K in the
C ring at solar elevation
They also slightly increase with solar elevation contrary to the currently accepted picture that this ring is thin enough
for the inter-particle shadowing effect to be negligible at these solar elevations.
The elevation of the observer is almost constant, but the solar elevation has increased by
A slight increase in the brightness temperature of the rings might then be expected in March 2000,
when solar rays enter deeper into Saturn's rings and heat up ring particles more efficiently.
But Saturn is also farther away from the Sun by 0.44 UA, which somewhat counterbalances this effect
so that comparable fluxes are expected in both data sets.
No East-West asymmetry were detected in the March 2000 observations (Fig. 2).
In July data, a 0.7 K (0.8%) asymmetry systematically appeared in each ring but it remained within the error bars
(0.3% for astrometry and 0.2% for
noise level, i.e. 0.3 K).
We cannot confirm the existence of a thermal gradient between East and West ansae in the C ring as seen in previous data.
|Figure 3: Azimuthal brightness temperature profiles of B and C rings in July 1999 and March 2000. Longitude is located at midnight in the middle of the shadow. The East (morning) ansa is at longitude , noon at and West (afternoon) ansa at longitude . Error bars represent uncertainty. Vertical dashed lines represent the shadow limits on the rings. The shadow extent does not change significantly between epochs at the longitudinal scale of this figure. The ellipsoid shape of Saturn was taken into account.|
|Open with DEXTER|
Azimuthal profiles of B and C ring are displayed in Fig. 3. Points were extracted every 10 degrees in the longitudinal direction inside the B ring and every 12 degrees in the C ring to be independent, given the limited spatial resolution of 1.4 . Because the A ring is most sensitive to IR flat-fielding and optical distortions, we were not able to extract an azimuthal profile with reasonable error bars. The Earth was near maximum elongation on the West ansa in July 1999 so that the entry into the shadow could be seen. In March 2000, the exit of the shadow could be observed despite the moderate East elongation of the Earth (Fig. 3). Data around noon were not plotted because they were contaminated with Saturn infrared emission. In the B ring, the temperature gradient between the exit from the shadow (longitude ) and the East ansa (longitude ) was in March 2000. It was about in July as the first five minutes of heating cannot be observed. Particles had recovered most of their initial temperature at East ansa. In the C ring, a maximum gradient of was observed between exit from the shadow and East ansa in March 2000 (Fig. 3). It fell down to in July 1999. The measurements at the exit of the shadow were certainly contaminated by the planet due to the limited spatial resolution.
A slight azimuthal asymmetry of could be seen on the West ansa (longitude -140 to -50) of the B ring in July, its amplitude was about in March 2000. It peaked at in July and at in March. After the West ansa (longitude ), a ground-based observer would be looking at the face of a spherical particle, which is illuminated both by direct sunlight and planet infrared radiation at the same time. If slowly spinning with low thermal inertia, the brightness temperature of the sphere will increase regularly with azimuth and peak just before entry into the shadow. This suggestion meant to explain the observed azimuthal asymmetry has to be more widely explored by including spherical spinning particles in our model. No significant azimuthal asymmetry other than the one due to the heating after the eclipse could be observed in the C ring.
The brightness temperature of a ring differs from the particle temperature due to interparticle shadowing and ring-filling factor. Aumann & Kieffer (1973) have related both temperatures with an anisotropy index Q, such that , to take the emission function of individual particles and the filling factor into account. This model implicitly assumes a monolayer structure of the ring. In the case of very low optical depth, Q=Q(particle) and Q=1 for very large optical depth. The anisotropy of the emission function of ring particles depends on their spin rate, axis, and thermal inertia.
Table 1: Definitions of variables.
In the present paper, the ring is also assumed to spread in a monolayer as described by Froidevaux (1981).
The brightness temperature is related to the particle temperature by:
A particle at distance a from Saturn center receives energy from different sources, i.e. the Sun, Saturn, and the nearby particles.
The energy balance, away from the shadow is:
|Figure 4: Absorbed fluxes of B and C rings as a function of ring longitude and heating sources for Bond albedo and f=2 at July 1999 epoch. (- ) F1, (- -) F1+F3, (- -) F1+F3+F4, (...) F1+F3+F4+infrared neighbours, (solid line) F1+F2+F3+F4 +infrared neighbours. Longitude is zero at midnight; West ansa (W) is located at and East ansa (E) at .|
|Open with DEXTER|
The inter-particle shadowing is represented in this model by the fractional area of the particle visible from elevation B' above the ring plane. It is dissociated from the anisotropy of the particle emission function, which is described by rotation factor f. Formulation of has been provided by Froidevaux for two viewing geometries of a monolayer of ring particles regularly distributed along a reseau and producing minimum and maximum shadowing factors. We prefered, however, to use his experimental curve of derived from laboratory measurements of the shadowed fractional area in random configurations of spherical balls, mentioned as the random shadowing case.
Ground-based observations of Saturn are limited to small phase angles. In this viewing geometry, only the day side of particles is observed. It corresponds to the warmest point of slowly rotating particles (assuming the thermal inertia is small), or to the average temperature for the faster ones, then isothermal (assuming their pole of rotation is not facing the Sun at any time). The anisotropy of the emission function is described here by rotation factor f. F=2 for a slow rotating particle, which emits mainly over one hemisphere, and f=4 for a fast rotating isothermal particle, which emits over sr. All particles are then supposed to have the same spinning behavior in this ring, an assumption which might be true in Saturn's rings only for the population of large particles. Centimeter-sized particles are gravitationally scattered by the larger ones and show a large dispersion in spin rate (Salo 1987; Richardson 1994). The total solid angle subtended by the eight closest neighbours is as the blocking area is . This assumes that the four closest particles are located at distance 2d and that the four left are at distance 2 . Average optical depths of the B and C rings are 1.4 and 0.1, respectively (Hanel et al. 1981; Spilker et al. 2002; Esposito et al. 1984).
Contribution of the different heating sources varies with the ring distances and optical depth (Fig. 4). The main heating source is, of course, the Sun (F1). The particles in the thiner C ring are warmer than the B ring because the absorbed flux is greater due to weaker inter-particle shadowing. The infrared contribution from the planet on ring particles (F3) is larger in the C ring as it is closer to the planet. Including infrared and visible contributions from the closest particles ( and F4) significantly raises the temperature in the B ring compared to the C ring, as particles are closer together in that thick ring. The infrared contribution of the neighbours is progressively decreasing in the shadow as particles are cooling together. The way this flux decreases in the shadow is unknown a priori as the surface temperature of particles has not been solved yet.
The contribution of Saturn's visible hemisphere F2 is larger in the C ring too. C ring particles absorb much more flux than B ring particles due to their lower distance from the planet and the limited mutual shadowing. This creates an azimuthal asymmetry in the rings which peaks at noon (Fig. 4). But the shape of this azimuthal asymmetry is model-dependent. Our model indeed computes the temperature of a flat ring element, which receives all the contributions at the same time, but this is not the case for a real spherical particle. At noon, F2 is absorbed by the face of the sphere, which is invisible to a ground-based observer. Then, in the case of a slow spinning particle, the present asymmetry due to F2 can never be seen in fact. If the particle is rapidly spinning in comparison with the orbital period and is of moderate thermal inertia, this energy input may be visible just after noon by the observer. At other longitudes, the F2 contribution becomes negligible. The flux F3 suffers the same problem but does not induce an azimuthal asymmetry here, as it is constant along the azimuth. Its contribution in the shadow is very important for both rings and far from negligible in the C ring. The F2 contribution of the Saturn visible hemisphere will be neglected from now on, as it is quite negligible and produces an artificial azimuthal asymmetry in flux (and then in temperature), which is not real. It is confirmed in the fit between data and model as the residuals are much worse when F2 is included. F3 will be included as its contribution is not negligible on average along the orbit. A new model which takes into account the spherical shape of particles, their spin rate and orientation should be built to bypass these drawbacks.
is the global physical property which controls both the temperature variations
at the surface of ring particles and the capacity of the sub-surface to store heat.
It mainly depends on thermal conductivity K, as volume density
and specific heat
vary less with the nature of the material.
The boundary condition for the heat flow at the surface is:
|Figure 5: Azimuthal variation of the surface temperature as a function of thermal inertia ( ). The particle is located at a=84 000 km from Saturn, , f=2 and . Calculations are made for the viewing geometry of the July 1999 epoch. In this case, and . All incident fluxes but F2 have been included.|
|Open with DEXTER|
The specific heat is the amount of heat necessary to increase the temperature of of solid material in 1 s by 1 K. It characterizes the ability of the material to store heat. Experimental studies by Klug et al. (1991) show that it increases only by a factor of about 1.7 between the Ih crystalline phase of water ice and its low density amorphous phase, at a ring temperature of 100 K. It also depends on temperature. For crystalline water ice, it follows the law where T is the temperature of the material (Klinger 1981). As temperature might significantly differ from B to C ring due to different inter-particle shadowing, we have introduced this temperature dependence in our model. However, instead of using the inside temperature T, the average surface temperature over a period is used and . Its value is then determined for an epoch of observation and a ring.
The thermal conductivity K is the amount of heat transfered in 1s due to a thermal gradient of 1 K/m in the material. In other words, it controls the heat flow, given a thermal gradient according to Fourier's law. K depends on the composition of the solid phase and more importantly on its structure. The more there are areas of contact between grains of the regolith, for example, the better the heat can propagate by conduction and the higher K. The thermal conductivity may then decrease because of fractures in a compacted regolith or because of porosity.
Particles are supposed to be much larger than the thermal skin depth, i.e. . Saturn's ring particles are known to be cm-to-m sized with a continuous power-law size distribution , where on average and with minimum and maximum sizes of about 3-30 cm and 10-20 m, respectively, in the C and B rings, respectively (French & Nicholson 2000, and references therein). Smaller dust particles are absent (Dones et al. 1993; Doyle et al. 1989; Cooke 1991). The validity of this assumption is checked a posteriori.
The diffusion Eq. (3) inside the particle is solved with a Crank-Nicolson algorithm. This finite difference implicit scheme, which is second order accurate in both space and time, has been chosen to insure stability and convergence for a wide range of Fourier numbers, even if it may be slower than any other simpler scheme. The boundary condition (4) is non-linear in and has been linearized. The surface temperature curve on longitude is re-calculated until the average emitted flux equals to the average absorbed flux. The tolerance on this criterion is 10-4. The relative accuracy of equality is about 10-5 for the values of the thermal inertia used in this paper.
|Figure 6: Best fit of the B ring data of Froidevaux et al. (1981) obtained with our eclipse model for f=2. Thermal inertia is and Bond albedo at 3. At the epoch , and the phase angle is on the East ansa. "Max'' and "Min'' show the temperature for the limit values and .|
|Open with DEXTER|
|Figure 7: Brightness temperature azimuthal profile of the B ring at range a=105 000 km ( ) in July 1999 and March 2000. Best fits (3) are represented in case of: (E) fit with east ansa data only, (E+W) fit with east and west ansae data.|
|Open with DEXTER|
|Figure 8: Best fits of brightness temperature azimuthal profile of the C ring at range a=84 000 km ( ) in July 1999 and March 200. See Fig. 7 for details.|
|Open with DEXTER|
Our eclipse model was fitted first to Froidevaux observations of the B ring (1981, Table 2). The thermal inertia was constrained to be for a Bond albedo at 3 (Fig. 6). It is close to the qualitative estimate made by Froidevaux of a thermal inertia comparable to that of the Galilean satellites. As can be seen in Fig. 6, the brightness temperature can decrease by more than 15 K in the planetary shadow when the thermal inertia is on the order of 5 to 20 . This is compatible with Pioneer measurements of a 10-15 K drop in the shadow for A and C rings (Froidevaux & Ingersoll 1980).
Two separate fits were run on the CAMIRAS data, for each ring at each epoch: the first one takes data of both ansae into account, the second one only fits data of East ansa. This is to account for the fact that the model fails to reproduce the azimuthal asymmetry on the West ansa of the B ring and not on the East ansa. Azimuthal temperature curves of the best fits for the B ring are shown in Fig. 7. The model fails to reproduce the observed ring temperatures with f=4, which means that the observed temperatures are too high for the particles to be almost isothermal emitters in the frame of our model. For f=2, B ring values for thermal inertia and Bond albedo are summarized in Fig. 9. The Bond albedo values found at both July 1999 and March 2000 epochs are very consistent. This means that the model is capable of treating the seasonal variations in temperature due to variations in the sun elevation and distance. Thermal inertias are found to be very low, about a few , i.e. about three orders of magnitude lower than the thermal inertia of solid water ice. Thermal inertia is significantly lower in March compared to July, but the residuals of the fit in March are very bad. The model does not seem to be able to reproduce the observed slope in temperature on the East ansa (Fig. 7). Given the residuals for the B ring in March, these solutions were rejected. On July 1999, the average Bond albedo of the B ring is and the average thermal inertia is .
Fits to the C ring data are better since there are no bright asymmetries on the West ansa that cannot be reproduced by the model (Fig. 8).
Also for this ring, no solution for the Bond albedo is found if f=4.
This, once again, rejects a ring of quasi-isothermal particles.
For f=2, results are summarized in Fig. 9.
Bond albedos are again very consistent between epochs, and are significantly lower than are those of B ring.
Thermal inertias are also very low in this ring and almost equal to that of the B ring.
They were marginally lower in March 2000 than in July 1999.
The average Bond albedo of the C ring is
and its average thermal inertia is
|Figure 9: Thermal inertias and Bond albedos of B and C rings. Solutions are found by using data from: (triangle) CAMIRAS, July 1999 with East ansa only; (square) CAMIRAS, July 1999 with both ansae; (cross) CAMIRAS, March 2000 with East ansa only; (full circle) CAMIRAS, March 2000 with both ansae; (losange) Froidevaux 1981(see Fig. 6). Errors bars are given at level.|
|Open with DEXTER|
The bolometric Bond albedoes are compatible with the most recent estimates obtained at visual wavelengths by Cuzzi & Estrada (1998) using Voyager images, i.e. 0.46 and 0.18 for B and C rings, respectively. They were obtained for highly anisotropic emitters exclusively, which is new in the case of the C ring. Given the small thermal inertia, particles have to be slowly spinning to emit by one hemisphere only (f=2). Froidevaux (1981) could obtain a fit between his monolayer-model, slowly spinning particles (f=2) and C ring data only for a bolometric Bond albedo of 0.55, which was too high compared to the expected value of 0.25 at epoch (Esposito et al. 1984, Table 6). A lower albedo value was possible for f=4, i.e. in case of very small particles, fast rotation or high thermal inertia. A large spin rate corresponds to a large velocity dispersion in the ring and a vertical thickness that is not compatible with the monolayer. The monolayer model was rejected for the C ring for this reason (Esposito et al. 1984). The C ring brightness temperature derived from CAMIRAS observations is about 7 K above the value expected by the model for . This shift in temperature by itself explains the low values of the Bond albedo obtained here. Monolayer, slowly spinning particles, and low albedo are no longer incompatible for the C ring. The monolayer structure for the C ring should be reconsidered then, especially as observations of waves in this ring have since suggested a vertical thickness of less than 1 m (Rosen 1989), i.e. a monolayer structure for a population whose effective size is about 2 m (French and Nicholson 2000).
Thermal inertias are found to be very low. Many other icy bodies in the solar system also exhibit comparable thermal inertias, yet slightly larger. Ground-based observations of Europa have given thermal inertia as low as 17 for low albedo regions ( , Blaney et al. 1999). The thermal inertia of higher albedo regions of Europa, of the other Galilean satellites or of saturnian satellites like Rhea, Dione or Tethys range between 30 and 70 (Spencer 1989, 1992). The thermal inertias of Chiron or Asbolus Centaurs have been very recently estimated at 10 or less (Fernandez et al. 2002; Groussin et al. 2000).
A low thermal inertia means that material on the surface of a particle cannot store heat or else it insulates the subsurface. For a compacted surface, the volumic density might be close to and its specific heat close to solid water ice value, i.e. or 845 for temperatures in the B and C ring, respectively. The low thermal inertia would then be due to an insulating regolith with low thermal conductivity K. For the average thermal inertias of B and C rings, K would be on the order of . This is five orders of magnitude smaller than the expected thermal conductivity for solid water ice. The standard formula given by Klinger (1981) for crystalline compact water ice is K0=567/T0, i.e. on average for B and C rings at our epochs. With such low thermal conductivities, the thermal skin depth is limited to 0.4 mm, given the mean motions of B and C ring particles. This is small compared to the typical size of Saturn's ring particles. The hypothesis we made to consider particles large compared to the thermal skin depth is reasonable then. Near-infrared spectroscopy suggests that Saturn's rings particles are covered with water ice frost with a typical grain size smaller than (Puetter & Russel 1977; Molina et al. 1992). More detailed modeling of ring spectra by Poulet et al. (2002) tends to confirm the predominence of such grain size even if some fraction of larger grains ( ) has to be present. The thermal skin depth might then be much larger than the typical regolith grain size. In these conditions the measured thermal inertia is the thermal inertia of the regolith and not that of the grain material. This depth is much larger than the electrical skin depth where is the wavelength of observation and nk the imaginary part of the refractive index of crystalline water ice. From laboratory experiments at 100 K (Bertie 1969), . Then . Therefore the observed temperature is actually the surface temperature at 20.5 .
Kouchi et al. (1992) have produced laboratory measurements of the thermal conductivity of amorphous water ice at temperatures around 130 K while slowly depositing a thin film on a cold metal substrate. The mean thermal conductivity of the deposit (mixing in fact amorphous and crystalline water ice) has been found to decrease with decreasing deposition rate, from 10-1 to 10-5 with rate from 10-5 to 10-9 . Most measurements are below 10-4 for deposition rate below 10-6 . This low thermal conductivity might not be due to the amorphous ice structure but simply to defects in the ice deposit (Kouchi et al. 1992; Bunch et al. 1998), like cracks that would prevent heat diffusion while keeping the porosity low. The porosity is indeed limited in their samples as the infrared absorptivity indicates a mean density of about 900 . These experimental results are very interesting, as they give a proper range for thermal conductivity and a plausible origin due to cracks. Cracks might indeed grow easily under the surfaces of particles which are eroded by inter-particle collisions. Fractures may be created at the impact point and ejecta from the collision be compressed on the surface, creating a compacted frosty regolith (Weidenschilling et al. 1984). Thermal stress during nights and days of the spinning particle (possibly few times 10 K-gradient in one orbit) or during the crossing of the planetary shadow may also create fractures on the surface, as would highly energetic particles which break molecules bonds on the surface and fragilize the frosty regolith cover (Borderies et al. 1984).
The porosity p of the regolith, if likely, is unknown. It may directly affect the heat transfer into the medium by reducing the available mass participating to heat conduction. Porosity p lowers the heat capacity by a factor (1-p). The effective thermal conductivity of this porous material is also affected as the area of contacts, through which heat flows, are reduced. In this case K=(1-p)K0, where K0 is the thermal conductivity of the regolith grains. For high porosity like p=0.9, the thermal inertia corresponds to a thermal conductivity . This is still three orders of magnitudes below the thermal conductivity of water ice. So, porosity, even high, cannot alone explain such low thermal inertias, as far as the theory of effective medium is valid for our case.
The thermal inertias are low, and particles have recovered most of their pre-eclipse temperature on the morning east ansa. Other azimuthal temperature variations are visible in the B ring on the West afternoon ansa, which cannot be reproduced by the present model. Also some systematic discrepancy appears between thermal inertias measured in July 1999 and March 2000. These may be due to some phase effect not handle by our model. A possible origin for such azimuthal temperature variations and apparent phase effect is under investigation (Ferrari & Leyrat 2005).
CAMIRAS data show that previously undetected azimuthal temperature asymmetries are present in the West ansa of the B ring, at a place where temperature is expected to be about constant due to the low thermal inertia of particles. Our work also shows that the C ring temperatures can be easily reproduced by a model that assumes it is a monolayer of slowly rotating particles with low albedo and low thermal inertia. Also, if distributed in a monolayer, B ring particles can hardly be quasi-isothermal.
In the next months, the CIRS spectrometer onboard the CASSINI spacecraft will send back unique data on the heating and cooling of A, B, and C rings with very diverse viewing geometries. It will undoubtedly provide new insights on the azimuthal temperature asymmetries and their variations with viewing angles. Direct observations of the shadow will also help to better constrain their thermal inertia and probably their rotational properties.
Many thanks to P. Masse and R. Jouan for their fundamental contribution to the success of this observing campaign. This work was supported by the Programme National de Planétologie and Centre National d'Études Spatiales.