Long term nonlinear propagation of uncertainties in perturbed geocentric dynamics using automatic domain splitting