Issue 
A&A
Volume 530, June 2011



Article Number  A134  
Number of page(s)  16  
Section  Planets and planetary systems  
DOI  https://doi.org/10.1051/00046361/201116738  
Published online  25 May 2011 
A study of asteroid polelatitude distribution based on an extended set of shape models derived by the lightcurve inversion method^{⋆}
^{1}
Astronomical Institute, Faculty of Mathematics and Physics, Charles University in Prague, V Holešovičkách 2, 18000 Prague, Czech Republic
email: hanus.home@gmail.com
^{2}
Palmer Divide Observatory, 17995 Bakers Farm Rd., Colorado Springs, CO 80908, USA
^{3}
4438 Organ Mesa Loop, Las Cruces, NM 88011, USA
^{4}
Goat Mountain Astronomical Research Station, 11355 Mount Johnson Court, Rancho Cucamonga, CA 91737, USA
^{5}
Kingsgrove, NSW, Australia
^{6}
Observatoire des Engarouines, 84570 MallemortduComtat, France
^{7}
Via M. Rosa 1, 00012 Colleverde di Guidonia, Rome, Italy
^{8}
Geneva Observatory, 1290 Sauverny, Switzerland
^{9}
Benoziyo Center for Astrophysics, The Weizmann Institute of Science, Rehovot 76100, Israel
^{10}
Astronomical Institute, Academy of Sciences of the Czech Republic, Fričova 1, 25165 Ondřejov, Czech Republic
^{11}
Severni 765, 50003 Hradec Kralove, Czech republic
^{12}
National Astronomical Observatory, Osawa 2211, Mitaka, 1818588 Tokyo, Japan
Received: 17 February 2011
Accepted: 13 April 2011
Context. In the past decade, more than one hundred asteroid models were derived using the lightcurve inversion method. Measured by the number of derived models, lightcurve inversion has become the leading method for asteroid shape determination.
Aims. Tens of thousands of sparseintime lightcurves from astrometric projects are publicly available. We investigate these data and use them in the lightcurve inversion method to derive new asteroid models. By having a greater number of models with known physical properties, we can gain a better insight into the nature of individual objects and into the whole asteroid population.
Methods. We use sparse photometry from selected observatories from the AstDyS database (Asteroids – Dynamic Site), either alone or in combination with dense lightcurves, to determine new asteroid models by the lightcurve inversion method. We investigate various correlations between several asteroid parameters and characteristics such as the rotational state and diameter or family membership. We focus on the distribution of ecliptic latitudes of pole directions. We create a synthetic uniform distribution of latitudes, compute the method bias, and compare the results with the distribution of known models. We also construct a model for the longterm evolution of spins.
Results. We present 80 new asteroid models derived from combined data sets where sparse photometry is taken from the AstDyS database and dense lightcurves are from the Uppsala Asteroid Photometric Catalogue (UAPC) and from several individual observers. For 18 asteroids, we present updated shape solutions based on new photometric data. For another 30 asteroids we present their partial models, i.e., an accurate period value and an estimate of the ecliptic latitude of the pole. The addition of new models increases the total number of models derived by the lightcurve inversion method to ~200. We also present a simple statistical analysis of physical properties of asteroids where we look for possible correlations between various physical parameters with an emphasis on the spin vector. We present the observed and debiased distributions of ecliptic latitudes with respect to different size ranges of asteroids as well as a simple theoretical model of the latitude distribution and then compare its predictions with the observed distributions. From this analysis we find that the latitude distribution of small asteroids (D < 30 km) is clustered towards ecliptic poles and can be explained by the YORP thermal effect while the latitude distribution of larger asteroids (D > 60 km) exhibits an evident excess of prograde rotators, probably of primordial origin.
Key words: minor planets, asteroids: general / thechniques: photometric / methods: numerical / methods: statistical
Tables 3–6 are available in electronic form at http://www.aanda.org
© ESO, 2011
1. Introduction
The lightcurve inversion method (LI) is a powerful tool that allows us to derive basic physical properties of asteroids (the rotational state and the shape) from their diskintegrated photometry (see Kaasalainen & Torppa 2001; Kaasalainen et al. 2001, 2002). This photometry can be denseintime, which typically consists of tens to a few hundreds of individual data points observed during one revolution. This is in contrast to sparseintime, where the typical separation of individual measurements is large compared to the rotation period. For sparse data, we usually have a few measurements per night, such as in the case of astrometric sky surveys. In the following text, we use the terms “dense lightcurves” and “sparse lightcurves”.
To obtain a unique spin and shape solution, we need a set of at least a few tens of dense lightcurves observed during at least three apparitions. Based on simulated photometric data and the survey cadence of the Panoramic Survey Telescope And Rapid Response System (PanSTARRS), Kaasalainen (2004) showed that we can also use only sparse data for the inversion technique. In this case, a unique model can be derived from more than about one hundred calibrated measurements observed during 3–5 years if the photometric accuracy is better than ~ 5% (Ďurech et al. 2005, 2007). Sparse data available so far are not that accurate. Nevertheless, for many asteroids with high lightcurve amplitudes, it is possible to derive their models from current sparse data (see Sect. 2 for more details). We can also combine sparse and dense data to derive models. First results from this approach were presented by Ďurech et al. (2009), where sparse data from the US Naval Observatory in Flagstaff (USNO) were used.
Currently (January 2011), there are 113 models of asteroids derived by the lightcurve inversion method; most of them are stored in the Database of Asteroid Models from Inversion Techniques (DAMIT^{1}, Ďurech et al. 2010). Most of these models were derived from dense lightcurves. Only 24 of them were computed from combined dense and sparse data (Ďurech et al. 2009). The AstDyS database (Asteroids – Dynamic Site^{2}), which contains data from astrometric projects, is another possible source of sparse data. However, most of the data are not accurate enough to be used for inversion alone. On the other hand, even noisy sparse data in combination with a few dense lightcurves can give us, in many cases, a unique solution (Ďurech et al. 2007). The aim of our work was to gather these data, keep only those that were useful, and then combine them with dense lightcurves in the lightcurve inversion method.
Dense data are best used to define the rotational period and constrain the period interval that must be searched during the model computation (see Sect. 3.1 for more details). On the other hand, sparse data usually cover a long time interval, typically over several apparitions, and carry information about brightness variations for different geometries, which constrains the pole directions.
A priori information about rotational periods of asteroids plays an important role in the process of model determination. When an approximate period is known, we search for the solution near this value (details in Sect. 3.1) and thus save considerable computational time. We use the latest update of the Minor Planet Lightcurve Database^{3} published by Warner et al. (2009) to check for previously derived periods. For many asteroids, there are only a few sparse lightcurves from different astrometric observatories available but no dense lightcurves. In these cases, we must scan the whole interval of expected period values (2–30 h). This approach is timeconsuming and there is no guarantee that the correct period will lie in the scanned interval.
The knowledge of rotational states of asteroids is fundamental for understanding the history of the Solar System, specifically the accretion of planets or the collisional processes. For example, it was presumed that due to collisional evolution, the spinvector distribution of main belt asteroids (MBAs) should be nearly isotropic, possibly with a small excess of prograde spins (Davis et al. 1989). Johansen & Lacerda (2010) performed a hydrodynamical simulation of the accretion of pebbles and rocks onto protoplanets and speculated that the trend of prograde rotators among the largest asteroids is primordial.
First statistical analyses of the spinvector distribution were presented by Magnusson (1986, 1990) and Drummond et al. (1988, 1991), later by Pravec et al. (2002) and Skoglöv & Erikson (2002). They all observed a lack of poles close to the ecliptic plane. Kryszczyńska et al. (2007) used more objects in the analysis, finding that the distribution was strongly anisotropic with a moderate excess of prograde spins in the limited size range from 100 to 150 km. Interpretation of this depopulation of poles close to the ecliptic plane is still unclear. Probable candidates are selection effects, the role of inclination, the YORP effect^{4} (Rubincam 2000; Vokrouhlický et al. 2003), or a combination of these. The YORP effect acts only on small bodies with D ≲ 40 km. Asteroids with these sizes have nonMaxwellian spin rate distribution (Pravec & Harris 2000) and is particularly evident for asteroids with D < 14 km (Warner et al. 2009). It is believed that the YORP effect is responsible for this trend since it can either spin up or spin down an irregularlyshaped asteroid on the timescale shorter than the typical time between collisions and also affects the obliquity of spin axes (Rubincam 2000; Bottke et al. 2006).
In the near Earth asteroids (NEAs) population, the latitude distribution of poles is different from that of MBAs (La Spina et al. 2004; Kryszczyńska et al. 2007), i.e., there is a significant excess of retrograde spins probably caused by the transport mechanism of MBAs to Earthcrossing space by gravitational resonances and the Yarkovsky effect^{5} (Morbidelli & Vokrouhlický 2003). There is no statistically significant clustering in the longitude of poles of either MBAs or NEAs (Kryszczyńska et al. 2007).
Comparison of estimated characteristics of residuals for 13 selected observatories: mode, “FWHM” and median.
As the number of asteroid models with known physical properties grows, we can have a better insight into the nature of individual objects and into the asteroid population as a whole. In Sect. 2, we describe available dense and sparse photometric data and the selection of astrometric observatories with quality sparse data. In Sect. 3, we present new asteroid models derived from combined photometric data sets or from sparse data alone, mentioning a few individual objects and define several procedures on how to test the reliability of new models. In Sect. 4, we present a statistical analysis of asteroid physical parameters that we derived using the lightcurve inversion method or adopted from different sources (proper elements from the AstDyS database, diameters from IRAS, ...). We also present results of a numerical simulation that allowed us to estimate the bias in pole directions of the lightcurve inversion method. Using these results, we then corrected the observed pole distributions for this effect. Finally, in order to explain the observed latitude distributions, we present a simple theoretical model of the latitude distribution in Sect. 5.
2. Photometric data
The main source of dense photometric lightcurves is the Uppsala Asteroid Photometric Catalogue (UAPC, Lagerkvist et al. 1987; Piironen et al. 2001), where the lightcurves for about 1 000 asteroids are stored. We also used data from several individual observers (Table 6).
Sparse photometry was first used in combination with dense data for lightcurve inversion by Ďurech et al. (2009). These sparse data were from the USNOFlagstaff station and had a typical photometric uncertainty of ~8–10%. Other sparse photometric measurements are produced by many astrometric surveys, but mostly as a byproduct. In most cases, asteroid magnitudes are given to only one decimal place, i.e., the accuracy is 0.1 mag at best. Whether or not this is sufficient for a unique shape determination for reasonable number of asteroids can be deduced from asteroids lightcurve amplitude distribution. We used lightcurve amplitude data for ~2500 asteroids from the Minor Planet Lightcurve Database (Warner et al. 2009) and found that the mean lightcurve amplitude is ~0.3 mag. For 19% of asteroids, the amplitude is ≥0.5 mag. This means that, in principle, photometry with an accuracy of ~0.1 mag carries sufficient information about rotational states and shapes for a significant number of asteroids.
Our goal was to find out which observatories produce photometry suitable for lightcurve inversion and to use these data for determining new asteroid models. Through to September 2009 (the time of the data download), data for more than 350 000 objects from almost 1500 observatories were archived on the AstDyS server. Some of the observatories contributed with only a few data points, while others contributed tens of thousands of photometric measurements (e.g., large sky surveys such as the Catalina Sky Survey, LONEOS, or Siding Spring Survey).
2.1. Data reduction
The quality of the sparse photometry archived on the AstDyS varies significantly. We investigated the photometry carefully by establishing criteria for its quality. Then, using those criteria, we choose only those data that were useful for inversion.
For each observatory, we extracted photometric data for the first 10 000 numbered asteroids if there were at least 30 data points for a single lightcurve. We then transformed this photometry to the standard format used in lightcurve inversion: we computed geometry of observation (astrocentric ecliptic coordinates of the Sun and the Earth), corrected for lighttime, normalized the brightness to the distance of 1 AU from the Sun and the Earth, and excluded clear outliers.
For further investigation, we selected 13 observatories that fulfilled the condition of having data for more than ~50 asteroids. This resulted in almost 30 000 sparse lightcurves for ~9000 asteroids. In the next step, we estimated mean uncertainties of individual observatories and, based on these uncertainties, we assigned a relative weight to the data from each observatory. In this process, we assumed that the brightness vs. solar phase angle relation can be fitted with a simple relation for each sparse lightcurve: (1)where α is the solar phase angle^{6} and a,b,c and d are free parameters. Then, we constructed a histogram of residuals (rms) for each observatory comparing actual data against the model given by Eq. (1). Four examples are plotted in Fig. 1. The dispersion is caused by observational uncertainties and by the amplitudes of the lightcurves. From these histograms, we estimated the “FWHM”^{7} values and the most frequent residual (the mode); median values of the residual distributions for each location (Table 1) were computed. Observatories with a high median or “FWHM” value (≳ 0.2 mag for both) are not suitable for the lightcurve inversion (e.g. observatory 691 in Fig. 1). Data from only seven observatories, listed in Table 1 with nonzero weights, had sufficient accuracy and so could be used for modeling. Based on the values of medians and “FWHMs”, we estimated a weight for the photometric data from each observatory relative to dense data, which has a unity weight (see Table 1). We assumed that the typical accuracy of dense lightcurves is ~0.02 mag.
Fig. 1 Four histograms of residuals comparing actual data against the model given by Eq. (1) of all sparse lightcurves belonging to the following observatories: 689 USNO, 691 Steward Observatory, 703 Catalina Sky Survey and E12 Siding Spring Survey. Number of bins is ~, where N is the total number of sparse lightcurves used for histogram construction. 
The USNO in Flagstaff (MPC code 689) and the Hipparcos satellite are clearly the best observatories with respect to photometric accuracy. Other observatories are less accurate but they still hold some information about rotational states and shapes of asteroids. For any given asteroid, we have typically 2–4 sparse lightcurves from different observatories covering the last ~10–15 years. Data from USNO were already used by Ďurech et al. (2009). We updated those data along with adding another 20–30% of new data from the years 2008 and 2009 if there was an apparition for the asteroid.
3. Results
3.1. Models of asteroids
Our lightcurve inversion (LI) method is based on the optimization of unknown parameters of the shape (modeled as a convex hull), the rotational state, and the scattering law (see Kaasalainen & Torppa 2001; Kaasalainen et al. 2001). The parameter space have many local minims. Since LI is a gradientbased method that converges to a local minimum near the initial choice of parameter values, it is critical to find the global minimum for the data set and then do the modeling. Finding a global minimum involves a systematic search through all relevant parameter values. Each model corresponding to a particular local minimum is characterized by a single value χ^{2}, which corresponds to the quality of the fit.
A unique solution is defined as follows: (i) the best period has at least 10% lower χ^{2} than all other periods in the scanned interval; (ii) for this period, there is only one pole solution with at least 10% lower χ^{2} than the others (with a possible ambiguity for λ ± 180°); and (iii) this solution fulfills our additional tests (see Sect. 3.3).
The most timeconsuming part of the lightcurve inversion method is scanning through all periods within a chosen interval, which we determined by using the periods and reliability codes given in the Minor Planet Lightcurve Database^{8}. For each asteroid, we scanned an interval centered at the reported period value P with a range of ± 1%, ± 5% and ± 20% of P for reliability codes 4, 3 and 2, respectively. Half and double periods were tested later on.
We combined relative lightcurves from the UAPC and from individual observers together with sparse data obtained from the AstDyS site to create a data set for each asteroid. This gave us data sets for ~2300 asteroids (in ~900 cases there were only sparse data available) to which we then applied the lightcurve inversion method and then we ran the additional tests described in Sect. 3.3. We derived 80 new unique models, 16 of which are based only on sparse data. Basic characteristics of these models are listed in Table 3. We estimated the uncertainty in the pole direction as ± 10–20° based on previous results with limited data sets. As might be expected, the uncertainty seems dependent on the number of dense and sparse photometric data. The longitude uncertainty increases for higher latitudes because meridians on a (λ, β)sphere are more dense with increasing latitude. These uncertainties are discussed in more detail in Sect. 4.3. The uncertainty of the rotational period depends on the time interval covered by observational data and is of the order of the last decimal place of period values P in Tables 3–5.
In some cases, we were able to determine a unique rotational period, but we had multiple pole solutions with similar ecliptic latitudes β. These models give us accurate period values and rough estimates of ecliptic latitudes β, which are also important parameters. In Table 4, we present results for 30 partial models, where β is the mean value for all different models if the dispersion is smaller than 50°. We defined a parameter Δ = β_{max} − β_{min}/2 as being the estimated uncertainty of β, where β_{max} and β_{min} are the extremal values within all β.
All new unique shape models are now included in DAMIT.
3.2. Comments to selected models
In DAMIT there are several solutions designated as “preliminary”. These models do not have a wellconstrained pole solution or are based on combined data sets. For 18 of those asteroids we derived updated model solutions based on additional photometric data (see Table 5). The difference between the old and new model for asteroid (1223) Neckar was significant. The new model has a slightly different period, but the pole directions and shapes are nearly similar to the old model. The current data suggest a period of P = 7.82401 h (previous value was P = 7.82123 h).
The asteroid (4483) Petofi was recently observed by Brian Warner. We derived a shape solution from three poor dense lightcurves and one sparse lightcurve from Catalina Sky Survey. Warner used these four lightcurves in combination with his new observations and also derived the shape model of Petofi (Warner 2011a). His period of P = 4.3330 h and pole direction (90°, 35°) are close to our solution of P = 4.33299 h and (107°, 40°).
The asteroid (832) Karin was also studied by Slivan & Molnar (2010); their solution with P = 18.352 h and pole (51° or 228°, 41°) confirms our results, see Table 3.
In past decades, occultations of stars by several asteroids were observed. These events give us additional information about the shape and can help resolve which mirror solutions of a model is the correct one. According to the recent work of Ďurech et al. (2011), asteroid occultation measurements prefer pole solutions of (122°, − 44°) for (10) Hygiea, (347°, 47°) for (152) Atala, (28°, − 72°) for (302) Clarissa, (223°, 67°) for (471) Papagena, and (296°, 41°) for (925) Alphonsina. Spin solutions preferred by asteroid occultations appear in bold font in Tables 3 and 5.
3.3. Models and method testing
We constructed five additional tests to be sure that the new models are reliable. We performed the first two tests for all models. For models derived only from sparse data, which are presented for the first time, we performed three additional tests:
Inertia tensor of the shape model.
The lightcurve inversion method we use assumes that asteroids are in a relaxed rotational state, which means that derived models should rotate around the axis with a maximum moment of inertia. For each derived shape, we computed principal moments using equations presented by Dobrovolskis (1996) and checked if the rotation axis was close to the principal axis of the maximum momentum of inertia. We rejected those models for which the angle between the spin axis and the axis with a maximum momentum of inertia was larger than ± 30°. However, this criterion is too strict for elongated models with similar sizes along the rotational axis and the axis that is both perpendicular to the rotational axis and is the minimal size of the model. In this case, the principal moments for these two axes are similar. Under these circumstances, the angle between the spin axis and the axis with a maximum momentum of inertia can be large even for realistic shapes and so we allowed these models to pass this test.
Half and doubleperiod models.
In cases where we have only a few dense lightcurves for a given asteroid, it is easy to confuse the correct rotational period with its half or double value. When the a priori period value in the Minor Planet Lightcurve Database is uncertain, which corresponds to a low reliability code (see Sect. 3.1), it is reasonable to check if the half and double period value give a better fit. If the period was in doubt, we searched for a solution also around 2P and P/2; if the χ^{2} was lower than 1.1χ^{2} of the solution with period P, we rejected the model as unreliable.
Reduction of the number of sparse photometric data.
In this test, we used only the sparse data sets for modeling. For 63 asteroids, this led to unique shape and spin state solutions (after performing the tests described above). For each solution, we randomly reduced the original amount of observed sparse data points to 90% and used these new limited data sets again in the lightcurve inversion. Our expectation was that we would not get a unique solution when using less data. This was true for five asteroids. These models, when using the full sparse data sets, are not necessarily wrong, but the amount of available data is probably just at the level when a unique solution can be derived. The important point of this test is that, for a given asteroid, we did not find two different but formally correct solutions when using the full versus reduced data sets. In Table 3, we present 16 models successfully derived only from sparse data that passed this test.
Models from sparse data vs. DAMIT.
Here, we used previously derived models based only on relative photometry which are stored in DAMIT. As can seen in the previous test, sparse data are sometimes sufficient to produce a unique model of an asteroid. In 16 cases, we were able to derive a model for an asteroid which was already included in DAMIT and thus a model based on entirely different photometric data sets is available. These two independent models can be then compared and should be similar. We obtained similar resulting models for all 16 asteroids.
Models of “mock” objects.
For each asteroid shape model derived only from sparse data, we created a set of ten “mock” objects of roughly the same appearance and spin state (see an example of such shape in Fig. 2). For these synthetic objects, we computed their photometric data using same epochs and geometries and with similar random noise level. These synthetic photometric data sets were then used in the lightcurve inversion method. A check that the original model using actual data is reliable is to be able to derive most of the models of the “mock” objects. The dispersions between the periods and pole directions of the “mock” objects represent the typical uncertainties of these parameters. For all studied asteroids, we were able to derive unique models for most of their “mock” objects. In all cases when we did not get a unique solution for the “mock” object, the best fit corresponded to the correct solution although other solutions could not be ruled out. The typical uncertainty in pole direction was ± 10° and, for the period, ~ 0.1 times the difference between the two local minimums as determined by the period value and the time span of the data.
Fig. 2 Asteroid (810) Atossa: shape model (left panel) and an example of a “mock” shape model (right panel). 
4. An analysis of periods, poles and sizes of asteroid models
Previous studies have been looking on spin state results using different techniques, e.g. amplitudeepoch, lightcurve inversion or radar methods. If there were multiple solutions for a given asteroid, the most probable one or simply the weighted mean was taken. However, this could cause systematic deviations. In our study, we used only results based on the lightcurve inversion method – i.e., unique and partial models presented in this work and models from the DAMIT. Our sample consists of 221 asteroid models: 80 new models, 18 updates for models from DAMIT, 30 new partial models, 84 models from DAMIT and 9 new models presented by Ďurech et al. (2011). Our sample consists of models for 206 MBAs, 10 NEAs, 3 Hungaria, 1 Trojan, and 1 Hilda and so a statistical study is only possible for the MBAs. In many cases there is an ambiguity in the pole direction since there are two, undistinguishable mirror solutions. For our statistical analysis we randomly chose only one.
In Fig. 3, we show (among other things to be discussed later) the relation between the proper semimajor axis and the proper eccentricity for asteroids in our sample and for all main belt asteroids. It is obvious that the positions of studied asteroids strongly correlate with the MBAs population and so derived models are not significantly biased with respect to orbits, e.g. they do not lie in the inner main belt. Several asteroids in Fig. 3 with semimajor axis a > 3.3 AU belong to the Cybele group, e.g. (121) Hermione.
Fig. 3 Relation between the proper semimajor axis and the proper eccentricity for asteroids in our sample and for the first 100 000 numbered asteroids for comparison. Main resonances are shown by dotted lines. Prograde rotators are plotted with blue circles and retrograde rotators with red squares. The horizontal lines represent for each asteroid its estimated past drift (i.e. where the asteroid came from) during the collisional lifetime computed with Eq. (7). Proper elements are from the AstDyS database. 
4.1. Pole distribution analysis
In the following study of spin axis directions, we did not use the Koronis family members because their spin states are correlated, i.e., their spin vectors are clustered towards two values of the obliquity (Slivan 2002). In Fig. 4a, we show the ecliptic latitude distribution of our MBA sample. As in all similar plots, the width of the latitude and longitude bins corresponds to equal surfaces on the (λ, β)sphere (bins are equidistant in sinβ for latitudes and in λ for longitudes). We confirmed the expectation that there is a lack of asteroids with latitudes close to the ecliptic plane. The latitude distribution is clearly not symmetric: about half of the retrograde rotators have latitudes in the bin (–53°, –90°). On the other hand, less than a third of asteroids with prograde spins are in the corresponding bin (53°, 90°). Moreover, the remaining prograde bins are more populated than the corresponding retrograde ones. From a detailed look at the plot we can see that there are up to 10% more prograde rotators among the MBAs.
In Fig. 5, we show the dependence of the ecliptic latitude β of the pole direction on the diameter D (most of the diameters used are based on the IRAS data (Tedesco et al. 2002) or occultations profiles and have an uncertainty ± 10%). Even for diameters D ≲ 50 km, the clustering of the latitudes towards higher absolute values, and conversely, the depletion close to the ecliptic plane is obvious and markedly so for D ≲ 30 km. Asteroids with larger diameters have a more isotropic distribution of latitudes and only a moderate excess of prograde. In Fig. 6, we plotted the latitude and longitude distributions of asteroids with respect to their diameters. Based on Fig. 5, we resolved three different size groups: 0–30 km, 30–60 km and > 60 km.
The latitude distribution for asteroids with D > 60 km (Fig. 6a) is close to the uniform distribution for latitudes lower than 11° and for larger latitudes it exhibits an evident excess of prograde rotators. This is in agreement with theoretical arguments presented by Davis et al. (1989) and recently by Johansen & Lacerda (2010). On the other hand, the latitude distribution for asteroids with D < 30 km (Fig. 6e) exhibits a strong depopulation of pole vectors close to the ecliptic plane (i.e. small absolute values of latitudes β). The few asteroids with small latitudes have diameters D > 25 km. The latitude distribution for asteroids with intermediate diameters of 30–60 km (Fig. 6c) is also somewhat clustered towards higher latitudes but the bins for small latitudes are more populated. Therefor, it is probably a transition region between the two distinct distributions.
Fig. 4 a) Shows the ecliptic latitude distribution of all MBAs in our sample except the Koronis family members. The width of bins in latitude β corresponds to similar surfaces on the (λ, β)sphere, so the bins are equidistant in sinβ. In b), the longitude distribution of all MBAs is plotted, again except the Koronis family members. The longitude bins are equidistant in λ. In the top right corners, there are the probability values of the χ^{2}tests (that the observed distributions are randomly drown from a uniform distribution, see Table 2). 
Fig. 5 Dependence of the ecliptic latitude β (plotted as sinβ) of the pole direction on the models diameter D. The gray box shows the gap of small latitudes for asteroids with D < 30 km. Members of the Koronis cluster are plotted with crosses. 
Fig. 6 Histograms showing the observed latitude and longitude distributions of MBAs (except Koronis cluster members) for different size ranges. a) Shows the latitude distribution for asteroids with diameters larger than 60 km, c) for asteroids with diameters in the range of 30–60 km, e) for asteroids with diameters smaller than 30 km, and similarly b), d) and f) for longitudes. The width of bins in the latitude β and longitude lambda λ corresponds to similar surfaces on the (λ, β)sphere, so the bins are equidistant in sinβ and λ. In the top right corners, there are the probability values of the χ^{2}tests (that the observed distributions are randomly drown from a uniform distribution, see Table 2). 
Fig. 7 Dependence of the ecliptic latitude β (plotted as sinβ) of the pole direction on the models spin rate, a) for asteroids with diameter D < 30 km and b) for asteroids with diameter D > 60 km. 
It is evident that the depopulation concerns mainly objects with diameters D ≲ 30 km (the distribution for the intermediate size sample shows that the limit is probably ~50 km). This size roughly corresponds to the value, when the YORP effect starts to act and hence it is a natural candidate for a physical explanation. It is known from previous studies (Pravec & Harris 2000; Rubincam 2000) that the YORP effect is significantly altering the periods and also spin vectors of these objects on a timescale shorter than the typical collisional lifetime of these objects (timescales are discussed in more details in Sect. 5). In Fig. 7, we show the relation between the spin rate and the latitude for the small (D < 30 km) and large (D > 60 km) groups of asteroids. In concert with YORP theory, the spin up and spin down and the simultaneous evolution of the latitudes towards higher absolute values are evident in the small asteroid sample.
Note that observed latitude and longitude distributions can be biased by the convex inversion method and, therefore, this bias should be taken into consideration. In general, models for asteroids with higher amplitudes are more often successfully derived than those asteroids with lower amplitudes. This is because, while the accuracy of the sparse photometry in both cases is roughly the same, the signaltonoise ratio is significantly better for higher amplitude lightcurves. If we assume two bodies with the same shapes and orbits and different ecliptic latitudes of the poles, the body with higher absolute value of ecliptic latitude usually has a higher amplitude. This effect is numerically investigated in Sect. 4.3.
Test of the hypothesis that the observed pole distributions and debiased latitude distribution are drown from uniform distributions (a χ^{2}test).
Skoglöv & Erikson (2002) discussed the role of orbital inclination in the observed distribution of latitudes. No indications of this effect were found in the asteroids sample of Kryszczyńska et al. (2007). We also did not find any indications of such correlation in our sample and so we conclude that orbital inclination does not affect the observed distribution of latitudes.
We are not aware of any other physical effects in the main belt that could explain the nonuniform observed latitude distribution of small asteroids (D < 30 km). Collisions are believed to produce uniform spin distributions and close encounters with planets are common only among NEAs.
There are many additional selection effects that influence the properties of derived models, e.g. the role of amplitude, orbit, time, accuracy and geometry of observations, among others. The significance of this bias is unknown and cannot be easily determined. The main problem here is that for almost every asteroid, the photometric data are from different observers with a different number of measurements, quality, and purpose. The significance of this effect can be determined only from a comparison of models derived from real and synthetic data of known properties. This may be possible in a few years when the photometry from the PanSTARRS is available, but not now. In the meantime, the role of the selection effects seems to be small and does not significantly affect, for example, the latitude distribution.
The longitude distributions of the MBAs are plotted in Figs. 4b and 6b, d, f. They are, contrary to the latitude distributions, without any statistically significant features and have very close to uniform distributions. The only exception are the asteroids with D < 30 km, but the excess appears to be just a random coincidence than a result of some physical process.
In all cases, we tested a hypothesis that the observed distribution of latitudes or longitudes is uniform (using a χ^{2}test^{9}). The computed chisquares and corresponding probabilities are listed in Table 2. Higher χ^{2}values and lower probabilities mean that the supposed hypothesis “the observed distribution is uniform” does not fit the observed data. If we assume a probability of 5% or lower as statistically significant, we can say that the latitude distributions for the whole sample and for asteroids with D < 30 km and 30 < D < 60 km do not agree with a uniform distribution. On the other hand, all longitude distributions are consistent with uniform distributions. Latitude distribution for the MBAs with diameters D > 60 km also disagrees with the uniform distribution; this is because of the excess of prograde rotators.
The overall view on the model positions within the main belt of asteroids, together with their estimated total drifts and the information about whether they are prograde or retrograde rotators (Fig. 3), shows behavior consistent with the Yarkovsky/YORP theory: there is an asymmetry of prograde and retrograde rotators near the main resonances and prograde asteroids drift outwards from the Sun and can reach the resonance only from the left. On the other hand, retrograde rotators drift the opposite direction and depopulate the zone left of the resonance because the resonance prevents the entry of new retrograde asteroids. This creates an excess of prograde rotators. The same mechanism works also in the zone right of the resonance, except, in this case, only an excess of retrograde asteroids is now created. This effect is obvious in the neighborhood of the ν_{6} and 3:1 resonances. The total drift from an asteroid’s original location during the collisional lifetime computed with Eq. (7) are inversely proportional to the size of the asteroids. Larger asteroids (D ≳ 50 km) do not drift significantly while smaller asteroids frequently drift ± 0.05 AU. Note that smallsized asteroids are found mainly in the inner or middle part of the main belt (due to selection effect; they have high albedos and/or are closer to the Earth). Asteroids drifting through the resonances during their collisional lifetime are interesting from the point of the dynamical evolution. These asteroids were either recently collisionallyaffected or their shape models are wrong. There seems to be three such models, two near the 3:1 resonance and one near the resonance 5:2.
4.2. The Koronis family members
The analysis of rotational state solutions for ten members of the Koronis asteroid family revealed a clustered distribution of their spin vectors (Slivan 2002; Slivan et al. 2003). This spin distribution was later explained by Vokrouhlický et al. (2003) as the result of the thermal torques and spinorbital resonances that modify the spin states over time of 2–3 Gyr. The modeling suggested an existence of two groups of asteroids: (a) lowobliquity retrograde objects with rotational periods P < 5 h or P > 13 h, and (b) prograde rotators with periods 4 < P < 7 h that became trapped in a spinorbit resonance with secular frequency s_{6} and thus have similar spin obliquities (42−51°) and also similar ecliptic longitudes in the range of (24−73°) and (204−259°). All ten members of the Koronis family studied by Slivan (2002) and Slivan et al. (2003) had the expected properties: periods for prograde rotators were shifted only to higher values of 7–10 h. Slivan et al. (2009) published spin state solutions for another four members of Koronis family. Only the solution for (253) Dresda was not in agreement with the theoretical expectation.
Here, we present three new models of asteroids belonging to the Koronis family: (832) Karin, (1482) Sebastiana, and (1635) Bohrmann, along with two partial models for (1350) Rosselina and (1389) Onnie. Only the spin state solutions for Bohrmann and Onnie fit the theoretical expectations. Rotational parameters for Karin (P = 18.3512 h, λ = 242°, β = 46°) are outside both groups. Asteroids Sebastiana and Rosselina are lowobliquity retrograde rotators, but their periods (10.49 h for Sebastiana and 8.14 h for Rosselina) are in the middle of the “wrong” range of P = 5−13 h. Karin is the namesake and largest member of a small and young (~5.8 My, Nesvorný & Bottke 2004) collisional family that is confined within the larger Koronis family. The spin state of Karin was thus likely affected during this catastrophic event and changed to a random state that disagrees with the clustered distribution.
We are not able to give a satisfactory explanation for the peculiar spin state solutions for Sebastiana and Rosselina. Nevertheless, we are aware of two possible scenarios: (i) the initial rotational state and shape did not allow being captured in the resonance or (ii) the objects were randomly reoriented by noncatastrophic collisions. The timescales of such collisions (given by Eq. (7)) are for Sebastiana τ_{reor} ~ 7.5 Gyr and for Rosselina τ_{reor} ~ 14.7 Gyr. This leads to the probability of a collision during the Koronis cluster lifetime (estimated to ~2.5 Gyr, Bottke et al. 2001) ~1/3 for Sebastiana and ~1/6 for Rosselina, respectively, which means that random collisional reorientation of the spin axis is likely for at least a few of 19 asteroids in the Koronis cluster with known spin state solutions (most of them have τ_{reor} ≲ 20 Gyr).
4.3. Biases of the LI method
We developed a numerical algorithm to estimate the selection effect of the lightcurve inversion method and used this approach to debias the observed distribution of asteroid’s pole directions. The algorithm was as follows:

1.
for a model with a known shape, we randomly generated a newpole direction (while the overall distribution of poles wasisotropic);

2.
for each shape with a new rotational state but with the period unchanged, we computed synthetic lightcurves for the same epochs as observed ones;

3.
to each data point i, we added the corresponding noise δ_{i} given by: (2)where is ith brightness observed and is ith brightness computed, both for the original model. This gave us synthetic lightcurve data equivalent to the original observed data, but for a new pole direction;

4.
finally, we performed a lightcurve inversion the same way as with the actual data and tried to derive a model;

5.
we repeated steps 1–4 for 50 random poles for each asteroid model.
In this simulation, we used 80 models derived from combined dense and sparse data sets and 89 models from the DAMIT.
For each successfully derived model we have the generated pole direction (λ_{gen}, β_{gen}) and period P_{gen}, and also the computed pole and period: λ_{comp}, β_{comp}, P_{comp}. We computed the angular differences between the generated and derived ecliptic latitude and longitude of the pole direction: β_{comp} − β_{gen} and (λ_{comp} − λ_{gen})cosβ_{gen} (the cosβ_{gen} factor is used for the correction of the different distances of meridians near the equator and poles). In Fig. 8, we show the histograms of these differences in a) ecliptic latitudes and b) ecliptic longitudes. We assumed that the histograms can be described by a normal distribution and we computed the mean and the standard deviation (μ; σ). We found values of (–0.2; 10.2) for latitudes and (–0.2; 5.2) for longitudes. The standard deviation σ is directly related to the typical uncertainty that we can expect in pole determination by the lightcurve inversion method, which is ~5°/cosβ in λ and ~10° in β.
In Fig. 9a, we constructed a histogram of the latitude distribution for all successful models. The bins in β were again equally spaced in sinβ. The latitude distribution of all generated models was not exactly uniform, the amount of latitudes in bins slightly differed. To remove this effect, we divided the latitude distribution of successfully derived models by the latitude distribution of all generated models normalized to unity. This correction was also applied to latitude distributions in Figs. 10a, c, e. It is obvious that the LI method is more efficient for asteroids with higher β. The amount of successfully derived models with β ~ 0° is about 30% lower than with β > 53°.
In Figs. 10a, c, e, we constructed the histograms of latitude distributions for successfully derived models and distinguished three size ranges. All three plots look very similar, except that with decreasing size, the ratio between models with β ~ 0° and β > 53° goes down. This ratio is ~75% for D > 60 km, ~ 65% for 30 < D < 60 km and ~60% for D < 30 km.
Fig. 8 Histograms of the angular differences between the generated and derived ecliptic latitudes and longitudes of the pole directions for all successfully derived synthetic models: a) β_{comp} − β_{gen}, where β_{gen} and β_{comp} are generated and computed ecliptic latitudes and similarly b) (λ_{comp} − λ_{gen})cosβ_{gen} for ecliptic longitudes. 
Fig. 9 a) Shows the simulated ecliptic latitude distribution of all successfully derived models. In b), the debiased observed latitude distribution is plotted. The bins are equidistant in sinβ. In the top right corner, there is the probability value of the χ^{2}test (that the observed distribution is drown from a uniform one, see Table 2). 
Fig. 10 Histograms showing the simulated and corrected (debiased) observed latitude distributions of MBAs for different size groups, excluding the Koronis cluster members. a) Shows the simulated latitude distribution for asteroids with the diameters larger than 60 km, c) for asteroids with the size in range of 30–60 km, e) for asteroids with diameters smaller than 30 km. Similarly b), d) and f) show the observed latitude distributions corrected by the bias of the LI method. The bins are equidistant in sinβ. In the top right corners, there are the probability values of the χ^{2}tests (that the observed distributions are drown from a uniform distribution, see Table 2). 
The histograms in Figs. 9a and 10a, c, e define the bias in latitude of the LI method and can be used for debiasing the observed latitude distributions presented in Figs. 4a and 6a, c, e. The debiased histograms of latitudes are plotted in Figs. 9b and 10b, d, f. The histograms changed only slightly and the conclusions from Sect. 4.1 are still valid, i.e., the latitude distribution differs significantly from a uniform distribution, and especially so for D < 30 km. The distribution of latitudes for asteroids with D > 60 has an evident excess of prograde rotators while the distribution for a subsample with 30 km < D < 60 km shows an enrichment of asteroids with large latitudes ( β > 53°). Other bins have similar populations. We also performed a χ^{2}–test in the same way as for the observed distributions (see Table 2, Cols. 6 and 7).
We did not find any significant correlation between the ecliptic longitude and the efficiency of the model determination.
5. A theoretical model of the latitude distribution
Fig. 11 Left panel: the distribution of the ecliptic latitudes β for the observed asteroids (thick line) and 100 synthetic samples (generated with a different random seed) at time t = 0 (thin lines). Right panel: steadystate synthetic latitude distributions at time t = 4 Gyr (thin lines), evolved by the YORP effect, collisions, mass shedding and spinorbit resonances. We also applied an observational bias to it. The steady state was reached already within ≃ 1 Gyr. 
In order to understand observations of mainbelt asteroids, namely the debiased distribution of their ecliptic latitudes β (Fig. 9b), we constructed a simple model for spin evolution that contains the following processes: (i) the YORP effect, i.e., torques arising from the emission of thermal radiation, (ii) random reorientations induced by noncatastrophic collisions, and (iii) mass shedding after a critical rotational frequency is reached.
On the other hand, we did not include gravitational torques of the Sun, spinorbital resonances, damping (dissipation of rotational energy), or tumbling. Even though individual asteroids may be substantially affected by these processes, our model is for a large statistical sample of asteroids and the effect on the overall latitude distribution is assumed to be only minor. For example, gravitational torques and spinorbital resonances usually cause large oscillations of β for prograderotating asteroids, but they remain bound to a certain interval (Vokrouhlický et al. 2006). Moreover, we tried to account for these (rather random) oscillations in our model as well (see below).
Our sample of 220 asteroids was the same as the observed sample discussed in Sect. 4. This means that the orbits and sizes correspond to real asteroids. The model for spin evolution was similar to that of Brož et al. (2011), where it was used for studies of the longterm evolution of asteroid families. We assumed the following relations for the rate of the angular velocity ω and the obliquity ϵ due to the YORP effect where f and gfunctions were given by Čapek & Vokrouhlický (2004) for a set of 200 shapes with mean radius R_{0} = 1 km, bulk density ρ_{0} = 2500 kg/m^{3}, located on a circular orbit with semimajor axis a_{0} = 2.5 AU. We assigned one of the artificial shapes (denoted by the index i) randomly to each individual asteroid^{10}. We had only to scale the f and gfunctions by a factor (5)where a, R, ρ_{bulk} are semimajor axis, radius, and density of the simulated body, respectively, and c_{YORP} is a free scaling parameter, which can account for an additional uncertainty of the YORP model. Because the values of f’s and g’s were computed for only a limited set of obliquities (with a step Δϵ = 30°) we used interpolation by Hermite polynomials (Hill 1982) of the data in Čapek & Vokrouhlický (2004) to obtain a smooth analytical functions for f_{i}(ϵ) and g_{i}(ϵ).
When the angular velocity approached a critical value (i.e., the gravity was equal to the centrifugal force) (6)we assumed a mass shedding event. We kept the orientation of the spin axis and the sense of rotation but reset the orbital period P = 2π/ω to a random value from the interval (P_{1},P_{2}) = (2.5,9) h. We also altered the assigned shape since any change of shape can produce a different YORP effect. We did not change the mass, however.
The differential Eqs. (3) and (4) were integrated numerically by a simple Euler integrator. The usual time step was Δt = 1000 yr. The time scale of the spin axis evolution for small bodies (D ≃ 10 km) is τ_{YORP} ≃ 500 Myr. After ≃ 3 times τ_{YORP} most of these bodies have spin axes perpendicular to the ecliptic.
We also included a MonteCarlo model for spin axis reorientations caused by collisions^{11}. We used an estimate of the time scale by Farinella et al. (1998)(7)where B = 84.5 kyr, β_{1} = 5/6, β_{2} = 4/3, D_{0} = 2 m and ω_{0} corresponds to period P = 5 h. These values are characteristic for the main belt. After a collision, we reset the spin axis periods to random values, using the interval h for the period. Since the time scale is τ_{reor} ≃ 3 Gyr for the smallest (D ≃ 5 km) bodies, reorientations are only of minor importance. However, note that the probability of the reorientation is enhanced when the YORP effect drives the angular velocity ω close to zero.
There were several free parameters in our model: the c_{YORP} parameter, thermal conductivity K, bulk density ρ_{bulk}, initial distribution of β and initial distribution of ω.
Our aim was to start with a simple β and ωdistribution, wait until a steady state was reached, and then compare the resulting synthetic to observed latitude distributions. We applied an observational bias derived in Sect. 4.3 to the synthetic distribution.
We partly accounted for spinorbital resonances acting on prograde asteroids by adding a sinusoidal oscillations to β with a random phase and an amplitude ≃ 40°, which are typically induced by resonances. This procedure naturally decreased the rightmost bin (sinβ = (0.8,1)) of the synthetic distribution and increased the next bin (sinβ = (0.6,0.8)).
We started with reasonable parameters of c_{YORP} = 0.33, K = 10^{2} W/K/m, ρ_{bulk} = 2500 kg/m^{3}, a Maxwellian distribution of ω, a uniform distribution of sinβ (i.e., an isotropic distribution of spin axes). We ran 100 such simulations with different random seeds. A steady state was reached within ≃ 1 Gyr. The resulting latitude distributions are shown in Figs. 11 and 12.
From these it can be seen that: (i) the observed distribution of β for small asteroids seems compatible with our model; the YORP effect is capable of creating such an uneven distribution and (ii) there is a discrepancy for large asteroids (especially in bins sinβ ∈ (− 1, − 0.8) and (0.2,0.4)), which can be explained as a preference for prograde rotators in the primordial population (see Davis et al. 1989; Johansen & Lacerda 2010). The results regarding the spin rates agree with Pravec et al. (2008), so we do not repeat the discussion here.
We also tested the sensitivity of our results with respect to the free parameters. The thermal conductivity did not seem important (we tested K = 10^{3} W/K/m). A simulation with c_{YORP} = 0.66, ρ = 1300 kg/m^{3}, and a uniform distribution of orbital periods P ∝ 1/ω produced almost the same resulting latitude distribution. Nevertheless, a value of c_{YORP} = 1.00 seems too high because the extreme bins of the βdistribution were overpopulated^{12}. There is only a weak dependence of our results on the period ranges that we used for resetting the orbital period after a massshedding event (P_{1},P_{2}) and collision (). As would be expected, the values P_{1},P_{2}, significantly affect the period distribution. The relatively weak dependence on the free parameters likely stems from the fact that we presume the steadystate. Even though the free parameters change, e.g., the strength of the YORP effect and the evolution of spins is slower and/or faster, after reaching a steady state, the basic characteristics of the latitude distribution remain similar. The observed βdistribution of small asteroids (D < 30 km) cannot be explained by our simulation without accounting the YORP effect.
Fig. 12 Steadystate synthetic latitude distributions for three different size ranges: D > 60 km (top), 30 km < D < 60 km (middle) and D < 30 km (bottom). The synthetic distributions are plotted by thin gray lines while the observed distribution by a thick black line. 
6. Conclusions
The results of this paper can be summarized as follows.
We used combined dense and sparse data to derive new asteroid shape models. We systematically gathered and processed all available sparse photometry from astrometric surveys and employed valuable data from seven observatories (see Table 1) in lightcurve inversion.
We derived 80 new unique models of asteroids, from which 16 are based only on sparse data. We also present 30 partial models with accurate rotational periods and estimated ecliptic latitudes of the pole directions and 18 updated solutions based on new data for asteroids already included in DAMIT.
In the future, quality sparse data sets will be produced by allsky surveys such as PanSTARRS, the Large Synoptic Survey Telescope (LSST), and the Gaia satellite. When these data are available, we will be able to apply the same methods in order to derive many more new asteroid models. These surveys will have one advantage over dense data: their selection effects (e.g., with respect to the orbit) will be known. This will allow us to make a more accurate analysis of the asteroid population.
As expected, the observed ecliptic longitude distribution of asteroid spin vector is independent of diameter and is compatible with a uniform distribution. Unlike the latitude distribution, the observed ecliptic longitude distribution is not significantly biased by the LI method. However, the effect of the LI bias is only minor and the global features of the observed latitude distribution do not change. The observed (and debiased) latitude distribution for asteroids with diameters D > 60 km shows an excess of prograde rotators in the latitude interval (11°, 90°). This excess is probably primordial. On the other hand, the latitude distributions for the entire sample and in particular for asteroids with D < 30 km, is strongly anisotropic.
The dynamical evolution of asteroid spins seems to be dominated by the YORP effect and also by collisions and mass shedding for asteroids with diameters D ≲ 30 km. We calculate that YORP (with a small contribution for the LI method’s bias) is capable of producing the observed depopulation of spin vectors for small asteroids.
We are not yet able to study small asteroids in individual families (small bodies at the outskirts of a family should have extreme spins); this is an aim of future work.
http://cfawww.harvard.edu/iau/lists/LightcurveDat.html, there is also an explanation and more details about the reliability codes too.
We did not use the convexhull shapes derived in this work for two reasons: (i) the two samples of shapes are believed to be statistically equivalent and it is thus not necessary to compute the YORP torques again; (ii) the YORP effect seems sensitive to smallscale surface structure (Scheeres & Mirrahimi 2007) which cannot be caught by our shape model. Nevertheless, the YORP torque remains of the same order, so the random assignment of shapes seems reasonable.
As an alternative hypothesis, we assumed the spin axis evolution without a YORP effect (c_{YORP} = 0). In this case, the initial βdistributions (Fig. 11, left panel, thin lines) do not change significantly in time.
Acknowledgments
The work of J.H. has been supported by the grant GA UK 134710 of the Grant agency of the Charles University, by the project SVV 261301 of the Charles University in Prague and by the grant GACR 205/08/H005 of the Czech grant agency. The work of J.H. and J.Ď. has been supported by the grant GACR 209/10/0537 of the Czech grant agency, the work of J.D. and M.B. by the Research Program MSM0021620860 of the Czech Ministry of Education and the work of M.B. has been also supported by the Grant Agency of the Czech Republic (grant 205/08/P196). The calculations were performed on the computational cluster Tiger at the Astronomical Institute of Charles University in Prague (http://sirrah.troja.mff.cuni.cz/tiger). We thank Brian D. Warner for significantly improving the linguistic quality of this paper and David Čapek for sending us the YORP effect data in an electronic form.
References
 Bottke, W. F., Vokrouhlický, D., Brož, M., Nesvorný, D., & Morbidelli, A. 2001, Science, 294, 1693 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Bottke, J. W. F., Vokrouhlický, D., Rubincam, D. P., & Nesvorný, D. 2006, Ann. Rev. Earth Planet. Sci., 34, 157 [Google Scholar]
 Brož, M., Vokrouhlický, D., Morbidelli, A., Nesvorný, D., & Bottke, W. F. 2011, MNRAS, in press [Google Scholar]
 Čapek, D., & Vokrouhlický, D. 2004, Icarus, 172, 526 [NASA ADS] [CrossRef] [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, 805 [Google Scholar]
 Dobrovolskis, A. R. 1996, Icarus, 124, 698 [NASA ADS] [CrossRef] [Google Scholar]
 Drummond, J. D., Weidenschilling, S. J., Chapman, C. R., & Davis, D. R. 1988, Icarus, 76, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Drummond, J. D., Weidenschilling, S. J., Chapman, C. R., & Davis, D. R. 1991, Icarus, 89, 44 [NASA ADS] [CrossRef] [Google Scholar]
 Ďurech, J., Grav, T., Jedicke, R., Denneau, L., & Kaasalainen, M. 2005, Earth Moon and Planets, 97, 179 [NASA ADS] [CrossRef] [Google Scholar]
 Ďurech, J., Scheirich, P., Kaasalainen, M., et al. 2007, IAU Symp. 236, ed. G. B. Valsecchi, D. Vokrouhlický, & A. Milani, 191 [Google Scholar]
 Ďurech, J., Kaasalainen, M., Warner, B. D., et al. 2009, A&A, 493, 291 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ďurech, J., Sidorin, V., & Kaasalainen, M. 2010, A&A, 513, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ďurech, J., Kaasalainen, M., Herald, D., et al. 2011, Icarus, accepted [Google Scholar]
 Farinella, P., Vokrouhlicky, D., & Hartmann, W. K. 1998, Icarus, 132, 378 [NASA ADS] [CrossRef] [Google Scholar]
 Hill, G. 1982, Publ. Dom. Astrophys. Obs. Victoria BC, 16, 67 [Google Scholar]
 Ito, T., & Yoshida, F. 2007, PASJ, 59, 269 [NASA ADS] [Google Scholar]
 Johansen, A., & Lacerda, P. 2010, MNRAS, 404, 475 [NASA ADS] [Google Scholar]
 Kaasalainen, M. 2004, A&A, 422, L39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kaasalainen, M., & Torppa, J. 2001, Icarus, 153, 24 [NASA ADS] [CrossRef] [Google Scholar]
 Kaasalainen, M., Torppa, J., & Muinonen, K. 2001, Icarus, 153, 37 [NASA ADS] [CrossRef] [Google Scholar]
 Kaasalainen, M., Mottola, S., & Fulchignoni, M. 2002, Asteroids III, 139 [Google Scholar]
 Kryszczyńska, A., La Spina, A., Paolicchi, P., et al. 2007, Icarus, 192, 223 [NASA ADS] [CrossRef] [Google Scholar]
 La Spina, A., Paolicchi, P., Kryszczyńska, A., & Pravec, P. 2004, Nature, 428, 400 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Lagerkvist, C., Barucci, M. A., Capria, M. T., et al. 1987, Asteroid photometric catalogue., ed. C.I. Lagerkvist, M. A. Barucci, M. T. Capria, M. Fulchignoni, L. Guerriero, E. Perozzi, & V. Zappala [Google Scholar]
 Magnusson, P. 1986, Icarus, 68, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Magnusson, P. 1990, Icarus, 85, 229 [NASA ADS] [CrossRef] [Google Scholar]
 Morbidelli, A., & Vokrouhlický, D. 2003, Icarus, 163, 120 [NASA ADS] [CrossRef] [Google Scholar]
 Nesvorný, D., & Bottke, W. F. 2004, Icarus, 170, 324 [NASA ADS] [CrossRef] [Google Scholar]
 Piironen, J., Lagerkvist, C., Torppa, J., Kaasalainen, M., & Warner, B. 2001, in BAAS, 33, 1562 [Google Scholar]
 Pilcher, F. 2008, Minor Planet Bulletin, 35, 135 [NASA ADS] [Google Scholar]
 Pilcher, F. 2009, Minor Planet Bulletin, 36, 133 [Google Scholar]
 Pilcher, F. 2010, Minor Planet Bulletin, 37, 167 [NASA ADS] [Google Scholar]
 Pilcher, F., & Cooney, J. W. R. 2006, Minor Planet Bulletin, 33, 13 [NASA ADS] [Google Scholar]
 Pilcher, F., & Jardine, D. 2009, Minor Planet Bulletin, 36, 52 [NASA ADS] [Google Scholar]
 Pilcher, F., & Oey, J. 2010, Minor Planet Bulletin, 37, 144 [NASA ADS] [Google Scholar]
 Pravec, P., & Harris, A. W. 2000, Icarus, 148, 12 [NASA ADS] [CrossRef] [Google Scholar]
 Pravec, P., Harris, A. W., & Michalowski, T. 2002, Asteroids III, 113 [Google Scholar]
 Pravec, P., Harris, A. W., Vokrouhlický, D., et al. 2008, Icarus, 197, 497 [NASA ADS] [CrossRef] [Google Scholar]
 Rubincam, D. P. 2000, Icarus, 148, 2 [NASA ADS] [CrossRef] [Google Scholar]
 Scheeres, D. J., & Mirrahimi, S. 2007, in BAAS, 38, 416 [Google Scholar]
 Skoglöv, E., & Erikson, A. 2002, Icarus, 160, 24 [NASA ADS] [CrossRef] [Google Scholar]
 Slivan, S. M. 2002, Nature, 419, 49 [Google Scholar]
 Slivan, S. M., & Molnar, L. A. 2010, in BAAS, 42, 1036 [Google Scholar]
 Slivan, S. M., Binzel, R. P., Crespo da Silva, L. D., et al. 2003, Icarus, 162, 285 [NASA ADS] [CrossRef] [Google Scholar]
 Slivan, S. M., Binzel, R. P., Kaasalainen, M., et al. 2009, Icarus, 200, 514 [NASA ADS] [CrossRef] [Google Scholar]
 Stephens, R. D. 2003, Minor Planet Bulletin, 30, 46 [NASA ADS] [Google Scholar]
 Stephens, R. D. 2004, Minor Planet Bulletin, 31, 40 [NASA ADS] [Google Scholar]
 Stephens, R. D. 2005a, Minor Planet Bulletin, 32, 82 [NASA ADS] [Google Scholar]
 Stephens, R. D. 2005b, Minor Planet Bulletin, 32, 27 [NASA ADS] [Google Scholar]
 Stephens, R. D. 2005c, Minor Planet Bulletin, 32, 2 [NASA ADS] [Google Scholar]
 Stephens, R. D., & Warner, B. D. 2004, Minor Planet Bulletin, 31, 3 [NASA ADS] [Google Scholar]
 Tedesco, E. F., Noah, P. V., Noah, M., & Price, S. D. 2002, AJ, 123, 1056 [Google Scholar]
 Vokrouhlický, D., Nesvorný, D., & Bottke, W. F. 2003, Nature, 425, 147 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Vokrouhlický, D., Nesvorný, D., & Bottke, W. F. 2006, Icarus, 184, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Warner, B. D. 2002, Minor Planet Bulletin, 29, 74 [NASA ADS] [Google Scholar]
 Warner, B. D. 2003, Minor Planet Bulletin, 30, 21 [NASA ADS] [Google Scholar]
 Warner, B. D. 2005a, Minor Planet Bulletin, 32, 54 [NASA ADS] [Google Scholar]
 Warner, B. D. 2005b, Minor Planet Bulletin, 32, 26 [NASA ADS] [Google Scholar]
 Warner, B. D. 2011a, Minor Planet Bulletin, 38, 25 [NASA ADS] [Google Scholar]
 Warner, B. D. 2011b, Minor Planet Bulletin, 38, 52 [NASA ADS] [Google Scholar]
 Warner, B. D. 2011c, Minor Planet Bulletin, 38, 63 [NASA ADS] [Google Scholar]
 Warner, B. D., Harris, A. W., & Pravec, P. 2009, Icarus, 202, 134 [NASA ADS] [CrossRef] [Google Scholar]
 Yoshida, F., Dermawan, B., Ito, T., et al. 2004, PASJ, 56, 1105 [NASA ADS] [Google Scholar]
Online material
List of new asteroid models derived from combined data sets or sparse data alone.
List of partial models derived from combined data sets.
List of improved asteroid models that were originally designated in DAMIT as “preliminary”.
Observations used for successful model determination that are not included in the UAPC.
All Tables
Comparison of estimated characteristics of residuals for 13 selected observatories: mode, “FWHM” and median.
Test of the hypothesis that the observed pole distributions and debiased latitude distribution are drown from uniform distributions (a χ^{2}test).
List of new asteroid models derived from combined data sets or sparse data alone.
List of improved asteroid models that were originally designated in DAMIT as “preliminary”.
Observations used for successful model determination that are not included in the UAPC.
All Figures
Fig. 1 Four histograms of residuals comparing actual data against the model given by Eq. (1) of all sparse lightcurves belonging to the following observatories: 689 USNO, 691 Steward Observatory, 703 Catalina Sky Survey and E12 Siding Spring Survey. Number of bins is ~, where N is the total number of sparse lightcurves used for histogram construction. 

In the text 
Fig. 2 Asteroid (810) Atossa: shape model (left panel) and an example of a “mock” shape model (right panel). 

In the text 
Fig. 3 Relation between the proper semimajor axis and the proper eccentricity for asteroids in our sample and for the first 100 000 numbered asteroids for comparison. Main resonances are shown by dotted lines. Prograde rotators are plotted with blue circles and retrograde rotators with red squares. The horizontal lines represent for each asteroid its estimated past drift (i.e. where the asteroid came from) during the collisional lifetime computed with Eq. (7). Proper elements are from the AstDyS database. 

In the text 
Fig. 4 a) Shows the ecliptic latitude distribution of all MBAs in our sample except the Koronis family members. The width of bins in latitude β corresponds to similar surfaces on the (λ, β)sphere, so the bins are equidistant in sinβ. In b), the longitude distribution of all MBAs is plotted, again except the Koronis family members. The longitude bins are equidistant in λ. In the top right corners, there are the probability values of the χ^{2}tests (that the observed distributions are randomly drown from a uniform distribution, see Table 2). 

In the text 
Fig. 5 Dependence of the ecliptic latitude β (plotted as sinβ) of the pole direction on the models diameter D. The gray box shows the gap of small latitudes for asteroids with D < 30 km. Members of the Koronis cluster are plotted with crosses. 

In the text 
Fig. 6 Histograms showing the observed latitude and longitude distributions of MBAs (except Koronis cluster members) for different size ranges. a) Shows the latitude distribution for asteroids with diameters larger than 60 km, c) for asteroids with diameters in the range of 30–60 km, e) for asteroids with diameters smaller than 30 km, and similarly b), d) and f) for longitudes. The width of bins in the latitude β and longitude lambda λ corresponds to similar surfaces on the (λ, β)sphere, so the bins are equidistant in sinβ and λ. In the top right corners, there are the probability values of the χ^{2}tests (that the observed distributions are randomly drown from a uniform distribution, see Table 2). 

In the text 
Fig. 7 Dependence of the ecliptic latitude β (plotted as sinβ) of the pole direction on the models spin rate, a) for asteroids with diameter D < 30 km and b) for asteroids with diameter D > 60 km. 

In the text 
Fig. 8 Histograms of the angular differences between the generated and derived ecliptic latitudes and longitudes of the pole directions for all successfully derived synthetic models: a) β_{comp} − β_{gen}, where β_{gen} and β_{comp} are generated and computed ecliptic latitudes and similarly b) (λ_{comp} − λ_{gen})cosβ_{gen} for ecliptic longitudes. 

In the text 
Fig. 9 a) Shows the simulated ecliptic latitude distribution of all successfully derived models. In b), the debiased observed latitude distribution is plotted. The bins are equidistant in sinβ. In the top right corner, there is the probability value of the χ^{2}test (that the observed distribution is drown from a uniform one, see Table 2). 

In the text 
Fig. 10 Histograms showing the simulated and corrected (debiased) observed latitude distributions of MBAs for different size groups, excluding the Koronis cluster members. a) Shows the simulated latitude distribution for asteroids with the diameters larger than 60 km, c) for asteroids with the size in range of 30–60 km, e) for asteroids with diameters smaller than 30 km. Similarly b), d) and f) show the observed latitude distributions corrected by the bias of the LI method. The bins are equidistant in sinβ. In the top right corners, there are the probability values of the χ^{2}tests (that the observed distributions are drown from a uniform distribution, see Table 2). 

In the text 
Fig. 11 Left panel: the distribution of the ecliptic latitudes β for the observed asteroids (thick line) and 100 synthetic samples (generated with a different random seed) at time t = 0 (thin lines). Right panel: steadystate synthetic latitude distributions at time t = 4 Gyr (thin lines), evolved by the YORP effect, collisions, mass shedding and spinorbit resonances. We also applied an observational bias to it. The steady state was reached already within ≃ 1 Gyr. 

In the text 
Fig. 12 Steadystate synthetic latitude distributions for three different size ranges: D > 60 km (top), 30 km < D < 60 km (middle) and D < 30 km (bottom). The synthetic distributions are plotted by thin gray lines while the observed distribution by a thick black line. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext 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 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.