(47171) 1999 TC36, A transneptunian triple

We present new analysis of HST images of (47171) 1999 TC₃₆ that confirm it as a triple system. Fits to the point-spread function (PSF) consistently show that the apparent primary is itself composed of two similar-sized components. The two central components, A1 and A2, can be consistently identified in each of nine epochs spread over 7 years of time. In each instance, the component separation, ranging from 0.023 ± 0.002 to 0.031 ± 0.003 arcsec, is roughly one half of the Hubble Space Telescope’s diffraction limit at 606 nm. The orbit of the central pair has a semi-major axis of a ∼ 867 km with a period of P ∼ 1.9 days. These orbital parameters yield a system mass that is consistent with M_sys = 12.75 ± 0.06 × 10¹⁸ kg derived from the orbit of the more distant secondary, component B. The diameters of the three components are . The relative sizes of these components are more similar than in any other known multiple in the Solar System. Taken together, the diameters and system mass yield a bulk density of . HST photometry shows that component B is variable with an amplitude of ?0.17 ± 0.05 magnitudes. Components A1 and A2 do not show variability larger than 0.08 ± 0.03 magnitudes approximately consistent with the orientation of the mutual orbit plane and tidally distorted equilibrium shapes. The system has high specific angular momentum of J/J′ = 0.93, comparable to most of the known transneptunian binaries.     

Energy balance of the Deep Impact experiment

We present results on the energy balance of the Deep Impact experiment based on analysis of 180 infrared spectra of the ejecta obtained by the Deep Impact spacecraft. We derive an output energy of 16.5 (+9.1/−4.1) GJ. With an input energy of 19.7 GJ, the error bars are large enough so that there may or may not be a balance between the kinetic energy of the impact and that of outflowing materials. Although possible, no other source of energy other than the impactor or the Sun is needed to explain the observations. Most of the energy (85%) goes into the hot plume in the first few seconds, which only represents a very small fraction (<0.01%) of the total ejected mass. The hot plume contains 190 (+263/−71) kg of H₂O, 1.6 ± 0.5 kg of CO₂, 8.2 (+11.3/3.1) kg of CO (assuming a CO/H₂O ratio of 4.3%), 27.9 (+25.0/−8.9) kg of organic material and 255 ± 128 kg of dust, while the ejecta contains ∼10⁷ kg of materials. About 12% of the energy goes into the ejecta (mostly water) and 3% to destroy the impactor. Volatiles species other than H₂O (CO₂, CO or organic molecules) contribute to <7% of the energy balance. In terms of physical processes, 68% of the energy is used to accelerate grains (kinetic energy), 16% to heat them, 6% to sublimate or melt them and 10% (upper limit) to break and compress dust and/or water ice aggregates into small micron size particles. For the hot plume, we derive a dust/H₂O ratio of 1.3 (+1.9/−1.0), a CO₂/H₂O ratio of 0.008 (+0.009/−0.006), an organics/H₂O ratio of 0.15 (+0.29/−0.11) and an organics/dust ratio of 0.11 (+0.30/−0.07). This composition refers to the impact site and is different from that of the bulk nucleus, consistent with the idea of layers of different composition in the nucleus sub-surface. Our results emphasize the importance of laboratory impact experiments to understand the physical processes involved at such a large scale.     

Photometry of Comet 9P/Tempel 1 during the 2004/2005 approach and the Deep Impact module impact

The results of the 9P/Tempel 1 CARA (Cometary Archive for Amateur Astronomers) observing campaign is presented. The main goal was to perform an extended survey of the comet as a support to the Deep Impact (DI) Mission. CCD R, I and narrowband aperture photometries were used to monitor the Afρ quantity. The observed behavior showed a peak of 310 cm 83 days before perihelion, but we argue that it can be distorted by the phase effect, too. The phase effect is roughly estimated around 0.0275 mag/degree, but we had no chance for direct determination because of the very similar geometry of the observed apparitions. The log-slope of Afρ was around −0.5 between about 180-100 days before the impact but evolved near the steady-state like 0 value by the impact time. The DI module impact caused about a 60% increase in the value of Afρ and a cloud feature in the coma profile which was observed just after the event. The expansion of the ejecta cloud was consistent with a fountain model with initial projected velocity of 0.2 km/s and β=0.73. Referring to a 25,000 km radius area centered on the nucleus, the total cross section of the ejected dust was 0.06 days after the impact, and 1.93 days after the impact (A is the dust albedo). Five days after the event no signs of the impact were detected, nor deviations from the expected activity referring both to the average pre-impact behavior and to the previous apparitions.     

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.     

Heavy ion escape from Mars, influence from solar wind conditions and crustal magnetic fields

We have used more than 4 years of Mars Express ion data to estimate the escape of heavy ions ( and ) from Mars. To take the limited field of view of the instrument into account, the data has been binned into spatial bins and angular bins to create average distribution functions for different positions in the near Mars space. The net escape flux for the studied low solar activity period, between May 2007 and May 2011, is 2.0 ± 0.2 × 10²⁴ s⁻¹. The escape has been calculated independently for four different quadrants in the Y_MSO − Z_MSO plane, south, dusk, north and dawn. Escape is highest from the northern and dusk quadrants, 0.6 ± 0.1 × 10²⁴ s⁻¹, and smallest from the south and dawn quadrants, 0.4 ± 0.1 × 10²⁴ s⁻¹. The flux ratio of molecular ( and ) to O⁺ ions is 0.9 ± 0.1, averaged over all quadrants. The flux difference between the north and south quadrants is statistically significant, and is presumed to be due to the presence of significant crustal magnetic fields in the southern hemisphere, reducing the outflow. The difference between the dawn and dusk quadrants is likely due to the magnetic tension associated with the nominal Parker angle spiral, which should lead to higher average magnetic tension on the dusk side. The escape increases during periods of high solar wind flux and during times when co-rotating interaction regions (CIR) affect Mars. In the latter case the increase is a factor 2.4-2.9 as compared to average conditions.     

Explosion of Comet 17P/Holmes as revealed by the Spitzer Space Telescope

An explosion on Comet 17P/Holmes occurred on 2007 October 23, projecting particulate debris of a wide range of sizes into the interplanetary medium. We observed the comet using the mid-Infrared Spectrograph (5-40 μm), on 2007 November 10 and 2008 February 27, and the imaging photometer (24 and 70 μm), on 2008 March 13, on board the Spitzer Space Telescope. The 2007 November 10 spectral mapping revealed spatially diffuse emission with detailed mineralogical features, primarily from small crystalline olivine grains. The 2008 February 27 spectra, and the central core of the 2007 November 10 spectral map, reveal nearly featureless spectra, due to much larger grains that were ejected from the nucleus more slowly. Optical images were obtained on multiple dates spanning 2007 October 27-2008 March 10 at the Holloway Comet Observatory and 1.5-m telescope at Palomar Observatory. The images and spectra can be segmented into three components: (1) a hemispherical shell fully 28′ on the sky in 2008 March, due to the fastest (262 m s⁻¹), smallest (2 μm) debris, with a mass ; (2) a ‘blob’ or ‘pseudonucleus’ offset from the true nucleus and subtending some 10′ on the sky, due to intermediate speed (93 m s⁻¹) and size (8 μm) particles, with a total mass ; and (3) a ‘core’ centered on the nucleus due to slower (9 m s⁻¹), larger (200 μm) ejecta, with a total mass . This decomposition of the mid-infrared observations can also explain the temporal evolution of the millimeter-wave flux. The orientation of the leading edge of the ejecta shell and the ejecta ‘blob,’ relative to the nucleus, do not change as the orientation of the Sun changes; instead, the configuration was imprinted by the orientation of the initial explosion. The distribution and speed of ejecta implies an explosion in a conical pattern directed approximately in the solar direction on the date of explosion. The kinetic energy of the ejecta >10²¹ erg is greater than the gravitational binding energy of the nucleus. We model the explosion as being due to crystallization and release of volatiles from interior amorphous ice within a subsurface cavity; once the pressure in the cavity exceeded the surface strength, the material above the cavity was propelled from the comet. The size of the cavity and the tensile strength of the upper layer of the nucleus are constrained by the observed properties of the ejecta; tensile strengths on >10 m scale must be greater than 10 kPa (or else the ejecta energy exceeds the binding energy of the nucleus) and they are plausibly 200 kPa. The appearance of the 2007 outburst is similar to that witnessed in 1892, but the 1892 explosion was less energetic by a factor of about 20.     

R-band light curve of Comet 9P/Tempel 1 during the Deep Impact event

We present high-speed CCD photometry of Comet 9P/Tempel 1 during the Deep Impact event on 2005 July 4 UT. Approximately 2 h and 50 min of R-band data were acquired at Mount Laguna Observatory with a temporal resolution of 5.5 s. The flux increased by 9% in the first minute after impact. This was followed by a more gradual two-part linear rise, with a change in slope at 9.2 min post-impact, at which time the rate of brightening increased from ∼ to ∼. An analysis of the light curve obtained with the guide camera on the United Kingdom Infrared Telescope and yields very similar results. These findings are mildly in disagreement with the 3-part linear rise found by Fernández et al. (2007) in that we do not find any evidence for a change at 4 min post-impact. We interpret the linear rise phase as due to solar illumination of the edge of an expanding optically thick dust ejecta plume. After approximately 20 min, the light curves begin to flatten out, perhaps coincident with the start of the transition to becoming optically thin. In the large apertures (>10^″) the light curve continues to gradually rise until the end of the observations. In smaller apertures, the light curves reach a peak at approximately 50 min, then decrease back towards the pre-impact flux level. The drop in flux in the smaller apertures may be caused by the ejecta expanding beyond the edge of the photometric aperture, and if so, we can use this timescale to infer an expansion velocity of ∼, consistent with previous published estimates.     

Mineralogy and thermal properties of V-type Asteroid 956 Elisa: Evidence for diogenitic material from the Spitzer IRS (5-35 μm) spectrum

We present the thermal infrared (5-35 μm) spectrum of 956 Elisa as measured by the Spitzer Infrared Spectrograph (“IRS”; Houck, J.R. et al. [2004]. Astrophys. J. Suppl. 154, 18-24) together with new groundbased lightcurve data and near-IR spectra. From the visible lightcurve photometry, we determine a rotation period of 16.494 ± 0.001 h, identify the rotational phase of the Spitzer observations, and estimate the visible absolute magnitude (H_V) at that rotational phase to be 12.58 ± 0.04. From radiometric analysis of the thermal flux spectrum, we find that at the time of observation 956 Elisa had a projected radius of 5.3 ± 0.4 km with a visible albedo p_V = 0.142 ± 0.022, significantly lower than that of the prototype V-type asteroid, 4 Vesta. (This corresponds to a radius of 5.2 ± 0.4 km at lightcurve mean.) Analysis with the standard thermal model (STM) results in a sub-solar temperature of 292.3 ±  2.8 K and beaming parameter η = 1.16 ± 0.05. Thermophysical modeling places a lower limit of on the thermal inertia of the asteroid’s surface layer (if the surface is very smooth) but more likely values fall between 30 and depending on the sense of rotation.The emissivity spectrum, calculated by dividing the measured thermal flux spectrum by the modeled thermal continuum, exhibits mineralogically interpretable spectral features within the 9-12 μm reststrahlen band, the 15-16.5 μm Si-O-Si stretching region, and the 16-25 μm reststrahlen region that are consistent with pyroxene of diogenitic composition: extant diogenitic pyroxenes fall within the narrow compositional range Wo_2±1En_74±2Fs_24±1. Spectral deconvolution of the 9-12 μm reststrahlen features indicates that up to ≈20% olivine may also be present, suggesting an olivine-diogenite-like mineralogy. The mid-IR spectrum is inconsistent with non-cumulate eucrite as the major component on the surface of 956 Elisa, although cumulate eucrite material may be present at abundances lower than that of the diogenite component.Analysis of new near-IR spectra of 956 Elisa with the Modified Gaussian Model (MGM; Sunshine, J.M., Pieters, C.M., Pratt, S.F. [1990]. J. Geophys. Res. 95 (May), 6955-6966) results in two pyroxene compositions: 75% magnesian low-Ca pyroxene and 25% high-Ca pyroxene. High-Ca pyroxene is not evident in the mid-IR data, but may belong to a component that is underrepresented in the mid-IR spectrum either because of its spatial distribution on the asteroid or because of its particle size. High-Ca pyroxenes that occur as exsolution lamellae may also be more evident spectrally in the NIR than in the mid-IR. In any case, we find that the mid-IR spectrum of 956 Elisa is dominated by emission from material of diogenite-like composition, which has very rarely been observed among asteroids.     

Micrometer-sized ice particles for planetary-science experiments - I. Preparation, critical rolling friction force, and specific surface energy

Coagulation models assume a higher sticking threshold for micrometer-sized ice particles than for micrometer-sized silicate particles. However, in contrast to silicates, laboratory investigations of the collision properties of micrometer-sized ice particles (in particular, of the most abundant H₂O-ice) have not been conducted yet. Thus, we used two different experimental methods to produce micrometer-sized H₂O-ice particles, i.e. by spraying H₂O droplets into liquid nitrogen and by spraying H₂O droplets into a cold nitrogen atmosphere. The mean particle radii of the ice particles produced with these experimental methods are (1.49 ± 0.79) μm and (1.45 ± 0.65) μm. Ice aggregates composed of the micrometer-sized ice particles are highly porous (volume filling factor: ? = 0.11 ± 0.01) or rather compact (volume filling factor: ? = 0.72 ± 0.04), depending on the method of production. Furthermore, the critical rolling friction force of F_Roll,ice = (114.8 ± 23.8) × 10⁻¹⁰ N was measured for micrometer-sized ice particles, which exceeds the critical rolling friction force of micrometer-sized SiO₂ particles . This result implies that the adhesive bonding between micrometer-sized ice particles is stronger than the bonding strength between SiO₂ particles. An estimation of the specific surface energy of micrometer-sized ice particles, derived from the measured critical rolling friction forces and the surface energy of micrometer-sized SiO₂ particles, results in γ_ice = 0.190 J m⁻².     

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

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

The distribution of water ice in the interior of Comet Tempel 1

The Deep Impact flyby spacecraft includes a 1.05 to 4.8 μm infrared (IR) spectrometer. Although ice was not observed on the surface in the impact region, strong absorptions near 3 μm due to water ice are detected in IR measurements of the ejecta from the impact event. Absorptions from water ice occur throughout the IR dataset beginning three seconds after impact through the end of observations, ∼45 min after impact. Spatially and temporally resolved IR spectra of the ejecta are analyzed in conjunction with laboratory impact experiments. The results imply an internal stratigraphy for Tempel 1 consisting of devolatilized materials transitioning to unaltered components at a depth of approximately one meter. At greater depths, which are thermally isolated from the surface, water ice is present. Up to depths of 10 to 20 m, the maximum depths excavated by the impact, these pristine materials consist of very fine grained (∼1±1 μm) water ice particles, which are free from refractory impurities.     

Asteroid age distributions determined by space weathering and collisional evolution models

We provide evidence of consistency between the dynamical evolution of main belt asteroids and their color evolution due to space weathering. The dynamical age of an asteroid’s surface (Bottke, W.F., Durda, D.D., Nesvorný, D., Jedicke, R., Morbidelli, A., Vokrouhlický, D., Levison, H. [2005]. Icarus 175 (1), 111-140; Nesvorný, D., Jedicke, R., Whiteley, R.J., Ivezi?, ?. [2005]. Icarus 173, 132-152) is the time since its last catastrophic disruption event which is a function of the object’s diameter. The age of an S-complex asteroid’s surface may also be determined from its color using a space weathering model (e.g. Willman, M., Jedicke, R., Moskovitz, N., Nesvorný, D., Vokrouhlický, D., Mothé-Diniz, T. [2010]. Icarus 208, 758-772; Jedicke, R., Nesvorný, D., Whiteley, R.J., Ivezi?, ?., Juri?, M. [2004]. Nature 429, 275-277; Willman, M., Jedicke, R., Nesvorny, D., Moskovitz, N., Ivezi?, ?., Fevig, R. [2008]. Icarus 195, 663-673. We used a sample of 95 S-complex asteroids from SMASS and obtained their absolute magnitudes and u, g, r, i, z filter magnitudes from SDSS. The absolute magnitudes yield a size-derived age distribution. The u, g, r, i, z filter magnitudes lead to the principal component color which yields a color-derived age distribution by inverting our color-age relationship, an enhanced version of the ‘dual τ’ space weathering model of Willman et al. (2010).We fit the size-age distribution to the enhanced dual τ model and found characteristic weathering and gardening times of τ_w = 2050 ± 80 Myr and respectively. The fit also suggests an initial principal component color of −0.05 ± 0.01 for fresh asteroid surface with a maximum possible change of the probable color due to weathering of ΔPC = 1.34 ± 0.04. Our predicted color of fresh asteroid surface matches the color of fresh ordinary chondritic surface of PC₁ = 0.17 ± 0.39.     

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

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

Investigations of the pre-Deep Impact morphology of the dust coma of Comet Tempel-1

The pre-Deep Impact images of Comet Tempel-1 obtained at the Indian Astronomical Observatory are used to investigate the morphology of the dust coma of the comet. We show that the trajectory of a cometary grain under the influence of solar radiation pressure is a reliable diagnostic to estimate its initial velocity. Four main active regions at mean latitudes +45° ± 5°(D), 0° ± 5° (E),−30° ± 5°(A) and−60° ± 5°(F) are found to explain the morphology of the dust coma in the ground-based and published images obtained by the High Resolution Instrument(HRI) cameras aboard the Deep Impact flyby spacecraft. From a χ² fit of the intensity distribution in the observed and the simulated images, we derive the fraction of the productivity of the active vents to the total dust emission of the comet to be 27%. Of this the southern source alone accounts for 19.8%. The grains are found to be ejected with a velocity distribution with an upper limit of 70 ± 7 m s⁻¹. However, the broad region ‘A’ appears to eject slower grains with an upper limit of 24 ± 2.5 m s⁻¹. This source, that is active throughout the cycle is likely to be driven by CO₂ sublimation. We compute the dependence of the percentage contribution of the southern source on the heliocentric distance and show that this ratio varies over the apparition and reaches a maximum at around 260 days before perihelion. The published images of the nucleus of Comet Tempel-1 show significant departure from sphericity. Therefore, the torque exerted by the enhanced activity of the southern region may be significant enough to produce changes in the rotational state of the nucleus before each perihelion passage.     

Gemini-N mid-IR observations of the dust properties of the ejecta excavated from Comet 9P/Tempel 1 during Deep Impact

We present mid-infrared spectra and images from the Gemini-N (+MICHELLE) observational campaign of Comet 9P/Tempel 1 before, during, and after its encounter with Deep Impact. We use our thermal grain model to probe the 10 μm properties of the dust grains in the coma of the comet. Before impact (3 July 2005 UT), and more than 24 h after impact (5, 16, and 28 July 2005 UT), the comet dust grains were composed mostly of amorphous olivine, and were relatively large (peak of the grain size distribution ). For the night of impact, we extract spectra by centering on the nucleus, and offset 1″ from the nucleus in the direction of the impact ejecta plume. We find small dust grains (∼0.2 μm) of a diverse mineralogy (amorphous olivine, amorphous pyroxene, amorphous carbon, and crystalline olivine) populating the ejecta. The submicron sized dust grains move faster than the other, larger grains (?0.7 μm), with amorphous olivine and amorphous carbon traveling together, and amorphous pyroxene and crystalline olivine dispersing at a similar rate. Deriving a velocity law from a time-of-flight analysis, we find that the material traveled with a velocity law scaled by and with a power of p=0.5. This velocity power-law requires a sustained release of grains for the duration of 45-60 min after impact. Since the mineral species are traveling at different speeds, and there was a sustained release of grains due to a possible “gas-plume,” we conclude that the different minerals did not originate from grain aggregates destroyed by the impact, but instead arise from an inhomogeneous nucleus.     

Mars magnetic field: Sources and models for a quarter of the southern hemisphere

The heavily-cratered southern hemisphere of Mars encompasses the planet’s strongest, most widespread magnetization. Our study concentrates on this magnetized region in the southern hemisphere within 40° of latitude 40°S, longitude 180°W. First we rotate the coordinates to position the center at −40°, 180° and treat these new latitudes and longitudes as if they were Cartesian coordinates. Then, using an ordinary two-dimensional Fourier analysis for downward continuation, the MGS (MAG/ER) magnetic field data at satellite mapping elevation of ∼400 km are extrapolated to 100 km, sources are estimated and used to model the fields. Quantitative comparison of the downward continued field with the aerobraking field for bins having angular deviation within ±30° gives correlation of .947, .868, and .769 for the components, respectively. This agreement of the fields may result from most of the power in the magnetization resting in wavelengths ∼400 km, with comparatively little at ∼100 km. Over this region, covering nearly an octant of the planet, just a dozen sources can account for 94% of the variance of the magnetic field at the surface. In these models for the field an obvious asymmetry in polarity exists, with majority of the sources being positive. The locations of strongest surface magnetization appear to be near - but not actually within - ancient multi-ringed basins. We test the likelihood of this association by comparing the observed sources found within and near basins for two alternative basin location scenarios with random distributions. For both alternatives we find the observed distributions to be low-probability occurrences. If contemporaneous, this would establish that Mars’ magnetic field extended to the time of impacts causing these basins.     

Exchange of ejecta between Telesto and Calypso: Tadpoles, horseshoes, and passing orbits

We have numerically integrated the orbits of ejecta from Telesto and Calypso, the two small Trojan companions of Saturn’s major satellite Tethys. Ejecta were launched with speeds comparable to or exceeding their parent’s escape velocity, consistent with impacts into regolith surfaces. We find that the fates of ejecta fall into several distinct categories, depending on both the speed and direction of launch.The slowest ejecta follow suborbital trajectories and re-impact their source moon in less than one day. Slightly faster debris barely escape their parent’s Hill sphere and are confined to tadpole orbits, librating about Tethys’ triangular Lagrange points L₄ (leading, near Telesto) or L₅ (trailing, near Calypso) with nearly the same orbital semi-major axis as Tethys, Telesto, and Calypso. These ejecta too eventually re-impact their source moon, but with a median lifetime of a few dozen years. Those which re-impact within the first 10 years or so have lifetimes near integer multiples of 348.6 days (half the tadpole period).Still faster debris with azimuthal velocity components ?10 m/s enter horseshoe orbits which enclose both L₄ and L₅ as well as L₃, but which avoid Tethys and its Hill sphere. These ejecta impact either Telesto or Calypso at comparable rates, with median lifetimes of several thousand years. However, they cannot reach Tethys itself; only the fastest ejecta, with azimuthal velocities ?40 m/s, achieve “passing orbits” which are able to encounter Tethys. Tethys accretes most of these ejecta within several years, but some 1% of them are scattered either inward to hit Enceladus or outward to strike Dione, over timescales on the order of a few hundred years.     

Using the youngest asteroid clusters to constrain the space weathering and gardening rate on S-complex asteroids

We have extended our earlier work on space weathering of the youngest S-complex asteroid families to include results from asteroid clusters with ages <10⁶ years and to newly identified asteroid pairs with ages <5 × 10⁵ years. We have identified three S-complex asteroid clusters amongst the set of clusters with ages in the range 10^5-6 years—(1270) Datura, (21509) Lucascavin and (16598) 1992 YC2. The average color of the objects in these clusters agrees with the color predicted by the space weathering model of Willman et al. (Willman, M., Jedicke, R., Nesvorný, D., Moskovitz, N., Ivezi?, Z., Fevig, R. [2008]. Icarus 195, 663-673). SDSS five-filter photometry of the members of the very young asteroid pairs with ages <10⁵ years was used to determine their taxonomic classification. Their types are consistent with the background population near each object. The average color of the S-complex pairs is PC₁ = 0.49 ± 0.03, over 5σ redder than predicted by Willman et al. (Willman, M., Jedicke, R., Nesvorný, D., Moskovitz, N., Ivezi?, Z., Fevig, R. [2008]. Icarus 195, 663-673). This may indicate that the most likely pair formation mechanism is a gentle separation due to YORP spin-up leaving much of the aged and reddened surface undisturbed. If this is the case then our color measurement allows us to set an upper limit of ∼64% on the amount of surface disturbed in the separation process. Using pre-existing color data and our new results for the youngest S-complex asteroid clusters we have extended our space weather model to explicitly include the effects of regolith gardening and fit separate weathering and gardening characteristic time scales of τ_w = 960 ± 160 Myr and τ_g = 2000 ± 290 Myr respectively. The first principal component color for fresh S-complex material is PC₁ = 0.37 ± 0.01 while the maximum amount of local reddening is ΔPC₁ = 0.33 ± 0.06. Our first-ever determination of the gardening time is in stark contrast to our calculated gardening time of τ_g ∼ 270 Myr based on main belt impact rates and reasonable assumptions about crater and ejecta blanket sizes. A possible resolution for the discrepancy is through a ‘honeycomb’ mechanism in which the surface regolith structure absorbs small impactors without producing significant ejecta. This mechanism could also account for the paucity of small craters on (433) Eros.     

Determination of a precise rotation period for the Deep Space 1 target, Comet 19P/Borrelly

The rotation period derived by Mueller and Samarasinha (Mueller, B.E.A., Samarasinha, N.H. [2002]. Earth Moon Planets 90, 463-471) of the Deep Space 1 (DS1) mission target, Comet 19P/Borrelly, using ground-based data from July 28 to August 1, 2000, is improved by two orders of magnitude. This precision is reached in a multistep process.Combining all available ground-based data in 2000 decreases the error by an order of magnitude. Next, assuming that the rotation period did not change between 2000 and 2001, constraints from the HST 2001 data (Weaver, H.A., Stern, S.A., Parker, J.Wm. [2003]. Astron. J. 126, 444-451) yield three possible rotation periods: P = 1.088 ± 0.003 days, P = 1.108 ± 0.002 days, and P = 1.135 ± 0.003 days, which are consistent with our initial derivation of P = 1.08 ± 0.04 days (Mueller, B.E.A., Samarasinha, N.H. [2002]. Earth Moon Planets 90, 463-471).These three periods are further refined and the error bars further improved by another order of magnitude by linking the combined ground-based data from 2000 to the nuclear orientation of Borrelly at the DS1 encounter in 2001 (see Table 2). Due to aliasing, there are seven possible rotation periods around P = 1.088 days, five possible periods around P = 1.108 days, and six possible periods around P = 1.135 days, with precisions of the order of 0.0002 days (≈17 s).     

Photometric analysis and disk-resolved thermal modeling of Comet 19P/Borrelly from Deep Space 1 data

Images returned from the Deep Space 1 (DS-1) spacecraft during its encounter with Comet 19P/Borrelly are used to study its disk-integrated and disk-resolved photometry and its thermal properties. A disk-integrated phase function was constructed from a combination of DS-1 images and ground-based observations, giving a geometric albedo of 0.072±0.020 and a phase slope of 0.043 mag deg⁻¹. The shape model of Borrelly [Kirk, R.L., Howington-Kraus, E., Soderblom, L.A., Giese, B., Oberst, J., 2004a. Icarus 167, 54-69] and the ephemerides of DS-1 were used to analyze the disk-resolved photometric data with Hapke's theoretical model. It was found that the surface of Borrelly displays large photometric heterogeneities in its photometric parameters. The single-scattering albedo, w, varies by a factor of 2.5 with an average of 0.057±0.009; the asymmetry factor, g, ranges from almost isotropic (−0.1) to strongly backscattering (−0.7) with an average of −0.43±0.07; the roughness parameter, , is less than 35° for most parts of surface but ranges up to 55° in some areas. Its average is 22°±5°. The observed 1-D temperature profile is modeled well by the standard thermal model (STM) for inactive regions and is found to be consistent with a very low thermal inertia. Water sublimation in the source region of the fan jet is observed to decrease the surface temperature from the STM predictions by 20-40 K. The source areas of two collimated jets could not be determined from either photometric model or thermal model. It is evident that the fan jet activity occurring on Borrelly's surface can be correlated to areas of relatively high albedo, weak backscattering, and high roughness.