Spitzer Space Telescope Observations of the Nucleus of Comet 103P/Hartley 2

We have used the Spitzer Space Telescope InfraRed Spectrograph (IRS) 22-μm peakup array to observe thermal emission from the nucleus and trail of comet 103P/Hartley 2, the target of NASA’s Deep Impact Extended Investigation (DIXI). The comet was observed on UT 2008 August 12 and 13, while 5.5 AU from the Sun. We obtained two 200 frame sets of photometric imaging over a 2.7 hr period. To within the errors of the measurement, we find no detection of any temporal variation between the two images. The comet showed extended emission beyond a point source in the form of a faint trail directed along the comet’s antivelocity vector. After modeling and removing the trail emission, a NEATM model for the nuclear emission with beaming parameter of 0.95 ± 0.20 indicates a small effective radius for the nucleus of 0.57 ± 0.08 km and low geometric albedo 0.028 ± 0.009 (1σ). With this nucleus size and a water production rate of 3 × 10^(28) molecules s^(-1) at perihelion, we estimate that ~100% of the surface area is actively emitting volatile material at perihelion. Reports of emission activity out to ~5 AU support our finding of a highly active nuclear surface. Compared to Deep Impact’s first target, comet 9P/Tempel 1, Hartley 2’s nucleus is one-fifth as wide (and about one-hundredth the mass) while producing a similar amount of outgassing at perihelion with about 13 times the active surface fraction. Unlike Tempel 1, comet Hartley 2 should be highly susceptible to jet driven spin-up torques, and so could be rotating at a much higher frequency. Since the amplitude of nongravitational forces are surprisingly similar for both comets, close to the ensemble average for ecliptic comets, we conclude that comet Hartley 2 must have a much more isotropic pattern of time-averaged outgassing from its nuclear surface. Barring a catastrophic breakup or major fragmentation event, the comet should be able to survive up to another 100 apparitions (~700 yr) at its current rate of mass loss

Water Ice and Dust in the Innermost Coma of Comet 103P/Hartley 2

On November 4th, 2010, the Deep Impact eXtended Investigation (DIXI) successfully encountered comet 103P/Hartley 2, when it was at a heliocentric distance of 1.06 AU. Spatially resolved near-IR spectra of comet Hartley 2 were acquired in the 1.05-4.83 micron wavelength range using the HRI-IR spectrometer. We present spectral maps of the inner ~10 kilometers of the coma collected 7 minutes and 23 minutes after closest approach. The extracted reflectance spectra include well-defined absorption bands near 1.5, 2.0, and 3.0 micron consistent in position, bandwidth, and shape with the presence of water ice grains. Using Hapke's radiative transfer model, we characterize the type of mixing (areal vs. intimate), relative abundance, grain size, and spatial distribution of water ice and refractories. Our modeling suggests that the dust, which dominates the innermost coma of Hartley 2 and is at a temperature of 300K, is thermally and physically decoupled from the fine-grained water ice particles, which are on the order of 1 micron in size. The strong correlation between the water ice, dust, and CO2 spatial distribution supports the concept that CO2 gas drags the water ice and dust grains from the nucleus. Once in the coma, the water ice begins subliming while the dust is in a constant outflow. The derived water ice scale-length is compatible with the lifetimes expected for 1-micron pure water ice grains at 1 AU, if velocities are near 0.5 m/s. Such velocities, about three order of magnitudes lower than the expansion velocities expected for isolated 1-micron water ice particles [Hanner, 1981; Whipple, 1951], suggest that the observed water ice grains are likely aggregates.Comment: 51 pages, 12 figures, accepted for publication in Icaru

Spitzer observations of the asteroid-comet transition object and potential spacecraft target 107P (4015) Wilson-Harrington

Context. Near-Earth asteroid-comet transition object 107P/ (4015) Wilson-Harrington is a possible target of the joint European Space Agency (ESA) and Japanese Aerospace Exploration Agency (JAXA) Marco Polo sample return mission. Physical studies of this object are relevant to this mission, and also to understanding its asteroidal or cometary nature. Aims. Our aim is to obtain significant new constraints on the surface thermal properties of this object. Methods. We present mid-infrared photometry in two filters (16 and 22 microns) obtained with NASA's Spitzer Space Telescope on February 12, 2007, and results from the application of the Near Earth Asteroid Thermal Model (NEATM).We obtained high S/N in two mid-IR bands allowing accurate measurements of its thermal emission. Results. We obtain a well constrained beaming parameter (eta = 1.39 +/- 0.26) and obtain a diameter and geometric albedo of D = 3.46 +/- 0.32 km, and pV = 0.059 +/- 0.011. We also obtain similar results when we apply this best-fitting thermal model to single-band mid-IR photometry reported by Campins et al. (1995), Kraemer et al. (2005) and Reach et al. (2007). Conclusions. The albedo of 4015 Wilson-Harrington is low, consistent with those of comet nuclei and primitive C-, P-, D-type asteorids. We establish a rough lower limit for the thermal inertia of W-H of 60 Jm^-2s^(-0.5)K^-1 when it is at r=1AU, which is slightly over the limit of 30 Jm^-2s^(-0.5)K-1 derived by Groussin et al. (2009) for the thermal inertia of the nucleus of comet 22P/Kopff.Comment: 4 pages, 1 figure and 3 tables. Paper accepted for publicatio

A Branch-Heterogeneous Model of Protein Evolution for Efficient Inference of Ancestral Sequences

International audienceMost models of nucleotide or amino acid substitution used in phylogenetic studies assume that the evolutionary process has been homogeneous across lineages and that composition of nucleotides or amino acids has remained the same throughout the tree. These oversimplified assumptions are refuted by the observation that compositional variability characterizes extant biological sequences. Branch-heterogeneous models of protein evolution that account for compositional variability have been developed, but are not yet in common use because of the large number of parameters required, leading to high computational costs and potential overparameterization. Here, we present a new branch-nonhomogeneous and nonstationary model of protein evolution that captures more accurately the high complexity of sequence evolution. This model, henceforth called Correspondence and likelihood analysis (COaLA), makes use of a correspondence analysis to reduce the number of parameters to be optimized through maximum likelihood, focusing on most of the compositional variation observed in the data. The model was thoroughly tested on both simulated and biological data sets to show its high performance in terms of data fitting and CPU time. COaLA efficiently estimates ancestral amino acid frequencies and sequences, making it relevant for studies aiming at reconstructing and resurrecting ancestral amino acid sequences. Finally, we applied COaLA on a concatenate of universal amino acid sequences to confirm previous results obtained with a nonhomogeneous Bayesian model regarding the early pattern of adaptation to optimal growth temperature, supporting the mesophilic nature of the Last Universal Common Ancestor

Simulations of the Population of Centaurs I: The Bulk Statistics

Large-scale simulations of the Centaur population are carried out. The evolution of 23328 particles based on the orbits of 32 well-known Centaurs is followed for up to 3 Myr in the forward and backward direction under the influence of the 4 massive planets. The objects exhibit a rich variety of dynamical behaviour with half-lives ranging from 540 kyr (1996 AR20) to 32 Myr (2000 FZ53). The mean half-life of the entire sample of Centaurs is 2.7 Myr. The data are analyzed using a classification scheme based on the controlling planets at perihelion and aphelion, previously given in Horner et al (2003). Transfer probabilities are computed and show the main dynamical pathways of the Centaur population. The total number of Centaurs with diameters larger than 1 km is estimated as roughly 44300, assuming an inward flux of one new short-period comet every 200 yrs. The flux into the Centaur region from the Edgeworth-Kuiper belt is estimated to be 1 new object every 125 yrs. Finally, the flux from the Centaur region to Earth-crossing orbits is 1 new Earth-crosser every 880 yrsComment: 15 pages, 2 figures, MNRAS in pres

Activity distribution of comet 67P/Churyumov-Gerasimenko from combined measurements of non-gravitational forces and torques

Aims. Understanding the activity is vital for deciphering the structure, formation, and evolution of comets. We investigate models of cometary activity by comparing them to the dynamics of 67P/Churyumov-Gerasimenko. Methods. We matched simple thermal models of water activity to the combined Rosetta datasets by fitting to the total outgassing rate and four components of the outgassing induced non-gravitational force and torque, with a final manual adjustment of the model parameters to additionally match the other two torque components. We parametrised the thermal model in terms of a distribution of relative activity over the surface of the comet, and attempted to link this to different terrain types. We also tested a more advanced thermal model based on a pebble structure. Results. We confirm a hemispherical dichotomy and non-linear water outgassing response to insolation. The southern hemisphere of the comet and consolidated terrain show enhanced activity relative to the northern hemisphere and dust-covered, unconsolidated terrain types, especially at perihelion. We further find that the non-gravitational torque is especially sensitive to the activity distribution, and to fit the pole-axis orientation in particular, activity must be concentrated (in excess of the already high activity in the southern hemisphere and consolidated terrain) around the south pole and on the body and neck of the comet over its head. This is the case for both the simple thermal model and the pebble-based model. Overall, our results show that water activity cannot be matched by a simple model of sublimating surface ice driven by the insolation alone, regardless of the surface distribution, and that both local spatial and temporal variations are needed to fit the data. Conclusions. Fully reconciling the Rosetta outgassing, torque, and acceleration data requires a thermal model that includes both diurnal and seasonal effects and also structure with depth (dust layers or ice within pebbles). This shows that cometary activity is complex. Nonetheless, non-gravitational dynamics provides a useful tool for distinguishing between different thermophysical models and aids our understanding

Narrow Dust Jets in a Diffuse Gas Coma: A Natural Product of Small Active Regions on Comets

Comets often display narrow dust jets but more diffuse gas comae when their eccentric orbits bring them into the inner solar system and sunlight sublimates the ice on the nucleus. Comets are also understood to have one or more active areas covering only a fraction of the total surface active with sublimating volatile ices. Calculations of the gas and dust distribution from a small active area on a comet's nucleus show that as the gas moves out radially into the vacuum of space it expands tangentially, filling much of the hemisphere centered on the active region. The dust dragged by the gas remains more concentrated over the active area. This explains some puzzling appearances of comets having collimated dust jets but more diffuse gaseous atmospheres. Our test case is 67P/Churyumov-Gerasimenko, the Rosetta mission target comet, whose activity is dominated by a single area covering only 4% of its surface.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/98552/1/0004-637X_749_1_29.pd

Exploring the surface properties of Transneptunian Objects and Centaurs with polarimetric FORS1/VLT observations

Polarization is a powerful remote-sensing method to investigate solar system bodies. It is an especially sensitive diagnostic tool to reveal physical properties of the bodies whose observational characteristics are governed by small scatterers (dust, regolith surfaces). For these objects, at small phase angles, a negative polarization is observed, i.e., the electric vector E oscillates predominantly in the scattering plane, contrary to what is typical for rather smooth homogeneous surfaces. The behavior of negative polarization with phase angle depends on the size, composition and packing of the scatterers. These characteristics can be unveiled by modelling the light scattering by the dust or regolith in terms of the coherent backscattering mechanism. We have investigated the surface properties of TNOs and Centaurs by means of polarimetric observations with FORS1 of the ESO VLT. TNOs Ixion and Quaoar, and Centaur Chiron show a negative polarization surge. The Centaur Chiron has the deepest polarization minimum (-1.5 - 1.4%). The two TNOs show differing polarization curves: for Ixion, the negative polarization increases rapidly with phase; for Quaoar, the polarization is relatively small (~ -0.6%), and nearly constant at the observed phase angles. For all three objects, modelling results suggest that the surface contains an areal mixture of at least two components with different single-scatterer albedos and photon mean-free paths.Comment: 11 pages, 7 postscript figures, accepted by A&A; astro-ph abstract has been replaced with a more complete on

The Highly Unusual Outgassing of Comet 103P/Hartley 2 from Narrowband Photometry and Imaging of the Coma

We report on photometry and imaging of Comet 103P/Hartley 2 obtained at Lowell Observatory from 1991 through 2011. We acquired photoelectric photometry on two nights in 1991, four nights in 1997/98, and 13 nights in 2010/11. We observed a strong secular decrease in water and all other observed species production in 2010/11 from the 1991 and 1997/98 levels. We see evidence for a strong asymmetry with respect to perihelion in the production rates of our usual bandpasses, with peak production occurring ~10 days post-perihelion and production rates considerably higher post-perihelion. The composition was "typical", in agreement with the findings of other investigators. We obtained imaging on 39 nights from 2010 July until 2011 January. We find that, after accounting for their varying parentage and lifetimes, the C2 and C3 coma morphology resemble the CN morphology we reported previously. These species exhibited an hourglass shape in October and November, and the morphology changed with rotation and evolved over time. The OH and NH coma morphology showed hints of an hourglass shape near the nucleus, but was also enhanced in the anti-sunward hemisphere. This tailward brightness enhancement did not vary significantly with rotation and evolved with the viewing geometry. We conclude that all five gas species likely originate from the same source regions on the nucleus, but that OH and NH were derived from small grains of water and ammonia ice that survived long enough to be affected by radiation pressure and driven in the anti-sunward direction. We detected the faint, sunward facing dust jet reported by other authors, and did not detect a corresponding gas feature. This jet varied little during a night but exhibited some variations from night to night, suggesting it is located near the total angular momentum vector.Comment: Accepted by Icarus; 20 pages of text (preprint style), 5 tables, 7 figure

Shape modeling technique KOALA validated by ESA Rosetta at (21) Lutetia

We present a comparison of our results from ground-based observations of asteroid (21) Lutetia with imaging data acquired during the flyby of the asteroid by the ESA Rosetta mission. This flyby provided a unique opportunity to evaluate and calibrate our method of determination of size, 3-D shape, and spin of an asteroid from ground-based observations. We present our 3-D shape-modeling technique KOALA which is based on multi-dataset inversion. We compare the results we obtained with KOALA, prior to the flyby, on asteroid (21) Lutetia with the high-spatial resolution images of the asteroid taken with the OSIRIS camera on-board the ESA Rosetta spacecraft, during its encounter with Lutetia. The spin axis determined with KOALA was found to be accurate to within two degrees, while the KOALA diameter determinations were within 2% of the Rosetta-derived values. The 3-D shape of the KOALA model is also confirmed by the spectacular visual agreement between both 3-D shape models (KOALA pre- and OSIRIS post-flyby). We found a typical deviation of only 2 km at local scales between the profiles from KOALA predictions and OSIRIS images, resulting in a volume uncertainty provided by KOALA better than 10%. Radiometric techniques for the interpretation of thermal infrared data also benefit greatly from the KOALA shape model: the absolute size and geometric albedo can be derived with high accuracy, and thermal properties, for example the thermal inertia, can be determined unambiguously. We consider this to be a validation of the KOALA method. Because space exploration will remain limited to only a few objects, KOALA stands as a powerful technique to study a much larger set of small bodies using Earth-based observations.Comment: 15 pages, 8 figures, 2 tables, accepted for publication in P&S