VRI imaging of comet 46P/Wirtanen

Imaging of comet 46P/Wirtanen was performed in the standard VRI filters on 10 and 11 November 1996 with the 1.0 m telescope of the South African Astronomical Observatory (SAAO). After proper processing and calibration, the images have been used to derive information about the isophote distribution, magnitude and colour indices. We observed a coma extending to at most ∼5×10³ km. We derive absolute magnitudes of the coma for the different filters and the colours V−R = 0.18±0.17 and R−I = 0.39±0.16. The Afρ value (related to the dust production rate) for the R filter is 10.5±0.6 cm. The radial profile (azimuthally averaged) of the surface brightness in the R filter has a slope of −0.9, which is consistent with a steady-state dust production. The results are compared with other observations in different periods.     

The long-term evolution and initial size of comets 46P/Wirtanen and 67P/Churyumov–Gerasimenko

We present a new method to study the long-term evolution of cometary nuclei in order to estimate their original size, and we consider the case of comets 46P/Wirtanen (hereafter 46P) and 67P/Churyumov–Gerasimenko (hereafter 67P). We calculate the past evolution of the orbital elements of both comets over 100 000 yr using a Bulirsch–Stoer integrator and over 450 000 yr using a Radau integrator, and we incorporate a realistic model of the erosion of their nucleus. Their long-term orbital evolution is prominently chaotic, resulting from several close encounters with planets, and this result is independent of the choice of the integrator and of the presence or not of non-gravitational forces. The dynamical lifetime of comet 46P is estimated at ∼133 000 yr and that of comet 67P at ∼105 000 yr. Our erosion model assumes a spherical nucleus composed of a macroscopic mixture of two thermally decoupled components, dust and pure water ice. Erosion strongly depends upon the active fraction and the density of the nucleus. It mainly takes place at heliocentric distances <4 au and lasts for only ∼7 per cent of the lifetime. Assuming a density of 300 kg m⁻³ and an average active fraction over time of 10 per cent, we find an initial radius of ∼1.3 km for 46P and ∼2.8 km for 67P. Upper limit are obtained assuming a density of 100 kg m⁻³ and an active fraction of 100 per cent, and amounts to 21 km for 46P and 25 km for 67P. Erosion acts as a rejuvenating process of the surface so that exposed materials on the surface may only contain very little quantities of primordial materials. However, materials located just under it (a few centimetres to metres) may still be much less evolved. We will apply this method to several other comets in the future.     

Energy balance and the gas flux from the surface of comet 46P/Wirtanen

Understanding the power balance at the surface of the nucleus is essential to study the chemical and physical evolution of a comet. Therefore, we present a detailed energy budget analysis for the surface of a model comet in the orbit of 46P/Wirtanen, target comet of the European space craft mission Rosetta, for a variety of parameters and assumptions. We will show that for a fast spinning Jupiter-family comet such as 46P/Wirtanen with a rotation period of about 6 h, a fast rotator approximation underestimates the effective energy input. This yields lower gas fluxes from the surface. For an 100% active, non-dust covered surface we obtain a water gas flux on the order of about 1.5×10²⁸ molecules s⁻¹ at perihelion, assuming a radius of 600 m. The calculated gas flux of water is within the order of measured values for comet 46P/Wirtanen. But our calculated values are maximum gas fluxes at noon—not averaged over one cometary day or taking the lesser insolation at the polar areas into account. Therefore, we conclude that either the radius of comet 46P/Wirtanen may be much larger than the accepted value of 600 m. A radius in the order of 2 km seems more likely to explain the measurements. Or, an other possibility could be that water-ice particles are blown off from the surface like dust particles. This may also increase the effective surface area of sublimation.     

Constraints on comet 46P/Wirtanen dust parameters provided by in-situ and ground-based observations

A dust environment working group was encouraged by ESA to provide coma dust environment models useful to plan the ROSETTA operations around the nucleus of short period comet 46P/Wirtanen. Among the many parameters describing the dust released from the nucleus surface, special care was devoted to the dust size distribution. Its present uncertainty makes all environment models sensitive, mainly, to which actual size distribution is adopted. In fact, it must be stressed that no other cometary dust parameter can be derived, such as dust loss rate or dust to gas ratio, if the size distribution remains undetermined. This paper will focus, therefore, on the available information on cometary dust size distributions, starting from the in situ measurement cornerstone provided by the GIOTTO-DIDSY results. Available ground-based observations are then reviewed, in order to disentangle the real sensitivity of them to this quantity; the size distribution is always embedded together with other dust parameters, and its influence on the published results is often forgotten.     

Lyman-alpha observations of comet 46 P/Wirtanen with swan on SOHO: H2O production rate near 1997 perihelion

The SWAN instrument on board SOHO is a Lyman-α photometer able to map the sky intensity with a resolution of 1°, primarily devoted to the study of the large scale distribution of solar wind from its imprints on the interplanetary sky background. In addition SWAN was extensively used to map the Lyman α emission of several comets since launch in December 1995. Here we report observations of Comet 46 P/Wirtanen near perihelion. From the recorded Lyman α intensity the H₂O production rate was derived for 45 observations from 21 December 1996–17 May 1997, with a peak of 1.6±0.4×10²⁸ mol/s just before perihelion. This should help to constrain the physical models of 46 P/Wirtanen for Rosetta mission planning purposes.     

The nucleus of Comet P/Arend-Rigaux

Infrared observations, including photometry, spectroscopy, and J, H, and K lightcurves, of Comet P/Arend-Rigaux, suspected of being a near-extinct comet nucleus, are consistent with the detection of a rotating nucleus of radius 5.1 ± 1.1 km and a faint coma. The visual geometric albedo of the nucleus is 0.02 ± 0.01, placing the nucleus among the darkest objects in the Solar System. The near-infrared spectrum exhibits a red slope with no clear evidence of ice. A nonvolatile dust mantle compositionally similar to D-type asteroids may explain the observations.     

The gas production of Comet 9P/Tempel 1 around the Deep Impact date

The target of the Deep Impact space mission (NASA), Comet 9P/Tempel 1, was observed from two nights before impact to eight nights after impact using the FORS spectrographs at the ESO VLT UT1 and UT2 telescopes. Low resolution optical long-slit spectra were obtained to study the evolution of the gas coma around the Deep Impact event. Following first results of this observing campaign on the CN and dust activity [Rauer, H., Weiler, M., Sterken, C., Jehin, E., Knollenberg, J., Hainaut, O., 2006. Astron. Astrophys. 459, 257-263], this work presents a study of the complete dataset on CN, C₂, C₃, and NH₂ activity of Comet 9P/Tempel 1. An extended impact gas cloud was observed moving radially outwards. No compositional differences between this impact cloud and the undisturbed coma were found as far as the observed radicals are concerned. The gas production rates before and well after impact indicate no change in the cometary activity on an intermediate time scale. Over the observing period, the activity of Comet 9P/Tempel 1 was found to be related to the rotation of the cometary nucleus. The rotational lightcurve for different gaseous species provides indications for compositional differences among different parts of the nucleus surface.     

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.     

The 'silent world' of Comet 15P/Finlay

Comet 15P/Finlay is unusual in that, contrary to ab initio expectations, it demonstrates no apparent linkage to any known meteor shower. Using data contained within the Electronic Atlas of Dynamical Evolutions of Short-Period Comets, we evaluate theoretical shower radiants for Comet 15P/Finlay, but find no evidence to link it to any meteoric anomalies in recorded antiquity. This result, however, must be tempered by the fact that any Comet 15P/Finlay-derived meteoroids will have a low, 16 km s⁻¹, encounter velocity with Earth's atmosphere. Typically, therefore, one would expect mostly faint meteors to be produced during an encounter with a Comet 15P/Finlay-derived meteoroid stream. We have conducted a D -criterion survey of meteoroid orbits derived from three southern hemisphere meteor radar surveys conducted during the 1960s, and again we find no evidence for any Comet 15P/Finlay-related activity. Numerical calculations following the orbital evolution of hypothetical meteoroids ejected from the comet, at each perihelion epoch since 1886, indicate that Jovian perturbations effectively 'drive' the meteoroids to orbits with nodal points beyond the Earth's orbit. The numerical calculations indicate that, even if Comet 15P/Finlay had been a copious emitter of meteoroids during the past 100 years, virtually none of them would have evolved into orbits capable of being sampled by the Earth. There are good observational data, however, to suggest that Comet 15P/Finlay is becoming a transitional comet–asteroid object, and that it has probably not been an efficient producer of meteoroids during the past several hundreds of years.     

The dust coma of Comet P/Giacobini-Zinner in the infrared

We present 1-20 micrometers photometry of P/Giacobini-Zinner obtained at the NASA Infrared Telescope Facility, during 1985 June-September (r = 1.57-1.03 AU). A broad, weak 10 micrometers silicate emission feature was detected on August 26.6; a similar weak emission feature could have been hidden in the broadband photometry on other dates. The total scattering and emitting cross section of dust in the inner coma was similar to that in other short-period comets, but a factor of 10 (r = 1.56 AU) to 100 (r = 1.03 AU) lower than the amount of dust in Comet Halley. The thermal emission continuum can be fit with models weighted toward either small or large absorbing grains. The dust production rate near perihelion was approximately 10(5) g/s (small-grain model) to approximately 10(6) g/s (large-grain model). The corresponding dust/gas mass ratio on August 26 was approximately 0.1-1. A silicate-rich heterogeneous grain model with an excess of large particles is compatible with the observed spectrum of Giacobini-Zinner on August 26. Thus, weak or absent silicate emission does not necessarily imply an absence of silicates in the dust, although the abundance of silicate particles < or = 1 micrometer radius must have been lower than in Comet Halley.     

The effect of internal inhomogeneity on the activity of comet nuclei - Application to Comet 67P/Churyumov-Gerasimenko

We present thermal evolution calculations of inhomogeneous asymmetric initial configurations of a spherical model of Comet 67P/Churyumov-Gerasimenko, using a fully 3-dimensional numerical code. The initial composition is amorphous H₂O ice and dust, in a “layered-pile” configuration, where layers differing in ice/dust ratio and thermal properties extend over a fraction of the surface area and about 10 m in depth and may overlap. We analyze the effect of one such layer, as well as the combined effect of many layers, randomly distributed. We find that internal inhomogeneities affect both the surface temperature and the activity pattern of the comet. In particular, they may lead to outbursts at large heliocentric distances and also to activity on the night-side of the nucleus. The rates of ablation and depths of dust mantle and crystalline ice outer layer as functions of longitude and latitude are shown to be affected as well.     

The excited spin state of Comet 2P/Encke

Ways to rationalize the different periods (e.g., 15.08 h, Luu and Jewitt, 1990, Icarus 86, 69-81; 11.01 h, Fernández et al., 2004, Icarus, in this issue; Lowry et al., 2003, Lunar Planet. Sci. XXXIV, Abstract 2056) seen in near aphelion R-band light curves of Comet 2P/Encke are explored. We show that the comet is usually active at aphelion and it's observed light curves contain signal from both the nucleus and an unresolved coma. The coma contribution to the observed brightness is generally found to dominate with the nucleus providing from 28 to 87% of the total brightness. The amplitude of the observed variations cannot be explained by the nucleus alone and are due to coma activity. We show that some seven periodicities exist in the observed light curves at various times and that this is likely the result of an active nucleus spinning in an excited spin state. The changing periodicities are probably due to changes in the relative strengths of the active areas. We work out possible excited states based on experience with model light curves and by using an analogy to light curve observations of Comet 1P/Halley for which the spin state has been separately determined from spacecraft observations. There is a possibility of a fully relaxed principal axis spin state (0.538 d⁻¹; P=44.6 h) but, because it provides a poorer fit to the observed periodicities than the best fit excited state together with the absence of a peak near 1.08 d⁻¹ (2fφ) in the frequency spectrum of the Fernández et al. (2000, Icarus 147, 145-160) thermal IR lightcurve, we consider it unlikely. Both SAM and LAM excited states are allowed by the underlying periodicities and additional information is needed to choose between these. Our choice of a low excitation SAM state, i.e., one in which the instantaneous spin axis nutates around the total angular momentum vector in a motion that is characterized by limited angular oscillations around the long axis, is based on Sekanina's (1988, Astron J. 95, 911-924, 1988, Astron. J. 96, 1455-1475) interpretation of the fan coma that this comet often displays. We argue that possible LAM states are excluded either because they are too difficult to excite or because they would be inconsistent with the formation of the observed fan morphology. Two possible SAM states emerge that provide good fits to the observed periodicities, one with a precessional frequency for the long axis about the total angular momentum vector of 1.614 d⁻¹ (P?=14.9 h) and an oscillation frequency around the long axis of 0.539 d⁻¹ (Pψ=44.5 h) and a second with a precessional frequency of 2.162 d⁻¹ (P?=11.1 h) combined with an oscillation around the long axis of 0.502 d⁻¹ (Pψ=47.8 h). While either solution is possible, the latter is, in a least squares sense, more likely to be the actual spin state. In both cases the direction of the total angular momentum vector (αM,δM[J2000]=198.6, −0.3 deg) is assumed to be defined by the evolving geometry and morphology of the coma (Sekanina, 1988, Astron J. 95, 911-924, 1988, Astron. J. 96, 1455-1475; Festou and Barale, 2000, Astron J. 119, 3119-3132). We discuss the possible locations of the primary active areas found by Sekanina (1988, Astron J. 95, 911-924, 1988, Astron. J. 96, 1455-1475) and, while they are at high cometographic latitudes, they do not have to be physically located close the region were the axis of maximum moment of inertia pierces the surface (i.e., at high cometocentric latitude). We offer a new interpretation of the 10.7 μm data by Fernández et al. (2000, Icarus 147, 145-160) which yields an axial ratio a/b=2.04. This, with the two SAM states that we have found, requires that b/c>1.18 or >1.09 implying a significant asymmetry in the shape of the elongated nucleus. For the observed fan morphology to be maintained, the true axial ratio b/c cannot be much larger than these limiting values otherwise the amplitude of the oscillation about the long axis becomes too large and the fan morphology would be destroyed. The precise phasing of the spin modes, i.e., the value of the Euler angles at a particular time, is not determinable from the current data set, but a set of well sampled thermal infrared observations of the nucleus covering many periods and a wide range of observing geometries could provide this information in the future as well as clearly distinguishing between the two excited spin states.     

Modelling of Shape Changes of the Nuclei of Comets C/1995 O1 Hale–Bopp and 46P/Wirtanen Caused by Water Ice Sublimation

The aim of this modelling work is to assess shape changes of cometary nuclei caused by sublimation of ices. The simplest possible model is assumed with the nucleus being initially spherical and its thermal conductivity being neglected. We have calculated the time-dependent sublimation flux versus cometographic latitude. If the rotation axis of the comet is inclined to the orbital plane, then sublimation leads to non-symmetrical changes of the nucleus shape. Calculations were performed for the nuclei of comets Hale–Bopp and Wirtanen.     

Models of P/Wirtanen nucleus: active regions versus non-active regions

From the current understanding we know that comet nuclei have heterogeneous compositions and complex structures. It is believed that cometary activity is the result of a combination of physical processes in the nucleus, like sublimation and recondensation of volatile ices, dust grains release, phase transition of water ice, depletion of the most volatile components in the outer layers and interior differentiation.The evolution of the comet depends on the sublimation of ices and the release of different gases and dust grains: the formation of a dust crust, the surface erosion and the development of the coma are related to the gas fluxes escaping from the nucleus. New observations, laboratory experiments and numerical simulations suggest that the gas and dust emissions are locally generated, in the so-called active regions. This localized activity is probably superimposed to the global nucleus activity. The differences between active and inactive regions can be attributed to differences in texture and refractory material content of the different areas.In this paper we present the results of numerical models of cometary nucleus evolution, developed in order to understand which are the processes leading to the formation of active and non-active regions on the cometary surface. The used numerical code solves the equations of heat transport and gas diffusion within a porous nucleus composed of different ices—such as water (the dominant constituent), CO₂, CO- and of dust grains embedded in the ice matrix.By varying the set of physical parameters describing the initial properties of comet P/Wirtanen, the different behaviour of the icy and dusty areas can be followed.Comet P/Wirtanen is the target of the international ROSETTA mission, the cornerstone ESA mission to a cometary nucleus. The successful design of ROSETTA requires some knowledge of comet status and activity: surface temperatures, amount of active and inactive surface areas, gas production rate and dust flux.     

The morphology and surface processes of Comet 19/P Borrelly

The flyby of the nucleus of the Comet 19P/Borrelly by the Deep Space 1 spacecraft produced the best views to date of the surface of these interesting objects. It transformed Borrelly from an astronomical object shrouded in coma of gas and dust into a geological object with complex surface processes and a rich history of erosion and landform evolution. Based on analysis of the highest resolution images, stereo images, photometry, and albedo we have mapped four major morphological units and four terrain features. The morphological units are named dark spots, mottled terrain, mesas, and smooth terrain. The features are named ridges, troughs, pits, and hills. In strong contrast to asteroids, unambiguous impact craters were not observed on Borrelly's surface. Because of the relatively short period of this comet, surface erosion by volatile sublimation is, in geologic terms, a very active process. The formation and the morphologies of units and features appear to be driven by differential rates of sublimation erosion. Erosional rates across the comet are probably controlled by solar energy input and the location of the subsolar point during perihelion. Differences in energy input may produce different varieties of sublimation erosional landforms. The terrains on Borrelly suggest that solar energy input could map directly into erosional processes and landforms.     

Charge-coupled device photometry of Comet P/Halley

Comet P/Halley has been observed during its approach to perihelion at heliocentric distancesR = 11.0 AU and R = 8.2 AU. No extended coma is seen and limits can be placed on the fraction of the total light contributed by coma. The brightness of the comet varies on a short time scale. The variations may be due to transient activity or to rotation of the irregular nucleus.     

Estimating the nucleus density of Comet 19P/Borrelly

The nucleus bulk density of Comet 19P/Borrelly has been estimated by modeling the sublimation-induced non-gravitational force acting upon the orbital motion, thereby reproducing the empirical perihelion advance (i.e., the shortening of the orbital period). The nucleus has been modeled as a prolate ellipsoid, covered by various surface activity maps which reproduce the observed water production rate. The theoretical water production rate of active areas has been obtained by applying a sophisticated thermophysical model. This model takes into account net sublimation of ice and thermal reradiation from the surface, solid state conductivity, sub-surface sublimation and recondensation, mass and heat transport by diffusing gas, layer absorption of solar energy, a full treatment of local time-dependent illumination conditions, and a detailed consideration of nucleus/coma interaction mechanisms. The outgassing properties of the modeled nucleus are physically consistent with the gas kinetic structure of the innermost coma since the molecular backflux and surface gas density required in the thermophysical model (as functions of the nucleus surface temperature and the sub-surface temperature profile) have been obtained from Direct Simulation Monte Carlo modeling of inelastic intermolecular collisions in the cometary Knudsen layer. The calculation of local normal forces acting on the nucleus due to outgassing has been made within the same framework—recoil and/or impact momentum transfer to the nucleus caused by sublimating molecules and by recondensing and/or scattered coma molecules is therefore evaluated in accordance with local nucleus/coma conditions. According to this model, the density is found to be 100-300 kg, depending on the applied spin axis orientation and surface activity map. This range can be narrowed down to 180-300 kg by also requiring that the empirical changes (per orbital revolution) of the argument of perihelion and the longitude of the ascending node are reproduced.     

Searching for the Deep Impact crater on Comet 9P/Tempel 1 using image processing techniques

In this work we attempt to obtain direct images of the crater associated with the impact of the Deep Impact impactor spacecraft on the nucleus of Comet 9P/Tempel 1 on July 4, 2005. The impact generated a large and bright ejecta cloud that hampers the clear view of the post-impact nucleus surface. We used image restoration techniques to enhance spatial resolution and contrast on a subset of selected post-impact high resolution images. No unambiguous evidence for the crater can be found; however, indirect evidence is consistent with a crater size in the 150-200 m range.