Long-term evolution of orbits about a precessing oblate planet: 3. A semianalytical and a purely numerical approach

Abstract

Construction of an accurate theory of orbits about a precessing and nutating oblate planet, in terms of osculating elements defined in a frame associated with the equator of date, was started in Efroimsky and Goldreich (2004) and Efroimsky (2004, 2005, 2006a, b) .H ere we continue this line of research by combining that analytical machinery with numerical tools. Our model includes three factors: the J2 of the planet, its nonuniform equinoctial pre- cession described by the Colombo formalism, and the gravitational pull of the Sun. This semianalytical and seminumerical theory, based on the Lagrange planetary equations for the Keplerian elements, is then applied to Deimos on very long time scales (up to 1 billion years). In parallel with the said semianalytical theory for the Keplerian elements defined in the co-precessing equatorial frame, we have also carried out a completely independent, purely numerical, integration in a quasi-inertial Cartesian frame. The results agree to within fractions of a percent, thus demonstrating the applicability of our semianalytical model over long timescales. Another goal of this work was to make an independent check of whether the equinoctial-precession variations predicted for a rigid Mars by the Colombo model could We use the term "precession" in its general meaning, which includes any change of the instantaneous spin axis. So generally defined precession embraces the entire spectrum of spin-axis variations—from the polar wander and nutations through the Chandler wobble through the equinoctial precession.