# Specific effects of large asteroids on the orbits of terrestrial planets and the ASETEP database

S. Aljbaae, J. Souchay

► To cite this version:

S. Aljbaae, J. Souchay. Specific effects of large asteroids on the orbits of terrestrial planets and the ASETEP database. Astronomy and Astrophysics - A&A, 2012, 540, pp.A21. 10.1051/0004-6361/201118564. hal-02480684

HAL Id: hal-02480684

https://hal.sorbonne-universite.fr/hal-02480684

Submitted on 2 Oct 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Specific effects of large asteroids on the orbits of terrestrial planets and the ASETEP database*

S. Aljbaae and J. Souchay

Observatoire de Paris, SYRTE/CNRS UMR8630, 61 avenue de l'Observatoire, 75014 Paris, France
e-mail: safwan.aljbaae; jean.souchay@obspm.fr

Received 2 December 2011 / Accepted 16 January 2012

## ABSTRACT

The necessity to take into account the perturbations caused by a large number of asteroids on the terrestrial planets is fundamental in the construction of modern numerical ephemeris on the solar system. Therefore about 300 of the largest asteroids were taken into account in recent ephemeris. Yet, the uncertainty on the mass values of the great majority of these asteroids constitutes a crucial and the main limit of accuracy of this ephemeris. Consequently, it is important to conduct a specific and detailed study of their individual effects especially on the terrestrial planets, which are far more affected than the giant planets. This was already done explicitly, but only for Mars and for only two orbital elements (a and $\lambda$). We aim both to confirm these previous results and to extend the study to all orbital elements and to the other three terrestrial planets (Mercury, Venus and the Earth), which are priori less affected by asteroid perturbations. Our methodology consists in several steps: we carried out precise computations of the orbital motions of the planets at short (100 y) and longer (1000 y) time scales with numerical integration. For that purpose we included the eight planets and also considered 43 of the most powerful asteroids. These were added to the numerical integrations once separately and once combined to determine their specific effects on the orbital elements of the Earth and the three other terrestrial planets. This procedure also allowed us to assess the spatial geocentric coordinates of the three terrestrial planets. We determined the signal that represents the effects by simple subtraction. Then we systematically analyzed this signal by FFT (fast Fourier transform), and finally we adjusted the signal with a set of sinusoidal components. We analyzed in detail the variations of the six orbital elements *a*, *e*, *i*, $\Omega$, $\tilde{\omega}$ and $\lambda$ of Mercury, Venus, the Earth-Moon barycenter (EMB) and Mars that are caused by the individual influences of the set of our 43 selected asteroids. We compared our results for Mars with the analytical ones on the semi-major axis and the longitude. The tow studies agree very well. All our results, consisting of 1032 different curves (43 asteroids $\times$ 4 planets $\times$ 6 orbital elements) and the related tables that provide the fitted Fourier and Poisson components are gathered the ASETEP database (asteroid effect on the terrestrial planets). Moreover, we include in this database the influence of our sample of 43 asteroids on three fundamental parameters: the distance and the bi-dimensional orientation vector $(\alpha, \delta)$ from the EMB to each of the other terrestrial planets. This database, which will be regularly updated by taking into account more asteroids with improved mass determinations, constitutes a precious tool for understanding specifically the influence of the large asteroids on the orbital motion of the terrestrial planets, and also for better understanding how modern ephemeris can be improved.

**Key words.** gravitation – minor planets, asteroids: general – ephemerides – astronomical databases: miscellaneous – astrometry – celestial mechanics

## 1. Introduction

There are two main ways of computing the ephemeris of celestial bodies. One is semi-analytical, which was used for the VSOP87 theory of the planets by Bretagnon et al (1987) and the ELP2000 solution for the orbital motion of the Moon (Chapront et al. 1988). In these successful and complex works, the coordinates of the moving bodies are given in the form of Fourier series of fundamental arguments as the mean longitudes of the planets or the Delaunay arguments. The interest of these analytical theories is that they enable us to understand the origin of the countless sinusoidal terms that originate in the perturbing function. But these theories can no longer reach the necessary accuracy to match numerical determinations of the motions that constitute the second way of computation.

Indeed, today the stability and rapidity of the numerical integrators allow an accurate construction of the ephemeris of the solar system's main bodies. As explained by Standish (2004), the accuracy of an ephemeris depends on three major factors. The first one is related to the equations of motion. These equations

might be imperfect for instance because they do not take all the gravitational effects into account or because the underlying analytical model contains some errors. The second factor concerns the integration of motion, which inevitably contains uncertainties that have to be estimated. The last factor is the exact determination of basic parameters such as the initial conditions (positions, velocities, masses etc.), which are given by the observations, and therefore contain errors that are conditioning the results of the integration. The two first factors are generally believed to be satisfactory, because the equations of motion look correct and complete, whereas the errors accompanying the integration algorithms have been proved to be negligible. The third factor is still subject to various investigations and numerous recent developments, in particular those coming from the recent missions to Mars (Standish 2004; Fienga et al. 2009). The JPL ephemeris are undoubtedly used the most worldwide. Since the construction of an initial release DE200 in 1979, a lot of other versions have been constructed, such as DE405 in 1997, and DE409 in 2003. This last version in particular uses very recent VLBI observations of the MGS (Mars Global Surveyor) in 2001–2003 and Odyssey in 2002–2003, which leads to an accuracy for the determination of the position of Mars at the level

<footnote>* Appendices A–C are available in electronic form at
http://www.aanda.org</footnote>
of a few milliarcseconds. The comparisons between the various versions of JPL ephemerides give a clear insight into the improvements. DE405 for instance, showed an improvement of an order of magnitude over DE200, especially for the inner planets, as explained by Standish (2004).

To evaluate the accuracy limits in the ephemeris of the four inner planets, Standish et al. (2002) considered the present-day knowledge of the asteroid masses determined by various methods, and modeled the distributions of those masses. Then they carried out a Monte-Carlo study, repeatedly adjusting the ephemeris to fit the observational data, by a suitable set of asteroid masses. the authors’ conclusion is that by combining accumulated data for more than two decades, the orbits become distorted both in right ascension and in declination by as much as 5 km or more. At the level of present-time ephemerides accuracy, the influence of asteroids can no longes be neglected, as in the case of DE409. Many asteroids have enough mass to produce a significant gravitational effect on the inner planets. Nevertheless, because of the large uncertainty in the determination of asteroid masses, it is difficult or impossible to evaluate their effect with sufficient accuracy. Yet no specific study of the influence of each asteroid has been conducted to date. A pioneering work was carried out by Williams (1984), giving a list of the most influential asteroids on the orbit of Mars over a short period of time, roughly 50 years. Recently, Mouret et al. (2009) identified and provided a list of the main perturbers of Mars, which could cause significant disturbances over the next 50 years as well. Their analytical analysis is based on a harmonic decomposition of the perturbations on Mars as a function of the mean longitudes of the planet and the asteroid, to the first order in asteroid mass. Still, these two last studies (Williams 1984; Mouret et al. 2009) can be considered as incomplete because they concentrate mainly on two orbital elements, the semi-major axis and longitude of Mars. Moreover, they only focus on this planet.

We aim to carry out an extensive study of the subject in the present paper. In Sect. 2, after determining the asteroidal mass and its accuracy, we select a set of 43 asteroids. We explain in detail in Sect. 3 the theoretical framework and our methodology with which we calculated the perturbing effect of any specific asteroid on a terrestrial planet. In Sect. 4, we determine the effects of our 43 selected asteroids on the six Mars orbital elements, while the effects on the other telluric planets, the Earth (replaced in the following by the EMB (Earth-Moon barycenter), Mercury and Venus are presented in Sect. 5. Thus we reconstruct the asteroidal perturbation as a combination of Fourier and Poisson components for each perturbing asteroid as well as each terrestrial planet and each orbital parameter. Moreover, as an important by-product of these computations, we calculate in Sect. 6 the influences of each asteroid on two fundamental astrometric parameters namely the distance from the EMB to the barycenter of the other terrestrial planets (Mercury, Venus and Mars) and their two-dimensional positioning (α, δ) in the celestial sphere. Finally we describe in Sect. 7 how we constructed our ASETEP (asteroidal effects on the terrestrial planets) database, which enables one to determine the influence of each asteroid on any of the four terrestrial planets precisely and in detail through their orbital elements and their geocentric positioning.

## 2. Mass determinations of the asteroids

Much information about the size, shape and pole direction of the leading asteroids has been obtained from analyses of light curves at different epochs, occultation, and/or images obtained through adaptive optics and speckle techniques. The mass is one of the

important characteristics, necessary for describing the perturbations on other bodies. There are two methods of estimating the mass: one is dynamical, the other is astrophysical.

The dynamical method is based on the analysis of perturbations of the studied asteroid on other celestial bodies. After a pioneering work by Williams (1984), studying how asteroid masses can be determined from their perturbations on Mars, Standish et al. (1989) showed that the measurements of Mars orbit perturbations can be conducted to estimate the masses of Ceres, Pallas and Vesta, with the intermediary of a least-squares adjustment of all ephemeris parameters, including the relevant asteroid masses. This adjustment was compiled on the basis of observational data encompassing the highly accurate Viking Lander ranging data. The accuracy of mass determination through this method highly depends on the precision of positional observations used and on the number of the observations and their distribution.

The astrophysical method is based on the data about the asteroid size and the taxonomic classification. This method assumes that the body has a spherical shape, a homogeneous structure, is single, etc. This method, in contrast to the dynamical one, has the advantage that it can be applied to many objects even if there is no close encounter with other bodies (Kochetova 2004), but is generally marred by a large uncertainty.

### 2.1. The mass of Ceres, Pallas and Vesta

The masses of the three largest asteroids, Ceres, Pallas and Vesta obtained by Standish et al. (1989) starting from the perturbations on Mars orbital motion were $(5.0 \pm 0.2) \times 10^{-10} M_S$, $(1.4 \pm 0.2) \times 10^{-10} M_S$ and $(1.5 \pm 0.3) \times 10^{-10} M_S$ respectively. They were updated by more recent studies, as we will show below. Hilton (1999) explained why the mass determination of Ceres is completely dependent on that of Pallas. The reason is the degeneracy between the two masses, which is in turn caused by the similarity of the semi-major axes of the two asteroids. In consequence, their mean motion is nearly the same, and their mutual synodic period is very long, around 2000 yr (Williams 1989). This amplifies the propagation of the errors on the mass of Pallas to the mass of Ceres, when the mass of the first asteroid is obtained starting from a series of close encounters with the second one. This explains why the mass estimations are highly dependent on old data with poor accuracy. Because of the high Pallas inclination and its high eccentricity, its close encounters with main-belt asteroids occur only rarely, and when they occur, the relative velocity is high so that mass determination is not very efficient. For Vesta, the mass was obtained starting from perturbations of both Ceres and 197 Arete. The values obtained by Hilton (1999) were $M_{Ceres} = 4.39(\pm0.04) \times 10^{-10} M_S$, $M_{Pallas} = 1.59(\pm0.05) \times 10^{-10} M_S$ and $M_{Vesta} = 1.69(\pm0.11)) \times 10^{-10} M_S$. Many new constraints for the masses of the three asteroids were obtained during one decade. For our computations, we have chosen the recent values determined by Pitjeva (2009) from data of the Martian landers and spacecraft by fitting the planet ephemeris to observations. Their values are significantly different from the estimations above, i.e. $4.712(\pm0.006) \times 10^{-10} M_S$ for Ceres, $1.027(\pm0.007) \times 10^{-10} M_S$ for Pallas, and $1.344(\pm0.003) \times 10^{-10} M_S$ for Vesta.

### 2.2. The mass of the other large asteroids

We limit our study to the perturbations produced by asteroids with relatively large and well-constrained masses. According to Mouret et al. (2009), we can rely on 36 asteroid masses derived
from direct measurements with a precision better than 10 per cent, while only 54 reach the 50 per cent mark. These statistics come from Mouret et al. (2007), Baer et al. (2008a) and Baer database<fnref n="1" />. All these masses were deduced from the gravitational pull exerted by the corresponding asteroid on other solar system objects, mainly smaller asteroids, and in special cases Mars or a spacecraft. In Table 1 we present the set of the 43 asteroids selected in this paper for which we will calculate the gravitational effects on the four terrestrial planets in the following, together with their mass determination and the relative uncertainty. Our choice was guided by two criteria: the asteroid must be big enough and its mass determination satisfactorily accurate. The data were taken from 12 different sources: Pitjeva et al. (2009), Fienga et al. (2009), Baer et al. (2008a), Baer et al. (2008b), Bange (1998), Marchis et al. (2008a), Marchis et al. (2008b), Kochetova (2004), Ivantsov (2008), Aslan et al. (2007), and Yeomans et al. (1997). Notice that in the near future, Gaia may significantly increase the number of favorable encounters to be analyzed for mass determination, as explained by Mouret et al. (2009).

### 3. Analytical expression of planet perturbations caused by an asteroid

After reviewing the classical theory of perturbations for a system of N planets orbiting the Sun, we show in this section how the individual effects of a given asteroid on one of these planets are determined analytically and explain our methodology of characterizing these effects with numerical integration.

#### 3.1. Theoretical principles

Consider N planets $P_1, ..., P_N$ orbiting around the Sun $P_0$, with mass $m_0$, according to Newton's law. The heliocentric position vector $\mathbf{r}_i$ of a given planet in the frame of this N+1 body problem is subject to the second order differential equation:

$$ \frac{d^2 \mathbf{r}_i}{dt^2} = \operatorname{grad}_{\mathbf{r}_i} \left( \frac{G(m_0 + m_i)}{r_i} + R_i \right), \quad (1) $$

where $i = 1, 2, ..., N$ and the parameters $m_i$ and $G$ are the mass of the ith. planet and the gravitational constant. $R_i$, usually called the disturbing function on the planet $P_i$, is written as follows:

$$ R_i = \sum_{j \neq i} Gm_j \left( \frac{1}{|\mathbf{r}_i - \mathbf{r}_j|} - \frac{\mathbf{r}_i \cdot \mathbf{r}_j}{r_j^3} \right). \quad (2) $$

This disturbing function causes a perturbation of the planet's osculating orbital elements; in principle each of these perturbations can be calculated analytically according to the classical Lagrange equations for the osculating elements $a, e, i, \Omega, \varpi$ and $\lambda$ (see for instance Murray & Dermott 1999) (for simplicity we did not take into account the index $i$ in $R_i$)

$$ \begin{aligned} \frac{da}{dt} &= \frac{2}{na} \frac{\partial R}{\partial \lambda} \\ \frac{de}{dt} &= -\frac{\sqrt{1-e^2}}{na^2e} \left(1 - \sqrt{1-e^2}\right) \frac{\partial R}{\partial \lambda} - \frac{\sqrt{1-e^2}}{na^2e} \frac{\partial R}{\partial \varpi} \end{aligned} $$

$$ \begin{aligned} \frac{d\lambda}{dt} &= n - \frac{2}{na} \frac{\partial R}{\partial a} + \frac{\sqrt{1 - e^2}(1 - \sqrt{1 - e^2})}{na^2e} \frac{\partial R}{\partial e} + \frac{\tan \frac{1}{2} I}{na^2 \sqrt{1 - e^2}} \frac{\partial R}{\partial I} \\ \frac{d\Omega}{dt} &= \frac{1}{na^2 \sqrt{1 - e^2} \sin I} \frac{\partial R}{\partial I} \\ \frac{d\varpi}{dt} &= \frac{\sqrt{1 - e^2}}{na^2e} \frac{\partial R}{\partial e} + \frac{\tan \frac{1}{2} I}{na^2 \sqrt{1 - e^2}} \frac{\partial R}{\partial I} \\ \frac{dI}{dt} &= \frac{-\tan \frac{1}{2} I}{na^2 \sqrt{1 - e^2}} \left( \frac{\partial R}{\partial \lambda} + \frac{\partial R}{\partial \varpi} \right) - \frac{1}{na^2 \sqrt{1 - e^2} \sin I} \frac{\partial R}{\partial \Omega}. \end{aligned} \quad (3) $$

We focus on the perturbation caused on a planet P (one of the $P_i$) by an asteroid P' of mass $m'$ added to the N+1 body system considered above. Let the position vector of the new body be denoted by $\mathbf{r}'$. Generally we define the perturbation $\Delta X$ of one orbital element X of the planet induced by the asteroid as the difference $\Delta X = \bar{X} - X$, where $\bar{X}$ stands for the perturbed element in the frame of the N+2 body problem (with the asteroid) and X for perturbed element in the frame of the N+1 body problem (without the asteroid). On an equivalent manner, the perturbing function $\Delta R$ can be written as

$$ \Delta R = \bar{R} - R = Gm' \left( \frac{1}{|\mathbf{r} - \mathbf{r}'|} - \frac{\mathbf{r} \cdot \mathbf{r}'}{|\mathbf{r}'|^3} \right). \quad (4) $$

We can write the derivative of the perturbation induced by the extra body P' on the elements of the planet in the following form (Kuchynka 2010):

$$ \begin{aligned} \frac{d\Delta a}{dt} &= \frac{2}{na} \frac{\partial \Delta R}{\partial \lambda} \\ \frac{d\Delta e}{dt} &= -\frac{\sqrt{1 - e^2}}{na^2e} \left(1 - \sqrt{1 - e^2}\right) \frac{\partial \Delta R}{\partial \lambda} - \frac{\sqrt{1 - e^2}}{na^2e} \frac{\partial \Delta R}{\partial \varpi} \\ \frac{d\Delta \lambda}{dt} &= \Delta n - \frac{2}{na} \frac{\partial \Delta R}{\partial a} + \frac{\sqrt{1 - e^2} \left(1 - \sqrt{1 - e^2}\right)}{na^2e} \frac{\partial \Delta R}{\partial e} \\ &\quad + \frac{\tan \frac{1}{2} I}{na^2 \sqrt{1 - e^2}} \frac{\partial \Delta R}{\partial I} \\ \frac{d\Delta \Omega}{dt} &= \frac{1}{na^2 \sqrt{1 - e^2} \sin I} \frac{\partial \Delta R}{\partial I} \\ \frac{d\Delta \varpi}{dt} &= \frac{\sqrt{1 - e^2}}{na^2e} \frac{\partial \Delta R}{\partial e} + \frac{\tan \frac{1}{2} I}{na^2 \sqrt{1 - e^2}} \frac{\partial \Delta R}{\partial I} \\ \frac{d\Delta I}{dt} &= \frac{-\tan \frac{1}{2} I}{na^2 \sqrt{1 - e^2}} \left( \frac{\partial \Delta R}{\partial \lambda} + \frac{\partial R}{\partial \varpi} \right) - \frac{1}{na^2 \sqrt{1 - e^2} \sin I} \frac{\partial \Delta R}{\partial \Omega}. \end{aligned} \quad (5) $$

The development of the perturbing function as a function of eccentricities and inclinations enables one to obtain the expressions of the perturbations as a function of the orbital elements both of the perturbed body (here the planet) and the perturbing body (the studied asteroid). This development was continued up to the fourth order by Brouwer & Clemence (1961). To develop the disturbing function $\Delta R$ given in Eq. (4), we distinguished between its direct part $R_D$ and its indirect part $R_E$. According to classical celestial mechanics (see for instance Laskar et al. 1995; or Murray & Dermott 1999) we have:

$$ \Delta R = \frac{Gm'}{a'} (R_D + \alpha R_E), \quad (6) $$

<footnote n="1">¹ URL: http://home.earthlink.net/~jimbaer1/astmass.htm</footnote>
With $\alpha = \frac{a'}{a}$, where $a'$ stands for the semi-major axis of the asteroid, and:

$$ R_D = \frac{a'}{|\mathbf{r} - \mathbf{r}'|} + (C_1^j + (e^2 + e'^2)C_2^j + (s^2 + s'^2)C_3^j)\cos[j\lambda - j\lambda'] + ee'C_4^j \cos[j\lambda' - j\lambda + \varpi' - \varpi] + ss'C_6^j \cos[j\lambda' - j\lambda + \Omega' - \Omega] + eC_7^j \cos[j\lambda' + (1-j)\lambda - \varpi] + e'C_8^j \cos[j\lambda' + (1-j)\lambda - \varpi'] + e^2C_9^j \cos[j\lambda' + (2-j)\lambda - 2\varpi] + ee'C_10^j \cos[j\lambda' + (2-j)\lambda - \varpi' - \varpi] + e'C_11^j \cos[j\lambda' + (2-j)\lambda - 2\varpi'] + e'C_12^j \cos[j\lambda' + (2-j)\lambda - 2\Omega] + ss'C_13^j \cos[j\lambda' + (2-j)\lambda - \Omega' - \Omega] + s'^2C_14^j \cos[j\lambda' + (2-j)\lambda - 2\Omega'] \quad (7) $$

$$ R_E = -\frac{a'^2}{a} \frac{\mathbf{r} \cdot \mathbf{r}'}{|\mathbf{r} - \mathbf{r}'|^3} = -\frac{r}{a} \left(\frac{a'}{|\mathbf{r}'|}\right)^2 \cos \psi = \begin{aligned}[t] & -\left(-1 + \frac{1}{2}e^2 + \frac{1}{2}e'^2 + s^2 + s'^2\right) \cos(\lambda' - \lambda) \\ & -ee' \cos(2\lambda' - 2\lambda - \varpi' + \varpi) - \frac{1}{2}e \cos(\lambda' - 2\lambda + \varpi) \\ & -2ss' \cos(\lambda' - \lambda - \Omega' + \Omega) + \frac{3}{2}e \cos(\lambda' - \varpi) \\ & -2e' \cos(2\lambda' - \lambda - \varpi') - \frac{3}{8}e^2 \cos(\lambda' - 3\lambda + 2\varpi) \\ & -\frac{1}{8}e^2 \cos(\lambda' + \lambda - 2\varpi) + 3ee' \cos(2\lambda - \varpi' - \varpi) \\ & -\frac{1}{8}e'^2 \cos(\lambda' + \lambda - 2\varpi') - s'^2 \cos(\lambda' + \lambda - 2\Omega') \\ & -s^2 \cos(\lambda' + \lambda - 2\Omega) + 2ss' \cos(\lambda' + \lambda - \Omega' - \Omega) \\ & -\frac{27}{8}e'^2 \cos(3\lambda' - \lambda - 2\varpi') \end{aligned} \quad (8) $$

where $e'$, $\lambda'$, $\varpi'$, $\Omega'$ stand for the orbital elements of the asteroid. $\psi$ is the angle between the position vectors $\mathbf{r}$ and $\mathbf{r}'$, $s = \sin \frac{1}{2}i$ $s' = \sin \frac{1}{2}i'$. The coefficients $C_i^j$ depend on the Laplace coefficients

$$ b_s^j(\alpha) = \frac{1}{\pi} \int_0^{2\pi} \frac{\cos j\rho \, d\rho}{(1 - 2\alpha \cos \rho + \alpha^2)^s} $$

The developments in (7) and (8) show that after integration of the variational Eqs. (5) we can expect that the variations of orbital elements be expressed as a series of linear combinations of $\lambda$ and $\lambda'$, taking into account that the parameters $\tilde{\omega}$, $\tilde{\omega}'$, $\tilde{\Omega}$ and $\tilde{\Omega}'$ have a very low frequency with respect both to $\lambda$ and $\lambda'$. Therefore, they can be considered as nearly constant during the interval of time considered (100 y or 1000 y). This explains why both Williams (1984) and Mouret et al. (2009) constructed their developments of $a_{\text{Mars}}$ and $\lambda_{\text{Mars}}$ only with the help of sinusoidal terms with high frequency, function of the arguments of the form $j\lambda + j'\lambda'$. Nevertheless this study is incomplete, because it does not take into account two long-periodic terms that can come from two different origins: first the resonant terms in (7) and (8) with very

<table><caption>**Table 1.** Mass estimation of the 43 asteroids selected in this paper, with their relative uncertainty.</caption><thead><tr><th></th><th>Name</th><th>M [M<sub>☉</sub>]</th><th>ΔM [%]</th><th>Authors</th></tr></thead><tbody><tr><td>1</td><td>Ceres</td><td>0.472 × 10<sup>-9</sup></td><td>0.127</td><td>Pitjeva et al. (2009)</td></tr><tr><td>2</td><td>Pallas</td><td>0.103 × 10<sup>-9</sup></td><td>0.682</td><td>Pitjeva et al. (2009)</td></tr><tr><td>3</td><td>Juno</td><td>0.750 × 10<sup>-11</sup></td><td>20.000</td><td>Fienga et al. (2009)</td></tr><tr><td>4</td><td>Vesta</td><td>0.134 × 10<sup>-9</sup></td><td>0.223</td><td>Pitjeva et al. (2009)</td></tr><tr><td>6</td><td>Hebe</td><td>0.759 × 10<sup>-11</sup></td><td>18.709</td><td>Baer et al. (2008a)</td></tr><tr><td>7</td><td>Iris</td><td>0.500 × 10<sup>-11</sup></td><td>8.621</td><td>Fienga et al. (2009)</td></tr><tr><td>8</td><td>Flora</td><td>0.426 × 10<sup>-11</sup></td><td>10.563</td><td>Baer et al. (2008b)</td></tr><tr><td>9</td><td>Metis</td><td>0.103 × 10<sup>-10</sup></td><td>23.301</td><td>Baer et al. (2008a)</td></tr><tr><td>10</td><td>Hygiea</td><td>0.454 × 10<sup>-10</sup></td><td>3.304</td><td>Baer et al. (2008a)</td></tr><tr><td>11</td><td>Parthenope</td><td>0.316 × 10<sup>-11</sup></td><td>1.899</td><td>Baer et al. (2008a)</td></tr><tr><td>13</td><td>Egeria</td><td>0.820 × 10<sup>-11</sup></td><td>19.512</td><td>Baer et al. (2008b)</td></tr><tr><td>14</td><td>Irene</td><td>0.413 × 10<sup>-11</sup></td><td>17.676</td><td>Baer et al. (2008b)</td></tr><tr><td>15</td><td>Eunomia</td><td>0.168 × 10<sup>-10</sup></td><td>4.762</td><td>Baer et al. (2008a)</td></tr><tr><td>16</td><td>Psyche</td><td>0.129 × 10<sup>-10</sup></td><td>13.178</td><td>Baer et al. (2008a)</td></tr><tr><td>17</td><td>Thetis</td><td>0.617 × 10<sup>-12</sup></td><td>10.373</td><td>Baer et al. (2008a)</td></tr><tr><td>18</td><td>Melpomene</td><td>0.151 × 10<sup>-11</sup></td><td>33.775</td><td>Baer et al. (2008b)</td></tr><tr><td>19</td><td>Fortuna</td><td>0.541 × 10<sup>-11</sup></td><td>14.048</td><td>Baer et al. (2008a)</td></tr><tr><td>20</td><td>Massalia</td><td>0.242 × 10<sup>-11</sup></td><td>16.942</td><td>Bange (1998)</td></tr><tr><td>21</td><td>Lutetia</td><td>0.129 × 10<sup>-11</sup></td><td>9.302</td><td>Baer et al. (2008b)</td></tr><tr><td>22</td><td>Kalliope</td><td>0.407 × 10<sup>-11</sup></td><td>2.457</td><td>Marchis et al. (2008b)</td></tr><tr><td>24</td><td>Themis</td><td>0.567 × 10<sup>-11</sup></td><td>37.919</td><td>Baer et al. (2008b)</td></tr><tr><td>28</td><td>Bellona</td><td>0.700 × 10<sup>-11</sup></td><td>28.571</td><td>Ivantsov (2008)</td></tr><tr><td>29</td><td>Amphitrite</td><td>0.100 × 10<sup>-10</sup></td><td>35.000</td><td>Baer et al. (2008a)</td></tr><tr><td>31</td><td>Euphrosyne</td><td>0.313 × 10<sup>-11</sup></td><td>18.850</td><td>Baer et al. (2008b)</td></tr><tr><td>45</td><td>Eugenia</td><td>0.900 × 10<sup>-11</sup></td><td>33.333</td><td>Ivantsov. (2008)</td></tr><tr><td>46</td><td>Hestia</td><td>0.109 × 10<sup>-10</sup></td><td>62.385</td><td>Bange (1998)</td></tr><tr><td>47</td><td>Aglaja</td><td>0.109 × 10<sup>-11</sup></td><td>39.450</td><td>Baer et al. (2008b)</td></tr><tr><td>48</td><td>Doris</td><td>0.610 × 10<sup>-11</sup></td><td>49.180</td><td>Kochetova (2004)</td></tr><tr><td>49</td><td>Pales</td><td>0.135 × 10<sup>-11</sup></td><td>18.519</td><td>Baer et al. (2008a)</td></tr><tr><td>52</td><td>Europa</td><td>0.976 × 10<sup>-11</sup></td><td>22.643</td><td>Baer et al. (2008a)</td></tr><tr><td>65</td><td>Cybele</td><td>0.759 × 10<sup>-11</sup></td><td>23.979</td><td>Baer et al. (2008a)</td></tr><tr><td>87</td><td>Sylvia</td><td>0.743 × 10<sup>-11</sup></td><td>0.404</td><td>Baer et al. (2008a)</td></tr><tr><td>88</td><td>Thisbe</td><td>0.572 × 10<sup>-11</sup></td><td>29.196</td><td>Baer et al. (2008a)</td></tr><tr><td>90</td><td>Antiope</td><td>0.414 × 10<sup>-12</sup></td><td>1.208</td><td>Baer et al. (2008a)</td></tr><tr><td>107</td><td>Camilla</td><td>0.563 × 10<sup>-11</sup></td><td>2.664</td><td>Marchis et al. (2008b)</td></tr><tr><td>111</td><td>Ate</td><td>0.840 × 10<sup>-10</sup></td><td>22.619</td><td>Ivantsov (2008)</td></tr><tr><td>121</td><td>Hermione</td><td>0.272 × 10<sup>-11</sup></td><td>5.515</td><td>Marchis et al. (2005)</td></tr><tr><td>130</td><td>Elektra</td><td>0.332 × 10<sup>-11</sup></td><td>6.024</td><td>Marchis et al. (2008a)</td></tr><tr><td>165</td><td>Loreley</td><td>0.160 × 10<sup>-10</sup></td><td>62.500</td><td>Aslan et al. (2007)</td></tr><tr><td>189</td><td>Phthia</td><td>0.187 × 10<sup>-13</sup></td><td>34.225</td><td>Baer et al. (2008a)</td></tr><tr><td>243</td><td>Ida</td><td>0.192 × 10<sup>-13</sup></td><td>4.688</td><td>Baer et al. (2008a)</td></tr><tr><td>253</td><td>Mathilde</td><td>0.519 × 10<sup>-13</sup></td><td>4.239</td><td>Yeomans et al. (1997)</td></tr><tr><td>283</td><td>Emma</td><td>0.694 × 10<sup>-12</sup></td><td>2.162</td><td>Marchis et al. (2008a)</td></tr></tbody></table>

low frequency; and second the long periodic terms in the disturbing function with arguments as a linear combination of the sole $\tilde{\omega}$, $\tilde{\omega}'$, $\tilde{\Omega}$ and $\tilde{\Omega}'$ not shown in (7) and (8).

The planetary orbital perturbations due to the asteroids were also discussed by Williams (1984) and Simon (1987), focusing exclusively on the planet Mars. Simon (1987) describes the analytical construction of programs to calculate the first and second derivatives of the Lagrange equations, and to express the theoretical perturbations due to an asteroid on the six ecliptical elements of a planet $\Delta\sigma_k$, where the variables $\sigma_{k=1,...,6}$ are: $a, \lambda$, the mean longitude (written as $\lambda = \lambda_0 + nt, n$ being the mean motion), $k = e \cos \varpi, h = e \sin \varpi, q = \sin \frac{1}{2} \cos \Omega, p = \sin \frac{1}{2} \sin \Omega$. The method is based on a harmonic decomposition with respect to the mean longitudes of the planet $\bar{\lambda}$ and the asteroid $\bar{\lambda}'$. Finally the perturbations can be written as

$$ \Delta\sigma_k = \sigma_k^{\text{sec}} t + \sum_{i,j} a_{ij}^k \cos(i\bar{\lambda} + j\bar{\lambda}') + b_{ij}^k \sin(i\bar{\lambda} + j\bar{\lambda}'), \quad (9) $$
as an extension of this work, Mouret et al. (2009) presented in a detailed manner the application of this decomposition for each perturbation of some selected asteroids on the Mars semi-major axis and longitude which were also the only two orbital elements studied by Williams (1984).

## 3.2. Methodology of studying asteroidal effects

For the studies mentioned above, the corresponding analytical developments of the perturbing function $\Delta R$ expressed in (6)–(8) from the orbital elements of the two bodies considered (the planet and the perturbing asteroid) are particularly cumbersome and complex, all the more so if we aim to include very long-periodic components, which are difficult to determine analytically, and could be treated as polynomial expansions. But it is possible to evaluate the perturbations in a straightforward manner from numerical integration. Therefore we computed the classical nine-body problem (the Sun and the eight planets), and then carried out the same integration with the sole addition of the perturbing asteroid (Ceres, Pallas, etc.). A simple subtraction of the two corresponding signals of one orbital parameter of the planet gives a clear and direct determination of the effect of the asteroid on this parameter.

As mentioned above, the objective of this work is to try to assess the specific influence of each of the 43 big selected asteroids on the motion of the terrestrial planets. For that goal we focused our efforts in the following complementary steps:

- we evaluated and then removed the quadratic variation of the differential signal, to isolate the periodic components;

- we performed the frequency analysis of the signal, using fast Fourier transform to determine the leading frequencies;

- we carried out a nonlinear regression in which the differential signal is modeled by the least-square method following an expression combining Fourier-type and Poisson-type components in the form

$$F(t) = \sum_{i=1}^{N} F_i(t) = \sum_{i=1}^{N} A_i \sin f_i t + B_i \cos f_i t + C_i t \sin f_i t + D_i t \cos f_i t. \quad (10)$$

We evaluated the individual influence of the gravitational action of each asteroid listed in Table 1 on the six orbital elements of each of the four terrestrial planets. Because of the number of asteroids considered (43), of orbital elements (6) and of the terrestrial planets (4), we obtained atotal of $43 \times 6 \times 4 = 1032$ curves, which are gathered together with the related tables of components $A_i$, $B_i$, $C_i$, $D_i$ in a database called ASETEP (asteroid effects on the terrestrial planets), as detailed in Sect. 7. In addition, we also calculated and included in this database the influences of each asteroid on the distance from the EMB to a given terrestrial planet (Mercury, Venus or Mars) and its orientation vector with respect to the EMB. This is of particular interest for the specific study of geocentric observations of the planets.

Because of the very huge number of results above, our purpose in the next section is only to present a few characteristic examples, and some generalizing tables, which aim to understand both the quality of our determinations and the impact of the gravitational effects of each asteroid on the orbital parameters of the terrestrial planets. Although these effects are included in the construction of modern ephemeris of planets, the individual signature of each asteroid is generally not explicitly given, and the effect of mass uncertainty is not fully characterized. One

<figure><img src="image_2.png" /><figcaption>Fig. 1. Effects of Ceres on $a_{Mars}$ for a 100 y time interval from JJ 2009.81; raw signal (black), fitting curve (red) and residuals (in green).</figcaption></figure>

of the interests of our study is to give complete empirical tables of the asteroidal effects on the orbital elements of a given planet, which could be easily re-adjusted (by simple proportionality) after revision of the asteroidal masses. Then we can classify the asteroidal effects by the range of amplitude or/and frequency.

## 4. Perturbation on the Mars orbital parameters

By using the methodology detailed in the previous section, we computed at first the effects of the asteroids on Mars orbital parameters with the help of a Runge-Kutta twelfth order integration, with a 0.01-day stepsize. The rectangular coordinates obtained with respect to the ecliptic and equinox of J2000.0 are directly converted into orbital parameters with respect to this reference plane and point.

### 4.1. Perturbation on $a_{Mars}$ and $\lambda_{Mars}$

As a basic and typical example of our computations, and to compare them with other authors, such as Mouret et al. (2009), we show in Figs. 1 and 2 the influence of the gravitational action of Ceres on the Mars semi-major axis and longitude. The time interval here is 1000 y starting from 2009 Sept. 26, i.e. between JD 2009.81 and JD 3009.81. For the clarity we plotted the figures for a 100 y time span. According to (10), we fitted the curves as a sum of Fourier (with coefficients $A_i$ and $B_i$) and Poisson components ($C_i$ and $D_i$).

In both curves, with peak-to-peak amplitudes of 400 m and 0.006" for $a_{Mars}$ and $\lambda_{Mars}$, respectively, we observe the good quality of the fit resulting from the least-squares-analysis: the fitting curve can hardly be distinguished from the original one. The flat curve represents the residuals obtained after subtracting a combination of 27 sinusoids for the semi-major axis from the original signal, and 10 sinusoids for the longitude. These sinusoids are listed in Tables 2 and 3 in the typical ASETEP database format. Notice that similar tables have been generated automatically for the 1032 cases mentioned above.

$t$ is given in Julian year: its value is $t = 2009.81$ at the origin and not $t = 0$, so that the value of the Poisson coefficients contribution can be quite high all along the interval of time explored. Nevertheless, this high value is fully compensated by the high value of some Fourier components. In other words, to best fit the signal according to (10) requires introducing two contributions that can be very high individually compared with the signal itself but which give the best fit when mixed together. Therefore
<table><caption>**Table 2.** Coefficients of Fourier and Poisson series for the effects of Ceres on $a_{\text{Mars}}$.</caption><thead><tr><th>Period (Day)</th><th>Period (Year)</th><th>SIN (m)</th><th>COS (m)</th><th>(m)</th><th>T SIN ×10<sup>-4</sup> (m)</th><th>T COS ×10<sup>-4</sup> (m)</th><th>×10<sup>-3</sup> (m)</th><th>Amplitude ×10<sup>1</sup> (m)</th><th>Asteroid</th></tr></thead><tbody><tr><td>27960.94</td><td>76.55</td><td>2.644</td><td>2.158</td><td>3.413</td><td>0.291</td><td>-3.054</td><td>0.432</td><td>0.563</td><td>1_Ceres</td></tr><tr><td>16191.36</td><td>44.33</td><td>-54.282</td><td>-53.165</td><td>75.981</td><td>67.182</td><td>62.329</td><td>8.815</td><td>5.757</td><td>1_Ceres</td></tr><tr><td>4886.42</td><td>13.38</td><td>-1.227</td><td>9.966</td><td>10.042</td><td>2.164</td><td>-13.844</td><td>1.958</td><td>0.691</td><td>1_Ceres</td></tr><tr><td>3754.72</td><td>10.28</td><td>15.277</td><td>57.810</td><td>59.794</td><td>-9.105</td><td>-0.557</td><td>0.079</td><td>5.954</td><td>1_Ceres</td></tr><tr><td>3047.01</td><td>8.34</td><td>-16.680</td><td>-0.495</td><td>16.687</td><td>7.171</td><td>0.221</td><td>0.031</td><td>1.542</td><td>1_Ceres</td></tr><tr><td>1877.01</td><td>5.14</td><td>24.877</td><td>16.008</td><td>29.582</td><td>-20.336</td><td>-10.365</td><td>1.466</td><td>2.495</td><td>1_Ceres</td></tr><tr><td>1251.31</td><td>3.43</td><td>8.691</td><td>-0.233</td><td>8.694</td><td>-10.205</td><td>0.628</td><td>0.089</td><td>0.671</td><td>1_Ceres</td></tr><tr><td>1161.51</td><td>3.18</td><td>6.924</td><td>-12.634</td><td>14.407</td><td>6.635</td><td>2.299</td><td>0.325</td><td>1.484</td><td>1_Ceres</td></tr><tr><td>887.06</td><td>2.43</td><td>3.562</td><td>-39.634</td><td>39.793</td><td>14.293</td><td>19.194</td><td>2.714</td><td>3.622</td><td>1_Ceres</td></tr><tr><td>717.51</td><td>1.96</td><td>-5.645</td><td>-11.862</td><td>13.136</td><td>9.587</td><td>5.867</td><td>0.830</td><td>1.149</td><td>1_Ceres</td></tr><tr><td>687.06</td><td>1.88</td><td>-3.017</td><td>-0.612</td><td>3.078</td><td>1.168</td><td>-1.459</td><td>0.206</td><td>0.304</td><td>1_Ceres</td></tr><tr><td>602.38</td><td>1.65</td><td>-4.163</td><td>-2.346</td><td>4.779</td><td>6.308</td><td>0.588</td><td>0.083</td><td>0.316</td><td>1_Ceres</td></tr><tr><td>580.76</td><td>1.59</td><td>-43.229</td><td>24.966</td><td>49.921</td><td>-13.036</td><td>-34.487</td><td>4.877</td><td>4.935</td><td>1_Ceres</td></tr><tr><td>502.95</td><td>1.38</td><td>-11.922</td><td>15.858</td><td>19.840</td><td>-7.375</td><td>-15.356</td><td>2.172</td><td>1.844</td><td>1_Ceres</td></tr><tr><td>443.53</td><td>1.21</td><td>-0.836</td><td>8.157</td><td>8.199</td><td>-5.500</td><td>-8.805</td><td>1.245</td><td>0.669</td><td>1_Ceres</td></tr><tr><td>408.23</td><td>1.12</td><td>11.695</td><td>-37.093</td><td>38.893</td><td>-56.454</td><td>184.861</td><td>26.143</td><td>1.956</td><td>1_Ceres</td></tr><tr><td>396.66</td><td>1.09</td><td>1.263</td><td>2.717</td><td>2.996</td><td>-3.399</td><td>-2.282</td><td>0.323</td><td>0.215</td><td>1_Ceres</td></tr><tr><td>387.17</td><td>1.06</td><td>19.798</td><td>0.132</td><td>19.798</td><td>-3.570</td><td>24.294</td><td>3.436</td><td>2.008</td><td>1_Ceres</td></tr><tr><td>350.98</td><td>0.96</td><td>9.686</td><td>-4.009</td><td>10.483</td><td>-0.729</td><td>12.473</td><td>1.764</td><td>0.962</td><td>1_Ceres</td></tr><tr><td>320.97</td><td>0.88</td><td>3.050</td><td>-3.512</td><td>4.652</td><td>0.594</td><td>5.911</td><td>0.836</td><td>0.402</td><td>1_Ceres</td></tr><tr><td>314.71</td><td>0.86</td><td>4.284</td><td>2.959</td><td>5.207</td><td>-5.783</td><td>3.139</td><td>0.444</td><td>0.473</td><td>1_Ceres</td></tr><tr><td>290.37</td><td>0.79</td><td>-7.373</td><td>-1.951</td><td>7.626</td><td>5.054</td><td>-12.518</td><td>1.770</td><td>0.819</td><td>1_Ceres</td></tr><tr><td>272.16</td><td>0.75</td><td>13.490</td><td>11.265</td><td>17.575</td><td>-66.768</td><td>-54.343</td><td>7.685</td><td>0.852</td><td>1_Ceres</td></tr></tbody><tfoot><tr><td colspan="10">AMPLITUDE Before = 409.294276 After = 57.641405</td></tr></tfoot></table>

**Notes.** Bias: $0.137559 \times 10^{+2}$; linear: $-0.377681 \times 10^{-2}$; $T^2$: $0.333290 \times 10^{-5}$. The next-to-last column indicates the peak-to-peak amplitude of the total mixed component for the 1000 y time interval considered.

<table><caption>**Table 3.** Coefficients of Fourier and Poisson series for the effects of Ceres on $\lambda_{\text{Mars}}$.</caption><thead><tr><th>period (Day)</th><th>period (Year)</th><th>SIN ×10<sup>-5</sup> (")</th><th>COS ×10<sup>-5</sup> (")</th><th>×10<sup>-4</sup> (")</th><th>T SIN ×10<sup>-8</sup> (")</th><th>T COS ×10<sup>-8</sup> (")</th><th>×10<sup>-8</sup> (")</th><th>Amplitude ×10<sup>-2</sup> (")</th><th>Asteroid</th></tr></thead><tbody><tr><td>27961.03</td><td>76.55</td><td>-10.158</td><td>13.521</td><td>1.691</td><td>0.841</td><td>0.666</td><td>0.942</td><td>0.021</td><td>1_Ceres</td></tr><tr><td>16190.53</td><td>44.33</td><td>149.608</td><td>-165.221</td><td>22.289</td><td>-14.986</td><td>19.171</td><td>27.112</td><td>0.177</td><td>1_Ceres</td></tr><tr><td>15442.23</td><td>42.28</td><td>-8.026</td><td>-1.417</td><td>0.815</td><td>2.600</td><td>0.650</td><td>0.919</td><td>0.001</td><td>1_Ceres</td></tr><tr><td>4886.18</td><td>13.38</td><td>-11.271</td><td>-0.854</td><td>1.130</td><td>1.674</td><td>0.200</td><td>0.282</td><td>0.007</td><td>1_Ceres</td></tr><tr><td>3754.68</td><td>10.28</td><td>-56.393</td><td>15.559</td><td>5.850</td><td>0.449</td><td>-0.867</td><td>1.227</td><td>0.057</td><td>1_Ceres</td></tr><tr><td>3046.73</td><td>8.34</td><td>-0.633</td><td>-5.797</td><td>0.583</td><td>-0.101</td><td>0.045</td><td>0.064</td><td>0.006</td><td>1_Ceres</td></tr><tr><td>1877.01</td><td>5.14</td><td>-7.474</td><td>12.606</td><td>1.465</td><td>0.362</td><td>-1.003</td><td>1.418</td><td>0.013</td><td>1_Ceres</td></tr><tr><td>1161.51</td><td>3.18</td><td>7.066</td><td>3.973</td><td>0.811</td><td>-0.223</td><td>0.274</td><td>0.387</td><td>0.008</td><td>1_Ceres</td></tr><tr><td>887.06</td><td>2.43</td><td>11.749</td><td>1.247</td><td>1.182</td><td>-0.559</td><td>0.397</td><td>0.562</td><td>0.011</td><td>1_Ceres</td></tr><tr><td>580.76</td><td>1.59</td><td>-5.492</td><td>-9.687</td><td>1.114</td><td>0.781</td><td>-0.210</td><td>0.297</td><td>0.011</td><td>1_Ceres</td></tr></tbody><tfoot><tr><td colspan="10">AMPLITUDE Before = 0.005746 After = 0.000604</td></tr></tfoot></table>

**Notes.** Bias: $0.221687$; linear: $-0.550293 \times 10^{-4}$; $T^2$: $-0.278178 \times 10^{-7}$. The next-to-last column indicates the peak-to-peak amplitude of the total mixed component for the 1000 y time interval considered.

we should not refer directly to the values of $A_i$, $B_i$, $C_i$ and $D_i$ to evaluate the amplitude of a component at a given frequency $f_i$ in (10) but to the peak-to-peak amplitude of $F_i(t)$ during the interval of time considered. This peak-to-peak amplitude is given in the next-to-last columns of Tables 2 and 3, where each component has been ranged in decreasing order of period. For information and to avoid too many data in the paper, the results of the effects of Ceres on $e_{\text{Mars}}$, $i_{\text{Mars}}$, $\Omega_{\text{Mars}}$ and $\varpi_{\text{Mars}}$ are given in Figs. A.1 to A.4 and Tables A.1 to A.4.

For comparison we also show in Figs. 3 and 4 the effects of Pallas and Vesta on the Mars semi-major axis for the same time interval as in Fig. 1. The peak-to-peak amplitudes reach roughly 200 m and 400 m, this last value is roughly the same as for Ceres, but the nature of the three signals is quite different

from each other. For the sake of clarity the tables characterizing the adjustment are presented in Tables A.5 and A.6.

## 4.2. Comparison with the analytical results in Mouret et al. (2009)

It is interesting and challenging to compare our numerical results with the analytical ones found by Mouret et al. (2009) when using developments as in (7) and (8). For that purpose it is necessary to choose the same time interval as these authors, that is to say 1000 y, because the developments in series are highly dependent on this time interval. In Tables 4 and 5 we show this comparison for the variation of the Mars semi-major axis and mean longitude due to the gravitational influence of Ceres, Pallas and
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Effects of Ceres on $\lambda_{\text{Mars}}$ for a 100 y time interval from JJ 2009.81; raw signal (in black), fitting curve (in red) and residuals (in green).</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Effects of Pallas on $a_{\text{Mars}}$ for a 100 y time interval from JJ 2009.81; raw signal (in black), fitting curve (in red) and residuals (in green).</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Effects of Vesta on $a_{\text{Mars}}$ for a 100 y time interval from JJ 2009.81; raw signal (in black), fitting curve (in red) and residuals (in green).</figcaption></figure>

Vesta. We can notice the remarkable general agreement between our numerical determination of the periods and amplitudes of the sinusoidal components that characterize the signal, and the analytical terms developed by Mouret et al. (2009), whose the arguments are a linear combination of the mean longitudes of Mars and the asteroid.

For the effect of Ceres on $a_{\text{Mars}}$ (Table 4) a small but significant discrepancy arises only for the leading term with period 44.33 y and amplitude 57.57 m, which seems to correspond to the component with the period 45.77 y and 60.59 m amplitude in the above paper. This discrepancy probably arises because we used a quadratic regression and removed it from the signal to isolate the periodic parts, which should affect the determination of low-frequency terms. We also notice components in our calculations that are not mentioned in the paper above, with significant amplitudes 19.56 m, 8.52 m, 5.63 m and 4.78 m and periods 1.12 y, 0.75 y, 76.55 y and 0.56 y. the reason for this could be that analytical developments, even at a low order of masses and eccentricities, are cumbersome with a mixing of many sinusoidal contributions involving the mean longitudes. Therefore some components may be missing.

For the effect of Ceres on $\lambda_{\text{Mars}}$, the comparison between our results and Mouret et al. (2009) still agrees very well for the majority of components both in amplitudes and periods. Differences similar to those noted for $a_{\text{Mars}}$ concern the leading term with period 44.33 y instead of 45.77 y and amplitude 1.77 mas (milliarcsecond).

The agreement is almost perfect for the influences of Pallas and Vesta, if we exclude large amplitudes without counterparts in the analytical developments, such as the component with amplitude 19.02 m and period 1.12 y for $a_{\text{Mars}}$ and the second-largest component with amplitude 0.333 mas and period 59.79 y for $\lambda_{\text{Mars}}$.

### 4.3. Perturbations on the other orbital parameters *e*, *i*, *Ω*, $\tilde{\omega}$

The gravitational effects of Ceres and other leading asteroids on Mars orbital parameters $e_{\text{Mars}}$, $i_{\text{Mars}}$, $\Omega_{\text{Mars}}$, $\tilde{\omega}_{\text{Mars}}$ have not been explicitly given either by Williams (1984) or by Mouret et al. (2009) although their relative amplitude is of the same order (and sometimes even larger) than the corresponding effects on $a_{\text{Mars}}$ and $\lambda_{\text{Mars}}$ described above.

Two important remarks about the curves plotted for Ceres, Pallas and Vesta in Figs. 5 to 8 need to be made. First of all, in contrast to what happens for the semi-major axes (see Fig. 1), long periodic effects that can be characterized by linear or quadratic expressions appear clearly, in particular for $i_{\text{Mars}}$ (Fig. 5) and $\Omega_{\text{Mars}}$ (Fig. 7), for which they dominate the periodic oscillations, even during the short time span considered here (100 y). Second, we can observe that although its mass is considerably smaller than that of Ceres (with a ratio roughly 1/3), the effect of Vesta on the Mars orbital elements can be significantly stronger, as is the case for $e_{\text{Mars}}$ (Fig. 6) and $\tilde{\omega}_{\text{Mars}}$ (Fig. 8). For the three asteroids considered here, the peak-to-peak amplitudes are on the order of $10^{-8}$ for $e_{\text{Mars}}$, 1 mas (milliarcsecond) for $i_{\text{Mars}}$, $0^{\circ}.05$ for $\Omega_{\text{Mars}}$ and $0^{\circ}.01$ for $\tilde{\omega}_{\text{Mars}}$, in the 100 y time interval considered.

In Figs. A.1–A.4 we show the effects of Ceres after the quadratic regression is removed. For each of the parameters, leading sinusoidal signatures appear clearly and the same fitting method as previously can be applied. Results of the decomposition are shown in Tables A.1–A.4.

### 4.4. Amplitude of the perturbations for each parameter and each asteroid

As we already mentioned in Sect. 3.2, we extended our study to the 1032 cases (43 asteroids × 6 planets × 6 orbital elements). The output of our computations is presented in Table 6: we
<table><caption>**Table 4.** Comparisons of the amplitude of the sinusoids found in this paper with the analytical ones of Mouret et al. (2009) of the effects of Ceres, Pallas and Vesta on $a_{Mars}$.</caption><thead><tr><th colspan="2">Our analysis</th><th colspan="4">Mouret et al. (2009)</th></tr><tr><th>Period (y)</th><th>Amplitude (m)</th><th>Period (y)</th><th>Amplitude (m)</th><th>Argument</th><th>iλ<sub>M</sub> + jλ<sub>c</sub></th></tr></thead><tbody><tr><td colspan="6" style="text-align: center;">Ceres</td></tr><tr><td>10.280</td><td>59.5370</td><td>10.25</td><td>59.390</td><td>1</td><td>-2</td></tr><tr><td>44.330</td><td>57.5675</td><td>45.77</td><td>60.587</td><td></td><td></td></tr><tr><td>1.590</td><td>49.3516</td><td>1.59</td><td>49.517</td><td>2</td><td>-2</td></tr><tr><td>2.429</td><td>36.2161</td><td>2.43</td><td>36.352</td><td>2</td><td>-3</td></tr><tr><td>5.139</td><td>24.9480</td><td>5.12</td><td>25.132</td><td>2</td><td>-4</td></tr><tr><td>1.118</td><td>19.5619</td><td></td><td></td><td></td><td></td></tr><tr><td>1.060</td><td>20.0837</td><td>1.06</td><td>19.897</td><td>3</td><td>-3</td></tr><tr><td>1.377</td><td>18.4410</td><td>1.38</td><td>18.550</td><td>3</td><td>-4</td></tr><tr><td>8.342</td><td>15.4159</td><td>8.37</td><td>15.708</td><td>1</td><td>-3</td></tr><tr><td>3.180</td><td>14.8437</td><td>3.18</td><td>14.900</td><td>1</td><td>-1</td></tr><tr><td>1.964</td><td>11.4949</td><td>1.96</td><td>11.369</td><td>3</td><td>-5</td></tr><tr><td>0.961</td><td>9.6248</td><td>0.96</td><td>9.694</td><td>4</td><td>-5</td></tr><tr><td>0.745</td><td>8.5178</td><td></td><td></td><td></td><td></td></tr><tr><td>0.795</td><td>8.1945</td><td>0.79</td><td>7.944</td><td>4</td><td>-4</td></tr><tr><td>13.378</td><td>6.9088</td><td>13.20</td><td>7.465</td><td>3</td><td>-7</td></tr><tr><td>3.426</td><td>6.7110</td><td>3.42</td><td>6.672</td><td>3</td><td>-6</td></tr><tr><td>1.214</td><td>6.6896</td><td>1.21</td><td>6.627</td><td>4</td><td>-6</td></tr><tr><td>76.553</td><td>5.6321</td><td></td><td></td><td></td><td></td></tr><tr><td>0.738</td><td>4.8935</td><td>0.74</td><td>4.952</td><td>5</td><td>-6</td></tr><tr><td>0.862</td><td>4.7347</td><td>0.86</td><td>4.682</td><td>3</td><td>-2</td></tr><tr><td>0.559</td><td>4.7800</td><td></td><td></td><td></td><td></td></tr><tr><td>0.879</td><td>4.0163</td><td>0.88</td><td>3.934</td><td>5</td><td>-7</td></tr><tr><td>1.649</td><td>3.1624</td><td>1.65</td><td>3.665</td><td>4</td><td>-7</td></tr><tr><td>0.678</td><td>3.4828</td><td>0.68</td><td>3.366</td><td>4</td><td>-3</td></tr><tr><td>1.881</td><td>3.0392</td><td>1.88</td><td>2.992</td><td>1</td><td>0</td></tr><tr><td colspan="6" style="text-align: center;">Pallas</td></tr><tr><td>49.058</td><td>73.5324</td><td>48.65</td><td>74.350</td><td>2</td><td>-5</td></tr><tr><td>1.118</td><td>19.0229</td><td></td><td></td><td></td><td></td></tr><tr><td>12.855</td><td>14.7998</td><td>12.89</td><td>14.750</td><td>3</td><td>-7</td></tr><tr><td>24.518</td><td>14.1215</td><td>24.32</td><td>13.883</td><td>4</td><td>-10</td></tr><tr><td>5.093</td><td>13.0090</td><td>5.10</td><td>13.224</td><td>2</td><td>-4</td></tr><tr><td>2.421</td><td>9.6624</td><td>2.42</td><td>9.619</td><td>2</td><td>-3</td></tr><tr><td>17.424</td><td>6.2069</td><td>17.54</td><td>7.121</td><td>5</td><td>-12</td></tr><tr><td>1.588</td><td>7.1841</td><td>1.59</td><td>6.537</td><td>2</td><td>-2</td></tr><tr><td>10.185</td><td>6.8003</td><td>10.19</td><td>6.508</td><td>1</td><td>-2</td></tr><tr><td>8.433</td><td>6.9895</td><td>8.43</td><td>6.358</td><td>1</td><td>-3</td></tr><tr><td>3.395</td><td>5.7901</td><td>3.40</td><td>5.789</td><td>3</td><td>-6</td></tr><tr><td>7.196</td><td>5.5975</td><td>7.18</td><td>5.505</td><td>3</td><td>-8</td></tr><tr><td>1.956</td><td>4.6407</td><td>1.96</td><td>4.638</td><td>3</td><td>-5</td></tr><tr><td>5.683</td><td>4.8822</td><td>5.69</td><td>4.593</td><td>4</td><td>-9</td></tr><tr><td>27.009</td><td>4.4647</td><td>27.44</td><td>3.964</td><td>7</td><td>-17</td></tr><tr><td>4.217</td><td>3.8491</td><td>4.21</td><td>3.934</td><td>2</td><td>-6</td></tr><tr><td>1.374</td><td>3.6767</td><td>1.37</td><td>3.590</td><td>3</td><td>-4</td></tr><tr><td>16.358</td><td>2.8954</td><td>16.22</td><td>3.366</td><td>6</td><td>-15</td></tr><tr><td>2.546</td><td>2.9761</td><td>2.55</td><td>2.857</td><td>4</td><td>-8</td></tr><tr><td>2.546</td><td>2.9761</td><td>1.64</td><td>2.423</td><td>4</td><td>-7</td></tr><tr><td>1.058</td><td>2.2976</td><td></td><td></td><td></td><td></td></tr><tr><td>1.210</td><td>2.1458</td><td>1.21</td><td>2.035</td><td>4</td><td>-6</td></tr><tr><td colspan="6" style="text-align: center;">Vesta</td></tr><tr><td>25.855</td><td>38.9394</td><td>25.93</td><td>38.746</td><td>2</td><td>-4</td></tr><tr><td>1.952</td><td>30.6734</td><td>1.95</td><td>30.368</td><td>2</td><td>-2</td></tr><tr><td>4.223</td><td>25.5473</td><td>4.22</td><td>24.983</td><td>2</td><td>-3</td></tr><tr><td>1.301</td><td>15.7973</td><td>1.30</td><td>15.708</td><td>3</td><td>-3</td></tr><tr><td>2.029</td><td>13.9884</td><td>2.03</td><td>13.254</td><td>3</td><td>-4</td></tr><tr><td>3.904</td><td>11.5442</td><td>3.90</td><td>11.325</td><td>1</td><td>-1</td></tr><tr><td>4.598</td><td>10.5150</td><td>4.59</td><td>9.829</td><td>3</td><td>-5</td></tr><tr><td>17.238</td><td>9.9542</td><td>17.29</td><td>9.140</td><td>3</td><td>-6</td></tr><tr><td>0.976</td><td>8.5505</td><td>0.98</td><td>8.198</td><td>4</td><td>-4</td></tr><tr><td>1.335</td><td>8.3356</td><td>1.33</td><td>8.078</td><td>4</td><td>-5</td></tr></tbody></table>

**Notes.** Coefficients are ranged in decreasing order of amplitude.

<table><caption>**Table 5.** Comparisons of the sinusoids found in this paper with those in Mouret et al. (2009) of the effects of Ceres, Pallas and Vesta on $\lambda_{Mars}$.</caption><thead><tr><th colspan="2">Our analysis</th><th colspan="4">Mouret et al. (2009)</th></tr><tr><th>Period (y)</th><th>Amplitude ×10<sup>-3</sup>"</th><th>Period (y)</th><th>Amplitude ×10<sup>-3</sup>"</th><th>Argument</th><th>iλ<sub>M</sub> + jλ<sub>c</sub></th></tr></thead><tbody><tr><td colspan="6" style="text-align: center;">Ceres</td></tr><tr><td>10.280</td><td>0.573</td><td>10.25</td><td>0.571</td><td>1</td><td>-2</td></tr><tr><td>44.330</td><td>1.768</td><td>45.77</td><td>1.914</td><td>2</td><td>-5</td></tr><tr><td>76.553</td><td>0.205</td><td></td><td></td><td></td><td></td></tr><tr><td>1.590</td><td>0.108</td><td>1.59</td><td>0.127</td><td>2</td><td>-4</td></tr><tr><td>2.429</td><td>0.108</td><td>2.43</td><td>0.109</td><td>2</td><td>-2</td></tr><tr><td>5.139</td><td>0.126</td><td>5.12</td><td>0.108</td><td>2</td><td>-3</td></tr><tr><td>8.342</td><td>0.058</td><td>8.37</td><td>0.060</td><td>1</td><td>-3</td></tr><tr><td>3.180</td><td>0.080</td><td>3.18</td><td>0.039</td><td>3</td><td>-4</td></tr><tr><td>13.378</td><td>0.069</td><td>13.20</td><td>0.025</td><td>3</td><td>-6</td></tr><tr><td colspan="6" style="text-align: center;">Pallas</td></tr><tr><td>49.058</td><td>2.521</td><td>48.65</td><td>2.537</td><td>2</td><td>-5</td></tr><tr><td>59.790</td><td>0.333</td><td></td><td></td><td></td><td></td></tr><tr><td>24.517</td><td>0.209</td><td>24.32</td><td>0.227</td><td>4</td><td>-10</td></tr><tr><td>12.855</td><td>0.148</td><td>12.89</td><td>0.153</td><td>3</td><td>-7</td></tr><tr><td>17.424</td><td>0.088</td><td>17.54</td><td>0.099</td><td>5</td><td>-12</td></tr><tr><td colspan="6" style="text-align: center;">Vesta</td></tr><tr><td>51.657</td><td>4.372</td><td></td><td></td><td></td><td></td></tr><tr><td>25.854</td><td>0.679</td><td>25.93</td><td>0.668</td><td>2</td><td>-4</td></tr><tr><td>4.223</td><td>0.112</td><td>4.22</td><td>0.110</td><td>2</td><td>-3</td></tr><tr><td>40.162</td><td>0.107</td><td></td><td></td><td></td><td></td></tr><tr><td>17.237</td><td>0.104</td><td>17.29</td><td>0.103</td><td>3</td><td>-6</td></tr><tr><td>1.118</td><td>0.096</td><td></td><td></td><td></td><td></td></tr><tr><td>15.779</td><td>0.083</td><td></td><td></td><td></td><td></td></tr><tr><td>1.952</td><td>0.078</td><td>1.95</td><td>0.078</td><td>2</td><td>-2</td></tr><tr><td>32.816</td><td>0.078</td><td></td><td></td><td></td><td></td></tr><tr><td>2.029</td><td>0.036</td><td>2.03</td><td>0.035</td><td>3</td><td>-4</td></tr><tr><td>1.301</td><td>0.032</td><td>1.30</td><td>0.031</td><td>3</td><td>-3</td></tr></tbody></table>

**Notes.** Coefficients are ranged in decreasing order of amplitude.

<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Effects of Ceres, Pallas and Vesta on $i_{Mars}$, for a 100 y time interval, starting from JD 2009.81.</figcaption></figure>

provide the peak-to-peak amplitudes of the effects of each of the 43 asteroids listed in Table 1 on each Mars orbital element. This is done for a 1000 y time interval, beginning in 2009. The amplitudes stand for the total signal including secular trends and periodic components. Moreover, to compare the amplitudes of the effects on $e_{Mars}$ to those on $a_{Mars}$ with the same scale, we multiplied these effects by the nominal value $a_0$ of the Mars semi-major axis. Thus we can remark that throughout, $a_0 \times \Delta e$ is significantly higher than $\Delta a$ with a ratio typically on the order of 2 or 3, but sometimes considerably more, as for Pallas and Vesta. The relatively large amplitude of the effects due to the asteroid
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Effects of Ceres, Pallas and Vesta on $e_{Mars}$, for a 100 y time interval starting from JD 2009.81.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Effects of Ceres, Pallas and Vesta on $\Omega_{Mars}$, for a 100 y time interval, starting from JD 2009.81.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Effects of Ceres, Pallas and Vesta on $\tilde{\omega}_{Mars}$, for a 100 y time interval, starting from JD 2009.81.</figcaption></figure>

111 Ate compared to those due to these two last asteroid is remarkable. This is because the mass of Ate is comparable to that of these two big objects. For Ceres, during the 1000 y interval of time considered, the peak-to-peak amplitudes are 416 m for $\Delta a$, 1300 m for $a_0 \times \Delta e$, 6.8 mas (millarcsec) for $\Delta i$, 0.410 for $\Delta\Omega$ and 0.197 for $\Delta\lambda$. These effects of the biggest asteroid are surpassed by those of Vesta for $\Delta a$ (546.7 m) and $a_0 \times \Delta e$ (4666 m), by Pallas and 111 Ate for $a_0 \times \Delta e$ (2025 m and 426 m).

As a direct by-product of our results, we present in Table 7 the uncertainty of each Mars orbital parameter due to the mass uncertainty of each selected asteroid. Indeed, according to (6), this uncertainty is directly dependent on the asteroid mass, it is

<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Effects of Ceres on $a_{EMB}$, for a 100 y time interval, from JD 2009.81. Residuals (in green) after removing the sinusoidal components from the least-squares analysis.</figcaption></figure>

sufficient for the purpose to multiply each number of Table 6 by the relative mass uncertainty $\Delta M/M$, whose value has been already given for each object in Table 1. Because of its relatively large mass (17.7% that of Ceres) and its very high mass uncertainty (22.6%), the asteroid 111 Ate is by far the most problematic one, causing uncertainties higher than 93 times that of Ceres for $a_{Mars}$, 197 times for $a_0\Delta e_{Mars}$, 44 times for $i_{Mars}$, 23 times for $\Omega_{Mars}$, 58 times for $\tilde{\omega}_{Mars}$ and 295 times for $\lambda_{Mars}$. Uncertainties on the same order can be observed for 6 Hebe (18.7% mass uncertainty), 7 Iris (8.6%), 9 Metis (23.3%) 19 Fortuna (14.0%), 20 Massalia (16.9%) and 46 Hestia (63%). We think that a table such as Table 7 shows in a very straightforward manner for which asteroid observational campaigns should be launched to minimize errors in the Mars ephemeris.

# 5. Perturbation on the EMB, Mercury and Venus

In this section we carry out the same kind of study as in the previous section for the terrestrial planets other than Mars. Instead of the Earth isolated, we consider the EMB. This facilitates our calculations noticeably. A specific study of the asteroidal effects on the Earth and the Moon taken individually is scheduled for a future paper.

## 5.1. Effects on EMB

As an example we show in Figs. 9–11 the effects of the gravitational action of Ceres, Pallas and Vesta on the EMB semi-major axis. They are characterized by very stable and regular leading oscillations with peak-to-peak amplitudes about 25 m for Ceres, 9 m for Pallas and 15 m for Vesta. The flat curves in each figure represent the residuals after subtracting the fit with the help of Fourier and Poisson components. We were able to reconstruct the curve of Ceres with as many as 50 frequencies. The leading components are shown in Table 8. We note that the amplitudes and periods of the leading oscillations are significantly shorter for the EMB than for Mars.

In Table 9 we present the peak-to-peak amplitudes of the effects of each asteroid on the orbital parameters of the EMB for a 1000 y time interval. They are generally, but no systematically smaller than the corresponding effect on Mars. This is because the mass of Mars is much smaller than the mass of our planet ($M_{Mars}/M_{Earth} \approx 0.107$) and its orbit is much closer to the asteroid main belt.
<table><caption>**Table 6.** Peak-to-peak amplitudes of the variations of the six Mars orbital elements for a 1000 y time interval starting from JD 2009.81.</caption><thead><tr><th rowspan="2">asteroid</th><th>Δa</th><th>Δe × a<sub>Mars</sub></th><th>Δi</th><th>ΔΩ</th><th>Δω</th><th>Δλ</th></tr><tr><th>m</th><th>m</th><th>×10<sup>-6</sup>"</th><th>×10<sup>-6</sup>"</th><th>×10<sup>-6</sup>"</th><th>×10<sup>-6</sup>"</th></tr></thead><tbody><tr><td>1 Ceres</td><td>415.976</td><td>1300.018</td><td>6764.552</td><td>409910.423</td><td>101767.790</td><td>197860.644</td></tr><tr><td>2 Pallas</td><td>212.853</td><td>2025.012</td><td>4517.624</td><td>110490.432</td><td>21951.167</td><td>143118.437</td></tr><tr><td>3 Juno</td><td>38.999</td><td>166.213</td><td>381.803</td><td>8282.676</td><td>1641.266</td><td>24833.521</td></tr><tr><td>4 Vesta</td><td>546.665</td><td>4665.960</td><td>5665.561</td><td>98770.967</td><td>50554.943</td><td>623437.866</td></tr><tr><td>6 Hebe</td><td>144.013</td><td>607.501</td><td>606.177</td><td>2122.292</td><td>6015.502</td><td>40524.178</td></tr><tr><td>7 Iris</td><td>67.551</td><td>589.408</td><td>77.461</td><td>11921.764</td><td>1707.940</td><td>3249.652</td></tr><tr><td>8 Flora</td><td>28.615</td><td>351.759</td><td>264.020</td><td>1527.953</td><td>1472.282</td><td>23161.386</td></tr><tr><td>9 Metis</td><td>68.308</td><td>608.089</td><td>96.367</td><td>10037.872</td><td>5999.291</td><td>6453.236</td></tr><tr><td>10 Hygiea</td><td>53.587</td><td>547.873</td><td>225.513</td><td>13260.778</td><td>5733.690</td><td>6918.631</td></tr><tr><td>11 Parthenope</td><td>8.210</td><td>65.473</td><td>76.711</td><td>197.072</td><td>688.702</td><td>2268.556</td></tr><tr><td>13 Egeria</td><td>14.479</td><td>57.479</td><td>77.631</td><td>15907.661</td><td>1666.410</td><td>6619.356</td></tr><tr><td>14 Irene</td><td>14.634</td><td>79.531</td><td>122.264</td><td>5390.858</td><td>1888.597</td><td>11809.004</td></tr><tr><td>15 Eunomia</td><td>41.963</td><td>416.526</td><td>593.602</td><td>16346.103</td><td>1348.323</td><td>48827.323</td></tr><tr><td>16 Psyche</td><td>10.560</td><td>135.172</td><td>101.183</td><td>2213.723</td><td>824.135</td><td>6357.709</td></tr><tr><td>17 Thetis</td><td>2.740</td><td>23.640</td><td>20.280</td><td>38.826</td><td>173.695</td><td>1944.757</td></tr><tr><td>18 Melpomene</td><td>6.828</td><td>57.461</td><td>138.647</td><td>1381.789</td><td>195.621</td><td>497.143</td></tr><tr><td>19 Fortuna</td><td>48.578</td><td>308.680</td><td>30.604</td><td>3799.559</td><td>1832.499</td><td>46311.253</td></tr><tr><td>20 Massalia</td><td>80.479</td><td>238.095</td><td>9.043</td><td>2019.272</td><td>2971.102</td><td>47763.291</td></tr><tr><td>21 Lutetia</td><td>6.698</td><td>48.954</td><td>10.871</td><td>207.539</td><td>516.701</td><td>7525.555</td></tr><tr><td>22 Kalliope</td><td>2.682</td><td>36.025</td><td>21.367</td><td>3708.021</td><td>401.749</td><td>1934.839</td></tr><tr><td>24 Themis</td><td>14.250</td><td>50.778</td><td>1.832</td><td>527.693</td><td>1237.710</td><td>16186.184</td></tr><tr><td>28 Bellona</td><td>16.140</td><td>74.403</td><td>214.068</td><td>1144.324</td><td>2820.137</td><td>7919.977</td></tr><tr><td>29 Amphitrite</td><td>17.267</td><td>165.764</td><td>224.757</td><td>2016.979</td><td>2826.971</td><td>12679.011</td></tr><tr><td>31 Euphrosyne</td><td>13.346</td><td>56.216</td><td>54.024</td><td>3639.611</td><td>482.373</td><td>6169.238</td></tr><tr><td>45 Eugenia</td><td>8.092</td><td>107.311</td><td>202.524</td><td>1652.153</td><td>1867.647</td><td>2507.753</td></tr><tr><td>46 Hestia</td><td>20.516</td><td>229.351</td><td>112.954</td><td>5941.046</td><td>1115.982</td><td>11354.487</td></tr><tr><td>47 Aglaja</td><td>0.743</td><td>4.292</td><td>10.767</td><td>84.131</td><td>29.633</td><td>59.064</td></tr><tr><td>48 Doris</td><td>4.005</td><td>28.678</td><td>58.280</td><td>2075.386</td><td>802.772</td><td>969.684</td></tr><tr><td>49 Pales</td><td>1.635</td><td>27.457</td><td>5.568</td><td>440.772</td><td>120.879</td><td>210.067</td></tr><tr><td>52 Europa</td><td>9.826</td><td>68.303</td><td>131.673</td><td>607.090</td><td>1828.279</td><td>13298.494</td></tr><tr><td>65 Cybele</td><td>3.061</td><td>33.727</td><td>37.098</td><td>520.461</td><td>303.235</td><td>2904.304</td></tr><tr><td>87 Sylvia</td><td>1.974</td><td>7.763</td><td>16.837</td><td>2297.924</td><td>229.907</td><td>419.793</td></tr><tr><td>88 Thisbe</td><td>7.671</td><td>48.380</td><td>54.682</td><td>3342.237</td><td>256.270</td><td>6600.855</td></tr><tr><td>90 Antiope</td><td>0.866</td><td>4.080</td><td>0.385</td><td>4.311</td><td>34.948</td><td>628.948</td></tr><tr><td>107 Camilla</td><td>3.072</td><td>6.013</td><td>58.238</td><td>1078.727</td><td>503.130</td><td>278.376</td></tr><tr><td>111 Ate</td><td>218.942</td><td>1425.543</td><td>1683.832</td><td>53568.836</td><td>33425.288</td><td>328096.755</td></tr><tr><td>121 Hermione</td><td>1.001</td><td>11.230</td><td>4.441</td><td>587.147</td><td>95.005</td><td>1318.865</td></tr><tr><td>130 Elektra</td><td>4.738</td><td>23.884</td><td>110.486</td><td>393.899</td><td>589.786</td><td>3224.094</td></tr><tr><td>165 Loreley</td><td>12.197</td><td>52.039</td><td>256.157</td><td>5936.699</td><td>937.680</td><td>901.986</td></tr><tr><td>189 Phthia</td><td>0.132</td><td>0.395</td><td>0.255</td><td>21.056</td><td>7.059</td><td>152.358</td></tr><tr><td>243 Ida</td><td>0.108</td><td>0.116</td><td>0.049</td><td>3.560</td><td>3.224</td><td>28.511</td></tr><tr><td>253 Mathilde</td><td>0.193</td><td>0.395</td><td>1.320</td><td>42.175</td><td>12.877</td><td>224.965</td></tr><tr><td>283 Emma</td><td>0.486</td><td>4.920</td><td>9.306</td><td>230.447</td><td>22.563</td><td>56.147</td></tr></tbody></table>

## 5.2. Effects on Mercury and Venus

The effects of the asteroids on the orbital parameters of the two internal planets Mercury and Venus are presented in Tables 10 and 11 also for a 1000 y time span. The peak-to-peak amplitudes for $a_{\text{Mercury}}$ and $a_{\text{Venus}}$ do not exceed 1 m and 10 m, i.e. roughly 400 times and 40 times less than for $a_{\text{Mars}}$. Nevertheless, the ratio is significantly lower for other parameters, also because of the predominance of long-term components, which can be characterized by quadratic or polynomial expressions in the interval of time considered.

# 6. Influences of asteroid on the EMB-Mars distance and Mars orientation

In view of the ground-based observations of Mars, the fundamental parameters are both the distance $d_{\text{Mars}}$ from the Earth geocenter (or the EMB in an alternative case considered here) to the planet geocenter, and the geocentric direction of the planet

generally given through the set of its geocentric coordinates ($\alpha_{\text{Mars}}$, $\delta_{\text{Mars}}$). Thus the gravitational effect of a given asteroid on these parameters comes directly from the combination of the influences of the asteroid on the orbital elements of the EMB and the planet, which we studied above. As for these influences, the uncertainty of the asteroid masses directly affects the accuracy of the determination of $d_{\text{Mars}}$, $\alpha_{\text{Mars}}$ and $\delta_{\text{Mars}}$. Moreover, the analytical characterization of the asteroidal effects on these three geocentric parameters is still much more complicated than for the heliocentric orbital parameters, because the determination of the vector from the EMB to planet, even assuming that the two bodies follow a Keplerian trajectory, is a much more complex function that involves Legendre polynomials, than for the heliocentric vector of the planet.

We recall that to facilitate the computations, we replaced the Moon and the Earth by the EMB, i.e. a single point whose mass is the sum of the masses of the two bodies, and which is located at their barycenter. We will calculate in a future paper the differential effects of our asteroids on the Earth to Moon distance,
<table><caption>**Table 7.** Uncertainty on the values of the six Mars orbital elements, for a 1000 y, because of the lack of accuracy in the measurements of the asteroid mass.</caption><thead><tr><th></th><th>Asteroid</th><th>Δa<br>[m]</th><th>Δe × a<sub>Mars</sub></th><th>Δi<br>×10<sup>-6</sup>"</th><th>ΔΩ<br>×10<sup>-6</sup>"</th><th>Δω<br>×10<sup>-6</sup>"</th><th>Δl<br>×10<sup>-6</sup>"</th></tr></thead><tbody><tr><td>1</td><td>Ceres</td><td>0.529</td><td>1.653</td><td>8.599</td><td>521.073</td><td>129.366</td><td>251.518</td></tr><tr><td>2</td><td>Pallas</td><td>1.451</td><td>13.802</td><td>30.792</td><td>753.100</td><td>149.618</td><td>975.491</td></tr><tr><td>3</td><td>Juno</td><td>7.800</td><td>33.243</td><td>76.361</td><td>1656.535</td><td>328.253</td><td>4966.704</td></tr><tr><td>4</td><td>Vesta</td><td>1.220</td><td>10.415</td><td>12.646</td><td>220.471</td><td>112.846</td><td>1391.602</td></tr><tr><td>6</td><td>Hebe</td><td>26.943</td><td>113.656</td><td>113.409</td><td>397.056</td><td>1125.430</td><td>7581.599</td></tr><tr><td>7</td><td>Iris</td><td>13.510</td><td>117.882</td><td>15.492</td><td>2384.353</td><td>341.588</td><td>649.930</td></tr><tr><td>8</td><td>Flora</td><td>3.023</td><td>37.158</td><td>27.890</td><td>161.404</td><td>155.523</td><td>2446.625</td></tr><tr><td>9</td><td>Metis</td><td>15.917</td><td>141.691</td><td>22.455</td><td>2338.922</td><td>1397.893</td><td>1503.667</td></tr><tr><td>10</td><td>Hygiea</td><td>1.771</td><td>18.102</td><td>7.451</td><td>438.131</td><td>189.439</td><td>228.589</td></tr><tr><td>11</td><td>Parthenope</td><td>0.156</td><td>1.243</td><td>1.457</td><td>3.742</td><td>13.077</td><td>43.074</td></tr><tr><td>13</td><td>Egeria</td><td>2.825</td><td>11.215</td><td>15.148</td><td>3103.934</td><td>325.153</td><td>1291.582</td></tr><tr><td>14</td><td>Irene</td><td>2.587</td><td>14.058</td><td>21.611</td><td>952.864</td><td>333.820</td><td>2087.306</td></tr><tr><td>15</td><td>Eunomia</td><td>1.998</td><td>19.835</td><td>28.267</td><td>778.386</td><td>64.206</td><td>2325.111</td></tr><tr><td>16</td><td>Psyche</td><td>1.392</td><td>17.813</td><td>13.334</td><td>291.731</td><td>108.607</td><td>837.838</td></tr><tr><td>17</td><td>Thetis</td><td>0.284</td><td>2.452</td><td>2.104</td><td>4.027</td><td>18.017</td><td>201.725</td></tr><tr><td>18</td><td>Melpomene</td><td>2.306</td><td>19.408</td><td>46.828</td><td>466.697</td><td>66.071</td><td>167.909</td></tr><tr><td>19</td><td>Fortuna</td><td>6.824</td><td>43.364</td><td>4.299</td><td>533.764</td><td>257.431</td><td>6505.832</td></tr><tr><td>20</td><td>Massalia</td><td>13.635</td><td>40.338</td><td>1.532</td><td>342.108</td><td>503.369</td><td>8092.128</td></tr><tr><td>21</td><td>Lutetia</td><td>0.623</td><td>4.554</td><td>1.011</td><td>19.306</td><td>48.065</td><td>700.052</td></tr><tr><td>22</td><td>Kalliope</td><td>0.066</td><td>0.885</td><td>0.525</td><td>91.106</td><td>9.871</td><td>47.539</td></tr><tr><td>24</td><td>Themis</td><td>5.404</td><td>19.255</td><td>0.695</td><td>200.095</td><td>469.326</td><td>6137.619</td></tr><tr><td>28</td><td>Bellona</td><td>4.612</td><td>21.258</td><td>61.162</td><td>326.950</td><td>805.754</td><td>2262.851</td></tr><tr><td>29</td><td>Amphitrite</td><td>6.044</td><td>58.018</td><td>78.665</td><td>705.943</td><td>989.440</td><td>4437.654</td></tr><tr><td>31</td><td>Euphrosyne</td><td>2.516</td><td>10.597</td><td>10.184</td><td>686.061</td><td>90.927</td><td>1162.892</td></tr><tr><td>45</td><td>Eugenia</td><td>2.698</td><td>35.770</td><td>67.508</td><td>550.718</td><td>622.549</td><td>835.918</td></tr><tr><td>46</td><td>Hestia</td><td>12.799</td><td>143.081</td><td>70.467</td><td>3706.341</td><td>696.209</td><td>7083.534</td></tr><tr><td>47</td><td>Aglaja</td><td>0.293</td><td>1.693</td><td>4.248</td><td>33.190</td><td>11.690</td><td>23.301</td></tr><tr><td>48</td><td>Doris</td><td>1.970</td><td>14.104</td><td>28.663</td><td>1020.682</td><td>394.806</td><td>476.894</td></tr><tr><td>49</td><td>Pales</td><td>0.303</td><td>5.085</td><td>1.031</td><td>81.624</td><td>22.385</td><td>38.901</td></tr><tr><td>52</td><td>Europa</td><td>2.225</td><td>15.466</td><td>29.815</td><td>137.466</td><td>413.986</td><td>3011.237</td></tr><tr><td>65</td><td>Cybele</td><td>0.734</td><td>8.087</td><td>8.896</td><td>124.801</td><td>72.713</td><td>696.421</td></tr><tr><td>87</td><td>Sylvia</td><td>0.008</td><td>0.031</td><td>0.068</td><td>9.277</td><td>0.928</td><td>1.695</td></tr><tr><td>88</td><td>Thisbe</td><td>2.240</td><td>14.125</td><td>15.965</td><td>975.793</td><td>74.820</td><td>1927.173</td></tr><tr><td>90</td><td>Antiope</td><td>0.010</td><td>0.049</td><td>0.005</td><td>0.052</td><td>0.422</td><td>7.596</td></tr><tr><td>107</td><td>Camilla</td><td>0.082</td><td>0.160</td><td>1.552</td><td>28.741</td><td>13.405</td><td>7.417</td></tr><tr><td>111</td><td>Ate</td><td>49.523</td><td>322.444</td><td>380.867</td><td>12116.761</td><td>7560.482</td><td>74212.364</td></tr><tr><td>121</td><td>Hermione</td><td>0.055</td><td>0.619</td><td>0.245</td><td>32.379</td><td>5.239</td><td>72.732</td></tr><tr><td>130</td><td>Elektra</td><td>0.285</td><td>1.439</td><td>6.656</td><td>23.729</td><td>35.529</td><td>194.223</td></tr><tr><td>165</td><td>Loreley</td><td>7.623</td><td>32.525</td><td>160.098</td><td>3710.437</td><td>586.050</td><td>563.742</td></tr><tr><td>189</td><td>Phthia</td><td>0.045</td><td>0.135</td><td>0.087</td><td>7.206</td><td>2.416</td><td>52.144</td></tr><tr><td>243</td><td>Ida</td><td>0.005</td><td>0.005</td><td>0.002</td><td>0.167</td><td>0.151</td><td>1.336</td></tr><tr><td>253</td><td>Mathilde</td><td>0.008</td><td>0.017</td><td>0.056</td><td>1.788</td><td>0.546</td><td>9.536</td></tr><tr><td>283</td><td>Emma</td><td>0.011</td><td>0.106</td><td>0.201</td><td>4.982</td><td>0.488</td><td>1.214</td></tr></tbody></table>

which likely has some interesting implications on lunar laser ranging (LLR) measurements. We calculated the geocentric coordinates of Mars by simply subtracting the heliocentric vector of the EMB from that of Mars. The effects of the leading asteroids Ceres, Pallas and Vesta on the EMB-Mars distance are plotted in Figs. 12–14. In any case the signal is dominated by a single wave with increasing amplitude and a period that slightly varies at each oscillation, which makes the adjustment through a least-squares analysis uneasy, indeed even inadequate. As expected, the period is clearly close to the synodic period of the planet, which changes slightly at each synodic cycle. The increase of amplitudes with respect to the time are undoubtedly caused by the secular effects of the asteroid on the EMB and Mars orbital parameters, characterized by quadratic trends as studied in Sect. 4.3. These secular effects gradually enhance the effects on the EMB-Mars vector. The peak-to-peak amplitudes of the effects of Ceres, Pallas and Vesta on the EMB-to-Mars distance are, 34 km, 31 km and 116 km respectively during the 100 y time interval considered.

In Figs. 15 and 16, we present the influence of Ceres on Mars geocentric right ascension and declination. Their peak-to-peak amplitudes reach 0.08 and 0.05 arcseconds during the 100 y time interval. The peaks correspond to Mars oppositions, for which the relatively close distance of the planet from the Earth enhances the angular effects of the perturbations. In parallel we can observe the asymmetrical nature of the two curves originating from the fact that the effects are proportional to the distance from the Earth to Mars, which varies roughly from 0.5 to 2.5 AU.

## 7. ASETEP database

### 7.1. Construction of the ASETEP database.

As a result of all our computations, we constructed the ASETEP (asteroid effects on the terrestrial planets) database, which contains at its present extent the curves showing the individual effects of each of the 43 asteroids listed in Table 1 on each of the 6 orbital elements for any of the 4 terrestrial planets (Mercury,
<table><caption>**Table 8.** Coefficients of Fourier and Poisson series for the effects of Ceres on $a_{\text{EMB}}$.</caption><thead><tr><th>Period (Day)</th><th>Period (Year)</th><th>SIN ×10<sup>-2</sup></th><th>COS ×10<sup>-2</sup></th><th>×10<sup>-1</sup></th><th>T SIN ×10<sup>-5</sup></th><th>T COS ×10<sup>-5</sup></th><th>×10<sup>-5</sup></th><th>Amplitude ×10<sup>0</sup></th><th>Asteroid</th></tr></thead><tbody><tr><td>4264.15</td><td>11.67</td><td>-5.833</td><td>-4.725</td><td>0.751</td><td>4.218</td><td>4.478</td><td>6.332</td><td>0.089</td><td>1_Ceres</td></tr><tr><td>2785.82</td><td>7.63</td><td>14.760</td><td>-17.093</td><td>2.258</td><td>-9.642</td><td>14.800</td><td>20.930</td><td>0.164</td><td>1_Ceres</td></tr><tr><td>1048.02</td><td>2.87</td><td>14.503</td><td>-23.187</td><td>2.735</td><td>-13.280</td><td>16.173</td><td>22.872</td><td>0.199</td><td>1_Ceres</td></tr><tr><td>645.80</td><td>1.77</td><td>23.772</td><td>1.851</td><td>2.384</td><td>4.417</td><td>7.592</td><td>10.736</td><td>0.389</td><td>1_Ceres</td></tr><tr><td>583.74</td><td>1.60</td><td>-4.586</td><td>80.143</td><td>8.027</td><td>2.322</td><td>-39.840</td><td>56.342</td><td>0.048</td><td>1_Ceres</td></tr><tr><td>466.62</td><td>1.28</td><td>-107.748</td><td>26.526</td><td>11.096</td><td>-3.096</td><td>-9.959</td><td>14.085</td><td>1.127</td><td>1_Ceres</td></tr><tr><td>399.51</td><td>1.09</td><td>-1.046</td><td>24.761</td><td>2.478</td><td>-0.533</td><td>-9.255</td><td>13.089</td><td>0.106</td><td>1_Ceres</td></tr><tr><td>365.29</td><td>1.00</td><td>2.167</td><td>4.551</td><td>0.504</td><td>1.556</td><td>-0.782</td><td>1.106</td><td>0.131</td><td>1_Ceres</td></tr><tr><td>337.45</td><td>0.92</td><td>0.631</td><td>7.730</td><td>0.776</td><td>-0.767</td><td>-3.066</td><td>4.336</td><td>0.080</td><td>1_Ceres</td></tr><tr><td>322.86</td><td>0.88</td><td>8.798</td><td>48.792</td><td>4.958</td><td>1.456</td><td>-5.255</td><td>7.431</td><td>0.358</td><td>1_Ceres</td></tr><tr><td>312.83</td><td>0.86</td><td>7.583</td><td>-2.635</td><td>0.803</td><td>-4.257</td><td>0.966</td><td>1.366</td><td>0.216</td><td>1_Ceres</td></tr><tr><td>299.95</td><td>0.82</td><td>-6.889</td><td>11.029</td><td>1.300</td><td>5.119</td><td>-5.229</td><td>7.395</td><td>0.294</td><td>1_Ceres</td></tr><tr><td>291.95</td><td>0.80</td><td>24.131</td><td>76.178</td><td>7.991</td><td>-12.129</td><td>-38.203</td><td>54.027</td><td>0.144</td><td>1_Ceres</td></tr><tr><td>281.05</td><td>0.77</td><td>9.514</td><td>-6.059</td><td>1.128</td><td>-2.537</td><td>2.629</td><td>3.718</td><td>0.399</td><td>1_Ceres</td></tr><tr><td>273.99</td><td>0.75</td><td>-51.020</td><td>-30.936</td><td>5.967</td><td>25.466</td><td>14.530</td><td>20.548</td><td>0.519</td><td>1_Ceres</td></tr></tbody><tfoot><tr><td colspan="10">AMPLITUDE Before= 22.943527 After= 0.353778</td></tr></tfoot></table>

**Notes.** Bias: $0.561802 \times 10^{+3}$; linear: $-0.549138$; $T^2$: $0.133527 \times 10^{-3}$. The next-to-last column indicates the peak-to-peak amplitude of the total mixed component for the 100 y time interval.

<figure><img src="image_11.png" /><figcaption>**Fig. 10.** Effects of Pallas on $a_{\text{EMB}}$, for a 100 y time interval, from JD 2009.81. Residuals (in green) after removing the sinusoidal components from the least-squares analysis.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Effects of Vesta on $a_{\text{EMB}}$, for a 100 y time interval, from JD 2009.81. Residuals (in green) after removing the sinusoidal components from the least-squares analysis.</figcaption></figure>

Venus, the EMB, Mars) both for a 100 y time interval and a 1000 y time interval from JD 2009.7. In parallel we deliver for each of the curves the corresponding tables giving the Fourier and Poisson components for fitting these curves. This ASETEP database will be subject to regular and continuous updates with an increasing number of new asteroids taken into account, if their mass is determined with a satisfactory accuracy, i.e. better than 65% uncertainty, as for the 43 asteroids listed in Table 1. One of the advantages of the ASETEP database is that the figures and tables already produced for a given asteroid can be updated immediately each time a new and more accurate determination of the mass of this asteroid is reached, by simple proportionnality.

## 7.2. Individual and combined asteroidal effects

In all computations we considered the effect of each asteroid taken individually, which neglects the effects from the combined gravitational attraction of the asteroids themselves. A priori these effects which result principally from the reciprocal gravitational interactions of the asteroids affecting their positions, are of second order because of the smallness of asteroid masses. They can be neglected at first approximation as long as

the interval of time of the study $\Delta t$ is short, as for $\Delta t = 100$ y. To obtain some information about the error caused by neglecting the effects due to the combination of asteroids, we carried out a test with numerical integration: it consists of studying the difference between two numerical computations: first we summed up the individual effects of the 43 asteroids considered here on the Mars orbital parameters, each of which was obtained in the frame of an n-body problem with $n = 10$ (Sun + 8 planets + asteroid); second we included all asteroids in the frame of an n-body problem with $n = 52$ (Sun + 8 planets + 43 asteroids).

In Fig. 17, we plot the combined effect of the 43 asteroids taken together on the Mars semi-major axis. The residuals of the two computation methods presented above do not exceed 2.5 m peak-to-peak, as can be seen in Fig. 18. This is by two orders less than the signal itself (650 m peak-to-peak), and the error caused by treating individual effects separately can be neglected in first approximation. This is still true for the other Mars parameters (see Appendix B), for which the errors are below $1.2 \times 10^{-11}$ for e (Figs. B.1 and B.2), 1.8 mas for i (Figs. B.3 and B.4), $10^{-2}$ mas for $\Omega$ (Figs. B.5 and B.6), $5 \times 10^{-3}$ mas for $\tilde{\omega}$ (Figs. B.7 and B.8).
<table><caption>**Table 9.** Peak-to-peak amplitudes of the variations of the six EMB orbital elements due to each given asteroid for a 1000 y time interval starting from JD 2009.81.</caption><thead><tr><th>Asteroid</th><th>Δa<br>m</th><th>Δe × a<sub>Earth</sub><br>m</th><th>Δi<br>×10<sup>-6</sup>"</th><th>ΔΩ<br>×10<sup>-6</sup>"</th><th>Δω<br>×10<sup>-6</sup>"</th><th>Δλ<br>×10<sup>-6</sup>"</th></tr></thead><tbody><tr><td>1 Ceres</td><td>25.443</td><td>543.236</td><td>4929.530</td><td>464926.356</td><td>6311.044</td><td>2882.002</td></tr><tr><td>2 Pallas</td><td>10.608</td><td>147.959</td><td>219.655</td><td>1138450.592</td><td>2759.495</td><td>26618.466</td></tr><tr><td>3 Juno</td><td>1.347</td><td>29.962</td><td>29.840</td><td>55190.496</td><td>1976.844</td><td>5736.100</td></tr><tr><td>4 Vesta</td><td>19.679</td><td>221.168</td><td>1729.805</td><td>303243.082</td><td>55222.264</td><td>40243.056</td></tr><tr><td>6 Hebe</td><td>1.575</td><td>44.829</td><td>128.567</td><td>66754.335</td><td>1169.070</td><td>2754.659</td></tr><tr><td>7 Iris</td><td>1.916</td><td>49.911</td><td>67.099</td><td>5392.495</td><td>2470.755</td><td>2850.537</td></tr><tr><td>8 Flora</td><td>1.141</td><td>36.128</td><td>61.188</td><td>13731.489</td><td>989.258</td><td>194.963</td></tr><tr><td>9 Metis</td><td>1.324</td><td>23.418</td><td>97.963</td><td>18171.585</td><td>2828.013</td><td>301.391</td></tr><tr><td>10 Hygiea</td><td>2.154</td><td>28.291</td><td>122.134</td><td>18072.496</td><td>6670.653</td><td>4666.845</td></tr><tr><td>11 Parthenope</td><td>0.402</td><td>8.018</td><td>19.419</td><td>6613.715</td><td>1071.967</td><td>1323.398</td></tr><tr><td>13 Egeria</td><td>0.647</td><td>6.346</td><td>114.239</td><td>52651.902</td><td>504.693</td><td>173.308</td></tr><tr><td>14 Irene</td><td>0.624</td><td>17.746</td><td>52.094</td><td>3431.527</td><td>414.283</td><td>1124.330</td></tr><tr><td>15 Eunomia</td><td>2.074</td><td>60.980</td><td>236.253</td><td>49917.699</td><td>1191.803</td><td>1711.442</td></tr><tr><td>16 Psyche</td><td>1.673</td><td>37.652</td><td>18.645</td><td>13866.105</td><td>2789.469</td><td>7791.016</td></tr><tr><td>17 Thetis</td><td>0.119</td><td>0.879</td><td>4.588</td><td>1577.995</td><td>306.011</td><td>204.078</td></tr><tr><td>18 Melpomene</td><td>0.631</td><td>16.064</td><td>17.729</td><td>13773.197</td><td>234.733</td><td>1303.507</td></tr><tr><td>19 Fortuna</td><td>1.107</td><td>30.447</td><td>8.395</td><td>5701.767</td><td>689.533</td><td>2320.502</td></tr><tr><td>20 Massalia</td><td>0.688</td><td>1.965</td><td>1.549</td><td>1315.434</td><td>922.438</td><td>953.087</td></tr><tr><td>21 Lutetia</td><td>0.235</td><td>6.677</td><td>7.069</td><td>685.747</td><td>598.547</td><td>132.662</td></tr><tr><td>22 Kalliope</td><td>0.381</td><td>3.543</td><td>41.715</td><td>8466.995</td><td>239.465</td><td>775.256</td></tr><tr><td>24 Themis</td><td>0.456</td><td>7.424</td><td>1.398</td><td>908.453</td><td>142.142</td><td>794.665</td></tr><tr><td>28 Bellona</td><td>0.675</td><td>10.883</td><td>40.815</td><td>25052.303</td><td>1082.842</td><td>801.938</td></tr><tr><td>29 Amphitrite</td><td>0.773</td><td>14.573</td><td>6.314</td><td>36621.458</td><td>577.161</td><td>632.069</td></tr><tr><td>31 Euphrosyne</td><td>0.318</td><td>2.090</td><td>23.015</td><td>18359.137</td><td>122.700</td><td>67.111</td></tr><tr><td>45 Eugenia</td><td>0.618</td><td>10.152</td><td>35.986</td><td>25097.892</td><td>1631.268</td><td>501.838</td></tr><tr><td>46 Hestia</td><td>3.892</td><td>78.350</td><td>3.219</td><td>18091.046</td><td>4244.686</td><td>10698.904</td></tr><tr><td>47 Aglaja</td><td>0.140</td><td>1.865</td><td>0.227</td><td>2219.030</td><td>223.737</td><td>84.944</td></tr><tr><td>48 Doris</td><td>0.412</td><td>1.419</td><td>1.141</td><td>12170.155</td><td>138.494</td><td>549.257</td></tr><tr><td>49 Pales</td><td>0.156</td><td>3.498</td><td>3.626</td><td>498.753</td><td>155.474</td><td>501.291</td></tr><tr><td>52 Europa</td><td>0.569</td><td>4.259</td><td>41.151</td><td>14775.024</td><td>553.168</td><td>2008.514</td></tr><tr><td>65 Cybele</td><td>0.349</td><td>1.907</td><td>8.621</td><td>5582.948</td><td>883.675</td><td>375.725</td></tr><tr><td>87 Sylvia</td><td>0.350</td><td>3.944</td><td>33.215</td><td>6532.738</td><td>334.762</td><td>847.921</td></tr><tr><td>88 Thisbe</td><td>0.704</td><td>11.842</td><td>32.942</td><td>2892.127</td><td>1809.126</td><td>1924.267</td></tr><tr><td>90 Antiope</td><td>0.090</td><td>0.626</td><td>0.569</td><td>128.930</td><td>74.697</td><td>95.295</td></tr><tr><td>107 Camilla</td><td>0.503</td><td>1.431</td><td>6.525</td><td>11211.884</td><td>73.744</td><td>476.956</td></tr><tr><td>111 Ate</td><td>9.223</td><td>76.054</td><td>449.106</td><td>151257.961</td><td>12292.524</td><td>28037.019</td></tr><tr><td>121 Hermione</td><td>0.192</td><td>2.596</td><td>9.092</td><td>1821.472</td><td>60.102</td><td>613.565</td></tr><tr><td>130 Elektra</td><td>0.277</td><td>4.675</td><td>29.562</td><td>17655.726</td><td>232.881</td><td>840.610</td></tr><tr><td>165 Loreley</td><td>1.280</td><td>5.862</td><td>102.851</td><td>29936.733</td><td>1867.874</td><td>974.511</td></tr><tr><td>189 Phthia</td><td>0.056</td><td>0.036</td><td>0.066</td><td>62.766</td><td>3.047</td><td>179.217</td></tr><tr><td>243 Ida</td><td>0.076</td><td>0.131</td><td>0.016</td><td>8.814</td><td>11.517</td><td>233.687</td></tr><tr><td>253 Mathilde</td><td>0.069</td><td>0.386</td><td>0.016</td><td>231.658</td><td>24.807</td><td>92.696</td></tr><tr><td>283 Emma</td><td>0.094</td><td>1.347</td><td>3.668</td><td>1093.080</td><td>54.125</td><td>181.180</td></tr></tbody></table>

<figure><img src="image_13.png" /><figcaption>**Fig. 12.** Effects of Ceres on the distance from the EMB to Mars for a 100 y time interval from JD 2009.81. The curve in green corresponds to the post-fit residuals.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Fig. 13.** Effects of Pallas on the distance from the EMB to Mars for a 100 y time interval from JD 2009.81. The curve in green corresponds to the post-fit residuals.</figcaption></figure>
<table><caption>**Table 10.** Peak-to-peak amplitudes of the variations of the six Mercury orbital elements due to each given asteroid for a 1000 y time interval starting from JD 2009.81.</caption><thead><tr><th>Asteroid</th><th>Δa<br>m</th><th>Δe × a<sub>Mercury</sub><br>m</th><th>Δi<br>×10<sup>-6</sup>"</th><th>ΔΩ<br>×10<sup>-6</sup>"</th><th>Δω<br>×10<sup>-6</sup>"</th><th>Δλ<br>×10<sup>-6</sup>"</th></tr></thead><tbody><tr><td>1 Ceres</td><td>1.027</td><td>17.040</td><td>506.486</td><td>1525.642</td><td>5274.314</td><td>1031.921</td></tr><tr><td>2 Pallas</td><td>0.236</td><td>21.632</td><td>546.181</td><td>3524.193</td><td>680.519</td><td>2027.054</td></tr><tr><td>3 Juno</td><td>0.031</td><td>0.226</td><td>23.731</td><td>194.155</td><td>53.392</td><td>1045.685</td></tr><tr><td>4 Vesta</td><td>0.444</td><td>1.697</td><td>312.317</td><td>1284.460</td><td>2605.141</td><td>9801.069</td></tr><tr><td>6 Hebe</td><td>0.035</td><td>1.712</td><td>38.624</td><td>149.068</td><td>85.198</td><td>742.517</td></tr><tr><td>7 Iris</td><td>0.050</td><td>0.311</td><td>3.638</td><td>190.978</td><td>91.108</td><td>1633.283</td></tr><tr><td>8 Flora</td><td>0.036</td><td>0.610</td><td>10.833</td><td>69.548</td><td>81.990</td><td>1386.485</td></tr><tr><td>9 Metis</td><td>0.032</td><td>0.097</td><td>6.998</td><td>53.990</td><td>161.183</td><td>1166.247</td></tr><tr><td>10 Hygiea</td><td>0.054</td><td>1.107</td><td>12.770</td><td>512.880</td><td>391.011</td><td>1769.981</td></tr><tr><td>11 Parthenope</td><td>0.020</td><td>0.202</td><td>5.388</td><td>50.541</td><td>51.114</td><td>965.908</td></tr><tr><td>13 Egeria</td><td>0.033</td><td>0.866</td><td>7.069</td><td>161.930</td><td>104.715</td><td>559.951</td></tr><tr><td>14 Irene</td><td>0.043</td><td>0.472</td><td>6.181</td><td>2.022</td><td>61.844</td><td>2266.552</td></tr><tr><td>15 Eunomia</td><td>0.046</td><td>1.757</td><td>38.713</td><td>391.935</td><td>218.771</td><td>398.350</td></tr><tr><td>16 Psyche</td><td>0.037</td><td>0.389</td><td>9.404</td><td>146.207</td><td>111.991</td><td>868.428</td></tr><tr><td>17 Thetis</td><td>0.021</td><td>0.015</td><td>1.276</td><td>9.532</td><td>10.570</td><td>1164.296</td></tr><tr><td>18 Melpomene</td><td>0.024</td><td>0.282</td><td>6.566</td><td>45.892</td><td>22.050</td><td>776.751</td></tr><tr><td>19 Fortuna</td><td>0.028</td><td>0.299</td><td>2.340</td><td>119.341</td><td>79.661</td><td>99.164</td></tr><tr><td>20 Massalia</td><td>0.021</td><td>0.275</td><td>0.760</td><td>48.928</td><td>36.882</td><td>163.080</td></tr><tr><td>21 Lutetia</td><td>0.023</td><td>0.194</td><td>0.909</td><td>16.061</td><td>22.953</td><td>380.720</td></tr><tr><td>22 Kalliope</td><td>0.020</td><td>0.055</td><td>2.041</td><td>34.978</td><td>37.148</td><td>904.473</td></tr><tr><td>24 Themis</td><td>0.035</td><td>0.325</td><td>0.393</td><td>43.222</td><td>44.341</td><td>2114.352</td></tr><tr><td>28 Bellona</td><td>0.031</td><td>0.613</td><td>15.834</td><td>93.611</td><td>61.361</td><td>745.187</td></tr><tr><td>29 Amphitrite</td><td>0.032</td><td>0.383</td><td>13.422</td><td>69.474</td><td>135.370</td><td>275.212</td></tr><tr><td>31 Euphrosyne</td><td>0.015</td><td>0.351</td><td>5.117</td><td>60.569</td><td>13.988</td><td>107.380</td></tr><tr><td>45 Eugenia</td><td>0.032</td><td>0.463</td><td>16.170</td><td>134.059</td><td>102.246</td><td>1300.468</td></tr><tr><td>46 Hestia</td><td>0.056</td><td>0.989</td><td>8.991</td><td>226.888</td><td>156.690</td><td>2138.406</td></tr><tr><td>47 Aglaja</td><td>0.019</td><td>0.084</td><td>0.771</td><td>6.255</td><td>11.698</td><td>333.500</td></tr><tr><td>48 Doris</td><td>0.024</td><td>0.331</td><td>5.801</td><td>81.630</td><td>39.975</td><td>1628.596</td></tr><tr><td>49 Pales</td><td>0.026</td><td>0.041</td><td>0.345</td><td>16.042</td><td>9.940</td><td>566.625</td></tr><tr><td>52 Europa</td><td>0.032</td><td>0.298</td><td>11.911</td><td>65.363</td><td>66.695</td><td>1956.630</td></tr><tr><td>65 Cybele</td><td>0.020</td><td>0.069</td><td>4.196</td><td>52.779</td><td>46.187</td><td>221.947</td></tr><tr><td>87 Sylvia</td><td>0.026</td><td>0.165</td><td>1.920</td><td>18.840</td><td>41.512</td><td>869.579</td></tr><tr><td>88 Thisbe</td><td>0.029</td><td>0.288</td><td>3.084</td><td>107.780</td><td>73.901</td><td>517.634</td></tr><tr><td>90 Antiope</td><td>0.029</td><td>0.026</td><td>0.075</td><td>2.631</td><td>3.275</td><td>544.525</td></tr><tr><td>107 Camilla</td><td>0.018</td><td>0.242</td><td>6.226</td><td>53.845</td><td>24.304</td><td>708.428</td></tr><tr><td>111 Ate</td><td>0.195</td><td>4.328</td><td>93.162</td><td>1380.844</td><td>1113.267</td><td>1479.395</td></tr><tr><td>121 Hermione</td><td>0.015</td><td>0.088</td><td>0.568</td><td>0.221</td><td>15.036</td><td>342.469</td></tr><tr><td>130 Elektra</td><td>0.017</td><td>0.498</td><td>11.278</td><td>34.221</td><td>10.621</td><td>593.566</td></tr><tr><td>165 Loreley</td><td>0.042</td><td>0.823</td><td>22.573</td><td>186.959</td><td>138.388</td><td>1484.663</td></tr><tr><td>189 Phthia</td><td>0.019</td><td>0.017</td><td>0.024</td><td>0.565</td><td>0.252</td><td>792.788</td></tr><tr><td>243 Ida</td><td>0.023</td><td>0.012</td><td>0.005</td><td>0.243</td><td>0.353</td><td>1287.098</td></tr><tr><td>253 Mathilde</td><td>0.011</td><td>0.011</td><td>0.099</td><td>1.381</td><td>0.617</td><td>456.920</td></tr><tr><td>283 Emma</td><td>0.012</td><td>0.043</td><td>0.732</td><td>8.355</td><td>5.916</td><td>591.482</td></tr></tbody></table>

<figure><img src="image_15.png" /><figcaption>**Fig. 14.** Effects of Vesta on the distance from the EMB to Mars for a 100 y time interval from JD 2009.81. The curve in green corresponds to the post-fit residuals.</figcaption></figure>

### 7.3. Use of ASETEP for ephemerides study

Most of these asteroids have a very large mass uncertainty, sometimes by far exceeding their nominal values, which is why it appears to be ill-advised to include them individually in the ephemerides computations. But they can be assigned a reasonable distribution of masses according to some models such as the statistical asteroid model (Tedesco et al. 2005). The problem of mass undetermination of small asteroids was there partly solved by replacing all asteroids except for the 303 most perturbing ones by a ring. Through Monte Carlo experiments, where asteroids are assigned random but reasonable masses, it is possible to calculate the corresponding perturbations from the selected set of asteroids (Fienga et al. 2009) An estimate of the mass of the ring by Kuchynka et al. (2008) was $M_{\text{ring}} = 1 \pm 0.3 \times 10^{-10} M_{\text{Sun}}$ for a ring at 3.14 AU. Yet compared this statistical evaluation, the uncertainty concerning the largest asteroids looks like a severe limitation for the ephemeris accuracy. For instance, the mass
<table><caption>**Table 11.** Peak-to-peak amplitudes of the variations of the six Venus orbital elements due to the action of each asteroid for a 1000 y time interval starting from JD 2009.81.</caption><thead><tr><th></th><th>Asteroid</th><th>Δa<br>m</th><th>Δe × a<sub>Venus</sub><br>m</th><th>Δi<br>×10<sup>-6</sup>"</th><th>ΔΩ<br>×10<sup>-6</sup>"</th><th>Δω<br>×10<sup>-6</sup>"</th><th>Δλ<br>×10<sup>-6</sup>"</th></tr></thead><tbody><tr><td>1</td><td>Ceres</td><td>8.873</td><td>104.443</td><td>48.038</td><td>31824.583</td><td>36061.152</td><td>70697.133</td></tr><tr><td>2</td><td>Pallas</td><td>3.622</td><td>42.975</td><td>1518.868</td><td>4481.157</td><td>2990.961</td><td>13040.604</td></tr><tr><td>3</td><td>Juno</td><td>0.520</td><td>12.100</td><td>67.545</td><td>237.562</td><td>844.492</td><td>749.814</td></tr><tr><td>4</td><td>Vesta</td><td>5.436</td><td>113.260</td><td>411.520</td><td>6658.590</td><td>35311.328</td><td>16328.430</td></tr><tr><td>6</td><td>Hebe</td><td>0.363</td><td>12.581</td><td>87.162</td><td>486.578</td><td>2490.318</td><td>662.926</td></tr><tr><td>7</td><td>Iris</td><td>0.647</td><td>17.653</td><td>1.688</td><td>902.456</td><td>2010.286</td><td>3453.414</td></tr><tr><td>8</td><td>Flora</td><td>0.225</td><td>10.603</td><td>18.079</td><td>143.247</td><td>567.063</td><td>1067.780</td></tr><tr><td>9</td><td>Metis</td><td>0.292</td><td>10.441</td><td>8.033</td><td>336.941</td><td>2196.530</td><td>174.656</td></tr><tr><td>10</td><td>Hygiea</td><td>0.655</td><td>14.605</td><td>16.991</td><td>2080.629</td><td>4350.355</td><td>1846.343</td></tr><tr><td>11</td><td>Parthenope</td><td>0.129</td><td>0.871</td><td>9.425</td><td>12.093</td><td>1109.277</td><td>555.197</td></tr><tr><td>13</td><td>Egeria</td><td>0.239</td><td>0.504</td><td>53.465</td><td>939.036</td><td>928.800</td><td>2454.464</td></tr><tr><td>14</td><td>Irene</td><td>0.179</td><td>4.529</td><td>2.946</td><td>292.024</td><td>649.127</td><td>210.247</td></tr><tr><td>15</td><td>Eunomia</td><td>0.435</td><td>17.318</td><td>69.851</td><td>2515.550</td><td>1056.317</td><td>2316.562</td></tr><tr><td>16</td><td>Psyche</td><td>0.681</td><td>8.146</td><td>18.972</td><td>245.119</td><td>1104.041</td><td>647.785</td></tr><tr><td>17</td><td>Thetis</td><td>0.034</td><td>0.526</td><td>2.198</td><td>4.257</td><td>214.398</td><td>83.850</td></tr><tr><td>18</td><td>Melpomene</td><td>0.191</td><td>4.177</td><td>16.546</td><td>5.985</td><td>669.946</td><td>832.947</td></tr><tr><td>19</td><td>Fortuna</td><td>0.195</td><td>8.576</td><td>6.152</td><td>369.171</td><td>512.983</td><td>1946.517</td></tr><tr><td>20</td><td>Massalia</td><td>0.176</td><td>1.709</td><td>1.578</td><td>146.463</td><td>923.534</td><td>881.744</td></tr><tr><td>21</td><td>Lutetia</td><td>0.085</td><td>0.941</td><td>0.244</td><td>7.793</td><td>692.551</td><td>295.897</td></tr><tr><td>22</td><td>Kalliope</td><td>0.113</td><td>1.642</td><td>6.392</td><td>315.607</td><td>171.884</td><td>98.980</td></tr><tr><td>24</td><td>Themis</td><td>0.135</td><td>1.459</td><td>0.592</td><td>110.140</td><td>457.939</td><td>793.104</td></tr><tr><td>28</td><td>Bellona</td><td>0.267</td><td>0.896</td><td>34.004</td><td>52.256</td><td>1521.361</td><td>1362.879</td></tr><tr><td>29</td><td>Amphitrite</td><td>0.216</td><td>5.020</td><td>43.100</td><td>337.808</td><td>252.462</td><td>1806.797</td></tr><tr><td>31</td><td>Euphrosyne</td><td>0.101</td><td>0.899</td><td>21.759</td><td>238.441</td><td>141.318</td><td>516.297</td></tr><tr><td>45</td><td>Eugenia</td><td>0.169</td><td>4.213</td><td>33.420</td><td>75.524</td><td>850.141</td><td>667.179</td></tr><tr><td>46</td><td>Hestia</td><td>0.652</td><td>13.456</td><td>21.361</td><td>588.879</td><td>3725.682</td><td>168.093</td></tr><tr><td>47</td><td>Aglaja</td><td>0.085</td><td>0.211</td><td>2.609</td><td>23.395</td><td>250.509</td><td>277.676</td></tr><tr><td>48</td><td>Doris</td><td>0.134</td><td>0.905</td><td>15.407</td><td>180.786</td><td>166.177</td><td>1191.363</td></tr><tr><td>49</td><td>Pales</td><td>0.083</td><td>1.194</td><td>0.484</td><td>65.283</td><td>55.044</td><td>315.793</td></tr><tr><td>52</td><td>Europa</td><td>0.153</td><td>0.666</td><td>21.136</td><td>112.815</td><td>824.430</td><td>911.623</td></tr><tr><td>65</td><td>Cybele</td><td>0.119</td><td>0.888</td><td>7.796</td><td>66.585</td><td>786.489</td><td>359.006</td></tr><tr><td>87</td><td>Sylvia</td><td>0.131</td><td>1.011</td><td>4.817</td><td>230.323</td><td>445.904</td><td>256.907</td></tr><tr><td>88</td><td>Thisbe</td><td>0.261</td><td>0.760</td><td>3.815</td><td>476.506</td><td>1825.573</td><td>703.813</td></tr><tr><td>90</td><td>Antiope</td><td>0.061</td><td>0.103</td><td>0.056</td><td>3.868</td><td>94.980</td><td>676.628</td></tr><tr><td>107</td><td>Camilla</td><td>0.154</td><td>0.265</td><td>14.943</td><td>65.902</td><td>140.830</td><td>170.971</td></tr><tr><td>111</td><td>Ate</td><td>2.328</td><td>15.164</td><td>198.173</td><td>6787.420</td><td>15143.916</td><td>10854.568</td></tr><tr><td>121</td><td>Hermione</td><td>0.064</td><td>0.781</td><td>1.315</td><td>48.973</td><td>144.406</td><td>429.713</td></tr><tr><td>130</td><td>Elektra</td><td>0.098</td><td>1.419</td><td>24.887</td><td>146.941</td><td>390.856</td><td>201.809</td></tr><tr><td>165</td><td>Loreley</td><td>0.361</td><td>1.349</td><td>41.093</td><td>1179.523</td><td>1587.762</td><td>431.086</td></tr><tr><td>189</td><td>Phthia</td><td>0.055</td><td>0.035</td><td>0.069</td><td>1.701</td><td>6.693</td><td>268.361</td></tr><tr><td>243</td><td>Ida</td><td>0.086</td><td>0.105</td><td>0.014</td><td>0.812</td><td>29.902</td><td>879.399</td></tr><tr><td>253</td><td>Mathilde</td><td>0.040</td><td>0.123</td><td>0.267</td><td>2.917</td><td>32.359</td><td>300.357</td></tr><tr><td>283</td><td>Emma</td><td>0.063</td><td>0.319</td><td>1.421</td><td>46.358</td><td>89.405</td><td>437.905</td></tr></tbody></table>

evaluation of asteroid 111 Ate is $0.84 \times 10^{-10} M_{\text{Sun}}$. This is very close to the evaluation of the mass of the ring above, and is accompanied by a 22% mass uncertainty (see Table 1). This is will cause a severe degradation of the accuracy in the Mars positioning in particular. To clearly demonstrate this degradation, we show in Fig. 19 the difference of the EMB-to-Mars distance from two types of planetary ephemeris DE405 (Standish 1998) and INPOP08 (Fienga et al. 2008). In the same figure we represent also the difference in the EMB-to-Mars distance due to the uncertainty on the sole mass of 111 Ate.

We point out the astonishing correspondence between the two curves in phase, frequency and amplitude. This suggests that the differences in the values of masses for this asteroid between the two ephemerides are highly responsible for their mutual discrepancy. For comparison we show in Figs. 20 to 22 the effects on the EMB-Mars distance of the mass uncertainty of 46 Hestia, 19 Fortuna, 1 Ceres, 2 Pallas, and 4 Vesta (Appendix C). We note that for these three last asteroids, which are the biggest ones, the

effect is negligible, because of the very good mass accuracy. In contrast, for the first two the effect becomes clearly significant after a few years with an amplitude of roughly 1/8th and 1/16th that of 111 Ate (see Fig. 19) for a 100 y time interval, that is to say, 2 km and 1 km instead of 16 km.

In this context, our ASETEP database is an interesting and useful tool for understanding in a direct manner the limitations of modern ephemerides of terrestrial planets caused by the uncertainty on asteroid masses. It enables one to clearly indicate which asteroid would have to be observed preferentially in order to limitate its degradation on ephemeris accuracy. We limited our study to a small sample of 43 asteroids for which the mass is known with relatively good accuracy, i.e. with less than 65% uncertainty, as we explained in Sect. 2. But we can easily extend this study to a much larger set of asteroids. To some extent our work is similar to that of Fienga et al. (2009), who constructed the INPOP08 planetary ephemeris: they also analyzed the individual effects of very many asteroids one by one,
<figure><img src="image_16.png" /><figcaption>**Fig. 15.** Effects of Ceres on $\alpha_{\text{Mars}}$ for a 100 y time interval from JD 2009.81. The curve in green corresponds to the post-fit residuals.</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>**Fig. 16.** Effects of Ceres on $\delta_{\text{Mars}}$ for a 100 y time interval from JD 2009.81. The curve in green corresponds to a post-fit residuals.</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>**Fig. 17.** Peak-to-peak amplitude of the combined effects of the 43 asteroids on $\alpha_{\text{Mars}}$. The curve in green represents the difference between two ways of computation: adding individual perturbations, or making a combined integration.</figcaption></figure>

but without carrying out a Fourier and Poisson series decomposition of each signal as we did, and without comparing this decomposition with specific studies (Williams 1984). Indeed, the great advantage of such a decomposition is to understand, at a given period (or frequency) range, which asteroid is dominant both in terms of amplitude (related to the mass *M*) and amplitude uncertainty (related to $\delta M/M$).

<figure><img src="image_19.png" /><figcaption>**Fig. 18.** Zoom of the green curve of Fig. 17, showing the difference between two ways of computation.</figcaption></figure>

<figure><img src="image_20.png" /><figcaption>**Fig. 19.** Uncertainty of the EMB-to-Mars distance due to the lack of accuracy in Ate's mass (red curve) and difference in the EMB-to-Mars distance between DE405 and INPOP08 (black curve).</figcaption></figure>

<figure><img src="image_21.png" /><figcaption>**Fig. 20.** Uncertainty of the EMB-to-Mars distance due to the lack of accuracy in Hestia's mass (red curve) and difference in the EMB-to-Mars distance between DE405 and INPOP08 (black curve).</figcaption></figure>

Another use of ASETEP is the straightforward listing of asteroidal effects. In Table 12, we show the effects of the asteroid mass uncertainty on fundamental parameters such as the determination of the distance from the EMB to Mars, Venus and Mercury. For each planet, the value of the uncertainty on the EMP-to-planet distance is given in decreasing order of amplitude. The most problematic asteroids for Mars are 111 Ate, with a dramatic 139.0 km uncertainty, followed by 46 Hestia
<figure><img src="image_22.png" /><figcaption>**Fig. 21.** Uncertainty of the EMB-to-Mars distance due to the lack of accuracy in Fortuna's mass (red curve) and difference in the EMB-to-Mars distance between DE405 and INPOP08 (black curve).</figcaption></figure>

<figure><img src="image_23.png" /><figcaption>**Fig. 22.** Uncertainty of the EMB-to-Mars distance due to the lack of accuracy in Ceres' mass (red curve) and difference in the EMB-to-Mars distance between DE405 and INPOP08 (black curve).</figcaption></figure>

(22.0 km), 20 Massilia (14.0 km), 19 Fortuna (12.1 km), and 3 Juno (10.5 km). In the case of Venus, 46 Hestia is dominant with a 7.2 km amplitude uncertainty, followed by 111 Ate (4.2 km), 3 Juno (1.4 km), 7 Iris (1.3 km) and 16 Psyche (1.2 km). The leading uncertainties for Mercury come from 111 Ate (4.1 km), 46 Hestia (3.8 km) and 165 Loreley (0.9 km). Additional information of ASETEP is the ratio of the amplitude of the uncertainty of Mars to that of Venus and Mercury, respectively. Because of the proximity of Mars to the asteroids this ratio would a priori be significantly higher by one or two orders. This is indeed what we can notice except for a few cases. The minimum and maximum values of the ratio are 0.15 (47 Aglaja) and 186.27 (19 Fortuna) respectively for Venus; they are 0.32 (47 Aglaja) and 224.54 (31 Euphrosyne) for Mercury.

In Fig. 23 we represent as a function of the mass of each of the 43 asteroids both the amplitude of its effect on the EMB-to-Mars distance and the amplitude of the uncertainty of determination of this effect. Clearly, if we exclude 111 Ate which offers by far the largest uncertainty, the following asteroids in decreasing amplitude of effects are not the largest ones but constitute a cluster with mass ranging between $3 \times 10^{-12} M_S$ and $3 \times 10^{-11} M_S$.

## 8. Conclusion and prospects

We have carried out detailed studies of the individual effects on the orbital parameters and geocentric positions of the terrestrial planets of a set of 43 of the most interesting asteroids,

<figure><img src="image_24.png" /><figcaption>**Fig. 23.** Peak-to-peak effects (black circles) and respective uncertainty (red circles) of the EMB-to-Mars distance for each asteroid mass.</figcaption></figure>

selected for their large mass and the relatively good accuracy of their mass determination, This continues previous research of Williams (1984) and Mouret et al. (2009) in a more systematic and extended way, These authors had restricted their results sole on the planet Mars and only two orbital parameters, i.e. $a_{Mars}$ and $\lambda_{Mars}$.

1. In contrast to these authors we performed a purely numerical study, based on the numerical integration of the n-body problem. To validate our method, we compared our results for the effects of Ceres, Pallas, Vesta and other leading asteroids on $a_{Mars}$ and $\lambda_{Mars}$ with Mouret et al. (2009). We characterized the signals by a set of Fourier and Poisson series, whose arguments and amplitudes fit the analytical tables of Mouret et al. (2009) very well, which were constructed as a Fourier series as a function of the mean longitudes of Mars and the asteroids. This both validates our methodology and confirms the analytical tables of Williams (1984) and Mouret et al. (2009).

2. We extended our investigations to the other parameters $e_{Mars}$, $i_{Mars}$, $\Omega_{Mars}$, $\tilde{\omega}_{Mars}$, and we summarized in a table the amplitude of the effect of each asteroid, and the uncertainty on the modelisation of this effect directly proportional to the relative mass uncertainty $\Delta M_{ast.} / M_{ast.}$.

3. Moreover, we studied the other terrestrial planets Mercury, Venus and the Earth (the Earth-Moon barycenter) in the same way, gathering our results in recapitulative tables.

4. We also evaluated the influences of each asteroid on the equatorial coordinates $\alpha$, $\delta$ and the geocentric distance $d$ of the three terrestrial planets Mercury, Venus and Mars, as well as the uncertainty on these parameters that are directly linked to the mass uncertainty of each asteroid. Our results will be compared, with those obtained recently by Kuchynka et al. (2010), in the case of common asteroids.

5. All results of our study are gathered in the ASETEP database (asteroid effects on the terrestrial planets), i.e. the tables of the decomposition of the signals with the leading components of the fit as a function of Fourier and Poisson series, and also the curves showing the asteroid effects and the residuals after subtraction of the fit. These tables are available only for a 1000 y time span. The ASTEP database can be accessed at http://hpiers.obspm.fr/icrs-pc/. As a first step of construction of ASETEP we restricted our study on
<table><caption>**Table 12.** Uncertainty of the distances Earth-Mercury, Earth-Venus, and Earth-Mars, due to the asteroids mass unaccuracy for a 1000 y time interval.</caption><thead><tr><th colspan="2">Mars</th><th colspan="3">Venus</th><th colspan="3">Mercury</th></tr><tr><th>Asteroide</th><th>EMB-Mars</th><th>Asteroide</th><th>EMB-Venus</th><th>Ratio</th><th>Asteroide</th><th>EMB-Mercury</th><th>Ratio</th></tr></thead><tbody><tr><td>111_Ate</td><td>139044.10</td><td>46_Hestia</td><td>7184.52</td><td>3.07</td><td>111_Ate</td><td>4073.74</td><td>34.13</td></tr><tr><td>46_Hestia</td><td>22055.13</td><td>111_Ate</td><td>4188.84</td><td>33.19</td><td>46_Hestia</td><td>3782.45</td><td>5.83</td></tr><tr><td>20_Massalia</td><td>14016.59</td><td>3_Juno</td><td>1399.89</td><td>7.55</td><td>165_Loreley</td><td>916.20</td><td>0.86</td></tr><tr><td>19_Fortuna</td><td>12090.17</td><td>7_Iris</td><td>1334.75</td><td>0.85</td><td>3_Juno</td><td>824.34</td><td>12.82</td></tr><tr><td>3_Juno</td><td>10566.12</td><td>16_Psyche</td><td>1181.23</td><td>2.55</td><td>16_Psyche</td><td>718.58</td><td>4.19</td></tr><tr><td>24_Themis</td><td>10531.30</td><td>165_Loreley</td><td>931.12</td><td>0.84</td><td>48_Doris</td><td>429.47</td><td>1.30</td></tr><tr><td>6_Hebe</td><td>9948.96</td><td>88_Thisbe</td><td>806.58</td><td>5.23</td><td>88_Thisbe</td><td>428.74</td><td>9.85</td></tr><tr><td>29_Amphitrite</td><td>7805.64</td><td>18_Melpomene</td><td>770.59</td><td>1.14</td><td>24_Themis</td><td>409.79</td><td>25.70</td></tr><tr><td>52_Europa</td><td>5856.20</td><td>28_Bellona</td><td>662.53</td><td>6.81</td><td>18_Melpomene</td><td>391.27</td><td>2.24</td></tr><tr><td>28_Bellona</td><td>4510.37</td><td>13_Egeria</td><td>468.42</td><td>5.07</td><td>6_Hebe</td><td>283.14</td><td>35.14</td></tr><tr><td>8_Flora</td><td>4421.21</td><td>29_Amphitrite</td><td>436.53</td><td>17.88</td><td>7_Iris</td><td>249.36</td><td>4.53</td></tr><tr><td>88_Thisbe</td><td>4221.34</td><td>6_Hebe</td><td>432.03</td><td>23.03</td><td>19_Fortuna</td><td>220.82</td><td>54.75</td></tr><tr><td>15_Eunomia</td><td>4199.23</td><td>48_Doris</td><td>335.68</td><td>1.66</td><td>189_Pthia</td><td>214.23</td><td>0.83</td></tr><tr><td>14_Irene</td><td>3522.35</td><td>20_Massalia</td><td>333.67</td><td>42.01</td><td>45_Eugenia</td><td>213.72</td><td>6.66</td></tr><tr><td>16_Psyche</td><td>3008.18</td><td>2_Pallas</td><td>288.46</td><td>4.93</td><td>29_Amphitrite</td><td>193.07</td><td>40.43</td></tr><tr><td>9_Metis</td><td>2772.05</td><td>52_Europa</td><td>270.33</td><td>21.66</td><td>9_Metis</td><td>182.47</td><td>15.19</td></tr><tr><td>4_Vesta</td><td>2382.66</td><td>10_Hygiea</td><td>235.19</td><td>1.22</td><td>14_Irene</td><td>180.04</td><td>19.56</td></tr><tr><td>13_Egeria</td><td>2376.39</td><td>14_Irene</td><td>197.95</td><td>17.79</td><td>2_Pallas</td><td>114.65</td><td>12.39</td></tr><tr><td>31_Euphrosyne</td><td>1995.00</td><td>65_Cybele</td><td>184.85</td><td>7.53</td><td>20_Massalia</td><td>104.96</td><td>133.55</td></tr><tr><td>45_Eugenia</td><td>1423.89</td><td>49_Pales</td><td>156.28</td><td>0.65</td><td>52_Europa</td><td>100.20</td><td>58.45</td></tr><tr><td>2_Pallas</td><td>1420.99</td><td>4_Vesta</td><td>138.08</td><td>17.26</td><td>8_Flora</td><td>98.61</td><td>44.84</td></tr><tr><td>65_Cybele</td><td>1391.56</td><td>47_Aglaja</td><td>133.58</td><td>0.15</td><td>65_Cybele</td><td>85.12</td><td>16.35</td></tr><tr><td>21_Lutetia</td><td>1211.49</td><td>45_Eugenia</td><td>129.20</td><td>11.02</td><td>28_Bellona</td><td>81.86</td><td>55.10</td></tr><tr><td>7_Iris</td><td>1128.66</td><td>8_Flora</td><td>103.45</td><td>42.74</td><td>17_Thetis</td><td>72.16</td><td>4.60</td></tr><tr><td>18_Melpomene</td><td>875.10</td><td>1_Ceres</td><td>98.58</td><td>4.94</td><td>10_Hygiea</td><td>70.72</td><td>4.05</td></tr><tr><td>165_Loreley</td><td>785.95</td><td>31_Euphrosyne</td><td>95.83</td><td>20.82</td><td>13_Egeria</td><td>63.56</td><td>37.39</td></tr><tr><td>48_Doris</td><td>557.76</td><td>189_Pthia</td><td>74.72</td><td>2.38</td><td>47_Aglaja</td><td>61.06</td><td>0.32</td></tr><tr><td>1_Ceres</td><td>486.84</td><td>19_Fortuna</td><td>64.90</td><td>186.27</td><td>15_Eunomia</td><td>55.59</td><td>75.54</td></tr><tr><td>130_Elektra</td><td>402.76</td><td>24_Themis</td><td>59.06</td><td>178.32</td><td>4_Vesta</td><td>46.54</td><td>51.19</td></tr><tr><td>17_Thetis</td><td>332.10</td><td>9_Metis</td><td>56.95</td><td>48.67</td><td>49_Pales</td><td>46.45</td><td>2.19</td></tr><tr><td>10_Hygiea</td><td>286.44</td><td>243_Ida</td><td>54.77</td><td>0.29</td><td>243_Ida</td><td>46.23</td><td>0.34</td></tr><tr><td>189_Pthia</td><td>177.82</td><td>130_Elektra</td><td>44.75</td><td>9.00</td><td>130_Elektra</td><td>16.80</td><td>23.98</td></tr><tr><td>121_Hermione</td><td>174.91</td><td>15_Eunomia</td><td>36.12</td><td>116.25</td><td>21_Lutetia</td><td>14.81</td><td>81.82</td></tr><tr><td>11_Parthenope</td><td>111.49</td><td>21_Lutetia</td><td>25.10</td><td>48.27</td><td>121_Hermione</td><td>14.04</td><td>12.45</td></tr><tr><td>22_Kalliope</td><td>107.91</td><td>22_Kalliope</td><td>19.96</td><td>5.41</td><td>253_Mathilde</td><td>12.39</td><td>1.14</td></tr><tr><td>49_Pales</td><td>101.64</td><td>11_Parthenope</td><td>16.64</td><td>6.70</td><td>283_Emma</td><td>10.60</td><td>0.49</td></tr><tr><td>47_Aglaja</td><td>19.69</td><td>253_Mathilde</td><td>15.27</td><td>0.93</td><td>31_Euphrosyne</td><td>8.88</td><td>224.54</td></tr><tr><td>243_Ida</td><td>15.80</td><td>17_Thetis</td><td>15.20</td><td>21.85</td><td>11_Parthenope</td><td>7.98</td><td>13.98</td></tr><tr><td>253_Mathilde</td><td>14.17</td><td>283_Emma</td><td>14.15</td><td>0.37</td><td>107_Camilla</td><td>7.17</td><td>1.93</td></tr><tr><td>107_Camilla</td><td>13.81</td><td>107_Camilla</td><td>13.15</td><td>1.05</td><td>22_Kalliope</td><td>6.90</td><td>15.64</td></tr><tr><td>90_Antiope</td><td>13.21</td><td>121_Hermione</td><td>11.46</td><td>15.26</td><td>90_Antiope</td><td>2.88</td><td>4.58</td></tr><tr><td>87_Sylvia</td><td>7.56</td><td>90_Antiope</td><td>7.92</td><td>1.67</td><td>1_Ceres</td><td>2.43</td><td>200.25</td></tr><tr><td>283_Emma</td><td>5.21</td><td>87_Sylvia</td><td>4.36</td><td>1.73</td><td>87_Sylvia</td><td>1.40</td><td>5.38</td></tr></tbody></table>

Notes. The third column for Mercury and Venus indicates the ratio of the effect with respect to the corresponding effect on Mars.

a limited set of 43 asteroids. But one of the advantages of this database is that it is generated automatically and can be updated each time a new and more accurate determination of the masses of selected asteroids is made or when a new asteroid is included. The ASETEP database enables one to characterize the influence of each asteroid on the orbital motion of the terrestrial planets in an accurate and detailed manner, and will guide the efforts toward improvements of the ephemeris. The case of 111 Ate is particularly significant, because we were able to show that the differences given by two types of ephemerides (DE405 and IPOP08) on the EMB-Mars distance is remarkably similar to the effect caused by the mass uncertainty on this asteroid, which suggests the choice of two different masses for 111 Ate in the ephemeris given above.

Acknowledgements. This work is supported for the first author by a Ph.D. fellowship from the Damsacus University, Syria.

References

Aslan, Z., Gumerov, R., Hudkova, L., et al. 2007, ASP Conf. Ser., 370

Baer, J., URL: http://home.earthlink.net/~jimbaer1/astmass.htm

Baer, J., & Chesley, R. S. 2008a, Cel. Mech. Dyn. Astr., 100, 27

Baer, J., Milani, A., Chesley, S., & Matson, R. D. 2008b, DPS meeting, 40, 493

Bange, J. 1998, A&A, 340, L1

Bange, J. F., & Bec-Borsenberger, A. 1997, ESA, 402, 169

Bretagnon, P., & Francou, G. 1987, A&A, 202, 309

Brouwer, D., & Clemence, G. M. 1961, Methods of Celestial Mechanics (New York: Academic Press)

Chapront-Touzé, M., & Chapront, J. 1988, A&A, 190, 342

Fienga, A., Manche, H., Laskar, J., & Gastineau, M. 2008, A&A, 477, 315
S. Aljbaae and J. Souchay: Effects of the large asteroids on terrestrial planets

Fienga, A., Laskar, J., Morley, T., et al. 2009, A&A, 507, 1675

Hilton, J. L. 1999, AJ, 117, 1077

Ivantsov, A. 2008, Planet. Space Sci., 56, 1857

Kochetova, O. M. 2004, Sol. Syst. Res., 38, 66

Kuchynka, P. 2010, Observatoire de Paris, Ph.D.

Kuchynka, P., Laskar, J., Fienga, A., Manche, H., & Somenzi, L. 2008, Journees Astrometry, Geodynamics and Astronomical Reference Systems, ed. M. Soffel, & N. Capitaine (Dresden)

Kuchynka, P., Laskar, J., Fienga, A., & Manche, H. 2010, A&A, 514, A96

Laskar, J., & Robutel, P. 1995, Cel. Mech. Dyn. Astron., 62, 193

Marchis, F., Descamps, P., Berthier, J., et al. 2008a, Icarus, 195, 295

Marchis, F., Descamps, P., Baek, M., et al. 2008b, Icarus, 196, 97

Mouret, S., Hestroffer, D., & Mignard. F. 2007, A&A, 472, 1017

Mouret, S., Simon, J. L., Mignard, F., & Hestroffer, D. 2009, A&A, 508, 479

Murray, C. D., & Dermott, S. F. 1999, Solar System Dynamics (Cambridge University Press)

Pitjeva, E. V., & Standish, E. M. 2009, Cel. Mech. Dyn. Astron., 103, 365

Simon, J. L. 1987, A&A, 175, 303

Standish, E. M. 1998, IOM, 312.F, 98

Standish, E. M. 2004, A&A, 417, 1165

Standish, E. M., & Fienga, A. 2002, A&A, 384, 322

Standish, E. M., & Hellings, R. W. 1989, Icarus, 80, 326

Tedesco, E. F., Cellino, A., & Zappala, V. 2005, AJ, 129, 2869

Williams, J. G. 1984, Icarus, 57, 1

Williams, J. G. 1989, Asteroid family identifications and proper elements, ed. R. P. Binzel, T. Gehrels, M. Matthews (Asteroids II, Univ. of Arizona Press), 1034

Yeomans, D. K., Barriot, J. P., Dunham, D. W., et al. 1997, Science, 278, 2106
Appendix A: Effects on the Mars orbital motion

<figure><img src="image_25.png" /><figcaption>**Fig. A.1.** Effects of Ceres on $e_{Mars}$ for a 100 y time interval from JJ 2009.81; raw signal (in black), fitting curve (in red) and residuals (in green).</figcaption></figure>

<figure><img src="image_26.png" /><figcaption>**Fig. A.2.** Effects of Ceres on $i_{Mars}$ for a 100 y time interval from JJ 2009.81; raw signal (in black), fitting curve (in red) and residuals (in green).</figcaption></figure>

<figure><img src="image_27.png" /><figcaption>**Fig. A.3.** Effects of Ceres on $\Omega_{Mars}$ for a 100 y time interval from JJ 2009.81; raw signal (in black), fitting curve (in red) and residuals (in green).</figcaption></figure>

<figure><img src="image_28.png" /><figcaption>**Fig. A.4.** Effects of Ceres on $\dot{\omega}_{Mars}$ for a 100 y time interval from JJ 2009.81; raw signal (in black), fitting curve (in red) and residuals (in green).</figcaption></figure>
<table><caption>**Table A.1.** Coefficients of Fourier and Poisson series for the effects of Ceres on $e_{Mars}$.</caption><thead><tr><th>period<br>(Day)</th><th>period<br>(Year)</th><th>SIN<br>×10<sup>-10</sup></th><th>COS<br>×10<sup>-10</sup></th><th>×10<sup>-10</sup></th><th>T SIN<br>×10<sup>-14</sup></th><th>T COS<br>×10<sup>-14</sup></th><th>×10<sup>-14</sup></th><th>Amplitude<br>×10<sup>-8</sup></th><th>Asteroid</th></tr></thead><tbody><tr><td>16800.43</td><td>46.00</td><td>71.042</td><td>-57.060</td><td>91.120</td><td>-336.820</td><td>219.769</td><td>310.801</td><td>0.145</td><td>l_Ceres</td></tr><tr><td>4862.41</td><td>13.31</td><td>0.650</td><td>4.085</td><td>4.137</td><td>1.708</td><td>-18.545</td><td>26.227</td><td>0.064</td><td>l_Ceres</td></tr><tr><td>3752.63</td><td>10.27</td><td>-5.746</td><td>-12.352</td><td>13.623</td><td>-14.963</td><td>26.378</td><td>37.305</td><td>0.127</td><td>l_Ceres</td></tr><tr><td>3046.97</td><td>8.34</td><td>2.753</td><td>-1.422</td><td>3.099</td><td>8.535</td><td>4.919</td><td>6.957</td><td>0.046</td><td>l_Ceres</td></tr><tr><td>1875.99</td><td>5.14</td><td>-4.305</td><td>-4.983</td><td>6.585</td><td>11.563</td><td>32.195</td><td>45.530</td><td>0.028</td><td>l_Ceres</td></tr><tr><td>1682.02</td><td>4.61</td><td>-4.362</td><td>-1.963</td><td>4.783</td><td>26.924</td><td>7.169</td><td>10.138</td><td>0.013</td><td>l_Ceres</td></tr><tr><td>1250.69</td><td>3.42</td><td>-1.101</td><td>7.132</td><td>7.217</td><td>6.247</td><td>-31.438</td><td>44.460</td><td>0.009</td><td>l_Ceres</td></tr><tr><td>1083.67</td><td>2.97</td><td>1.155</td><td>-0.487</td><td>1.254</td><td>-7.228</td><td>0.406</td><td>0.574</td><td>0.007</td><td>l_Ceres</td></tr><tr><td>886.91</td><td>2.43</td><td>-3.654</td><td>0.482</td><td>3.686</td><td>22.982</td><td>4.080</td><td>5.769</td><td>0.018</td><td>l_Ceres</td></tr><tr><td>840.95</td><td>2.30</td><td>-0.706</td><td>-0.207</td><td>0.736</td><td>0.831</td><td>-0.538</td><td>0.761</td><td>0.007</td><td>l_Ceres</td></tr><tr><td>717.34</td><td>1.96</td><td>1.234</td><td>1.552</td><td>1.983</td><td>-2.558</td><td>-8.717</td><td>12.328</td><td>0.009</td><td>l_Ceres</td></tr><tr><td>687.04</td><td>1.88</td><td>0.063</td><td>-0.479</td><td>0.483</td><td>-3.217</td><td>2.391</td><td>3.381</td><td>0.006</td><td>l_Ceres</td></tr><tr><td>502.91</td><td>1.38</td><td>1.721</td><td>-0.327</td><td>1.752</td><td>-8.277</td><td>-1.340</td><td>1.894</td><td>0.006</td><td>l_Ceres</td></tr></tbody><tfoot><tr><td colspan="10">DISPERSION Before = 0.437036 ×10<sup>-8</sup> After = 0.449648 ×10<sup>-9</sup></td></tr></tfoot></table>

**Notes.** Bias: $-0.147732 \times 10^{-5}$; linear: $0.767680 \times 10^{-9}$; $T^2$: $-0.246106 \times 10^{-13}$. The next-to-last column indicates the peak-to-peak amplitude of the total mixed component for the 100 y time interval.

<table><caption>**Table A.2.** Coefficients of Fourier and Poisson series for the effects of Ceres on $i_{Mars}$.</caption><thead><tr><th>Period<br>(Day)</th><th>Period<br>(Year)</th><th>SIN<br>×10<sup>-5</sup></th><th>COS<br>×10<sup>-5</sup></th><th>×10<sup>-5</sup></th><th>T SIN<br>×10<sup>-9</sup></th><th>T COS<br>×10<sup>-9</sup></th><th>×10<sup>-9</sup></th><th>Amplitude<br>×10<sup>-4</sup></th><th>Asteroid</th></tr></thead><tbody><tr><td>16704.15</td><td>45.73</td><td>-30.404</td><td>-41.976</td><td>51.830</td><td>128.455</td><td>201.120</td><td>284.427</td><td>0.354</td><td>l_Ceres</td></tr><tr><td>4871.57</td><td>13.34</td><td>-3.481</td><td>-0.198</td><td>3.486</td><td>16.174</td><td>2.520</td><td>3.564</td><td>0.081</td><td>l_Ceres</td></tr><tr><td>3745.29</td><td>10.25</td><td>2.366</td><td>2.832</td><td>3.691</td><td>-6.410</td><td>-15.360</td><td>21.722</td><td>0.143</td><td>l_Ceres</td></tr><tr><td>3043.72</td><td>8.33</td><td>0.387</td><td>0.930</td><td>1.007</td><td>-1.532</td><td>-9.749</td><td>13.788</td><td>0.111</td><td>l_Ceres</td></tr><tr><td>1878.41</td><td>5.14</td><td>0.518</td><td>0.763</td><td>0.922</td><td>-1.876</td><td>-5.321</td><td>7.525</td><td>0.049</td><td>l_Ceres</td></tr><tr><td>1683.67</td><td>4.61</td><td>1.111</td><td>0.683</td><td>1.304</td><td>-7.355</td><td>-4.239</td><td>5.995</td><td>0.052</td><td>l_Ceres</td></tr><tr><td>1251.11</td><td>3.43</td><td>-1.761</td><td>0.037</td><td>1.762</td><td>8.127</td><td>0.333</td><td>0.471</td><td>0.026</td><td>l_Ceres</td></tr><tr><td>1161.41</td><td>3.18</td><td>1.597</td><td>0.156</td><td>1.604</td><td>-5.001</td><td>-0.533</td><td>0.754</td><td>0.067</td><td>l_Ceres</td></tr><tr><td>1083.47</td><td>2.97</td><td>0.253</td><td>0.348</td><td>0.430</td><td>-2.531</td><td>-1.831</td><td>2.590</td><td>0.027</td><td>l_Ceres</td></tr><tr><td>887.20</td><td>2.43</td><td>0.772</td><td>-1.393</td><td>1.593</td><td>-2.531</td><td>7.513</td><td>10.625</td><td>0.043</td><td>l_Ceres</td></tr><tr><td>840.94</td><td>2.30</td><td>-0.626</td><td>0.521</td><td>0.814</td><td>0.989</td><td>-3.728</td><td>5.273</td><td>0.050</td><td>l_Ceres</td></tr><tr><td>687.10</td><td>1.88</td><td>-0.399</td><td>0.843</td><td>0.932</td><td>1.609</td><td>-5.043</td><td>7.132</td><td>0.024</td><td>l_Ceres</td></tr><tr><td>580.74</td><td>1.59</td><td>-0.469</td><td>-0.210</td><td>0.514</td><td>0.789</td><td>0.457</td><td>0.647</td><td>0.038</td><td>l_Ceres</td></tr><tr><td>560.56</td><td>1.53</td><td>-0.124</td><td>-0.154</td><td>0.198</td><td>1.263</td><td>1.183</td><td>1.673</td><td>0.017</td><td>l_Ceres</td></tr><tr><td>503.00</td><td>1.38</td><td>-0.861</td><td>0.430</td><td>0.963</td><td>3.812</td><td>-2.672</td><td>3.779</td><td>0.020</td><td>l_Ceres</td></tr><tr><td>487.77</td><td>1.34</td><td>0.452</td><td>-0.048</td><td>0.454</td><td>-1.446</td><td>1.033</td><td>1.462</td><td>0.023</td><td>l_Ceres</td></tr><tr><td>387.16</td><td>1.06</td><td>-0.661</td><td>-0.235</td><td>0.701</td><td>3.708</td><td>1.657</td><td>2.344</td><td>0.017</td><td>l_Ceres</td></tr></tbody><tfoot><tr><td colspan="10">AMPLITUDE Before = 0.000024 After = 0.237359 ×10<sup>-4</sup></td></tr></tfoot></table>

**Notes.** Bias: $0.258949 \times 10^{-1}$; linear: $-0.251401 \times 10^{-4}$; $T^2$: $0.600592 \times 10^{-8}$. The next-to-last column indicates the peak-to-peak amplitude of the total mixed component for the 100 y time interval.
<table><caption>**Table A.3.** Coefficients of Fourier and Poisson series for the effects of Ceres on $\Omega_{\text{Mars}}$.</caption><thead><tr><th>Period (Day)</th><th>Period (Year)</th><th>SIN ×10<sup>-4</sup></th><th>COS ×10<sup>-4</sup></th><th>×10<sup>-4</sup></th><th>T SIN ×10<sup>-7</sup></th><th>T COS ×10<sup>-7</sup></th><th>×10<sup>-7</sup></th><th>Amplitude ×10<sup>-2</sup></th><th>Asteroid</th></tr></thead><tbody><tr><td>17237.30</td><td>47.19</td><td>-254.868</td><td>163.129</td><td>302.603</td><td>120.532</td><td>-84.520</td><td>119.529</td><td>0.126</td><td>1_Ceres</td></tr><tr><td>10945.54</td><td>29.97</td><td>3.333</td><td>13.201</td><td>13.616</td><td>-1.550</td><td>-6.445</td><td>9.114</td><td>0.116</td><td>1_Ceres</td></tr><tr><td>7586.87</td><td>20.77</td><td>-3.312</td><td>-4.891</td><td>5.907</td><td>1.250</td><td>2.352</td><td>3.326</td><td>0.054</td><td>1_Ceres</td></tr><tr><td>4921.78</td><td>13.48</td><td>7.651</td><td>-0.338</td><td>7.658</td><td>-3.232</td><td>0.853</td><td>1.207</td><td>0.036</td><td>1_Ceres</td></tr><tr><td>3754.45</td><td>10.28</td><td>1.482</td><td>-1.680</td><td>2.241</td><td>2.192</td><td>0.170</td><td>0.241</td><td>0.078</td><td>1_Ceres</td></tr><tr><td>3045.26</td><td>8.34</td><td>0.149</td><td>5.993</td><td>5.995</td><td>-2.205</td><td>-4.246</td><td>6.005</td><td>0.056</td><td>1_Ceres</td></tr><tr><td>1877.72</td><td>5.14</td><td>-6.977</td><td>1.328</td><td>7.102</td><td>4.562</td><td>-0.518</td><td>0.733</td><td>0.030</td><td>1_Ceres</td></tr><tr><td>1681.32</td><td>4.60</td><td>-0.203</td><td>-3.210</td><td>3.216</td><td>-1.143</td><td>1.428</td><td>2.019</td><td>0.029</td><td>1_Ceres</td></tr><tr><td>1521.01</td><td>4.16</td><td>3.786</td><td>-0.861</td><td>3.882</td><td>-1.957</td><td>0.150</td><td>0.212</td><td>0.006</td><td>1_Ceres</td></tr><tr><td>1252.35</td><td>3.43</td><td>-0.665</td><td>0.300</td><td>0.729</td><td>0.284</td><td>0.253</td><td>0.357</td><td>0.014</td><td>1_Ceres</td></tr><tr><td>1161.48</td><td>3.18</td><td>-2.464</td><td>1.778</td><td>3.038</td><td>0.017</td><td>-1.536</td><td>2.172</td><td>0.031</td><td>1_Ceres</td></tr><tr><td>1083.21</td><td>2.97</td><td>-1.423</td><td>-0.137</td><td>1.429</td><td>1.120</td><td>0.322</td><td>0.456</td><td>0.010</td><td>1_Ceres</td></tr><tr><td>887.02</td><td>2.43</td><td>5.332</td><td>-2.557</td><td>5.913</td><td>-3.223</td><td>1.151</td><td>1.628</td><td>0.018</td><td>1_Ceres</td></tr><tr><td>840.68</td><td>2.30</td><td>-0.150</td><td>-0.303</td><td>0.338</td><td>0.904</td><td>0.412</td><td>0.582</td><td>0.019</td><td>1_Ceres</td></tr><tr><td>717.62</td><td>1.96</td><td>3.271</td><td>1.048</td><td>3.435</td><td>-1.827</td><td>-0.739</td><td>1.045</td><td>0.009</td><td>1_Ceres</td></tr><tr><td>686.93</td><td>1.88</td><td>0.479</td><td>-5.096</td><td>5.119</td><td>0.187</td><td>2.734</td><td>3.866</td><td>0.012</td><td>1_Ceres</td></tr><tr><td>580.74</td><td>1.59</td><td>-3.886</td><td>-1.628</td><td>4.214</td><td>2.118</td><td>1.406</td><td>1.989</td><td>0.016</td><td>1_Ceres</td></tr><tr><td>560.52</td><td>1.53</td><td>0.764</td><td>1.457</td><td>1.645</td><td>-0.390</td><td>-0.971</td><td>1.374</td><td>0.006</td><td>1_Ceres</td></tr><tr><td>502.94</td><td>1.38</td><td>-1.693</td><td>-1.005</td><td>1.969</td><td>1.027</td><td>0.679</td><td>0.960</td><td>0.007</td><td>1_Ceres</td></tr><tr><td>487.69</td><td>1.34</td><td>-2.055</td><td>1.376</td><td>2.473</td><td>0.739</td><td>-0.930</td><td>1.316</td><td>0.008</td><td>1_Ceres</td></tr><tr><td>431.64</td><td>1.18</td><td>2.216</td><td>3.012</td><td>3.739</td><td>-1.231</td><td>-1.633</td><td>2.309</td><td>0.006</td><td>1_Ceres</td></tr><tr><td>387.16</td><td>1.06</td><td>0.068</td><td>3.096</td><td>3.097</td><td>-0.005</td><td>-1.781</td><td>2.519</td><td>0.006</td><td>1_Ceres</td></tr><tr><td>343.47</td><td>0.94</td><td>-1.495</td><td>-3.287</td><td>3.611</td><td>0.833</td><td>1.882</td><td>2.661</td><td>0.008</td><td>1_Ceres</td></tr></tbody><tfoot><tr><td>AMPLITUDE</td><td colspan="4">Before = 0.005407</td><td colspan="4">After = 0.000974</td><td></td></tr></tfoot></table>

**Notes.** Bias: $-0.811114$; linear: $0.427426 \times 10^{-3}$; $T^2$: $-0.121763 \times 10^{-7}$. The next-to-last column indicates the peak-to-peak amplitude of the total mixed component for the 100 y time interval.

<table><caption>**Table A.4.** Coefficients of Fourier and Poisson series for the effects of Ceres on $\tilde{\omega}_{\text{Mars}}$.</caption><thead><tr><th>Period (Day)</th><th>Period (Year)</th><th>SIN ×10<sup>-4</sup></th><th>COS ×10<sup>-4</sup></th><th>×10<sup>-4</sup></th><th>T SIN ×10<sup>-7</sup></th><th>T COS ×10<sup>-7</sup></th><th>×10<sup>-7</sup></th><th>Amplitude ×10<sup>-2</sup></th><th>Asteroid</th></tr></thead><tbody><tr><td>17112.22</td><td>46.85</td><td>-156.276</td><td>611.164</td><td>630.828</td><td>70.223</td><td>-290.043</td><td>410.183</td><td>0.311</td><td>1_Ceres</td></tr><tr><td>10842.14</td><td>29.68</td><td>276.449</td><td>3.355</td><td>276.470</td><td>-134.342</td><td>-2.937</td><td>4.153</td><td>0.378</td><td>1_Ceres</td></tr><tr><td>4924.23</td><td>13.48</td><td>0.168</td><td>-17.541</td><td>17.542</td><td>0.160</td><td>9.637</td><td>13.629</td><td>0.115</td><td>1_Ceres</td></tr><tr><td>3755.58</td><td>10.28</td><td>25.200</td><td>-2.573</td><td>25.331</td><td>-0.378</td><td>2.019</td><td>2.855</td><td>0.273</td><td>1_Ceres</td></tr><tr><td>3046.85</td><td>8.34</td><td>-7.000</td><td>-5.383</td><td>8.831</td><td>2.778</td><td>-2.097</td><td>2.965</td><td>0.099</td><td>1_Ceres</td></tr><tr><td>1877.78</td><td>5.14</td><td>-2.572</td><td>-2.221</td><td>3.398</td><td>3.756</td><td>1.278</td><td>1.808</td><td>0.060</td><td>1_Ceres</td></tr><tr><td>1681.96</td><td>4.60</td><td>-2.301</td><td>12.463</td><td>12.674</td><td>0.668</td><td>-7.302</td><td>10.327</td><td>0.031</td><td>1_Ceres</td></tr><tr><td>1251.92</td><td>3.43</td><td>-5.802</td><td>-7.270</td><td>9.301</td><td>3.454</td><td>3.746</td><td>5.298</td><td>0.014</td><td>1_Ceres</td></tr><tr><td>1083.69</td><td>2.97</td><td>6.023</td><td>-8.188</td><td>10.165</td><td>-3.391</td><td>4.274</td><td>6.044</td><td>0.012</td><td>1_Ceres</td></tr><tr><td>887.10</td><td>2.43</td><td>6.554</td><td>1.009</td><td>6.632</td><td>-4.692</td><td>-1.180</td><td>1.669</td><td>0.039</td><td>1_Ceres</td></tr><tr><td>840.86</td><td>2.30</td><td>0.293</td><td>0.606</td><td>0.673</td><td>-0.137</td><td>0.308</td><td>0.435</td><td>0.012</td><td>1_Ceres</td></tr><tr><td>717.65</td><td>1.96</td><td>3.002</td><td>2.539</td><td>3.931</td><td>-1.857</td><td>-1.823</td><td>2.578</td><td>0.017</td><td>1_Ceres</td></tr><tr><td>687.08</td><td>1.88</td><td>4.640</td><td>3.120</td><td>5.591</td><td>-2.526</td><td>-1.043</td><td>1.474</td><td>0.012</td><td>1_Ceres</td></tr><tr><td>502.96</td><td>1.38</td><td>0.513</td><td>-1.741</td><td>1.815</td><td>0.007</td><td>1.204</td><td>1.703</td><td>0.010</td><td>1_Ceres</td></tr></tbody><tfoot><tr><td>AMPLITUDE</td><td colspan="4">Before = 0.009966</td><td colspan="4">After = 0.000944</td><td></td></tr></tfoot></table>

**Notes.** Bias: $-0.605169 \times 10^{+2}$; linear: $-0.121834 \times 10^{-1}$; $T^2$: $0.156335 \times 10^{-5}$. The next-to-last column indicates the peak-to-peak amplitude of the total mixed component for the 100 y time interval.
<table><caption>**Table A.5.** Coefficients of Fourier and Poisson series for the effects of Pallas on $a_{\text{Mars}}$.</caption><thead><tr><th>Period (Day)</th><th>Period (Year)</th><th>SIN ×10<sup>0</sup></th><th>COS ×10<sup>0</sup></th><th>×10<sup>0</sup></th><th>T SIN ×10<sup>-3</sup></th><th>T COS ×10<sup>-3</sup></th><th>×10<sup>-3</sup></th><th>Amplitude ×10<sup>1</sup></th><th>Asteroid</th></tr></thead><tbody><tr><td>21839.28</td><td>59.79</td><td>-10.566</td><td>1.282</td><td>10.644</td><td>3.470</td><td>-0.149</td><td>0.210</td><td>1.031</td><td>2_Pallas</td></tr><tr><td>17918.54</td><td>49.06</td><td>-107.297</td><td>3.679</td><td>107.360</td><td>16.781</td><td>-2.048</td><td>2.896</td><td>7.353</td><td>2_Pallas</td></tr><tr><td>9865.09</td><td>27.01</td><td>-9.001</td><td>4.336</td><td>9.991</td><td>2.672</td><td>-0.949</td><td>1.342</td><td>0.446</td><td>2_Pallas</td></tr><tr><td>8955.04</td><td>24.52</td><td>12.925</td><td>19.794</td><td>23.640</td><td>-3.318</td><td>-3.703</td><td>5.236</td><td>1.412</td><td>2_Pallas</td></tr><tr><td>6364.17</td><td>17.42</td><td>5.432</td><td>9.561</td><td>10.996</td><td>-0.971</td><td>-1.935</td><td>2.736</td><td>0.621</td><td>2_Pallas</td></tr><tr><td>5974.81</td><td>16.36</td><td>4.910</td><td>-2.340</td><td>5.439</td><td>-0.825</td><td>1.045</td><td>1.478</td><td>0.290</td><td>2_Pallas</td></tr><tr><td>4695.27</td><td>12.85</td><td>13.748</td><td>-17.279</td><td>22.081</td><td>-2.510</td><td>2.891</td><td>4.088</td><td>1.480</td><td>2_Pallas</td></tr><tr><td>3720.02</td><td>10.18</td><td>-3.842</td><td>-4.707</td><td>6.075</td><td>-0.432</td><td>0.013</td><td>0.018</td><td>0.680</td><td>2_Pallas</td></tr><tr><td>3080.06</td><td>8.43</td><td>-2.749</td><td>-5.035</td><td>5.736</td><td>0.396</td><td>-0.555</td><td>0.785</td><td>0.699</td><td>2_Pallas</td></tr><tr><td>2628.38</td><td>7.20</td><td>-8.639</td><td>-1.328</td><td>8.740</td><td>1.655</td><td>0.137</td><td>0.194</td><td>0.560</td><td>2_Pallas</td></tr><tr><td>2075.66</td><td>5.68</td><td>1.660</td><td>7.042</td><td>7.236</td><td>-0.231</td><td>-1.369</td><td>1.937</td><td>0.488</td><td>2_Pallas</td></tr><tr><td>1860.14</td><td>5.09</td><td>15.833</td><td>-9.731</td><td>18.584</td><td>-2.247</td><td>1.721</td><td>2.434</td><td>1.301</td><td>2_Pallas</td></tr><tr><td>1540.23</td><td>4.22</td><td>3.956</td><td>-4.351</td><td>5.880</td><td>-0.578</td><td>0.830</td><td>1.173</td><td>0.385</td><td>2_Pallas</td></tr><tr><td>1240.07</td><td>3.40</td><td>0.060</td><td>9.006</td><td>9.006</td><td>-0.144</td><td>-1.609</td><td>2.275</td><td>0.579</td><td>2_Pallas</td></tr><tr><td>930.06</td><td>2.55</td><td>-3.330</td><td>-2.179</td><td>3.980</td><td>0.558</td><td>0.278</td><td>0.393</td><td>0.298</td><td>2_Pallas</td></tr><tr><td>884.18</td><td>2.42</td><td>-8.246</td><td>9.304</td><td>12.432</td><td>0.721</td><td>-1.233</td><td>1.744</td><td>0.966</td><td>2_Pallas</td></tr><tr><td>714.39</td><td>1.96</td><td>-2.652</td><td>-6.204</td><td>6.747</td><td>0.530</td><td>0.901</td><td>1.274</td><td>0.464</td><td>2_Pallas</td></tr><tr><td>579.91</td><td>1.59</td><td>1.465</td><td>-5.380</td><td>5.576</td><td>0.679</td><td>-0.289</td><td>0.408</td><td>0.718</td><td>2_Pallas</td></tr><tr><td>501.71</td><td>1.37</td><td>3.219</td><td>3.088</td><td>4.461</td><td>-0.590</td><td>-0.005</td><td>0.008</td><td>0.368</td><td>2_Pallas</td></tr><tr><td>442.08</td><td>1.21</td><td>-2.850</td><td>0.973</td><td>3.011</td><td>0.416</td><td>-0.174</td><td>0.246</td><td>0.215</td><td>2_Pallas</td></tr><tr><td>408.23</td><td>1.12</td><td>-10.042</td><td>36.276</td><td>37.640</td><td>4.874</td><td>-18.183</td><td>25.715</td><td>1.902</td><td>2_Pallas</td></tr><tr><td>386.59</td><td>1.06</td><td>-0.348</td><td>0.504</td><td>0.613</td><td>-0.644</td><td>-0.204</td><td>0.288</td><td>0.230</td><td>2_Pallas</td></tr><tr><td>272.16</td><td>0.75</td><td>-15.594</td><td>-2.638</td><td>15.815</td><td>7.718</td><td>1.194</td><td>1.688</td><td>0.770</td><td>2_Pallas</td></tr></tbody><tfoot><tr><td colspan="10">AMPLITUDE Before = 214.335797 After = 72.316238</td></tr></tfoot></table>

**Notes.** Bias: $0.398536 \times 10^{+4}$; linear: $-0.326098 \times 10^{+1}$; $T^2$: $0.651096 \times 10^{-3}$. The next-to-last column indicates the peak-to-peak amplitude of the total mixed component for the 1000 y time interval considered.

<table><caption>**Table A.6.** Coefficients of Fourier and Poisson series for the effects of Vesta on $a_{\text{mars}}$.</caption><thead><tr><th>Period (Day)</th><th>Period (Year)</th><th>SIN ×10<sup>0</sup></th><th>COS ×10<sup>0</sup></th><th>×10<sup>0</sup></th><th>T SIN ×10<sup>-3</sup></th><th>T COS ×10<sup>-3</sup></th><th>×10<sup>-3</sup></th><th>Amplitude ×10<sup>2</sup></th><th>Asteroid</th></tr></thead><tbody><tr><td>18867.98</td><td>51.66</td><td>83.834</td><td>-107.509</td><td>136.331</td><td>-3.213</td><td>4.449</td><td>6.291</td><td>1.267</td><td>4_Vesta</td></tr><tr><td>9443.39</td><td>25.85</td><td>17.064</td><td>-41.235</td><td>44.626</td><td>-0.487</td><td>2.848</td><td>4.027</td><td>0.389</td><td>4_Vesta</td></tr><tr><td>6296.23</td><td>17.24</td><td>3.447</td><td>-4.872</td><td>5.968</td><td>1.009</td><td>-0.857</td><td>1.211</td><td>0.100</td><td>4_Vesta</td></tr><tr><td>1679.56</td><td>4.60</td><td>0.997</td><td>12.655</td><td>12.695</td><td>0.599</td><td>-1.360</td><td>1.923</td><td>0.105</td><td>4_Vesta</td></tr><tr><td>1542.42</td><td>4.22</td><td>0.669</td><td>28.454</td><td>28.462</td><td>-0.482</td><td>-1.519</td><td>2.149</td><td>0.255</td><td>4_Vesta</td></tr><tr><td>1425.93</td><td>3.90</td><td>-1.302</td><td>10.131</td><td>10.215</td><td>-0.384</td><td>0.443</td><td>0.627</td><td>0.115</td><td>4_Vesta</td></tr><tr><td>740.95</td><td>2.03</td><td>-0.607</td><td>13.687</td><td>13.700</td><td>-0.711</td><td>-0.238</td><td>0.337</td><td>0.140</td><td>4_Vesta</td></tr><tr><td>712.97</td><td>1.95</td><td>-8.158</td><td>30.122</td><td>31.208</td><td>-1.160</td><td>-0.688</td><td>0.972</td><td>0.307</td><td>4_Vesta</td></tr><tr><td>487.59</td><td>1.33</td><td>-1.036</td><td>10.484</td><td>10.535</td><td>-1.058</td><td>-1.360</td><td>1.924</td><td>0.083</td><td>4_Vesta</td></tr><tr><td>475.32</td><td>1.30</td><td>-5.351</td><td>15.041</td><td>15.964</td><td>-1.000</td><td>-0.548</td><td>0.775</td><td>0.158</td><td>4_Vesta</td></tr><tr><td>408.23</td><td>1.12</td><td>30.965</td><td>-107.157</td><td>111.541</td><td>-14.906</td><td>53.474</td><td>75.624</td><td>0.555</td><td>4_Vesta</td></tr><tr><td>356.49</td><td>0.98</td><td>-2.836</td><td>7.467</td><td>7.988</td><td>-1.067</td><td>-0.451</td><td>0.638</td><td>0.086</td><td>4_Vesta</td></tr><tr><td>272.16</td><td>0.75</td><td>40.158</td><td>30.629</td><td>50.506</td><td>-19.865</td><td>-14.742</td><td>20.848</td><td>0.240</td><td>4_Vesta</td></tr></tbody><tfoot><tr><td colspan="10">AMPLITUDE Before = 545.913598 After = 146.660553</td></tr></tfoot></table>

**Notes.** Bias: $0.152498 \times 10^{+3}$; linear: $-0.209840 \times 10^{-1}$; $T^2$: $0.445801 \times 10^{-5}$. The next-to-last column indicates the peak-to-peak amplitude of the total mixed component for the 1000 y time interval considered.
## Appendix B: individual and combined asteroidal effects

<figure><img src="image_29.png" /></figure>

**Fig. B.1.** Combined effects of the 43 asteroids on $e_{Mars}$. The curve in green represents the difference between two ways of computation: adding individual perturbations, or making a combined integration.

<figure><img src="image_30.png" /></figure>

**Fig. B.2.** Zoom of the green curve of Fig. B.1, showing the difference between two ways of computation.

<figure><img src="image_31.png" /></figure>

**Fig. B.3.** Combined effects of the 43 asteroids on $i_{Mars}$. The curve in green represents the difference between two ways of computation: adding individual perturbations, or making a combined integration.

<figure><img src="image_32.png" /></figure>

**Fig. B.4.** Zoom of the green curve of Fig. B.3, showing the difference between two ways of computation.

<figure><img src="image_33.png" /></figure>

**Fig. B.5.** Combined effects of the 43 asteroids on $\Omega_{Mars}$. The curve in green represents the difference between two ways of computation: adding individual perturbations, or making a combined integration.

<figure><img src="image_34.png" /></figure>

**Fig. B.6.** Zoom of the green curve of Fig. B.5, showing the difference between two ways of computation.
<figure><img src="image_35.png" /></figure>

**Fig. B.7.** Combined effects of the 43 asteroids on $\tilde{\omega}_{\text{Mars}}$. The curve in green represents the difference between two ways of computation: adding individual perturbations, or making a combined integration.

<figure><img src="image_36.png" /></figure>

**Fig. B.8.** Zoom of the green curve of Fig. B.7, showing the difference between two ways of computation.

<figure><img src="image_37.png" /></figure>

**Fig. B.9.** Combined effects of the 43 asteroids on $\lambda_{\text{Mars}}$. The curve in green represents the difference between two ways of computation: adding individual perturbations, or making a combined integration.

<figure><img src="image_38.png" /></figure>

**Fig. B.10.** Zoom of the green curve of Fig. B.9, showing the difference between two ways of computation.

## Appendix C: Use of ASTEP for ephemerides study

<figure><img src="image_39.png" /></figure>

**Fig. C.1.** Uncertainty of the EMB-to-Mars distance due to the lack of accuracy in Pallas' mass (red curve) and difference in the EMB-to-Mars distance between DE405 and INPOP08 (black curve).

<figure><img src="image_40.png" /></figure>

**Fig. C.2.** Uncertainty of the EMB-to-Mars distance due to the lack of accuracy in Vesta's mass (red curve) and difference in the EMB-to-Mars distance between DE405 and INPOP08 (black curve).