Semianalytic Integration of High-Altitude Orbits under Lunisolar Effects

The long-term effect of lunisolar perturbations on high-altitude orbits is studied after a double averaging procedure that removes both the mean anomaly of the satellite and that of the moon. Lunisolar effects acting on high-altitude orbits are comparable in magnitude to the Earth’s oblateness perturbation. Hence, their accurate modeling does not allow for the usual truncation of the expansion of the third-body disturbing function up to the second degree. Using canonical perturbation theory, the averaging is carried out up to the order where second-order terms in the Earth oblateness coefficient are apparent. This truncation order forces to take into account up to the fifth degree in the expansion of the lunar disturbing function. The small values of the moon’s orbital eccentricity and inclination with respect to the ecliptic allow for some simplification. Nevertheless, as far as the averaging is carried out in closed form of the satellite’s orbit eccentricity, it is not restricted to low-eccentricity orbits.