Precovery of near-Earth asteroids by a citizen-science project of the Spanish Virtual Observatory

This article describes a citizen-science project conducted by the Spanish Virtual Observatory (SVO) to improve the orbits of near-Earth asteroids (NEAs) using data from astronomical archives. The list of NEAs maintained at the Minor Planet Center (MPC) is checked daily to identify new objects or changes in the orbital parameters of already catalogued objects. Using NEODyS we compute the position and magnitude of these objects at the observing epochs of the 938 046 images comprising the Eigth Data Release of the Sloan Digitised Sky Survey (SDSS). If the object lies within the image boundaries and the magnitude is brighter than the limiting magnitude, then the associated image is visually inspected by the project's collaborators (the citizens) to confirm or discard the presence of the NEA. If confirmed, accurate coordinates and, sometimes, magnitudes are submitted to the MPC. Using this methodology, 3,226 registered users have made during the first fifteen months of the project more than 167,000 measurements which have improved the orbital elements of 551 NEAs (6% of the total number of this type of asteroids). Even more remarkable is the fact that these results have been obtained at zero cost to telescope time as NEAs were serendipitously observed while the survey was being carried out. This demonstrates the enormous scientific potential hidden in astronomical archives. The great reception of the project as well as the results obtained makes it a valuable and reliable tool for improving the orbital parameters of near-Earth asteroids.Comment: 9 pages, 5 figures. Accepted in Astron. Nach

Prediction of transits of solar system objects in Kepler/K2 images: An extension of the Virtual Observatory service SkyBoT

All the fields of the extended space mission Kepler/K2 are located within the ecliptic. Many solar system objects thus cross the K2 stellar masks on a regular basis. We aim at providing to the entire community a simple tool to search and identify solar system objects serendipitously observed by Kepler. The SkyBoT service hosted at IMCCE provides a Virtual Observatory (VO) compliant cone-search that lists all solar system objects present within a field of view at a given epoch. To generate such a list in a timely manner, ephemerides are pre-computed, updated weekly, and stored in a relational database to ensure a fast access. The SkyBoT Web service can now be used with Kepler. Solar system objects within a small (few arcminutes) field of view are identified and listed in less than 10 sec. Generating object data for the entire K2 field of view (14{\deg}) takes about a minute. This extension of the SkyBot service opens new possibilities with respect to mining K2 data for solar system science, as well as removing solar system objects from stellar photometric time-series

Asteroid Models from Multiple Data Sources

In the past decade, hundreds of asteroid shape models have been derived using the lightcurve inversion method. At the same time, a new framework of 3-D shape modeling based on the combined analysis of widely different data sources such as optical lightcurves, disk-resolved images, stellar occultation timings, mid-infrared thermal radiometry, optical interferometry, and radar delay-Doppler data, has been developed. This multi-data approach allows the determination of most of the physical and surface properties of asteroids in a single, coherent inversion, with spectacular results. We review the main results of asteroid lightcurve inversion and also recent advances in multi-data modeling. We show that models based on remote sensing data were confirmed by spacecraft encounters with asteroids, and we discuss how the multiplication of highly detailed 3-D models will help to refine our general knowledge of the asteroid population. The physical and surface properties of asteroids, i.e., their spin, 3-D shape, density, thermal inertia, surface roughness, are among the least known of all asteroid properties. Apart for the albedo and diameter, we have access to the whole picture for only a few hundreds of asteroids. These quantities are nevertheless very important to understand as they affect the non-gravitational Yarkovsky effect responsible for meteorite delivery to Earth, or the bulk composition and internal structure of asteroids.Comment: chapter that will appear in a Space Science Series book Asteroids I

Physical and dynamical properties of the main belt triple asteroid (87) Sylvia

We present the analysis of high angular resolution observations of the triple Asteroid (87) Sylvia collected with three 8-10 m class telescopes (Keck, VLT, Gemini North) and the Hubble Space Telescope. The moons' mutual orbits were derived individually using a purely Keplerian model. We computed the position of Romulus, the outer moon of the system, at the epoch of a recent stellar occultation which was successfully observed at less than 15 km from our predicted position, within the uncertainty of our model. The occultation data revealed that the Moon, with a surface-area equivalent diameter Ds=23.1$\pm$0.7km, is strongly elongated (axes ratio of 2.7$\pm$0.32.7$\pm$0.3), significantly more than single asteroids of similar size in the main-belt. We concluded that its shape is probably affected by the tides from the primary. A new shape model of the primary was calculated combining adaptive-optics observations with this occultation and 40 archived light-curves recorded since 1978. The difference between the J2=0.024-0.009+0.016 derived from the 3-D shape model assuming an homogeneous distribution of mass for the volume equivalent diameter Dv=273$\pm$10km primary and the null J2 implied by the Keplerian orbits suggests a non-homogeneous mass distribution in the asteroid's interior

SsODNet: The Solar system Open Database Network

The sample of Solar system objects has dramatically increased over the last decade. The amount of measured properties (e.g., diameter, taxonomy, rotation period, thermal inertia) has grown even faster. However, this wealth of information is spread over a myriad of articles, under many different designations per object. We provide a solution to the identification of Solar system objects from any of their multiple names or designations. We also compile and rationalize their properties to provide an easy access to them. We aim to continuously update the database as new measurements become available. We built a Web Service, SsODNet, that offers four access points, each corresponding to an identified necessity in the community: name resolution (quaero), compilation of a large corpus of properties (datacloud), determination of the best estimate among compiled values (ssoCard), and statistical description of the population (ssoBFT). The SsODNet interfaces are fully operational and freely accessible to everyone. The name resolver quaero translates any of the ~5.3 million designations of objects into their current official designation. The datacloud compiles about 105 million parameters (osculating and proper elements, pair and family membership, diameter, albedo, mass, density, rotation period, spin coordinates, phase function parameters, colors, taxonomy, thermal inertia, and Yarkovsky drift) from over 3,000 articles (and growing). For each of the known asteroids and dwarf planets (~1.2 million), a ssoCard providing a single best-estimate for each parameter is available. The SsODNet service provides these resources in a fraction of second upon query. Finally, the large ssoBFT table compiles all the best-estimates in a single table for population-wide studies

Near-infrared spatially resolved spectroscopy of (136108) Haumea's multiple system

The transneptunian region of the solar system is populated by a wide variety of icy bodies showing great diversity. The dwarf planet (136108) Haumea is among the largest TNOs and displays a highly elongated shape and hosts two moons, covered with crystalline water ice like Hamuea. Haumea is also the largest member of the sole TNO family known to date. A catastrophic collision is likely responsible for its unique characteristics. We report here on the analysis of a new set of observations of Haumea obtained with SINFONI at the ESO VLT. Combined with previous data, and using light-curve measurements in the optical and far infrared, we carry out a rotationally resolved spectroscopic study of the surface of Haumea. We describe the physical characteristics of the crystalline water ice present on the surface of Haumea for both regions, in and out of the Dark Red Spot (DRS), and analyze the differences obtained for each individual spectrum. The presence of crystalline water ice is confirmed over more than half of the surface of Haumea. Our measurements of the average spectral slope confirm the redder characteristic of the spot region. Detailed analysis of the crystalline water-ice absorption bands do not show significant differences between the DRS and the remaining part of the surface. We also present the results of applying Hapke modeling to our data set. The best spectral fit is obtained with a mixture of crystalline water ice (grain sizes smaller than 60 micron) with a few percent of amorphous carbon. Improvements to the fit are obtained by adding ~10% of amorphous water ice. Additionally, we used the IFU-reconstructed images to measure the relative astrometric position of the largest satellite Hi`iaka and determine its orbital elements. An orbital solution was computed with our genetic-based algorithm GENOID and our results are in full agreement with recent results.Comment: Accepted for publication in A&

The triaxial ellipsoid dimensions, rotational pole, and bulk density of ESA Rosetta target asteroid (21) Lutetia

We seek the best size estimates of the asteroid (21) Lutetia, the direction of its spin axis, and its bulk density, assuming its shape is well described by a smooth featureless triaxial ellipsoid, and to evaluate the deviations from this assumption. Methods. We derive these quantities from the outlines of the asteroid in 307 images of its resolved apparent disk obtained with adaptive optics (AO) at Keck II and VLT, and combine these with recent mass determinations to estimate a bulk density. Our best triaxial ellipsoid diameters for Lutetia, based on our AO images alone, are a x b x c = 132 x 101 x 93 km, with uncertainties of 4 x 3 x 13 km including estimated systematics, with a rotational pole within 5 deg. of ECJ2000 [long,lat] = [45, -7], or EQJ2000 [RA, DEC] = [44, +9]. The AO model fit itself has internal precisions of 1 x 1 x 8 km, but it is evident, both from this model derived from limited viewing aspects and the radius vector model given in a companion paper, that Lutetia has significant departures from an idealized ellipsoid. In particular, the long axis may be overestimated from the AO images alone by about 10 km. Therefore, we combine the best aspects of the radius vector and ellipsoid model into a hybrid ellipsoid model, as our final result, of 124 +/- 5 x 101 +/- 4 x 93 +/- 13 km that can be used to estimate volumes, sizes, and projected areas. The adopted pole position is within 5 deg. of [long, lat] = [52, -6] or[RA DEC] = [52, +12]. Using two separately determined masses and the volume of our hybrid model, we estimate a density of 3.5 +/- 1.1 or 4.3 +/- 0.8 g cm-3 . From the density evidence alone, we argue that this favors an enstatite-chondrite composition, although other compositions are formally allowed at the extremes (low-porosity CV/CO carbonaceous chondrite or high-porosity metallic). We discuss this in the context of other evidence.Comment: 9 pages, 8 figures, 5 tables, submitted to Astronomy and Astrophysic

The Resolved Asteroid Program - Size, shape, and pole of (52) Europa

With the adaptive optics (AO) system on the 10 m Keck-II telescope, we acquired a high quality set of 84 images at 14 epochs of asteroid (52) Europa on 2005 January 20. The epochs covered its rotation period and, by following its changing shape and orientation on the plane of sky, we obtained its triaxial ellipsoid dimensions and spin pole location. An independent determination from images at three epochs obtained in 2007 is in good agreement with these results. By combining these two data sets, along with a single epoch data set obtained in 2003, we have derived a global fit for (52) Europa of diameters (379x330x249) +/- (16x8x10) km, yielding a volume-equivalent spherical-diameter of 315 +/- 7 km, and a rotational pole within 7 deg of [RA; Dec] = [257,+12] in an Equatorial J2000 reference frame (ECJ2000: 255,+35). Using the average of all mass determinations available forEuropa, we derive a density of 1.5 +/- 0.4, typical of C-type asteroids. Comparing our images with the shape model of Michalowski et al. (A&A 416, 2004), derived from optical lightcurves, illustrates excellent agreement, although several edge features visible in the images are not rendered by the model. We therefore derived a complete 3-D description of Europa's shape using the KOALA algorithm by combining our imaging epochs with 4 stellar occultations and 49 lightcurves. We use this 3-D shape model to assess these departures from ellipsoidal shape. Flat facets (possible giant craters) appear to be less distinct on (52) Europa than on other C-types that have been imaged in detail. We show that fewer giant craters, or smaller craters, is consistent with its expected impact history. Overall, asteroid (52) Europa is still well modeled as a smooth triaxial ellipsoid with dimensions constrained by observations obtained over several apparitions.Comment: Accepted for publication in Icaru

Comparison of machine learning algorithms used to classify the asteroids observed by all-sky surveys

Context. Multifilter photometry from large sky surveys is commonly used to assign asteroid taxonomic types and study various problems in planetary science. To maximize the science output of those surveys, it is important to use methods that best link the spectro-photometric measurements to asteroid taxonomy. Aims. We aim to determine which machine learning methods are the most suitable for the taxonomic classification for various sky surveys. Methods. We utilized five machine learning supervised classifiers: logistic regression, naive Bayes, support vector machines (SVMs), gradient boosting, and MultiLayer Perceptrons (MLPs). Those methods were found to reproduce the Bus-DeMeo taxonomy at various rates depending on the set of filters used by each survey. We report several evaluation metrics for a comprehensive comparison (prediction accuracy, balanced accuracy, F1 score, and the Matthews correlation coefficient) for 11 surveys and space missions. Results. Among the methods analyzed, multilayer perception and gradient boosting achieved the highest accuracy and naive Bayes achieved the lowest accuracy in taxonomic prediction across all surveys. We found that selecting the right machine learning algorithm can improve the success rate by a factor of >2. The best balanced accuracy (similar to 85% for a taxonomic type prediction) was found for the Visible and Infrared Survey telescope for Astronomy (VISTA) and the ESA Euclid mission surveys where broadband filters best map the 1 mu m and 2 mu m olivine and pyroxene absorption bands. Conclusions. To achieve the highest accuracy in the taxonomic type prediction based on multifilter photometric measurements, we recommend the use of gradient boosting and MLP optimized for each survey. This can improve the overall success rate even when compared with naive Bayes. A merger of different datasets can further boost the prediction accuracy. For the combination of the Legacy Survey of Space and Time and VISTA survey, we achieved 90% for the taxonomic type prediction.Peer reviewe

Origin of volatiles in the main belt

We propose a scenario for the formation of the main belt in which asteroids incorporated icy particles formed in the outer solar nebula. We calculate the composition of icy planetesimals formed beyond a heliocentric distance of 5 au in the nebula by assuming that the abundances of all elements, in particular that of oxygen, are solar. As a result, we show that ices formed in the outer solar nebula are composed of a mix of clathrate hydrates, hydrates formed above 50 K and pure condensates produced at lower temperatures. We then consider the inward migration of solids initially produced in the outer solar nebula and show that a significant fraction may have drifted to the current position of the main belt without encountering temperature and pressure conditions high enough to vaporize the ices they contain. We propose that, through the detection and identification of initially buried ices revealed by recent impacts on the surfaces of asteroids, it could be possible to infer the thermodynamic conditions that were present within the solar nebula during the accretion of these bodies, and during the inward migration of icy planetesimals. We also investigate the potential influence that the incorporation of ices in asteroids may have on their porosities and densities. In particular, we show how the presence of ices reduces the value of the bulk density of a given body, and consequently modifies its macroporosity from that which would be expected from a given taxonomic typ