The size-frequency distribution of dormant Jupiter family comets

We estimate the total number and the slope of the size-frequency distribution (SFD) of dormant Jupiter family comets (JFCs) by fitting a one-parameter model to the known population. We first select 61 near-Earth objects (NEOs) that are likely to be dormant JFCs because their orbits are dynamically coupled to Jupiter [Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J., Levison, H.F., Michel, P., Metcalfe, T.S., 2002a. Icarus 156, 399-433]. Then, from the numerical simulations of Levison and Duncan [1997. Icarus 127, 13-32], we construct an orbit distribution model for JFCs in the NEO orbital element space. We assume an orbit-independent SFD for all JFCs, the slope of which is our unique free parameter. Finally, we compute observational biases for dormant JFCs using a calibrated NEO survey simulator [Jedicke, R., Morbidelli, A., Spahr, T., Petit, J., Bottke, W.F., 2003. Icarus 161, 17-33]. By fitting the biased model to the data, we estimate that there are ∼75 dormant JFCs with H<18 in the NEO region and that the slope of their cumulative SFD is −1.5±0.3. Our slope for the SFD of dormant JFCs is very close to that of active JFCs as determined by Weissman and Lowry [2003. Lunar Planet. Sci. 34. Abstract 2003]. Thus, we argue that when JFCs fade they are likely to become dormant rather than to disrupt and that the fate of faded comets is size-independent. Our results imply that the size distribution of the JFC progenitors—the scattered disk trans-neptunian population—either (i) has a similar and shallow SFD or (i^′) is slightly steeper and physical processes acting on the comets in a size-dependent manner creates the shallower active comet SFD. Our measured slope, typical of collisionally evolved populations with a size-dependent impact strength [Benz, W., Asphaug, E., 1999. Icarus 142, 5-20], suggests that scattered disk bodies reached collisional equilibrium inside the protoplanetary disk prior to their removal from the planetary region.     

The population of faint Jupiter family comets near the Earth

We study the population of faint Jupiter family comets (JFCs) that approach the Earth (perihelion distances q<1.3 AU) by applying a debiasing technique to the observed sample. We found for the debiased cumulative luminosity function (CLF) of absolute total magnitudes H₁₀ a bimodal distribution in which brighter comets (H₁₀?9) follow a linear relation with a steep slope α=0.65±0.14, while fainter comets follow a much shallower slope α=0.25±0.06 down to H₁₀∼18. The slope can be pushed up to α=0.35±0.09 if a second break in the H₁₀ distribution to a much shallower slope is introduced at H₁₀∼16. We estimate a population of about 10³ faint JFCs with q<1.3 AU and 10<H₁₀<15 (radii ∼0.1-0.5 km). The shallowness of the CLF for faint near-Earth JFCs may be explained either as: (i) the source population (the scattered disk) has an equally very shallow distribution in the considered size range, or (ii) the distribution is flattened by the disintegration of small objects before that they have a chance of being observed. The fact that the slope of the magnitude distribution of the faint active JFCs is very similar to that found for a sample of dormant JFCs candidates suggests that for a surviving (i.e., not disintegrated) object, the probability of becoming dormant versus keeping some activity is roughly size independent.     

Plasma environment of Jupiter family comets

As any comet nears the Sun, gas sublimes from the nucleus taking dust with it. Jupiter family comets are no exception. The neutral gas becomes ionized, and the interaction of a comet with the solar wind starts with ion pickup. This key process is also important in other solar system contexts wherever neutral particles become ionized and injected into a flowing plasma such as at Mars, Venus, Io, Titan and interstellar neutrals in the solar wind. At comets, ion pickup removes momentum and energy from the solar wind and puts it into cometary particles, which are then thermalised via plasma waves. Here we review what comets have shown us about how this process operates, and briefly look at how this can be applied in other contexts. We review the processes of pitch angle and energy scattering of the pickup ions, and the boundaries and regions in the comet-solar wind interaction. We use in-situ measurements from the four comets visited to date by spacecraft carrying plasma instrumentation: 21P/Giacobini-Zinner, 1P/Halley, 26P/Grigg-Skjellerup and 19P/Borrelly, to illustrate the process in action. While, of these, comet Halley is not a Jupiter class comet, it has told us the most about cometary plasma environments. The other comets, which are from the Jupiter family, give an interesting comparison as they have lower gas production rates and less-developed interactions. We examine the prospects for Rosetta at comet Churyumov-Gerasimenko, another Jupiter family comet where a wide range of gas production rates will be studied.     

Hilda asteroids among Jupiter family comets

Hilda asteroids and comets are similar from the compositional point of view. The D-taxonomic class prevailing among Hildas has all the characteristics found in cometary spectra. Jupiter Family Comets (JFCs) coming from the trans-neptunian region are under the gravitational control of Jupiter, making them a dynamically unstable population with a mean dynamical lifetime of 10⁴ to ¹⁰⁵ years. In contrast, Hilda asteroids residing in the 3:2 mean motion resonance with Jupiter are a very stable population. But once they escape from the resonance, they are dynamically controlled by Jupiter, and in this sense their behavior resembles that of JFC. We performed a numerical simulation to analyze the dynamical evolution that Hildas follow after escaping from the resonance, and their contribution to the JFC population. We found that 8% of the particles leaving the resonance end up impacting Jupiter. 98.7% of the escaped Hildas live at least 1000 years as a JFC, with a mean lifetime of 1.4×¹⁰⁶ years. In particular, escaped Hildas stay mainly in the region of perihelion distances greater than 2.5 AU. On the other hand, the number of escaped Hildas reaching the inner Solar System (q<2.5 AU) is negligible. So, there are almost no Hilda asteroids among the NEO population. We also analyzed the possibility that the Shoemaker-Levy 9 were an escaped Hilda asteroid. In this case, it would be possible to give stronger constraints to its pre-capture orbital elements.     

Meteor streams associated with Jupiter family comets

Meteor showers have been observed for a considerable time, and the cause, meteoroids from a meteoroid stream ablating in the Earth's atmosphere, has also been understood for centuries. The connection between meteoroid streams and comets was also established 150 years ago. Since that time our ability both to understand the physics and to numerically model the situation has steadily increased. We will review the current state of knowledge. However, just as there are differences between the behaviour of long period comets, Halley family comets and Jupiter family comets, so also differences exist between the associated meteoroid streams. Streams associated with Jupiter family comets show much more variety in their behaviour, driven by the gravitational perturbations from Jupiter. The more interesting showers associated with Jupiter family comets will be discussed individually.     

Optical properties of dust from Jupiter family comets

To try to define specific physical properties of the dust of Jupiter-family comets (JFCs), we compare the light scattered by them. Amongst the more than 1000 JFCs, less than 200 are numbered, 40 of them being rather bright. In the present work we use data from the latter. In situ observations of three nuclei show low albedo surfaces. The albedo of the dust particles in the coma is low, with generally a red colour. The A(α)fρ product is a measure of cometary activity and secular changes. Images of different regions (jets and fans) give indications on the nucleus rotation and position of the emitting areas, as compared to the position of the rotation axis. Differences in physical properties between the particles in different regions are pointed out by differences in the linear polarization of the scattered light and by spectral variations in brightness and polarization. Jupiter family comets are considered as dust-poor comets. Tails and trails’ studies give an estimation of the size distribution of the particles. However the dust production rates depend on the largest particles (up to centimetre size), which are mainly observed in the trails where large dark compact particles are found. These dark particles are also responsible for the high polarization in the inner most coma of some comets. The meaning, in terms of physical properties, of the linear polarization is discussed through different examples such as 2P/Encke, 9P/Tempel 1 or the fragments of 73P/Schwassmann-Wachmann 3. Cometary outbursts and splitting events show that the properties of the dust ejected from the interior of the nucleus are similar to the ones of more active comets (new or with larger semi-major axis).     

Nuclear magnitudes and the size distribution of Jupiter family comets

We present a new catalog of absolute nuclear magnitudes of Jupiter family (JF) comets, which is an updated version of our previous catalog [Tancredi, G., Fernández, J.A., Rickman, H., Licandro, J., 2000. Astron. Astrophys. Suppl. Ser. 146, 73-90]. From the new catalog we find a linear cumulative luminosity function (CLF) of slope 0.54±0.05 for JF comets with q?2.5 AU. By considering this CLF combined with the few measured geometric albedos with their respective uncertainties, and assuming a canonical albedo of 0.035±0.012 for those comets with undetermined albedos, we derive a cumulative size distribution that follows a power-law of index −2.7±0.3. The slope is similar to that derived from some theoretical collisional models and from some populations of Solar System bodies like the trans-neptunian objects. We also discuss and compare our size distribution with those by other authors that have recently appeared in the literature. Some striking differences in the computed slopes are explained in terms of biases in the studied samples, the different weights given to the brightest members of the samples, and discrepancies in the values of a few absolute nuclear magnitudes. We also compute sizes and fractions of active surface area of JF comets from their estimated absolute nuclear magnitudes and their water production rates. With the outgassing model that we use, about 60% of the computed fractions f of active surface area are found to be smaller than 0.2, with one case (28P/Neujmin 1) of no more than 0.001, which suggests that JF comets may transit through stages of very low activity, or even dormancy. There is an indication that JF comets with radii RN?3 km have active fractions f?0.01, which might be due to the rapid formation of insulating dust mantles on larger nuclei.     

On the population, physical decay and orbital distribution of Jupiter family comets: Numerical simulations

We study the Jupiter family comet (JFC) population assumed to come from the Scattered Disk and transferred to the Jupiter’s zone through gravitational interactions with the Jovian planets. We shall define as JFCs those with orbital periods and Tisserand parameters in the range 2<T?3.1, while those comets coming from the same source, but that do not fulfill the previous criteria (mainly because they have periods ) will be called ‘non-JFCs’. We performed a series of numerical simulations of fictitious comets with a purely dynamical model and also with a more complete dynamical-physical model that includes besides nongravitational forces, sublimation and splitting mechanisms. With the dynamical model, we obtain a poor match between the computed distributions of orbital elements and the observed ones. However with the inclusion of physical effects in the complete model we are able to obtain good fits to observations. The best fits are attained with four splitting models with a relative weak dependence on q, and a mass loss in every splitting event that is less when the frequency is high and vice versa. The mean lifetime of JFCs with radii and is found to be of about 150-200 revolutions (∼. The total population of JFCs with radii within Jupiter’s zone is found to be of 450±50. Yet, the population of non-JFCs with radii in Jupiter-crossing orbits may be ∼4 times greater, thus leading to a whole population of JFCs + non-JFCs of ∼2250±250. Most of these comets have perihelia close to Jupiter’s orbit. On the other hand, very few non-JFCs reach the Earth’s vicinity (perihelion distances ) which gives additional support to the idea that JFCs and Halley-type comets have different dynamical origins. Our model allows us to define the zones of the orbital element space in which we would expect to find a large number of JFCs. This is the first time, to our knowledge, that a physico-dynamical model is presented that includes sublimation and different splitting laws. Our work helps to understand the role played by these erosion effects in the distribution of the orbital elements and lifetimes of JFCs.     

Activity of comets: Gas transport in the near-surface porous layers of a cometary nucleus

The gas transport through non-volatile random porous media is investigated numerically. We extend our previous research of the transport of molecules inside the uppermost layer of a cometary surface ( [Skorov and Rickman, 1995] and [Skorov et al., 2001]). We assess the validity of the simplified capillary model and its assumptions to simulate the gas flux trough the porous dust mantle as it has been applied in cometary physics. A microphysical computational model for molecular transport in random porous media formed by packed spheres is presented. The main transport characteristics such as the mean free path distribution and the permeability are calculated for a wide range of model parameters and compared with those obtained by more idealized models. The focus in this comparison is on limitations inherent in the capillary model. Finally a practical way is suggested to adjust the algebraic Clausing formula taking into consideration the nonlinear dependence of permeability on layer porosity. The retrieved dependence allows us to accurately calculate the permeability of layers whose thickness and porosity vary in the range of values expected for the near-surface regions of a cometary nucleus.     

Stochasticity of two comets in resonance with Jupiter

The orbital stochasticity of comets P/Ciffréo (1985 XVI) and P/Maury (1985 VI), at the present time near the 5/3 and 4/3 resonances with Jupiter, is investigated using Lyapunov Characteristic Indicators. First results indicate a strong stochastic behaviour for the two comets, mainly induced by encounters with Jupiter, which looks roughly like the behaviour of the group of comets in 1/1 resonance with Jupiter.     

Coma morphology of Jupiter-family comets

We present a general review of cometary coma morphology, with specific reference to how it used in studies of Jupiter-family comets. We introduce the most common features that are seen in gas and dust observations, and summarize of how they are used to infer the properties of the nucleus and coma. We also expand the discussion to cover other topics relating to morphology, including the general shape of the coma (characterized by radial gradient profiles) and spatial maps of the color, albedo and polarization of the dust. We address the pros and cons of the different approaches used in the interpretation and analysis of the features. Finally, we review the results obtained for specific comets and compare the Jupiter-family comets to those from other classes.     

The internal structure of Jupiter family cometary nuclei from Deep Impact observations: The “talps” or “layered pile” model

We consider the hypothesis that the layering observed on the surface of Comet 9P/Tempel 1 from the Deep Impact spacecraft and identified on other comet nuclei imaged by spacecraft (i.e., 19P/Borrelly and 81P/Wild 2) is ubiquitous on Jupiter family cometary nuclei and is an essential element of their internal structure. The observational characteristics of the layers on 9P/Tempel 1 are detailed and considered in the context of current theories of the accumulation and dynamical evolution of cometary nuclei. The works of Donn [Donn, B.D., 1990. Astron. Astrophys. 235, 441-446], Sirono and Greenberg [Sirono, S.-I., Greenberg, J.M., 2000. Icarus 145, 230-238] and the experiments of Wurm et al. [Wurm, G., Paraskov, G., Krauss, O., 2005. Icarus 178, 253-263] on the collision physics of porous aggregate bodies are used as basis for a conceptual model of the formation of layers. Our hypothesis is found to have implications for the place of origin of the JFCs and their subsequent dynamical history. Models of fragmentation and rubble pile building in the Kuiper belt in a period of collisional activity (e.g., [Kenyon, S.J., Luu, J.X., 1998. Astron. J. 115, 2136-2160; 1999a. Astron. J. 118, 1101-1119; 1999b. Astrophys. J. 526, 465-470; Farinella, P., Davis, D.R., Stern, S.A., 2000. In: Mannings, V., Boss, A.P., Russell, S.S. (Eds.), Protostars and Planets IV. Univ. of Arizona Press, Tucson, pp. 1255-1282; Durda, D.D., Stern, S.J., 2000. Icarus 145, 220-229]) following the formation of Neptune appear to be in conflict with the observed properties of the layers and irreconcilable with the hypothesis. Long-term residence in the scattered disk [Duncan, M.J., Levison, H.F., 1997. Science 276, 1670-1672; Duncan, M., Levison, H., Dones, L., 2004. In: Festou, M., Keller, H.U., Weaver, H.A. (Eds.), Comets II. Univ. of Arizona Press, Tucson, pp. 193-204] and/or a change in fragmentation outcome modeling may explain the long-term persistence of primordial layers. In any event, the existence of layers places constraints on the environment seen by the population of objects from which the Jupiter family comets originated. If correct, our hypothesis implies that the nuclei of Jupiter family comets are primordial remnants of the early agglomeration phase and that the physical structure of their interiors, except for the possible effects of compositional phase changes, is largely as it was when they were formed. We propose a new model for the interiors of Jupiter family cometary nuclei, called the talps or “layered pile” model, in which the interior consists of a core overlain by a pile of randomly stacked layers. We discuss how several cometary characteristics—layers, surface texture, indications of flow, compositional inhomogeneity, low bulk density low strength, propensity to split, etc., might be explained in terms of this model. Finally, we make some observational predictions and suggest goals for future space observations of these objects.     

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 origin of the CN radical in comets: A review from observations and models

The origin of CN radicals in comets is not completely understood so far. We present a study of CN and HCN production rates and CN Haser scale lengths showing that: (1) at heliocentric distances larger than 3 AU, CN radicals could be entirely produced by HCN photolysis; (2) closer to the Sun, for a fraction of comets CN production rates are higher than HCN ones whereas (3) in the others, CN distribution cannot be explained by the HCN photolysis although CN and HCN production rates seem to be similar. Thus, when the comets are closer than 3 AU to the Sun, an additional process to the HCN photolysis seems to be required to explain the CN density in some comets.The photolysis of HC₃N or C₂N₂ could explain the CN origin. But the HC₃N production rate is probably too low to reproduce CN density profile, even if uncertainties on its photolysis leave the place for all possible conclusions. The presence of C₂N₂ in comets is a reliable hypothesis to explain the CN origin; thus, its detection is a challenging issue. Since C₂N₂ is very difficult to detect from ground-based observations, only in situ measurements or space observations could determine the contribution of this compound in the CN origin.Another hypothesis is a direct production of CN radicals by the photo- or thermal degradation of complex refractory organic compounds present on cometary grains. This process could explain the spatial profile of CN inside jets and the discrepancy noted in the isotopic ratio ¹⁴N/¹⁵N between CN and HCN. Laboratory studies of the thermal and UV-induced degradation of solid nitrogenated compounds are required to model and validate this hypothesis.     

Short-period Jupiter family comets after Stardust

The NASA Stardust mission has provided for laboratory study an extensive data set of cometary dust of known provenance (from comet 81P/Wild 2) yielding detailed insights into the composition of the comet. Combined with the results of data from other missions to short-period Jupiter family comets (JFC), this has greatly deepened the understanding of such objects. If depressions on the surface of comet 81P/Wild 2 are all taken as evidence of impact cratering, their number suggests a long occupancy in the outer region of the Solar System. The dust from comet 81P/Wild 2 has been shown to be heavily deficient in pre-Solar grains and rich in materials formed at high temperatures in the inner Solar System. Although it is too early to know if this is typical of JFC, it does argue for rapid and thorough mixing of materials in the disk on timescales related to comet formation, and may also suggest outward migration of small icy bodies after their formation. Thus, instead of providing mainly new knowledge of the pre-Solar materials expected to be rich in comets, Stardust and comet 81P/Wild 2 have instead focussed attention on large-scale transport processes during the critical period when cometary parent bodies were forming in the early Solar System.     

Structure and evolution of the Jupiter family

Both physical and dynamical issues are important in order to judge the origin and evolution of the Jupiter family of short-period comets. The steady-state condition for maintaining this structure at its present size by captures from the classical Oort cloud is reviewed on the basis of recent results concerning the absolute number of Jupiter family comets as a function of perihelion distance as well as the coupled physical and dynamical evolutions that evidently occur. Like in most earlier investigations, a clear shortage is found in the classical Oort cloud source. The shortage seems, however, less extreme than sometimes assumed. Monte Carlo simulations are envisaged as a way to shed further light on the fate of Jupiter family comets.     

Appearance of layered structures in numerical simulations of polydisperse bodies accretion: Application to cometary nuclei

A model for the aggregation of size distribution of cometesimals (Gaussian or power law) into cometary nuclei is developed. Upon disruption induced by collisions, sticking and evolution of the tensile strength and density of the cometesimals by sintering processes are taken into account. The resulting cometary nuclei present specific internal structures that have been quantified to allow the comparison with observational constraints and future in situ observations and cometary nucleus sounding with the CONSERT radar on-board the Rosetta mission. A parameter called the homogeneity exponent, μ, determines different aggregation regimes. Fractal aggregates are formed for μ < 0.4. Radial variations in tensile strength appear for 0.4 < μ < 0.6 and vanish for larger values of μ. The initial size distribution (following a Gaussian or power law) of aggregating cometesimals does not influence strongly these values but can change the extent of corresponding layers. If the layering observed on the surface of some cometary nuclei occurs often and originates from primordial structures, this constrains the velocity distribution of aggregating bodies to follow v∝m^-0.25, while a differential size distribution following a power law with exponent between −2 and −3 should result for large bodies, in agreement with current estimations of the size distributions. Such a layered structure would lead to more cohesive, dense and less porous material located near the center of mass of the nucleus predicting an increase of bulk density of comet nuclei with their erosion state.     

Radar observations of 8P/Tuttle: A contact-binary comet

Arecibo radar imagery of Comet 8P/Tuttle reveals a 10-km-long nucleus with a highly bifurcated shape consistent with a contact binary. A separate echo component was also detected from large (>cm-size), slow-moving grains of the type expected to contribute to the Ursid meteor stream.     

The capture of comets by Jupiter

The capture of comets with parabolic orbits by Jupiter is investigated. The influence of the gravitational force of the Sun on the cometary orbit during the passage of Jupiter's sphere of influence is taken into account. A comparison of the present results with previous calculations demonstrate the importance of the solar perturbations.It is also shown that captures of comets with parabolic orbits and repeated close passages to Jupiter cannot explain all of the observed cometary orbits found in the family of Jupiter.     

A ∼25 K temperature of formation for the submicron ice grains which formed comets

Following the observations of ice grains in cometary comae and their size distributions, we reexamined experimentally our previous conclusion that the ice grains which agglomerated to form comet nuclei were formed at ∼25 K. The suggestion of a ∼25 K formation temperature was confirmed experimentally. Moreover, we suggest that these ice grains had to be of submicron size.