The effect of the Earth’s oblateness on the Moon’s physical libration in latitude

The Moon’s physical libration in latitude generated by gravitational forces caused by the Earth’s oblateness has been examined by a vector analytical method. Libration oscillations are described by a close set of five linear inhomogeneous differential equations, the dispersion equation has five roots, one of which is zero. A complete solution is obtained. It is revealed that the Earth’s oblateness: a) has little effect on the instantaneous axis of Moon’s rotation, but causes an oscillatory rotation of the body of the Moon with an amplitude of 0.072″ and pulsation period of 16.88 Julian years; b) causes small nutations of poles of the orbit and of the ecliptic along tight spirals, which occupy a disk with a cut in a center and with radius of 0.072″. Perturbations caused by the spherical Earth generate: a) physical librations in latitude with an amplitude of 34.275″; b) nutational motion for centers of small spiral nutations of orbit (ecliptic) pole over ellipses with semi-major axes of 113.850″ (85.158″) and the first pole rotates round the second one along a circle with radius of 28.691″; c) nutation of the Moon’s celestial pole over an ellipse with a semi-major axis of 45.04″ and with an axes ratio of about 0.004 with a period of T = 27.212 days. The principal ellipse’s axis is directed tangentially with respect to the precession circumference, along which the celestial pole moves nonuniformly nearly in one dimension. In contrast to the accepted concept, the latitude does not change while the Moon’s poles of rotation move. The dynamical reason for the inclination of the Moon’s mean equator with respect to the ecliptic is oblateness of the body of the Moon.     

Physical librations due to the third and fourth degree harmonics of the lunar gravity potential

The existence of third and fourth harmonics of the lunar gravity potential gives rise to sizable lunar physical librations. Using one recent set of potential estimates, the following effects are noted: the mean sub-Earth point is displaced from the earthward principal moment of inertia axis by 168″; the inclination of the lunar equator to the ecliptic is decreased by 14″.5; and a six year period libration in longitude, with amplitude 13″.1, is induced.     

Effect of the Earth’s compression on the physical libration of the Moon

The effect of the Earth??s compression on the physical libration of the Moon is studied using a new vector method. The moment of gravitational forces exerted on the Moon by the oblate Earth is derived considering second order harmonics. The terms in the expression for this moment are arranged according to their order of magnitude. The contribution due to a spherically symmetric Earth proves to be greater by a factor of 1.34 × 10⁶ than a typical term allowing for the oblateness. A linearized Euler system of equations to describe the Moon??s rotation with allowance for external gravitational forces is given. A full solution of the differential equation describing the Moon??s libration in longitude is derived. This solution includes both arbitrary and forced oscillation harmonics that we studied earlier (perturbations due to a spherically symmetric Earth and the Sun) and new harmonics due to the Earth??s compression. We posed and solved the problem of spinorbital motion considering the orientation of the Earth??s rotation axis with regard to the axes of inertia of the Moon when it is at a random point in its orbit. The rotation axes of the Earth and the Moon are shown to become coplanar with each other when the orbiting Moon has an ecliptic longitude of L _? = 90° or L _? = 270°. The famous Cassini??s laws describing the motion of the Moon are supplemented by the rule for coplanarity when proper rotations in the Earth-Moon system are taken into account. When we consider the effect of the Earth??s compression on the Moon??s libration in longitude, a harmonic with an amplitude of 0.03?? and period of T ₈ = 9.300 Julian years appears. This amplitude exceeds the most noticeable harmonic due to the Sun by a factor of nearly 2.7. The effect of the Earth??s compression on the variation in spin angular velocity of the Moon proves to be negligible.     

Orientation of the moon by numerical integration

The differential equations of rotational motion of the Moon are solved by numerical integration methods. Euler's dynamical equations transformed to a convenient form are treated by techniques analogous to ordinary orbit determination procedures. The proposed method is fully consistent with the ephemeris of the Moon and can utilize a variety of observational material for the solution of the selected parameters. The parameters are grouped into three distinct groups, namely:

--The physical libration angles of the Moon and their time rates at an arbitrary initial epoch.

--Physical constants featuring the principal moments of ineria of the Moon.

--Parameters associated with the particular observational material being used.

Examples are given of comparison between the proposed method and Eckhardt's 1970 model of the physical librations of the Moon. The merits of the new method are discussed in the light of conventional data sources like Earth-based or satellite-based photography as well as newly available data types like Laser ranging to retroreflectors on the Moon.     

Effects of a physical librations of the moon caused by a liquid core,and determination of the fourth mode of a free libration

An analytical theory of lunar physical librations based on its two-layer model consisting of a non-spherical solid mantle and ellipsoidal liquid core is developed. The Moon moves on a high-precision orbit in the gravitational field of the Earth and other celestial bodies. The defined fourth mode of a free libration is caused by the influence of the liquid core, with a long period of 205.7 yr, with amplitude S = 0″0395 and with an initial phase Π₀ = ?134° (for the initial epoch 2000.0). Estimates of dynamic (meridional) oblatenesses of a liquid core of the Moon have been estimated: ?_D = 4.42 × 10^?4, μ_D = 2.83 × 10^?4 (?_D + μ_D = 7.24 × 10^?4). These results have been obtained as a result of comparison of the developed analytical theory of physical librations of the Moon with the empirical theory of librations of the Moon constructed on the basis of laser observations.     

Planetary perturbations on the Libration of the Moon

A theory of the libration of the Moon, completely analytical with respect to the harmonic coefficients of the lunar gravity field, was recently built (Moons, 1982). The Lie transforms method was used to reduce the Hamiltonian of the main problem of the libration of the Moon and to produce the usual libration series p₁, p₂ and . This main problem takes into account the perturbations due to the Sun and the Earth on the rotation of a rigid Moon about its center of mass. In complement to this theory, we have now computed the planetary effects on the libration, the planetary terms being added to the mean Hamiltonian of the main problem before a last elimination of the angles. For the main problem, as well as for the planetary perturbations, the motion of the center of mass of the Moon is described by the ELP 2000 solution (Chapront and Chapront-Touze, 1983).     

A nonlinear analysis of the moon's physical libration in longitude

The Euler equations for the forced physical librations of the Moon have already been solved by using a digital computer to perform the semi-literal mathematical manipulations. Very near resonance, the computer solution for the physical libration in longitude is complemented by the solution of the appropriate Duffing equation with a dissipation term. Because of its apparent proximity to a resonant frequency, the term whose argument is 2 - twice the mean angular distance of the Moon's perigee from the ascending node of its orbit - is especially important. Its phase, which soon should be measurable, is related to the Moon's anelasticity. The term's frequency, in units of the sidereal month, increases as the semi-major axis of the Moon's orbit about the Earth increases. Using the Moon's mechanical ellipticity of Koziel and the rate of increase of the semi-major axis of MacDonald, it is estimated that the 2 term will cross the resonant frequency in 130 million years and, if the rate of energy dissipation is sufficiently low, a transient libration will be induced.     

Long-period forcing of Mercury's libration in longitude

Planetary perturbations of Mercury's orbit lead to forced librations in longitude in addition to the 88-day forced libration induced by Mercury's orbital motion. The forced librations are a combination of many periods, but 5.93 and 5.66 years dominate. These two periods result from the perturbations by Jupiter and Venus respectively, and they lead to a 125-year modulation of the libration amplitude corresponding to the beat frequency. Other periods are also identified with Jupiter and Venus perturbations as well as with those of the Earth, and these and other periods in the perturbations cause several arc second fluctuations in the libration extremes. The maxima of these extremes are about 30″ above and the minima about 7″ above the superposed ∼^60″, 88-day libration during the 125-year modulation. Knowledge of the nature of the long-period forced librations is important for the interpretation of the details of Mercury's rotation state to be obtained from radar and spacecraft observations. We show that the measurement of the 88-day libration amplitude for the purposes of determining Mercury's core properties is not compromised by the additional librations, because of the latter's small amplitude and long period. If the free libration in longitude has an amplitude that is large compared with that of the forced libration, its ∼10-year period will dominate the libration spectrum with the 88-day forced libration and the long-period librations from the orbital perturbations superposed. If the free libration has an amplitude that is comparable to those of the long-period forced libration, it will be revealed by erratic amplitude, period and phase on the likely time span of a series of observations. However, a significant free libration component is not expected because of relatively rapid damping.     

Lunar and terrestrial tidal effects on the Moon's rotational motion

An accurate model of the rotation of the Moon, constructed by numerical integration, has been presented in a previous paper. All direct perturbations capable of producing at least 10^–4 seconds of arc on the Moon's rotational motion have been included, and the physical librations resulting from planetary effects and Earth-Moon figure-figure interactions have been presented. The present study deals with the Moon's physical librations resulting from the non-rigidities of the Moon and the Earth. The effects of the Moon's elasticity and of a lunar phase lag are analyzed. Physical librations due to lunar tides and those due to terrestrial tides are presented and described.     

On the possibility of creating a radio-astronomical inertial coordinate system based on the measurement of arcs between QSOs with the aid of VLBI

The aim of this paper has been to discuss a possibility of the formation of an inertial system of coordinates, precise within 0″.001, with the aid of the measures of mutual angular distances of quasars by the VLBI technique used in radio-astronomy. Such a system depends on neither the rotation, nor revolution of the Earth around the Sun — a fact which permits a separation of the astronomical and geophysical aspects of the problem. The proposed system should permit us to resolve some general astrometric problems — such as of the precession and nutation of the Earth as well as of the motion of the terrestrial axis of rotation within the Earth, or fluctuations in the angular veoocity of terrestrial rotation.     

Lunar libration tables

The selenographic direction cosines of the Earth and of the pole of the ecliptic are developed in trigonometric series; the time dependence of each term enters through the Delaunay arguments. These and other series are tabulated for different libration parameters to provide means for calculating at any epoch the lunar librations and their partial derivatives with respect to the parameters.     

Sun–Earth L_{1} and L_{2} to Moon transfers exploiting natural dynamics

This paper examines the design of transfers from the Sun–Earth libration orbits, at the \(L_{1}\) and \(L_{2}\) points, towards the Moon using natural dynamics in order to assess the feasibility of future disposal or lifetime extension operations. With an eye to the probably small quantity of propellant left when its operational life has ended, the spacecraft leaves the libration point orbit on an unstable invariant manifold to bring itself closer to the Earth and Moon. The total trajectory is modeled in the coupled circular restricted three-body problem, and some preliminary study of the use of solar radiation pressure is also provided. The concept of survivability and event maps is introduced to obtain suitable conditions that can be targeted such that the spacecraft impacts, or is weakly captured by, the Moon. Weak capture at the Moon is studied by method of these maps. Some results for planar Lyapunov orbits at \(L_{1}\) and \(L_{2}\) are given, as well as some results for the operational orbit of SOHO.     

The problem of the Eulerian oscillations: A weakness of numerical versus analytical methods

The application of computer technology has permitted more and more problems of dynamical astronomy to be solved more easily, quickly, and accurately. In this area, numerical integration is often very efficient, and sometimes essential. There is often, however, a temptation to choose numerical integration simply because it is the easiest way to attack the problem. Sometimes this works to the detriment of a satisfactory understanding of the physics of the problem under study. It is particularly the case for the free, or Eulerian, oscillations. The forces that create such a motion are not of gravitational origin and are not even conservative. The theory can only specify the frequencies of oscillation, not their amplitudes nor phases. The case is complicated when the free oscillations interact with gravitationally-forced oscillations, a situation that is almost inevitable, since nothing is isolated in the Universe. The first author has particularly studied this problem in the case of the rotation of the Moon, and published the first credible determinations of the lunar free libration. In this kind of problem, the observations have to be used and care must to be taken to create no spurious free librations in the results by using numerical integrations to describe the other related motions. A differential correction of the starting conditions to fit the observations does not necessarily give any valid information on the real free oscillation contained in the data. An analytical model is necessary, if the goal of the research is to understand the origins and characteristics of an Eulerian oscillation in such a system.     

Selenographic control

Evaluation of selenographic data obtained with use of different observational means require the formulation of rigorous algorithms connecting the systems of coordinates, which the various methods have been referred to. The lunar principal axes of inertia are suggested as most appropriate for reference in lunar mapping and selenographic coordinate catalogues. The connection between the instantaneous axis of lunar rotation (involved in laser ranging, radar studies, astronomical observations from the surface of the Moon and VLBI observations of ALSEPs), the ecliptic system of coordinates (which in reductions of observations was considered as fixed in space), the Cassini mean selenographic coordinates (to which physical libration measures were referred), the lunar principal axes of inertia and the invariable plane of the solar system is discussed.On leave from the University of Manchester, England.Lunar Science Institute Contribution No. 138.Communication presented at the Conference on Lunar Dynamics and Observational Coordinate Systems, Held January 15–17, 1973, at the Lunar Science Institute, Houston, Tex., U.S.A.     

The vector approach to the problem of physical libration of the Moon: the linearized problem

A flexible and informative vector approach to the problem of physical libration of the rigid Moon has been developed in which three Euler differential equations are supplemented by 12 kinematic ones. A linearized system of equations can be split into an even and odd systems with respect to the reflection in the plane of the lunar equator, and rotational oscillations of the Moon are presented by superposition of librations in longitude and latitude. The former is described by three equations and consists of unrestricted oscillations with a period of T ₁ = 2.878 Julian years (amplitude of 1.855″) and forced oscillations with periods of T ₂ = 27.201 days (15.304″), one stellar year (0.008″), half a year (0.115″), and the third of a year (0.0003″) (five harmonics altogether). A zero frequency solution has also been obtained. The effect of the Sun on these oscillations is two orders of magnitude less than that of the Earth. The libration in latitude is presented by five equations and, at pertrubations from the Earth, is described by two harmonics of unrestricted oscillations (T ₅ ≈ 74.180 Julian years, T ₆ ≈ 27.347 days) and one harmonic of forced oscillations (T ₃ = 27.212 days). The motion of the true pole is presented by the same harmonics, with the maximum deviation from the Cassini pole being 45.3″. The fifth (zero) frequency yields a stationary solution with a conic precession of the rotation axis (previously unknown). The third Cassini law has been proved. The amplitudes of unrestricted oscillations have been determined from comparison with observations. For the ratio $ \frac{{\sin I}} {{\sin \left( {I + i} \right)}} \approx 0.2311 $ \frac{{\sin I}} {{\sin \left( {I + i} \right)}} \approx 0.2311 , the theory gives 0.2319, which confirms the adequacy of the approach. Some statements of the previous theory are revised. Poinsot’s method is shown to be irrelevant in describing librations of the Moon. The Moon does not have free (Euler) oscillations; it has oscillations with a period of T ₅ ≈ 74.180 Julian years rather than T ≈ 148.167 Julian years.     

The vector approach to the problem of physical libration of the Moon. II. The nonlinear problem

By the new vector method in a nonlinear setting, a physical libration of the Moon is studied. Using the decomposition method on small parameters we derive the closed system of nine differential equations with terms of the first and second order of smallness. The conclusion is drawn that in the nonlinear case a connection between the librations in a longitude and latitude, though feeble, nevertheless exists; therefore, the physical libration already is impossible to subdivide into independent from each other forms of oscillations, as usually can be done. In the linear approach, ten characteristic frequencies and two special invariants of the problem are found. It is proved that, taking into account nonlinear terms, the invariants are periodic functions of time. Therefore, the stationary solution with zero frequency, formally supposing in the linear theory a resonance, in the nonlinear approach gains only small (proportional to e) periodic oscillations. Near to zero frequency of a resonance there is no, and solution of the nonlinear equations of physical libration is stable. The given nonlinear solution slightly modifies the previously unknown conical precession of the Moon’s spin axis. The character of nonlinear solutions near the basic forcing frequency Ω₁, where in the linear approach there are beats, is carefully studied. The average method on fast variables is obtained by the linear system of differential equations with almost periodic coefficients, which describe the evolution of these coefficients in a nonlinear problem. From this follows that the nonlinear components only slightly modify the specified beats; the interior period T ≈ 16.53 days appears 411 times less than the exterior one T ≈ 18.61 Julian years. In particular, with such a period the angle between ecliptic plane and Moon orbit plane also varies. Resonances, on which other researches earlier insisted, are not discovered. As a whole, the nonlinear analysis essentially improves and supplements a linear picture of the physical libration.     

On the 3/1 planar and circular resonant problem

The simplest model of a resonant problem of second order is the planar and circular case. Simplification like this is very old and for 3/1 resonance, several authors have studied this problem with different purposes. In this work, we test this model for the available asteroids, by applying Hori's perturbation method. Explicit solutions of the intermediate orbit are obtained. In the plane of two constants of the problem, all types of motion are described. By testing the model, it is shown that, in general, one can confirm results of numerical integrations indicating libration for a few number of asteroids and circulation for most of them. However, agreement in numerical values for amplitude and period of librations seems to be not possible mainly if Jupiter's eccentricity is neglected. On the other hand, even though there might be some physical reasons determining that only asteroids with high eccentricity may librate, it is shown that, from mathematical point of view, libration may occur even in the case of small eccentricities provided that some relations are satisfied.     

Determining parameters of Moon’s orbital and rotational motion from LLR observations using GRAIL and IERS-recommended models

The aim of this work is to combine the model of orbital and rotational motion of the Moon developed for DE430 with up-to-date astronomical, geodynamical, and geo- and selenophysical models. The parameters of the orbit and physical libration are determined in this work from lunar laser ranging (LLR) observations made at different observatories in 1970–2013. Parameters of other models are taken from solutions that were obtained independently from LLR. A new implementation of the DE430 lunar model, including the liquid core equations, was done within the EPM ephemeris. The postfit residuals of LLR observations make evident that the terrestrial models and solutions recommended by the IERS Conventions are compatible with the lunar theory. That includes: EGM2008 gravitational potential with conventional corrections and variations from solid and ocean tides; displacement of stations due to solid and ocean loading tides; and precession-nutation model. Usage of these models in the solution for LLR observations has allowed us to reduce the number of parameters to be fit. The fixed model of tidal variations of the geopotential has resulted in a lesser value of Moon’s extra eccentricity rate, as compared to the original DE430 model with two fit parameters. A mixed model of lunar gravitational potential was used, with some coefficients determined from LLR observations, and other taken from the GL660b solution obtained from the GRAIL spacecraft mission. Solutions obtain accurate positions for the ranging stations and the five retroreflectors. Station motion is derived for sites with long data spans. Dissipation is detected at the lunar fluid core-solid mantle boundary demonstrating that a fluid core is present. Tidal dissipation is strong at both Earth and Moon. Consequently, the lunar semimajor axis is expanding by 38.20 mm/yr, the tidal acceleration in mean longitude is \(-25.90 {{}^{\prime \prime }}/\mathrm{cy}^2\), and the eccentricity is increasing by \(1.48\times 10^{-11}\) each year.     

Passing through resonance: The excitation and dissipation of the lunar free libration in longitude

The acceleration of the mean lunar longitude has a small effect on the periods of most terms in a Fourier expansion of the longitude. There are several planetary perturbation terms that have small amplitudes, but whose periods are close to the resonant period of the lunar libration in longitude. Some of these terms are moving toward resonance, some are moving away from resonance, and the periods of those terms that do not include the Delaunay variables in their arguments are not moving. Because of its acceleration of longitude, the Moon is receding from the Earth, so the magnitude of the restoring torque that the Earth exerts on the rotating Moon is gradually attenuating; thus resonance itself is moving, but at a much slower rate than the periods of the accelerating planetary perturbations. There are five planetary perturbation terms from the ELP-2000 Ephemeris (with amplitudes of 0.00001 or greater) that have passed through resonance in the past two million years. One of them is of special interest because it appears to be the excitation source of a supposed free libration in longitude that has been detected by the lunar laser ranging experiment. The amplitude of the term is only 0.00021 but it could be the source of the 1 amplitude free libration term if the viscoelastic properties of the Moon are similar to those of the Earth.     

Relation of a contracting earth to the apparent accelerations of the sun and moon

The tidal theory of the evolution of the lunar orbit has remained inconsistent with the observational values of the apparent secular accelerations of the Sun and Moon since it was first developed by Jeffreys in 1920. Allowance for a changing moment of inertia of the Earth enables the discrepancy to be completely removed if a decrease is occurring at a rate of just about the amount already required by the phase-change theory of the nature of the terrestrial core. The agreement of the resulting theory with the latest determinations of the lunar acceleration increases confidence in the phase-change hypothesis. On the other hand the theory renders it most unlikely that a changing constant of gravitation will prove necessary to account for the observations. On the present theory of itself the Moon would have been extremely close to the Earth only about 10⁹ yr ago which suggests that some additional process may at times have influenced the lunar orbit.