We propose a method to account for the Earth oblateness effect in preliminary orbit determination of satellites in low orbits with radar observations. This method is an improvement of the one described in (Gronchi et al 2015), which uses a pure Keplerian dynamical model. Since the effect of the Earth oblateness is strong at low altitudes, its inclusion in the model can sensibly improve the initial orbit, giving a better starting guess for differential corrections and increasing the chances to obtain their convergence. The input set consists of two tracks of radar observations, each one composed of at least 4 observations taken during the same pass of the satellite. A single observation gives the topocentric position of the satellite, where the range is very accurate, while the line of sight direction is poorly determined. From these data we can compute by a polynomial fit the values of the range and range rate at the mean epochs of the two tracks. In order to obtain a preliminary orbit we wish to compute the angular velocities, that is the rate of change of the line of sight. In the same spirit of (Gronchi et al 2015), we also wish to correct the values of the angular measurements, so that they fit the selected dynamical model if the same holds for the radial distance and velocity. The selected model is a perturbed Keplerian dynamics, where the only perturbation included is the secular effect of the $J_2$ term of the geopotential. The proposed algorithm models this problem with 8 equations in 8 unknowns.

The modern optical telescopes produce a huge number of asteroid observations, that are grouped into very short arcs (VSAs), each containing a few observations of the same object in one single night. To decide whether two VSAs, collected in different nights, refer to the same observed object we can attempt to compute an orbit with the observations of both arcs: this is called the linkage problem. Since the number of orbit computations to be performed is very large, we need efficient methods of orbit determination. Using the first integrals of Kepler's motion we can write algebraic equations for the linkage problem, which can be put in polynomial form. The equations introduced in (Gronchi et al. 2015) can be reduced to a univariate polynomial of degree 9: the unknown is the topocentric distance $\rho$ of the observed body at the mean epoch of one of the VSAs. Using elimination theory we show an optimal property of this polynomial: it has the least degree among the univariate polynomials in the same variable that are consequence of the algebraic conservation laws and are obtained without squaring operations, that can be used to bring these algebraic equations in polynomial form. In this paper we also introduce a procedure to join three VSAs belonging to different nights: from the conservation of angular momentum at the three mean epochs of the VSAs, we obtain a univariate polynomial equation of degree 8 in the topocentric distance $\rho_2$ at the intermediate epoch. This algorithm has the same computational complexity as the classical method by Gauss, but uses more information, therefore we expect that it can produce more accurate results. For both methods, linking two and three VSAs, we also discuss how to select the solutions, making use of the full two-body dynamics, and show some numerical tests comparing the results with the ones obtained by Gauss' method.

We report on our search for genetically related asteroids amongst the near-Earth object (NEO) population - families of NEOs akin to the well known main belt asteroid families. We used the technique proposed by Fu et al. (2005) supplemented with a detailed analysis of the statistical significance of the detected clusters. Their significance was assessed by comparison to identical searches performed on 1,000 'fuzzy-real' NEO orbit distribution models that we developed for this purpose. The family-free 'fuzzy-real' NEO models maintain both the micro and macro distribution of 5 orbital elements (ignoring the mean anomaly). Three clusters were identified that contain four or more NEOs but none of them are statistically significant at \geq 3{\sigma}. The most statistically significant cluster at the \sim 2{\sigma} level contains 4 objects with H < 20 and all members have long observational arcs and concomitant good orbital elements. Despite the low statistical significance we performed several other tests on the cluster to determine if it is likely a genetic family. The tests included examining the cluster's taxonomy, size-frequency distribution, consistency with a family-forming event during tidal disruption in a close approach to Mars, and whether it is detectable in a proper element cluster search. None of these tests exclude the possibility that the cluster is a family but neither do they confirm the hypothesis. We conclude that we have not identified any NEO families.

Given a set of astrometric observations of the same object, the problem of orbit determination is to compute the orbit and to assess its uncertainty and reliability. For the next generation surveys, with much larger number density of observed objects, new algorithms or substantial revisions of the classical ones are needed. The problem has three main steps, preliminary orbit, least squares orbit, and quality control. The classical theory of preliminary orbits was incomplete: the consequences of the topocentric correction had not been fully studied. We show that it is possible to account for this correction, possibly with an increase in the number of preliminary solutions, without impairing the overall orbit determination performance. We have developed modified least squares orbit determination algorithms that can be used to improve the reliability of the procedure. We have tested the complete procedure on two simulations with number densities comparable to that expected from the next generation surveys such as Pan-STARRS and LSST. To control the problem of false identifications we have introduced a quality control on the fit residuals based on an array of metrics and a procedure to remove duplications and contradictions in the output. The results confirm that large sets of discoveries can be obtained with good quality orbits and very high success rate losing only 0.6 to 1.3% of objects and a false identification rate in the range 0.02 to 0.06%.