We present a review of the problem of asteroid shape and spin reconstruction from generalised projections; i.e., from lightcurves, disk-resolved images, occultation silhouettes, radar range-Doppler data, and interferometry. The aim of this text is to summarize all important mathematical facts and proofs related to this inverse problem, to describe their implications to observers and modellers, and to provide the reader with all relevant references.

Time-resolved photometry of asteroids can be used for shape and spin reconstruction. If the number of measurements per asteroid is not sufficient to create a model, the whole data set can be used to reconstruct the distribution of shape elongations and pole latitudes in the population. This is done by reconstructing amplitudes of lightcurves that are estimated from dispersion of points observed at (assumed) constant aspect angle. Here, we formulate orbit-averaged approach where the observable is the orbit-averaged dispersion of brightness.

With an estimated diameter in the 320 to 350 km range, (704) Interamnia is the fifth largest main belt asteroid and one of the few bodies that fills the gap in size between the four largest bodies with $D$ > 400 km (Ceres, Vesta, Pallas and Hygiea) and the numerous smaller bodies with $D$ $\lesssim$ 200 km. However, despite its large size, little is known about the shape and spin state of Interamnia and, therefore, about its bulk composition and past collisional evolution. We aimed to test at what size and mass the shape of a small body departs from a nearly ellipsoidal equilibrium shape (as observed in the case of the four largest asteroids) to an irregular shape as routinely observed in the case of smaller ($D$ $\lesssim$ 200 km) bodies. We observed Interamnia as part of our ESO VLT/SPHERE large program (ID: 199.C-0074) at thirteen different epochs. In addition, several new optical lightcurves were recorded. These data, along with stellar occultation data from the literature, were fed to the All-Data Asteroid Modeling (ADAM) algorithm to reconstruct the 3D-shape model of Interamnia and to determine its spin state. Interamnia's volume-equivalent diameter of 332 $\pm$ 6 km implies a bulk density of $\rho$=1.98 $\pm$ 0.68 gcm$^{-3}$ , which suggests that Interamnia - like Ceres and Hygiea - contains a high fraction of water ice, consistent with the paucity of apparent craters. Our observations reveal a shape that can be well approximated by an ellipsoid, and that is compatible with a fluid hydrostatic equilibrium at the 2 $\sigma$ level. The rather regular shape of Interamnia implies that the size and mass limit, under which the shapes of minor bodies with a high amount of water ice in the subsurface become irregular, has to be searched among smaller ($D$ $\lesssim$ 300km) less massive ($m$ $\lesssim$ 3x10$^{19}$ kg) bodies.

Asteroid (7) Iris is an ideal target for disk-resolved imaging owing to its brightness (V$\sim$7-8) and large angular size of 0.33 arcsec during its apparitions. Iris is believed to belong to the category of large unfragmented asteroids that avoided internal differentiation, implying that its current shape and topography may record the first few 100 Myr of the solar system's collisional evolution. We recovered information about the shape and surface topography of Iris from disk-resolved VLT/SPHERE/ZIMPOL images acquired in the frame of our ESO large program. We used the All-Data Asteroid Modeling (ADAM) shape reconstruction algorithm to model the 3D shape of Iris, using optical disk-integrated data and disk-resolved images from SPHERE as inputs. We analyzed the SPHERE images to infer the asteroid's global shape and the morphology of its main craters. We present the volume-equivalent diameter D$_{{\rm eq}}$=214$\pm$5 km, and bulk density $\rho$=2.7$\pm$0.3 g cm$^{-3}$ of Iris. Its shape appears to be consistent with that of an oblate spheroid with a large equatorial excavation. We identified eight putative surface features 20--40 km in diameter detected at several epochs, which we interpret as impact craters. Craters on Iris have depth-to-diameter ratios that are similar to those of analogous 10 km craters on Vesta. The bulk density of Iris is consistent with that of its meteoritic analog, namely LL ordinary chondrites. Considering the absence of a collisional family related to Iris and the number of large craters on its surface, we suggest that its equatorial depression may be the remnant of an ancient (at least 3 Gyr) impact. Iris's shape further opens the possibility that large planetesimals formed as almost perfect oblate spheroids. Finally, we attribute the difference in crater morphology between Iris and Vesta to their different surface gravities.

Context. Over the past decades, several interplanetary missions have studied small bodies in situ, leading to major advances in our understanding of their geological and geophysical properties. These missions, however, have had a limited number of targets. Among them, the NASA Dawn mission has characterised in detail the topography and albedo variegation across the surface of asteroid (4) Vesta down to a spatial resolution of approximately 20 m/pixel scale. Aims. Here our aim was to determine how much topographic and albedo information can be retrieved from the ground with VLT/SPHERE in the case of Vesta, having a former space mission (Dawn) providing us with the ground truth that can be used as a benchmark. Methods. We observed Vesta with VLT/SPHERE/ZIMPOL as part of our ESO large programme (ID 199.C-0074) at six different epochs, and deconvolved the collected images with a parametric point spread function (PSF). We then compared our images with synthetic views of Vesta generated from the 3D shape model of the Dawn mission, on which we projected Vesta's albedo information. Results. We show that the deconvolution of the VLT/SPHERE images with a parametric PSF allows the retrieval of the main topographic and albedo features present across the surface of Vesta down to a spatial resolution of approximately 20-30 km. Contour extraction shows an accuracy of approximately 1 pixel (3.6 mas). The present study provides the very first quantitative estimate of the accuracy of ground-based adaptive-optics imaging observations of asteroid surfaces. Conclusions. In the case of Vesta, the upcoming generation of 30-40m telescopes (ELT, TMT, GMT) should in principle be able to resolve all of the main features present across its surface, including the troughs and the north-south crater dichotomy, provided that they operate at the diffraction limit.

Context. CM-like asteroids (Ch and Cgh classes) are a major population within the broader C-complex, encompassing about 10% of the mass of the main asteroid belt. Their internal structure has been predicted to be homogeneous, based on their compositional similarity as inferred from spectroscopy (Vernazza et al., 2016, AJ 152, 154) and numerical modeling of their early thermal evolution (Bland & Travis, 2017, Sci. Adv. 3, e1602514). Aims. Here we aim to test this hypothesis by deriving the density of the CM-like asteroid (41) Daphne from detailed modeling of its shape and the orbit of its small satellite. Methods. We observed Daphne and its satellite within our imaging survey with the Very Large Telescope extreme adaptive-optics SPHERE/ZIMPOL camera (ID 199.C-0074, PI P. Vernazza) and complemented this data set with earlier Keck/NIRC2 and VLT/NACO observations. We analyzed the dynamics of the satellite with our Genoid meta-heuristic algorithm. Combining our high-angular resolution images with optical lightcurves and stellar occultations, we determine the spin period, orientation, and 3-D shape, using our ADAM shape modeling algorithm. Results. The satellite orbits Daphne on an equatorial, quasi-circular, prograde orbit, like the satellites of many other large main-belt asteroids. The shape model of Daphne reveals several large flat areas that could be large impact craters. The mass determined from this orbit combined with the volume computed from the shape model implies a density for Daphne of 1.77+/-0.26 g/cm3 (3 {\sigma}). This density is consistent with a primordial CM-like homogeneous internal structure with some level of macroporosity (~17%). Conclusions. Based on our analysis of the density of Daphne and 75 other Ch/Cgh-type asteroids gathered from the literature, we conclude that the primordial internal structure of the CM parent bodies was homogeneous.

Asteroid (16) Psyche is the target of the NASA Psyche mission. It is considered one of the few main-belt bodies that could be an exposed proto-planetary metallic core and that would thus be related to iron meteorites. Such an association is however challenged by both its near- and mid-infrared spectral properties and the reported estimates of its density. Here, we aim to refine the density of (16) Psyche to set further constraints on its bulk composition and determine its potential meteoritic analog. We observed (16) Psyche with ESO VLT/SPHERE/ZIMPOL as part of our large program (ID 199.C-0074). We used the high angular resolution of these observations to refine Psyche's three-dimensional (3D) shape model and subsequently its density when combined with the most recent mass estimates. In addition, we searched for potential companions around the asteroid. We derived a bulk density of 3.99\,$\pm$\,0.26\,g$\cdot$cm$^{-3}$ for Psyche. While such density is incompatible at the 3-sigma level with any iron meteorites ($\sim$7.8\,g$\cdot$cm$^{-3}$), it appears fully consistent with that of stony-iron meteorites such as mesosiderites (density $\sim$4.25\,$\cdot$cm$^{-3}$). In addition, we found no satellite in our images and set an upper limit on the diameter of any non-detected satellite of 1460\,$\pm$\,200}\,m at 150\,km from Psyche (0.2\%\,$\times$\,R$_{Hill}$, the Hill radius) and 800\,$\pm$\,200\,m at 2,000\,km (3\%\,$\times$\,$R_{Hill}$). Considering that the visible and near-infrared spectral properties of mesosiderites are similar to those of Psyche, there is merit to a long-published initial hypothesis that Psyche could be a plausible candidate parent body for mesosiderites.

This study aims at advancing mathematical and computational techniques for reconstructing the interior structure of a near earth asteroid (NEA) via Computed Radar Tomography (CRT). We introduce a far-field model for full-wave CRT and validate it numerically for sparse and limited-angle data and an orbiting distance of about 5 km. As a potential future application of the proposed method, we consider the Deep Interior Scanning CUbeSat (DISCUS) concept in which two small spacecraft carrying a bistatic radar localize internal macropores inside a rubble pile asteroid.

Many asteroid databases with lightcurve brightness measurements (e.g. WISE, Pan-STARRS1) contain enormous amounts of data for asteroid shape and spin modelling. While lightcurve inversion is not plausible for individual targets with scarce data, it is possible for large populations with thousands of asteroids, where the distributions of the shape and spin characteristics of the populations are obtainable. We aim to introduce a software implementation of a method that computes the joint shape elongation p and spin latitude beta distributions for a population, with the brightness observations given in an asteroid database. Other main goals are to include a method for performing validity checks of the algorithm, and a tool for a statistical comparison of populations. The LEADER software package read the brightness measurement data for a user-defined subpopulation from a given database. The observations were used to compute estimates of the brightness variations of the population members. A cumulative distribution function (CDF) was constructed of these estimates. A superposition of known analytical basis functions yielded this CDF as a function of the (shape, spin) distribution. The joint distribution can be reconstructed by solving a linear constrained inverse problem. To test the validity of the method, the algorithm can be run with synthetic asteroid models, where the shape and spin characteristics are known, and by using the geometries taken from the examined database. LEADER is a fast and robust software package for solving shape and spin distributions for large populations. There are major differences in the quality and coverage of measurements depending on the database used, so synthetic simulations are always necessary before a database can be reliably used. We show examples of differences in the results when switching to another database.

Context. A lot of photometric data is produced by surveys such as Pan-STARRS, LONEOS, WISE or Catalina. These data are a rich source of information about the physical properties of asteroids. There are several possible approaches for utilizing these data. Lightcurve inversion is a typical method that works with individual asteroids. Our approach in this paper is statistical when we focused on large groups of asteroids like dynamical families and taxonomic classes, and the data were not sufficient for individual models. Aims. Our aim was to study the distributions of shape elongation $b/a$ and the spin axis latitude $\beta$ for various subpopulations of asteroids and to compare our results, based on Pan-STARRS1 survey, with statistics previously done using different photometric data (Lowell database, WISE data). Methods. We use the LEADER algorithm to compare the $b/a$ and $\beta$ distributions for different subpopulations of asteroids. The algorithm creates a cumulative distributive function (CDF) of observed brightness variations, and computes the $b/a$ and $\beta$ distributions using analytical basis functions that yield the observed CDF. A variant of LEADER is used to solve the joint distributions for synthetic populations to test the validity of the method. Results. When comparing distributions of shape elongation for groups of asteroids with different diameters $D$, we found that there are no differences for $D < 25$ km. We also constructed distributions for asteroids with different rotation periods and revealed that the fastest rotators with $P = 0 - 4$ h are more spheroidal than the population with $P = 4 - 8$ h.

We present shape models and volume estimates of twenty asteroids based on relative photometry and adaptive optics images. We discuss error estimation and the effects of myopic deconvolution on shape solutions. For further analysis of the information capacities of data sources, we also present and discuss ambiguity and uniqueness results for the reconstruction of nonconvex shapes from photometry.

Context. Many databases on asteroid brightnesses (e.g. ALCDEF, WISE) are potential sources for extensive asteroid shape and spin modelling. Individual lightcurve inversion models require several apparitions and hundreds of data points per target. However, we can analyse the coarse shape and spin distributions over populations of at least thousands of targets even if there are only a few points and one apparition per asteroid. This is done by examining the distribution of the brightness variations observed within the chosen population. Aims. Brightness variation has been proposed as a population-scale rather than individual-target observable in two studies so far. We aim to examine this approach rigorously to establish its theoretical validity, degree of ill-posedness, and practical applicability. Methods. We model the observed brightness variation of a target population by considering its cumulative distribution function (CDF) caused by the joint distribution function of two fundamental shape and spin indicators. These are the shape elongation and the spin latitude of a simple ellipsoidal model. The main advantage of the model is that we can derive analytical basis functions that yield the observed CDF as a function of the shape and spin distribution. The inverse problem can be treated linearly. Even though the inaccuracy of the model is considerable, databases of thousands of targets should yield some information on the distribution. Results. We establish the theoretical soundness and the typical accuracy limits of the approach both analytically and numerically. Using simulations, we derive a practical estimate of the model distribution in the (shape, spin)-plane. We show that databases such as Wide-field Infrared Survey Explorer (WISE) yield coarse but robust estimates of this distribution, and as an example compare various asteroid families with each other.

Disk-integrated photometric data of asteroids do not contain accurate information on shape details or size scale. Additional data such as disk-resolved images or stellar occultation measurements further constrain asteroid shapes and allow size estimates. We aim to use all available disk-resolved images of about forty asteroids obtained by the Near-InfraRed Camera (Nirc2) mounted on the W.M. Keck II telescope together with the disk-integrated photometry and stellar occultation measurements to determine their volumes. We can then use the volume, in combination with the known mass, to derive the bulk density. We download and process all asteroid disk-resolved images obtained by the Nirc2 that are available in the Keck Observatory Archive (KOA). We combine optical disk-integrated data and stellar occultation profiles with the disk-resolved images and use the All-Data Asteroid Modeling (ADAM) algorithm for the shape and size modeling. Our approach provides constraints on the expected uncertainty in the volume and size as well. We present shape models and volume for 41 asteroids. For 35 asteroids, the knowledge of their mass estimates from the literature allowed us to derive their bulk densities. We clearly see a trend of lower bulk densities for primitive objects (C-complex) than for S-complex asteroids. The range of densities in the X-complex is large, suggesting various compositions. Moreover, we identified a few objects with rather peculiar bulk densities, which is likely a hint of their poor mass estimates. Asteroid masses determined from the Gaia astrometric observations should further refine most of the density estimates.

Asteroid (130) Elektra belongs to one of the six known triple asteroids in the main belt, so its mass has been reliably determined. We aim to use all available disk-resolved images of (130) Elektra obtained by the SPHERE instrument at VLT and by the Nirc2 of the Keck telescope together with the disk-integrated photometry to determine its shape model and its size. The volume can be then used in combination with the known mass to derive the bulk density of the primary. We apply the All-Data Asteroid Modeling (ADAM) algorithm to the optical disk-integrated data, 2 disk-resolved images obtained by the SPHERE instrument and 13 disk-resolved images from the Nirc2 of the Keck telescope, and derive the shape model and size of Elektra. We present the shape model, volume-equivalent diameter (199$\pm$7 km) and bulk density (1.60$\pm$0.13 g cm$^{-3}$) of the C-type asteroid Elektra.

Context: Large all-sky surveys provide us with a lot of photometric data that are sparse in time (typically few measurements per night) and can be potentially used for the determination of shapes and rotational states of asteroids. The method which is generally used to derive these parameters is the lightcurve inversion. However, for most asteroids their sparse data are not accurate enough to derive a unique model and the lightcurve inversion method is thus not very efficient. Aims: To fully utilize photometry sparse in time, we developed a new simplified model and applied it on the data from the Lowell photometric database. Our aim was to derive spin axis orientations and shape elongations of asteroids and to reconstruct distributions of these parameters for selected subpopulations to find if there are some differences. Methods: We model asteroids as geometrically scattering triaxial ellipsoids. Observed values of mean brightness and the dispersion of brightness are compared with computed values obtained from the parameters of the model -- ecliptical longitude $\lambda$ and latitude $\beta$ of the pole and the ratios $a/b$, $b/c$ of axes of the ellipsoid. These parameters are optimized to get the best agreement with the observation. Results: We found that the distribution of $\lambda$ for main-belt asteroids is not uniform (in agreement with findings of Bowell et al., 2014, M&PS, 49, 95) and is dependent on the inclination of orbit. Surprisingly, the non-uniformity of $\lambda$ distribution is larger for asteroids residing on low-inclination orbits. We also studied distributions of $a/b$ for several groups of asteroids and found that small asteroids ($D<25\,$km) are on average more elongated than large ones.

This paper concerns electromagnetic 3D subsurface imaging in connection with sparsity of signal sources. We explored an imaging approach that can be implemented in situations that allow obtaining a large amount of data over a surface or a set of orbits but at the same time require sparsity of the signal sources. Characteristic to such a tomography scenario is that it necessitates the inversion technique to be genuinely three-dimensional: For example, slicing is not possible due to the low number of sources. Here, we primarily focused on astrophysical subsurface exploration purposes. As an example target of our numerical experiments we used a synthetic small planetary object containing three inclusions, e.g. voids, of the size of the wavelength. A tetrahedral arrangement of source positions was used, it being the simplest symmetric point configuration in 3D. Our results suggest that somewhat reliable inversion results can be produced within the present a priori assumptions, if the data can be recorded at a specific resolution. This is valuable early-stage knowledge especially for design of future planetary missions in which the payload needs to be minimized, and potentially also for the development of other lightweight subsurface inspection systems.

The goal of this paper is to advance mathematical and computational methodology for orbiter-to-orbiter radio tomography of small planetary objects. In this study, an advanced full waveform forward model is coupled with a total variation based inversion technique. We use a satellite formation model in which a single unit receives a signal that is transmitted by one or more transponder satellites. Numerical results for a 2D domain are presented.

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.

We introduce ADAM, the All-Data Asteroid Modelling algorithm. ADAM is simple and universal since it handles all disk-resolved data types (adaptive optics or other images, interferometry, and range-Doppler radar data) in a uniform manner via the 2D Fourier transform, enabling fast convergence in model optimization. The resolved data can be combined with disk-integrated data (photometry). In the reconstruction process, the difference between each data type is only a few code lines defining the particular generalized projection from 3D onto a 2D image plane. Occultation timings can be included as sparse silhouettes, and thermal infrared data are efficiently handled with an approximate algorithm that is sufficient in practice due to the dominance of the high-contrast (boundary) pixels over the low-contrast (interior) ones. This is of particular importance to the raw ALMA data that can be directly handled by ADAM without having to construct the standard image. We study the reliability of the inversion by using the independent shape supports of function series and control-point surfaces. When other data are lacking, one can carry out fast nonconvex lightcurve-only inversion, but any shape models resulting from it should only be taken as illustrative global-scale ones.

This study concerned conebeam travel-time tomography. The focus was on a sparse distribution of signal sources that can be necessary in a challenging in situ environment such as in asteroid tomography. The goal was to approximate the minimum number of source positions needed for robust detection of refractive anomalies, e.g., voids within an asteroid or a casting defects in concrete. Experimental ultrasonic data were recorded utilizing as a target a 150 mm plastic cast cube containing three stones with diameter between 22 and 41 mm. A signal frequency of 55 kHz (35 mm wavelength) was used. Source counts from one to six were tested for different placements. Based on our statistical inversion approach and analysis of the results, three or four sources were found to lead to reliable inversion. The source configurations investigated were also ranked according to their performance. Our results can be used, for example, in the planning of planetary missions as well as in material testing.

We present a simple choice of integration variables that can be used to exploit the near-integrable character of problems in celestial mechanics. The approach is based on the well-known principle of variation of parameters: instead of orbital elements, we use the phase-space coordinates the object would have at a given point in its (Keplerian) orbit if the perturbing forces were removed. This formulation is suitable for almost any numerical integrator; thus, multistep schemes are easy to build, stepsize can be adjusted, and dissipative forces are allowed. Compared with traditional non-symplectic N-body integrators, the approach often offers increase in speed or accuracy if perturbations are small.

We seek to characterize giant-planet systems by their gravitational scattering properties. We do this to a given system by integrating it numerically along with a large number of hypothetical small bodies that are initially in eccentric habitable zone (HZ)-crossing orbits. Our analysis produces a single number, the escape rate, which represents the rate at which the small-body flux is perturbed away by the giant planets into orbits that no longer pose a threat to terrestrial planets inside the HZ. Obtaining the escape rate this way is similar to computing the largest Liapunov exponent as the exponential rate of divergence of two nearby orbits. For a terrestrial planet inside the HZ, the escape rate value quantifies the "protective" effect that the studied giant-planet system offers. Therefore, escape rates could provide information on whether certain giant-planet configurations produce a more desirable environment for life than the others. We present some computed escape rates on selected planetary systems, focusing on effects of varying the masses and semi-major axes of the giant planets. In the case of our Solar System we find rather surprisingly that Jupiter, in its current orbit, may provide a minimal amount of protection to the Earth.

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.

The spin state of small asteroids can change on a long timescale by the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect, the net torque that arises from anisotropically scattered sunlight and proper thermal radiation from an irregularly-shaped asteroid. The secular change in the rotation period caused by the YORP effect can be detected by analysis of asteroid photometric lightcurves. We analyzed photometric lightcurves of near-Earth asteroids (1865) Cerberus, (2100) Ra-Shalom, and (3103) Eger with the aim to detect possible deviations from the constant rotation caused by the YORP effect. We carried out new photometric observations of the three asteroids, combined the new lightcurves with archived data, and used the lightcurve inversion method to model the asteroid shape, pole direction, and rotation rate. The YORP effect was modeled as a linear change in the rotation rate in time d\omega /dt. Values of d\omega/ dt derived from observations were compared with the values predicted by theory. We derived physical models for all three asteroids. We had to model Eger as a nonconvex body because the convex model failed to fit the lightcurves observed at high phase angles. We probably detected the acceleration of the rotation rate of Eger d\omega / dt = (1.4 +/- 0.6) x 10^{-8} rad/d (3\sigma error), which corresponds to a decrease in the rotation period by 4.2 ms/yr. The photometry of Cerberus and Ra-Shalom was consistent with a constant-period model, and no secular change in the spin rate was detected. We could only constrain maximum values of |d\omega / dt| < 8 x 10^{-9} rad/d for Cerberus, and |d\omega / dt| < 3 x 10^{-8} rad/d for Ra-Shalom.

We compare various approaches to find the most efficient method for the practical computation of the lightcurves (integrated brightnesses) of irregularly shaped bodies such as asteroids at arbitrary viewing and illumination geometries. For convex models, this reduces to the problem of the numerical computation of an integral over a simply defined part of the unit sphere. We introduce a fast method, based on Lebedev quadratures, which is optimal for both lightcurve simulation and inversion in the sense that it is the simplest and fastest widely applicable procedure for accuracy levels corresponding to typical data noise. The method requires no tessellation of the surface into a polyhedral approximation. At the accuracy level of 0.01 mag, it is up to an order of magnitude faster than polyhedral sums that are usually applied to this problem, and even faster at higher accuracies. This approach can also be used in other similar cases that can be modelled on the unit sphere. The method is easily implemented in lightcurve inversion by a simple alteration of the standard algorithm/software.

We present a comparison of our results from ground-based observations of asteroid (21) Lutetia with imaging data acquired during the flyby of the asteroid by the ESA Rosetta mission. This flyby provided a unique opportunity to evaluate and calibrate our method of determination of size, 3-D shape, and spin of an asteroid from ground-based observations. We present our 3-D shape-modeling technique KOALA which is based on multi-dataset inversion. We compare the results we obtained with KOALA, prior to the flyby, on asteroid (21) Lutetia with the high-spatial resolution images of the asteroid taken with the OSIRIS camera on-board the ESA Rosetta spacecraft, during its encounter with Lutetia. The spin axis determined with KOALA was found to be accurate to within two degrees, while the KOALA diameter determinations were within 2% of the Rosetta-derived values. The 3-D shape of the KOALA model is also confirmed by the spectacular visual agreement between both 3-D shape models (KOALA pre- and OSIRIS post-flyby). We found a typical deviation of only 2 km at local scales between the profiles from KOALA predictions and OSIRIS images, resulting in a volume uncertainty provided by KOALA better than 10%. Radiometric techniques for the interpretation of thermal infrared data also benefit greatly from the KOALA shape model: the absolute size and geometric albedo can be derived with high accuracy, and thermal properties, for example the thermal inertia, can be determined unambiguously. We consider this to be a validation of the KOALA method. Because space exploration will remain limited to only a few objects, KOALA stands as a powerful technique to study a much larger set of small bodies using Earth-based observations.

Asteroid sizes can be directly measured by observing occultations of stars by asteroids. When there are enough observations across the path of the shadow, the asteroid's projected silhouette can be reconstructed. Asteroid shape models derived from photometry by the lightcurve inversion method enable us to predict the orientation of an asteroid for the time of occultation. By scaling the shape model to fit the occultation chords, we can determine the asteroid size with a relative accuracy of typically ~ 10%. We combine shape and spin state models of 44 asteroids (14 of them are new or updated models) with the available occultation data to derive asteroid effective diameters. In many cases, occultations allow us to reject one of two possible pole solutions that were derived from photometry. We show that by combining results obtained from lightcurve inversion with occultation timings, we can obtain unique physical models of asteroids.

Aims. We determine the physical properties (spin state and shape) of asteroid (21) Lutetia, target of the ESA Rosetta mission, to help in preparing for observations during the flyby on 2010 July 10 by predicting the orientation of Lutetia as seen from Rosetta. Methods. We use our novel KOALA inversion algorithm to determine the physical properties of asteroids from a combination of optical lightcurves, disk-resolved images, and stellar occultations, although the latter are not available for (21) Lutetia. Results. We find the spin axis of (21) Lutetia to lie within 5 degrees of ({\lambda} = 52 deg., {\beta} = -6 deg.) in Ecliptic J2000 reference frame (equatorial {\alpha} = 52 deg., {\delta} = +12 deg.), and determine an improved sidereal period of 8.168 270 \pm 0.000 001 h. This pole solution implies the southern hemisphere of Lutetia will be in "seasonal" shadow at the time of the flyby. The apparent cross-section of Lutetia is triangular as seen "pole-on" and more rectangular as seen "equator-on". The best-fit model suggests the presence of several concavities. The largest of these is close to the north pole and may be associated with large impacts.

We show that the boundary curves (profiles) in $\R^2$ of the generalized projections of a body in $\R^3$ uniquely determine a large class of shapes, and that sparse profile data, combined with projection volume (brightness) data, can be used to reconstruct the shape and the spin state of a body. We also present an optimal strategy for the relative weighting of the data modes in the inverse problem, and derive the maximum compatibility estimate (MCE) that corresponds to the maximum likelihood or maximum a posteriori estimates in the case of a single data mode. MCE is not explicitly dependent on the noise levels, scale factors or numbers of data points of the complementary data modes, and can be determined without the mode weight parameters. We present a solution method well suitable for adaptive optics images in particular, and discuss various choices of regularization functions.

We acquired and analyzed adaptive-optics imaging observations of asteroid (2) Pallas from Keck II and the Very Large Telescope taken during four Pallas oppositions between 2003 and 2007, with spatial resolution spanning 32-88 km (image scales 13-20 km/pix). We improve our determination of the size, shape, and pole by a novel method that combines our AO data with 51 visual light-curves spanning 34 years of observations as well as occultation data. The shape model of Pallas derived here reproduces well both the projected shape of Pallas on the sky and light-curve behavior at all the epochs considered. We resolved the pole ambiguity and found the spin-vector coordinates to be within 5 deg. of [long, lat] = [30 deg., -16 deg.] in the ECJ2000.0 reference frame, indicating a high obliquity of ~84 deg., leading to high seasonal contrast. The best triaxial-ellipsoid fit returns radii of a=275 km, b= 258 km, and c= 238 km. From the mass of Pallas determined by gravitational perturbation on other minor bodies [(1.2 +/- 0.3) x 10-10 Solar Masses], we derive a density of 3.4 +/- 0.9 g.cm-3 significantly different from the density of C-type (1) Ceres of 2.2 +/- 0.1 g.cm-3. Considering the spectral similarities of Pallas and Ceres at visible and near-infrared wavelengths, this may point to fundamental differences in the interior composition or structure of these two bodies. We define a planetocentric longitude system for Pallas, following IAU guidelines. We also present the first albedo maps of Pallas covering ~80% of the surface in K-band. These maps reveal features with diameters in the 70-180 km range and an albedo contrast of about 6% wrt the mean surface albedo.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We report on the results of a six-month photometric study of the main-belt binary C-type asteroid 121 Hermione, performed during its 2007 opposition. We took advantage of the rare observational opportunity afforded by one of the annual equinoxes of Hermione occurring close to its opposition in June 2007. The equinox provides an edge-on aspect for an Earth-based observer, which is well suited to a thorough study of Hermione's physical characteristics. The catalog of observations carried out with small telescopes is presented in this work, together with new adaptive optics (AO) imaging obtained between 2005 and 2008 with the Yepun 8-m VLT telescope and the 10-m Keck telescope. The most striking result is confirmation that Hermione is a bifurcated and elongated body, as suggested by Marchis et al., (2005). A new effective diameter of 187 +/- 6 km was calculated from the combination of AO, photometric and thermal observations. The new diameter is some 10% smaller than the hitherto accepted radiometric diameter based on IRAS data. The reason for the discrepancy is that IRAS viewed the system almost pole-on. New thermal observations with the Spitzer Space Telescope agree with the diameter derived from AO and lightcurve observations. On the basis of the new AO astrometric observations of the small 32-km diameter satellite we have refined the orbit solution and derived a new value of the bulk density of Hermione of 1.4 +0.5/-0.2 g cm-3. We infer a macroscopic porosity of ~33 +5/-20%.

We obtain full information on the orbital parameters by combining radial velocity and astrometric measurements by means of Bayesian inference. We sample the parameter probability densities of orbital model parameters with a Markov chain Monte Carlo (McMC) method in simulated observational scenarios to test the detectability of planets with orbital periods longer than the observational timelines. We show that, when fitting model parameters simultaneously to measurements from both sources, it is possible to extract much more information from the measurements than when using either source alone. We demonstrate this by studying the orbit of recently found extra-solar planet HD 154345 b.

We study the inverse problem of deducing the dynamical characteristics (such as the potential field) of large systems from kinematic observations. We show that, for a class of steady-state systems, the solution is unique even with fragmentary data, dark matter, or selection (bias) functions. Using spherically symmetric models for simulations, we investigate solution convergence and the roles of data noise and regularization in the inverse problem. We also present a method, analogous to tomography, for comparing the observed data with a model probability distribution such that the latter can be determined.

In 2003, we initiated a long-term adaptive optics campaign to study the orbit of various main-belt asteroidal systems. Here we present a consistent solution for the mutual orbits of four binary systems: 22 Kalliope, 45 Eugenia, 107 Camilla and 762 Pulcova. With the exception of 45 Eugenia, we did not detect any additional satellites around these systems although we have the capability of detecting a loosely-bound fragment (located at 1/4 x RHill) that is ~40 times smaller in diameter than the primary. The common characteristic of these mutual orbits is that they are roughly circular. Three of these binary systems belong to a C-"group" taxonomic class. Our estimates of their bulk densities are consistently lower (~1 g/cm3) than their associated meteorite analogs, suggesting an interior porosity of 30-50% (taking CI-CO meteorites as analogs). 22 Kalliope, a W-type asteroid, has a significantly higher bulk density of ~3 g/cm3, derived based on IRAS radiometric size measurement. We compare the characteristics of these orbits in the light of tidal-effect evolution.

In 2007, the M-type binary asteroid 22 Kalliope reached one of its annual equinoxes. As a consequence, the orbit plane of its small moon, Linus, was aligned closely to the Sun's line of sight, giving rise to a mutual eclipse season. A dedicated international campaign of photometric observations, based on amateur-professional collaboration, was organized and coordinated by the IMCCE in order to catch several of these events. The set of the compiled observations is released in this work. We developed a relevant model of these events, including a topographic shape model of Kalliope refined in the present work, the orbit solution of Linus as well as the photometric effect of the shadow of one component falling on the other. By fitting this model to the only two full recorded events, we derived a new estimation of the equivalent diameter of Kalliope of 166.2+/-2.8km, 8% smaller than its IRAS diameter. As to the diameter of Linus, considered as purely spherical, it is estimated to 28+/-2 km. This substantial "shortening" of Kalliope gives a bulk density of 3.35+/-0.33g/cm3, significantly higher than past determinations but more consistent with its taxonomic type. Some constraints can be inferred on the composition.

We obtained N- and Q-band observations of the Apollo-type asteroid 25143 Itokawa during its close Earth approach in July 2004 with TIMMI2 at the ESO 3.6 m telescope. Our photometric measurement, in combination with already published data, allowed us to derive a radiometric effective diameter of 0.32+/-0.03 km and an albedo of 0.19 +0.11/-0.03 through a thermophysical model. This effective diameter corresponds to a slightly asymmetrical and flattened ellipsoid of the approximate size of 520(+/-50) x 270(+/-30) x 230(+/-20) m, based on the Kaasalainen et al. (2005) shape model. Our studies show that the thermal observations lead to size estimates which are about 15% smaller than the radar results (Ostro et al. 2005), slightly outside the stated radar uncertainties of +/-10%. We determined a rather high thermal inertia of 750 Jm-2s-0.5K-1. This is an indication for a bare rock dominated surface, a thick dust regolith can be excluded as well as a metallic surface. From our data we constructed a 10.0 micrometer thermal lightcurve which is nicely matched in amplitude and phase by the shape and spin vector solution in combination with our TPM description. The assumed S-type bulk density in combination with radiometric size lead to a total mass estimate of 4.5 +2.0/-1.8 x 10^10 kg.