Structure and density of cometary nuclei

Abstract— Understanding the nature of the cometary nucleus remains one of the major problems in solar system science. Whipple's (1950) icy conglomerate model has been very successful at explaining a range of cometary phenomena, including the source of cometary activity and the nongravitational orbital motion of the nuclei. However, the internal structure of the nuclei is still largely unknown. We review herein the evidence for cometary nuclei as fluffy aggregates or primordial rubble piles, as first proposed by Donn et al. (1985) and Weissman (1986). These models assume that cometary nuclei are weakly bonded aggregations of smaller, icy‐conglomerate planetesimals, possibly held together only by self‐gravity. Evidence for this model comes from studies of the accretion and subsequent evolution of material in the solar nebula, from observations of disrupted comets, and in particular comet Shoemaker‐Levy 9, from measurements of the ensemble rotational properties of observed cometary nuclei, and from recent spacecraft missions to comets. Although the evidence for rubble pile nuclei is growing, the eventual answer to this question will likely not come until we can place a spacecraft in orbit around a cometary nucleus and study it in detail over many months to years. ESA's Rosetta mission, now en route to comet 67P/Churyumov‐Gerasimenko, will provide that opportunity.     

Shape and obliquity effects on the thermal evolution of the Rosetta target 67P/Churyumov-Gerasimenko cometary nucleus

This work is dedicated to the application to 67P/Churyumov-Gerasimenko of a new quasi-3D approach for non-spherically shaped comet nuclei with the aim to interpret the current activity of the comet in terms of initial characteristics and to predict shape and internal stratification evolution of the nucleus. The model is applied to differently shaped nuclei taking into account the characteristics of Comet 67P/Churyumov-Gerasimenko deduced from observations. We focus our attention on the combined effects that shapes and obliquity have on the comet surface and sub-surface evolution. We discuss the results in terms of activity, local dust mantle formation and disruption, erosion of the surface and internal stratigraphy.The results show that differently shaped nuclei can have different internal structures leading to different activity patterns and behaviors. Our calculations have shown that local variations in the dust and gas fluxes can be induced by the nucleus shape. The distribution of “active” areas on Comet 67P/Churyumov-Gerasimenko is different because of different shapes, reflecting the illumination conditions on the surface. These shapes can influence the structure of the inner coma, but the coma far away from the nucleus is only marginally affected by the nucleus shape. However, different comet behaviors can arise from differently shaped comet nuclei, especially in terms of local activity, surface and sub-surface characteristics and properties. The water flux local distribution is the most influenced by the shape as it is directly linked to the illumination. Irregular shapes have large shadowing effects that can result in activity patterns on the comet surface.The effects of different pole directions are discussed to see the relations with the nucleus activity and internal structure. It is shown that the orientation of the rotation axis plays a strong role on the surface evolution of 67P/Churyumov-Gerasimenko, determining seasonal effects on the fluxes. The activity of the comet changes greatly with the nucleus obliquity leading to pre-post-perihelion differences in the activity and seasonal effects. The effects of the dust deposition and crust formation on the cometary activity have also been simulated and are discussed with respect to 67P/Churyumov-Gerasimenko observations. The dust mantling is also strongly obliquity dependent, with different surface distributions of the dust-covered regions according to the different comet pole orientations. Finally, we show that our model can reproduce the fluxes behavior near perihelion in terms of amplitude and asymmetry, and we estimate 20% of the illuminated surface to be active.     

Probing the internal structure of the nuclei of comets

There is no direct evidence about the internal structure of cometary nuclei, which are mostly hidden by their gas and dust comae, and have not yet been orbited by any spacecraft. Their densities are low, typically of about 400 kg m⁻³ for 9P/Tempel 1 (that was impacted by the Deep Impact probe) and 67P/Churyumov-Gerasimenko (that is the target of the Rosetta mission). Such low densities are in favour of a high macro-porosity, or a high micro-porosity, or both. Observations of disruption or splitting of nuclei indeed suggest that some huge sub-nuclei or some meter-sized fragments could be the building blocks of comets. Analysis, from in-situ measurements and from remote light scattering observations, of the structure of the dust particles, which significantly consist of fluffy aggregates of submicron-sized grains, could be in favour of a fractal structure. However, the presence of huge icy grains in the innermost coma, and of flat layers on the surface of 9P/Tempel 1, are clues to the complexity of these objects, which have suffered drastic erosion phenomena on their elongated orbits. It is expected that the Rosetta mission will provide a fair understanding of the structure of the deep interior of the nucleus of 67P/Churyumov-Gerasimenko, thanks to the on-board CONSERT experiment.     

The nuclei of Jupiter family comets: A critical review of our present knowledge

Of the currently over 300 identified Jupiter family comets (JFCs), we have estimated nucleus sizes and shapes for fewer than 70 and have detailed nucleus observations arising from spacecraft fly-bys for just 3: 19P/Borrelly (Deep Space 1), 81P/Wild 2 (Stardust), and 9P/Tempel 1 (Deep Impact). These observations reveal similarities but also significant diversity. In this review, we make a critical assessment of our knowledge of JFC nuclei and suggest a priority list for observations of the nucleus of the JFC, 67P/Churyumov-Gerasimenko, the Rosetta target comet.     

Radar observations of Comet 2P/Encke during the 2003 apparition

The nucleus of Comet 2P/Encke was detected with the Arecibo radar during the close approach of November, 2003, making this the first comet to yield radar detections at two different apparitions. Although the measured radar cross section of 1.0 km² was close to that obtained in 1980, the Doppler bandwidth was nearly four times larger. Most of this bandwidth difference can simply be attributed to a different observing aspect relative to the spin axis proposed by Sekanina [1988, Astron. J. 95, 911] and Festou and Barale [2000, Astron. J. 119, 3119]. Comparison of the 2003 Doppler bandwidth with infrared-based size estimates supports an 11-h dominant rotation period and excludes slower 15- and 22-h periods that have also been suggested. If one assumes a short-axis-mode rotation with an 11-h period, then the Doppler bandwidth indicates that the nucleus is an oblong object with a long-axis dimension of 9 km. The estimated radar albedo of 0.05 is similar to that measured for C/IRAS-Araki-Alcock, providing further evidence that comet nuclei have relatively low surface densities of ∼0.5-1.0 g cm⁻³. No broadband echo component was detected from large coma grains despite predictions, based on optical/infrared models, that such a component might be detectable.     

Radar observations of Comet P/2005 JQ5 (Catalina)

Arecibo radar observations of Comet P/2005 JQ5 (Catalina) have produced the first delay-Doppler images of a comet nucleus and the first radar detection of large-grain ejection from a Jupiter-family comet. The nucleus is small (1.4 km diameter), rough, and rapidly rotating. The large (>cm) grains have low velocities (∼1 m/s) and a low production rate.     

Rosetta—one comet rendezvous and two asteroid fly-bys

One of the two planetary cornerstone missions of the European Space Agency is the Rosetta mission to comet 67P/Churyumov-Gerasimenko. Rosetta is a rendezvous mission with a comet nucleus, which combines an Orbiter with a Lander. It will monitor the evolution of the comet nucleus and the coma as a function of increasing and decreasing solar flux input along the comet’s pre- and post-perihelion orbit. Different instrumentations will be used in parallel, from multi-wavelength spectrometry to in-situ measurements of coma and nucleus composition and physical properties. Rosetta will go in orbit around the nucleus of its target comet 67P/Churyumov-Gerasimenko, when it is still far from the Sun and accompany the comet along its way to perihelion and beyond. In addition the Rosetta Lander Philae will land on the nucleus surface, before the comet is too active to permit such a landing (i.e. at around r = 3 AU) and examine the surface and subsurface composition of the comet nucleus as well as its physical properties.     

Comet 67P/Churyumov-Gerasimenko: Modeling of orientation and structure

The comet 67P/Churyumov-Gerasimenko is the current target of the mission Rosetta, initially planned to investigate comet 46P/Wirtanen. These two comets have similar orbits, except the distance to the Sun at perihelion, but different orbital histories and different masses. Thus, structures of the nuclei can be significantly different. The evolution of comet Wirtanen was simulated by several authors, while comet Churyumov-Gerasimenko became an object of high interest only recently and is not well investigated. In the present work we simulate the evolution of the nucleus, down to tens of meters below the surface, using an extended version of the model previously applied for comet Wirtanen [Kossacki et al., 1999. Comet 46P/Wirtanen: evolution of the subsurface layer. Icarus 142, 202-218.]. The model includes strengthening of the nucleus due to sintering of the ice grains. Simulations are performed for different latitudes, accounting for the evolution of the orbit and for changes of the nucleus orientation, as well as diurnal and seasonal changes of insolation. The calculated loss of water vapor from the comet is integrated over the nucleus surface and is compared with the observational data. We have found, that the sublimation through the dust mantle can be large enough to reproduce the profile of the total water production as a function of time from perihelion. The required dependence of thickness of the dust layer on latitude qualitatively matches present distribution of the absorbed solar flux. The non-gravitational acceleration in the comet motion together with the simulated sublimation flux are used in order to estimate the mass and the bulk density of the nucleus.     

A fully 3-dimensional thermal model of a comet nucleus

A 3-D numerical model of comet nuclei is presented. An implicit numerical scheme was developed for the thermal evolution of a spherical nucleus composed of a mixture of ice and dust. The model was tested against analytical solutions, simplified numerical solutions, and 1-D thermal evolution codes. The 3-D code was applied to comet 67P/Churyumov-Gerasimenko; surface temperature maps and the internal thermal structure was obtained as function of depth, longitude and hour angle. The effect of the spin axis tilt on the surface temperature distribution was studied in detail. It was found that for small tilt angles, relatively low temperatures may prevail on near-pole areas, despite lateral heat conduction. A high-resolution run for a comet model of 67P/Churyumov-Gerasimenko with low tilt angle, allowing for crystallization of amorphous ice, showed that the amorphous/crystalline ice boundary varies significantly with depth as a function of cometary latitude.     

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.     

Evolution of Comet Nucleus Rotation

The secular evolution of comet nucleus rotation states subject to outgassing torques is studied. The dynamical model assumes that the nucleus inertia ellipsoid is axially symmetric. The outgassing torques acting on the surface are modeled using standard cometary activity formulae. The general rotational equations of motion are derived and separately averaged over the fast rotational dynamics terms and the comet orbit. Special cases where the averaging assumptions cannot be applied are evaluated separately. The modification of the comet orbit due to comet outgassing is neglected. Resulting from this analysis is a system of secular differential equations that describes the dynamics of the comet nucleus angular momentum and rotation state. We find that the qualitative secular evolution of the rotation state is controlled by a single parameter that combines parameters related to the comet orbit and parameters related to the nucleus surface geometry and activity. From this solution, we find qualitatively different evolutionary paths for comet nuclei whose entire surface is active, as compared to nuclei with only a single active region. For surface activity models between these extremes, we show that certain evolutionary paths are more likely than others. Additionally, our solution indicates that a comet nucleus' rotational angular momentum will tend to increase over time, potentially contributing to the observed phenomenon of comet nucleus splitting.     

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.     

短周期彗星表面水冰升华分布研究

彗星是太阳系遗留的原始星子,研究彗星彗核的演化对理解太阳系其他天体的形成和演化历史具有重要意义.在太阳的辐射作用下,彗星携带的挥发性成分会发生升华,并带动尘埃运动,造成彗核物质的损失.因此,彗核的升华活动对其表面形貌甚至整体形状演化都会产生影响.从IAU (International Astronomical Union) MPC (Minor Planet Center)获取轨道数据,并考虑了彗核的自转以及进动,利用MONET (Mass lossdriven shape evolution model)形状演化模型对短周期彗星做数值模拟,计算得到了短周期彗星1P/Halley、9P/Tempel 1、 19P/Borrelly、 67P/C-G (Churyumov-Gerasimenko)、 81P/Wild 2和103P/Hartley 2在一个轨道周期内的太阳辐射能量以及表面侵蚀深度的分布,结合其动力学参数讨论了自转、进动和公转等特性对其表面水冰升华分布的影响以及造成南北侵蚀差异的可能性.     

Outgassing-induced effects in the rotational state of comet 67P/Churyumov—Gerasimenko during the Rosetta mission

The new target of the Rosetta mission is comet 67P/Churyumov-Gerasimenko (hereafter 67P/C-G). In order to support the planning of the mission, in particular the strategy during the mapping and landing phases, we have performed numerical simulations of the rotational evolution of a comet in the orbit of 67P/C-G. In these simulations, the currently known observational constraints have been taken into account and a large set of initial conditions were considered. For most of the simulations, we observe that the sublimation-induced torques produce significant changes in the rotational parameters of a 67P/C-G-like comet. Typical rates of change for the spin period from the rendezvous up to the end of the nominal mission range from 0.001 to depending on different circumstances as described in the text. At perihelion, rates of change of the orientation of the angular momentum vector amount to about 0.01-. These simulations suggest that a specific strategy should be defined in order to monitor likely variations of the rotational parameters. As an example we show a possible optimized schedule for observations with the OSIRIS instrument to determine the rotational parameters of comet 67P/C-G and their possible evolution.     

Challenging a Paradigm: Do We Need Active and Inactive Areas to Account for Near-Nuclear Jet Activity?

We briefly describe an advanced 3D gas dynamical model developed for the simulation of theenvironment of active cometary nuclei. The model canhandle realistic nucleus shapes and alternative physical models for the gas and dust production mechanism.The inner gas coma structure is computed by solving self-consistently(a) near to the surface the Boltzman Equation(b) outside of it, Euler or Navier-Stokes equations.The dust distribution is computed from multifluid ``zero-temperature' Euler equations,extrapolated with the help of a Keplerian fountain model.The evolution of the coma during the nucleus orbital and spin motion,is computed as a succession of quasi-steady solutions. Earlier versions of the model using simple,``paedagogic' nuclei have demonstrated that the surface orographyand the surface inhomogeneity contribute similarly to structuring the near-nucleusgas and dust coma,casting a shadow on the automatic attribution of such structures to ``active areas'.The model was recently applied to comet P/Halley, for whichthe nucleus shape is available. In the companion paper of this volume,we show that most near-nucleus dust structuresobserved during the 1986 Halley flybys are reproduced, assuming that the nucleus is strictly homogeneous. Here, we investigate the effect of shape perturbations and homogeneityperturbations. We show that the near nucleus gas coma structure is robust vis-a-vissuch effects. In particular, a random distribution of active and inactive areaswould not affect considerably this structure, suggesting that such areas,even if present, could not be easily identified on images of the coma.     

Evolution of depressions on Comet 67P/Churyumov-Gerasimenko: Role of ice metamorphism

This work is intended to investigate the influence of temperature-dependent metamorphism of ice on the shape of small depressions in the surface of cometary nuclei. We are mainly interested in the role of initial cohesivity of a nucleus. For this purpose we simulate sublimation of ice from the facets of initially cylindrical depressions in ice of different initial structure. The simulations account for the diurnal and orbital changes of insolation and its dependence on the current shape of the depressions. Our model includes heat transport in the cometary material and metamorphism of ice. We present the results obtained for the nucleus of the Comet 67P/Churyumov-Gerasimenko, target of the ESA cornerstone mission Rosetta.     

Modeling of meteoroid streams: The velocity of ejection of meteoroids from comets (a review)

An analytical review of the models of ejection of meteoroids from cometary nuclei is presented. Different formulas for the ejection velocity of meteoroids and the corresponding parameters are discussed and compared with the use of comet Halley and the Geminids meteoroid stream as examples. The ejection velocities obtained from observations of the dust trails of comets are discussed, and the values for comets 2P/Encke, 4P/Faye, 17P/Holmes, 22P/Kopff, and 67P/Churyumov-Gerasimenko are compared to the velocities yielded by Whipple’s model. The uncertainty intervals of the results are estimated.     

Comparison of some characteristics of comets 1P/Halley and 67P/Churyumov–Gerasimenko from the <Emphasis Type="Italic">Vega</Emphasis> and <Emphasis Type="Italic">Rosetta</Emphasis> mission data

On March 6 and 9, 1986, for the first time in the history of science, the Russian spacecraft Vega-1 and -2 approached the nucleus of comet 1P/Halley and flew by at a small distance. A while later, on March 14, 1986, the Giotto spacecraft (European Space Agency (ESA)) followed them. Together with the Japanese spacecraft Suisei (Japan Aerospace Exploration Agency (JAXA)), they obtained spaceborne investigations of cometary nuclei. Direct studies of cometary bodies that bear traces of the Solar System formation were continued in the next missions to comets. Starting from 2014 and up to 2016 September, the Rosetta spacecraft (ESA), being in a low orbit around the nucleus of comet 67P/Churyumov–Gerasimenko, has performed extremely sophisticated investigations of this comet. Here, we compare some results of these missions. The paper is based on the reports presented at the memorial conference dedicated to the 30th anniversary of the Vega mission, which took place at the Space Research Institute of the Russian Academy of Sciences in March, 2016, and does not pretend to comprehensively cover the problems of cometary physics.     

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.