Near-Earth asteroid (410777) 2009 FD is a potentially hazardous asteroid with potential impacts on Earth at the end of the 22nd century. The astrometry collected during the 2019 apparition provides information on the trajectory of (410777) by constraining the Yarkovsky effect, which is the main source of uncertainty for future predictions, and informing the impact hazard assessment. We included the Yarkovsky effect in the force model and estimated its magnitude from the fit to the (410777) optical and radar astrometric data. We performed the (410777) hazard assessment over 200 years by using two independent approaches: the NEODyS group adopted a generalisation of the Line Of Variations method in a 7-dimensional space, whereas the JPL team resorted to the Multilayer Clustered Sampling technique. We obtain a 4-sigma detection of the Yarkovsky effect acting on (410777), which corresponds to a semimajor axis drift of (3.6 +/- 0.9) x 10^{-3} au/Myr. As for the hazard results of both teams, the main impact possibility in 2185 is ruled out and the only remaining one is in 2190, but with a probability smaller than 10^{-8}.

The completeness limit is a key quantity to measure the reliability of an impact monitoring system. It is the impact probability threshold above which every virtual impactor has to be detected. A goal of this paper is to increase the completeness without increasing the computational load. We propose a new method to sample the Line Of Variations (LOV) with respect to the previously one used in NEODyS. The step-size of the sampling is not uniform in the LOV parameter, since the probability of each LOV segment between consecutive points is kept constant. Moreover, the sampling interval has been extended to the larger interval [-5,5] in the LOV parameter and a new decomposition scheme in sub-returns is provided to deal with the problem of duplicated points in the same return. The impact monitoring system CLOMON-2 has been upgraded with all these features, resulting in a decrease of the generic completeness limit IP* by a factor $\simeq 4$ and in an increase of the computational load by a factor $\simeq 2$. Moreover, we statistically investigate the completeness actually reached by the system with two different methods: a direct comparison with the results of the independent system Sentry at JPL and an empirical power-law to model the number of virtual impactors as a function of the impact probability. We found empirically that the number of detected virtual impactors with $IP>IP^*$ appears to grow according to a power-law, proportional to $IP^{-2/3}$. An analytical model explaining this power-law is currently an open problem, but we think it is related to the way the number of virtual impactors within a time $t_{rel}$ from the first observed close approach accumulates, for which we prove the existence of a power-law proportional to $t_{rel}^3$. The power-law allows us to detect a loss of efficiency in the virtual impactors search for impact probabilities near the generic completeness limit.

We continue our investigation of the bulk properties of asteroid dynamical families identified using only asteroid proper elements (Milani et al. 2014) to provide plausible collisional interpretations. We focus on cratering families consisting of a substantial parent body and many small fragments. We propose a quantitative definition of cratering families based on the fraction in volume of the fragments with respect to the parent body; fragmentation families are above this empirical boundary. We assess the compositional homogeneity of the families and their shape in proper element space by computing the differences of the proper elements of the fragments with respect to the ones of the major body, looking for anomalous asymmetries produced either by post-formation dynamical evolution, or by multiple collisional/cratering events, or by a failure of the Hierarchical Clustering Method (HCM) for family identification. We identified a total of 25 dynamical families with more than 100 members ranging from moderate to heavy cratering. For three families (4, 15 and 283) we confirm the occurrence of two separate cratering events. For family 20, we propose a double collision origin, not previously identified. In four cases (31, 480, 163 and 179) we performed a dedicated search for dynamical resonant transport mechanisms that could have substantially changed the shape of the family. By using a new synthetic method for computation of secular frequencies, we found possible solutions for families 31, 480, and 163, but not for family 179, for which we propose a new interpretation, based on a secular resonance contaminating this family: the family of 179 should be split into two separate clusters, one containing (179) itself and the other, family (9506) Telramund, of fragmentation type, for which we have computed an age.

The ESA Euclid mission is a space telescope that will survey ~15,000 square degrees of the sky, primarily to study the distant universe (constraining cosmological parameters through the lensing of galaxies). It is also expected to observe ~150,000 Solar System Objects (SSOs), primarily in poorly understood high inclination populations, as it will mostly avoid +/-15 degrees from the ecliptic plane. With a launch date of 2022 and a 6 year survey, Euclid and LSST will operate at the same time, and have complementary capabilities. We propose a LSST mini-survey to coordinate quasi-simultaneous observations between these two powerful observatories, when possible, with the primary aim of greatly improving the orbits of SSOs discovered by these facilities. As Euclid will operate from a halo orbit around the Sun-Earth L2 Lagrangian point, there will be significant parallax between observations from Earth and Euclid (0.01 AU). This means that simultaneous observations will give an independent distance measurement to SSOs, giving additional constraints on orbits compared to single Euclid visits.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

We present an updated set of near-Earth asteroids with a Yarkovsky-related semimajor axis drift detected from the orbital fit to the astrometry. We find 87 reliable detections after filtering for the signal-to-noise ratio of the Yarkovsky drift estimate and making sure the estimate is compatible with the physical properties of the analyzed object. Furthermore, we find a list of 24 marginally significant detections, for which future astrometry could result in a Yarkovsky detection. A further outcome of the filtering procedure is a list of detections that we consider spurious because unrealistic or not explicable with the Yarkovsky effect. Among the smallest asteroids of our sample, we determined four detections of solar radiation pressure, in addition to the Yarkovsky effect. As the data volume increases in the near future, our goal is to develop methods to generate very long lists of asteroids with reliably detected Yarkovsky effect, with limited amounts of case by case specific adjustments. Furthermore, we discuss the improvements this work could bring to impact monitoring. In particular, we exhibit two asteroids for which the adoption of a non-gravitational model is needed to make reliable impact predictions.

The Gaia spacecraft of the European Space Agency (ESA) has been securing observations of solar system objects (SSOs) since the beginning of its operations. Gaia Data Release 2 (DR2) contains the observations of a selected sample of 14,099 SSOs. These asteroids have been already identified and have been numbered by the Minor Planet Center. Positions are provided for each Gaia observation at CCD level. As additional information, the apparent brightness of SSOs in the unfiltered G band is also provided for selected observations. We explain the processing of SSO data, and describe the criteria we used to select the sample published in Gaia DR2. We then explore the data set to assess its quality. To exploit the epoch astrometry of asteroids in Gaia DR2 it is necessary to take into account the unusual properties of the uncertainty, as the position information is nearly one-dimensional. When this aspect is handled appropriately, an orbit fit can be obtained with post-fit residuals that are overall consistent with the a-priori error model that was used to define individual values of the astrometric uncertainty. The distribution of residuals allowed us to identify possible contaminants in the data set. Photometry in the G band was compared to computed values from reference asteroid shapes and to the flux registered at the corresponding epochs by the red and blue photometers (RP and BP). The overall astrometric performance is close to the expectations, with an optimal range of brightness G~12-17. In this range, the typical transit-level accuracy is well below 1 mas. For fainter asteroids, the growing photon noise deteriorates the performance. Asteroids brighter than G~12 are affected by a lower performance of the processing of their signals. The dramatic improvement brought by Gaia DR2 astrometry of SSOs is demonstrated by preliminary tests on the detection of subtle non-gravitational effects.

The second Gaia data release (DR2), contains very precise astrometric and photometric properties for more than one billion sources, astrophysical parameters for dozens of millions, radial velocities for millions, variability information for half a million of stellar sources and orbits for thousands of solar system objects. Before the Catalogue publication, these data have undergone dedicated validation processes. The goal of this paper is to describe the validation results in terms of completeness, accuracy and precision of the various Gaia DR2 data. The validation processes include a systematic analysis of the Catalogue content to detect anomalies, either individual errors or statistical properties, using statistical analysis, and comparisons to external data or to models. Although the astrometric, photometric and spectroscopic data are of unprecedented quality and quantity, it is shown that the data cannot be used without a dedicated attention to the limitations described here, in the Catalogue documentation and in accompanying papers. A particular emphasis is put on the caveats for the statistical use of the data in scientific exploitation.

Short-arc orbit determination is crucial when an asteroid is first discovered. In these cases usually the observations are so few that the differential correction procedure may not converge. We have developed an initial orbit computation method, based on the systematic ranging, an orbit determination techniques which systematically explores a raster in the topocentric range and range-rate space region inside the admissible region. We obtain a fully rigorous computation of the probability for the asteroid that could impact the Earth within few days from the discovery, without any a priori assumption. We test our method on the two past impactors 2008 TC3 and 2014 AA, on some very well known cases, and on two particular objects observed by the ESA Gaia mission.

The Gaia Data Release 1 (GDR1) is a first, important step on the path of evolution of astrometric accuracy towards a much improved situation. Although asteroids are not present in GDR1, this intermediate release already impacts asteroid astrometry. Our goal is to investigate how the GDR1 can change the approach to a few typical problems, including the determination of orbits from short-arc astrometry, the exploitation of stellar occultations, and the impact risk assessment. We employ optimised asteroid orbit determination tools, and study the resulting orbit accuracy and post-fit residuals. For this goal, we use selected ground-based asteroid astrometry, and occultation events observed in the past. All measurements are calibrated by using GDR1 stars. We show that, by adopting GDR1, very short measurement arcs can already provide interesting orbital solutions, capable of correctly identifying Near Earth Asteroids (NEAs) and providing a much more accurate risk rating. We also demonstrate that occultations, previously used to derive asteroid size and shapes, now reach a new level of accuracy at which they can be fruitfully used to obtain astrometry at the level of accuracy of Gaia star positions.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

In this work we have estimated 10 collisional ages of 9 families for which for different reasons our previous attempts failed. In general, these are difficult cases that required dedicated effort, such as a new family classifications for asteroids in mean motion resonances, as well as a revision of the classification inside the $3/2$ resonance. Of the families locked in mean motion resonances, we succeeded in determining ages of the families of (1911) Schubart and of the "super-Hilda" family, assuming this is actually a severely eroded original family of (153) Hilda. In the Trojan region we found families with almost no Yarkovsky evolution, for which we could compute only physically implausible ages. Hence, we interpreted their modest dispersions of proper eccentricities and inclinations as implying that the Trojan asteroid families are fossil families, frozen at their proper elements determined by the original ejection velocity field. We have found a new family, among the Griquas locked in the 2/1 resonance with Jupiter: (11097) 1994 UD1. We have estimated the ages of 6 families affected by secular resonances: families of (5) Astraea, (25) Phocaea, (283) Emma, (363) Padua, (686) Gersuind, and (945) Barcelona. By using a numerical calibration method, we have shown that the secular resonances do not affect significanly the secular change of proper a. For the family of (145) Adeona we could estimate the age only after removal of a number of assumed interlopers. With the present paper we have concluded the series dedicated to the determination of asteroid ages with a uniform method. We computed the ages for a total of 57 families with $>100$ members. There remain families too small at present to provide reliable estimates, as well as some complex families (221, 135, 298) which may have more ages than we could currently estimate.

Orbit determination is possible for a chaotic orbit of a dynamical system, given a finite set of observations, provided the initial conditions are at the central time. In a simple discrete model, the standard map, we tackle the problem of chaotic orbit determination when observations extend beyond the predictability horizon. If the orbit is hyperbolic, a shadowing orbit is computed by the least squares orbit determination. We test both the convergence of the orbit determination iterative procedure and the behaviour of the uncertainties as a function of the maximum number $n$ of map iterations observed. When the initial conditions belong to a chaotic orbit, the orbit determination is made impossible by numerical instability beyond a computability horizon, which can be approximately predicted by a simple formula. Moreover, the uncertainty of the results is sharply increased if a dynamical parameter is added to the initial conditions as parameter to be estimated. The uncertainty of the dynamical parameter decreases like $n^a$ with $a<0$ but not large (of the order of unity). If only the initial conditions are estimated, their uncertainty decreases exponentially with $n$. If they belong to a non-chaotic orbit the computational horizon is much larger, if it exists at all, and the decrease of the uncertainty is polynomial in all parameters, like $n^a$ with $a\simeq 1/2$. The Shadowing Lemma does not dictate what the asymptotic behaviour of the uncertainties should be. These phenomena have significant implications, which remain to be studied, in practical problems of orbit determination involving chaos, such as the chaotic rotation state of a celestial body and a chaotic orbit of a planet-crossing asteroid undergoing many close approaches.

A new family classification, based on a catalog of proper elements with $\sim 384,000$ numbered asteroids and on new methods is available. For the $45$ dynamical families with $>250$ members identified in this classification, we present an attempt to obtain statistically significant ages: we succeeded in computing ages for $37$ collisional families. We used a rigorous method, including a least squares fit of the two sides of a V-shape plot in the proper semimajor axis, inverse diameter plane to determine the corresponding slopes, an advanced error model for the uncertainties of asteroid diameters, an iterative outlier rejection scheme and quality control. The best available Yarkovsky measurement was used to estimate a calibration of the Yarkovsky effect for each family. The results are presented separately for the families originated in fragmentation or cratering events, for the young, compact families and for the truncated, one-sided families. For all the computed ages the corresponding uncertainties are provided. We found 2 cases where two separate dynamical families form together a single V-shape with compatible slopes, thus indicating a single collisional event. We have also found 3 examples of dynamical families containing multiple collisional families, plus a dubious case. We have found 2 cases of families containing a conspicuous subfamily, such that it is possible to measure the slope of a distinct V-shape, thus the age of the secondary collision. We also provide data on the central gaps appearing in some families. The ages computed in this paper are obtained with a single and uniform methodology, thus the ages of different families can be compared, providing a first example of collisional chronology of the asteroid main belt.

Asteroid 2009 FD could impact Earth between 2185 and 2196. The long term propagation to the possible impacts and the intervening planetary encounters make 2009 FD one of the most challenging asteroids in terms of hazard assessment. To compute accurate impact probabilities we model the Yarkovsky effect by using the available physical characterization of 2009 FD and general properties of the Near Earth Asteroid population. We perform the hazard assessment with two independent methods: the first method is a generalization of the standard impact monitoring algorithms in use by NEODyS and Sentry, while the second one is based on a Monte Carlo approach. Both methods generate orbital samples in a 7 dimensional space that includes orbital elements and the parameter characterizing the Yarkovsky effect. The highest impact probability is $2.7 \times 10^{-3}$ for an impact during the 2185 Earth encounter. Impacts after 2185 corresponding to resonant returns are possible, the most relevant being in 2190 with a probability of $3 \times 10^{-4}$. Both numerical methods can be used in the future to handle similar cases. The structure of resonant returns and the list of the possible keyholes on the Target Plane of the scattering encounter in 2185 can be predicted by an analytic theory.

The target asteroid of the OSIRIS-REx asteroid sample return mission, (101955) Bennu (formerly 1999 RQ$_{36}$), is a half-kilometer near-Earth asteroid with an extraordinarily well constrained orbit. An extensive data set of optical astrometry from 1999--2013 and high-quality radar delay measurements to Bennu in 1999, 2005, and 2011 reveal the action of the Yarkovsky effect, with a mean semimajor axis drift rate $da/dt = (-19.0 \pm 0.1)\times 10^{-4}$ au/Myr or $284\pm 1.5\;\rm{m/yr}$. The accuracy of this result depends critically on the fidelity of the observational and dynamical model. As an example, neglecting the relativistic perturbations of the Earth during close approaches affects the orbit with $3\sigma$ significance in $da/dt$. The orbital deviations from purely gravitational dynamics allow us to deduce the acceleration of the Yarkovsky effect, while the known physical characterization of Bennu allows us to independently model the force due to thermal emissions. The combination of these two analyses yields a bulk density of $\rho = 1260\pm70\,\rm{kg/m^3}$, which indicates a macroporosity in the range $40\pm10$% for the bulk densities of likely analog meteorites, suggesting a rubble-pile internal structure. The associated mass estimate is $(7.8\pm0.9)\times 10^{10}\, \rm{kg}$ and $GM = 5.2\pm0.6\,\rm{m^3/s^2}$. Bennu's Earth close approaches are deterministic over the interval 1654--2135, beyond which the predictions are statistical in nature. In particular, the 2135 close approach is likely within the lunar distance and leads to strong scattering and therefore numerous potential impacts in subsequent years, from 2175--2196. The highest individual impact probability is $9.5\times 10^{-5}$ in 2196, and the cumulative impact probability is $3.7\times 10^{-4}$, leading to a cumulative Palermo Scale of -1.70.

We study the time evolution of the impact probability for synthetic but realistic impacting and close approaching asteroids detected in a simulated all-sky survey. We use the impact probability to calculate the impact warning time as the time interval between when an object reaches a Palermo Scale value of -2 and when it impacts Earth. A simple argument shows that warning time ~ D^x with the exponent in the range [1.0,1.5] and our derived value was $x=1.3\pm0.1$. The low-precision astrometry from the single simulated all-sky survey could require many days or weeks to establish an imminent impact for asteroids >100m diameter that are discovered far from Earth. Most close approaching asteroids are quickly identified as not being impactors but a size-dependent percentage, even for those >50m diameter, have a persistent impact probability of >10^{-6} on the day of closest approach. Thus, a single all-sky survey can be of tremendous value in identifying Earth impacting and close approaching asteroids in advance of their closest approach but it can not solve the problem on its own: high-precision astrometry from other optical or radar systems is necessary to rapidly establish an object as an impactor or close approacher. We show that the parallax afforded by surveying the sky from two sites is only of benefit for a small fraction of the smallest objects detected within a couple days before impact: probably not enough to justify the increased operating costs of a 2-site survey. Finally, the survey cadence within a fixed time span is relatively unimportant to the impact probability calculation. We tested three different reasonable cadences and found that one provided ~10x higher (better) value for the impact probability on the discovery night for the smallest (10m diameter) objects but the consequences on the overall impact probability calculation is negligible.

The number of asteroids with accurately determined orbits increases fast. The catalogs of asteroid physical observations have also increased, although the number of objects is still smaller than in the orbital catalogs. We developed a new approach to the asteroid family classification by combining the Hierarchical Clustering Method (HCM) with a method to add new members to existing families. This procedure makes use of the much larger amount of information contained in the proper elements catalogs, with respect to classifications using also physical observations for a smaller number of asteroids. Our work is based on the large catalog of the high accuracy synthetic proper elements (available from AstDyS). We first identify a number of core families; to these we attribute the next layer of smaller objects. Then, we remove all the family members from the catalog, and reapply the HCM to the rest. This gives both halo families which extend the core families and new independent families, consisting mainly of small asteroids. These two cases are discriminated by another step of attribution of new members and by merging intersecting families. By using information from absolute magnitudes, we take advantage of the larger size range in some families to analyze their shape in the proper semimajor axis vs. inverse diameter plane. This leads to a new method to estimate the family age (or ages). The results from the previous steps are then analyzed, using also auxiliary information on physical properties including WISE albedos and SDSS color indexes. This allows to solve some difficult cases of families overlapping in the proper elements space but generated by different collisional events. We analyze some examples of cratering families (Massalia, Vesta, Eunomia) which show internal structures, interpreted as multiple collisions. We also discuss why Ceres has no family.

We seek evidence of the Yarkovsky effect among Near Earth Asteroids (NEAs) by measuring the Yarkovsky-related orbital drift from the orbital fit. To prevent the occurrence of unreliable detections we employ a high precision dynamical model, including the Newtonian attraction of 16 massive asteroids and the planetary relativistic terms, and a suitable astrometric data treatment. We find 21 NEAs whose orbital fits show a measurable orbital drift with a signal to noise ratio (SNR) greater than 3. The best determination is for asteroid (101955) 1999 RQ36, resulting in the recovery of one radar apparition and an orbit improvement by two orders of magnitude. In addition, we find 16 cases with a lower SNR that, despite being less reliable, are good candidates for becoming stronger detections in the future. In some cases it is possible to constrain physical quantities otherwise unknown by means of the detected orbital drift. Furthermore, the distribution of the detected orbital drifts shows an excess of retrograde rotators that can be connected to the delivery mechanism from the most important NEA feeding resonances and allows us to infer the distribution for NEAs obliquity. We discuss the implications of the Yarkovsky effect for impact predictions. In particular, for asteroid (29075) 1950 DA our results favor a retrograde rotation that would rule out an impact in 2880.