Issue |
A&A
Volume 510, February 2010
|
|
---|---|---|
Article Number | A43 | |
Number of page(s) | 12 | |
Section | Planets and planetary systems | |
DOI | https://doi.org/10.1051/0004-6361/200913322 | |
Published online | 05 February 2010 |
SDSS-based taxonomic classification and orbital distribution of main belt asteroids
J. M. Carvano1 - P. H. Hasselmann1,2 - D. Lazzaro1 - T. Mothé-Diniz2
1 - Observatório Nacional (COAA), rua Gal. José Cristino 77,
São Cristóvão, CEP 20921-400 Rio de Janeiro RJ, Brazil
2 -
Universidade Federal do Rio de Janeiro/Observatório do Valongo,
Lad.Pedro Antônio, 43, 20080-090 Rio de Janeiro, Brazil
Received 18 September 2009 / Accepted 4 November 2009
Abstract
Aims. The present paper aims to derive a new classification
scheme for SDSS MOC asteroid colors that is compatible with previous
taxonomies based on spectroscopic data. The distribution of these can
give important clues to the formation and evolution of this region of
the Solar System, as well as to locate candidates with mineralogically
interesting spectra for detailed observations.
Methods. The methodology is based on the large database SDSS
MOC4. Templates of the main taxonomic classes are derived and then used
to classify the asteroid observations in the SDSS MOC4. The derived
taxonomic scheme is compatible with the Bus taxonomy and is suitable to
the peculiarities of the SDSS observations, in particular, the low
spectral resolution.
Results. Density maps of the seven classes defined by the method
reproduce classical results for the background which is mainly
dominated by the
class in the inner belt and by the
and the
classes beyond 2.8 AU. It also shows new structures, such as the fact that the
and
seem evenly distributed in the inner belt while in the outer belt the
class increase in density only at the location of asteroid families. Although their overall distribution is similar, the
class seems less frequently associated with large families than do the
class asteroids. Although only clustering around the Vesta family, the
class asteroidsnappear scattered all around the main belt. Besides the
lack satisfactory explanations of most of the highlighted features,
they may provide strong constraints on the models of the formation and
evolution of the Solar System.
Key words: minor planets, asteroids: general - methods: miscellaneous
1 Introduction
Since the early works in asteroid taxonomy, there has been a tendency to establish firm relations between the derived taxonomic classes and mineralogy. At that time, the only analytical tool available to researchers were direct comparison between the visible reflectance spectra of asteroids and those of meteorite and mineral samples obtained in laboratory. Indeed, by using such methodology McCord et al. (1970) were able to establish a very robust relation between the asteroid (4) Vesta and the howardite, eucrite, and diogenite (HED) meteorites. Possibly inspired by this accomplishment, in the early taxonomies (Chapman et al. 1975; Bowell et al. 1978) the classes were denominated by letters according to the mineralogical interpretation of their observational properties. In this way, the C class was associated to carbonaceous chondrites, the S class to stony-iron meteorites, the E class to enstatite achondrites, and the M class to the metallic iron-nickel meteorites. Later taxonomies (Tholen 1984; Bus & Binzel 2002a; Barucci et al. 1987; Tedesco et al. 1989) kept most of the class denominations. With the increase in the number of observations available and the improvement in classification algorithms, the number of classes increased between the first and the second generation of taxonomies, culminating in the Bus taxonomy (Bus & Binzel 2002a), which defines 26 classes.
The direct mineralogical characterization of the asteroids, however,
has proved to be much harder to achieve. It is severely complicated by
(Gaffey et al. 2002): 1) the existence of several factors
(surface texture, observation geometry, space weathering effects) that
are not directly related to mineralogy but may affects the observed
reflectance spectra of asteroids; 2) the fact that not all
mineralogical assemblages exhibit diagnostic features in their
spectra; and 3) most of the diagnostic spectral bands of the most
abundant minerals tend to be broad and to overlap heavily. Despite the
difficulties, considerable progress has been made in this area, some
of which have shown how uncertain the mineralogical inferences made
from taxonomy can be. A good example here is the discovery that some
E- and M-class asteroids show hydration bands around
(Jones et al. 1990; Rivkin et al. 1995,2000). The
presence of this band indicate that the asteroids did not suffer
severe heating during their evolution, which is clearly inconsistent
with the thermal history of enstatite achondrites and nickel-iron
meteorites. In the face of the evolution of our understanding of
the mineralogy of asteroids, it is becoming evident that taxonomy
alone is a poor tool for infering composition.
These limitations do not mean, however, that asteroid taxonomies are of no use. The fact that some classes have been, with minor modifications, systematically recovered by independent authors using diverse data sets and different methodologies suggests that they convey very real information about the diversity of spectral properties of asteroids. Taxonomic classification is paramount to selecting targets for detailed investigation. Also, given the scarcity of reliable mineralogic characterization of asteroids, the use of taxonomic classes is still the best available way to study the large-scale distribution of compositions in the main asteroid belt. And, although taxonomy alone is not sufficient to make meaningful mineralogical inferences, the analysis of the classes distribution in the main belt can, at least, highlight regions and objects that deserve further investigation.
Given these two main applications of taxonomy - target selection and statistical studies - it is fair to say that the usefulness of a given taxonomic scheme is proportional to the number of objects it potentially can classify. This depends on several factors, among which are the number of suitable observations that are available at any given time. The Bus taxonomy represented an increase of about an order of magnitude over early efforts. Here taxonomies based on data of large photometric surveys, like the Sloan Digital Sky Survey (SDSS), offer a clear advantage, since they can provide observations for much more objects than what can be reasonably obtained from spectroscopy. Several works have been using taxonomy classification inferred from the SDSS colors, but so far without formal linkage to previous taxonomies. Such work, however, is needed to provide a robust link to previous studies that use those taxonomies. Also, it will be useful as a basis of comparison to the next logical step, which is to derive an independent taxonomic system based on the SDSS data, which may reveal new putative spectral classes that could have eluded previous efforts based on smaller datasets.
In the present work we define a classification algorithm based on the SDSS colors. It seeks to be compatible with previous taxonomic schemes based on spectroscopic data. The goal here is not to define a new taxonomy, but merely to provide a concise method of assigning taxonomic classes to SDSS observations of asteroids that are compatible with the limitations of the dataset. This scheme is then applied to the asteroids listed in the fourth release of the SDSS Moving Object Catalogue (MOC4) and to meteorite spectra from the RELAB database. These results are used to analyze the distributions of taxonomies in the main asteroid belt. We also use the large number of asteroids that were observed multiple times in the MOC4 to analyze the issue of taxonomic variations among observations. Since there are very few definite conclusions about mineralogy that can be drawn based solely on taxonomy, the ultimate goal here is to identify interesting targets for NIR spectroscopic follow-up observations and to study the distribution of taxonomies in families and in the main belt as a whole.
2 Methodology
2.1 Parameterization of the observations
The definition of a taxonomic scheme depends on the choice of the
observational parameters and of the methodology used for separating the
observations among classes.
The Sloan Digital Sky Survey uses the five filters that define the
u'g'r'i'z' system. These filters have central wavelengths equal to,
respectively, 0.354, 0.477, 0.6230, 0.7630 and
.
This filter set was designed for minimum superposition and
together cover a spectral interval from 0.3 to
(Fukugita et al. 1996), which is compatible with the
spectral interval
used in most of the previous taxonomic classifications. To obtain
information on the reflectance of the asteroids, it is necessary first
to remove the solar contribution in their observed magnitudes. This is
usually done by subtracting the solar colors from
those obtained from the observed asteroid magnitudes, which
implicitly results also on the normalization of the spectral
reflectance distribution of the observation at the band used as the
reference for calculating the colors. The result of this operation can
be called reflectance colors:
![]() |
(1) |
where



Normalization was also a necessary step for reflectance measurements obtained through visible spectroscopy, since i) the spectra were not usually flux-calibrated and ii) the observed fluxes depend both on the reflectance of the surface (as a function of wavelength) and on the total illuminated area of the asteroid that was visible at the instant of the observation, which is usually unknown. The latter reason still applies to the flux-calibrated SDSS magnitudes.
After removal of the solar contribution, these reflectance colors can be used directly to forge a taxonomic scheme (Binzel et al. 2007) or can first be converted to reflectance spectra (Roig & Gil-Hutton 2006; Misra & Bus 2008), which is more directly comparable to what is obtained through low-resolution spectroscopy. In either case, each observation ends up parameterized by four quantities (since the reflectance at the reference filter becomes redundant) whose values depend on the choice of the normalization point.
The choice of the normalization wavelength is usually done so as to avoid
regions with spectral features, in order to minimize the effects of
normalization when comparing spectra. A common choice for low-resolution
spectra is to normalize it at the central wavelength of the Johnson V band
(
), while most authors working with SDSS asteroid observations
choose to normalize at the r' band (whose central wavelength is
). It has been recently realized, however, (Carvano et al., in preparation)
that normalization can affect spectral comparison even if made at these
relatively featureless spectral regions. To avoid this potential
pitfall we parameterize each observation by its reflectance color gradients,
which we define as
![]() |
(2) |
These four quantities, which are independent of the normalization wavelength, represent the spectral inclination of the logarithm of the reflectance between contiguous filters. Therefore, a positive value indicates an increase in the reflectance, a negative indicates decrease in the reflectance and zero indicates that the reflectance do not varies between the filters in consideration. The reflectance colors for the u' filter can be readily calculated from the SDSS magnitudes as
![]() |
(3) |
where


2.2 Classification algorithm
Once the data parameterization is defined there are several possible paths to a supervised classification scheme. Binzel et al. (2007) use previously classified asteroids to define templates for each class and then classify the remaining observations based on the distance to the templates of each class, calculated with a given metric; Roig & Gil-Hutton (2006) perform a principal component analysis on the SDSS reflectance spectra and then use previously classified asteroids to define the locus of each class in principal component space. Misra & Bus (2008) use previously classified asteroids and SDSS observations classified by an human operator to train a neural network that, in turn, classify the remaining SDSS observations. We believe, however, that the existing schemes do not contemplate an important characteristics of the SDSS asteroid observations: that the classification can be very sensitive to magnitude uncertainties when compared to classifications made from spectroscopy.
This issue arises from the low spectral resolution of the SDSS observations. Although each observation is flux-calibrated, hence better corrected for atmospheric extinction and systematic effects than most asteroid spectra obtained in surveys, the fact that each reflectance spectra is sampled at only 5 wavelengths means that the uncertainties in the measurements can substantially change the shape of the resulting spectra. This induces a fuzziness into the classification procedure that, we believe, must be made explicit by the classification algorithm to allow informed decisions by anyone that uses the taxonomy for target selection and/or statistical studies.
To satisfy this requirement we devised the following scheme: initially, limits
on the values of the four reflectance color gradients are calculated for each
class, using observations of previously classified asteroids, meteorite
spectra, or synthetic spectra (see Sect. 2.3). To classify one
SDSS observation of an asteroid, we calculated the probability distribution of
each reflectance color gradient as a normal distribution with mean equal to
and standard deviation equal to the uncertainty of the
measurement. The probability that the measurement in a given spectral
interval is within the limits of a class is then the area of the density
distribution that falls within those limits. The probability that an
observation is compatible with a class is then the product of the
probabilities at each spectral range. Because of the overlapping of the
limits of the classes, the probabilities of the observation being
compatible with
different classes will not in general sum up to unity. The best
classification for the observation is chosen as the one with at least 60% of
the sum of the probabilities for all possible classes. The final
classification is then presented as a letter for the class plus the
corresponding probability. If no class satisfy
this requirement (meaning that the probabilities for the observation to be
compatible with two or more classes are comparable), all classes with
probabilities over 30% of the total probability are considered
possible and a multiple letter classification is issued, with the
corresponding probability being the average value for the possible classes.
The order of the letters in multiple classification has no
relation with the probability of each class. If an observation fails to be
assigned to any class, it is labeled as unclassified (``U'').
One consequence of this classification scheme is that observations with high observational uncertainties in one or more filters will be classified with low probabilities. These classifications will then be easily recognized as problematic, but could still be useful in some contexts.
Other distinctive characteristic of the SDSS MOC4 is that it includes several asteroids with multiple observations, which increases the possibility of having asteroids with multiple classifications. Since a unique classification for each asteroid is needed for statistical studies, whenever there is variation in the taxonomic class with different observations, we adopt for the asteroid the most frequent taxonomy, with probability equal to the higher value for that class. If two classes are equally frequent, a multiple letter classification is issued, with probability equal to the mean of the higher probabilities for each class, and a flag is set to indicate that this situation has happened.
2.3 Definition of the classes
A final step is to choose which classes of the spectroscopic taxonomies
to represent and how to define
their limits. Another problem here is the nomenclature: the norm in
the literature is for new taxonomic schemes to use the same designations for
their classes that were used in previous taxonomies, even though they usually
differ in the way they are defined. Here we append a ``p'' to the
corresponding designation in the Bus taxonomy, to indicate clearly that this
taxonomy is obtained from photometric data. We decided to define 9 classes,
based on their current
mineralogical interpretation and the ability of the SDSS data to represent
them. The adopted classes are ,
,
,
,
,
,
,
,
and
.
Except for the
class, all others are intended to be directly comparable to the corresponding
classes in Bus taxonomy. It should be noted that, although
this classification scheme was designed to be compatible with the Bus
taxonomy, such compatibility is limited by the characteristics of the
dataset.
The first five are characterized by an absorption band around
,
associated with olivine/pyroxene assemblages. On the SDSS data,
this band causes as a drop in reflectance at the z' filter. The Vclass characteristic spectrum is distinctive by having the deepest
band among observed asteroid spectra. This class is firmly
associated with HED meteorites (McCord et al. 1970), and is
observed mostly in and around the Vesta family
(Alvarez-Candal et al. 2006; Binzel & Xu 1993; Burbine et al. 2001; Florczak et al. 2002)
in the inner main belt, with a few observed in the intermediate
(Nesvorný et al. 2008; Binzel et al. 2006; Roig et al. 2008; Moskovitz et al. 2008)
and outer belt (Duffard & Roig 2009; Lazzaro et al. 2000). The O
and Q classes are characterized by a not very deep
band and
by a flatter spectrum shortwards the i' band than the V class (with
the Q class being more reddened in the region than the O class). Both
have been associated with ordinary chondrite meteorites
(McFadden et al. 1984; Binzel et al. 1993; Bell & Keil 1988; McFadden et al. 1985),
since their visible spectra is similar to spectra of this type of
meteorite (Burbine & Binzel 2002).
Direct comparison and mineralogical inferences based on NIR spectra
tend to confirm this association for the Q class
(Binzel et al. 2001), but the situation is not as clear for the
O class
(Chapman 1996a; Burbine & Binzel 2002; Binzel et al. 2001; DeMeo et al. 2009). Both
classes are rare among asteroids, with Q-class asteroids observed
mostly among the NEO population
(Binzel et al. 2004; Mothé-Diniz & Nesvorný 2008). The S class, on the
other hand, is the most common for asteroids in the inner and
intermediate main belt (Mothé-Diniz et al. 2003). This class harbors
a good deal of spectral diversity, and in the Bus taxonomy
(Bus & Binzel 2002a) it was deemed a ``complex'' composed of 6
subtypes. Spectra of asteroids with this classification tend to be
redder than the Q-class shortward of the i' band, present a
shallower
band, and, in general, are not well-matched by
most meteorite spectra. Mineralogical inferences made based on the
position and area of the 1 and
olivine/pyroxene and
pyroxene bands suggest that a considerable mineralogical diversity can
be found among the asteroids with this classification, including
ordinary chondrites, mixtures of metal and basalts and
olivine-dominated assemblages (Gaffey et al. 1993). Further
evidence from laboratory experiments
(Sasaki et al. 2001; Hapke 2001) and spacecraft
observations (McFadden et al. 1985) suggest that at least part of
the S-class asteroids may have a superficial composition compatible
with ordinary chondrite material affected by space weathering
processes.
Finally, A-class asteroids present extremely reddened spectra
shortwards the i' band and a moderate to deep
band. These
asteroids are rare in the main belt, and their characteristic spectrum
is not well-matched to spectra of meteorites
(Burbine 2000). Analysis from NIR spectra shows that
A-class asteroids have compositions dominated by olivine, with varying
fosterite content. These extremes may represent mantle material from
a completely differentiated parent body or material similar in
composition to R chondrites (Sunshine et al. 2007).
The D, X, and C classes show no clear
band, and their
spectra longwards of the g' band is mostly featureless. Their main
defining characteristic is their spectral gradient in this region,
with the D class presenting a very steep gradient, the X class a
moderate gradient, and the C class a flat to slightly bluish gradient
(Bus & Binzel 2002a). Given the absence of mineralogically
diagnostic bands in most of the cases, mineralogical inferences for
these classes rely mostly on comparison with meteorite spectra and
synthetic spectra generated from mixing models. The D-type asteroids
are generally believed to be mixtures of ices, anhydrous silicates,
and organics. Their characteristic spectra is mostly not seem among
meteorites with the exception of a slight similarity with some iron
meteorites (Carvano et al. 2003) and the Tagish Lake meteroite
(Hiroi et al. 2003), which is an ungrouped carbonaceous
chondrite. Also, the D-class are the predominant class in the far end
of the outer main belt and in the Trojan region, and their spectra are
similar to spectra of part of the TNOs and comet nuclei
(Burbine 2000). The X and C classes, on the other hand,
show some spectral diversity and were also deemed ``complexes'' in the
Bus taxonomy, each subdivided into several subclasses. Comparison to
meteorites (Burbine 2000) suggest that their possible
mineralogies include all carbonaceous chondrites subtypes,
shock-darkened ordinary chondrites, enstatite chondrites, iron
meteorites, and some achondrites. Also, spectra similar to both
classes can be seen among TNOs and comet nuclei
(Jewitt 2002), so that volatile and organic rich
compositions are also a possibility. Notable exceptions to this
general indetermination is the Ch Bus subclass, characterized by a
wide absorption band around
that is associated to water
alteration in mafic phyllosilicates and that is present on some CM
carbonaceous chondrites (Vilas & Gaffey 1989).
Finally, the
class, as defined in this work, differs
from the other because it is not directly comparable to classes
defined in previous works. Rather, it is defined by the junction of
the K, L, and Ld classes from Bus taxonomy, and for this reason, it is
represented here by a different typesetting. These three classes have
spectra that, shortwards of the i' band, vary from moderately (K) to
extremely reddened (Ld), all presenting a flattening of the spectra
longward of that band. This last characteristic may stem from the
presence of a very shallow
band, and in this sense it
represents a transition between the classes with and without the band.
2.4 Defining the limits of the classes
A first step is to select the objects that have spectroscopic classification in the Bus system and that have also been observed by the SDSS. We use here the objects observed by the SMASSII (Bus & Binzel 2002b) and S3OS2 surveys (Lazzaro et al. 2004). Since we want these objects to define each class, it is important to exclude object whose classification is known to change between observations. The issue of spectral variations among the two spectroscopic surveys is discussed in Carvano et al. (2003) and in Lazzaro et al. (2004). Using the compilation of taxonomic classification of asteroids (Neese 2006), we looked for objects that have classifications in the classes that form the classes we adopt here. If an object was observed in both spectroscopic surveys, it was only used if its classifications according to both fall in the same class. We then examined the distribution of the uncertainties in the color gradients of the whole sample: observations where the errors were greater than the 3rd quartile for any of the filters were discarded at this stage. The result of this search is shown in Table 1. In this table, the third and fourth columns the number in parenthesis refer to the number of objects and their observations that remained after removing the observations with uncertainties above the adopted upper limit.
Table 1: Number of objects with unambiguous spectral classification that were observed by the SDSS.
We note that no A-, Q-, and O-type asteroids with unambiguous
taxonomy have been observed by the SDSS. Q- and O- types are thought
to represent the unweathered end-members of the population of
asteroids with ordinary chondrite composition. It is then reasonable
to use meteorite spectra to obtain representative SDSS reflectances
for these classes. To do so, we used the templates of the Q and O
classes in Bus taxonomy to search the RELAB database for OC meteorites
with a spectrum similar to each class, with the methodology described
in Mothé-Diniz & Carvano (2005). SDSS reflectances for each meteorite
were then obtained by convolution of the meteorite spectra with the
SDSS band passes for each filter. This approach is not possible for
A-type asteroids, since they have no clear analogs among the meteorite
collection. However, it is well established that A-type asteroids are
olivine-dominated assemblages (Sunshine et al. 2007), with
compositions ranging from magnesium-rich to iron-rich. Therefore, to
obtain reflectances in the SDSS range we adopted the following
approach: for each of the 16 unambiguous A-type visible spectra, we
used the implementation of the Modified Gaussian Model
(Sunshine et al. 1990) described in Mothé-Diniz et al. (2008) to
fit olivine bands. Since the visible spectrum alone is not sufficient
to allow a unique fit, the three olivine bands were fixed to a given
fosterite composition by both the band center calibration and the
band-depth ratio calibration (Sunshine et al. 2007) . For each
spectrum, we then produced five synthetic spectra by fitting with
compositions
,
,
,
,
and
.
These synthetic spectra were then convolved with the SDSS
band passes to yield the SDSS reflectance spectra for the A-class.
To check whether the reflectance gradient ranges of the template spectra of each class are able to be used as classification criteria we inspected for superpositions. We did that by calculating the kernel density plot for each class at each spectral range. This is done in the following fashion: the density of probability of a given observation of an object of a given class having a reflectance gradient at a given value was calculated as a normal distribution, with the mean equal to the reflectance gradient of the observation and the variance equal to the uncertainty in the reflectance gradient. The density distribution of the class in that range was then given by the sum of the densities of all observations of asteroids of that class. To use the proposed scheme for classification, it was necessary that each class presents little superposition with all the other in at least one spectral interval. However, since the classes represent a continuous distribution, some degree superposition is unavoidable in some instances.
It is found that there is considerable overlapping between the ,
and
classes, in particular between the
and
the
classes and between the
and the
classes. Upon visual inspection of the spectra, we see that
several asteroids classified as S-type present in some observations
SDSS reflectance spectra, which are closer to templates of the L and K
classes, with a flattening of the reflectance towards
.
Conversely, many SDSS observations of asteroids in the
class show spectra more similar to S-type asteroids,
with a decrease in the reflectance towards
.
Similarly, the
overlapping between the
and
classes happens
mostly because of the Ld-class asteroids in the latter. Both Ld- and
D-type spectra have steep visible spectral gradients, but for the
Ld-class there is a flattening of the reflectance towards
,
whereas the reflectance of the D-type asteroids shows a continued
increase in that range.
A solution is to seting strict limits for reflectance gradients of
the ,
,
and
classes templates in the (z-i)
range. A natural choice of theses limits are the values of the
reflectance gradient where the density distribution of the classes
intercept each other. Observations of asteroids classified as the
,
and
classes by spectroscopic surveys that
fall outside these limits are not used to define the classes from SDSS
reflectance color gradients.
The resulting densities are then used to define the final limits of the
classes. The limits of each class are, at this step, given by the
region in each spectral range that contains a
probability of containing
a template of that class. Figure 1 shows schematically
the limits of the classes in reflectance gradient space. The log reflectance
spectra corresponding to the peak of the distributions for each class is shown
in Fig. 2.
![]() |
Figure 1: Final reflectance color gradient limits for the classes. |
Open with DEXTER |
![]() |
Figure 2: Median log reflectance spectra for the classes, shifted vertically for clarity. |
Open with DEXTER |
![]() |
Figure 3: Distribution of the SDSS taxonomic classes for observations of the asteroids previously classified based on spectroscopy. The solid filling represents the number of observations with uncertainties above the threshold set in Sect. 2.4. |
Open with DEXTER |
3 Results
3.1 Comparison with spectroscopic classifications
It is instructive to apply this classification scheme first to the sample of asteroids with unambiguous classification. Here we use the whole sample, with no exclusions based on the uncertainties. Figure 3 shows the distribution of the classification made from the SDSS observations of the objects with spectroscopic classification for each class.
In this figure, observations for which the higher probability was lower than 5% were considered to be unclassified (U). Also, the amount of observations in each class made from observations considered bad are shown separately.
In general, visual inspection of the spectra shows that they are
consistent with the assigned class. Exceptions are the spectroscopic C
asteroids that have observations classified as
or
,
which
would fit in the
or
classes better. Most of the
observations that failed to be classified have high
uncertainties. Visual inspection of the observations assigned to the U
class that have uncertainties below the threshold adopted show that
all are definitively ``strange''.
3.2 Application to meteorite spectra
![]() |
Figure 4: Distribution of the SDSS taxonomic classes for meteorite spectra from the RELAB database. |
Open with DEXTER |
![]() |
Figure 5: Distribution of the SDSS taxonomic classes for meteorite spectra from the RELAB database. |
Open with DEXTER |
Direct comparison between meteorite and asteroid spectra is hampered by the effects of space weathering processes acting on the surface of the asteroids. This effect, together with the selectiveness of the transport mechanisms of meteoroid are currently understood as the main culprits for the observed mismatch between the spectral properties of asteroids and meteorites (Gaffey et al. 2002). Even so, such an exercise can provide some insight to the mineralogical interpretation of the taxonomic classes. To this end, we calculated the SDSS color (see Sect. 2.4) and applied the classification algorithm derived in the previous sessions to a sample of bidirectional reflectance spectra of several meteorite classes from the RELAB spectroscopic database (Pieters & Hiroi 2004). For each spectral interval, we used the median uncertainties from the sample of asteroids with previous spectroscopic classification that were observed by the SDSS. This essentially tells us what would be the taxonomic classification of asteroids covered by unweathered material of each meteorite type.
The results are shown as histograms in Figs. 4 and 5. From these figures, it is clear that none of the meteorite
types is spectroscopically homogeneous, with specimens falling into different
SDSS taxonomic classes. There is also many of spectra that fail to
be classified. Visual inspection reveals that the main culprit for most of the
unclassified meteorite spectra is the steep decrease in the reflectance
towards the UV that is commonly seen in meteorite spectra but rare among
asteroids (Britt et al. 1992; Carvano et al. 2003; Burbine 2000). This characteristic is also responsible for the
classification of some carbonaceous chondrite and ordinary chondrite spectra
in the
class.
As expected, most ordinary chondrites fall into the
or
classes, but not exclusively: the spectra of some H and L chondrites
with unusually deep
band end up classified as
-class. Also, a number of spectra with unusually swallow
band are assigned to the
or
classes. These are mostly
spectra of shock-darkened ordinary chondrites, also known as black
chondrites (Britt & Pieters 1989). Also as expected, the
majority of HED spectra that end up classified were assigned to the
class, but a few spectra with shallower
fell into the
and
classes. Besides the UV drop in reflectance, the large
number of HED spectra that could not be classified is also a
consequence of a
that is deeper than what is seen among
V-class asteroid spectra in Bus taxonomy. Most spectra of carbonaceous
chondrites subtypes fall into the
class, although a notable
exception are the CO carbonaceous chondrites that show more affinity
for the
,
and
classes. Also
predominantely classified as
-class are the spectra of
mesosiderites, pallasites and iron meteorites. The
-class spectra
are seen among CI and CM and (just one case) L ordinary chondrite and
also among iron meteorites. It is interesting to note that, apart from
2 L chondrites, 2 urellites and 1 HED spectra, spectra compatible with
the
class, are not seen among meteorites, which agrees with the
current view that the S-type spectra are a consequence of space
weathering processes acting on ordinary chondrite (Chapman 1996b)
and possibly other mineralogies.
3.3 Application to the MOC4
We proceed then to apply the classification scheme to all observations of asteroids in the MOC4. Here we do not exclude any observation a priori, but only flag as bad the ones with uncertainties above the threshold. The MOC4 contains 218 996 observations of 104 687 asteroids. Out of these, a total of 107 466 observations of 63 468 asteroids are classified in this fashion. Out of the 111 530 observations that failed to be classified, 109 559 had one or more filters with uncertainties above the threshold, while 87 002 observations flagged as bad were still classified with probability greater than 5%. Table 2 shows the number of asteroids and observations classified into each of the adopted classes and also into some of the most common ambiguous classifications. As in Table 1, the number inside parenthesis refer to the asteroids and observations with uncertainties below the threshold - for asteroids this is the number of objects whose best classified observation have uncertaities above the threshold.
Table 2: Number of asteroids and observations classified.
Table 3: Objects with taxonomic variation between observations.
3.3.1 Taxonomic variations
Color variation among observations of asteroids in the SDSS have been analyzed by (Szabó et al. 2004) for a sample of 7531 asteroids in the third release of the MOC. The authors conclude that the observed variations are real and caused by surface inhomogeneities in the asteroids, rather than by observational artifacts or random errors. Here we use a much bigger sample to analyze color variations that lead to different taxonomic classifications between observations.
A total of 22 019 asteroids have two or more observations with
successful classifications. Of these, 14 962 show taxonomic
variations. Table 3 shows the number of asteroids
associated with each combination of taxonomies found, for the
combinations that are seen on 100 asteroids or more. As should be
expected, most of the variations seen are between classes with
adjacent limits in one or more spectral intervals. This is the case of
the most frequent transitions, from
to
,
from
to
and from
to
,
and the transitions from
to
and from
to
,
from
to
and
to
also fall into this category. Such variations can
be caused by random errors for the objects with spectra falling close
to the classes limits, but they can also reflect real surface
variations, caused either by variations in the prevailing mineralogy
or variations in the prevailing degree of space weathering maturity.
Alternatively, some of the variations can also be caused by
differences in the solar phase angle among observations and
phase-reddening effects.
More interesting are the transitions between more dissimilar classes,
to
,
to
,
to
,
to
and
to
.
These transitions imply that those asteroids at some
time presented spectra dominated by olivine/pyroxene bands that become
suppressed at some other point of their orbit and/or rotational
aspect. Possible mechanisms for accomplishing this include extreme
space weathering, contamination by metal, or other spectrally dark
phase, the coexistence of different mineralogies in the same body or
phase reddening. Again, it is not possible to rule out random error as
the cause of the taxonomic variation in some instances.
![]() |
Figure 6: Density map of all numbered asteroids. The location of several families are marked. |
Open with DEXTER |
3.3.2 Distribution of taxonomies in the main belt
The distribution of taxonomies in the main belt can be highlighted
through the analysis of the density maps for each class. To construct
such maps we simply divide the main belt in bins of osculating
semi-major axis and inclination and then count the number of objects
of a given class inside each bin. Here we adopted bins with sides of
in semi-major axis and
in inclination. The
usefulness of this kind of plot for analyzing the distribution of
asteroids can be appreciated from Figure 6, which
shows a density map made from all asteroids listed in the ASTORB file
(Bowell 2009) for semi-major axes between 1.8 and
and inclinations up to
.
The patchy nature of the
asteroid distribution in the main belt is clearly seen, as is the
localized increase in density in the vicinity of the major asteroid
families. The loci in osculating element space of the members of
the families of Merxia, Veritas, Nysa, Adeona, Padua, Koronis,
Maria/Renate, Mineva, Eunomia, Vesta, Erigone, Dora, Agnia, Eos,
Hygiea, Baptistina, Meliboea, Hoffmeister, Massalia, and Themis (as
defined in Mothé-Diniz et al. 2005) are delimited by closed lines
in Fig. 6.
Figure 7 shows the density distributions for the classes
,
,
,
,
,
,
and
,
for all asteroids with classification probability
.
In general,
limiting the sample based on the uncertaities of the measurements
produce qualitatively similar results to what is shown in
Fig. 7.
![]() |
Figure 7: Density maps for each photometric class. |
Open with DEXTER |
This figure summarizes several results for the distribution of the taxonomies in the main belt that have been already mentioned in the literature, but also reveals some features that have not been previously discussed. The most interesting features can be summarized as follows.
- -
- The
and, particularly, the
classes seem to dominate the background beyond
until the outer border of the main belt. In this region the
-class asteroids mostly appear in the vicinities of the Koronis and Eos families. This predominance of the
and
classes in the background becomes even more explicit at inclinations greater than
.
- -
- The
-class asteroids, on the other hand, dominate the region before
, showing a considerable increase in density even in regions not associated with dynamical families. In this region the increase in the density of
-class asteroids only occur in regions associated to dynamical families; however, unlike the
-class in the outer belt that are rare and patchy outside families, the background
- and
-class asteroids in the inner belt seem evenly distributed.
- -
- Although their overall distribution is similar, the
-class seems less frequently associated with large families than
-class asteroids. The most notable exception occurs in the vicinity of the Eos family.
- -
-class asteroids dominate beyond
. In the Trojan region, in particular, only
-class asteroids are observed with inclinations greater than
. In the main belt,
-class asteroids seem evenly distributed, with two notable exceptions: they are not seen in the outer belt with inclinations greater than
, and they seem to cluster in the region of the Eos family.
- -
- The distribution of the
-,
-, and
-classes seem to follow the same pattern as the
-type asteroids, with two noticeable exceptions:
-class asteroids show a clustering in the vicinity of the Vesta family that is not evident in the
-class asteroids; the
-class asteroids show a clustering in the vicinity of the Eos family that is clearly more intense than what is seen for the
-class asteroids and also seem to cluster in the vicinity of the Meliboea family, while no clear clustering of
-class asteroids can be seen in this region.
- -
- The
-class asteroid only show a significant clustering in the vicinity of the Vesta family, but are scattered in the intermediate and outer belts.


4 Conclusions
Here we have developed a taxonomic scheme that is compatible with the Bus taxonomy and that is suitable to the peculiarities of the SDSS observations. In particular, it explicitly takes into account the uncertainty in the classification introduced by the low spectral resolution of the SDSS data by assigning probabilities to each classification. In this way, it is possible to maximize the number of classified objects while allowing the putative users to select the confidence level of the classification that is adequate to their particular needs. The application of the scheme to the fourth release of the SDSS Moving Object Catalogue yields the classification of 107 466 observations of 63 468 asteroids with a probability over 5%.
Applying of this taxonomy to meteorite spectra showed that objects with a given mineralogy can be spread through different taxonomic classes. In turn, each taxonomic class is, in general, compatible with several different meteorite classes. This exercise again shows that even in the absence of processes such as space weathering, the information concerning the mineralogy of an asteroid that can be obtained only by taxonomy is very limited.
We also used asteroids with multiple observations to verify the issue of taxonomic variability. It was found that color variations that are sufficient to lead to varying taxonomic classifications between observations of the same body are fairly common. Most of these variations may be the result of random errors and the low spectral resolution of the SDSS data, but some of the most extreme could be real and deserve further attention. Finally, the large number of classified objects allow a much more detailed and robust view of the distribution of taxonomic classes in the main belt. Most of the features highlighted in the previous section lack satisfactory explanations. Although some might be related to relatively recent events, like family formation, some others might contain clues to processes that took place in the early stages of proto-planetary disk and may provide valuable constraints to the models of formation and evolution of the Solar System. Again, the taxonomic classification provided in this work can be instrumental in these studies.
AcknowledgementsThe authors acknowledge FAPERJ and CNPq for several grants and fellowships.
References
- Alvarez-Candal, A., Duffard, R., Lazzaro, D., & Michtchenko, T. 2006, A&A, 459, 969 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Barucci, M. A., Capria, M. T., Coradini, A., & Fulchignoni, M. 1987, Icarus, 72, 304 [NASA ADS] [CrossRef] [Google Scholar]
- Bell, J. E., & Keil, K. 1988, in Lunar and Planetary Science Conference Proceedings, ed. G. Ryder, 18, 573 [Google Scholar]
- Binzel, R. P., & Xu, S. 1993, Science, 260, 186 [NASA ADS] [CrossRef] [Google Scholar]
- Binzel, R. P., Xu, S., Bus, S. J., et al. 1993, Science, 262, 1541 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Binzel, R. P., Harris, A. W., Bus, S. J., & Burbine, T. H. 2001, Icarus, 151, 139 [NASA ADS] [CrossRef] [Google Scholar]
- Binzel, R. P., Rivkin, A. S., Stuart, J. S., et al. 2004, Icarus, 170, 259 [NASA ADS] [CrossRef] [Google Scholar]
- Binzel, R. P., Masi, G., & Foglia, S. 2006, in BAAS, 38, 627 [Google Scholar]
- Binzel, R. P., Masi, G., Foglia, S., et al. 2007, in Lunar and Planetary Institute Science Conference Abstracts, 38, 1851 [Google Scholar]
- Bowell, E. 2009, VizieR Online Data Catalog, 1, 2001 [Google Scholar]
- Bowell, E., Chapman, C. R., Gradie, J. C., & Zellner, B. 1978, Icarus, 35, 313 [NASA ADS] [CrossRef] [Google Scholar]
- Britt, D. T., & Pieters, C. M. 1989, Meteoritics, 24, 255 [NASA ADS] [Google Scholar]
- Britt, D. T., Tholen, D. J., Bell, J. F., & Pieters, C. M. 1992, Icarus, 99, 153 [NASA ADS] [CrossRef] [Google Scholar]
- Burbine, T. H. 2000, Ph.D. Thesis, MIT [Google Scholar]
- Burbine, T. H., & Binzel, R. P. 2002, Icarus, 159, 468 [NASA ADS] [CrossRef] [Google Scholar]
- Burbine, T. H., Buchanan, P. C., Binzel, R. P., et al. 2001, Meteoritics and Planetary Science, 36, 761 [NASA ADS] [CrossRef] [Google Scholar]
- Bus, S. J. 1999, Ph.D. Thesis, AA Massachusetts Institute of Technology [Google Scholar]
- Bus, S. J., & Binzel, R. P. 2002a, Icarus, 158, 146 [NASA ADS] [CrossRef] [Google Scholar]
- Bus, S. J., & Binzel, R. P. 2002b, Icarus, 158, 106 [CrossRef] [Google Scholar]
- Carvano, J. M., Mothé-Diniz, T., & Lazzaro, D. 2003, Icarus, 161, 356 [NASA ADS] [CrossRef] [Google Scholar]
- Cellino, A., Bus, S. J., Doressoundiram, A., & Lazzaro, D. 2002, Asteroids III, 633 [Google Scholar]
- Chapman, C. R. 1996a, Meteoritics and Planetary Science, 31, 699 [Google Scholar]
- Chapman, C. R. 1996b, Meteoritics and Planetary Science, 31, 699 [Google Scholar]
- Chapman, C. R., Morrison, D., & Zellner, B. 1975, Icarus, 25, 104 [NASA ADS] [CrossRef] [Google Scholar]
- DeMeo, F. E., Binzel, R. P., Slivan, S. M., & Bus, S. J. 2009, Icarus, 202, 160 [NASA ADS] [CrossRef] [Google Scholar]
- Duffard, R., & Roig, F. 2009, Planet. Space Sci., 57, 229 [NASA ADS] [CrossRef] [Google Scholar]
- Florczak, M., Lazzaro, D., & Duffard, R. 2002, Icarus, 159, 178 [NASA ADS] [CrossRef] [Google Scholar]
- Fukugita, M., Ichikawa, T., Gunn, J. E., et al. 1996, ApJ, 111, 1748 [Google Scholar]
- Gaffey, M. J., Burbine, T. H., Piatek, J. L., et al. 1993, Icarus, 106, 573 [NASA ADS] [CrossRef] [Google Scholar]
- Gaffey, M. J., Cloutis, E. A., Kelley, M. S., & Reed, K. L. 2002, Asteroids III, 183 [Google Scholar]
- Gradie, J. C., Chapman, C. R., & Tedesco, E. F. 1989, in Asteroids II, ed. R. P. Binzel, T. Gehrels, & M. S. Matthews, 316 [Google Scholar]
- Hapke, B. 2001, J. Geophys. Res., 106, 10039 [NASA ADS] [CrossRef] [Google Scholar]
- Hiroi, T., Kanno, A., Nakamura, R., et al. 2003, in Lunar and Planetary Inst. Technical Report, Lunar and Planetary Institute Science Conference Abstracts, ed. S. Mackwell & E. Stansbery, 34, 1425 [Google Scholar]
- Ivezic, Z., Tabachnik, S., Rafikov, R., et al. 2001, ApJ, 122, 2749 [Google Scholar]
- Jewitt, D. C. 2002, AJ, 123, 1039 [NASA ADS] [CrossRef] [Google Scholar]
- Jones, T. D., Lebofsky, L. A., Lewis, J. S., & Marley, M. S. 1990, Icarus, 88, 172 [NASA ADS] [CrossRef] [Google Scholar]
- Lazzaro, D., Michtchenko, T., Carvano, J. M., et al. 2000, Science, 288, 2033 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Lazzaro, D., Angeli, C. A., Carvano, J. M., et al. 2004, Icarus, 172, 179 [CrossRef] [Google Scholar]
- McCord, T. B., Adams, J. B., & Johnson, T. V. 1970, Science, 168, 1445 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- McFadden, L. A., Gaffey, M. J., & McCord, T. B. 1984, Icarus, 59, 25 [NASA ADS] [CrossRef] [Google Scholar]
- McFadden, L. A., Gaffey, M. J., & McCord, T. B. 1985, Science, 229, 160 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Misra, A., & Bus, S. J. 2008, in BAAS, 40, 40, 508 [Google Scholar]
- Moskovitz, N. A., Jedicke, R., Gaidos, E., et al. 2008, Icarus, 198, 77 [NASA ADS] [CrossRef] [Google Scholar]
- Mothé-Diniz, T., & Carvano, J. M. 2005, A&A, 442, 727 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mothé-Diniz, T., & Nesvorný, D. 2008, A&A, 486, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mothé-Diniz, T., Carvano, J. M. Á., & Lazzaro, D. 2003, Icarus, 162, 10 [NASA ADS] [CrossRef] [Google Scholar]
- Mothé-Diniz, T., Roig, F., & Carvano, J. M. 2005, Icarus, 174, 54 [NASA ADS] [CrossRef] [Google Scholar]
- Mothé-Diniz, T., Carvano, J. M., Bus, S. J., & Burbine, T. H. 2008, Icarus, 195, 277 [NASA ADS] [CrossRef] [Google Scholar]
- Neese, C. 2006, NASA Planetary Data System, EAR-A-5-DDR-TAXONOMY-V5.0, 48 [Google Scholar]
- Nesvorný, D., Roig, F., Gladman, B., et al. 2008, Icarus, 193, 85 [NASA ADS] [CrossRef] [Google Scholar]
- Parker, A., Ivezic, Z., Juric, M., et al. 2008, Icarus, 198, 138 [NASA ADS] [CrossRef] [Google Scholar]
- O'Brien, D. P., Morbidelli, A., & Bottke, W. F. 2007, Icarus, 191, 434 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- Pieters, C. M., & Hiroi, T. 2004, in Lunar and Planetary Inst. Technical Report, Lunar and Planetary Institute Science Conference Abstracts, ed. S. Mackwell, & E. Stansbery, 35, 1720 [Google Scholar]
- Rivkin, A. S., Howell, E. S., Britt, D. T., et al. 1995, Icarus, 117, 90 [NASA ADS] [CrossRef] [Google Scholar]
- Rivkin, A. S., Howell, E. S., Lebofsky, L. A., & Britt, D. T. 2000, Icarus, 145, 351 [NASA ADS] [CrossRef] [Google Scholar]
- Roig, F., & Gil-Hutton, R. 2006, Icarus, 183, 411 [NASA ADS] [CrossRef] [Google Scholar]
- Roig, F., Nesvorný, D., Gil-Hutton, R., & Lazzaro, D. 2008, Icarus, 194, 125 [NASA ADS] [CrossRef] [Google Scholar]
- Sasaki, S., Nakamura, K., Hamabe, Y., et al. 2001, Nature, 410, 555 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Sunshine, J. M., Pieters, C. M., & Pratt, S. F. 1990, J. Geophys. Res., 95, 6955 [Google Scholar]
- Sunshine, J. M., Bus, S. J., Corrigan, C. M., & Burbine, T. H. 2007, Meteoritics and Planetary Science, 42, 155 [NASA ADS] [CrossRef] [Google Scholar]
- Szabó, G. M., Ivezic, Z., Juric, M., & Kiss, L. L. 2004, MNRAS, 348, 987 [NASA ADS] [CrossRef] [Google Scholar]
- Tedesco, E. F., Williams, J. G., Matson, D. L., et al. 1989, AJ, 97, 580 [NASA ADS] [CrossRef] [Google Scholar]
- Tholen, D. J. 1984, Ph.D. Thesis, AA Arizona Univ., Tucson [Google Scholar]
- Vilas, F., & Gaffey, M. J. 1989, Science, 246, 790 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
All Tables
Table 1: Number of objects with unambiguous spectral classification that were observed by the SDSS.
Table 2: Number of asteroids and observations classified.
Table 3: Objects with taxonomic variation between observations.
All Figures
![]() |
Figure 1: Final reflectance color gradient limits for the classes. |
Open with DEXTER | |
In the text |
![]() |
Figure 2: Median log reflectance spectra for the classes, shifted vertically for clarity. |
Open with DEXTER | |
In the text |
![]() |
Figure 3: Distribution of the SDSS taxonomic classes for observations of the asteroids previously classified based on spectroscopy. The solid filling represents the number of observations with uncertainties above the threshold set in Sect. 2.4. |
Open with DEXTER | |
In the text |
![]() |
Figure 4: Distribution of the SDSS taxonomic classes for meteorite spectra from the RELAB database. |
Open with DEXTER | |
In the text |
![]() |
Figure 5: Distribution of the SDSS taxonomic classes for meteorite spectra from the RELAB database. |
Open with DEXTER | |
In the text |
![]() |
Figure 6: Density map of all numbered asteroids. The location of several families are marked. |
Open with DEXTER | |
In the text |
![]() |
Figure 7: Density maps for each photometric class. |
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.