Radar observations of Comet 2P/Encke during the 2003 apparition

The nucleus of Comet 2P/Encke was detected with the Arecibo radar during the close approach of November, 2003, making this the first comet to yield radar detections at two different apparitions. Although the measured radar cross section of 1.0 km² was close to that obtained in 1980, the Doppler bandwidth was nearly four times larger. Most of this bandwidth difference can simply be attributed to a different observing aspect relative to the spin axis proposed by Sekanina [1988, Astron. J. 95, 911] and Festou and Barale [2000, Astron. J. 119, 3119]. Comparison of the 2003 Doppler bandwidth with infrared-based size estimates supports an 11-h dominant rotation period and excludes slower 15- and 22-h periods that have also been suggested. If one assumes a short-axis-mode rotation with an 11-h period, then the Doppler bandwidth indicates that the nucleus is an oblong object with a long-axis dimension of 9 km. The estimated radar albedo of 0.05 is similar to that measured for C/IRAS-Araki-Alcock, providing further evidence that comet nuclei have relatively low surface densities of ∼0.5-1.0 g cm⁻³. No broadband echo component was detected from large coma grains despite predictions, based on optical/infrared models, that such a component might be detectable.     

Rotation and color properties of the nucleus of Comet 2P/Encke

We present results from CCD observations of Comet 2P/Encke acquired at Steward Observatory's 2.3 m Bok Telescope on Kitt Peak. The observations were carried out in October 2002 when the comet was near aphelion. Rotational lightcurves in B-, V-, and R-filters were acquired over two nights of observations, and analysed to study the physical and color properties of the nucleus. The average apparent R-filter magnitude across both nights corresponds to a mean effective radius of 3.95±0.06 km, and this value is similar to that found for the V- and B-filters. Taking the observed brightness range, we obtain a/b?1.44±0.06 for the semi-axial ratio of Encke's nucleus. Applying the axial ratio to the R-filter photometry gives nucleus semi-axes of [3.60±0.09]×[5.20±0.13] km, using the empirically-derived albedo and phase coefficient. No coma or tail was seen despite deep imaging of the comet, and flux limits from potential unresolved coma do not exceed a few percent of the total measured flux, for standard coma models. This is consistent with many other published data sets taken when the comet was near aphelion. Our data includes the first detailed time series multi-color measurements of a cometary nucleus, and significant color variations were seen on October 3, though not repeated on October 4. The average color indices across both nights are: (V−R)=0.39±0.06 and (B−V)=0.73±0.06 (). We analysed the R-filter time-series photometry using the method of Harris et al. [Harris, A.W., Young, J.W., Bowell, E., Martin, L.J., Millis, R.L., Poutanen, M., Scaltriti, F., Zappala, V., Schober, H.J., Debehogne, H., Zeigler, K.W., 1989. Icarus 77, 171-186] to constrain the rotation period of the comet's nucleus, and find that a period of ∼11.45 h will satisfy the data, however the errors bars are large. We have successfully linked our data with the September 2002 data from Fernández et al. [Fernández, Y.R., Lowry, S.C., Weissman, P.R., Mueller, B.E.A., Samarasinha, N.H., Belton, M.J.S., Meech, K.J., 2005. Icarus 175, 194-214]—taken just 2-3 weeks before the current data set—and we show that a rotation period of just over 11 h works extremely well for the combined data set. The resulting best-fit period is 11.083±0.003 h, consistent with the Fernández et al. value.     

Near-nucleus photometry of comets using archived NEAT data

Though optimized to discover and track fast moving Near-Earth Objects (NEOs), the Near-Earth Asteroid Tracking (NEAT) survey dataset can be mined to obtain information on the comet population observed serendipitously during the asteroid survey. We have completed analysis of over 400 CCD images of comets obtained during the autonomous operations of two 1.2-m telescopes: the first on the summit of Haleakala on the Hawaiian island of Maui and the second on Palomar Mountain in southern California. Photometric calibrations of each frame were derived using background catalog stars and the near-nucleus comet photometry measured. We measured dust production and normalized magnitudes for the coma and nucleus in order to explore cometary activity and comet size-frequency distributions. Our data over an approximately two-year time frame (2001 August-2003 February) include 52 comets: 12 periodic, 19 numbered, and 21 non-periodic, obtained over a wide range of viewing geometries and helio/geocentric distances. Nuclear magnitudes were estimated for a subset of comets observed. We found that for low-activity comets (Afρ<100 cm) our model gave reasonable estimates for nuclear size and magnitude. The slope of the cumulative luminosity function of our sample of low-activity comets was 0.33 ± 0.04, consistent with the slope we measured for the Jupiter-family cometary nuclei collected by Fernández et al. [Fernández, J.A., Tancredi, G., Rickman, H., Licandro, J., 1999. Astron. Astrophys. 392, 327-340] of 0.38 ± 0.02. Our slopes of the cumulative size distribution α=1.50±0.08 agree well with the slopes measured by Whitman et al. [Whitman, K., Morbidelli, A., Jedicke, R., 2006. Icarus 183, 101-114], Meech et al. [Meech, K.J., Hainaut, O.R., Marsden, B.G., 2004. Icarus 170, 463-491], Lowry et al. [Lowry, S.C., Fitzsimmons, A., Collander-Brown, S., 2003. Astron. Astrophys. 397, 329-343], and Weissman and Lowry [Weissman, P.R., Lowry, S.C., 2003. Lunar Planet. Sci. 34. Abstract 34].     

Nucleus properties of Comet 67P/Churyumov-Gerasimenko estimated from non-gravitational force modeling

By considering model comet nuclei with a wide range of sizes, prolate ellipsoidal shapes, spin axis orientations, and surface activity patterns, constraints have been placed on the nucleus properties of the primary Rosetta target, Comet 67P/Churyumov-Gerasimenko. This is done by requiring that the model bodies simultaneously reproduce the empirical nucleus rotational lightcurve, the water production rate as function of time, and non-gravitational changes (per apparition) of the orbital period (ΔP), longitude of perihelion (Δ?), and longitude of the ascending node (ΔΩ). Two different thermophysical models are used in order to calculate the water production rate and non-gravitational force vector due to nucleus outgassing of the model objects. By requiring that the nominal water production rate measurements are reproduced as well as possible, we find that the semi-major axis of the nucleus is close to 2.5 km, the nucleus axis ratio is approximately 1.4, while the spin axis argument is either 60°±15° or 240°±15°. The spin axis obliquity can only be preliminarily constrained, indicating retrograde rotation for the first argument value, and prograde rotation for the second suggested spin axis argument. A nucleus bulk density in the range 100-370 kg m⁻³ is found for the nominal ΔP, while an upper limit of 500 kg m⁻³ can be placed if the uncertainty in ΔP is considered. Both considered thermophysical models yield the same spin axis, size, shape, and density estimates. Alternatively, if calculated water production rates within an envelope around the measured data are considered, it is no longer possible to constrain the size, shape, and spin axis orientation of the nucleus, but an upper limit on the nucleus bulk density of 600 kg m⁻³ is suggested.     

Prometheus and Pandora: masses and orbital positions during the Cassini tour

Hubble Space Telescope (HST) images of Prometheus and Pandora show longitude discrepancies of about 20° with respect to the Voyager ephemerides, with an abrupt change in mean motion at the end of 2000 (French et al., 2003, Icarus 162, 143-170; French and McGhee, 2003, Bull. Am. Astron. Soc. 34, 06.07). These discrepancies are anti-correlated and arise from chaotic interactions between the two moons, occurring at interval of 6.2 yr, when their apses are anti-aligned (Goldreich and Rappaport, 2003a, Icarus 162, 391-399). This behavior is attributed to the overlap of four 121:118 apse-type mean motion resonances (Goldreich and Rappaport, 2003b, Icarus 166, 320-327). We study the Prometheus-Pandora system using a Radau-type integrator taking into account Saturn's oblateness up to and including terms in J₆, plus the effects of the major satellites. We first confirm the chaotic behavior of Prometheus and Pandora. By fitting the numerical integrations to the HST data (French et al., 2003, Icarus 162, 143-170; French and McGhee, 2003, Bull. Am. Astron. Soc. 34, 06.07), we derive the satellite masses. The resulting GM values (with their standard 3-σ errors) for Prometheus and Pandora are respectively and . Using the nominal shape of the two moons (Thomas, 1989, Icarus 77, 248-274), we derive Prometheus and Pandora's densities, 0.40^+0.03_−0.07 and 0.49^+0.05_−0.09 g cm⁻³, respectively. Our numerical fits also enable us to constrain the time of the latest apse anti-alignment in 2000. Finally, using our fit, we predict the orbital positions of the two satellites during the Cassini tour, and provide a lower limit of the uncertainties due to chaos. These uncertainties amount to about 0.2° in mean longitude at the arrival of the Cassini spacecraft in July 2004, and to about 3° in 2008, at the end of the nominal tour.     

Spin rate distribution of small asteroids

The spin rate distribution of main belt/Mars crossing (MB/MC) asteroids with diameters 3-15 km is uniform in the range from f=1 to 9.5 d⁻¹, and there is an excess of slow rotators with f<1 d⁻¹. The observed distribution appears to be controlled by the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect. The magnitude of the excess of slow rotators is related to the residence time of slowed down asteroids in the excess and the rate of spin rate change outside the excess. We estimated a median YORP spin rate change of ≈0.022 d⁻¹/Myr for asteroids in our sample (i.e., a median time in which the spin rate changes by 1 d⁻¹ is ≈45 Myr), thus the residence time of slowed down asteroids in the excess is ≈110 Myr. The spin rate distribution of near-Earth asteroids (NEAs) with sizes in the range 0.2-3 km (∼5 times smaller in median diameter than the MB/MC asteroids sample) shows a similar excess of slow rotators, but there is also a concentration of NEAs at fast spin rates with f=9-10 d⁻¹. The concentration at fast spin rates is correlated with a narrower distribution of spin rates of primaries of binary systems among NEAs; the difference may be due to the apparently more evolved population of binaries among MB/MC asteroids.     

The nucleus of comet C/1983 H1 IRAS-Araki-Alcock

We present a synthetic analysis of all available infrared (2-20 μm) and radio (1.3-6.1 cm) observations of comet C/1983 H1 IRAS-Araki-Alcock performed during its close approach to Earth in May 1983. We implement a model based on a spherical nucleus with a macroscopic mosaic of small and numerous active and inactive regions, and take into account the strong phase effect in the calculations of the thermal flux (often neglected in past interpretations). The orientation of the spin axis is assumed to be that determined by Sekanina [1988. Astron. J. 95, 1876-1894]. Additional constraints coming from visible photometry, measurements of the water production rate and the temporal variations of the cometary activity are introduced. We derive an equivalent nucleus radius of 3.4±0.5 km, consistent with a geometric albedo of 0.04 ±0.01 and a phase coefficient in the visible, and an active fraction of 2.9 ±1.9%. Although the nucleus is probably elongated as found in the past (Sekanina, 1988), we show that the relevant measurements were likely contaminated by the contribution of a variable coma.     

Simultaneous visible and near-infrared time resolved observations of the outer Solar System object (29981) 1999 TD10

The outer Solar System object (29981) 1999 TD₁₀ was observed simultaneously in the R, and J and H bands in September 2001, and in B, V, R, and I in October 2002. We derive B−V=0.80±0.05 mag, V−R=0.48±0.05 mag, R−I=0.44±0.05 mag, R−J=1.24±0.05 mag, and J−H=0.61±0.07 mag. Combining our data with the data from Rousselot et al. (2003, Astron. Astrophys. 407, 1139) we derive a synodic period of 15.382±0.001 hr in agreement with the period from Rousselot et al. Our observations at the same time, with better S/N and seeing, show no evidence of a coma, contrary to the claim by Choi et al. (2003, Icarus 165, 101).     

Astrometry and dynamics of Anthe (S/2007 S 4), a new satellite of Saturn

We describe the astrometry and dynamics of Anthe (S/2007 S 4), a new satellite of Saturn discovered in images obtained using the Imaging Science Subsystem (ISS) of the Cassini spacecraft. Included are details of 63 observations, of which 28 were obtained with Cassini's narrow-angle camera (NAC) and 35 using its wide-angle camera (WAC), covering an observation time-span of approximately 3 years. We estimate the diameter of Anthe to be ∼1.8 km. Orbit modeling based on a numerical integration of the full equations of motion fitted to the observations show that Anthe is in a first-order 11:10 mean motion resonance with Mimas. Two resonant arguments are librating: ?₁=11λ^′−10λ−? and ?₂=11λ^′−10λ−?^′−Ω^′+Ω, where λ, ? and Ω refer to the mean longitude, longitude of pericenter and longitude of ascending node of Mimas and Anthe, with the primed quantities corresponding to Anthe. These resonances cause periodic variations in the orbital elements. The semi-major axis varies by ±26 km over a 913-day period. Anthe is also close to a second-order eccentricity-type mean motion resonant relationship of the form 77:75 with Methone. Since Methone is also in a first-order resonance with Mimas [Spitale, J.N., Jacobson, R.A., Porco, C.C., Owen, W.M., 2006. Astron. J. 132, 692-710], an additional indirect perturbation exists between Methone and Anthe via Mimas. Neither effect is detectable in the orbit fitting and the short-term dynamical evolution of Anthe is dominated by the Mimas-Anthe resonances alone. The expected modulation effect from the Mimas-Tethys 4:2 inclination resonance is also insignificant over this time period. By including Cassini ISS observations of Mimas in the numerical integration fit, we estimate the GM of Mimas to be , consistent with Jacobson et al. [Jacobson, R.A., Spitale, J., Porco, C.C., Owen, W.M., 2006. Astron. J. 132, 711-713].     

Structure of self-gravity wakes in Saturn's A ring as measured by Cassini CIRS

The CIRS infrared spectrometer onboard the Cassini spacecraft has scanned Saturn's A ring azimuthally from several viewing angles since its orbit insertion in 2004. A quadrupolar asymmetry has been detected in this ring at spacecraft elevations ranging between 16° to 37°. Its fractional amplitude decreases from 22% to 8% from 20° to 37° elevations. The patterns observed in two almost complete azimuthal scans at elevations 20° and 36° strongly favor the self-gravity wakes as the origin of the asymmetry. The elliptical, infinite cylinder model of Hedman et al. [Hedman, M.M., Nicholson, P.D., Salo, H., Wallis, B.D., Buratti, B.J., Baines, K.H., Brown, R.H., Clark, R.N., 2007. Astron. J. 133, 2624-2629] can reproduce the CIRS observations well. Such wakes are found to have an average height-to-spacing ratio H/λ=0.1607±0.0002, a width-over-spacing W/λ=0.3833±0.0008. Gaps between wakes, which are filled with particles, have an optical depth τG=0.1231±0.0005. The wakes mean pitch angle ΦW is 70.70°±0.07°, relative to the radial direction. The comparison of ground-based visible data with CIRS observations constrains the A ring to be a monolayer. For a surface mass density of 40 g cm⁻² [Tiscarino, M.S., Burns, J.A., Nicholson, P.D., Hedman, M.M., Porco, C.C., 2007. Icarus 189, 14-34], the expected spacing of wakes is λ≈60 m. Their height and width would then be H≈10 m and W≈24 m, values that match the maximum size of particles in this ring as determined from ground-based stellar occultations [French, R.G., Nicholson, P.D., 2000. Icarus 145, 502-523].     

The proximity of Mercury's spin to Cassini state 1 from adiabatic invariance

In determining Mercury's core structure from its rotational properties, the value of the normalized moment of inertia, C/MR², from the location of Cassini 1 is crucial. If Mercury's spin axis occupies Cassini state 1, its position defines the location of the state, where the axis is fixed in the frame precessing with the orbit. Although tidal and core-mantle dissipation drive the spin to the Cassini state with a time scale O(¹⁰⁵) years, the spin might still be displaced from the Cassini state if the variations in the orbital elements induced by planetary perturbations, which change the position of the Cassini state, cause the spin to lag behind as it attempts to follow the state. After being brought to the state by dissipative processes, the spin axis is expected to follow the Cassini state for orbit variations with time scales long compared to the 1000 year precession period of the spin about the Cassini state because the solid angle swept out by the spin axis as it precesses is an adiabatic invariant. Short period variations in the orbital elements of small amplitude should cause displacements that are commensurate with the amplitudes of the short period terms. The exception would be if there are forcing terms in the perturbations that are nearly resonant with the 1000 year precession period. The precision of the radar and eventual spacecraft measurements of the position of Mercury's spin axis warrants a check on the likely proximity of the spin axis to the Cassini state. How confident should we be that the spin axis position defines the Cassini state sufficiently well for a precise determination of C/MR²? By following simultaneously the spin position and the Cassini state position during long time scale orbital variations over past 3 million years [Quinn, T.R., Tremaine, S., Duncan, M., 1991. Astron. J. 101, 2287-2305] and short time scale variations for 20,000 years [JPL Ephemeris DE 408; Standish, E.M., private communication, 2005], we show that the spin axis will remain within one arcsec of the Cassini state after it is brought there by dissipative torques. In this process the spin is located in the orbit frame of reference, which in turn is referenced to the inertial ecliptic plane of J2000. There are no perturbations with periods resonant with the precession period that could cause large separations. We thus expect Mercury's spin to occupy Cassini state 1 well within the uncertainties for both radar and spacecraft measurements, with correspondingly tight constraints on C/MR² and the extent of Mercury's molten core. Two unlikely caveats for this conclusion are: (1) an excitation of a free spin precession by an unknown mechanism or (2) a displacement by a dissipative core mantle interaction that exceeds the measurement uncertainties.     

Analyses of dust coma morphology of Comet Hyakutake (1996 B2) near perigee: outburst behavior, jet motion, source region locations, and nucleus pole orientation

We present inner-coma dust imaging of Comet Hyakutake (1996 B2) obtained on 11 consecutive nights in late March 1996, an interval including a major outburst and the comet’s closest approach to Earth. The evolution of the outburst morphology is followed, along with the motion along the tail of several outburst fragments. Two spiral dust jets—a primary jet, along with a much weaker secondary jet—are visible throughout the interval and are produced by two source regions on a rotating nucleus. These are examined as a function of rotational phase and viewing geometry, with their appearance changing from a nearly face-on view on March 18 to side-on by March 28. The dust outflow velocity as a function of distance from the nucleus is derived, with the dust continuing to accelerate to a distance of 4000 km or more and reaching an average outflow velocity of 0.38 km s⁻¹ between 3000 and 8000 km. We present details of our Monte Carlo modeling of the jets and our methodology of fitting the model to the images. The modeling yields the pole orientation of the nucleus, with an obliquity of approximately 108°, corresponding to an RA of 13^h41^m and a Dec of −1.1°. For an assumed spherical nucleus, the primary active region is centered at approximately −66° latitude, has a radius of about 56°, and therefore covers about 22% of the surface. The source of the secondary jet is at a latitude of −28°, has a radius of about 16°, and is located at a longitude nearly 180° away from the primary source. Estimated uncertainties for the pole orientation and the source locations and sizes are each about 3°. This solution for the nucleus orientation and source locations explains the strong asymmetry in measured production rates before and after perihelion in radio observations (Biver et al., 1999, Astron. J. 118, 1850-1872). The modeling also tightly constrains the sidereal rotation period as 0.2618 ± 0.0001 day, completely consistent with the expected +0.0003 day difference from the observed solar rotation period of 0.2614 ± 0.0004 day determined by Schleicher and Osip (2002, Icarus 159, 210-233), given the pole orientation and position of the comet in its orbit.     

Formation of Jupiter using opacities based on detailed grain physics

Numerical simulations, based on the core-nucleated accretion model, are presented for the formation of Jupiter at 5.2 AU in three primordial disks with three different assumed values of the surface density of solid particles. The grain opacities in the envelope of the protoplanet are computed using a detailed model that includes settling and coagulation of grains and that incorporates a recalculation of the grain size distribution at each point in time and space. We generally find lower opacities than the 2% of interstellar values used in previous calculations (Hubickyj, O., Bodenheimer, P., Lissauer, J.J. [2005]. Icarus 179, 415-431; Lissauer, J.J., Hubickyj, O., D’Angelo, G., Bodenheimer, P. [2009]. Icarus 199, 338-350). These lower opacities result in more rapid heat loss from and more rapid contraction of the protoplanetary envelope. For a given surface density of solids, the new calculations result in a substantial speedup in formation time as compared with those previous calculations. Formation times are calculated to be 1.0, 1.9, and 4.0 Myr, and solid core masses are found to be 16.8, 8.9, and 4.7 M_⊕, for solid surface densities, σ, of 10, 6, and 4 g cm⁻², respectively. For σ = 10 and σ = 6 g cm⁻², respectively, these formation times are reduced by more than 50% and more than 80% compared with those in a previously published calculation with the old approximation to the opacity.     

HCN and CN in Comet 2P/Encke: Models of the non-isotropic, rotation-modulated coma and CN parent life time

Axisymmetric models of the outgassing of a cometary nucleus have been constructed. Such models can be used to describe a nucleus with a single active region. The models may include a solar zenith angle dependence of the outgassing. They retrieve the outgassing flux at distances from the nucleus where collisions between molecules are unimportant, as function of the angle with respect to the outgassing axis. The observed emissions must be optically thin. Furthermore the models assume that the outflow speed at large distance from the nucleus does not depend on direction. The value of the outflow speed is retrieved. The models are applied to CN images and HCN spectra of Comet 2P/Encke, obtained nearly simultaneously in November 2003 with the 2 m optical telescope on Mount Rozhen, Bulgaria, and with the 10 m Heinrich Hertz Submillimeter Telescope on Mount Graham, Arizona, USA. According to Sekanina (1988), Astron. J. 95, 911-924, at that time a single outgassing source was active. Input parameters to the models like the rotation period of the nucleus and a small correction to Sekanina’s rotation axis are determined from a simpler jet position angle model. The rotation is prograde with a sideric period of 11.056 ± 0.024 h, in agreement with literature values. The best fit model has an outflow speed of 0.95 ± 0.04 km s⁻¹. The same value has been derived from the corkscrew appearing in the CN images. The location of the outgassing axis is at colatitude δ_a = 7.4° ± 2.9° and longitude λ_a = 235° ± 17° (a definition of zero longitude is provided). Comet Encke’s outgassing corresponds approximately to the longitudinally averaged solar input on a spherical nucleus (i.e. very likely comes from deeper layers) but with some deficiency of outgassing at mid-latitudes and non-zero outgassing from the dark polar cap. The presence of gas flow from the dark polar cap is explained as evidence of gas flow across the terminator. The models rely mostly on the CN images. The HCN spectra are more noisy. They provide information how to determine the best fit outflow velocity and the sense of rotation. The model HCN spectra are distinctly non-Gaussian. Within error limits they are consistent with the observations. Models based solely on the HCN spectra are also presented but, because of the lower quality of the data and the unfavorable observing geometry, yield inferior results. As a by-product we determine the CN parent life time from our CN observations. The solar EUV and Lyα radiation field at the time of our observations is taken into account.     

2006 Fragmentation of Comet 73P/Schwassmann-Wachmann 3B observed with Subaru/Suprime-Cam

The fragmentation of the split Comet 73P/Schwassmann-Wachmann 3 B was observed with the prime-focus camera Suprime-Cam attached to the Subaru 8.2-m telescope. The fragmentation revealed dozens of miniature comets [Fuse, T., Yamamoto, N., Kinoshita, D., Furusawa, H., Watanabe, J., 2007. Publ. Astron. Soc. Jpn. 59 (2), 381-386]. We analyzed the Subaru/Suprime-Cam images, detecting no fewer than 154 mini-comets, mostly extending to the southwest. Three were close to the projected orbit of fragment B. We applied synchrone-syndyne analysis, modified for rocket effect analysis, to the mini-fragment spatial distribution. We found that most of these mini-comets were ejected from fragment B by an outburst occurring around 1 April 2006, and three fragments on the leading side of nucleus B could have been released sunward on the previous return. Several fragments might have been released by successive outbursts around 24 April and 2 May 2006. The ratio of the rocket force to solar gravity was 7-23 times larger than that exerted on fragment B. No significant color variation was found. The mean color index, V-R = 0.50 ± 0.07, was slightly redder than that of the Sun and similar to that of the largest fragment, C, which suggests that these mini-fragments were detected mainly through sunlight reflected by dust particles and materials on the nuclei. We examined the surface brightness profiles of all detected fragments and estimated the sizes of 154 fragments. We found that the radius of these mini-fragments was in the 5- to 108-m range (equivalent size of Tunguska impactor). The power-law index of the differential size distribution was q = −3.34 ± 0.05. Based on this size distribution, we found that about 1-10% of the mass of fragment B was lost in the April 2006 outbursts. Modeling the cometary fragment dynamics [Desvoivres, E., Klinger, J., Levasseur-Regourd, A.C., Lecacheux, J., Jorda, L., Enzian, A., Colas, F., Frappa, E., Laques, P., 1999. Mon. Not. Roy. Astron. Soc. 303 (4), 826-834; Desvoivres, E., Klinger, J., Levasseur-Regourd, A.C., Jones, G.H., 2000. Icarus 144, 172-181] revealed that it is likely that mini-fragments smaller than ∼10-20 m could be depleted in water ice and become inactive, implying that decameter-sized comet fragments could survive against melting and remain as near-Earth objects. We attempted to detect the dust trail, which was clearly found in infrared wavelengths by Spitzer. No brightness enhancement brighter than 30.0 mag arcsec⁻² (3σ) was detected in the orbit of fragment B.     

Gravitational scattering within the Himalia group of jovian prograde irregular satellites

We test the hypothesis that gravitational scattering within the Himalia prograde irregular satellite group is responsible for the large velocity dispersion identified by Nesvorný et al. (2003, Astron. J. 126, 398-429). We carry out numerical simulations of the dynamical evolution of this group under several scenarios for Himalia's mass for ¹⁰⁸ yr. We find that modification of the satellite orbits is significant, 10-15% of their semimajor axes differences but <5% of the eccentricities. Additionally, the inclination of the satellite Lysithea may have been modified by a secular resonance with Himalia. We scale the observed orbital element diffusion to the age of the Solar System using a power law approximation. The projected changes in a and e can reduce the observable ejection speeds of Elara and Lysithea with respect to Himalia under the 100 m s⁻¹ mark, rendering them consistent with hydrocode simulation results (Michel et al., 2002, Icarus 160, 10-23). The dispersion magnitude required to migrate Elara to its present orbit implies that Himalia's size has been underestimated by Cassini (Porco et al., 2003, Science 299, 1541-1547) and/or that its density is significantly higher than that of extensively fractured or “rubble-pile” bodies. In addition, Lysithea could have been created as recently as the last 500 Myr of the Solar System's history. Simple calculations indicate that such diffusion would also be evident in groups of satellites genetically related to Phoebe at Saturn, Caliban and Sycorax at Uranus and Nereid at Neptune.     

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.     

Multi-wavelength observations of Asteroid 2100 Ra-Shalom

We observed near-Earth asteroid (NEA) 2100 Ra-Shalom over a six-year period, obtaining rotationally resolved spectra in the visible, near-infrared, thermal-infrared, and radar wavelengths. We find that Ra-Shalom has an effective diameter of Deff=2.3±0.2 km, rotation period P=19.793±0.001 h, visual albedo pv=0.13±0.03, radar albedo , and polarization ratio μc=0.25±0.04. We used our radar observations to generate a three-dimensional shape model which shows several structural features of interest. Based on our thermal observations, Ra-Shalom has a high thermal inertia of ∼10³ J m⁻² s^−0.5 K⁻¹, consistent with a coarse or rocky surface and the inferences of others [Harris, A.W., Davies, J.K., Green, S.F., 1998. Icarus 135, 441-450; Delbo, M., Harris, A.W., Binzel, R.P., Pravec, P., Davies, J.K., 2003. Icarus 166, 116-130]. Our spectral data indicate that Ra-Shalom is a K-class asteroid and we find excellent agreement between our spectra and laboratory spectra of the CV3 meteorite Grosnaja. Our spectra show rotation-dependent variations consistent with global variations in grain size. Our radar observations show rotation-dependent variations in radar albedo consistent with global variations in the thickness of a relatively thin regolith.     

The Existence of a Planet beyond 50 AU and the Orbital Distribution of the Classical Edgeworth-Kuiper-Belt Objects

We study the effects of a Mars-like planetoid with a semimajor axis at about ∼60 AU orbiting embedded in the primordial Edgeworth-Kuiper belt (EKB). The origin of such an object can be explained in the framework of our current understanding of the origin of the outer Solar System, and a scenario for the orbital transport mechanism to its present location is given. The existence of such an object would produce a gap in the EKB distribution with an edge at about 50 AU, which seems to be in agreement with the most recent observations. No object at low eccentricity with semimajor axis beyond 50 AU has been detected so far, even though the present observing capabilities would allow an eventual detection (B. Gladman et al. 1998, Astron. J.116, 2042-2054; D. Jewitt et al. 1998, Astron. J.115, 2125-2135; E. I. Chiang and M. E. Brown 1999, Astron. J.118, 1411-1422; R. L. Allen et al. 2000, Astrophys. J.549, 241-244; C. A. Trujillo et al. 2001, Astron. J.122, 457-473; B. Gladman et al. 2001, Astron. J.122, 1051-1066; C. A. Trujillo and M. E. Brown 2001, Astrophys. J.554, 95-98). Finally, ranges for the magnitude and proper motion of the proposed object are given.     

Moonlet induced wakes in planetary rings: Analytical model including eccentric orbits of moon and ring particles

Saturn’s rings host two known moons, Pan and Daphnis, which are massive enough to clear circumferential gaps in the ring around their orbits. Both moons create wake patterns at the gap edges by gravitational deflection of the ring material (Cuzzi, J.N., Scargle, J.D. [1985]. Astrophys. J. 292, 276-290; Showalter, M.R., Cuzzi, J.N., Marouf, E.A., Esposito, L.W. [1986]. Icarus 66, 297-323). New Cassini observations revealed that these wavy edges deviate from the sinusoidal waveform, which one would expect from a theory that assumes a circular orbit of the perturbing moon and neglects particle interactions. Resonant perturbations of the edges by moons outside the ring system, as well as an eccentric orbit of the embedded moon, may partly explain this behavior (Porco, C.C., and 34 colleagues [2005]. Science 307, 1226-1236; Tiscareno, M.S., Burns, J.A., Hedman, M.M., Spitale, J.N., Porco, C.C., Murray, C.D., and the Cassini Imaging team [2005]. Bull. Am. Astron. Soc. 37, 767; Weiss, J.W., Porco, C.C., Tiscareno, M.S., Burns, J.A., Dones, L. [2005]. Bull. Am. Astron. Soc. 37, 767; Weiss, J.W., Porco, C.C., Tiscareno, M.S. [2009]. Astron. J. 138, 272-286). Here we present an extended non-collisional streamline model which accounts for both effects. We describe the resulting variations of the density structure and the modification of the nonlinearity parameter q. Furthermore, an estimate is given for the applicability of the model. We use the streamwire model introduced by Stewart (Stewart, G.R. [1991]. Icarus 94, 436-450) to plot the perturbed ring density at the gap edges.We apply our model to the Keeler gap edges undulated by Daphnis and to a faint ringlet in the Encke gap close to the orbit of Pan. The modulations of the latter ringlet, induced by the perturbations of Pan (Burns, J.A., Hedman, M.M., Tiscareno, M.S., Nicholson, P.D., Streetman, B.J., Colwell, J.E., Showalter, M.R., Murray, C.D., Cuzzi, J.N., Porco, C.C., and the Cassini ISS team [2005]. Bull. Am. Astron. Soc. 37, 766), can be well described by our analytical model. Our analysis yields a Hill radius of Pan of 17.5 km, which is 9% smaller than the value presented by Porco (Porco, C.C., and 34 colleagues [2005]. Science 307, 1226-1236), but fits well to the radial semi-axis of Pan of 17.4 km. This supports the idea that Pan has filled its Hill sphere with accreted material (Porco, C.C., Thomas, P.C., Weiss, J.W., Richardson, D.C. [2007]. Science 318, 1602-1607). A numerical solution of a streamline is used to estimate the parameters of the Daphnis-Keeler gap system, since the close proximity of the gap edge to the moon induces strong perturbations, not allowing an application of the analytic streamline model. We obtain a Hill radius of 5.1 km for Daphnis, an inner edge variation of 8 km, and an eccentricity for Daphnis of 1.5 × 10⁻⁵. The latter two quantities deviate by a factor of two from values gained by direct observations (Jacobson, R.A., Spitale, J., Porco, C.C., Beurle, K., Cooper, N.J., Evans, M.W., Murray, C.D. [2008]. Astron. J. 135, 261-263; Tiscareno, M.S., Burns, J.A., Hedman, M.M., Spitale, J.N., Porco, C.C., Murray, C.D., and the Cassini Imaging team [2005]. Bull. Am. Astron. Soc. 37, 767), which might be attributed to the neglect of particle interactions and vertical motion in our model.