Current approaches to uncertainty propagation in astrodynamics mainly refer to linearized models or Monte Carlo simulations. Naive linear methods fail in nonlinear dynamics, whereas Monte Carlo simulations tend to be computationally intensive. Differential algebra has already proven to be an efficient compromise by replacing thousands of pointwise integrations of Monte Carlo runs with the fast evaluation of the arbitrary order Taylor expansion of the flow of the dynamics. However, the current implementation of the DA-based high-order uncertainty propagator fails when the non-linearities of the dynamics prohibit good convergence of the Taylor expansion in one or more directions. We solve this issue by introducing automatic domain splitting. During propagation, the polynomial expansion of the current state is split into two polynomials whenever its truncation error reaches a predefined threshold. The resulting set of polynomials accurately tracks uncertainties, even B Alexander Wittig