Abstract The variation equation of the satellite orbit elements under both J2 and drag perturbations was derived based on the Gauss perturbation equation.