Deep Space 1 photometry of the nucleus of Comet 19P/Borrelly

The NASA-JPL Deep Space 1 Mission (DS1) encountered the short-period Jupiter-family Comet 19P/Borrelly on September 22, 2001, about 8 days after perihelion. DS1's payload contained a remote-sensing package called MICAS (Miniature Integrated Camera Spectrometer) that included a 1024 square CCD and a near IR spectrometer with ∼12 nm resolution. Prior to its closest approach of 2171 km, the remote-sensing package on the spacecraft obtained 25 CCD images of the comet and 45 near-IR spectra (L. Soderblom et al., 2002, Science 296, 1087-1091). These images provided the first close-up view of a comet's nucleus sufficiently unobscured to perform quantitative photometric studies. At closest approach, corresponding to a resolution of 47 meters per pixel, the intensity of the coma was less than 1% of that of the nucleus. An unprecedented range of high solar phase angles (52-89 degrees), viewing geometries that are in general attainable only when a comet is active, enabled the first quantitative and disk resolved modeling of surface photometric physical parameters, including the single particle phase function and macroscopic roughness. The disk-integrated geometric albedo of Borrelly's nucleus is 0.029±0.006, comparable to the dark hemisphere of Iapetus, the lowest albedo C-type asteroids, and the uranian rings. The Bond albedo, 0.009±0.002, is lower than that of any Solar System object measured. Such a low value may enhance the heating of the nucleus and sublimation of volatiles, which in turn causes the albedo to decrease even further. A map of normal reflectance of Borrelly shows variations far greater than those seen on asteroids. The two main terrain types, smooth and mottled, exhibit mean normal reflectances of 0.03 and 0.022. The physical photometric parameters of Borrelly's nucleus are typical of other small dark bodies, particularly asteroids, except preliminary modeling results indicate its regolith may be substantially fluffier. The nucleus exhibits significant variations in macroscopic roughness, with the oldest, darkest terrain being slightly smoother. This result suggests the infilling of low-lying areas with dust and particles that have not been able to leave the comet. The surface of the comet is backscattering, but there are significant variations in the single particle phase function. One region exhibits a flat particle phase function between solar phase angles of 50° and 75° (like cometary dust and unlike planetary surfaces), suggesting that its regolith is controlled by native dust rather than by meteoritic bombardment.     

Deep Space 1 at comet 19P/Borrelly: Magnetic field and plasma observations

On September 22, 2001 the Deep Space 1 spacecraft performed a flyby at comet 19P/Borrelly at a solar distance of 1.36 AU leading the Earth by 74^° in longitude. The spacecraft-comet distance at closest approach was 2171 km. The bow shock had a magnetic compression ratio of 2.5 at a distance of 147 100 km from the nucleus. Deep Space 1 first entered the sheath region essentially from the north polar region. Fluctuations from the cometary ion pickup were present throughout the sheath region and even well upstream of the shock, as expected. The magnetic field pileup region had a peak field strength of 83 nT and was shown to be consistent with a pressure equal to the solar wind ram pressure. The peak field location was offset from the time of closest approach. It is uncertain whether this is a spatial or temporal variation. Draping of magnetic fields around the nucleus was sought, but evidence for this was not apparent in the data. A possible explanation is that the interplanetary solar wind was composed of turbulent short-scale fields, and thus the fields were not symmetric about the point of closest approach. During the flyby phase there were in general few intervals of ACE data where there were large scale Parker spiral fields. With the addition of plasma data, the shock properties are investigated. The characteristics of magnetic draping, pileup and fluctuations are explored. These comet 19P/Borrelly results are contrasted with other cometary flyby results.     

Near-infrared spectroscopy of 3:1 Kirkwood Gap Asteroids: (3760) Poutanen and (974) Lioba

A mineralogical assessment of 3:1 Kirkwood Gap Asteroids (3760) Poutanen and (974) Lioba, has been carried out from spectral data obtained May 21, 2008 (974) and May 25, 2009 (3760) UT using the NASA Infrared Telescope Facility. Results indicate (3760) Poutanen as a spinel-bearing asteroid. The spinel contained within this asteroid appears to be embedded in calcium aluminum inclusions (CAIs) much like that of the CV3 chondrite Allende. The preservation of the CAIs implies a lack of igneous processing and places further constraints on the heating that occurred within the early Solar System. This analysis suggests (3760) Poutanen originated in the general formation zone of the CV3 meteorites, but in a sub-region depleted in iron-rich matrix material and enriched in pyroxene-bearing clasts. Analysis of (974) Lioba indicates it is an S(IV) asteroid with a surface silicate assemblage consisting of olivine and orthopyroxene, (Fs_21(±5)Wo_8(±3)), consistent with measured band centers, (0.96 ± 0.01 μm; 1.95 ± 0.02 μm), and the band area ratio (0.59 ± 0.15). The location, delivery efficiency ((24.5%) Farinella, P., Gunczi, R., Froeschlé, Ch., Froeschlé, C. [1993]. Icarus 101, 174-187), calculated mineralogy, and BAR strongly suggest (974) Lioba is a plausible parent body for the LL- or L-ordinary chondrites.     

Size, albedo, and taxonomic type of potential spacecraft target Asteroid (10302) 1989 ML

The Amor-type near-Earth Asteroid (10302) 1989 ML has an “Earth-like” orbit (period 1.44 yr, eccentricity 0.14, inclination 4.4°), therefore the energy required to reach it from the Earth is relatively small making it a very attractive target for rendezvous missions. We have observed 1989 ML in the thermal-infrared using the Spitzer Space Telescope, and compared these data with optical and near-infrared observations. The Spitzer results imply a diameter of 0.28±0.05 km and a geometric albedo of 0.37±0.15; together with the reflectance spectrum they are consistent with the relatively rare E classification.     

Visible and near infrared spectroscopic investigation of E-type asteroids, including 2867 Steins, a target of the Rosetta mission

We present the results of a visible spectroscopic survey of igneous asteroids belonging to the small and intriguing E-class, including 2867 Steins, a target of the Rosetta mission. The survey was carried out at the 3.5 m Telescopio Nazionale Galileo (TNG), and at the 3.5 m New Technology Telescope (NTT) of the European Southern Observatory. We obtained new visible spectra for eighteen E-type asteroids, and near infrared spectra for eight of them. We confirm the presence of three different mineralogies in the small E-type populations. We classify each object in the E[I], E[II] or E[III] subgroups [Gaffey, M.J., Kelley, M.S., 2004. Lunar Planet. Sci. XXXV. Abstract 1812] on the basis of the spectral behavior and of the eventual presence of absorption features attributed to sulfides (such the 0.49 μm band, on E[II]), or to iron bearing silicates (0.9 μm band, on E[III]). We suggest that some asteroids (i.e. 64 Angelina, 317 Roxane, and 434 Hungaria), which show different spectral behavior comparing our data with those available in literature, have an inhomogeneous surface composition. 2867 Steins, a target of the Rosetta mission, shows a spectral behavior typical of the E[II] subgroup, as already suggested by Barucci et al. [Barucci, M.A., Fulchignoni, M., Fornasier, S., Dotto, E., Vernazza, P., Birlan, M., Binzel, R.P., Carvano, J., Merlin, F., Barbieri, C., Belskaya, I., 2005. Astron. Astrophys. 430, 313-317] and Fornasier et al. [Fornasier, S., Marzari, F., Dotto, E., Barucci, M.A., Migliorini, A., 2007. Astron. Astrophys. 474, 29-32]. Litva and 1990 TN1, initially classified as E-types, show a visible and near infrared behavior consistent with the olivine rich A-class asteroids, while 5806 Archieroy, also supposed to belong to the E-class, has a spectral behavior consistent with the S(V) classification following the Gaffey et al. [Gaffey, M.J., Burbine, T.H., Piatek, J.L., Reed, K.L., Chaky, D.A., Bell, J.F., Brown, R.H., 1993. Icarus 106, 573-602] classification scheme. To fully investigate the E-type population, we enlarged our sample including 6 E-type asteroids spectra available in literature, resulting in a total sample of 21 objects. The analysis of the spectral slope for the 3 different E-type subgroups versus the orbital elements show that E[III] members have the lowest mean spectral slope value inside the whole sample, and that they are located between 2.2-2.7 AU in low inclination orbits. E[II] members has the highest spectral slope inside the sample, half of them are located in the Hungaria region, 2 are NEA and 2 (64 Angelina and 2867 Steins), are in the main belt. A similar distribution is found for the 5 featureless E[I] members, located mainly in the Hungaria region (3 members), one in the middle main belt while one is a NEA (2004 VD17). Finally, for the five E-type asteroids observed both in the visible and near infrared range, plus 2867 Steins, we attempt to model their surface composition using linear geographical mixtures of no more than 3 components, selected from aubrite meteorites and correlated minerals. In particular we suggest that the aubrite Peña Blanca might have the E[III] Asteroid 317 Roxane as parent body, and that the aubrite ALH78113 might be related to the E[II] subgroup asteroids.     

Deep Impact: Optical spectroscopy and photometry obtained at MIRA

We present spectroscopic and photometric observations, spanning the optical UV to the far red, before, during, and after the NASA Deep Impact event of July 4, 2005. The inner 2000 km of the pre- and post-impact coma was about 0.3 magnitude redder in B-R than in the outer coma. The pre-impact spectrum was a faint reflected solar spectrum dominated by molecular emissions extending >40,000 km from the nucleus. The post-impact light curve in R and I showed a rapid rise consistent with an expanding optically thick cloud during the first 18 min after impact. During the next 8 min the cloud became optically thin. Sixty minutes after impact the impact R-band flux reached a plateau at , the comet brightening by a factor of ∼4.3 above its pre-impact value observed in a 15″ aperture. The mean expansion velocity of the grains during the first 49 min was . The spectrum became dominated by scattered sunlight during the first hour after impact. The volume scattering function (VSF) observed 32 min after impact shows strong reddening. At 49 min, however, the VSF shows an additional twofold increase in the blue but only a 20% increase at 5500 Å. Post-impact spectra and R-I photometry showed rapid reddening. The particle size distribution, dominated by 1-2.5 μm particles shortly after impact, changed dramatically during the first hour due to sublimation of water-ice particles of this size. On the night following impact the comet was still substantially brighter than before impact, but R-I had returned to its pre-impact value. B-R remained significantly redder. The ejecta 25 h after impact was fan-shaped subtending ∼180° roughly symmetrical about position angle 225°. The mean expansion velocity 90° from the direction to the Sun was .     

The space infrared telescope for cosmology and astrophysics: SPICA A joint mission between JAXA and ESA

The Space Infrared telescope for Cosmology and Astrophysics (SPICA) is planned to be the next space astronomy mission observing in the infrared. The mission is planned to be launched in 2017 and will feature a 3.5 m telescope cooled to <5 K through the use of mechanical coolers. These coolers will also cool the focal plane instruments thus avoiding the use of consumables and giving the mission a long lifetime. SPICA’s large, cold aperture will provide a two order of magnitude sensitivity advantage over current far infrared facilities (>30 microns wavelength). We describe the scientific advances that will be made possible by this large increase in sensitivity and give details of the mission, spacecraft and focal plane conceptual design.

A note concerning the 2:1 and the 3:2 resonances in the asteroid belt

The dynamical behavior of asteroids inside the 2:1 and 3:2 commensurabilities with Jupiter presents a challenge. Indeed most of the studies, either analytical or numerical, point out that the two resonances have a very similar dynamical behavior. In spite of that, the 3:2 resonance, a little outside the main belt, hosts a family of asteroids, called the Hildas, while the 2:1, inside the main belt, is associated to a gap (the Hecuba gap) in the distribution of asteroids.In his search for a dynamical explanation for the Hecuba gap, Wisdom (1987) pointed out the existence of orbits starting with low eccentricity and inclination inside the 2:1 commensurability and going to high eccentricity, and thus to possible encounters with Mars. It has been shown later (Henrard et al.), that these orbits were following a path from the low eccentric belt of secondary resonances to the high eccentric domain of secular resonances. This path crosses a bridge, at moderate inclination and large amplitude of libration, between the two chaotic domains associated with these resonances.The 3:2 resonance being similar in many respects to the 2:1 resonance, one may wonder whether it contains also such a path. Indeed we have found that it exists and is very similar to the 2:1 one. This is the object of the present paper.     

Ultraviolet spectroscopy of Comet 9P/Tempel 1 with Alice/Rosetta during the Deep Impact encounter

We report on spectroscopic observations of periodic Comet 9P/Tempel 1 by the Alice ultraviolet spectrograph on the Rosetta spacecraft in conjunction with NASA's Deep Impact mission. Our objectives were to measure an increase in atomic and molecular emissions produced by the excavation of volatile sub-surface material. We unambiguously detected atomic oxygen emission from the quiescent coma but no enhancement at the 10% (1-σ) level following the impact. We derive a quiescent H₂O production rate of 9×¹⁰²⁷ molecules s⁻¹ with an estimated uncertainty of ∼30%. Our upper limits to the volatiles produced by the impact are consistent with other estimates.     

Comet Hale-Bopp (C/1995 O1): UVSTAR-FUV Spectroscopy from the Space Shuttle

We report here on unique post-perihelion (2.3 AU) measurements of Comet Hale-Bopp in the FUV-range (950–1250 Å) by means of the UVSTAR spectrometer from the space shuttle with the main purpose of searching for argon and other FUV emitters. New methods for separating the strong airglow emission at shuttle altitudes are here discussed in detail. Due to our low resolution (15 Å) and S/N ratio the possible rocket-borne detection of argon near perihelion (0.9 AU) could not be confirmed. New species as N₂ are suspected but difficult to separate from the strong airglow emission at shuttle altitudes. From the Lyα brightness (1.30± 0.08 kRy) a water production rate Q = 5.9 ± 0.4 × 10²⁹ molecules s^?1 could be derived and compared with other post-perihelion observations.     

The shape, topography, and geology of Tempel 1 from Deep Impact observations

Deep Impact images of the nucleus of Comet Tempel 1 reveal pervasive layering, possible impact craters, flows with smooth upper surfaces, and erosional stripping of material. There are at least 3 layers 50-200 m thick that appear to extend deep into the nucleus, and several layers 1-20 m thick that parallel the surface and are being eroded laterally. Circular depressions show geographical variation in their forms and suggest differences in erosion rates or style over scales >1 km. The stratigraphic arrangement of these features suggests that the comet experienced substantial periods of little erosion. Smooth surfaces trending downslope suggest some form of eruption of materials from this highly porous object. The Deep Impact images show that the nucleus of Tempel 1 cannot be modeled simply as either an onion-layer or rubble pile structure.     

Imaging polarimetry of Comet 9P/Tempel before and after the Deep Impact

The NASA's Deep Impact mission was the first impact experiment to a cometary nucleus. The target of the mission was Comet 9P/Tempel, one of the Jupiter family comets. The impact was performed on July 4th, 2005. Imaging polarimetric observations were carried out by Polarimetric Imager for COmets (PICO) mounted on the Lulin One-meter Telescope (LOT) at Lulin Observatory, Taiwan. Intensity and linear polarization degree maps were obtained on July 3-5, 2005. Impact ejecta plume was clearly recognized in the enhanced intensity map. Furthermore, arc-shaped region of high polarization was recognized in the polarization map. Dust grains in this region had larger expansion velocity than the grains which provided the brightest area in the intensity map. comparing our results with the MIR spectroscopy obtained by Subaru Telescope we conclude that very small carbonaceous grains might be responsible for the region of high polarization.     

The Maria asteroid family: Genetic relationships and a plausible source of mesosiderites near the 3:1 Kirkwood Gap

We present a mineralogical assessment of 12 Maria family asteroids, using near-infrared spectral data obtained over the years 2000-2009 combined with visible spectral data (when available) to cover the spectral interval of 0.4-2.5 μm. Our analysis indicates the Maria asteroid family, which is located adjacent to the chaotic region of the 3:1 Kirkwood Gap, appears to be a true genetic family composed of assemblages analogous to mesosiderite-type meteorites. Dynamical models by Farinella et al. (Farinella, P., Gunczi, R., Froeschlé, Ch., Froeschlé, C., [1993]. Icarus 101, 174-187) predict this region should supply meteoroids into Earth-crossing orbits. Thus, the Maria family is a plausible source of some or all of the mesosiderites in our meteorite collections. These individual asteroids were most likely once part of a larger parent object that was broken apart and dispersed. One of the Maria dynamical family members investigated, ((695) Bella), was found to be unrelated to the genetic Maria family members. The parameters of (695) Bella indicate an H-chondrite assemblage, and that Bella may be a sister or daughter of Asteroid (6) Hebe.     

Spin vectors in the Koronis family: II. Additional clustered spins, and one stray

We recorded 101 new rotation lightcurves of five Koronis family members, and then combined the new observations with previous data to determine the objects' sidereal rotation periods, spin vector orientations, and model shape solutions. The observing program was tailored specifically for spin vector analyses by determining single-apparition Lumme–Bowell solar phase coefficients, and by measuring synodic rotation periods precisely enough to unambiguously count the rotations between two consecutive oppositions, which is a prerequisite for identifying the correct sidereal period. The new data make possible first pole and shape determinations for (263) Dresda, (462) Eriphyla, and (1289) Kutaïssi, and they improve the models for (277) Elvira and (534) Nassovia, two objects previously studied by Slivan et al. [Slivan, S.M., Binzel, R.P., Crespo da Silva, L.D., Kaasalainen, M., Lyndaker, M.M., Kr?o, M., 2003. Icarus 162, 285–307]. Our results increase the number of Koronis family spin vectors reported in the literature to fourteen, a sample which now includes the seven largest family members. The spin properties of Eriphyla (rotation period , spin vector obliquity ε=51°) and Kutaïssi (P=3.62 h, ε=165°) are consistent with the markedly nonrandom distribution reported by Slivan [Slivan, S.M., 2002. Nature 419, 49–51], and explained by Vokrouhlický et al. [Vokrouhlický, D., Nesvorný, D., Bottke, W.F., 2003. Nature 425, 147–151] as the result of the effects of thermal “YORP” torques combined with solar and planetary gravitational torques. Dresda (P=16.81 h, ε=16°) is the first prograde Koronis member whose spin obliquity and spin rate significantly differ from the clustered spin properties previously found for other prograde Koronis members; nevertheless, its spin vector is consistent with several of the spin evolution possibilities that were identified in the YORP modeling.     

PEGASE, an infrared interferometer to study stellar environments and low mass companions around nearby stars

PEGASE is a mission dedicated to the exploration of the environment (including habitable zone) of young and solar-type stars (particularly those in the DARWIN catalogue) and the observation of low mass companions around nearby stars. It is a space interferometer project composed of three free flying spacecraft, respectively featuring two 40 cm siderostats and a beam combiner working in the visible and near infrared. It has been proposed to ESA as an answer to the first “Cosmic Vision” call for proposals, as an M mission. The concept also enables full-scale demonstration of space nulling interferometry operation for DARWIN.

A Heliospheric Imager for Deep Space: Lessons Learned from Helios, SMEI, and STEREO

The zodiacal-light photometers on the twin Helios spacecraft, the Solar Mass Ejection Imager (SMEI) on the Coriolis spacecraft, and the Heliospheric Imagers (HIs) on the Solar-TErrestrial RElations Observatory (STEREO) twin spacecraft all point the way to optimizing future remote-sensing Thomson-scattering observations from deep space. Such data could be provided by wide-angle viewing instruments on Solar Orbiter, Solar Probe, or other deep-space probes. Here, we present instrument specifications required for a successful heliospheric imager, and the measurements and data-processing steps that make the best use of this remote-sensing system. When this type of instrument is properly designed and calibrated, its data are capable of determining zodiacal-dust properties, and of three-dimensional reconstructions of heliospheric electron density over large volumes of the inner heliosphere. Such systems can measure fundamental properties of the inner heliospheric plasma, provide context for the in-situ monitors on board spacecraft, and enable physics-based analyses of this important segment of the Sun-spacecraft connection.     

Numerical investigation on the orbital evolution of 1:2 resonant asteroids between Jupiter and Saturn

According to some investigations (Lecar and Franklin, 1973; Franklin et al., 1989; Soper et al., 1990) asteroids cannot remain for along time between Jupiter and Saturn. But as it is well known there is a near 5:2 commensurability between Jupiter and Saturn. So there might be a possibility that asteroids between Jupiter and Saturn could be trapped in a resonant relation.In order to investigate this possibility, the changes of orbital elements of an asteroid whose initial value of semi-major axis corresponds to that of a 1:2 resonant orbit were investigated by means of a double precision Cowell method. The integration routine was kindly supplied by Dr Yoshikawa.We considered first a planar restricted problem of three bodies, Sun-Jupiter-Asteroid, then a four body problem, Sun-Jupiter-Asteroid-Saturn. When integrating the equations of motion, short periodic terms were not eliminated and in the second test the interactions between Jupiter and Saturn were retained. Whether a close approach occured or not was not investigated. In every case a _j = 5.20, a _s = 9.54 and a = 8.26 were adopted as initial values of the semi-major axis of Jupiter, Saturn and Asteroid respectively.     

Shape, size and multiplicity of main-belt asteroids: I. Keck Adaptive Optics survey

This paper presents results from a high spatial resolution survey of 33 main-belt asteroids with diameters >40 km using the Keck II Adaptive Optics (AO) facility. Five of these (45 Eugenia, 87 Sylvia, 107 Camilla, 121 Hermione, 130 Elektra) were confirmed to have satellite. Assuming the same albedo as the primary, these moonlets are relatively small (∼5% of the primary size) suggesting that they are fragments captured after a disruptive collision of a parent body or captured ejecta due to an impact. For each asteroid, we have estimated the minimum size of a moonlet that can positively detected within the Hill sphere of the system by estimating and modeling a 2-σ detection profile: in average on the data set, a moonlet located at 2/100×RHill (1/4×RHill) with a diameter larger than 6 km (4 km) would have been unambiguously seen. The apparent size and shape of each asteroid was estimated after deconvolution using a new algorithm called AIDA. The mean diameter for the majority of asteroids is in good agreement with IRAS radiometric measurements, though for asteroids with a D<200 km, it is underestimated on average by 6-8%. Most asteroids had a size ratio that was very close to those determined by lightcurve measurements. One observation of 104 Klymene suggests it has a bifurcated shape. The bi-lobed shape of 121 Hermione described in Marchis et al. [Marchis, F., Hestroffer, D., Descamps, P., Berthier, J., Laver, C., de Pater, I., 2005c. Icarus 178, 450-464] was confirmed after deconvolution. The ratio of contact binaries in our survey, which is limited to asteroids larger than 40 km, is surprisingly high (∼6%), suggesting that a non-single configuration is common in the main-belt. Several asteroids have been analyzed with lightcurve inversions. We compared lightcurve inversion models for plane-of-sky predictions with the observed images (9 Metis, 52 Europa, 87 Sylvia, 130 Elektra, 192 Nausikaa, and 423 Diotima, 511 Davida). The AO images allowed us to determine a unique photometric mirror pole solution, which is normally ambiguous for asteroids moving close to the plane of the ecliptic (e.g., 192 Nausikaa and 52 Europa). The photometric inversion models agree well with the AO images, thus confirming the validity of both the lightcurve inversion method and the AO image reduction technique.     

A ballistics analysis of the Deep Impact ejecta plume: Determining Comet Tempel 1's gravity, mass, and density

In July of 2005, the Deep Impact mission collided a 366 kg impactor with the nucleus of Comet 9P/Tempel 1, at a closing speed of 10.2 km s⁻¹. In this work, we develop a first-order, three-dimensional, forward model of the ejecta plume behavior resulting from this cratering event, and then adjust the model parameters to match the flyby-spacecraft observations of the actual ejecta plume, image by image. This modeling exercise indicates Deep Impact to have been a reasonably “well-behaved” oblique impact, in which the impactor-spacecraft apparently struck a small, westward-facing slope of roughly 1/3-1/2 the size of the final crater produced (determined from initial ejecta plume geometry), and possessing an effective strength of not more than . The resulting ejecta plume followed well-established scaling relationships for cratering in a medium-to-high porosity target, consistent with a transient crater of not more than 85-140 m diameter, formed in not more than 250-550 s, for the case of (gravity-dominated cratering); and not less than 22-26 m diameter, formed in not less than 1-3 s, for the case of (strength-dominated cratering). At , an upper limit to the total ejected mass of 1.8×¹⁰⁷ kg (1.5-2.2×¹⁰⁷ kg) is consistent with measurements made via long-range remote sensing, after taking into account that 90% of this mass would have stayed close to the surface and then landed within 45 min of the impact. However, at , a lower limit to the total ejected mass of 2.3×¹⁰⁵ kg (1.5-2.9×¹⁰⁵ kg) is also consistent with these measurements. The expansion rate of the ejecta plume imaged during the look-back phase of observations leads to an estimate of the comet's mean surface gravity of (0.17-0.90 mm s⁻²), which corresponds to a comet mass of mt=4.5×¹⁰¹³ kg (2.3-12.0×¹⁰¹³ kg) and a bulk density of (200-1000 kg m⁻³), where the large high-end error is due to uncertainties in the magnitude of coma gas pressure effects on the ejecta particles in flight.