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.     

The dust trail of Comet 67P/Churyumov-Gerasimenko

We report the detection of Comet 67P/Churyumov-Gerasimenko's dust trail and nucleus in 24 μm Spitzer Space Telescope images taken February 2004. The dust trail is not found in optical Palomar images taken June 2003. Both the optical and infrared images show a distinct neck-line tail structure, offset from the projected orbit of the comet. We compare our observations to simulated images using a Monte Carlo approach and a dynamical model for comet dust. We estimate the trail to be at least one orbit old (6.6 years) and consist of particles of size ?100 μm. The neck-line is composed of similar sized particles, but younger in age. Together, our observations and simulations suggest grains 100 μm and larger in size dominate the total mass ejected from the comet. The radiometric effective radius of the nucleus is 1.87±0.08 km, derived from the Spitzer observation. The Rosetta spacecraft is expected to arrive at and orbit this comet in 2014. Assuming the trail is comprised solely of 1 mm radius grains, we compute a low probability (∼¹⁰⁻³) of a trail grain impacting with Rosetta during approach and orbit insertion.     

Meteor outbursts from long-period comet dust trails

Long-period comets have narrow one-revolution old dust trails that can cause meteor outbursts when encountered by Earth. To facilitate observing campaigns that will characterize and perhaps help find Earth-threatening, long-period comets from their trace of meteoric debris, we use past accounts of outbursts from 14 different showers to calculate the future dust trail positions near Earth’s orbit. We also examine known near-Earth, long-period comets and identify five potential new showers, which can be utilized to learn more about these objects. We demonstrate that it is the one-revolution trail that is responsible for meteor outbursts. A method that calculates in what year these showers are likely to return and at what hour is presented. The calculations improve on earlier approximate methods that used the Sun’s reflex motion to gauge the trail motion relative to Earth’s orbit.     

Solar flare track densities in interplanetary dust particles: The determination of an asteroidal versus cometary source of the zodiacal dust cloud

Dust particles that are larger than 1 μm, when injected into the Solar System from comets and asteroids, will spiral into the Sun due to the Poynting-Robertson effect. During the process of spiraling in, such dust particles accumulate solar flare tracks in their component minerals. The accumulated track density for a given dust grain is a function of the duration of its space exposure and its distance from the Sun. Using a computer model, it was determined that the expected track density distributions from grains produced by comets are very different from those produced by asteroids. Individual asteroids produce populations of particles that arrive at 1 AU with scaled track density distributions containing “spikes,” while comets supply particles with a flatter and wider distribution of track densities. Particles with track densities above 3 × 10⁷ (sϱA/v) tracks/cm² have probably been exposed to solar flare tracks prior to injection into the interplanetary medium and are therefore likely to be asteroidal. Particles with track densities below 0.7 × 10⁷(sϱA/v) tracks/cm² must be derived from comets or Earth-crossing asteroids. Earth-crossing asteroids are not responsible for all the dust collected at 1 AU since they cannot produce the large track densities observed in some of the interplanetary dust particles collected in the stratosphere. The track densities observed in the stratospheric dust fall within the predicted range, but there is at present an insufficient number of carefully determined densities to make strong statements about the sources of the present dust population.     

Physical properties of asteroid dust bands and their sources

Disruptive collisions in the main belt can liberate fragments from parent bodies ranging in size from several micrometers to tens of kilometers in diameter. These debris bodies group at initially similar orbital locations. Most asteroid-sized fragments remain at these locations and are presently observed as asteroid families. Small debris particles are quickly removed by Poynting-Robertson drag or comminution but their populations are replenished in the source locations by collisional cascade. Observations from the Infrared Astronomical Satellite (IRAS) showed that particles from particular families have thermal radiation signatures that appear as band pairs of infrared emission at roughly constant latitudes both above and below the Solar System plane. Here we apply a new physical model capable of linking the IRAS dust bands to families with characteristic inclinations. We use our results to constrain the physical properties of IRAS dust bands and their source families. Our results indicate that two prominent IRAS bands at inclinations ≈2.1° and ≈9.3° are byproducts of recent asteroid disruption events. The former is associated with a disruption of a ≈30-km asteroid occurring 5.8 Myr ago; this event gave birth to the Karin family. The latter came from the breakup of a large >100-km-diameter asteroid 8.3 Myr ago that produced the Veritas family. Using an N-body code, we tracked the dynamical evolution of ≈10⁶ particles, 1 μm to 1 cm in diameter, from both families. We then used these results in a Monte Carlo code to determine how small particles from each population undergo collisional evolution. By computing the thermal emission of particles, we were able to compare our results with IRAS observations. Our best-fit model results suggest the Karin and Veritas family particles contribute by 5-9% in 10-60-μm wavelengths to the zodiacal cloud's brightness within 50° latitudes around the ecliptic, and by 9-15% within 10° latitudes. The high brightness of the zodiacal cloud at large latitudes suggests that it is mainly produced by particles with higher inclinations than what would be expected for asteroidal particles produced by sources in the main belt. From these results, we infer that asteroidal dust represents a smaller fraction of the zodiacal cloud than previously thought. We estimate that the total mass accreted by the Earth in Karin and Veritas particles with diameters 20-400 μm is ≈15,000-20,000 tons per year (assuming 2 g cm⁻³ particles density). This is ≈30-50% of the terrestrial accretion rate of cosmic material measured by the Long Duration Exposure Facility. We hypothesize that up to ≈50% of our collected interplanetary dust particles and micrometeorites may be made up of particle species from the Veritas and Karin families. The Karin family IDPs should be about as abundant as Veritas family IDPs though this ratio may change if the contribution of third, near-ecliptic source is significant. Other sources of dust and/or large impact speeds must be invoked to explain the remaining ≈50-70%. The disproportional contribution of Karin/Veritas particles to the zodiacal cloud (only 5-9%) and to the terrestrial accretion rate (30-50%) suggests that the effects of gravitational focusing by the Earth enhance the accretion rate of Karin/Veritas particles relative to those in the background zodiacal cloud. From this result and from the latitudinal brightness of the zodiacal cloud, we infer that the zodiacal cloud emission may be dominated by high-speed cometary particles, while the terrestrial impactor flux contains a major contribution from asteroidal sources. Collisions and Poynting-Robertson drift produce the size-frequency distribution (SFD) of Karin and Veritas particles that becomes increasingly steeper closer to the Sun. At 1 AU, the SFD is relatively shallow for small particle diameters D (differential slope exponent of particles with D?100 μm is ≈2.2-2.5) and steep for D?100 μm. Most of the mass at 1 AU, as well as most of the cross-sectional area, is contributed by particles with D≈100-200 μm. Similar result has been found previously for the SFD of the zodiacal cloud particles at 1 AU. The fact that the SFD of Karin/Veritas particles is similar to that of the zodiacal cloud suggests that similar processes shaped these particle populations. We estimate that there are ≈5×¹⁰²⁴ Karin and ≈10²⁵ Veritas family particles with D>30 μm in the Solar System today. The IRAS observation of the dust bands may be satisfactorily modeled using ‘averaged’ SFDs that are constant with semimajor axis. These SFDs are best described by a broken power-law function with differential power index α≈2.1-2.4 for D?100 μm and by α?3.5 for 100 μm?D?1 cm. The total cross-sectional surface area of Veritas particles is a factor of ≈2 larger than the surface area of the particles producing the inner dust bands. The total volumes in Karin and Veritas family particles with 1 μm<D<1 cm correspond to D=11 km and D=14 km asteroids with equivalent masses ≈1.5×¹⁰¹⁸ g and ≈3.0×¹⁰¹⁸ g, respectively (assuming 2 g cm⁻³ bulk density). If the size-frequency and radial distribution of particles in the zodiacal cloud were similar to those in the asteroid dust bands, we estimate that the zodiacal cloud represents ∼3×¹⁰¹⁹ g of material (in particles with 1 μm<D<1 cm) at ±10° around the ecliptic and perhaps as much as ∼¹⁰²⁰ g in total. The later number corresponds to about a 23-km-radius sphere with 2 g cm⁻³ density.     

Ensemble properties of comets in the Sloan Digital Sky Survey

We present the ensemble properties of 31 comets (27 resolved and 4 unresolved) observed by the Sloan Digital Sky Survey (SDSS). This sample of comets represents about 1 comet per 10 million SDSS photometric objects. Five-band (u, g, r, i, z) photometry is used to determine the comets’ colors, sizes, surface brightness profiles, and rates of dust production in terms of the Afρ formalism. We find that the cumulative luminosity function for the Jupiter Family Comets in our sample is well fit by a power law of the form N(<H) ∝ 10^{(0.49±0.05)H} for H < 18, with evidence of a much shallower fit N(<H) ∝ 10^{(0.19±0.03)H} for the faint (14.5 < H < 18) comets. The resolved comets show an extremely narrow distribution of colors (0.57 ± 0.05 in g ? r for example), which are statistically indistinguishable from that of the Jupiter Trojans. Further, there is no evidence of correlation between color and physical, dynamical, or observational parameters for the observed comets.     

Optical properties of dust from Jupiter family comets

To try to define specific physical properties of the dust of Jupiter-family comets (JFCs), we compare the light scattered by them. Amongst the more than 1000 JFCs, less than 200 are numbered, 40 of them being rather bright. In the present work we use data from the latter. In situ observations of three nuclei show low albedo surfaces. The albedo of the dust particles in the coma is low, with generally a red colour. The A(α)fρ product is a measure of cometary activity and secular changes. Images of different regions (jets and fans) give indications on the nucleus rotation and position of the emitting areas, as compared to the position of the rotation axis. Differences in physical properties between the particles in different regions are pointed out by differences in the linear polarization of the scattered light and by spectral variations in brightness and polarization. Jupiter family comets are considered as dust-poor comets. Tails and trails’ studies give an estimation of the size distribution of the particles. However the dust production rates depend on the largest particles (up to centimetre size), which are mainly observed in the trails where large dark compact particles are found. These dark particles are also responsible for the high polarization in the inner most coma of some comets. The meaning, in terms of physical properties, of the linear polarization is discussed through different examples such as 2P/Encke, 9P/Tempel 1 or the fragments of 73P/Schwassmann-Wachmann 3. Cometary outbursts and splitting events show that the properties of the dust ejected from the interior of the nucleus are similar to the ones of more active comets (new or with larger semi-major axis).     

Oxygen isotopic composition of an enstatite ribbon of probable cometary origin

Filamentary enstatite crystals are found in interplanetary dust particles (IDPs) of likely cometary origin but are very rare or absent in meteorites. Crystallographic characteristics of filamentary enstatites indicate that they condensed directly from vapor. We measured the O isotopic composition of an enstatite ribbon from a giant cluster IDP to be δ¹⁸O = 25 ± 55, δ¹⁷O = ?19 ± 129, ?¹⁷O = ?32 ± 134 (2σ errors), which is inconsistent at the 2σ level with the composition of the Sun inferred from the Genesis solar wind measurements. The particle's O isotopic composition, consistent with the terrestrial composition, implies that it condensed from a gas of nonsolar O isotopic composition, possibly as a result of vaporization of disk region enriched in ¹⁶O‐depleted solids. The relative scarcity of filamentary enstatite in asteroids compared to comets implies either that this crystal condensed from dust vaporized in situ in the outer solar system where comets formed or it condensed in the inner solar system and was subsequently transported outward to the comet‐forming region.     

Collisional balance of the meteoritic complex

Taking into account meteoroid measurements by in situ experiments, zodiacal light observations, and oblique angle hypervelocity impact studies, it is found that the observed size distributions of lunar microcraters usually do not represent the interplanetary meteoroid flux for particles with masses ?10^?10g. From the steepest observed lunar crater size distribution a “lunar flux” is derived which is up to 2 orders of magnitude higher than the interplanetary flux at the smallest particle masses. New models of the “lunar” and “interplanetary” meteoroid fluxes are presented. The spatial mass density of interplanetary meteoritic material at 1 AU is ～10^?16g/m³. A large fraction of this mass is in particles of 10^?6 to 10^?4 g. A detailed analysis of the effects of mutual collisions (i.e., destruction of meteoroids and production of fragment particles) and of radiation pressure has been performed which yielded a new picture of the balance of the meteoritic complex. It has been found that the collisional lifetime at 1 AU is shortest (～10⁴years) for meteoroids of 10^?4 to 1 g mass. For particles with masses m > 10^?5g, Poynting-Robertson lifetimes are considerably larger than collisional lifetimes. The collisional destruction rate of meteoroids with masses $\text{m ? 10}{}^{\mathrm{?3}}\text{g}$ is about 10 times larger than the rate of collisional production of fragment particles in the same mass range. About 9 tons/sec of these “meteor-sized” (m > 10^?5g) particles are lost inside 1 AU due to collisions and have to be replenished by other sources, e.g., comets. Under steady-state conditions, most of these large particles are “young”; i.e., they have not been fragmented by collisions and their initial orbits are not altered much by radiation pressure drag. Many more micrometeoroids of masses m ? 10^?5g are generated by collisions from more massive particles than are destroyed by collisions. The net collisional production rate of intermediate-sized particles 10^?10g ? m ? 10^?5g is found to be about 16 times larger at 1 AU than the Poynting-Robertson loss rate. The total Poynting-Robertson loss rate inside 1 AU is only about 0.26 tons/sec. The smallest fragment particles (m ? 10^?10g) will be largely injected into hyperbolic trajectories under the influence of radiation pressure (β meteoroids). These particles provide the most effecient loss mechanism from the meteoritic complex. When it is assumed that meteoroids fragment similarly to experimental impact studies with basalt, then it is found that interplanetary meteoroids in the mass range 10^?10g ? m ? 10^?5g cannot be in temporal balance under collisions and Poynting-Robertson drag but their spatial density is presently increasing with time.     

Properties of the dust cloud caused by the Deep Impact experiment

We present an analysis of the observations of the Deep Impact event performed by the OSIRIS narrow angle camera aboard the Rosetta spacecraft over two weeks, in an effort to characterize the cometary dust grains ejected from the nucleus of Comet 9P/Tempel 1. We adopt a Monte Carlo approach to generate calibrated synthetic images, and a linear combination of them is fitted to the calibrated images so as to determine the physical parameters of the dust cloud. Our model considers spherical olivine particles with a density of 3780 kg m⁻³. It incorporates constraints on the direction of the cone of emission coming from additional images obtained at Pic du Midi observatory, and constraints on the dust terminal velocities coming from the physics of the impact. We find that the slope of the differential dust size distribution of grains with radii <20 μm (β>0.008) is 3.1±0.3, a value typical of cometary dust tails. This shows that there is no evidence in our data for an enhancement in sub-micron particles in the ejecta compared to the typical dust distribution of active comets. We estimate the mass of particles with radii <1.4 μm (β>0.14) to be 1.5±0.2×¹⁰⁵ kg. These particles represent more than 80% of the cross-section of the observed dust cloud. The mass carried by larger particles depends whether the gas significantly increases the kinetic energy of the grains in the inner coma; it lies in the range 1-14×¹⁰⁶ kg for particles with radii <100 μm (β>0.002). We obtain the distribution of terminal velocities reached by the dust after the dust-gas interaction which is very well constrained between 10 and 600 m s⁻¹. It is characterized by Gaussian with a maximum at about 190 m s⁻¹ and a width at half maximum of 150 m s⁻¹.     

A tale of two very different comets: ISO and MSX measurements of dust emission from 126P/IRAS (1996) and 2P/Encke (1997)

We present the characteristics of the dust comae of two comets, 126P/IRAS, a member of the Halley family (a near-isotropic comet), and 2P/Encke, an ecliptic comet. We have primarily used mid- and far-infrared data obtained by the ISOPHOT instrument aboard the Infrared Space Observatory (ISO) in 1996 and 1997, and mid-infrared data obtained by the SPIRIT III instrument aboard the Midcourse Space Experiment (MSX) in 1996. We find that the dust grains emitted by the two comets have markedly different thermal and physical properties. P/IRAS's dust grain size distribution appears to be similar to that of fellow family member 1P/Halley, with grains smaller than 5 microns dominating by surface area, whereas P/Encke emits a much higher fraction of big (20 μm and higher) grains, with the grain mass distribution being similar to that which is inferred for the interplanetary dust population. P/Encke's dearth of micron-scale grains accounts for its visible-wavelength classification as a “gassy” comet. These conclusions are based on analyses of both imaging and spectrophotometry of the two comets; this combination provides a powerful way to constrain cometary dust properties. Specifically, P/IRAS was observed preperihelion while 1.71 AU from the Sun, and seen to have a 15-arcmin long mid-infrared dust tail pointing in the antisolar direction. No sunward spike was seen despite the vantage point being nearly in the comet's orbital plane. The tail's total mass at the time was about 8×10⁹ kg. The spectral energy distribution (SED) is best fit by a modified greybody with temperature T=265±15 K and emissivity ε proportional to a steep power law in wavelength λ: ε∝λ^−α, where α=0.50±0.20(2σ). This temperature is elevated with respect to the expected equilibrium temperature for this heliocentric distance. The dust mass loss rate was between 150-600 kg/s (95% confidence), the dust-to-gas mass loss ratio was about 3.3, and the albedo of the dust was 0.15±0.03. Carbonaceous material is depleted in the comet's dust by a factor of 2-3, paralleling the C₂ depletion in P/IRAS's gas coma. P/Encke, on the other hand, observed while 1.17 AU from the Sun, had an SED that is best fit by a Planck function with T=270±15 K and no emissivity falloff. The dust mass loss rate was 70-280 kg/s (95% confidence), the dust-to-gas mass loss ratio was about 2.3, and the albedo of the dust was about 0.06±0.02. These conclusions are consistent with the strongly curved dust tail and bright dust trail seen by Reach et al. (2000; Icarus 148, 80) in their ISO 12-μm imaging of P/Encke. The observed differences in the P/IRAS and P/Encke dust are most likely due to the less evolved and insolated state of the P/IRAS nuclear surface. If the dust emission behavior of P/Encke is typical of other ecliptic comets, then comets are the major supplier of the interplanetary dust cloud.     

Comparison of X-rays from Comets LINEAR (C/1999 S4) and McNaught-Hartley (C/1999 T1)

The Chandra X-ray Observatory (CXO) observations of Comets McNaught-Hartley (MH) and LINEAR S4 (S4) have been processed in the same way to compare X-rays from those comets. The X-ray isophotes are crescent-like in S4 and more circular in MH because of the different phase angles (98° and 44°, respectively). The peak X-ray brightness is greater in S4 than that in MH by a factor of 1.5 and smaller by a factor of 1.7 after the correction for heliocentric distance. The X-ray luminosities of MH and S4 are equal to 8.6 and 1.4×10¹⁵ erg s⁻¹ inside the apertures of ρ=1.5 and 0.5×10⁴ km, respectively. (Brightness is 20% of the peak value at these ρ.) Efficiencies of X-ray excitation corrected to the solar wind flow are similar and equal to 4.3×10⁻¹⁴ erg AU^3/2 in both comets. This confirms the solar wind excitation of X-rays in comets. Spectra of the comets were extracted with a special care of the background correction and using an energy-dependent spectral resolution code. The MH spectrum consists of ten emissions instead of nine emissions in the previously published spectrum. The new emission at 307 eV fills in a strong minimum in the previous spectrum and removes the major difference between that spectrum and the synthetic spectrum. This emission is assigned to the C⁺⁴ and Mg⁺⁹ lines. The positions of the other emissions and their identification are similar to those in the previous spectrum. The S4 spectrum consists of eight emissions, and four emissions are the same as in MH. The line identification is given. Ion ratios in the solar wind have been extracted from the spectra. O⁺⁸/O⁺⁷ is equal to 0.29±0.04 and 0.14±0.02 in MH and S4, and this difference correlates with the higher solar wind speed in S4. Ne⁺⁹/O⁺⁷ is (15±6)×10⁻³ and (19±7)×10⁻³, and these are the first data on Ne⁺⁹ in the solar wind. C⁺⁶/O⁺⁷ is 0.7±0.2 in both MH and S4. X-ray spectroscopy of comets may be used as a diagnostic tool to study the solar wind composition.     

Very small bodies in the solar system: The Earth's atmosphere as detector

Particles of mass less than about 1 gm are a minor fraction of the total matter impinging on the Earth averaged over millennia time scales. However, these particles dominate during a single particular year and produce the most obvious evidence of incoming extra-terrestrial matter in the form of ablation trails in the atmosphere which are visible at night as meteors.Observations of meteors give astronomical information on the composition, structure, and cometary associations of the particles. The composition is deduced from optical spectra of meteors, whilst telescopic studies of the trails during formation give information on the physical structure of the particles. Any cometary associations are deduced from measurement of meteor orbits determined photographically, using television, or by radar.Meteors occur in the atmosphere at heights from about 70 to 120 km. Optical observations are restricted to night-time and usually under conditions of low moonlight. A typical television based detector can record +8^M meteors with a sporadic rate of 15–20 per hour and velocities accurate to about 3%. The luminosity of the trail is strongly dependent on the velocity of the meteoroid (to about the third power).Radar observations of meteors are unrestricted by weather or time of day, and can readily detect meteors at least two orders of magnitude smaller in mass than those detectable optically. Again the observations are heavily biased toward the higher velocities as the electron line density varies approximately asV ^3.5. However, the higher the velocity of the meteoroid the greater the height of the meteor trail, and the reduced probability of radar detection due to rapid diffusion of the trail. Thus radar observations tend to select meteors in the intermediate velocity range 30–40 km s^–1.     

CN, C2 radicals and dust in Comets Shoemaker-Levy 1991 T2 and P/deVico 1995 S1

The results of the multiaperture photometry of Comet Shoemaker-Levy 1991 T2 in the pre-perihelion and P/deVico in the post-perihelion period with the narrowband CN, C₂ and Blue Continuum (BC) IHW filters are presented. A Haser model of the molecular coma was used for the determination of the parent and daughter scale-lengths and production rates of the radicals. The comets showed some substantial differences between their parent scale-lengths. The CN parent scale-length (at 1.0 AU) was 16×¹⁰³ km for Comet Shoemaker-Levy and 39×¹⁰³ for P/deVico, the C₂ parent scale-lengths were respectively 29×¹⁰³ and 54×¹⁰³ km. Such divergences could be interpreted in the frame of different scenarios of emission of cometary parents, either from a nucleus or from a volume source. The daughter scale-lengths for these comets were quite similar, namely: 306×¹⁰³ and 318×¹⁰³ km for CN and 69×¹⁰³ and 66×¹⁰³ km for C₂. We determined the Afρ parameter for apertures of different radii. A Monte Carlo model of the dust coma was used to obtain the dust ejection velocity. It was of the order of 0.1 km s⁻¹ for both comets. The power index of the distribution of the β-parameter of dust particles (ratio of light pressure to the solar gravitation) was of the order of 3 for C/Shoemaker-Levy and close to 2 for P/deVico. The dependence on heliocentric distance (rh) of the radical and dust production rates for P/deVico in the range of 0.7-1.0 AU was described by the power law function with a power index equal to: 5.55±0.14 for CN, 5.70±0.24 for C₂ and 5.22±0.19 for dust. Relative abundances of the dynamically new Comet Shoemaker-Levy and short-period P/deVico were quite similar with an enhancement of C₂ comparing with standard values taken from A'Hearn et al. (1995).     

Are There Many Inactive Jupiter-Family Comets among the Near-Earth Asteroid Population?

We analyze the dynamical evolution of Jupiter-family (JF) comets and near-Earth asteroids (NEAs) with aphelion distances Q>3.5 AU, paying special attention to the problem of mixing of both populations, such that inactive comets may be disguised as NEAs. From numerical integrations for 2×10⁶ years we find that the half lifetime (where the lifetime is defined against hyperbolic ejection or collision with the Sun or the planets) of near-Earth JF comets (perihelion distances q<1.3 AU) is about 1.5×10⁵ years but that they spend only a small fraction of this time (∼ a few 10³ years) with q<1.3 AU. From numerical integrations for 5×10⁶ years we find that the half lifetime of NEAs in “cometary” orbits (defined as those with aphelion distances Q>4.5 AU, i.e., that approach or cross Jupiter's orbit) is 4.2×10⁵ years, i.e., about three times longer than that for near-Earth JF comets. We also analyze the problem of decoupling JF comets from Jupiter to produce Encke-type comets. To this end we simulate the dynamical evolution of the sample of observed JF comets with the inclusion of nongravitational forces. While decoupling occurs very seldom when a purely gravitational motion is considered, the action of nongravitational forces (as strong as or greater than those acting on Encke) can produce a few Enckes. Furthermore, a few JF comets are transferred to low-eccentricity orbits entirely within the main asteroid belt (Q<4 AU and q>2 AU). The population of NEAs in cometary orbits is found to be adequately replenished with NEAs of smaller Q's diffusing outward, from which we can set an upper limit of ∼20% for the putative component of deactivated JF comets needed to maintain such a population in steady state. From this analysis, the upper limit for the average time that a JF comet in near-Earth orbit can spend as a dormant, asteroid-looking body can be estimated to be about 40% of the time spent as an active comet. More likely, JF comets in near-Earth orbits will disintegrate once (or shortly after) they end their active phases.     

Cometary dust properties retrieved from polarization observations: Application to C/1995 O1 Hale-Bopp and 1P/Halley

The analysis of the polarized light scattered by cometary dust particles provides information on the physical properties of the solid component of cometary comae for C/1995 O1 Hale-Bopp and 1P/Halley. A model of light scattering by a size distribution of aggregates of up to 256 submicron-sized grains (spherical or spheroidal) mixed with single spheroidal particles has been developed, with its parameters adjusted to fit the phase angle and wavelength dependence of the polarization observations. The particles are built of two materials: a non-absorbing silicates-type material and a more absorbing organic-type material. The model reproduces accurately the inversion angle and the positive branch of the polarization phase curves from the visible to the near-infrared spectral domains. A negative branch of the polarization phase curves appears in our model, although the negative branch is not deep enough to reproduce accurately the observations. Significant differences are shown between the two comets, with dominance of small grains in the coma of Comet C/1995 O1 Hale-Bopp, well fitted by a distribution of the volume-equivalent diameter, a, following a^−3.0 with a lower cutoff around 0.20 μm and an upper cutoff of at least 40 μm. For 1P/Halley, the size distribution follows a^−2.8 with a lower cutoff around 0.26 μm and an upper cutoff of about 38 μm. The relative amount of organic-type particles is larger for 1P/Halley while the amount of aggregates, significant for both comets, is larger for C/1995 O1 Hale-Bopp.     

Evolution of comet orbits under the perturbing influence of the giant planets and nearby stars

The orbital evolution of 500 hypothetical comets during 10⁹ years is studied numerically. It is assumed that the birthplace of such comets was the region of Uranus and Neptune from where they were deflected into very elongated orbits by perturbations of these planets. Then, we adopted the following initial orbital elements: perihelion distances between 20 and 30 AU, inclinations to the ecliptic plane smaller than 20°, and semimajor axes from 5 × 10³ to 5 × 10⁴ AU. Gravitational perturbations by the four giant planets and by hypothetical stars passing at distances from the Sun smaller than 5 × 10⁵ AU are considered. During the simulation, somewhat more than 50% of the comets were lost from the solar system due to planetary or stellar perturbations. The survivors were removed from the planetary region and left as members of what is generally known as the cometary cloud. At the end of the studied period, the semimajor axes of the surviving comets tend to be concentrated in the interval 2 × 10⁴ < a < 3 × 10⁴ AU. The orbital planes of the comets with initial $\text{a ≧ 3 × 10}{}^4\text{AU}$ acquired a complete randomization while the others still maintain a slight predominance of direct orbits. In addition, comet orbits with final a < 6 × 10⁴AU preserve high eccentricities with an average value greater than 0.8 Most “new” comets from the sample entering the region interior to Jupiter's orbit had already registered earlier passages through the planetary region. By scaling up the rate of paritions of hypothetical new comets with the observed one, the number of members of the cometary cloud is estimated to be about 7 × 10¹⁰ and the conclusion is drawn that Uranus and Neptune had to remove a number of comets ten times greater.     

Hilda asteroids among Jupiter family comets

Hilda asteroids and comets are similar from the compositional point of view. The D-taxonomic class prevailing among Hildas has all the characteristics found in cometary spectra. Jupiter Family Comets (JFCs) coming from the trans-neptunian region are under the gravitational control of Jupiter, making them a dynamically unstable population with a mean dynamical lifetime of 10⁴ to ¹⁰⁵ years. In contrast, Hilda asteroids residing in the 3:2 mean motion resonance with Jupiter are a very stable population. But once they escape from the resonance, they are dynamically controlled by Jupiter, and in this sense their behavior resembles that of JFC. We performed a numerical simulation to analyze the dynamical evolution that Hildas follow after escaping from the resonance, and their contribution to the JFC population. We found that 8% of the particles leaving the resonance end up impacting Jupiter. 98.7% of the escaped Hildas live at least 1000 years as a JFC, with a mean lifetime of 1.4×¹⁰⁶ years. In particular, escaped Hildas stay mainly in the region of perihelion distances greater than 2.5 AU. On the other hand, the number of escaped Hildas reaching the inner Solar System (q<2.5 AU) is negligible. So, there are almost no Hilda asteroids among the NEO population. We also analyzed the possibility that the Shoemaker-Levy 9 were an escaped Hilda asteroid. In this case, it would be possible to give stronger constraints to its pre-capture orbital elements.     

Dust Trails of 8P/Tuttle and the Unusual Outbursts of the Ursid Shower

We calculate the position of dust trails from comet 8P/Tuttle, in an effort to explain unusual Ursid meteor shower outbursts that were seen when the comet was near aphelion. Comet 8P/Tuttle is a Halley-type comet in a 13.6-year orbit, passing just outside of Earth's orbit. We find that the meteoroids tend to be trapped in the 12:14 mean motion resonance with Jupiter, while the comet librates in a slightly shorter period orbit around the 13:15 resonance. It takes 6 centuries to decrease the perihelion of the meteoroid orbits enough to intersect Earth's orbit, during which time the meteoroids and comet separate in mean anomaly by 6 years, thus explaining the 6-year lag between the comet's return and Ursid outbursts. The resonances also prevent dispersion along the comet orbit and limit viewing to only one year in each return. We identified past dust trail encounters with dust trails from 1392 (Dec. 1945) and 1378 (Dec. 1986) and predicted another outburst on 2000 December 22 at around 7:29 and 8:35 UT, respectively, from dust trails dating to the 1405 and 1392 returns. This event was observed from California using video and photographic techniques. At the same time, five Global-MS-Net stations in Finland, Japan, and Belgium counted meteors using forward meteor scatter. The outburst peaked at 8:06±07 UT, December 22, at zenith hourly rate ∼90 per hour, and the Ursid rates were above half peak intensity during 4.2 h. We find that most Ursid orbits do scatter around the anticipated positions, confirming the link with comet 8P/Tuttle and the epoch of ejection. The 1405 and 1392 dust trails appear to have contributed similar amounts to the activity profile. Some orbits provide a hint of much older debris being present as well. This work is the strongest evidence yet for the relevance of mean motion resonances in Halley-type comet dust trail evolution.     

Mechanisms accounting for variations in the proportions of carbonaceous and ordinary chondrites in different mass ranges

The number ratio of carbonaceous to ordinary chondrites (the CC/OC ratio) varies with mass. It is very high (≳90) in small mass ranges (10⁻⁸ to 10⁻¹² kg) among interplanetary dust particles and micrometeorites; it is moderately high (~5 to 30) for 1 to 10 m size fireball meteoroids (with estimated masses between ~10³ and ~10⁶ kg). In the range of most normal-sized meteorite falls (0.01–20 kg), the ratio is low (0.04–0.05); the ratio increases at greater mass ranges: at ≥200 kg, the ratio is 0.09; at ≥500 kg, the ratio is 0.20. The CC/OC ratio also increases from 0.05 to 0.16 for small meteorite finds (10⁻³ to 10⁻⁴ kg). High CC/OC ratios at low and high mass ranges are due to the predominance of CC material in the outer solar system. Small particles from this region spiral into the inner solar system typically in ≤10⁶ years due to Poynting–Robertson drag. Meter-sized meteoroids in this region are affected by Yarkovsky forces, pushing them into resonances where they are efficiently transferred to the inner solar system. Normal-sized meteorites are derived from centimeter-to-decimeter-sized meteoroids that have sluggish drift rates (i.e., they are less affected by the seasonal Yarkovsky effect) compared to larger bodies. Consequently, the centimeter-to-decimeter-sized meteoroids spend more time in interplanetary space (where they are subject to collisions) than larger objects. The greater friability of carbonaceous chondrites relative to ordinary chondrites tends to winnow the carbonaceous chondrites out in this size/mass range during their long interplanetary sojourn, thereby decreasing the CC/OC ratio.