Radar and infrared observations of binary near-Earth Asteroid 2002 CE26

We observed near-Earth Asteroid (NEA) 2002 CE26 in August and September 2004 using the Arecibo S-band (2380-MHz, 12.6-cm) radar and NASA's Infrared Telescope Facility (IRTF). Shape models obtained based on inversion of our delay-Doppler images show the asteroid to be 3.5±0.4 km in diameter and spheroidal; our corresponding nominal estimates of its visual and radar albedos are 0.07 and 0.24, respectively. Our IRTF spectrum shows the asteroid to be C-class with no evidence of hydration. Thermal models from the IRTF data provide a size and visual albedo consistent with the radar-derived estimate. We estimate the spin-pole to be within a few tens of degrees of λ=317°, β=−20°. Our radar observations reveal a secondary approximately 0.3 km in diameter, giving this binary one of the largest size differentials of any known NEA. The secondary is in a near-circular orbit with period 15.6±0.1 h and a semi-major axis of 4.7±0.2 km. Estimates of the binary orbital pole and secondary rotation rate are consistent with the secondary being in a spin-locked equatorial orbit. The orbit corresponds to a primary mass of M=1.95±0.25×¹⁰¹³ kg, leading to a primary bulk density of , one of the lowest values yet measured for a main-belt or near-Earth asteroid.     

Photometric study of Centaur (60558) 2000 EC98 and trans-neptunian object (55637) 2002 UX25 at different phase angles

We present photometric observations of Centaur (60558) 2000 EC₉₈ and trans-neptunian object (55637) 2002 UX₂₅ at different phase angles and with different filters (mainly R but also V and B for some data). Results for 2000 EC₉₈ are: (i) a rotation period of 26.802±0.042 h if a double-peaked lightcurve is assumed, (ii) a lightcurve amplitude of 0.24±0.06 for the R band, (iii) a phase curve with H=9.03±0.01 and G=−0.39±0.08 (R filter) and H=9.55±0.04 and G=−0.50±0.35 (V filter) or a slope of (R filter) and 0.22±0.06 (V filter), (iv) the color indices B-V=0.76±0.15 and V-R=0.51±0.09 (for α=0.1-0.5°) and 0.55±0.08 (for α=1.4-1.5°). The rotation period is amongst the longest ever measured for Centaurs and TNOs. We also show that our photometry was not contaminated by any cometary activity down to magnitude ?27/arcsec². For 2002 UX₂₅ the results are: (i) a rotation period of 14.382±0.001 h or 16.782±0.003 h (if a double-peaked lightcurve is assumed) (ii) a lightcurve amplitude of 0.21±0.06 for the R band (and the 16.782 h period), (iii) a phase curve with H=3.32±0.01 and G=+0.16±0.18 or a slope of (R filter), (iv) the color indices B-V=1.12±0.26 and V-R=0.61±0.12. The phase curve reveals also a possible very narrow and bright opposition surge. Because such a narrow surge appears only for one point it needs to be confirmed.     

Radar observations and a physical model of Asteroid 1580 Betulia

We report Arecibo (2380-MHz, 13-cm) observations of Asteroid 1580 Betulia in May-June 2002. We combine these continuous-wave Doppler spectra and delay-Doppler images with optical lightcurves from the 1976 and 1989 apparitions in order to estimate Betulia's shape and spin vector. We confirm the spin vector solution of Kaasalainen et al. [Kaasalainen, M., and 21 colleagues, 2004. Icarus 167, 178-196], with sidereal period P=6.13836 h and ecliptic pole direction (λ,β)=(136°,+22°), and obtain a model that resembles the Kaasalainen et al. convex-definite shape reconstruction but is dominated by a prominent concavity in the southern hemisphere. We find that Betulia has a maximum breadth of 6.59±0.66 km and an effective diameter of 5.39±0.54 km. These dimensions are in accord with reanalyzed polarimetric and radar data from the 1970s. Our effective diameter is 15% larger than the best radiometric estimate of Harris et al. [Harris, A.W., Mueller, M., Delbó, M., Bus, S.J., 2005. Icarus 179, 95-108], but this difference is much smaller than the size differences between past models. Considering orbits of test particles around Betulia, we find that this asteroid's unusual shape results in six equilibrium points close to its equatorial plane rather than the usual four points; two of these six points represent stable synchronous orbits while four are unstable. Betulia's close planetary encounters can be predicted for over four thousand years into the future.     

Optical and thermal infrared observations of six near-Earth asteroids in 2002

We present thermal infrared photometry and spectrophotometry of six Near-Earth Asteroids (NEAs) using the 3.8 m United Kingdom Infrared Telescope (UKIRT) together with quasi-simultaneous optical observations of five NEAs taken at the 1.0 m Jacobus Kapteyn Telescope (JKT). For Asteroid (6455) 1992 HE we derive a rotational period P=2.736±0.002 h, and an absolute visual magnitude H=14.32±0.24. For Asteroid 2002 HK₁₂ we derive . The Standard Thermal Model (STM), the Fast Rotating Model (FRM) and the Near-Earth Asteroid Thermal Model (NEATM) have been fitted to the measured fluxes to derive albedos and effective diameters. The derived geometric albedos and effective diameters are (6455) 1992 HE: pv=0.26±0.08, Deff=3.55±0.5 km; 1999 HF₁: pv=0.18±0.07, ; 2000 ED₁₀₄: pv=0.18±0.05, Deff=1.21±0.2 km; 2002 HK₁₂: , Deff=0.62±0.2 km; 2002 NX₁₈: pv=0.031±0.009, Deff=2.24±0.3 km; 2002 QE₁₅: , Deff=1.94±0.4 km. The limitations of using the NEATM to observe NEAs at high phase angles are discussed.     

Physical and compositional studies of Comet 81P/Wild 2 at multiple apparitions

We present analyses and results from both narrowband photometry and CCD imaging of Comet 81P/Wild 2 from multiple apparitions, obtained in support of the Stardust mission. These data include photometric measurements from 12 days before the encounter and imaging from 3 days after. Using narrowband photometry from the different apparitions, we analyzed the dust and gas production rates as a function of heliocentric distance, finding a substantial seasonal effect where the production of OH, NH, and dust peaks 11-12 weeks before perihelion. The CN, C₂, and C₃ production show no such asymmetry, suggesting that there may be heterogeneities among different sources on the nucleus. The water production peaked at a level of approximately in 1997. A comparison of the relative abundances of minor gas species places Wild 2 in the “depleted” category in the A'Hearn et al. (1995, Icarus 118, 223) taxonomic classifications. Continuum measurements at multiple wavelengths indicate that the comet has a low dust-to-gas ratio, with moderately reddened dust. In our images we see a dust tail, an anti-tail and two well-defined jets. The primary jet, which persists for several months and is roughly aligned with the spin axis, has a source latitude >+75°, while the secondary jet is located on the opposite hemisphere between −37° and −62°. We used the apparent position angle of the primary jet to determine the pole orientation, α=281±5°, δ=+13±7°, and surmise that the nucleus is likely in a state of simple rotation. The primary source is continuously illuminated when Wild 2 is inbound and turns away from the Sun at about the time that the comet reaches perihelion, explaining the seasonal effects in the production rates. We measured lightcurves on several observing runs but saw no significant modulation, so no constraints can be set on the rotation rate. Images at different wavelengths show that the jets have the same colors as the dust in other regions in the coma and tail, indicating that the grain properties are similar throughout the coma. Radial profiles of the coma were measured in various directions on a number of different observing runs, and we discuss the findings from these measurements. Finally, we compare our results with other published data and attempt to predict future times at which observations should be obtained to help constrain additional properties.     

A global solution for the Mars static and seasonal gravity, Mars orientation, Phobos and Deimos masses, and Mars ephemeris

With the collection of six years of MGS tracking data and three years of Mars Odyssey tracking data, there has been a continual improvement in the JPL Mars gravity field determination. This includes the measurement of the seasonal changes in the gravity coefficients (e.g., , , , , , ) caused by the mass exchange between the polar ice caps and atmosphere. This paper describes the latest gravity field MGS95J to degree and order 95. The improvement comes from additional tracking data and the adoption of a more complete Mars orientation model with nutation, instead of the IAU 2000 model. Free wobble of the Mars' spin axis, i.e. polar motion, has been constrained to be less than 10 mas by looking at the temporal history of and . A strong annual signature is observed in , and this is a mixture of polar motion and ice mass redistribution. The Love number solution with a subset of Odyssey tracking data is consistent with the previous liquid outer core determination from MGS tracking data [Yoder et al., 2003. Science 300, 299-303], giving a combined solution of k₂=0.152±0.009 using MGS and Odyssey tracking data. The solutions for the masses of the Mars' moons show consistency between MGS, Odyssey, and Viking data sets; Phobos GM=(7.16±0.005)×¹⁰⁻⁴ km³/s² and Deimos GM=(0.98±0.07)×¹⁰⁻⁴ km³/s². Average MGS orbit errors, determined from differences in the overlaps of orbit solutions, have been reduced to 10-cm in the radial direction and 1.5 m along the spacecraft velocity and normal to the orbit plane. Hence, the ranging to the MGS and Odyssey spacecraft has resulted in position measurements of the Mars system center-of-mass relative to the Earth to an accuracy of one meter, greatly reducing the Mars ephemeris errors by several orders of magnitude, and providing mass estimates for Asteroids 1 Ceres, 2 Pallas, 3 Juno, 4 Vesta, and 324 Bamberga.     

Non-constant rotation period of Comet C/2001 K5 (LINEAR)

The article presents results of CCD photometry in R-band of a dynamically new Comet C/2001 K5 (LINEAR), obtained at a heliocentric distance of about 5.6 AU, after the perihelion passage. Being so distant from the Sun, this comet was extremely active (Afρ close to 2000 cm), exhibiting quite well developed dust coma and tail. During the observations, general photometric behavior of the comet with heliocentric distance r was well described by the 2.5nlog(r) function with coefficient n=5. The radial profiles of the coma were found to be undulated, with mean slope of the dependence between cometary magnitude and 2.5log of aperture radius (at comet distance) equal to . The light curve of Comet LINEAR exhibited short-term variability which we attributed to cyclic changes of dust emission, induced by nucleus rotation. Model computations by some authors have revealed that active comets can change their spin status quite substantially even during a single orbital revolution. Thus, attempting to search for a rotation frequency, we have modified the classical PDM approach by including the spin acceleration term. Such DynamicalPDM (DPDM) method revealed the most reliable solution for the frequency f₀=0.019048±0.000013 h⁻¹ and its first time-derivative (index “zero” denotes reference to the mid time of the whole observing run), indicating a rapid spin-down of the nucleus. These parameters are equivalent to the rotation period of 52.499±0.036 h and its relative increment of 0.02729±0.00013. We present the most probable evolution of the rotation frequency of Comet LINEAR, based on the results of periodicity analysis and a simple, almost parameter independent, dynamical model of nucleus rotation. It is also shown that the DPDM may be an effective tool for determination of a nucleus radius, which provided us with the value of 1.53±0.25 km for Comet LINEAR.     

Arecibo radar observations of Phobos and Deimos

In November 2005, we observed the moons of Mars using the Arecibo 2380-MHz (13-cm) radar, obtaining a result for the OC radar albedo of Phobos (0.056±0.014) consistent with its previously reported radar albedo and implying an upper bound on its near-surface bulk density of . We detected Deimos by radar for the first time, finding its OC radar albedo to be 0.021±0.006, implying an upper bound on its near-surface density of , consistent with a high-porosity regolith. We briefly discuss reasons for these low radar albedos, Deimos' being possibly the lowest of any Solar System body yet observed by radar.     

Mars high resolution gravity fields from MRO, Mars seasonal gravity, and other dynamical parameters

With 2 years of tracking data collection from the MRO spacecraft, there is noticeable improvement in the high frequency portion of the spherical harmonic Mars gravity field. The new JPL Mars gravity fields, MRO110B and MRO110B2, show resolution near degree 90. Additional years of MGS and Mars Odyssey tracking data result in improvement for the seasonal gravity changes which compares well to global circulation models and Odyssey neutron data and Mars rotation and precession (). Once atmospheric dust is accounted for in the spacecraft solar pressure model, solutions for Mars solar tide are consistent between data sets and show slightly larger values (k₂ = 0.164 ± 0.009, after correction for atmospheric tide) compared to previous results, further constraining core models. An additional 4 years of Mars range data improves the Mars ephemeris, determines 21 asteroid masses and bounds solar mass loss (dGM_Sun/dt < 1.6 × 10⁻¹³ GM_Sun year⁻¹).     

A strong vortex in Saturn's South Pole

Saturn's southern pole was observed at high resolution by the Cassini Imaging Science Subsystem (ISS) during the spacecraft insertion orbit in July 2004. Cloud tracking of individual features on images taken at a wavelength of 938 nm reveal the existence of a strong polar vortex enclosed by a jet with maximum speed of relative to System III rotation frame, and peak at 87 °S planetographic latitude. Radiative transfer models of the reflected light, based on the Cassini images complemented by Hubble Space Telescope images from March 2004, indicate that the aerosol particles in the vortex are structured vertically in three detached layers. We find two hazes and one dense cloud distributed in altitude between ∼500 mbar (top of the dense cloud) and few mbar (top of the stratospheric haze), spanning a vertical altitude range of ∼200 km. The vortex area coincides with a thermal hot spot recently reported, indicating that winds decrease with altitude above polar clouds.     

On the vertical wind shear of Saturn's Equatorial Jet at cloud level

In this work we analyze and compare the vertical cloud structure of Saturn's Equatorial Zone in two different epochs: the first one close to the Voyagers flybys (1979-1981) and the second one in 2004, when the Cassini spacecraft entered its orbit around the planet. Our goal is to retrieve the altitude of cloud features used as zonal wind tracers in both epochs. We reanalyze three different sets of photometrically calibrated published data: ground-based in 1979, Voyager 2 PPS and ISS observations in 1981, and we analyze a new set of Hubble Space Telescope images for 2004. For all situations we reproduced the observed reflectivity by means of a similar vertical model with three layers. The results indicate the presence of a changing tropospheric haze in 1979-1981 (Ptop∼100 mbar, τ∼10) and in 2004 (Ptop∼50 mbar, τ∼15) where the tracers are embedded. According to this model the Voyager 2 ISS images locate cloud tracers moving with zonal velocities of 455 to 465 (±2) m/s at a pressure level of 360 ± 140 mbar. For HST observations, our previous works had showed cloud tracers moving with zonal wind speeds of 280±10 m/s at a pressure level of about 50±10 mbar. All these values are calculated in the same region (3°±2° N). This speed difference, if interpreted as a vertical wind shear, requires a change of per scale height, two times greater than that estimated from temperature observations. We also perform an initial guess on Cassini ISS vertical sounding levels, retrieving values compatible with HST ones and Cassini CIRS derived vertical wind shear, but not with Voyager wind measurements. We conclude that the wind speed velocity differences measured between 1979-1981 and 2004 cannot be explained as a wind shear effect alone and demand dynamical processes.     

Interferometric millimeter observations of water vapor on Mars and comparison with Mars Express measurements

We present interferometric mapping of the 225.9-GHz HDO and 203.4-GHz lines on Mars obtained with the IRAM Plateau de Bure facility (PdBI). The observations were performed during martian year 28 (MY28), at L_s=320.3° for the HDO line, and at L_s=324.3° for the line. The HDO line is detected at the eastern (morning) and western (evening) limbs of the northern hemisphere, corresponding to a water column density in the range 3-6 pr.-μm. The line is not detected, which is compatible with the column densities derived from the HDO line. Quasi-simultaneous far infrared measurements obtained by the Planetary Fourier Spectrometer (PFS) onboard the Mars Express spacecraft confirm our PdBI results, yielding a 5±1 pr.-μm meridionally constant water column abundance.Such a low water abundance during the southern mid-autumn of MY28 does not correspond to the standard martian climatology as observed during the previous years. It was however already retrieved from near-infrared observations performed by the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) onboard the Mars Reconnaissance Orbiter spacecraft [Smith, M.D., Wolff, M.J., Clancy, R.T., Murchie, S.L. 2009. CRISM observations of water vapor and carbon monoxide. J. Geophys. Res. 114, doi: 10.1029/2008JE003288]. Our observations thus confirm that the planet-encircling dust storm that occurred during MY28 significantly affected the martian water cycle. Our observations also demonstrate the usefulness of interferometric submillimeter observations to survey the martian water cycle from ground-based facilities.     

Interplanetary dust detected by the Cassini CDA Chemical Analyser

During its cruise phase, prior to encountering Jupiter, the Cosmic Dust Analyser (CDA) onboard the Cassini spacecraft returned time of flight mass spectra (TOF MS) of two interplanetary dust particles. Both particles were found to be iron-rich, with possible traces of hydrogen, carbon, nickel, chromium, manganese, titanium, vanadium and minor silicates. Carbon, hydrogen, oxygen and potassium are also present as possible contaminants of the impact target of CDA. Silicates and magnesium do not feature predominantly in the spectra; this is surprising considering the expected dominance of silicate-rich minerals in interplanetary dust particles. The particle masses are and . The corresponding radii ranges for the particles, assuming densities from 7874-2500 kg m⁻³ are 0.7-4 μm and 2.6-6.8 μm, respectively. With the same density assumptions the β values (ratio of radiation pressure to gravitational force) are estimated as 0.027-0.21 and 0.016-0.06 respectively, allowing possible orbits to be calculated. The resulting orbits are bound and prograde with semi-major axes, eccentricities and inclinations in the region of 0.3-1.26 AU, 0.4-1.0 and 0-60° for the first particle and 0.8-2.5 AU, 0.2-0.9 and 0-30° for the second. The more probable orbits within these ranges indicate that the first particle is in an Aten-like orbit, whilst the second particle is in an Apollo-like orbit, despite both grains having very similar, predominantly metallic compositions. Other possible orbital solutions for both particles encompass orbits which more closely resemble those of Jupiter-family comets.     

The orbit, mass, and albedo of transneptunian binary (66652) 1999 RZ253

We have observed (66652) 1999 RZ₂₅₃ with the Hubble Space Telescope at seven separate epochs and have fit an orbit to the observed relative positions of this binary. Two orbital solutions have been identified that differ primarily in the inclination of the orbit plane. The best fit corresponds to an orbital period, days, semimajor axis a=4660±170 km and orbital eccentricity e=0.460±0.013 corresponding to a system mass m=3.7±0.4×¹⁰¹⁸ kg. For a density of the albedo at 477 nm is p₄₇₇=0.12±0.01, significantly higher than has been commonly assumed for objects in the Kuiper belt. Multicolor, multiepoch photometry shows this pair to have colors typical for the Kuiper belt with a spectral gradient of 0.35 per 100 nm in the range between 475 and 775 nm. Photometric variations at the four epochs we observed were as large as 12±3% but the sampling is insufficient to confirm the existence of a lightcurve.     

Direct measurement of the size, shape, and pole of 511 Davida with Keck AO in a single night

Using the high-quality data set of 165 images taken at 11 epochs over the 5.13 h rotation of the large C-type Asteroid 511 Davida, we find the dimensions of its triaxial ellipsoid model to be 357±2×294±2×231±50 km. The images were acquired with the adaptive optics system on the 10 m Keck II telescope on December 27, 2002. The a and b diameters are much better determined than previously estimated from speckle interferometry and indirect measurements, and our mean diameter, (abc)^1/3=289±21 km, is 19% below previous estimates. We find the pole to lie within 2° of [RA=295°; Dec=0°] or in Ecliptic coordinates [λ=297°; β=+21°], a significant improvement to the pole direction. Otherwise, previous determinations of the axial ratios agree with our new results. These observations illustrate that our technique of finding the dimensions and pole of an asteroid from its changing projected size and shape is very powerful because it can be done in essentially one night as opposed to decades of lightcurves. Average departures of 3% (5 km) of the asteroid's mean radius from a smooth outline are detected, with at least two local positive-relief features and at least one flat facet showing approximately 15 km deviations from the reference best-fit ellipsoid. The facet is reminiscent of large global-scale craters on Asteroid 253 Mathilde (also a C-type) when seen edge-on in close-up images from the NEAR mission flyby. We show that giant craters (up to 150 km diameter, the size of the largest facets seen on Davida) can be expected from the impactor size distribution, without likelihood of catastrophic disruption of Davida.     

New near-aphelion light curves of Comet 2P/Encke

We present new, near-aphelion, time series of photometry of Comet 2P/Encke in Cousins-R band. With these light curves we find that the dominant, synodic rotational periodicity is either P₀=11.079±0.009 h or 2P₀=22.158±0.012 h. This is in contrast to data from the 1980s published by others that are consistent with 15.08- and 22.6-h periods. Those periods do not satisfy our phased light curves, and also the 1980s data are not easily reconciled with our periods. This could be due to P/Encke having non-principal axis rotation or due to a drift in the rotation period caused by outgassing torques. We observed the comet at five epochs: July, August, September, and October 2001, and September 2002, and the comet was at times intrinsically brighter than expected for a bare nucleus, due to an apparent contribution from an unresolved coma. Three-quarters of the data were obtained in the second and fifth epochs, and we analyzed these two time series using both the phase-dispersion minimization and “WindowCLEAN” techniques. At both epochs and with both techniques strong periodicities were found near frequencies and . By then using visual inspection of the phased light curves to corroborate these frequencies, and by using the data from the other three epochs to properly align light curve features, we were able to derive P₀ and 2P₀ as the only solutions that satisfy all our observations. The periodicity due to f₁ is clearly seen in our data, but we cannot tell from our data alone whether it is a manifestation of the nucleus's shape, non-principal axis rotation, or both.     

Dark red debris from three short-period comets: 2P/Encke, 22P/Kopff, and 65P/Gunn

We present observations of the extended dust structures near the orbits of three short-period comets: 2P/Encke, 22P/Kopff, and 65P/Gunn. The dust trails were originally discovered by the Infrared Astronomical Satellite (IRAS). Our observations were made using wide-field optical CCD cameras on the University of Hawaii 2.24-m telescope, the Canada-France-Hawaii 3.6-m telescope, and the Kiso 1.05-m Schmidt telescope. We compared the observed images with models and found that the extended structures seen around 2P/Encke and 22P/Kopff before perihelion passage were most likely “dust trails,” whereas images taken after perihelion passage show a high contamination by recently released particles (i.e., particles in Neck-Line structures are visible). We could not confirm the existence of a dust trail from 65P/Gunn within the field of view of the camera used. The effective sizes of the particles responsible for the scattered light were estimated at 1-100 mm (2P/Encke), 1-10 mm (22P/Kopff), and 100 μm-1 mm (65P/Gunn), respectively, which is consistent with previous studies of dust trails made with infrared space telescopes and optical telescopes. We evaluated the mass loss rates of these comets, averaged over their orbits, as reaching (2P/Encke), (22P/Kopff), and (65P/Gunn). These values are consistent with previous work. Therefore, the total amount of material ejected from these three comets is , which would contribute a considerable fraction of the lost within 1 AU that needs to be replaced if the zodiacal cloud is to be maintained in a steady state. We also found that the particles in the dust structures are significantly redder than the Sun and the zodiacal light, and might be redder than the average short-period comet nuclei. Specifically, the reflectivity gradients of 2P/Encke, 22P/Kopff, and 65P/Gunn are 13±7 (% ¹⁰³ Å⁻¹), 20±5 (% ¹⁰³ Å⁻¹), and 15±4 (% ¹⁰³ Å⁻¹), respectively. We examined the change in color with distance from the nucleus. No clear correlation was detected for 2P/Encke or 22P/Kopff to an accuracy of 3-11%, while the 65P/Gunn tail did show color variation, becoming redder with increasing distance from the nucleus. This dark red material, consisting of particles of sand-cobble size, has marginally escaped from the nuclei and will evolve into finer-grained interplanetary dust particles after subsequent collisions.     

Mass and interior of Enceladus from Cassini data analysis

Gravity results are available from radio Doppler data acquired by the Deep Space Network during the encounter of the Cassini spacecraft with Enceladus in February 2005. We report the mass of Enceladus to be (1.0798±0.0016)×¹⁰²⁰ kg, which implies a density of . For a core made of hydrated silicates with a density of 2500 kg m⁻³ the core radius is ∼190 km and the quadrupole moment C₂₂∼1.4×¹⁰⁻³. If Enceladus is in hydrostatic equilibrium, the larger than previously anticipated density implies that the recently proposed secondary spin-orbit resonance cannot be present. Therefore, the source of endogenic activity of Enceladus remains unexplained.