Bounded motion design in the Earth zonal problem using differential algebra based normal form methods