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.     

Carbon Monoxide on Jupiter: Evidence for Both Internal and External Sources

Thirteen lines of the CO band near 4.7 μm have been observed on a jovian hot spot at a resolution of 0.045 cm⁻¹. The measured line profiles indicate that the CO mole fraction is 1.0±0.2 ppb around the 6-bar level and is larger in the upper troposphere and/or stratosphere. An external source of CO providing an abundance of 4⁺³₋₂×10¹⁶ molecules cm⁻² is implied by the observations in addition to the amount deposited at high altitude by the Shoemaker-Levy 9 collision. From a simple diffusion model, we estimate that the CO production rate is (1.5-10)×10⁶ molecules cm⁻² s⁻¹ assuming an eddy diffusion coefficient around the tropopause between 300 and 1500 cm² s⁻¹. Precipitation of oxygen atoms from the jovian magnetosphere or photochemistry of water vapor from meteoroidal material can only provide a negligible contribution to this amount. A significant fraction of the CO in Jupiter's upper atmosphere may be formed by shock chemistry due to the infall of kilometer- to subkilometer-size Jupiter family comets. Using the impact rate from Levison et al. (2000, Icarus143, 415-420) rescaled by Bottke et al. (2002, Icarus156, 399-433), this source can provide the observed stratospheric CO only if the eddy diffusion coefficient around the tropopause is 100-300 cm² s⁻¹. Higher values, ∼700 cm² s⁻¹, would require an impact rate larger by a factor of 5-10, which cannot be excluded considering uncertainties in the distribution of Jupiter family comets. Such a large rate is indeed consistent with the observed cratering record of the Galilean satellites (Zahnle et al. 1998, Icarus136, 202-222). On the other hand, the ∼1 ppb concentration in the lower troposphere requires an internal source. Revisiting the disequilibrium chemistry of CO in Jupiter, we conclude that rapid vertical mixing can provide the required amount of CO at ∼6 bar for a global oxygen abundance of 0.2-9 times the solar value considering the uncertainties in the convective mixing rate and in the chemical constants.     

Keck Adaptive Optics Images of Uranus and Its Rings

We present adaptive optic images of Uranus obtained with the 10-m W. M. Keck II telescope in June 2000, at wavelengths between 1 and 2.4 μm. The angular resolution of the images is ∼0.06-0.09″. We identified eight small cloud features on Uranus's disk, four of which were in the northern hemisphere. The latter features are ∼1000-2000 km in extent and located in the upper troposphere, above the methane cloud, at pressures between 0.5 and 1 bar. Our data have been combined with HST data by Hammel et al. (2001, Icarus153, 229-235); the combination of Keck and HST data allowed derivation of an accurate wind velocity profile. Our images further show Uranus's entire ring system: the asymmetric ? ring, as well as the three groups of inner rings (outward from Uranus): the rings 6+5+4, α+β, and the η+γ+δ rings. We derived the equivalent I/F width and ring particle reflectivity for each group of rings. Typical particle albedos are ∼0.04-0.05, in good agreement with HST data at 0.9 μm.     

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.     

A Prominent Apparition of Neptune's South Polar Feature

During the last week of June 2001, a bright apparition of Neptune's South Polar Feature (SPF) at 70°S was observed to develop and decay in less than 30 hours, displaying contrast of ∼2.5 at 619 nm. Assuming that the same SPF was observed on two consecutive rotations of Neptune, the feature moved eastward at 3.2±1.8° hr⁻¹ (130±80 m s⁻¹). The SPF made no obvious appearances during eight other Hubble Space Telescope (HST) observations of Neptune between July 2000 and June 2001, although there was a faint feature at 70°S in one image in October 2000. A prominent SPF was present in near-IR Keck Telescope images in August 2000. Bright SPFs are seen on ∼10% of the HST images of Neptune obtained since 1994, and a fainter SPF is visible on another ∼10%. An SPF bright enough to be visible at HST resolution was present around half the time during the last week of Voyager's approach to Neptune in August 1989, with one prominent brightening, suggesting that the SPF is less visible now than in 1989.     

No Hexagonal Wave around Saturn's Southern Pole

Hubble Space Telescope (HST) images (1997-2002) do not show a hexagonal wave in the south pole that is a counterpart of the northern one (D. A. Godfrey 1988, Icarus76, 335-356). However, a polar jet similar to that in which the north polar hexagonal wave is embedded has been discovered in the southern hemisphere. The HST images also reveal the existence of a small polar “cap” about 2000 km in diameter that marks the rotational pole of the planet, as well as unexpected changes in the nearby cloud structure on a timescale of years.     

Changes in Jupiter’s zonal velocity between 1979 and 2008

We show that the peak velocity of Jupiter’s visible-cloud-level zonal winds near 24°N (planetographic) increased from 2000 to 2008. This increase was the only change in the zonal velocity from 2000 to 2008 for latitudes between ±70° that was statistically significant and not obviously associated with visible weather. We present the first automated retrieval of fast (∼130 m s⁻¹) zonal velocities at 8°N planetographic latitude, and show that some previous retrievals incorrectly found slower zonal winds because the eastward drift of the dark projections (associated with 5-μm hot spots) “fooled” the retrieval algorithms.We determined the zonal velocity in 2000 from Cassini images from NASA’s Planetary Data System using a global method similar to previous longitude-shifting correlation methods used by others, and a new local method based on the longitudinal average of the two-dimensional velocity field. We obtained global velocities from images acquired in May 2008 with the Wide Field Planetary Camera 2 (WFPC2) on the Hubble Space Telescope (HST). Longer-term variability of the zonal winds is based on comparisons with published velocities based on 1979 Voyager 2 and 1995-1998 HST images. Fluctuations in the zonal wind speeds on the order of 10 m s⁻¹ on timescales ranging from weeks to months were found in the 1979 Voyager 2 and the 1995-1998 HST velocities. In data separated by 10 h, we find that the east-west velocity uncertainty due to longitudinal fluctuations are nearly 10 m s⁻¹, so velocity fluctuations of 10 m s⁻¹ may occur on timescales that are even smaller than 10 h. Fluctuations across such a wide range of timescales limit the accuracy of zonal wind measurements. The concept of an average zonal velocity may be ill-posed, and defining a “temporal mean” zonal velocity as the average of several zonal velocity fields spanning months or years may not be physically meaningful.At 8°N, we use our global method to find peak zonal velocities of ∼110 m s⁻¹ in 2000 and ∼130 m s⁻¹ in 2008. Zonal velocities from 2000 Cassini data produced by our local and global methods agree everywhere, except in the vicinity of 8°N. There, the local algorithm shows that the east-west velocity has large variations in longitude; vast regions exceed ∼140 m s⁻¹. Our global algorithm, and all of the velocity-extraction algorithms used in previously-published studies, found the east-west drift velocities of the visible dark projections, rather than the true zonal velocity at the visible-cloud level. Therefore, the apparent increase in zonal winds between 2000 and 2008 at 8°N is not a true change in zonal velocity.At 7.3°N, the Galileo probe found zonal velocities of 170 m s⁻¹ at the 3-bar level. If the true zonal velocity at the visible-cloud level at this latitude is ∼140 m s⁻¹ rather than ∼105 m s⁻¹, then the vertical zonal wind shear is much less than the currently accepted value.     

Jupiter's Cloud Structure as Constrained by Galileo Probe and HST Observations

The Galileo Probe sampled Jupiter's atmosphere at the edge of a 5-μm hot spot, where it found very little cloud opacity above the 700 mb level. Only τ=1-2 at λ=0.5 μm was inferred from Net Flux Radiometer observations (Sromovsky et al. 1998, J. Geophys. Res.103, 22,929-22,977), in seeming conflict with Chanover et al. (1997, Icarus128, 294-305) who inferred τ=6-8 above the 700 mb level (at λ∼0.9 μm) from 893-nm and 953-nm WFPC2 observations of a group of hot spots. Postulating a heterogeneous cloud structure is one way to resolve the conflict. We obtained a more satisfying resolution by reinterpretation of the HST observations with Probe-compatible assumptions about the vertical distribution of cloud particles. Assuming a physically thin upper (putative NH₃) cloud with adjustable optical depth and effective pressure (p_eff<440 mb) and a physically thin midlevel (putative NH₄SH) cloud with adjustable optical depth but a fixed pressure of 1.2 bars, we are able to fit WPFC2 observations with probe-consistent opacities in hot spot regions. With the same cloud pressures, but higher middle cloud opacities, we are even able to fit the visibly bright regions. Little variability is seen in the upper cloud. Best fits to October 1995 WFPC2 observations in dark regions (5-μm hot spots) yielded τ_upper=1.3-1.9 at 0.9 μm and p_eff=240 mb−270 mb, while in visibly bright regions between hot spots we obtained τ_upper=1.6-2.2 and p_eff=250 mb−290 mb. May 1996 observations yielded slightly higher values of τ_upper (1.8-2.3 and 2.0-2.8) and p_eff (250 mb−310 mb and 265 mb−320 mb). We found that the most important variable parameter is the opacity of the middle cloud, which ra nged from τ=1, 2 in dark regions, to τ=8-30 in bright regions. From limb darkening characteristics, we inferred a wavelength-dependent haze opacity ranging from 0.2±0.05 at 660 nm to 0.35±0.05 at 953 nm, and an effective haze pressure near 120 mb. We did not find it necessary to use low single scattering albedos that require effective imaginary indices, that are several orders of magnitude larger than the values of the main putative cloud components.     

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

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

An Extended CO Source around Comet 29P/Schwassmann-Wachmann 1

Comet 29P/Schwassmann-Wachmann 1 has been studied during seven days in August 1998 with the SEST submillimeter telescope at ESO, La Silla, Chile. The CO (J=2−1) emission at 230 GHz was mapped by directing the telescope beam at the nucleus and six off-nucleus positions. The CO line profiles exhibit the blue- and redshifted components previously observed by various observers. The strength of the observed lines does not decrease with projected distance to the nucleus as expected if CO molecules were coming from the nucleus only. Instead, the line area is nearly constant throughout the map. This can be explained if CO molecules are being released from both the sunlit side of the nucleus and CO-bearing particles distributed in a shell-like cloud. The extended source must consist of icy grains globally moving toward the Sun at ∼50 m s⁻¹ released ∼30 days before the observations were made. The nuclear and extended sources produce (7±1)×10²⁷ and 2.4×10²⁸ molecules s⁻¹, respectively. Our 1996 observations of the comet (Festou, M., M. Gunnarsson, A. Winnberg, H. Rickman, and G. Tancredi 2001. Icarus150, 140-150) were reexamined using this new two-source model. In this case, the nuclear and extended CO sources produced 10±1×10²⁷ and 2.9×10²⁸ CO molecules s⁻¹, respectively. It is not necessary to postulate night side outgassing, but a large quantity of solid grains has to be expelled into the coma.     

Dynamics of cloud features on Uranus

Near-infrared adaptive optics imaging of Uranus by the Keck 2 telescope during 2003 and 2004 has revealed numerous discrete cloud features, 70 of which were used to extend the zonal wind profile of Uranus up to 60° N. We confirmed the presence of a north-south asymmetry in the circulation [Karkoschka, E., 1998. Science 280, 570-572], and improved its characterization. We found no clear indication of long term change in wind speed between 1986 and 2004, although results of Hammel et al. [Hammel, H.B., Rages, K., Lockwood, G.W., Karkoschka, E., de Pater, I., 2001. Icarus 153, 229-235] based on 2001 HST and Keck observations average ∼10 m/s less westward than earlier and later results, and 2003 observations by Hammel et al. [Hammel, H.B., de Pater, I., Gibbard, S., Lockwood, G.W., Rages, K., 2005. Icarus 175, 534-545] show increased wind speeds near 45° N, which we do not see in our 2003-2004 observations. We observed a wide range of lifetimes for discrete cloud features: some features evolve within ∼1 h, many have persisted at least one month, and one feature near 34° S (termed S34) seems to have persisted for nearly two decades, a conclusion derived with the help of Voyager 2 and HST observations. S34 oscillates in latitude between 32° S and 36.5° S, with a period of ∼1000 days, which may be a result of a non-barotropic Rossby wave. It also varied its longitudinal drift rate between −20°/day and −31°/day in approximate accord with the latitudinal gradient in the zonal wind profile, exhibiting behavior similar to that of the DS2 feature observed on Neptune [Sromovsky, L.A., Limaye, S.S., Fry, P.M., 1993. Icarus 105, 110-141]. S34 also exhibits a superimposed rapid oscillation with an amplitude of 0.57° in latitude and period of 0.7 days, which is approximately consistent with an inertial oscillation.     

Brightness power spectral distribution and waves in Jupiter's upper cloud and hazes

In this work we analyze the spatial structure of Jupiter's cloud reflectivity field in order to determine brightness periodicities and power spectra characteristics together with their relationship with Jupiter's dynamics and turbulence. The research is based on images obtained in the near-infrared (∼950 nm), blue (∼430 nm) and near-ultraviolet (∼260 nm) wavelengths with the Hubble Space Telescope in 1995 and the Cassini spacecraft Imaging Science Subsystem in 2000. Zonal reflectivity scans were analyzed by means of spatial periodograms and power spectra. The periodograms have been used to search for waves as a function of latitude. We present the values of the dominant wavenumbers for latitude bands between 32° N and 42° S. The brightness power spectra analysis has been performed in the meridional and zonal directions. The meridional analysis of albedo profiles are close to a k⁻⁵ law similarly to the wind profiles at blue and infrared wavelengths, although results differ from that in the ultraviolet. The zonal albedo analysis results in two distributions characterized by different slopes. In the near infrared and blue wavelengths, average spectral slopes are n₁=−1.3±0.4 for shorter wavenumbers (k<80), and n₂=−2.5±0.7 for greater wavenumbers, whereas for the ultraviolet n₁=−1.9±0.4 and n₂=−0.7±0.4, possibly showing a different dynamical regime. We find a turning point in the spectra between both regimes at wavenumber k∼80 (corresponding to L∼1000 km) for all wavelengths.     

Measurement of Impact Ejecta from Regolith Targets in Oblique Impacts

This paper reports the results of experiments on projectile impact into regolith targets at various impact angles. Copper projectiles of 240 mg are accelerated to 197 to 272 m s⁻¹ using an electromagnetic gun. The ejecta are detected by thin Al foil targets as secondary targets, and the resulting holes on the foil are measured to derive the spatial distribution of the ejecta. The ejecta that penetrated the foil are concentrated toward the downrange azimuths of impacting projectiles in oblique impacts. In order to investigate the ejecta velocity distribution, the nondimensional volume of ejecta with velocities higher than a given value is calculated from the spatial distribution. In the case of the vertical impact of the projectile, most ejecta have velocities lower than 24% of the projectile speed (∼50 m s⁻¹), and there are only several ejecta with velocities higher than 72 m s⁻¹. This result confirms the existence of an upper limit to the ejection velocity in the ejecta velocity distribution (Hartmann cutoff velocity) (W. K. Hartmann, 1985, Icarus63, 69-98). On the other hand, it is found that, in the oblique impacts, there are a large number of ejecta with velocities higher than the Hartmann cutoff velocity. The relative quantity of ejecta above the Hartmann cutoff velocity increases as the projectile impact angle decreases. Taking these results with the results of S. Yamamoto and A. M. Nakamura (1997, Icarus128, 160-170) from impact experiments using an impact angle of 30°, it can be concluded that the ejecta from these regolith targets exhibit a bimodal velocity distribution. Below a few tens of m s⁻¹, we see the expected velocity distribution of ejecta, but above this velocity we see a separate group of high-velocity ejecta.     

Fresh Ammonia Ice Clouds in Jupiter: I. Spectroscopic Identification, Spatial Distribution, and Dynamical Implications

We report the first spectroscopic detection of discrete ammonia ice clouds in the atmosphere of Jupiter, as discovered utilizing the Galileo Near-Infrared Mapping Spectrometer (NIMS). Spectrally identifiable ammonia clouds (SIACs) cover less than 1% of the globe, as measured in complete global imagery obtained in September 1996 during Galileo's second orbit. More than half of the most spectrally prominent SIACs reside within a small latitudinal band, extending from 2° to 7° N latitude, just south of the 5-μm hot spots. The most prominent of these are spatially correlated with nearby 5-μm-bright hot spots lying 1.5°-3.0° of latitude to the north: they reside over a small range of relative longitudes on the eastward side of hot spots, about 37% of the longitudinal distance to the next hot spot to the east. This strong correlation between the positions of hot spots and the most prominent equatorial SIACs suggests that they are linked by a common planetary wave. Good agreement is demonstrated between regions of condensation predicted by the Rossby wave model of A. J. Friedson and G. S. Orton (1999, Bull. Am. Astron. Assoc31, 1155-1156) and the observed longitudinal positions of fresh ammonia clouds relative to 5-μm hot spots. Consistency is also demonstrated between (1) the lifetime of particles as determined by the wave phase speed and cloud width and (2) the sedimentation time for 10-μm radius particles consistent with previously reported ammonia particle size by T. Y. Brooke et al. (1998, Icarus136, 1-13). A young age (<two days) for most SIAC cloud particles is indicated. To the south, the most prominent SIACs are located to the northwest of the Great Red Spot, in a region where a westward flow of jovian air, diverted approximately 10° of latitude northward by the Great Red Spot, encounters a large eastward flow. SIACs have been observed repeatedly by NIMS at this location during Galileo's first four years in Jupiter orbit. It is speculated that due to the three-dimensional interactions of these flows, relatively large amounts of ammonia gas are steadily transported from the sub-cloud troposphere (below the ∼600-mbar level) to the high troposphere, nearly continuously forming fresh ammonia ice clouds to the northwest of the Great Red Spot.     

Detection of the Forbidden SO aΔ→XΣ Rovibronic Transition on Io at 1.7 μm

We report the discovery of the forbidden electronic a¹Δ→X³Σ⁻ transition of the SO radical on Io at 1.7 μm with the W. M. Keck II telescope on 24 September 1999 (UT), while the satellite was eclipsed by Jupiter. The shape of the SO emission band suggests a rotational temperature of ∼1000 K; i.e., the gas is extremely hot. We interpret the observed emission rate of ∼2×10²⁷ photons s⁻¹ to be caused by SO molecules in the excited a¹Δ state being directly ejected from the vent at a thermodynamic quenching temperature of ∼1500 K, assuming a SO/SO₂ abundance ratio of ∼0.1 and a total venting rate of ∼10³¹ molecules s⁻¹ (Strobel and Wolven 2001, Astrophys. Space Sci. 277, 1-17). The shape of our complete (1.6-2.5 μm) spectrum suggests that the volcano Loki contains a small (∼2 km²) hot spot at 960±12 K, as well as a larger (∼50 km²) area at 640±5 K.     

Improved measurement of Asteroid (4) Vesta’s rotational axis orientation

We report an improved measurement of the rotational axis orientation of Asteroid (4) Vesta. By analyzing and combining all previous measurements using a limb-fitting technique from ground/HST data collected from 1983 to 2006, we derive a pole solution of (RA = 304.5°, Dec = 41.5°). Images of Vesta acquired with the Wide Field Camera 3 onboard the Hubble Space Telescope (HST) in February 2010 are combined with images from the Wide Field Planetary Camera 2 on HST obtained in 1994, 1996, and 2007 at similar spatial resolution and wavelengths to perform new measurements. Control point stereogrammetry returns a pole solution of (305.1°, 43.4°). An alternate method tracks surface features and fits their projected paths with ellipses to determine a great circle containing the pole for each HST observation. Combined, the four great circles yield a pole solution of (309.3°, 41.9°). These three solutions obtained with almost independent methods are within 3.5° of each other, suggesting a robust solution. Combining the results from all three techniques, we propose an improved value of the rotational axis of Vesta as RA = 305.8° ± 3.1°, Dec = 41.4° ± 1.5° (1-σ error). This new solution changes from (301°, 41°) reported by Thomas et al. (Thomas, P.C., Binzel, R.P., Gaffey, M.J., Zellner, B.H., Storrs, A.D., Wells, E. [1997a]. Icarus 128, 88-94) by 3.6°, and from (306°, 38°) reported by Drummond and Christou (Drummond, J.D., Christou, J. [2008]. Icarus 197, 480-496) by 3.4°. It changes the obliquity of Vesta by up to ∼3°, but increases the Sun-centered RA of Vesta at equinox by ∼8°, and postpones the date of equinox by ∼35 days. The change of the pole position is less than the resolution of all previous images of Vesta, and should not change the main science conclusions of previous research about Vesta.     

Near-Infrared Spectrophotometry of Phobos and Deimos

We have observed the leading and trailing hemispheres of Phobos from 1.65 to 3.5 μm and Deimos from 1.65 to 3.12 μm near opposition. We find the trailing hemisphere of Phobos to be brighter than its leading hemisphere by 0.24±0.06 magnitude at 1.65 μm and brighter than Deimos by 0.98±0.07 magnitude at 1.65 μm. We see no difference larger than observational uncertainties in spectral slope between the leading and trailing hemispheres when the spectra are normalized to 1.65 μm. We find no 3-μm absorption feature due to hydrated minerals on either hemisphere to a level of ∼5-10% on Phobos and ∼20% on Deimos. When the infrared data are joined to visible and near-IR data obtained by previous workers, our data suggest the leading (Stickney-dominated) side of Phobos is best matched by T-class asteroids. The spectral slope of the trailing side of Phobos and leading side of Deimos are bracketed by the D-class asteroids. The best laboratory spectral matches to these parts of Phobos are mature lunar soils and heated carbonaceous chondrites. The lack of 3-μm absorption features on either side of Phobos argues against the presence of a large interior reservoir of water ice according to current models of Phobos' interior (F. P. Fanale and J. R. Salvail 1989, Geophys. Res. Lett.16, 287-290; Icarus88, 380-395).     

The nature of Neptune’s increasing brightness: evidence for a seasonal response

Hubble Space Telescope (HST) observations in August 2002 show that Neptune’s disk-averaged reflectivity increased significantly since 1996, by 3.2 ± 0.3% at 467 nm, 5.6 ± 0.6% at 673 nm, and 40 ± 4% in the 850-1000 nm band, which mainly results from dramatic brightness increases in restricted latitude bands. When 467-nm HST observations from 1994 to 2002 are added to the 472-nm ground-based results of Lockwood and Thompson (2002, Icarus 56, 37-51), the combined disk-averaged variation from 1972 to 2002 is consistent with a simple seasonal model having a hemispheric response delay relative to solar forcing of ∼30 years (∼73% of a full season).     

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.