Issue |
A&A
Volume 521, October 2010
|
|
---|---|---|
Article Number | A23 | |
Number of page(s) | 5 | |
Section | Planets and planetary systems | |
DOI | https://doi.org/10.1051/0004-6361/201014884 | |
Published online | 15 October 2010 |
Impactor flux and cratering on the Pluto-Charon system
G. C. de Elía - R. P. Di Sisto - A. Brunini
Facultad de Ciencias Astronómicas y Geofísicas, Universidad Nacional de La Plata and Instituto de Astrofísica de La Plata, CCT La Plata-CONICET-UNLP Paseo del Bosque S/N (1900), La Plata, Argentina
Received 28 April 2010 / Accepted 2 July 2010
Abstract
Aims. We study the impactor flux and cratering on
Pluto and Charon caused by the collisional evolution of Plutinos.
Plutinos are trans-Neptunian objects located at 39.5 AU, in the 3:2 mean motion
resonance with Neptune.
Methods. We develop a statistical code that includes
catastrophic collisions and cratering events, and takes into account
the stability and instability zones of the 3:2 mean motion resonance
with Neptune. Our numerical algorithm proposes different initial
populations that account for the uncertainty in the size distribution
of Plutinos at small sizes.
Results. Depending on the initial population, our
results indicate the following. The number of D
> 1 km Plutinos streaking Pluto over 3.5 Gyr
is between 1271 and 5552. For Charon, the number of D
> 1 km Plutino impactors is between 354 and
1545. The number of D > 1 km
craters on Pluto produced by Plutinos in the past 3.5 Gyr is between
43 076 and 113 879. For Charon, the number of D
> 1 km craters is between 20 351 and
50 688. On the other hand, the largest Plutino impactor onto
Pluto has a diameter of between 17 and 23 km, which produces a crater
with a diameter of
31-39 km.
In the same way, the largest Plutino impactor onto Charon has a
diameter of between
10
and 15 km, which produces a crater with a diameter of
24-33 km.
Finally, we test the dependence of results on the number of Pluto-sized
objects in the Plutino population. If two Pluto-sized objects are
assumed in the 3:2 Neptune resonance, the total number of Plutino
impactors onto both Pluto and Charon with diameters D
> 1 km is a factor of
1.6-1.8 larger than that obtained considering
only one Pluto-sized object in this resonant region.
Conclusions. Given the structure of the
trans-Neptunian region, with its dynamically different populations, it
is necessary to study in detail the contribution of all the potential
sources of impactors onto the Pluto-Charon system, to determine the
main contributor and the whole production of craters. Then, we will be
able to contrast those studies with observations, which will help us to
understand the geological processes and history of the surface of those
worlds.
Key words: methods: numerical - Kuiper belt: general - planets and satellites: individual: Pluto
1 Introduction
Pluto and Charon are members of a vast population of icy bodies beyond
Neptune and constitute the first discovered binary trans-Neptunian
object (TNO).
In the trans-Neptunian region, there are four dynamical classes (Chiang
et al. 2007):
the classical objects with semimajor axes a
greater than 42 AU
and low eccentricity orbits, the scattered disk objects (SDOs) with
perihelion distances of q > 30 AU
and high eccentricities, the resonant objects in mean motion resonances
with Neptune, and the centaurs with perihelion distances of q
< 30 AU.
The 3:2 mean motion resonance with Neptune, located at 39.5 AU,
is the most densely populated one in the trans-Neptunian region. The
residents of this resonant region are usually called Plutinos because
of the analogy of their orbits with that of Pluto, which is its most
representative member. Apart from Pluto and
its largest moon Charon, the Minor Planet Center (MPC) database
contains
200 Plutino
candidates.
Pluto and Charon have been exposed to impacts with minor bodies as all the objects in our Solar System. Cratering is one of the most important processes determining the morphology of the surface of a Solar System object. The understanding and quantification of the impactor source population onto an object and the observation of the object surface help us to understand the dynamical and physical history of both the impactor population and the target.
All of the detailed knowledge of the surface composition of Pluto and Charon has been obtained from telescopic observations of the spectrum of sunlight reflected from their surfaces. Pluto's reflectance spectrum shows absorption bands of methane ice, and an absorption band that could be related to the presence of CO and nitriles. The surface of Charon can be modeled by pure water ice darkened by a spectrally neutral continuum absorber (Protopapa et al. 2008).
However, it will not be until 2015 that we will have a real idea of the morphology of the surfaces of Pluto and Charon, with the fly-by of the Nasa's New Horizon Pluto-Kuiper belt mission to Pluto system. The New Horizons mission is the first one to the Pluto system and the Kuiper belt, and was launched on 19 January 2006 on a Jupiter Gravity Assist trajectory toward the Pluto system for a 14 July 2015 closest approach. It will study the Pluto system over a 5-month period beginning in early 2015 in particular providing measurements of cratering records.
The comparison of the predicted theoretical crater production from a given source with the observed surface of Pluto and Charon may account for the geological processes acting on the surface of the objects and if cratering collisions onto Pluto and Charon are an important surface modification process.
It is then very important to study all the possible sources of crater production on Pluto and Charon to characterize the total crater production and compare it with observations.
The possible main contributors to the impactor flux on Pluto and Charon would be on the one hand all type of comets. This topic has been addressed by Weissman & Stern (1994), Durda & Stern (2000), and Zahnle et al. (2003). In general, comets have eccentric orbits that can cross Pluto's orbit when they enter the planetary region from their source, either the trans-Neptunian zone or the Oort cloud, or when they leave the inner Solar System because of the perturbations of the planets.
Weissman & Stern (1994)
estimated current impact rates of comets on Pluto and Charon. They
showed that cratering on both bodies is dominated by Kuiper belt and
inner Oort cloud comets.
Durda & Stern (2000)
then calculated collision rates in the Kuiper belt and Centaur region
by means of a numerical model. They estimated that the flux of Kuiper
belt projectiles onto Pluto and Charon is 3-5 times that of Weissman &
Stern (1994).
Zahnle et al. (2003)
studied the cratering rates for the moons of the jovian planets and
Pluto produced mainly by ecliptic comets, obtaining results consistent
with the previous estimates.
We later analyze these studies in a discussion section.
On the other hand, Plutinos may be the other important source of impactors on Pluto and Charon. Plutinos shared the same dynamical conditions as Pluto and Charon, they are all located in the same 3:2 mean motion resonance with Neptune so they have a certain collision probability. de Elía et al. (2008) simulated the collisional evolution of Plutinos and determined collisional rates among these objects. This work allowed us to study the impactor flux on Pluto and Charon caused by the collisional evolution of Plutinos. In this paper, we evaluate the contribution of Plutinos to the impactor flux and cratering on Pluto and Charon and also determine whether the Plutino population can be considered a primary source of impactors on the Pluto-Charon system.
2 The full model
To simulate the collisional and dynamical evolution of the Plutino population, we use the statistical code developed by de Elía et al. (2008). This algorithm considers catastrophic collisions and cratering events, and takes into account the main dynamical characteristic associated with the 3:2 mean motion resonance with Neptune. In the following, we give a brief description of the initial populations, the collisional parameters and the main dynamical considerations used in our model.
2.1 Initial populations
From Kenyon et al. (2008),
the cumulative size distribution of the resonant population of the
trans-Neptunian region displays a break at a diameter D
near 40-80 km. For larger resonant objects, the population
seems to have a shallow
size distribution with a cumulative power-law index of 3. From
this, the general form of the cumulative initial population used in our
model to study the
collisional and dynamical evolution of the Plutinos can be written as
where C2 adopts a value of

- initial population 1, with a cumulative power-law index p
of 3.0 for
km;
- initial population 2, with a cumulative power-law index p
of 2.7 for
km;
- initial population 3, with a cumulative power-law index p
of 2.4 for
km.
2.2 Collisional parameters
We adopt constant values of the intrinsic collision probability
and
the mean impact velocity
for Plutinos derived by Dell'Oro et al. (2001).
Based on a sample of 46 Plutinos, these authors computed
values of
and
of
km-2 yr-1
and 1.44
0.71 km s-1,
respectively.
As for the impact strength, O'Brien & Greenberg (2005) showed
that the general shape of the final
evolved asteroid population is determined primarily by the impact
energy required for dispersal QD,
but variations in
the shattering impact specific energy QS
and the inelasticity parameter
can affect this final population
even if QD
is held the same. According to these arguments, we decide to choose a
combination of the parameters
QS and
that yield the QD
law from Benz & Asphaug (1999)
for icy bodies at 3 km s-1.
de Elía & Brunini (2007) analyzed the dependence of their numerical simulations on the shattering impact specific energy QS. According to this work, the smallest gaps between QS and QD curves lead to the smallest wave amplitudes in the size distribution of the final evolved population as well as to the highest ejection rates of collisional fragments. Moreover, that study also indicates that the formation of families is more effective for the simulations with a small gap between QS and QD laws. Following these arguments, we decide to use two QS laws, QS,1, and QS,2, with a small and a large gap with respect to the QD law from Benz & Asphaug (1999) for icy bodies at 3 km s-1, respectively. The results discussed in this work are those obtained using the QS,1 law. In Sect. 3.4, we develop numerical simulations using the QS,2 law to test the dependence of our results on this collisional parameter. Figure 1 shows the two QS laws used in our simulations and the QD law from Benz & Asphaug (1999) for icy bodies at 3 km s-1.
![]() |
Figure 1: Impact strength. The dashed lines represent the two different QS laws used in our simulations. The QD law from Benz & Asphaug (1999) for icy bodies at 3 km s-1 is plotted as a solid line. |
Open with DEXTER |
Once the QS
law is specified, we fit the inelasticity parameter
to obtain the Benz &
Asphaug (1999)
QD
law.
According to O'Brien & Greenberg (2005), we
express the parameter
as
where






2.3 Dynamical considerations
To study the orbital space occupied by the Plutino population, we develop a numerical integration of 197 Plutino candidates extracted from the Minor Planet Center database with semimajor axes between 39 and 40 AU. These objects are assumed to be massless particles subject to the gravitational field of the Sun (including the masses of the terrestrial planets) and the perturbations of the four giant planets. The simulation is performed with the simplectic code EVORB from Fernández et al. (2002). The evolution of the test particles is followed for 107 years, which is a timescale greater than any secular period found in this resonance (Morbidelli 1997). From this, we build maps of the distribution of Plutinos in the orbital element planes (a, e) and (a, i), which allows us to determine the main stability regions of the 3:2 Neptune resonance. These maps are used to assign a characteristic orbit for every colliding Plutino and to specify the final fates of the different fragments generated in the collisional evolution. A detailed description of this procedure can be found in de Elía et al. (2008).
3 Results
The previously described collisional code allows us to calculate the collisional rates of Plutinos onto Pluto and Charon. From this, it is possible to calculate the impactor flux of Plutinos of different sizes on the Pluto-Charon system. Using a suitable expression, we can calculate the crater diameters produced by Plutinos on Pluto and Charon. We present here our main results concerning the impactor flux and cratering onto Pluto and Charon due to the collisional evolution of the Plutino population.
![]() |
Figure 2: Cumulative number of Plutino impactors onto Pluto a) and Charon b) over the past 3.5 Gyr as a function of impactor diameter. |
Open with DEXTER |
3.1 Impactor flux onto Pluto and Charon
Figure 2
a) shows the cumulative number of Plutino impacts onto Pluto over the
past 3.5 Gyr as a function of impactor diameter, obtained
from the three different initial populations defined in
Sect. 2.1.
Moreover, Table 1
summarizes some of our results for the impactor flux onto Pluto
produced by the collisional evolution of Plutinos. From this, the
number of D > 1 km Plutinos
striking Pluto over 3.5 Gyr is between 1271 and 5552,
while the largest Plutino expected to have impacted Pluto during the
past 3.5 Gyr had a diameter of
17-23 km,
depending on the initial size distribution.
On the other hand, Fig. 2b) shows the
cumulative number of Plutino impacts onto Charon over the past
3.5 Gyr as a function of impactor diameter.
Moreover, results about the impactor flux onto Charon due to the
collisional evolution of Plutinos are summarized in Table 2.
From this, the number of D > 1 km
Plutinos striking Charon over 3.5 Gyr is
between 354 and 1545, while the largest Plutino
expected to have impacted Charon during the past 3.5 Gyr had a
diameter of
10-15 km,
depending on the initial size distribution.
Table 1: Collisional evolution results about Pluto.
Table 2: Collisional evolution results about Charon.
3.2 Cratering on Pluto and Charon
Using the impactor fluxes previously estimated, it is possible to
calculate the number of craters on Pluto and Charon over the past
3.5 Gyr as a function of crater diameter.
From Holsapple (1993),
the diameter
of a hemispherical crater is given by
where







![]() |
(4) |
where U is the hyperbolic encounter velocity, which is equal to the mean impact velocity given in Sect. 2.2, and

Figure 3
a) shows the cumulative number of craters on Pluto produced by Plutinos
over 3.5 Gyr as a function of crater diameter.
From this, the number of D > 1 km
craters on Pluto is between 43 076 and 113 879, while
the largest crater diameter is between
31 and
39 km, depending on the initial size distribution.
In the same way, Fig. 3 b) shows the
cumulative number of craters on Charon produced by Plutinos over
3.5 Gyr as a function of crater diameter.
From this, the number of D > 1 km
craters on Charon is between 20 351 and 50 688, while
the largest crater has a diameter of
24-33 km,
depending on the initial size distribution.
![]() |
Figure 3: Cumulative number of craters on Pluto a) and Charon b) produced by Plutinos over the past 3.5 Gyr as a function of crater diameter. |
Open with DEXTER |
Table 3:
Values of the density ,
the surface gravity g, and the encounter velocity
used for Pluto and Charon.
3.3 Dependence of results on the number of Pluto-sized objects
Kenyon et al. (2008)
suggested that the resonant trans-Neptunian population contains 0.01-0.05
in objects
with
-40 km. From these
estimates, it is possible to infer the existence of five Pluto-sized
objects in the whole resonant population. If these five Pluto-sized
objects
were all in the 3:2 Neptune resonance, we would have an upper limit for
the large objects in this resonance.
Brown (2008),
based on the completeness of the current surveys, argued that two or
three
more large KBOs are likely awaiting discovery. Since the existence of
five Pluto-sized objects in the 3:2 Neptune resonance would seem to be
an overestimation, we assume an upper limit of two Pluto-sized objects
in this resonant region to analyze the dependence of our simulations on
the number of large objects in this resonance.
If two Pluto-sized objects are assumed in the 3:2 Neptune
resonance, the total number of Plutino impactors onto both Pluto and
Charon with diameters D > 1 km is a
factor of 1.6-1.8 larger
than that obtained considering only one Pluto-sized object in this
resonant region.
3.4 Dependence of results on the shattering impact specific energy QS
To test the dependence of our results on the shattering impact specific energy QS, we carry out numerical simulations using the QS,1 and QS,2 laws, which show a small and a large gap with respect to the QD law from Benz & Asphaug (1999) for icy bodies at 3 km s-1, respectively. We find that the QS,2 law produces a size distribution with a wave amplitude larger than that of the one obtained using the QS,1 law as was noticed by de Elía & Brunini (2007). However, the whole impactor flux is almost preserved using both laws.
4 Discussion
Durda & Stern (2000) developed a model of collision rates in the present-day Kuiper belt and Centaur region consisting of a static, multizone, and multi-size-bin collision rate model that calculates instantaneous collision rates on Kuiper belt objects, Centaurs, and the Pluto-Charon system.
The model of Durda & Stern (2000) defines a
disk in the 30-50 AU zone whose surface mass density follows a
profile of the form
,
with
.
The disk adopts an average eccentricity
and assumes an equilibrium condition where the average inclination
.
The values adopted for
are 0.0256
and 0.2048, which are representative of
the values observed in the Kuiper belt at that time.
Once the disk properties are specified, the disk is binned into a
series of radially concentric tori 1 AU in width in each of
which the size distribution of the population is represented by a
two-component power law of the form
,
where
b = -3 for D <
10 km and b = -4.5 for larger bodies.
Durda & Stern (2000)
carried out their runs assuming the existence of
-1.4
105
objects with a radius r > 50 km
between 30 and 50 AU, according to the observational
evidence of Jewitt et al. (1998) and
Gladman et al. (1998).
Moreover, they suggested that about 40% of the total
population are in or near the 3:2 mean motion resonance with Neptune.
From this, the flux of r = 1 km
Kuiper belt projectiles onto Pluto and Charon over 3.5 Gyr is
found to be approximately 8.9
103
and 1.1
103,
respectively.
However, the observational advances that have occurred in the
past ten years have revealed a more complex dynamical structure of the
trans-Neptunian region. In particular, in the region
between 30 and 50 AU, there are very different
dynamical classes of TNOs, such as the classical, resonant, and
scattered-disk populations. They show different dynamical features and
ranges of orbital elements as well-defined semimajor axis zones,
eccentricities as high as 0.45
and inclinations as high as
.
Moreover, current data provide clear evidence of differences in the
mass and size distribution parameters among these dynamical classes
(Kenyon et al. 2008).
The results by Durda & Stern (2000)
are obtained for objects in that region, but they do not account for
the different dynamical populations. In particular, while the Plutino
population is included quantitatively in Durda & Stern's (2000) model, they
did not account for the dynamical properties of this resonant
population. Their results should therefore be taken with caution.
5 Conclusions
We have presented an analysis of the collisional and dynamical evolution of the Plutinos. Assuming the existence of one Pluto-sized object in the 3:2 Neptune resonance, our main results are the following.
Our results depend strongly on the initial size distribution
of Plutinos. Depending on the initial population, we have found that
the number of
D > 1 km Plutinos that have collided
with Pluto in the past 3.5 Gyr is between 1271 and 5552. For
Charon, the number of D > 1 km
Plutino impactors is between 354 and 1545. Moreover, the
number of D > 1 km craters on Pluto
produced by Plutinos in the past 3.5 Gyr is between
43 076 and 113 879.For Charon, the number of D
> 1 km craters is between 20 351 and
50 688.
On the other hand, the largest Plutino impactor onto Pluto has a
diameter of between
17 and
23 km, which produces a crater with a diameter of
31-39 km.
The largest Plutino impactor onto Charon has a diameter of between
10 and
15 km, which produces a crater with a diameter of
24-33 km.
We have tested the dependence of our results on the number of
Pluto-sized objects in the Plutino population. The number of Plutino
impactors on Pluto and Charon obtained by considering two Pluto-sized
objects is a factor of 1.6-1.8
larger than that obtained considering only one Pluto-sized object in
this resonant region.
Using two different QS laws, with a small and a large gap with respect to the QD law, we have obtained almost the same impactor flux onto Pluto and Charon.
The complex structure of the trans-Neptunian region, with its dynamically different populations, requires a detailed study of the contribution of all the potential sources of impactor on the Pluto-Charon system, to obtain the main contributor and the whole production of craters. When Nasa's New Horizon Pluto-Kuiper belt mission flies by Pluto system, we will obtain images of the surface, and cratering records will be able to be measured. We will then be able to compare the theoretical studies of the production of craters with the observations of the New Horizon mission, which will help us to understand the geological processes and history of the surfaces of those worlds. In future work, we will calculate the contribution of the different dynamical classes of the trans-Neptunian region to the flux of projectiles onto Pluto and Charon. From this, it would be possible to specify the primary source of impactors on these bodies and to determine whether the Plutino population can be considered a significant source.
References
- Benz, W., & Asphaug, E. 1999, Icarus, 142, 5 [NASA ADS] [CrossRef] [Google Scholar]
- Brown, M. E. 2008, in The Solar System Beyond Neptune, ed. M. A. Barucci, H. Boehnhardt, D. P. Cruikshank, & A. Morbidelli (Tucson, USA: University of Arizona Press, 335 [Google Scholar]
- Chiang, E., Lithwick, Y., Murray-Clay, R., et al. 2007, in Protostars and Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil (Tucson, USA: University of Arizona Press), 895 [Google Scholar]
- Davis, D. R., Weidenschilling, S. J., Farinella, P., Paolicchi, P., & Binzel, R. P. 1989, in Asteroids II, ed. R. P. Binzel, T. Gehrels, & M. S. Matthews (Tucson, USA: University of Arizona Press), 805 [Google Scholar]
- de Elía, G. C., & Brunini, A. 2007, A&A, 475, 375 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- de Elía, G. C., Brunini, A., & Di Sisto, R. P. 2008, A&A, 490, 835 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Dell'Oro, A., Marzari, F., Paolicchi, P., & Vanzani, V. 2001, A&A, 366, 1053 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Durda, D. D., & Stern, S. A. 2000, Icarus, 145, 220 [Google Scholar]
- Fernández, J. A., Gallardo, T., & Brunini, A. 2002, Icarus, 159, 358 [NASA ADS] [CrossRef] [Google Scholar]
- Gladman, B., Kavelaars, J. J., Nicholson, P. D., Loredo, T. J., & Burns, J. A. 1998, AJ, 116, 2042 [NASA ADS] [CrossRef] [Google Scholar]
- Holsapple, K. A. 1993, Annu. Rev. Earth Planet. Sci., 21, 333 [CrossRef] [Google Scholar]
- Jewitt, D., Luu, J. X., & Trujillo, C. 1998, AJ, 115, 2125 [NASA ADS] [CrossRef] [Google Scholar]
- Kenyon, S. J., Bromley, B. C., O'Brien, D. P., & Davis, D. R. 2008, in The Solar System Beyond Neptune, ed. M. A. Barucci, H. Boehnhardt, D. P. Cruikshank, & A. Morbidelli (Tucson, USA: University of Arizona Press), 293 [Google Scholar]
- Morbidelli, A. 1997, Icarus, 127, 1 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- O'Brien, D. P., & Greenberg, R. 2005, Icarus, 178, 179 [NASA ADS] [CrossRef] [Google Scholar]
- Protopapa, S., Boehnhardt, H., Herbst, T. M., et al. 2008, A&A, 490, 365 [Google Scholar]
- Weissman, P. R., & Stern, S. A. 1994, Icarus, 111, 378 [NASA ADS] [CrossRef] [Google Scholar]
- Zahnle, K., Schenk, P., Levison, H., & Dones, L. 2003, Icarus, 163, 263 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Table 1: Collisional evolution results about Pluto.
Table 2: Collisional evolution results about Charon.
Table 3:
Values of the density ,
the surface gravity g, and the encounter velocity
used for Pluto and Charon.
All Figures
![]() |
Figure 1: Impact strength. The dashed lines represent the two different QS laws used in our simulations. The QD law from Benz & Asphaug (1999) for icy bodies at 3 km s-1 is plotted as a solid line. |
Open with DEXTER | |
In the text |
![]() |
Figure 2: Cumulative number of Plutino impactors onto Pluto a) and Charon b) over the past 3.5 Gyr as a function of impactor diameter. |
Open with DEXTER | |
In the text |
![]() |
Figure 3: Cumulative number of craters on Pluto a) and Charon b) produced by Plutinos over the past 3.5 Gyr as a function of crater diameter. |
Open with DEXTER | |
In the text |
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.