Thermal model of water and CO activity of Comet C/1995 O1 (Hale-Bopp)

An investigation of the activity of Comet C/1995 O1 (Hale-Bopp) with a thermophysical nucleus model that does not rely on the existence of amorphous ice is presented. Our approach incorporates recent observations allowing to constrain important parameters that control cometary activity. The model accounts for heat conduction, heat advection, gas diffusion, sublimation, and condensation in a porous ice-dust matrix with moving boundaries. Erosion due to surface sublimation of water ice leads to a moving boundary. The movement of the boundary is modeled by applying a temperature remapping technique which allows us to account for the loss in the internal energy of the eroded surface material. These kind of problems are commonly referred to as Stefan problems. The model takes into account the diurnal rotation of the nucleus and seasonal effects due to the strong obliquity of Hale-Bopp as reported by Jorda et al. (Jorda, L., Rembor, K., Lecacheux, J., Colom, P., Colas, F., Frappa, E., Lara, L.M. [1997]. Earth Moon Planets 77, 167-180). Only bulk sublimation of water and CO ice are considered without further assumptions such as amorphous ices with certain amount of occluded CO gas. Confined and localized activity patterns are investigated following the reports of Lederer and Campins (Lederer, S.M., Campins, H. [2002]. Earth Moon Planets 90, 381-389) about the chemical heterogeneity of Hale-Bopp and of Bockelée-Morvan et al. (Bockelée-Morvan, D., Henry, F., Biver, N., Boissier, J., Colom, P., Crovisier, J., Despois, D., Moreno, R., Wink, J. [2009]. Astron. Astrophys. 505, 825-843) about a strong CO source at a latitude of 20°. The best fit to the observations of Biver et al. (Biver, N. et al. [2002]. Earth Moon Planets 90, 5-14) is obtained with a low thermal conductivity of 0.01 W m⁻¹ K⁻¹. This is in agreement with recent results of the Deep Impact mission to 9P/Tempel 1 (Groussin, O., A’Hearn, M.F., Li, J.-Y., Thomas, P.C., Sunshine, J.M., Lisse, C.M., Meech, K.J., Farnham, T.L., Feaga, L.M., Delamere, W.A. [2007]. Icarus 187, 16-25) and with previous thermal simulations (Kührt, E. [1999]. Space Sci. Rev. 90, 75-82). The water production curve matches the production rates well from −4 AU pre-perihelion to the outgoing leg while the model does not reproduce so well the water production beyond 4 AU pre-perihelion. The CO production curve is a good fit to the measurements of Biver et al. (2002) over the whole measured heliocentric range from −7 AU pre- to 15 AU post-perihelion.     

The opposition and tilt effects of Saturn’s rings from HST observations

The two major factors contributing to the opposition brightening of Saturn’s rings are (i) the intrinsic brightening of particles due to coherent backscattering and/or shadow hiding on their surfaces, and (ii) the reduced interparticle shadowing when the solar phase angle α → 0°. We utilize the extensive set of Hubble Space Telescope observations (Cuzzi, J.N., French, R.G., Dones, L. [2002]. Icarus 158, 199–223) for different elevation angles B and wavelengths λ to disentangle these contributions. We assume that the intrinsic contribution is independent of B, so that any B dependence of the phase curves is due to interparticle shadowing, which must also act similarly for all λ’s. Our study complements that of Poulet et al. (Poulet, F., Cuzzi, J.N., French, R.G., Dones, L. [2002]. Icarus 158, 224), who used a subset of data for a single B ∼ 10°, and the French et al. (French, R.G., Verbiscer, A., Salo, H., McGhee, C.A., Dones, L. [2007b] PASP 119, 623–642) study for the B ∼ 23° data set that included exact opposition. We construct a grid of dynamical/photometric simulation models, with the method of Salo and Karjalainen (Salo and Karjalainen [2003]. Icarus 164, 428–460), and use these simulations to fit the elevation-dependent part of opposition brightening. Eliminating the modeled interparticle component yields the intrinsic contribution to the opposition effect: for the B and A rings it is almost entirely due to coherent backscattering; for the C ring, an intraparticle shadow hiding contribution may also be present.Based on our simulations, the width of the interparticle shadowing effect is roughly proportional to B. This follows from the observation that as B decreases, the scattering is primarily from the rarefied low filling factor upper ring layers, whereas at larger B’s the dense inner parts are visible. Vertical segregation of particle sizes further enhances this effect. The elevation angle dependence of interparticle shadowing also explains most of the B ring tilt effect (the increase of brightness with elevation). From comparison of the magnitude of the tilt effect at different filters, we show that multiple scattering can account for at most a 10% brightness increase as B → 26°, whereas the remaining 20% brightening is due to a variable degree of interparticle shadowing. The negative tilt effect of the middle A ring is well explained by the the same self-gravity wake models that account for the observed A ring azimuthal brightness asymmetry (Salo, H., Karjalainen, R., French, R.G. [2004]. Icarus 170, 70–90; French, R.G., Salo, H., McGhee, C.A., Dones, L. [2007]. Icarus 189, 493–522).     

New Cassini RADAR results for Saturn’s icy satellites

Cassini radar tracks on Saturn’s icy satellites through the end of the Prime Mission in 2008 have increased the number of radar albedo estimates from 10 (Ostro et al., 2006) to 73. The measurements sample diverse subradar locations (and for Dione, Rhea, and Iapetus almost always use beamwidths less than half the target angular diameters), thereby constraining the satellites’ global radar albedo distributions. The echoes result predominantly from volume scattering, and their strength is thus strongly sensitive to ice purity and regolith maturity. The combination of the Cassini data set and Arecibo 13-cm observations of Enceladus, Tethys, Dione, Rhea (Black et al., 2007), and Iapetus (Black et al., 2004) discloses an unexpectedly complex pattern of 13-to-2-cm wavelength dependence. The 13-cm albedos are generally smaller than 2-cm albedos and lack the correlation seen between 2-cm and optical geometric albedos. Enceladus and Iapetus are the most interesting cases. We infer from hemispheric albedo variations that the E-ring has a prominent effect on the 13-cm radar “lightcurve”. The uppermost trailing-side regolith is too fresh for meteoroid bombardment to have developed larger-scale heterogeneities that would be necessary to elevate the 13-cm radar albedo, whereas all of Enceladus is clean and mature enough for the 2-cm albedo to be uniformly high. For, Iapetus, the 2-cm albedo is strongly correlated with optical albedo: low for the optically dark, leading-side material and high for the optically bright, trailing-side material. However, Iapetus’ 13-cm albedo values show no significant albedo dichotomy and are several times lower than 2-cm values, being indistinguishable from the weighted mean of 13-cm albedos for main-belt asteroids, 0.15 ± 0.10. The leading side’s optically dark contaminant must be present to depths of at least one to several decimeters, so 2-cm albedos can mimic the optical dichotomy; however, it does not have to extend any deeper than that. The fact that both hemispheres of Iapetus look Asteroid-like at 13 cm means that coherent backscattering itself is not nearly as effective as it is at 2 cm. Since Iapetus’ entire surface is mature regolith, the wavelength dependence must involve composition, not structure. Either the composition is a function of depth everywhere (with electrical loss much greater at depths greater than a decimeter or two), or the intrinsic electrical loss of some pervasive constituent is much higher at 13 cm than at 2 cm. Ammonia is a candidate for such a contaminant. If ammonia’s electrical properties do not depend on frequency, and if ammonia is globally much less abundant within the upper one or two decimeters than at greater depths, then coherent backscattering would effectively be shut down at 13 cm, explaining the Asteroid-like 13-cm albedo.     

On Titan’s Xanadu region

A large, circular marking ∼1800 km across is seen in near-infrared images of Titan. The feature is centered at 10°S, 120°W on Titan, encompasses much of Titan’s western Xanadu region, and has an off-center, quasi-circular, inner margin about 700 km across, with lobate outer margins extending 200-500 km from the inner margin. On the feature’s southern flank is Tui Regio, an area that has very high reflectivity at 5 μm, and is hypothesized to exhibit geologically recent cryovolcanic flows (Barnes, J.W. et al. [2006]. Geophys. Res. Lett. 33), similar to flows seen in Hotei Regio, a cryovolcanic area whose morphology may be controlled by pre-existing, crustal fractures resulting from an ancient impact (Soderblom, L.A. et al. [2009]. Icarus, 204). The spectral reflectivity of the large, circular feature is quite different than that of its surroundings, making it compositionally distinct, and radar measurements of its topography, brightness temperature and volume scattering also suggest that the feature is quite distinct from its surroundings. These and several other lines of evidence, in addition to the feature’s morphology, suggest that it may occupy the site of an ancient impact.     

Properties of the Hermean regolith: IV. Photometric parameters of Mercury and the Moon contrasted with Hapke modelling

A comparison of the photometric properties of Mercury and the Moon is performed, based on their integral phase curves and disk-resolved image data of Mercury obtained with the Swedish Vacuum Solar Telescope. Proper absolute calibration of integral V-band magnitude observations reveals that the near-side of the Moon is 10-15% brighter than average Mercury, and 0-5% brighter for the “bolometric” wavelength range 400-1000 nm. As shown, this is supported by recent estimates of their geometric albedos. Hapke photometric parameters of their surfaces are derived from identical approaches, allowing a contrasting study between their surface properties to be performed. Compared to the average near-side Moon, Mercury has a slightly lower single-scattering albedo, an opposition surge with smaller width and of marginally smaller amplitude, and a somewhat smoother surface with similar porosity. The width of the lobes of the single-particle scattering function are smaller for Mercury, and the backward scattering anisotropy is stronger. In terms of the double Henyey-Greenstein b-c parameter plot, the scattering properties of an average particle on Mercury is closer to the properties of lunar maria than highlands, indicating a higher density of internal scatterers than that of lunar particles. The photometric roughness of Mercury is well constrained by the recent study of Mallama et al. (2002, Icarus 155, 253-264) to a value of about 8°, suggesting that the surfaces sampled by the highest phase angle observations (Borealis, Susei, and Sobkou Planitia) are lunar mare-like in their textural properties. However, Mariner 10 disk brightness profiles obtained at intermediate phase angles indicate a surface roughness of about twice this value. The photometric parameters of the Moon are more difficult to constrain due to limited phase angle coverage, but the best Hapke fits are provided by rather small surface roughnesses. Better-calibrated, multiple-wavelength observations of the integral and disk-resolved brightnesses of both bodies, and obtained at higher phase angle values in the case of the Moon, are urgently needed to arrive at a more consistent picture of the contrasting light scattering properties of their surfaces.     

Laboratory Study of the Bidirectional Reflectance of Powdered Surfaces: On the Asymmetry Parameter of Asteroid Photometric Data

We performed laboratory experiments on the bidirectional reflectance of powdered surfaces of dunite, graphite, and Allende (CV3) and Gao (H5) meteorites. The particle size of each sample varied from tens to hundreds of micrometers. The absolute bidirectional reflectance was determined for two ranges of phase angle: from 2° to 80° at 0° angle of incidence and from 2° to 155° at 75° angle of incidence. The phase angle was incremented every 1° in between 2° and 5° and every 10° between 10° and 150°. We determined the Hapke parameters and found that the values of the asymmetry parameter retrieved for most of the samples were positive when the coverage of the phase angle was wide, 2° to 155°, although the values derived from remote sensing instruments on asteroid flyby missions have been negative or nearly zero. Among our sample surfaces, only those of graphite, Allende, and Gao with particle sizes of hundreds micrometers show negative or nearly zero values. The single-scattering albedos determined for the Gao samples are comparable to the values of S-class asteroids, while those for the Allende samples are much larger than those of C-class asteroids.     

A new look at photometry of the Moon

We use ROLO photometry (Kieffer, H.H., Stone, T.C. [2005]. Astron. J. 129, 2887-2901) to characterize the before and after full Moon radiance variation for a typical highlands site and a typical mare site. Focusing on the phase angle range 45° < α < 50°, we test two different physical models, macroscopic roughness and multiple scattering between regolith particles, for their ability to quantitatively reproduce the measured radiance difference. Our method for estimating the rms slope angle is unique and model-independent in the sense that the measured radiance factor I/F at small incidence angles (high Sun) is used as an estimate of I/F for zero roughness regolith. The roughness is determined from the change in I/F at larger incidence angles. We determine the roughness for 23 wavelengths from 350 to 939 nm. There is no significant wavelength dependence. The average rms slope angle is 22.2° ± 1.3° for the mare site and 34.1° ± 2.6° for the highland site. These large slopes, which are similar to previous “photometric roughness” estimates, require that sub-mm scale “micro-topography” dominates roughness measurements based on photometry, consistent with the conclusions of Helfenstein and Shepard (Helfenstein, P., Shepard, M.K. [1999]. Icarus 141, 107-131). We then tested an alternative and very different model for the before and after full Moon I/F variation: multiple scattering within a flat layer of realistic regolith particles. This model consists of a log normal size distribution of spheres that match the measured distribution of particles in a typical mature lunar soil 72141,1 (McKay, D.S., Fruland, R.M., Heiken, G.H. [1974]. Proc. Lunar Sci. Conf. 5, Geochim. Cosmochim. Acta 1 (5), 887-906). The model particles have a complex index of refraction 1.65-0.003i, where 1.65 is typical of impact-generated lunar glasses. Of the four model parameters, three were fixed at values determined from Apollo lunar soils: the mean radius and width of the log normal size distribution and the real part of the refraction index. We used FORTRAN programs from Mishchenko et al. (Mishchenko, M.I., Dlugach, J.M., Yanovitskij, E.G., Zakharova, N.T. [1999]. J. Quant. Spectrosc. Radiat. Trans. 63, 409-432; Mishchenko, M.I., Travis, L.D., Lacis, A.A. [2002]. Scattering, Absorption and Emission of Light by Small Particles. Cambridge Univ. Press, New York. <http://www.giss.nasa.gov/staff/mmishchenko/books.html>) to calculate the scattering matrix and solve the radiative transfer equation for I/F. The mean single scattering albedo is ω = 0.808, the asymmetry parameter is 〈cos Θ〉 = 0.77 and the phase function is very strongly peaked in both the forward and backward scattering directions. The fit to the observations for the highland site is excellent and multiply scattered photons contribute ?80% of I/F. We conclude that either model, roughness or multiple scattering, can match the observations, but that the strongly anisotropic phase functions of realistic particles require rigorous calculation of many orders of scattering or spurious photometric roughness estimates are guaranteed. Our multiple scattering calculation is the first to combine: (1) a regolith model matched to the measured particle size distribution and index of refraction of the lunar soil, (2) a rigorous calculation of the particle phase function and solution of the radiative transfer equation, and (3) application to lunar photometry with absolute radiance calibration.     

Saturn’s wayward shepherds: the peregrinations of Prometheus and Pandora

Saturn’s narrow F ring is flanked by two nearby small satellites, Prometheus and Pandora, discovered in Voyager images taken in 1980 and 1981 (Synnott et al., 1983, Icarus 53, 156-158). Observations with the Hubble Space Telescope (HST) during the ring plane crossings (RPX) of 1995 led to the unexpected finding that Prometheus was ∼19° behind its predicted orbital longitude, based on the Synnott et al. (1983) Voyager ephemeris (Bosh and Rivkin, 1996 Science 272, 518-521; Nicholson et al., 1996, Science 272, 509-515). Whereas Pandora was at its predicted location in August 1995, McGhee (2000, Ph.D. thesis, Cornell University) found from the May and November 1995 RPX data that Pandora also deviates from the Synnott et al. (1983) Voyager ephemeris. Using archival HST data from 1994, previously unexamined RPX images, and a large series of targeted WFPC2 observations between 1996 and 2002, we have determined highly accurate sky-plane positions for Prometheus, Pandora, and nine other satellites found in our images. We compare the Prometheus and Pandora measurements to the predictions of substantially revised and improved ephemerides for the two satellites based on an extensive analysis of a large set of Voyager images (Murray et al., 2000, Bull. Am. Astron. Soc. 32, 1090; Evans, 2001 Ph.D. thesis, Queen Mary College). From December 1994 to December 2000, Prometheus’ orbital longitude lag was changing by −0.71° year⁻¹ relative to the new Voyager ephemeris. In contrast, Pandora is ahead of the revised Voyager prediction. From 1994 to 2000, its longitude offset changed by +0.44° year⁻¹, showing in addition an ∼585 day oscillatory component with amplitude Δλ_CR0 = 0.65 ± 0.07° whose phase matches the expected perturbation due to the nearby 3:2 corotation resonance with Mimas, modulated by the 71-year libration in the longitude of Mimas due to its 4:2 resonance with Tethys. We determine orbital elements for freely precessing equatorial orbits from fits to the 1994-2000 HST observations, from which we conclude that Prometheus’ semimajor axis was 0.31 km larger, and Pandora’s was 0.20 km smaller, than during the Voyager epoch. Subsequent observations in 2001-2002 reveal a new twist in the meanderings of these satellites: Prometheus’ mean motion changed suddenly by an additional −0.77° year⁻¹, equivalent to a further increase in semimajor axis of 0.33 km, at the same time that Pandora’s mean motion changed by +0.92° year⁻¹, corresponding to a change of −0.42 km in its semimajor axis. There is an apparent anticorrelation of the motions of these two moons seen in the 2001-2002 observations, as well as over the 20-year interval since the Voyager epoch. This suggests a common origin for their wanderings, perhaps through direct exchange of energy between the satellites as the result of resonances, possibly involving the F ring.     

Photometric analysis of Eros from NEAR data

A photometric model of (433) Eros at wavelengths from 450 to 1050 nm is constructed using the combination of the images from the multispectral imager (MSI) obtained during the one-year long orbital phase of the NEAR mission, ground-based lightcurves from earlier observations, and our theoretical forward modeling simulations coupled with the NEAR shape model. The single scattering albedo is found to be 0.33±0.03 at 550 nm, which is smaller than past findings by 30%. The amplitude and width of the opposition effect are 1.4±0.1 and 0.010±0.004 from ground based lightcurves. It is confirmed that the asymmetry factor of the single-particle phase function and the surface roughness parameter do not depend on wavelength from 450 to 1050 nm, and their values are estimated to be −0.25±0.02 and 28°±3°, respectively, comparable with the earlier measurements from the NEAR NIS data. The geometric albedo and the Bond albedo at 550 nm are calculated to be 0.23 and 0.093, respectively, which make Eros less reflective than previous models, but still slightly more reflective than average S-type asteroids. The lower albedos of Eros are more consistent with our forward modeling simulations, as well as with its spectrum. Eros is a typical S-type asteroid like (951) Gaspra and (243) Ida, and has similar surface regolith properties. Combining the single-scattering albedo with the olivine composition of ordinary chondrites, taking into account space weathering darkening, we constrain the grain size of the regolith particles on Eros to a range of 50 to 100 μm.     

Compositional heterogeneity of Asteroid 4 Vesta’s southern hemisphere: Implications for the Dawn mission

High signal-to-noise, rotationally-resolved spectra of Asteroid 4 Vesta’s southern hemisphere from the 2007 opposition were used to constrain its compositional and mineralogical variations. The spectra were rotationally-phased using closely timed HST observations of Vesta by Li et al. (Li, J.-Y., McFadden, L.A., Thomas, P.C., Mutchler, M.J., Parker, J.Wm., Young, E.F., Russell, C.T., Sykes, M.V., Schmidt, B.E. [2010]. Icarus 208, 238–251). The average surface of Vesta’s southern hemisphere is analogous to a howardite or polymict eucrite assemblage similar to the northern hemisphere, although the band parameters are distinctly shifted towards the diogenite zone on the Band–Band plot. A few distinct compositional units were detected and they might be related to albedo features detected by Hubble Space Telescope (Li et al., 2010). We have identified two compositionally distinct regions overlaying the background surface. The first unit is a polymict eucrite and/or low-Ca eucrite compositional unit at 143° longitude that border the eucrite zone on the Band–Band plot and the second is a diogenite unit at 159°. While we did not detect any distinct olivine units as suggested by Gaffey (Gaffey, M.J. [1997]. Icarus 127, 130–157), we cannot rule out the possibility of smaller olivine-rich units that are below the detection limit of the instrumentation we used. Based on the analysis and the limitations of the data, we do not suggest that Vesta’s surface is olivine-free. Mean pyroxene chemistry estimates for both hemispheres broadly agree with one another (to within one-sigma) with the northern hemisphere ferrosilite (Fs) and wollastonite (Wo) values being slightly higher than southern hemisphere.     

The jovian rings: new results derived from Cassini, Galileo, Voyager, and Earth-based observations

Cassini's Imaging Science Subsystem (ISS) instrument took nearly 1200 images of the Jupiter ring system during the spacecraft's 6-month encounter with Jupiter (Porco et al., 2003, Science 299, 1541-1547). These observations constitute the most complete data set of the ring taken by a single instrument, both in phase angle (0.5°-120° at seven angles) and wavelength (0.45-0.93 μm through eight filters). The main ring was detected in all targeted exposures; the halo and gossamer rings were too faint to be detected above the planet's stray light. The optical depth and radial profile of the main ring are consistent with previous observations. No broad asymmetries within the ring were seen; we did identify possible hints of 1000 km-scale azimuthal clumps within the ring. Cassini observations taken within 0.02° of the ring plane place an upper limit on the ring's full thickness of 80 km at a phase angle of 64°. We have combined the Cassini ISS and VIMS (Visible and Infrared Mapping Spectrometer) observations with those from Voyager, HST (Hubble Space Telescope), Keck, Galileo, Palomar, and IRTF (Infrared Telescope Facility). We have fit the entire suite of data using a photometric model that includes microscopic silicate dust grains as well as larger, long-lived ‘parent bodies’ that engender this dust. Our best-fit model to all the data indicates an optical depth of small particles of τ_s=4.7×10⁻⁶ and large bodies τ_l=1.3×10⁻⁶. The dust's cross-sectional area peaks near 15 μm. The data are fit significantly better using non-spherical rather than spherical dust grains. The parent bodies themselves must be very red from 0.4-2.5 μm, and may have absorption features near 0.8 and 2.2 μm.     

Photometric behavior of 20000 Varuna at very small solar phase angles

Six nights of R-band CCD observations of the classical Kuiper Belt Object (KBO) 20000 Varuna (2000 WR106) were obtained at the Palomar Mountain 60- and 200-in telescopes. The observations were scheduled to take advantage of a particularly favorable apparition which allowed us to sample down to extremely small solar phase angle (α=0.036°). After rotational lightcurve subtraction, we found that the KBO exhibited a strong opposition surge of ∼0.1 mag at phase angles α<0.1°. We modeled our composite solar phase curve of Varuna using both H-G parameterization and Hapke theory and concluded that similar opposition surges may be wide spread among KBOs and that the regolith of Varuna may be significantly more porous than a typical main-belt C-type asteroid. Wide-spread opposition surges lead to higher albedos than derived assuming linear phase behavior: on the whole KBOs may be brighter than previously assumed.     

Radar Observations of Asteroid 3908 Nyx

We report Doppler-only (cw) radar observations of basaltic near-Earth asteroid 3908 Nyx obtained at Arecibo and Goldstone in September and October of 1988. The circular polarization ratio of 0.75±0.03 exceeds ∼90% of those reported among radar-detected near-Earth asteroids and it implies an extremely rough near-surface at centimeter-to-decimeter spatial scales. Echo power spectra over narrow longitudinal intervals show a central dip indicative of at least one significant concavity. Inversion of cw spectra yields two statistically indistinguishable shape models that have similar shapes and dimensions but pole directions that differ by ∼100°. We adopt one as our working model and explore its implications. It has an effective diameter of 1.0±0.15 km and radar and visual geometric albedos of 0.15±0.075 and 0.16^+0.08_−0.05. The visual albedo supports the interpretation by D. P. Cruikshank et al. (1991, Icarus89, 1-13) that Nyx has a thermal inertia consistent with that of bare rock. The model is irregular, modestly asymmetric, and topographically rugged.     

The size distribution of Jupiter's main ring from Galileo imaging and spectroscopy

Galileo's Solid State Imaging experiment (SSI) obtained 36 visible wavelength images of Jupiter's ring system during the nominal mission (Ockert-Bell et al., 1999, Icarus 138, 188-213) and another 21 during the extended mission. The Near Infrared Mapping Spectrometer (NIMS) recorded an observation of Jupiter's main ring during orbit C3 at wavelengths from 0.7 to 5.2 μm; a second observation was attempted during orbit E4. We analyze the high phase angle NIMS and SSI observations to constrain the size distribution of the main ring's micron-sized dust population. This portion of the population is best constrained at high phase angles, as the light scattering behavior of small dust grains dominates at these geometries and contributions from larger ring particles are negligible. High phase angle images of the main ring obtained by the Voyager spacecraft covered phase angles between 173.8° and 176.9° (Showalter et al., 1987, Icarus 69, 458-498). Galileo images extend this range up to 178.6°. We model the Galileo phase curve and the ring spectra from the C3 NIMS ring observation as the combination of two power law distributions. Our analysis of the main ring phase curve and the NIMS spectra suggests the size distribution of the smallest ring particles is a power law with an index of 2.0±0.3 below a size of ∼15 μm that transitions to a power law with an index of 5.0±1.5 at larger sizes. This combined power law distribution, or “broken power law” distribution, yields a better fit to the NIMS data than do the power law distributions that have previously been fit to the Voyager imaging data (Showalter et al., 1987, Icarus 69, 458-498). The broken power law distribution reconciles the results of Showalter et al. (1987, Icarus 69, 458-498) and McMuldroch et al. (2000, Icarus 146, 1-11), who also analyzed the NIMS data, and can be considered as an obvious extension of a simple power law. This more complex size distribution could indicate that ring particle production rates and/or lifetimes vary with size and may relate to the physical processes that control their evolution. The significant near arm/far arm asymmetry reported elsewhere (see Showalter et al., 1987, Icarus 69, 458-498; Ockert-Bell et al., 1999, Icarus 138, 188-213) persists in the data even after the main ring is isolated in the SSI images. However, the sense of the asymmetry seen in Galileo images differs from that seen in Voyager images. We interpret this asymmetry as a broad-scale, azimuthal brightness variation. No consistent association with the magnetic field of Jupiter has been observed. It is possible that these longitudinal variations may be similar to the random brightness fluctuations observed in Saturn's F ring by Voyager (Smith et al., 1982, Science 215, 504-537) and during the 1995 ring plane crossings (Nicholson et al., 1996, Science 272, 509-515; Bosh and Rivkin, 1996, Science 272, 518-521; Poulet et al., 2000, Icarus 144, 135-148). Stochastic events may thus play a significant role in the evolution of the jovian main ring.     

Search for Steins’ surface inhomogeneities from OSIRIS Rosetta images

We investigate the possible presence of heterogeneous surface features on asteroid (2867) Steins, using the G-mode multivariate statistical method (Coradini et al., 1977) applied to Rosetta/OSIRIS images. We analyze both NAC and WAC images obtained near around the closest approach that occurred on September 5th, 2008, through different filters centered on wavelengths ranging from 295 to 986 nm. The shape of Steins is modeled as a polyhedron of almost 58 000 facets. Photometric corrections were performed using Hapke's (2002) model to compensate for the variable illuminations conditions at the surface. The G-mode classification method was performed on all visible and illuminated facets, i.e. in a region limited to [−50°,+60°] in latitude and [−40°,+90°] in longitude, that represents almost 30% of the total surface. The analyzed set of facets does not show any significant difference in the reflected light content, suggesting no surface inhomogeneities larger than 4% at the 95% confidence level.     

Neptune’s cloud and haze variations 1994-2008 from 500 HST-WFPC2 images

The analysis of all suitable images taken of Neptune with the Wide Field Planetary Camera 2 on the Hubble Space Telescope between 1994 and 2008 revealed the following results. The activity of discrete cloud features located near Neptune’s tropopause remained roughly constant within each year but changed significantly on the time scale of ∼5 years. Discrete clouds covered 1% of the disk on average, but more than 2% in 2002. The other ∼99% of the disk probed Neptune’s hazes at lower altitudes. At red and near-infrared wavelengths, two dark bands around −70° and 10° latitude were perfectly steady and originated in the upper two scale heights of the troposphere, either by decreased haze opacity or by an increased methane relative humidity. At blue wavelengths, a dark band between −60° and −30° latitude was most obvious during the early years, caused by dark aerosols below the 3-bar level with single scattering albedos reduced by ∼0.04, and this contrast was constant between 410 and 630 nm wavelength. The dark band decayed exponentially with a time constant of 5 ± 1 years, which can be explained by settling of the dark aerosols at a rate of 1 bar pressure difference per year. The other latitudes brightened with the same time constant but lower amplitudes. The only exception was a darkening event in the 15-30° latitude region between 1994 and 1996, which coincides with two dark spots observed in the same region during the same time period, the only dark spots seen since Voyager. The dark aerosols had a similar latitudinal distribution as the discrete clouds near the tropopause, although both were separated by four scale heights. Photometric analysis revealed a phase coefficient of 0.0028 ± 0.0010 mag/deg for the 0-2° phase-angle range observable from Earth. Neptune’s sub-Earth latitude varied by less than 3° throughout the observation period providing a data set with almost constant viewing geometry. The trends observed up to 2008 continued into 2010 based on images taken with the Wide Field Camera 3.     

About mineral composition of geologic units in the northern hemisphere of Vesta

In the present paper we seek to understand the geologic diversity of units in the northern hemisphere of Vesta using HST observations (Binzel et al., 1997). First, we compare colors R(0.673 μm)/R(0.953 μm) and R(0.673 μm)/R(1.042 μm) of Vesta’s units with those of V-type asteroids (vestoids) as well as howardite, eucrite, and diogenite meteorites (HEDs). This comparative analysis showed that: (i) on the color-color plot, regions on Vesta are clustered whereas vestoids and HEDs cover a wide range in color; (ii) very few vestoids or HEDs fall into Vesta’s color region. This implies that Vesta’s units are more homogenous than most vestoids and HEDs examined here and material of the units are slightly different from that of vestoids and HEDs. Assuming reasonable choice of end-member materials, an optical model (Shkuratov et al., 1999) was used to simulate intimate mixtures of particles at the surface of Vesta’s units. Simulation of albedo, colors, and four-point spectra of Vesta’s units reveals that the rock-forming material is nearly equal for all units and has HED-like composition. Diversity of the units depends on the minor constituents such as chromite and a neutral phase. The western units contain more chromite and neutral phase than the eastern, consequently albedo of the western units is lower and their four-point spectra are flatter. Olivine and feldspar are also needed to give the best fit for the calculated and observed albedos and colors of Vesta’s units, but being in minor amount in Vesta’s rocks they play a secondary role in contributing to the optical properties of the units. Questions about the proportions of HED-like rock and the constituent called neutral phase remain open. Spectrophotometric studies of Vesta with both higher spatial and spectral resolution as expected from NASA’s Dawn mission are needed for resolving these problems.     

On Hapke photometric model predictions on reflectance of closely packed particulate surfaces

In a recent paper Hapke et al. (Hapke, B., Shepard, M., Nelson, R., Smythe, W., Piatek, J. [2009]. Icarus 199, 210-218) performed bi-directional reflectance measurements on closely-packed particulate surfaces of micrometer-sized particles and compared these with both the Hapke IMSA photometric model, and a numerical radiative transfer algorithm, the MDYZ (Mishchenko, M., Dlugach, J., Yanovitskij, E., Zakharova, N. [1999a]. J. Quant. Spectrosc. Radiat. Trans. 63, 409-432). To account for the effects of close packing, Hapke et al. applied a diffraction truncation scheme to remove the diffraction spike and supplied the renormalized single scattering phase function to the IMSA. They found that the IMSA prediction is a better match with measurement than that of MDYZ. In this work we demonstrate that the diffraction truncation procedure outlined by Hapke et al. contains an error. By following Hapke et al.’s intended truncation scheme, we have found that the IMSA model is not sufficiently anisotropic to describe the reflectance pattern of measurements on surface reflectance of closely packed large spherical particles.     

The first detection of propane on Saturn

We report the first detection of propane, C₃H₈, in Saturn's stratosphere. Observations taken on September 8, 2002 UT at NASA's IRTF using TEXES, show multiple emission lines due to the 748 cm⁻¹ν₂₁ band of C₃H₈. Using a line-by-line radiative transfer code, we are able to fit the data by scaling the propane vertical mixing ratio profile from the photochemical model of Moses et al. [2000. Icarus 143, 244-298]. Multiplicative factors of 0.7 and 0.65 are required to fit the −20° and −80° planetocentric latitude spectra. The resultant profiles are characterized by a 5 mbar mixing ratio of 2.7±0.8×¹⁰⁻⁸ at −20° and at −80° latitude. These results suggest that the time scale for meridional circulation lies between the net photochemical lifetimes of C₂H₂ and C₃H₈, ≈30-600 years.