Evection in latitude

Next, let us consider terms in the solution of the lunar equations of motions that depend linearly on the lunar inclination, $I$.

According to Equation (11.155),

$\displaystyle c_4 = -\frac{3}{2}\,x_4\,z_1.$ (11.305)

It follows from Equations (11.209), (11.213), and (11.217) that

$\displaystyle z_4 = \frac{3}{8}\,z_1 + \left(\frac{91-12\,g}{64}\right)\,z_1\,m.$ (11.306)

According to Equation (11.156),

$\displaystyle c_5 = \frac{3}{2}\,x_4\,z_1.$ (11.307)

It follows from Equations (11.209), (11.214), and (11.217) that

$\displaystyle z_5 =\frac{3}{16}\,z_1.$ (11.308)

According to Equation (11.152),

$\displaystyle c_1 = -\frac{3}{2}\,x_4\,z_4.$ (11.309)

It follows from Equation (11.217), as well as the previous expression for $z_4$, that

$\displaystyle c_1 = \frac{9}{16}\,z_1.$ (11.310)

The preceding expression can be combined with Equation (11.208) to give

$\displaystyle \left({\mit\Omega}_1^{\,2} -1 -\frac{3}{2}\,m^{\,2} + \frac{9}{16}\,m^{\,3}\right)z_1=0.$ (11.311)

The non-trivial solution of this equation is such that

$\displaystyle {\mit\Delta}\equiv {\mit\Omega}_1^{\,2} -1 -\frac{3}{2}\,m^{\,2} + \frac{9}{16}\,m^{\,3}=0.$ (11.312)

It follows from Equation (11.210) that

$\displaystyle {\mit\Delta} = \left(\frac{-3+4\,g}{2}\right)m^{\,2} +\frac{9}{16}\,m^{\,3} + {\cal O}(m^{\,4}).$ (11.313)

Thus, setting ${\mit\Delta}=0$, we get

$\displaystyle g = \frac{3}{4}-\frac{9}{32}\,m+{\cal O}(m^{\,2}).$ (11.314)

The arbitrary parameter $z_1$ is chosen such that the parameter $d_1$, appearing in Equation (11.124), is the same as in the undisturbed motion. Thus, making use of Equation (11.184),

$\displaystyle z_1 = 1-\frac{1}{6}\,m^{\,2}.$ (11.315)

Hence, Equations (11.306) and (11.308) reduce to

$\displaystyle z_4$ $\displaystyle = \frac{3}{8}+\frac{41}{32}\,m,$ (11.316)
$\displaystyle z_5$ $\displaystyle = \frac{3}{16}.$ (11.317)

It follows from Equations (11.122)–(11.124), (11.184), (11.187), and (11.188), as well as the previous expressions for $z_1$, $z_4$, and $z_5$, that the net perturbation to the lunar orbit due to terms in the solution of the lunar equations of motion that depend linearly on $I$ is

$\displaystyle \delta R$ $\displaystyle =0$ (11.318)
$\displaystyle \delta \lambda$ $\displaystyle = 0,$ (11.319)
$\displaystyle \delta \beta$ $\displaystyle =I\,\sin F+\left(\frac{3}{8}\,m\,I +\frac{25}{32}\,m^{\,2}\,I\right)\sin(2\,D-F) + \frac{11}{16}\,m^{\,2}\,I\,\sin(2\,D+F).$ (11.320)

The previous expressions are accurate to ${\cal O}(m^{\,2}\,I)$.

The first term on the right-hand side of expression (11.320) is Keplerian in origin (i.e., it is independent of the perturbing influence of the Sun). The second term, which is known as evection in latitude, is due to the combined action of the Sun and the inclination of the lunar orbit to the ecliptic. Evection in latitude can be thought of as causing a slight increase in the inclination of the lunar orbit at the times of the first and last quarter moons (i.e., when $D=90^\circ$ and $D=270^\circ$), and a slight decrease at the times of the new moon and the full moon (i.e., when $D=0^\circ$ and $D=180^\circ$). (See Section 11.18, Exercise 4.) Evection in latitude generates a perturbation in the lunar ecliptic latitude that oscillates sinusoidally with a period of 32.3 days, and has an amplitude (calculated with $m=0.07480$ and $I= 0.09008$) of 602 arc seconds. As before, the oscillation period is in good agreement with observations, but the amplitude is somewhat inaccurate [it should be 624 arc seconds (Chapront-TouzĂ© and Chapront 1988)] due to the omission of higher order (in $m$ and $I$) contributions.