Conjunction assessment requires knowledge of the uncertainty in the predicted orbit. Errors in the atmospheric density are a major source of error in the prediction of low Earth orbits. Therefore, accurate estimation of the density and quantification of the uncertainty in the density is required. Most atmospheric density models, however, do not provide an estimate of the uncertainty in the density. In this work, we present a new approach to quantify uncertainties in the density and to include these for calculating the probability of collision Pc. For this, we employ a recently developed dynamic reduced-order density model that enables efficient prediction of the thermospheric density. First, the model is used to obtain accurate estimates of the density and of the uncertainty in the estimates. Second, the density uncertainties are propagated forward simultaneously with orbit propagation to include the density uncertainties for Pc calculation. For this, we account for the effect of cross-correlation in position uncertainties due to density errors on the Pc. Finally, the effect of density uncertainties and cross-correlation on the Pc is assessed. The presented approach provides the distinctive capability to quantify the uncertainty in atmospheric density and to include this uncertainty for conjunction assessment while taking into account the dependence of the density errors on location and time. In addition, the results show that it is important to consider the effect of cross-correlation on the Pc, because ignoring this effect can result in severe underestimation of the collision probability.

Inaccurate estimates of the thermospheric density are a major source of error in low Earth orbit prediction. To improve orbit prediction, real-time density estimation is required. In this work, we develop a reduced-order dynamic model for the thermospheric density by computing the main spatial modes of the atmosphere and deriving a linear model for the dynamics. The model is then used to estimate the density using two-line element (TLE) data by simultaneously estimating the reduced-order modes and the orbits and ballistic coefficients of several objects using an unscented Kalman filter. Accurate density estimation using the TLEs of 17 objects is demonstrated and validated against CHAMP and GRACE accelerometer-derived densities. Finally, the use of the model for density forecasting is shown.

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.