Abstract A numerical solution to the Earth’s precession in a relativistic framework for a long time span is presented here. We obtain the motion of the solar system in the Barycentric Celestial Reference System by numerical integration with a symplectic integrator. Special Newtonian corrections accounting for tidal dissipation are included in the force model. The part representing Earth’s rotation is calculated in the Geocentric Celestial Reference System by integrating the post-Newtonian equations of motion published by Klioner et al. All the main relativistic effects are included following Klioner et al. In particular, we consider several relativistic reference systems with corresponding time scales, scaled constants and parameters. Approximate expressions for Earth’s precession in the interval ±1 Myr around J2000.0 are provided. In the interval ±2000 years around J2000.0, the difference compared to the P03 precession theory is only several arcseconds and the results are consistent with other long-term precession theories.
Keywords astrometry — ephemeris — post-Newtonian
It accepts original submissions from all over the world and is internationally published and distributed by IOP