The end-of-life disposal of Galileo satellites is needed to avoid collisions with operational spacecraft and to prevent the generation of space debris. Either disposal in stable graveyard orbits or disposal into the atmosphere exploiting eccentricity growth caused by lunisolar resonances are possible. However, there is a concern about the predictability of MEO orbits because of possible chaotic behaviour caused by the overlap of resonances. In this work, we investigate if Galileo disposal orbits are predictable and if safe disposal is possible under initial uncertainties. For this, we employ finite-time Lyapunov exponents (FTLE) and sensitivity analysis and compare these two methods regarding their practicality for analysing the predictability of disposal orbits. The results show that FTLE are not suitable for determining if an orbit behaves chaotically or not on the time scale of interest. Sensitivity analysis, on the other hand, can be used to quantify the effect of uncertainties on the orbital evolution and to determine if safe disposal is possible. In addition, we show that reliable re-entry disposal is feasible with a limited \Delta V budget. However, safe disposal into a graveyard orbit is not always feasible considering uncertainties in the disposal manoeuvre and dynamical model when the available \Delta V for disposal is low.

The propagation and Poincar\'e mapping of perturbed Keplerian motion is a key topic in celestial mechanics and astrodynamics, e.g. to study the stability of orbits or design bounded relative trajectories. The high-order transfer map (HOTM) method enables efficient mapping of perturbed Keplerian orbits over many revolutions. For this, the method uses the high-order Taylor expansion of a Poincar\'e or stroboscopic map, which is accurate close to the expansion point. In this paper, we investigate the performance of the HOTM method using different element sets for building the high-order map. The element sets investigated are the classical orbital elements, modified equinoctial elements, Hill variables, cylindrical coordinates and Deprit's ideal elements. The performances of the different coordinate sets are tested by comparing the accuracy and efficiency of mapping low-Earth and highly-elliptical orbits perturbed by $J_2$ with numerical propagation. The accuracy of HOTM depends strongly on the choice of elements and type of orbit. A new set of elements is introduced that enables extremely accurate mapping of the state, even for high eccentricities and higher-order zonal perturbations. Finally, the high-order map is shown to be very useful for the determination and study of fixed points and centre manifolds of Poincar\'e maps.

A differential algebra based importance sampling method for uncertainty propagation and impact probability computation on the first resonant returns of Near Earth Objects is presented in this paper. Starting from the results of an orbit determination process, we use a differential algebra based automatic domain pruning to estimate resonances and automatically propagate in time the regions of the initial uncertainty set that include the resonant return of interest. The result is a list of polynomial state vectors, each mapping specific regions of the uncertainty set from the observation epoch to the resonant return. Then, we employ a Monte Carlo importance sampling technique on the generated subsets for impact probability computation. We assess the performance of the proposed approach on the case of asteroid (99942) Apophis. A sensitivity analysis on the main parameters of the technique is carried out, providing guidelines for their selection. We finally compare the results of the proposed method to standard and advanced orbital sampling techniques.