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.

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.

We study the local shadow of the Schwarzschild black hole with a quadrupole distortion and the influence of the external gravitational field on the photon dynamics. The external matter sources modify the light ring structure and lead to the appearance of multiple shadow images. In the case of negative quadrupole moments we identify the most prominent mechanism causing multiple shadow formation. Furthermore, we obtain a condition under which this mechanism can be realized. This condition depends on the quadrupole moment, but also on the position of the observer and the celestial sphere.

In a recent letter, arXiv:1509.00021, it was shown that the lensing of light around rotating boson stars and Kerr black holes with scalar hair can exhibit chaotic patterns. Since no separation of variables is known (or expected) for geodesic motion on these backgrounds, we examine the 2D effective potentials for photon trajectories, to obtain a deeper understanding of this phenomenon. We find that the emergence of stable light rings on the background spacetimes, allows the formation of "pockets" in one of the effective potentials, for open sets of impact parameters, leading to an effective trapping of some trajectories, dubbed quasi-bound orbits. We conclude that pocket formation induces chaotic scattering, although not all chaotic orbits are associated to pockets. These and other features are illustrated in a gallery of examples, obtained with a new ray-tracing code, PYHOLE, which includes tools for a simple, simultaneous visualization of the effective potential together with the spacetime trajectory, for any given point in a lensing image. An analysis of photon orbits allows us to further establish a positive correlation between photon orbits in chaotic regions and those with more than one turning point in the radial direction; we recall that the latter is not possible around Kerr black holes. Moreover, we observe that the existence of several light rings around a horizon (several fundamental orbits, including a stable one), is a central ingredient for the existence of multiple shadows of a single hairy black hole. We also exhibit the lensing and shadows by Kerr black holes with scalar hair, observed away from the equatorial plane, obtained with PYHOLE.