Issue 
A&A
Volume 587, March 2016



Article Number  A48  
Number of page(s)  6  
Section  Planets and planetary systems  
DOI  https://doi.org/10.1051/00046361/201527573  
Published online  16 February 2016 
Asteroid models from the Lowell photometric database^{⋆}
^{1}
Astronomical Institute, Faculty of Mathematics and Physics,
Charles University in Prague, V
Holešovičkách 2, 180 00
Prague 8, Czech Republic
email:
durech@sirrah.troja.mff.cuni.cz
^{2}
Centre National d’Études Spatiales, 2 place Maurice
Quentin, 75039
Paris Cedex 01,
France
^{3}
Laboratoire Lagrange, UMR 7293, Université de la Côte d’Azur,
CNRS, Observatoire de la Côte d’Azur, Bd de l’Observatoire,
CS 34229, 06304
Nice Cedex 04,
France
^{4}
Astronomical Observatory Institute, Faculty of Physics, A.
Mickiewicz University, Słoneczna
36, 60286
Poznań,
Poland
Received: 16 October 2015
Accepted: 18 December 2015
Context. Information about shapes and spin states of individual asteroids is important for the study of the whole asteroid population. For asteroids from the main belt, most of the shape models available now have been reconstructed from diskintegrated photometry by the lightcurve inversion method.
Aims. We want to significantly enlarge the current sample (~350) of available asteroid models.
Methods. We use the lightcurve inversion method to derive new shape models and spin states of asteroids from the sparseintime photometry compiled in the Lowell Photometric Database. To speed up the timeconsuming process of scanning the period parameter space through the use of convex shape models, we use the distributed computing project Asteroids@home, running on the Berkeley Open Infrastructure for Network Computing (BOINC) platform. This way, the periodsearch interval is divided into hundreds of smaller intervals. These intervals are scanned separately by different volunteers and then joined together. We also use an alternative, faster, approach when searching the bestfit period by using a model of triaxial ellipsoid. By this, we can independently confirm periods found with convex models and also find rotation periods for some of those asteroids for which the convexmodel approach gives too many solutions.
Results. From the analysis of Lowell photometric data of the first 100 000 numbered asteroids, we derived 328 new models. This almost doubles the number of available models. We tested the reliability of our results by comparing models that were derived from purely Lowell data with those based on dense lightcurves, and we found that the rate of falsepositive solutions is very low. We also present updated plots of the distribution of spin obliquities and pole ecliptic longitudes that confirm previous findings about a nonuniform distribution of spin axes. However, the models reconstructed from noisy sparse data are heavily biased towards more elongated bodies with high lightcurve amplitudes.
Conclusions. The Lowell Photometric Database is a rich and reliable source of information about the spin states of asteroids. We expect hundreds of other asteroid models for asteroids with numbers larger than 100 000 to be derivable from this data set. More models will be able to be reconstructed when Lowell data are merged with other photometry.
Key words: minor planets, asteroids: general / methods: data analysis / techniques: photometric
Tables 1 and 2 are only available at the CDS via anonymous ftp to cdsarc.ustrasbg.fr (130.79.128.5) or via http://cdsarc.ustrasbg.fr/vizbin/qcat?J/A+A/587/A48
© ESO, 2016
1. Introduction
Large allsky surveys like Catalina, PanSTARRS, etc. image the sky every night to discover new asteroids and detect those that are potentially hazardous. The main output of these surveys is a steadily increasing number of asteroids with known orbits. Apart from astrometry that is used for orbit computation, these surveys also produce photometry of asteroids. This photometry contains, in principle, information about asteroid rotation, shape, and surface properties. However, because of its poor quality (when compared with a dedicated photometric measurements of a single asteroid) the signal corresponding to asteroid’s rotation is usually drowned in noise and systematic errors. However, there have been recent attempts to use sparseintime photometry to reconstruct the shape of asteroids. Kaasalainen (2004) has shown that sparse photometry can be used to solve the lightcurve inversion problem and further simulations confirm this (Ďurech et al. 2005, 2007). Afterwards, real sparse data were used either alone or in combination with dense lightcurves and new asteroid models were derived (Ďurech et al. 2009; Cellino et al. 2009; Hanuš et al. 2011, 2013c). The aim of these efforts was to derive new unique models of asteroids, i.e., their sidereal rotation periods, shapes, and direction of spin axis.
Another approach to utilize sparse data was to look for changes in the mean brightness as a function of the aspect angle, which led to estimations of spinaxis longitudes for more than 350 000 asteroids (Bowell et al. 2014) from the socalled Lowell Observatory photometric database (Oszkiewicz et al. 2011).
In this paper, we show that the Lowell photometric data set can also be used for solving the full inversion problem. By processing Lowell photometry for the first 100 000 numbered asteroids, we derived new shapes and spin states for 328 asteroids, which almost doubles the number of asteroids for which the photometrybased physical model is known.
We describe the data, the inversion method, and the reliability tests in Sect. 2, the results in Sect. 3, and we conclude in Sect. 4.
2. Method
The lightcurve inversion method of Kaasalainen et al. (2001) that we applied was reviewed by Kaasalainen et al. (2002) and more recently by Ďurech et al. (2016a). We used the same implementation of the method as Hanuš et al. (2011), where the reader is referred to for details. Here we describe only the general approach and the details specific for our work.
2.1. Data
As the data source, we used the Lowell Observatory photometric database (Bowell et al. 2014). This is photometry provided to Minor Planet Centre (MPC) by 11 of the largest surveys that were recalibrated in the Vband using the accurate photometry of the Sloan Digital Sky Survey. Details about the data reduction and calibration can be found in Oszkiewicz et al. (2011). The data are available for about ~330 000 asteroids. Typically, there are several hundreds of photometric points for each asteroid. The length of the observing interval is ~10–15 yr. The largest amount of data is for the lownumbered asteroids and decreases with increasing asteroid numbers. For example, the average number of data points is ~480 for asteroids with number < 10 000 and ~45 for those > 300 000. The accuracy of the data is around 0.10−0.20 mag.
For each asteroid and epoch of observation, we computed the asteroidcentric vectors towards the Sun and the Earth in the Cartesian ecliptic coordinate frame – these were needed to compute the illumination and viewing geometry in the inversion code.
2.2. Convex models
To derive asteroid models from the optical data, we used the lightcurve inversion method of Kaasalainen & Torppa (2001) and Kaasalainen et al. (2001), the same way as Hanuš et al. (2011). Essentially, we searched for the bestfit model by densely scanning the rotation period parameter space. We decided to search in the interval of 2−100 h. The lower limit roughly corresponds to the observed rotation limit of asteroids larger than ~150 m (Pravec et al. 2002), the upper limit was set arbitrarily to cover most of the rotation periods for asteroids determined so far. For each trial period, we started with ten initial pole directions that were isotropically distributed on a sphere. This turned out to be enough not to miss any local minimum in the pole parameter space. In each period run, we recorded the period and χ^{2} value that correspond to the best fit. Then we looked for the global minimum of χ^{2} on the whole period interval and tested the uniqueness and stability of this globally best solution (see details in Sect. 2.6).
For a typical data set, the number of trial periods is 200 000−300 000, which takes about a month on one CPU. Because the number of asteroids we wanted to process was ~100 000, the only way to finish the computations in a reasonable time was to use tens of thousands of CPUs. For this task, we used the distributed computing project Asteroids@home^{1}.
2.3. Asteroids@home
Asteroids@home is a volunteerbased computing project built on the Berkeley Open Infrastructure for Network Computing (BOINC) platform. Because the scanning of the period parameter space is the socalled embarrassingly parallel problem, we divided the whole interval of 2−100 h into smaller intervals (typically hundreds), which were searched individually on the computers of volunteers connected to the project. The units sent to volunteers had about the same CPUtime demand. Results from volunteers were sent back to the BOINC server and validated. When all units belonging to one particular asteroid were ready, they were connected and the global minimum was found. The technical details of the project are described in Ďurech et al. (2015)
2.4. Ellipsoids
To find the rotation period in sparse data, we also used an alternative approach that was based on the triaxial ellipsoid shape model and a geometrical lightscattering model (Kaasalainen & Ďurech 2007). Its advantage is that it is much faster than using convex shapes because the brightness can be computed analytically (it is proportional to the illuminated projected area, Ostro & Connelly 1984). On top of that, contrary to the convex modelling, all shape models automatically fulfill the physical condition of rotating along the principal axis with the largest momentum of inertia. The accuracy of this simplified model is sufficient to reveal the correct rotation period as a significant minimum of χ^{2} in the period parameter space. That period is then used as a start point for the convex inversion for the final model. In many cases when the convex models gives many equally good solutions with different periods, this method provides a unique and correct rotation period.
2.5. Restricted period interval
As mentioned above, the interval for period search was 2–100 h. However, for many asteroids, their rotation period is known from observations of their lightcurves. The largest database of asteroid rotation periods is the Lightcurve Asteroid Database (LCDB) compiled by Warner et al. (2009) and regularly updated^{2}. If we take information about the rotation period as an a priori constraint, we can narrow the interval of possible periods and significantly shrink the parameter space. For this purpose, we used only reliable period determinations from LCDB with quality codes U equal to 3, 3, or 2+. However, even for these quality codes, the LCDB period can be wrong (for examples see Marciniak et al. 2015) resulting in a wrong shape model. For quality codes 3 and 3, we restricted the search interval to P ± 0.05P, where P was the rotation period reported in LCDB. Similarly for U equal to 2+, we restricted the search interval to P ± 0.1P. We applied this approach to both convex and ellipsoidbased period search.
2.6. Tests
For each periodogram, there is formally one best model that corresponds to the period with the lowest χ^{2}. However, the global minimum in χ^{2} has to be significantly deeper than all other local minima to be considered as a reliable solution, rather than just a random fluctuation. We could not use formal statistical tools to decide whether the lowest χ^{2} value is statistically significant or not, because the data were also affected by systematic errors. Instead, to select only robust models, we set up several tests, which each model had to pass to be considered as a reliable model.

1.
The lowest χ^{2} corresponding to the rotation period P_{min} is at least 5% lower than all other χ^{2} values for periods outside the P_{min}± interval, where ΔT is the time span of observations (Kaasalainen 2004). The value of 5% was chosen such that the number of unique models was as large as possible while keeping the number of false positive solutions very low (~1%). The comparison was done with respect to models in DAMIT (see Sect. 3.1).

2.
When using convex models for scanning the period parameter space, we ran the period search for two resolutions of the convex model – the degree and order of the harmonics series expansion that parametrized the shape was three or six. The periods P_{min} corresponding to these two resolutions had to agree within their errors (and both had to pass the test nr. 1).

3.
Because we realized that P_{min} ≳ 20 h often produced false positive solutions, we accepted only models with P_{min} shorter than 20 h (when there was no information about the rotation period from LCDB).

4.
For a given P_{min}, there are no more than two distinct (farther than 30° apart) pole solutions with χ^{2} at least 5% deeper than other poles.

5.
Because of the geometry limited close to the ecliptic plane, two models that have the same pole latitudes β and pole longitudes λ that are different by 180° provide the same fit to diskintegrated data, and they cannot be distinguished from each other (Kaasalainen & Lamberg 2006). Therefore we accepted only such solutions that fulfilled the condition that if there were two pole directions (λ_{1},β_{1}) and (λ_{2},β_{2}), the difference in ecliptic latitudes  β_{1} − β_{2}  has to be less than 50° and the difference in ecliptic longitudes mod (  λ_{1} − λ_{2}  ,180°) has to be larger than 120°.

6.
The ratio of the moment of inertia along the principal axis to that along the actual rotation axis should be less than 1.1. Otherwise the model is too elongated along the direction of the rotation axis and it is not considered a realistic shape.

7.
For each asteroid that passed the above test, we created a bootstrapped lightcurve data set by randomly selecting the same number of observations from the original data set. This new data set was processed the same way as the original one (using either convex shapes or ellipsoids for the period search) and the model was considered stable only if the bestfit period P_{min} from the bootstrapped data agreed with that from the original data.

8.
We also visually checked all shape models, periodograms, and fits to the data to be sure that the shape model looked realistic and that there were no clear problems with the data and residuals. In some rare cases we rejected models that formally fitted the data, passed all the test, but were unrealistically elongated or flat.
Fig. 1 Histogram of differences between the pole directions of models derived from Lowell data and those archived in DAMIT. 
3. Results
3.1. Comparison with independent models
From all ~600 models that successfully passed the tests described in Sect. 2.6, some were already modeled from other photometric data and the models were stored in the Database of Asteroid Models from Inversion Techniques (DAMIT^{3}, Ďurech et al. 2010). For this subset, we could compare our results from an inversion of Lowell data with independent models (assumed to be reliable) from DAMIT. In total, there were 279 models in DAMIT for comparison. For these models, we computed the difference between the DAMIT and Lowell rotation periods and also the difference between the pole directions. Out of this set, almost all (275 models) have the same rotation periods (within the uncertainties) and the pole differences < 50° of arc. The histogram of pole differences between DAMIT and our models is shown in Fig. 1. Although there are some asteroids for which we got differences as large as ~40−50°, the mean value is 15° and the median 13°, which can be interpreted as a typical error in the pole determination that was based on Lowell data, assuming that the poles from DAMIT have smaller errors (typically 5−10°).
As an example of the difference between shape models, we show results for asteroid (63) Ausonia. In Fig. 2, we compare our shape model, which we derived from Lowell sparse photometry, with that obtained by inversion of dense lightcurves (Torppa et al. 2003). In general, the shapes derived from sparse photometry are more angular than those derived from dense lightcurves and often have artificial sharp edges.
The four asteroids (5) Astraea, (367) Amicita, (540) Rosamunde, and (4954) Eric, for which we got different solutions to DAMIT, are discussed below. We also discuss the five asteroids – (1753) Mieke, (2425) Shenzen, (6166) Univsima, (11958) Galiani, and (12753) Povenmire – for which there is no model in DAMIT, but the period we derived from the Lowell data does not agree with the data in LCDB.
Fig. 2 Comparison between the shape model of (63) Ausonia reconstructed from Lowell sparse data (top) and from dense lightcurves (bottom). 
(5) Astraea.
From Lowell data, we got two pole directions (λ,β) = (121°, −20°) and (296°, −15°), the former being about 60° away from the DAMIT model of Hanuš et al. (2013b) with the pole (126°,40°). The DAMIT model agrees with the adaptive optics data as well as with the occultation silhouette from 2008 and it is not clear why there is so large a difference in the pole direction, while the rotation periods are the same and the number of Lowell photometric points is also large (447 points).
(367) Amicita.
The model derived from Lowell data has two pole solutions (17°, −52°) and (194°, −45°) and rotation period of 5.05578 h, while the DAMIT model of Hanuš et al. (2011) has prograde rotation with poles of (21°,32°) and (203°,38°), with a significantly different period of 5.05502 h. However, the DAMIT model is based on sparse data from US Naval Observatory and Catalina and only two pieces of lightcurve by Wisniewski et al. (1997) and it might not be correct.
(540) Rosamunde.
Although the periodogram obtained with the convex model approach shows a minimum for 9.34780 h – the same as the DAMIT model of Hanuš et al. (2013a) – this minimum was not deep enough to pass the test nr. 1. However, the secondbest minimum for a convex model at 7.82166 h appeared as the best solution for the ellipsoid approach and passed all tests leading to a wrong model.
(4954) Eric.
The DAMIT model of Hanuš et al. (2013c) has a pole direction of (86°, −55°), which is almost exactly opposite to our value of (261°,70°). Moreover, even the rotation periods are different by about 0.0003 h, which is more than the uncertainty interval.
(1753) Mieke.
The rotation period of 8.9 h was determined by Lagerkvist (1978) from two (1.5 and 5 h) noisy lightcurves. Given the quality of the data, this period is not in contradiction with our value of 10.19942 h.
(2425) Shenzen.
The rotation period of 14.715 ± 0.012 h was determined by Hawkins & Ditteon (2008). Our value of 9.83818 h is close to 2/3 of their. In the periodogram, there is no significant minimum around 14.7 h.
(6166) Univsima.
The lightcurve is published online in the database of R. Behrend^{4}. However, the period of 9.6 h is based on only 12 points, which covers about half of the reported period, so we think that this preliminary result is not in contradiction with our period of ~11.4 h.
(11958) Galiani.
This asteroid was observed by Clark (2014), who determined the period 9.8013 ± 0.0023 h, which does not agree with our value of 8.24720 h. The reason is not clear, because the data of Clark (2014) seem to fit this period correctly. We do not see any significant minimum in χ^{2} near 9.8 h in the periodogram.
(12753) Povenmire.
The period of 12.854 h reported in the LCDB is based on the observations of Gary (2004). However, according to the same author^{5}, the correct rotation period that is based on observations from 2010 is 17.5752 ± 0.0008 h, which agrees with our value.
In summary, the frequency of false positive solutions that pass all reliability tests seems to be sufficiently low, around a few percent. However, the sample of models in DAMIT that we use for comparison is itself biased against lowamplitude longperiod asteroids (Marciniak et al. 2015), so the real number of false positive solutions might be higher.
3.2. New models
After applying all the tests described in Sect. 2.6, we selected only those asteroids with no model in DAMIT for publication. These are listed in Tables 1 (models from full interval 2−100 h) and 2 (models derived from a restricted period interval). The tables list the pole direction(s) (one or two models), the sidereal rotation period (with uncertainty corresponding to the order of the last decimal place). The C/E code means the method by which P_{min} was found – convex models (C) or ellipsoids (E). In some cases, both methods independently gave the same value of P_{min} (then CE code). All new shape models and the photometric data are available in DAMIT.
For some of these asteroids, Hanuš et al. (2016) obtained independent models by applying the same lightcurve inversion method on sparse data, which they combined with dense lightcurves. These asteroids (not yet published in DAMIT) are marked by asterisk in the Tables 1 and 2. For all of them (56 in total), our rotation periods agree with those of Hanuš et al. (2016) within their uncertainties and pole directions differ by 10–20 degrees on average. By way of comparison, this is a similar result to the DAMIT models (Sect. 3.1, Fig. 1) and it independently confirms the reliability of our models based on only Lowell data.
Fig. 3 Distribution of pole obliquity ε as a function of size for 575 mainbelt and 13 nearEarth asteroids. 
Fig. 4 Histograms of the distribution of pole obliquities ε for asteroids with diameters < 30 km and > 60 km, respectively. 
3.3. Statistics of pole directions
Together with models from DAMIT, we now have a sample of shape models for 717 asteroids (685 MBAs, 13 NEAs, 10 Marscrossers, 7 Hungarias, 1 Hilda, and 1 Trojan). The statistical analysis of the pole distribution confirms the previous findings. Namely, the distribution of spin directions is not isotropic (Kryszczyńska et al. 2007). Moreover, the distribution of pole obliquities (an angle between the spin vector and the normal to the orbital plane) depends on the size of an asteroid. We plot the dependence of obliquity on the size in Fig. 3 for mainbelt (MBAs) and nearEarth (NEAs) asteroids. There is a clear trend of smaller asteroids clustering towards extreme values of obliquity. This was explained by Hanuš et al. (2011) as YORPinduced evolution of spins (Hanuš et al. 2013c). The distribution of obliquities is not symmetric around 90° (Fig. 4). As noticed by Hanuš et al. (2013c), the retrograde rotators are more concentrated to −90°, probably because prograde rotators are affected by resonances. For larger asteroids, there is an excess of prograde rotators that might be primordial (Kryszczyńska et al. 2007; Johansen & Lacerda 2010).
However, the current sample of asteroid models is far from being representative of the whole asteroid population. Because the period search in sparse data is strongly dependent on the lightcurve amplitude – the larger the amplitude the easier is to detect the correct rotation period in noisy data – more elongated asteroids are reconstructed more easily than spherical ones. That is why almost all the asteroids listed in Tables 1 and 2 have large amplidudes of ≳ 0.3 mag. The lightcurve inversion (based mostly or exclusively on sparse data) is also less efficient for asteroids with poles close to the ecliptic plane because, during some apparitions, we observe them almost poleon, thus with very small amplitudes. This bias in the method was estimated to be of the order of several tens percent (Hanuš et al. 2011). A much higher discrepancy (factor 3–4) in the successfully recovered pole directions between poles closeto and perpendicularto the ecliptic was found by SantanaRos et al. (2015). But even such a large selection effect cannot fully explain the significant “gap” for obliquities between 60–120°. To clearly show how the unbiased distribution of pole obliquities looks like, we would have to carry out an extensive simulation on a synthetic population with realistic systematic and random errors to see the bias that is induced by the method, shape, and geometry. This sort of simulation would be more computationally demanding than processing real data from the Lowell database. Therefore, we postpone this investigation for a future paper.
For nearEarth asteroids, the excess of retrograde rotators can be explained by the Yarkovskyinduced delivery mechanism from the main belt through resonances (La Spina et al. 2004), although the number of NEA models in our sample is too small for any reliable statistics.
We also see a clear deviation from a uniform distribution of pole longitudes in Fig. 5. Because of ambiguity in pole direction (often there are two solutions with similar latitudes and the difference in longitudes of about 180°), we plotted the distribution modulo 180°. The histogram shows an excess of longitudes around 50−100°. This was already announced by Bowell et al. (2014), who processed the Lowell data set using a different approach, estimated spinaxis longitudes for more than 350 000 asteroids, and revealed an excess of longitudes at 30−110° and a paucity at 120−180°. The explanation of this phenomenon remains unclear.
Fig. 5 Histograms of the distribution of pole longitude λ for 685 mainbelt asteroids. 
4. Conclusions
The new models presented in this paper significantly enlarge the sample of asteroids for which their spin axis direction and approximate shape are known. Because these models are based on a limited number of data points, the shapes have to be interpreted as only approximations of the real shapes of asteroids. Also the pole directions need to be refined with more data if one is interested in a particular asteroid. However, as an ensemble, the models can be used in future statistical studies of asteroid spins, for example.
We believe that this is only the beginning of a mass production of shape and spin models from sparse photometry. Although the number of models derivable from the Lowell Observatory photometric database is small compared to the total number of asteroids, the potential of Lowell photometry consists in its combination with other data. Even a priori information about the rotation period shrinks the parameter space that has to be scanned, and a local minimum in a large parameters space becomes a global minimum on a restricted interval. Of course, the reliability of this type of model depends critically on the reliability of the period. Lowell photometry can be combined with dense lightcurves that constrain the rotation period. This way, models for about 250 asteroids were derived recently by Hanuš et al. (2016), some of which confirm the models presented in this paper. The database of asteroid rotation periods has been increased dramatically by Waszczak et al. (2015) – their data can also be combined with Lowell photometry, and we expect that other hundreds of models will be reconstructed from this data set. Another promising approach is the combination of sparse photometry with data from the Widefield Infrared Survey Explorer (WISE) mission (Wright et al. 2010). Although WISE data were observed in midinfrared wavelengths, Ďurech et al. (2016b) showed that thermally emitted flux can be treated as reflected light to derive the correct rotation period and the shape and spin model. This opens up a new possibility, because both Lowell and WISE data are available for tens of thousands of asteroids.
In general, the combination of more data sources is always better than using them separately. By using Lowell photometry with dense lightcurves, WISE data, photometry from Gaia, etc., the number of available models will increase and the statistical studies of spin and shape distribution will become more robust, being based on larger sets of models. Nevertheless, any inference based on the models derived from lightcurves (and sparse lightcurves in particular) has to take into account that the sample of models is biased against more spherical shapes with low lightcurve amplitudes and poles near the plane of ecliptic.
Acknowledgments
We would not be able to process data for hundreds of thousands of asteroids without the help of tens of thousand of volunteers who joined the Asteroids@home BOINC project and provided computing resources from their computers. We greatly appreciate their contribution. The work of J.Ď. was supported by the grant 1504816S of the Czech Science Foundation. J.H. greatly appreciates the CNES postdoctoral fellowship program. J.H. was supported by the project under the contract 11BS56008 (SHOCKS) of the French Agence National de la Recherche (ANR). D.O. was supported by the grant NCN 2012/S/ST9/00022 Polish National Science Center. We thank the referee A. Marciniak for providing constructive comments that improved the contents of this paper.
References
 Bowell, E., Oszkiewicz, D. A., Wasserman, L. H., et al. 2014, Meteor. Planet. Sci., 49, 95 [NASA ADS] [CrossRef] [Google Scholar]
 Cellino, A., Hestroffer, D., Tanga, P., Mottola, S., & Dell’Oro, A. 2009, A&A, 506, 935 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Clark, M. 2014, Minor Planet Bull., 41, 178 [NASA ADS] [Google Scholar]
 Ďurech, J., Grav, T., Jedicke, R., Kaasalainen, M., & Denneau, L. 2005, Earth Moon Planets, 97, 179 [Google Scholar]
 Ďurech, J., Scheirich, P., Kaasalainen, M., et al. 2007, in Near Earth Objects, our Celestial Neighbors: Opportunity and Risk, eds. A. Milani, G. B. Valsecchi, & D. Vokrouhlický (Cambridge: Cambridge University Press), 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., Hanuš, J., & Vančo, R. 2015, Astron. Comp., 13, 80 [Google Scholar]
 Ďurech, J., Carry, B., Delbo, M., Kaasalainen, M., & Viikinkoski, M. 2016a, in Asteroids IV, eds. P. Michel, F. DeMeo, & W. Bottke (Tucson: University of Arizona Press), in press [Google Scholar]
 Ďurech, J., Hanuš, J., AliLagoa, V., Delbo, M., & Oszkiewicz, D. 2016b, in Proc. IAU Symp. 318, eds. S. Chesley, R. Jedicke, A. Morbidelli, & D. Farnocchia (Cambridge: Cambridge University Press) [Google Scholar]
 Gary, B. L. 2004, Minor Planet Bull., 31, 56 [NASA ADS] [Google Scholar]
 Hanuš, J., Ďurech, J., Brož, M., et al. 2011, A&A, 530, A134 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hanuš, J., Brož, M., Durech, J., et al. 2013a, A&A, 559, A134 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hanuš, J., Marchis, F., & Ďurech, J. 2013b, Icarus, 226, 1045 [NASA ADS] [CrossRef] [Google Scholar]
 Hanuš, J., Ďurech, J., Brož, M., et al. 2013c, A&A, 551, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hanuš, J., Ďurech, J., Oszkiewicz, D. A., et al. 2016, A&A, 586, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hawkins, S., & Ditteon, R. 2008, Minor Planet Bull., 35, 1 [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., & Ďurech, J. 2007, in Near Earth Objects, our Celestial Neighbors: Opportunity and Risk, eds. A. Milani, G. B. Valsecchi, & D. Vokrouhlický (Cambridge: Cambridge University Press), 151 [Google Scholar]
 Kaasalainen, M., & Lamberg, L. 2006, Inverse Problems, 22, 749 [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., & Fulchignomi, M. 2002, in Asteroids III, eds. W. F. Bottke, A. Cellino, P. Paolicchi, & R. P. Binzel (Tucson: University of Arizona Press), 139 [Google Scholar]
 Kryszczyńska, A., La Spina, A., Paolicchi, P., et al. 2007, Icarus, 192, 223 [NASA ADS] [CrossRef] [Google Scholar]
 LaSpina, A., Paolicchi, P., Kryszczyńska, A., & Pravec, P. 2004, Nature, 428, 400 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Lagerkvist, C.I. 1978, A&AS, 31, 361 [NASA ADS] [Google Scholar]
 Marciniak, A., Pilcher, F., Oszkiewicz, D., et al. 2015, Planet. Space Sci., 118, 256 [NASA ADS] [CrossRef] [Google Scholar]
 Ostro, S. J., & Connelly, R. 1984, Icarus, 57, 443 [NASA ADS] [CrossRef] [Google Scholar]
 Oszkiewicz, D., Muinonen, K., Bowell, E., et al. 2011, J. Quant. Spectr. Rad. Transf., 112, 1919 [NASA ADS] [CrossRef] [Google Scholar]
 Pravec, P., Harris, A. W., & Michałowski, T. 2002, in Asteroids III, eds. W. F. Bottke, A. Cellino, P. Paolicchi, & R. P. Binzel (Tucson: University of Arizona Press), 113 [Google Scholar]
 SantanaRos, T., Bartczak, P., Michałowski, T., Tanga, P., & Cellino, A. 2015, MNRAS, 450, 333 [NASA ADS] [CrossRef] [Google Scholar]
 Torppa, J., Kaasalainen, M., Michalowski, T., et al. 2003, Icarus, 164, 346 [NASA ADS] [CrossRef] [Google Scholar]
 Warner, B. D., Harris, A. W., & Pravec, P. 2009, Icarus, 202, 134 [NASA ADS] [CrossRef] [Google Scholar]
 Waszczak, A., Chang, C.K., Ofek, E. O., et al. 2015, AJ, 150, 75 [NASA ADS] [CrossRef] [Google Scholar]
 Wisniewski, W. Z., Michalowski, T. M., Harris, A. W., & McMillan, R. S. 1997, Icarus, 126, 395 [NASA ADS] [CrossRef] [Google Scholar]
 Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
Fig. 1 Histogram of differences between the pole directions of models derived from Lowell data and those archived in DAMIT. 

In the text 
Fig. 2 Comparison between the shape model of (63) Ausonia reconstructed from Lowell sparse data (top) and from dense lightcurves (bottom). 

In the text 
Fig. 3 Distribution of pole obliquity ε as a function of size for 575 mainbelt and 13 nearEarth asteroids. 

In the text 
Fig. 4 Histograms of the distribution of pole obliquities ε for asteroids with diameters < 30 km and > 60 km, respectively. 

In the text 
Fig. 5 Histograms of the distribution of pole longitude λ for 685 mainbelt asteroids. 

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.