Secular evolution of asteroid orbits

Let us now consider the perturbing influence of the planets on the orbit of an asteroid. Because asteroids have much smaller masses than planets, it is reasonable to suppose that the perturbing influence of the asteroid on the planetary orbits is negligible. Let the asteroid have the standard osculating orbital elements $a$, $\skew{5}\bar{\lambda}_0$, $e$, $I$, $\varpi $, ${\mit \Omega }$, and the alternative elements $h=e\,\sin\varpi$, $k=e\,\cos\varpi$, $p=\sin I\,\sin {\mit\Omega}$, and $q=\sin I\,\cos{\mit\Omega}$. Thus, the mean orbital angular velocity of the asteroid is $n=(G\,M/a^{\,3})^{1/2}$, where $M$ is the solar mass. Likewise, let the eight planets have the standard osculating orbital elements $a_i$, $\skew{5}\bar{\lambda}_{0\,i}$, $e_i$, $I_i$, $\varpi_i$, ${\mit\Omega}_i$, and the alternative elements $h_i=e_i\,\sin\varpi_i$, $k_i=e_i\,\cos\varpi_i$, $p_i=\sin I_i\,\sin{\mit\Omega}_i$, and $q_i=\sin I_i\,\cos{\mit\Omega}_i$, for $i=1,8$. It is helpful to define the following parameters:

\begin{displaymath}\alpha_{i} =\left\{
\begin{array}{ccc}
a/a_i&\mbox{\hspace{1cm}}&a_i>a\\ [0.5ex]
a_i/a&&a_i<a
\end{array}\right.,\end{displaymath} (10.85)

and

\begin{displaymath}\bar{\alpha}_{i} =\left\{
\begin{array}{ccc}
a/a_i&\mbox{\hspace{1cm}}&a_i>a\\ [0.5ex]
1&&a_i<a
\end{array}\right.,\end{displaymath} (10.86)

as well as

$\displaystyle \epsilon_i = \frac{m_i}{M}.$ (10.87)

Figure: 10.6 The osculating eccentricity plotted against the sine of the osculating inclination (relative to J2000 equinox and ecliptic) for the orbits of the first 100,000 numbered asteroids at MJD 55400. Raw data from JPL Small-Body Database.
\includegraphics[height=4in]{Chapter09/fig9_06.eps}

By analogy with the analysis in the previous section, the secular terms in the disturbing function of the asteroid, generated by the perturbing influence of the planets, cause the asteroid's osculating orbital elements to evolve in time as

$\displaystyle \frac{dh}{dt}$ $\displaystyle = A\,k + \sum_{i=1,8} A_i\,k_i,$ (10.88)
$\displaystyle \frac{dk}{dt}$ $\displaystyle = -A\,h -\sum_{i=1,8}A_i\,h_i,$ (10.89)
$\displaystyle \frac{dp}{dt}$ $\displaystyle = B\,q + \sum_{i=1,8} B_i\,q_i,$ (10.90)
$\displaystyle \frac{dq}{dt}$ $\displaystyle = -B\,p -\sum_{i=1,8}B_i\,p_i,$ (10.91)

where

$\displaystyle A$ $\displaystyle = \sum_{i=1,8} \frac{n}{4}\,\epsilon_i\,\alpha_i\,\bar{\alpha}_i\,b_{3/2}^{(1)}(\alpha_i),$ (10.92)
$\displaystyle A_i$ $\displaystyle = - \frac{n}{4}\,\epsilon_i\,\alpha_i\,\bar{\alpha}_i\,b_{3/2}^{(2)}(\alpha_i),$ (10.93)
$\displaystyle B$ $\displaystyle =- \sum_{i=1,8} \frac{n}{4}\,\epsilon_i\,\alpha_i\,\bar{\alpha}_i\,b_{3/2}^{(1)}(\alpha_i),$ (10.94)
$\displaystyle B_i$ $\displaystyle = \frac{n}{4}\,\epsilon_i\,\alpha_i\,\bar{\alpha}_i\,b_{3/2}^{(1)}(\alpha_i).$ (10.95)

However, as we have already seen, the planetary osculating elements themselves evolve in time as

$\displaystyle h_i$ $\displaystyle =\sum_{l=1,8} e_{il}\,\sin(g_l\,t+\beta_l),$ (10.96)
$\displaystyle k_i$ $\displaystyle =\sum_{l=1,8} e_{il}\,\cos(g_l\,t+\beta_l),$ (10.97)
$\displaystyle p_i$ $\displaystyle =\sum_{l=1,8} I_{il}\,\sin(f_l\,t+\gamma_l),$ (10.98)
$\displaystyle q_i$ $\displaystyle =\sum_{l=1,8} I_{il}\,\cos(f_l\,t+\gamma_l).$ (10.99)

Equations (10.88)–(10.91) can be solved to give

$\displaystyle h(t)$ $\displaystyle = e_{\rm free}\,\sin(A\,t+\beta_{\rm free})+ h_{\rm forced}(t),$ (10.100)
$\displaystyle k(t)$ $\displaystyle = e_{\rm free}\,\cos(A\,t+\beta_{\rm free})+ k_{\rm forced}(t),$ (10.101)
$\displaystyle p(t)$ $\displaystyle = \sin I_{\rm free}\,\sin(B\,t+\gamma_{\rm free})+ p_{\rm forced}(t),$ (10.102)
$\displaystyle q(t)$ $\displaystyle = \sin I_{\rm free}\,\cos(B\,t+\gamma_{\rm free})+ q_{\rm forced}(t),$ (10.103)

where

$\displaystyle h_{\rm forced}$ $\displaystyle = - \sum_{l=1,8} \frac{\nu_l}{A-g_l}\,\sin(g_l\,t+\beta_l),$ (10.104)
$\displaystyle k_{\rm forced}$ $\displaystyle = - \sum_{l=1,8} \frac{\nu_l}{A-g_l}\,\cos(g_l\,t+\beta_l),$ (10.105)
$\displaystyle p_{\rm forced}$ $\displaystyle = - \sum_{l=1,8} \frac{\mu_l}{B-f_l}\,\sin(f_l\,t+\gamma_l),$ (10.106)
$\displaystyle q_{\rm forced}$ $\displaystyle = - \sum_{l=1,8} \frac{\mu_l}{B-f_l}\,\cos(f_l\,t+\gamma_l),$ (10.107)

as well as

$\displaystyle \nu_l$ $\displaystyle = \sum_{i=1,8} A_i\,e_{il},$ (10.108)
$\displaystyle \mu_l$ $\displaystyle =\sum_{i=1,8} B_i\,I_{il}.$ (10.109)

The parameters, $e_{\rm free}$ and $I_{\rm free}$ appearing in Equations (10.100)–(10.103) are the eccentricity and inclination, respectively, that the asteroid orbit would possess were it not for the perturbing influence of the planets. These parameters are usually called the free, or proper, eccentricity and inclination, respectively. Roughly speaking, the planetary perturbations cause the osculating eccentricity, $e=(h^{\,2}+k^{\,2})^{1/2}$, and inclination, $I=\sin^{-1}([p^{\,2}+q^{\,2}]^{1/2})$, to oscillate about the corresponding free quantities, $e_{\rm free}$ and $I_{\rm free}$, respectively.

Figure: 10.7 The free eccentricity plotted against the sine of the free inclination (relative to J2000 equinox and ecliptic) for the orbits of the first 100,000 numbered asteroids at MJD 55400. The free orbital elements are determined from standard Laplace-Lagrange secular evolution theory. The most prominent Hirayama families are labeled. Raw data from JPL Small-Body Database.
\includegraphics[height=4in]{Chapter09/fig9_07.eps}

Figure 10.6 shows the osculating eccentricity plotted against the sine of the osculating inclination for the orbits of the first 100,000 numbered asteroids (asteroids are numbered in order of their discovery). No particular patten is apparent. Figures 10.7 and 10.8 show the free eccentricity plotted against the sine of the free inclination for the same 100,000 orbits. In Figure 10.7, the free orbital elements are determined from standard Laplace-Lagrange secular evolution theory, whereas in Figure 10.8 they are determined from Brouwer and van Woerkom's refinement of this theory. It can be seen that many of the points representing the asteroid orbits have condensed into clumps. These clumps, which are somewhat clearer in Figure 10.8 than in Figure 10.7, are known as Hirayama families after their discoverer, the Japanese astronomer Kiyotsugu Hirayama (1874-1943). It is thought that the asteroids making up a given family had a common origin; most likely due to the break up of some much larger body (Bertotti et al. 2003). As a consequence of this origin, the asteroids originally had similar orbital elements. However, as time progressed, these elements were jumbled by the perturbing influence of the planets. Thus, only when this influence is removed does the commonality of the orbits becomes apparent. Hirayama families are named after their largest member. The most prominent families are the (4) Vesta, (15) Eunomia, (24) Themis, (44) Nysa, (158) Koronis, (221) Eos, and (1272) Gefion families. (The number in brackets is that of the corresponding asteroid.)

Figure: 10.8 The free eccentricity plotted against the sine of the free inclination (relative to J2000 equinox and ecliptic) for the orbits of the first 100,000 numbered asteroids at MJD 55400. The free orbital elements are determined from Brouwer and van Woerkom's improved secular evolution theory. The most prominent Hirayama families are labeled. Raw data from JPL Small-Body Database.
\includegraphics[height=4in]{Chapter09/fig9_08.eps}