The aim of this work is to analyze the orbital evolution of the mean eccentricity given by the Two-Line Elements (TLE) set of the Molniya satellites constellation. The approach is bottom-up, aiming at a synergy between the observed dynamics and the mathematical modeling. Being the focus the long-term evolution of the eccentricity, the dynamical model adopted is a doubly-averaged formulation of the third-body perturbation due to Sun and Moon, coupled with the oblateness effect on the orientation of the satellite. The numerical evolution of the eccentricity, obtained by a two-degree-of-freedom model assuming different orders in the series expansion of the third-body effect, is compared against the mean evolution given by the TLE. The results show that, despite being highly elliptical orbits, the second order expansion catches extremely well the behavior. Also, the lunisolar effect turns out to be non-negligible for the behavior of the longitude of the ascending node and the argument of pericenter. The role of chaos in the timespan considered is also addressed. Finally, a frequency series analysis is proposed to show the main contributions that can be detected from the observational data.

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.

We use a secular model to describe the non-resonant dynamics of trans-Neptunian objects in the presence of an external ten-earth-mass perturber. The secular dynamics is analogous to an "eccentric Kozai mechanism" but with both an inner component (the four giant planets) and an outer one (the eccentric distant perturber). By the means of Poincar\'e sections, the cases of a non-inclined or inclined outer planet are successively studied, making the connection with previous works. In the inclined case, the problem is reduced to two degrees of freedom by assuming a non-precessing argument of perihelion for the perturbing body. The size of the perturbation is typically ruled by the semi-major axis of the small body: we show that the classic integrable picture is still valid below about 70 AU, but it is progressively destroyed when we get closer to the external perturber. In particular, for a>150 AU, large-amplitude orbital flips become possible, and for a>200 AU, the Kozai libration islands are totally submerged by the chaotic sea. Numerous resonance relations are highlighted. The most large and persistent ones are associated to apsidal alignments or anti-alignments with the orbit of the distant perturber.

We use a secular representation to describe the long-term dynamics of transneptunian objects in mean-motion resonance with Neptune. The model applied is thoroughly described in Saillenfest et al. (2016). The parameter space is systematically explored, showing that the secular trajectories depend little on the resonance order. High-amplitude oscillations of the perihelion distance are reported and localised in the space of the orbital parameters. In particular, we show that a large perihelion distance is not a sufficient criterion to declare that an object is detached from the planets. Such a mechanism, though, is found unable to explain the orbits of Sedna or 2012VP113, which are insufficiently inclined (considering their high perihelion distance) to be possibly driven by such a resonant dynamics. The secular representation highlights the existence of a high-perihelion accumulation zone due to resonances of type 1:k with Neptune. That region is found to be located roughly at semi-major axes in [100;300] AU, perihelion distances in [50;70] AU and inclinations [30;50]{\deg}. In addition to the flux of objects directly coming from the Scattered Disc, numerical simulations show that the Oort Cloud is also a substantial source for such objects. Naturally, as that mechanism relies on fragile captures in high-order resonances, our conclusions break down in the case of a significant external perturber. The detection of such a reservoir could thus be an observational constraint to probe the external Solar System.

Two semi-analytical one-degree-of-freedom secular models are presented for the motion of small bodies beyond Neptune. A special attention is given to trajectories entirely exterior to the planetary orbits. The first one is the well-known non-resonant model of Kozai (1962) adapted to the transneptunian region. Contrary to previous papers, the dynamics is fully characterized with respect to the fixed parameters. A maximum perihelion excursion possible of 16.4 AU is determined. The second model handles the occurrence of a mean-motion resonance with one of the planets. In that case, the one-degree-of-freedom integrable approximation is obtained by postulating the adiabatic invariance, and is much more general and accurate than previous secular models found in the literature. It brings out in a plain way the possibility of perihelion oscillations with a very high amplitude. Such a model could thus be used in future studies to deeper explore that kind of motion. For complex resonant orbits (especially of type 1:k), a segmented secular description is necessary since the trajectories are only "integrable by parts". The two models are applied to the Solar System but the notations are kept general so that it could be used for any quasi-circular and coplanar planetary system.

Juno is a NASA mission launched in 2011 with the goal of studying Jupiter. The probe will arrive to the planet in 2016 and will be placed for one year in a polar high-eccentric orbit to study the composition of the planet, the gravity and the magnetic field. The Italian Space Agency (ASI) provided the radio science instrument KaT (Ka-Band Translator) used for the gravity experiment, which has the goal of studying the Jupiter's deep structure by mapping the planet's gravity: such instrument takes advantage of synergies with a similar tool in development for BepiColombo, the ESA cornerstone mission to Mercury. The Celestial Mechanics Group of the University of Pisa, being part of the Juno Italian team, is developing an orbit determination and parameters estimation software for processing the real data independently from NASA software ODP. This paper has a twofold goal: first, to tell about the development of this software highlighting the models used, second, to perform a sensitivity analysis on the parameters of interest to the mission.

The radioscience experiment is one of the on board experiment of the Mercury ESA mission BepiColombo that will be launched in 2014. The goals of the experiment are to determine the gravity field of Mercury and its rotation state, to determine the orbit of Mercury, to constrain the possible theories of gravitation (for example by determining the post-Newtonian (PN) parameters), to provide the spacecraft position for geodesy experiments and to contribute to planetary ephemerides improvement. This is possible thanks to a new technology which allows to reach great accuracies in the observables range and range rate; it is well known that a similar level of accuracy requires studying a suitable model taking into account numerous relativistic effects. In this paper we deal with the modelling of the space-time coordinate transformations needed for the light-time computations and the numerical methods adopted to avoid rounding-off errors in such computations.