Lunar theory


Lunar theory attempts to account for the motions of the Moon. There are many small variations in the Moon's motion, and many attempts have been made to account for them. After centuries of being problematic, lunar motion is now modeled to a very high degree of accuracy.
Lunar theory includes:
Lunar theory has a history of over 2000 years of investigation. Its more modern developments have been used over the last three centuries for fundamental scientific and technological purposes, and are still being used in that way.

Applications

Applications of lunar theory have included the following:
The Moon has been observed for millennia. Over these ages, various levels of care and precision have been possible, according to the techniques of observation available at any time. There is a correspondingly long history of lunar theories: it stretches from the times of the Babylonian and Greek astronomers, down to modern lunar laser ranging.
Among notable astronomers and mathematicians down the ages, whose names are associated with lunar theories, are:
; Babylonian/Chaldean:
; Greek/Hellenistic:
; Arab:
; European, 16th to early 20th centuries:
; North American, 19th to early 20th centuries:
Other notable mathematicians and mathematical astronomers also made significant contributions.
The history can be considered to fall into three parts: from ancient times to Newton; the period of classical physics; and modern developments.

Ancient times to Newton

Babylon

Of Babylonian astronomy, practically nothing was known to historians of science before the 1880s. Surviving ancient writings of Pliny had made bare mention of three astronomical schools in Mesopotamia – at Babylon, Uruk, and 'Hipparenum'. But definite modern knowledge of any details only began when Joseph Epping deciphered cuneiform texts on clay tablets from a Babylonian archive: In these texts he identified an ephemeris of positions of the Moon. Since then, knowledge of the subject, still fragmentary, has had to be built up by painstaking analysis of deciphered texts, mainly in numerical form, on tablets from Babylon and Uruk.
To the Babylonian astronomer Kidinnu has been attributed the invention of what is now called “System B” for predicting the position of the moon, taking account that the moon continually changes its speed along its path relative to the background of fixed stars. This system involved calculating daily stepwise changes of lunar speed, up or down, with a minimum and a maximum approximately each month. The basis of these systems appears to have been arithmetical rather than geometrical, but they did approximately account for the main lunar inequality now known as the equation of the center.
The Babylonians kept very accurate records for hundreds of years of new moons and eclipses. Some time between the years 500 BC and 400 BC they identified and began to use the 19 year cyclic relation between lunar months and solar years now known as the Metonic cycle.
This helped them built up a numerical theory of the main irregularities in the Moon's motion, reaching remarkably good estimates for the periods of the three most prominent features of the Moon's motion:
The Babylonian estimate for the synodic month was adopted for the greater part of two millennia by Hipparchus, Ptolemy, and medieval writers.

Greece and Hellenistic Egypt

Thereafter, from Hipparchus and Ptolemy in the Bithynian and Ptolemaic epochs down to the time of Newton's work in the seventeenth century, lunar theories were composed mainly with the help of geometrical ideas, inspired more or less directly by long series of positional observations of the moon. Prominent in these geometrical lunar theories were combinations of circular motions – applications of the theory of epicycles.

Hipparchus

, whose works are mostly lost and known mainly from quotations by other authors, assumed that the Moon moved in a circle inclined at 5° to the ecliptic, rotating in a retrograde direction once in 18 years. The circle acted as a deferent, carrying an epicycle along which the Moon was assumed to move in a retrograde direction. The center of the epicycle moved at a rate corresponding to the mean change in Moon's longitude, while the period of the Moon around the epicycle was an anomalistic month. This epicycle approximately provided for what was later recognized as the elliptical inequality, the equation of the center, and its size approximated to an equation of the center of about 5° 1'. This figure is much smaller than the modern value: but it is close to the difference between the modern coefficients of the equation of the center and that of the evection: the difference is accounted for by the fact that the ancient measurements were taken at times of eclipses, and the effect of the [|evection] was at that time unknown and overlooked. For further information see also separate article Evection.

Ptolemy

's work the Almagest had wide and long-lasting acceptance and influence for over a millennium. He gave a geometrical lunar theory that improved on that of Hipparchus by providing for a second inequality of the Moon's motion, using a device that made the apparent apogee oscillate a little – prosneusis of the epicycle. This second inequality or second anomaly accounted rather approximately, not only for the equation of the center, but also for what became known as the evection. But this theory, applied to its logical conclusion, would make the distance of the Moon appear to vary by a factor of about 2, which is clearly not seen in reality. This defect of the Ptolemaic theory led to proposed replacements by Ibn al-Shatir in the 14th century and by Copernicus in the 16th century.

Ibn al-Shatir and Copernicus

Significant advances in lunar theory were made by the Arab astronomer, Ibn al-Shatir . Drawing on the observation that the distance to the Moon did not change as drastically as required by Ptolemy's lunar model, he produced a new lunar model that replaced Ptolemy's crank mechanism with a double epicycle model that reduced the computed range of distances of the Moon from the Earth. A similar lunar theory, developed some 150 years later by the Renaissance astronomer Nicolaus Copernicus, had the same advantage concerning the lunar distances.

Tycho Brahe, Johannes Kepler, and Jeremiah Horrocks

and Johannes Kepler refined the Ptolemaic lunar theory, but did not overcome its central defect of giving a poor account of the variations in the Moon's distance, apparent diameter and parallax. Their work added to the lunar theory three substantial further discoveries.
  1. The nodes and the inclination of the lunar orbital plane both appear to librate, with a monthly or semi-annual period.
  2. The lunar longitude has a twice-monthly [|Variation], by which the Moon moves faster than expected at new and full moon, and slower than expected at the quarters.
  3. There is also an annual effect, by which the lunar motion slows down a little in January and speeds up a little in July: the annual equation.
The refinements of Brahe and Kepler were recognized by their immediate successors as improvements, but their seventeenth-century successors tried numerous alternative geometrical configurations for the lunar motions to improve matters further. A notable success was achieved by Jeremiah Horrocks, who proposed a scheme involving an approximate 6 monthly libration in the position of the lunar apogee and also in the size of the elliptical eccentricity. This scheme had the great merit of giving a more realistic description of the changes in distance, diameter and parallax of the Moon.

Newton

A first gravitational period for lunar theory started with the work of Newton. He was the first to define the problem of the perturbed motion of the Moon in recognisably modern terms. His groundbreaking work is shown for example in the Principia in all versions including the first edition published in 1687.

Solar perturbation of lunar motion

Newton identified how to evaluate the perturbing effect on the relative motion of the Earth and Moon, arising from their gravity towards the Sun, in Book 1, Proposition 66, and in Book 3, Proposition 25. The starting-point for this approach is Corollary VI to the laws of motion. This shows that if the external accelerative forces from some massive body happens to act equally and in parallel on some different other bodies considered, then those bodies would be affected equally, and in that case their motions would continue as if there were no such external accelerative forces at all. It is only in the case that the external forces are different in size or in direction in their accelerative effects on the different bodies considered, that consequent effects are appreciable on the relative motions of the latter bodies.
Thus Newton concluded that it is only the difference between the Sun's accelerative attraction on the Moon and the Sun's attraction on the Earth that perturbs the motion of the Moon relative to the Earth.
Newton then in effect used vector decomposition of forces, to carry out this analysis. In Book 1, Proposition 66 and in Book 3, Proposition 25, he showed by a geometrical construction, starting from the total gravitational attraction of the Sun on the Earth, and of the Sun on the Moon, the difference that represents the perturbing effect on the motion of the Moon relative to the Earth. In summary, line LS in Newton's diagram as shown below represents the size and direction of the perturbing acceleration acting on the Moon in the Moon's current position P.
Shown here is Newton's diagram from the first Latin edition of the Principia. Here he introduced his analysis of perturbing accelerations on the Moon in the Sun-Earth-Moon system. Q represents the Sun, S the Earth, and P the Moon.
Parts of this diagram represent distances, other parts gravitational accelerations. In a dual significance, SQ represents the Earth-Sun distance, and then it also represents the size and direction of the Earth-Sun gravitational acceleration. Other distances in the diagram are then in proportion to distance SQ. Other attractions are in proportion to attraction SQ.
The Sun's attractions are SQ and LQ. The size of LQ is drawn so that the ratio of attractions LQ:SQ is the inverse square of the ratio of distances PQ:SQ. The Earth's attraction on the Moon acts along direction PS..
After showing solar attractions LQ on the Moon and SQ on the Earth, on the same scale, Newton then makes a vector decomposition of LQ into components LM and MQ. Then he identifies the perturbing acceleration on the Moon as the difference of this from SQ. SQ and MQ are parallel to each other, so SQ can be directly subtracted from MQ, leaving MS. The resulting difference, after subtracting SQ from LQ, is therefore the vector sum of LM and MS: these add up to a perturbing acceleration LS.
Later Newton identified another resolution of the perturbing acceleration LM+MS = LS, into orthogonal components: a transverse component parallel to LE, and a radial component, effectively ES.
Newton's diagrammatic scheme, since his time, has been re-presented in other and perhaps visually clearer ways. Shown here is a vector presentation indicating, for two different positions, P1 and P2, of the Moon in its orbit around the Earth, the respective vectors LS1 and LS2 for the perturbing acceleration due to the Sun. The Moon's position at P1 is fairly close to what it was at P in Newton's diagram; corresponding perturbation LS1 is like Newton's LS in size and direction. At another position P2, the Moon is farther away from the Sun than the Earth is, the Sun's attraction LQ2 on the Moon is weaker than the Sun's attraction SQ=SQ2 on the Earth, and then the resulting perturbation LS2 points obliquely away from the Sun.
Constructions like those in Newton's diagram can be repeated for many different positions of the Moon in its orbit. For each position, the result is a perturbation vector like LS1 or LS2 in the second diagram. Shown here is an often-presented form of the diagram that summarises sizes and directions of the perturbation vectors for many different positions of the Moon in its orbit. Each small arrow is a perturbation vector like LS, applicable to the Moon in the particular position around the orbit from which the arrow begins. The perturbations on the Moon when it is nearly in line along the Earth-Sun axis, i.e. near new or full moon, point outwards, away from the Earth. When the Moon-Earth line is 90° from the Earth-Sun axis they point inwards, towards the Earth, with a size that is only half the maximum size of the axial perturbations.
Newton also showed that the same pattern of perturbation applies, not only to the Moon, in its relation to the Earth as disturbed by the Sun, but also to other particles more generally in their relation to the solid Earth as disturbed by the Sun ; for example different portions of the tidal waters at the Earth's surface. The study of the common pattern of these perturbing accelerations grew out of Newton's initial study of the perturbations of the Moon, which he also applied to the forces moving tidal waters. Nowadays this common pattern itself has become often known as a tidal force whether it is being applied to the disturbances of the motions of the Moon, or of the Earth's tidal waters – or of the motions of any other object that suffers perturbations of analogous pattern.
After introducing his diagram 'to find the force of the Sun to perturb the Moon' in Book 3, Proposition 25, Newton developed a first approximation to the solar perturbing force, showing in further detail how its components vary as the Moon follows its monthly path around the Earth. He also took the first steps in investigating how the perturbing force shows its effects by producing irregularities in the lunar motions.
For a selected few of the lunar inequalities, Newton showed in some quantitative detail how they arise from the solar perturbing force.
Much of this lunar work of Newton's was done in the 1680s, and the extent and accuracy of his first steps in the gravitational analysis was limited by several factors, including his own choice to develop and present the work in what was, on the whole, a difficult geometrical way, and by the limited accuracy and uncertainty of many astronomical measurements in his time.

Classical gravitational period after Newton

The main aim of Newton's successors, from Leonhard Euler, Alexis Clairaut and Jean d'Alembert in the mid-eighteenth century, down to E.W. Brown in the late nineteenth and early twentieth century, was to account completely and much more precisely for the moon's motions on the basis of Newton's laws, i.e. the laws of motion and of universal gravitation by attractions inversely proportional to the squares of the distances between the attracting bodies. They also wished to put the inverse-square law of gravitation to the test, and for a time in the 1740s it was seriously doubted, on account of what was then thought to be a large discrepancy between the Newton-theoretical and the observed rates in the motion of the lunar apogee. However Clairaut showed shortly afterwards that at least the major cause of the discrepancy lay not in the lunar theory based on Newton's laws, but in excessive approximations that he and others had relied on to evaluate it.
Most of the improvements in theory after Newton were made in algebraic form: they involved voluminous and highly laborious amounts of infinitesimal calculus and trigonometry. It also remained necessary, for completing the theories of this period, to refer to observational measurements.

Results of the theories

The lunar theorists used many different mathematical approaches to analyse the gravitational problem. Not surprisingly, their results tended to converge. From the time of the earliest gravitational analysts among Newton's successors, Euler, Clairaut and d'Alembert, it was recognized that nearly all of the main lunar perturbations could be expressed in terms of just a few angular arguments and coefficients. These can be represented by:
From these basic parameters, just four basic differential angular arguments are enough to express, in their different combinations, nearly all of the most significant perturbations of the lunar motions. They are given here with their conventional symbols due to Delaunay; they are sometimes known as the Delaunay arguments:
This work culminated into Brown's lunar theory and Tables of the Motion of the Moon. These were used in the American Ephemeris and Nautical Almanac until 1968, and in a modified form until 1984.

Largest or named lunar inequalities

Several of the largest lunar perturbations in longitude have been named. In terms of the differential arguments, they can be expressed in the following way, with coefficients rounded to the nearest second of arc :

Equation of the center

The analysts of the mid-18th century expressed the perturbations of the Moon's position in longitude using about 25-30 trigonometrical terms. However, work in the nineteenth and twentieth century led to very different formulations of the theory so these terms are no longer current. The number of terms needed to express the Moon's position with the accuracy sought at the beginning of the twentieth century was over 1400; and the number of terms needed to emulate the accuracy of modern numerical integrations based on laser-ranging observations is in the tens of thousands: there is no limit to the increase in number of terms needed as requirements of accuracy increase.

Modern developments

Digital computers and lunar laser ranging

Since the Second World War and especially since the 1960s, lunar theory has been further developed in a somewhat different way. This has been stimulated in two ways: on the one hand, by the use of automatic digital computation, and on the other hand, by modern observational data-types, with greatly increased accuracy and precision.
Wallace John Eckert, a student of Ernest William Brown who worked at IBM, used the experimental digital computers developed there after the Second World War for computation of astronomical ephemerides. One of the projects was to put Brown's lunar theory into the machine and evaluate the expressions directly. Another project was something entirely new: a numerical integration of the equations of motion for the Sun and the four major planets. This became feasible only after electronic digital computers became available. Eventually this led to the Jet Propulsion Laboratory Development Ephemeris series.
In the meantime, Brown's theory was improved with better constants and the introduction of Ephemeris Time and the removal of some empirical corrections associated with this. This led to the Improved Lunar Ephemeris, which, with some minor successive improvements, was used in the astronomical almanacs from 1960 through 1983 and was used in lunar landing missions.
The most significant improvement of position observations of the Moon have been the Lunar Laser Ranging measurements, obtained using Earth-bound lasers and special retroreflectors placed on the surface of the Moon. The time-of-flight of a pulse of laser light to one of the retroreflectors and back gives a measure of the Moon's distance at that time. The first of five retroreflectors that are operational today was taken to the Moon in the Apollo 11 spacecraft in July 1969 and placed in a suitable position on the Moon's surface by Neil Armstrong.
Its precision is still being extended further by the Apache Point Observatory Lunar Laser-ranging Operation, established in 2005.

Numerical integrations, relativity, tides, librations

The lunar theory, as developed numerically to fine precision using these modern measures, is based on a larger range of considerations than the classical theories: It takes account not only of gravitational forces but also of many tidal and geophysical effects and a greatly extended theory of lunar libration. Like many other scientific fields this one has now developed so as to be based on the work of large teams and institutions. An institution notably taking one of the leading parts in these developments has been the Jet Propulsion Laboratory at California Institute of Technology; and names particularly associated with the transition, from the early 1970s onwards, from classical lunar theories and ephemerides towards the modern state of the science include those of J. Derral Mulholland and J.G. Williams, and for the linked development of solar system ephemerides E. Myles Standish.
Since the 1970s, the Jet Propulsion Laboratory has produced a series of numerically integrated Development Ephemerides, incorporating Lunar Ephemerides. Planetary and lunar ephemerides DE200/LE200 were used in the official Astronomical Almanac ephemerides for 1984–2002, and ephemerides DE405/LE405, of further improved accuracy and precision, have been in use as from the issue for 2003.

Analytical developments

In parallel with these developments, a new class of analytical lunar theory has also been developed in recent years, notably the Ephemeride Lunaire Parisienne by Jean Chapront and Michelle Chapront-Touzé from the Bureau des Longitudes. Using computer-assisted algebra, the analytical developments have been taken further than previously could be done by the classical analysts working manually. Also, some of these new analytical theories have been fitted to the numerical ephemerides previously developed at JPL as mentioned above. The main aims of these recent analytical theories, in contrast to the aims of the classical theories of past centuries, have not been to generate improved positional data for current dates; rather, their aims have included the study of further aspects of the motion, such as long-term properties, which may not so easily be apparent from the modern numerical theories themselves.