New Trends in the Development of the Lunar Physical Libration Theory

A review of the modern state of the lunar libration theory is presented. A significant progress in the lunar investigation is achieved due to the simultaneous processing of results of the satellite Doppler tracing and of the lunar laser ranging. The data evidencing existence of a small iron core in the Moon are discussed. In this connection, the further development of the theory of rotation of the Moon presents the study of internal structure and dynamics of a lunar body. A model of a two-layer Moon can have a very advanced application to explain some observed phenomena and to be as a first approach in the modelling of internal processes determining the lunar rotation.     

Analytical theory of the libration of the Moon

This paper presents a new theory of the libration of the Moon, completely analytical with respect to the harmonic coefficients of the lunar gravity field. This field is represented through its third degree harmonics for the torque due to the Earth (second degree for the torque due to the Sun).The orbital motion of the Moon is described by the ELP 2000 solution (Chapront-Touzé, 1980) of the main problem of lunar theory.the physical libration variables are obtained as Poisson series and comparisons with the results of Eckhardt (Eckhardt, 1981) and Migus (Migus, 1980) are presented.     

On the analytic lunar and solar perturbations of a near Earth satellite

The disturbing function of the Moon (Sun) is expanded as a sum of products of two harmonic functions, one depending on the position of the satellite and the other on the position of the Moon (Sun). The harmonic functions depending on the position of the perturbing body are developed into trigonometric series with the ecliptic elementsl, l′, F, D and Γ of the lunar theory which are nearly linear with respect to time. Perturbation of elements are in the form of trigonometric series with the ecliptic lunar elements and the equatorial elements ω and Ω of the satellite so that analytic integration is simple and the results accurate over a long period of time.     

On the implications of an olivine dominated upper mantle on the development of a moon of fission origin

In a series of previous papers, a petrological model for the Moon has been developed based on the assumption that the Moon is a globe of differentiated terrestrial mantle material which fissioned from the Earth. One of the major constraints which this model matches is the hypothesis that the lunar upper mantle is dominated by pyroxene. However, it has been recently shown that olivine is most probably the major constituent of the lunar upper mantle and that, at least that part of the Moon has a composition which is very similar to that of pyrolite - the proposed composition of the Earth's mantle. As a result of this new model constraint, the previously proposed differentiation scheme for a Moon of fission origin is reviewed and found to be inadequate, despite modification, for explaining the near pyrolite composition of the lunar upper mantle. As a result, a solidification sequence, which has been proposed to explain the rhythmic banding in terrestrial ultra mafic complexes, is investigated and found to be able to account for the high olivine content of the upper mantle, assuming a pyrolite composition for the Moon.     

Precession/nutation solution consistent with the general planetary theory

In the present paper the equations of the translatory motion of the major planets and the Moon and the Poisson equations of the Earth’s rotation in Euler parameters are reduced to the secular system describing the evolution of the planetary and lunar orbits (independent of the Earth’s rotation) and the evolution of the Earth’s rotation (depending on the planetary and lunar evolution). Hence, the theory of the Earth’s rotation is presented by means of the series in powers of the evolutionary variables with quasi-periodic coefficients.     

The great-circle pattern of large circular maria: product of an earth-moon encounter

The circular maria - Orientale, Imbrium, Serenitatis, Crisium, Smythii, and Tsiolkovsky -lie nearly on a lunar great circle. This pattern can be considered the result of a very close, non-capture encounter between Moon and Earth early in solar-system history. Of critical importance in analyzing the effects of such an encounter is the position of the weightlessness limit of the Earth-Moon System which is located at about 1.63R _e, measured from the center of Earth to center of Moon. Within this weightlessness limit, material can be pulled from the lunar surface and interior by Earth's gravity and either escape from the Moon or be redistributed onto the lunar surface. In the case of an encounter with a non-spinning Moon, backfalling materials would be distributed along a lunar great circle. However, if the Moon is rotating during the encounter, the backfall pattern will deviate from the great circle, the amount depending on the rate and direction of spin. Such a close encounter model may be related to the pattern of circular maria if materials departing from the source region are visualized as spheroids of molten lunar upper mantle basalt. These spheroids, then, would impact onto the lunar surface to form a pattern of lava lakes. Radiometric dates from mare rocks are consistent with this model of mare formation if the older mare rock dates are considered to date the encounter and younger dates are considered to date subsequent volcanic eruptions on a structurally weakened Moon.     

Insight-building models for lunar range and range rate

The analysis of range or Doppler data between sites on the Earth and Moon requires an accurate computation of the lunar orbit and detailed models of the orientation of the Earth and Moon. Models constructed to understand range and range rate can lack detail, but if they include the largest lunar orbit variations, tracking stations on a rotating Earth, and lunar sites on a synchronously rotating Moon, then they will display the largest effects for orbit elements, Earth orientation, tracking station locations, and lunar site coordinates. The range and range rate are expanded into periodic series. To understand accurate solutions, the largest periodic terms that are sensitive to various solution parameters indicate the sensitivity of data to solution parameters and the time spans needed for their determination. Conclusions include: cylindrical coordinates work well for sites on the rapidly rotating Earth, but Cartesian coordinates are more natural for the synchronously rotating Moon since the series for the three coordinate projections are distinct. For range and range rate data, daily, semimonthly, monthly, and longer periods are present. For Doppler data, the daily periods may be stronger and more useful than the long periods, particularly for terms associated with the terrestrial tracking station. Doppler data do not determine the lander coordinate toward the Earth well. Observational strategies for range and Doppler data are not identical. For all data types, one wishes a variety of hour angles, lunar declinations, times of month, and longer periods. A long span of high-quality range data can improve the lunar orbit, orientation of the Earth’s equator, and physical librations. The locations of new lunar sites or new tracking stations can be determined from shorter spans of data.     

Lunar cratering asymmetries with high lunar orbital obliquity and inclination of the Moon

Accurate estimation of cratering asymmetry on the Moon is crucial for understanding Moon evolution history. Early studies of cratering asymmetry have omitted the contributions of high lunar obliquity and inclination. Here, we include lunar obliquity and inclination as new controlling variables to derive the cratering rate spatial variation as a function of longitude and latitude. With examining the influence of lunar obliquity and inclination on the asteroids population encountered by the Moon, we then have derived general formulas of the cratering rate spatial variation based on the crater scaling law. Our formulas with addition of lunar obliquity and inclination can reproduce the lunar cratering rate asymmetry at the current Earth-Moon distance and predict the apex/ant-apex ratio and the pole/equator ratio of this lunar cratering rate to be 1.36 and 0.87, respectively. The apex/ant-apex ratio is decreasing as the obliquity and inclination increasing. Combining with the evolution of lunar obliquity and inclination, our model shows that the apex/ant-apex ratio does not monotonically decrease with Earth-Moon distance and hence the influences of obliquity and inclination are not negligible on evolution of apex/ant-apex ratio. This model is generalizable to other planets and moons, especially for different spin-orbit resonances.     

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.     

The improvement of the lunar ephemeris

At present the fundamental lunar ephemeris is based on Brown's theory of the motion of the Moon with improvements based on the bypassing of Brown's Tables, the removal of the great empirical term, the substitution of the relevant constants of the IAU system of astronomical constants and the retransformation of Brown's series in rectangular coordinates to spherical coordinates. Even so this ephemeris does not represent adequately the recent range and range-rate radio observations, and it will be inadequate for use in the analysis of laser observations of corner reflectors on the Moon. Numerical integrations for these purposes have already been made at the Jet Propulsion Laboratory, but improved theoretical developments are also required; new solutions of the main problem are in hand elsewhere. Work at H.M. Nautical Almanac Office is aimed at obtaining improved values of the constants of the lunar orbit by a rediscussion of occultation observations made since 1943 and at the redevelopment of the series for the planetary perturbations using more precise theories of the motion of the Sun and planets. The techniques and preliminary results of exploratory numerical integrations were briefly described.Presented at the Conference on Celestial Mechanics, Oberwolfach, Germany, 17–23 August, 1969.     

On constructing the general Earth’s rotation theory

In the present paper the equations of the orbital motion of the major planets and the Moon and the equations of the three–axial rigid Earth’s rotation in Euler parameters are reduced to the secular system describing the evolution of the planetary and lunar orbits (independent of the Earth’s rotation) and the evolution of the Earth’s rotation (depending on the planetary and lunar evolution). Hence, the theory of the Earth’s rotation can be presented by means of the series in powers of the evolutionary variables with quasi-periodic coefficients with respect to the planetary–lunar mean longitudes. This form of the Earth’s rotation problem is compatible with the general planetary theory involving the separation of the short–period and long–period variables and avoiding the appearance of the non–physical secular terms.     

Detectability of lunar X-rays

We have examined several theories that imply the generation of X-rays by the Moon. The X-ray fluxes to be expected at the top of the Earth's atmosphere are estimated and compared. For example, we find that an X-ray flux is to be expected when the Moon is full and K_p high, as a consequence of the configuration of the auroral electrons in space deduced from the long tail model of the magnetosphere. The X-ray photons are caused by energetic electrons in the tail that bombard the lunar surface. Alternatively, Gold has suggested that lunar X-rays are produced by the bombardment by solar-wind electrons; this results in a lunar phase dependence that is different from the long tail model. The background is discussed and we conclude that the lunar X-ray flux may be detectable. Experiments of this kind may provide useful tools for investigating the models.     

Effects of physical librations of the Moon on the orbital elements of a lunar satellite

Physical librations of the Moon are small cyclic perturbations with periods of one month and longer, and amplitudes of 100 arc seconds or less. These cause the selenographic axes fixed in the true Moon to have a different orientation than similar axes fixed in the mean Moon.Physical librations have two types of effects of present interest. If the orbital elements of a lunar satellite are referred to selenographic axes in the true Moon as it rotates and librates, then the librations cause changes in the orientation angles (node, inclination, and periapsis argument of the satellite) large enough that long-period perturbation theory cannot be used without compensation for such geometrical effects. As a second effect, the gravitational potential of the Moon is actually wobbled in inertial space, a condition not included in the potential expression used in perturbation theory.This paper gives data on the magnitude of the physical librations, the geometrical effects on the orbital elements and the equivalent changes in the coefficients in the potential. It is shown that geometrical effects can be accommodated either by using an inertial axes system or by compensating for the lunar librations and precession when the selenographic axes are used. Further, it is shown that physical effects are small and negligible for all but the most exacting endeavors.     

Comparison of low-energy lunar transfer trajectories to invariant manifolds

In this study, transfer trajectories from the Earth to the Moon that encounter the Moon at various flight path angles are examined, and lunar approach trajectories are compared to the invariant manifolds of selected unstable orbits in the circular restricted three-body problem. Previous work focused on lunar impact and landing trajectories encountering the Moon normal to the surface, and this research extends the problem with different flight path angles in three dimensions. The lunar landing geometry for a range of Jacobi constants is computed, and approaches to the Moon via invariant manifolds from unstable orbits are analyzed for different energy levels.     

Two-dimensional guiding-center model of the solar wind-moon interaction

This paper presents a continued study of the two-dimensional guiding-center model of the solar wind interaction with the Moon. The characteristics theory and the computational method are discussed. The magnetic permeability of plasma is (1 + /2)^–1 in the solar wind flow upstream of the Moon, and it changes to 1 in the void region of the lunar wake. The gradual change of the magnetic permeability in the penumbral region from the interplanetary condition to the void condition is explained as the source of field perturbations in the lunar wake. Perturbations of the magnetic field propagate as magnetoacoustic waves in a frame of reference moving with the plasma flow. Computer solutions were obtained to show that (i) the two principal perturbations of the magnetic field in the lunar wake (the umbral increase and the penumbral decrease) are confined to a region bounded by a Mach cone tangent to the lunar body, and (ii) the penumbral increases occur outside the lunar Mach cone. Computer solutions are also used to identify the source of field perturbations and to simulate the solar wind-moon interaction under varying interplanetary conditions.     

A theory for the interpretation of lunar surface magnetometer data

The solution to the problem of the motion of the Moon relative to spatial irregularities in the interplanetary magnetic field is found. The lunar electrical conductivity is modeled by a two-layer conductivity profile. For the interaction of the Moon with the corotating sector structure of the interplanetary magnetic field it is found that the magnetic field in the lunar shell is the superposition of an oscillatory uniform field, an oscillatory dipole field and anoscillatory field that is toroidal about the axis of the motional electric field. With various lunar conductivity models and the theory of this paper, lunar surface magnetometer data can be quantitatively interpreted to yield information on the conductivity and consequently the temperature of the lunar core.Presently visiting the Max-Planck-Institut für Physik und Astrophysik, München, Germany.     

A new theory for the formation of the maria and cayley type lunar regions

A new liquefaction theory for the origin of the flat marial and Cayley areas on the lunar surface is described. It is supposed that the flat terrain in these areas resulted from periods in the development of the Moon when these regions, although not liquid, had a sufficiently low viscosity for the surfaces to relax more or less completely to a level form. To account for this low viscosity a model is developed in which, within these regions and for relatively short periods in the early history of the Moon, preferentially high temperatures were maintained close to the lunar surface. The paper examines in some detail the possibility that these high temperatures may have resulted from instabilities in the lunar heat flow pattern caused by the presence of a surface layer of very low thermal conductivity produced by the debris of early meteorite impacts.A comparison is made between current models for the formation of the lunar surface and the theory here proposed: the advantages of the latter are enumerated and discussed.Normally at Queen Mary College, University of London, England.     

Independent selenocentric system of coordinates by large-scale star-calibrated lunar photography

The selenocentric coordinates of 263 craters and one mountain are given. The scale and orientation of the selenocentric coordinates system is determined by star-calibrated lunar photographs obtained with the long-focus horizontal telescope. The origin of the system is defined by the theory of the motion of the Moon. The physical libration is taken into account for the reduction of the coordinates system to the principal axes of the lunar inertia. The reference stars coordinates are reduced to the FK4 system. The position of the lunar mass centre is determined by LURE-2 theory, the rotation parameters are taken from the Migus tables (Migus, 1980).     

On the systems of selenographic coordinates,their determination and terminology

Attention is drawn to the absence in literature of the precise definitions of selenographic and celestial selenocentric coordinate systems. In certain cases inaccuracies in the formulation of the first Cassini law occur. This is due to the fact that the principal directions dealt with in the theory of lunar rotation are being constantly confused. A clear-cut definition of the principal coordinate systems concerned with the lunar rotation is given. It is indicated that there is no necessity in a special astronomical time service on the Moon. Since the future expeditions to the Moon will be able to keep terrestrial time, the problem of the hour angle is simply solved by the Formula (11).     

Accretion process of the moon

Recent geochemical and geophysical data suggest that the initial temperature of the Moon was strongly peaked toward the lunar surface. To explain such an initial temperature distribution, a simple model of accretion process of the Moon is presented. The model assumes that the Moon was formed from the accumulation of the solid particles or gases in the isolated, closed cloud. Two equations are derived to calculate the accretion rate and surface temperature of the accreting Moon. Numerical calculations are made for a wide range of the parameters particle concentration and particle velocity in the cloud. A limited set of the parameters gives the initial temperature profiles as required by geochemical and geophysical data. These models of the proto-moon cloud indicate that the lunar outershell, about 400 km thick, was partially or completely molten just after the accretion of the Moon and that the Moon should have been formed in a period shorter than 1000 yr. If the Moon formed at a position nearer to the Earth than its present one, the Moon might have been formed in a period of less than one year.On leave from Geophysical Institute, University of Tokyo.Contribution No. 2104, Division of Geological and Planetary Sciences, California Institute of Technology.