Tidal acceleration

Last updated
A picture of Earth and the Moon from Mars. The presence of the Moon (which has about 1/81 the mass of Earth), is slowing Earth's rotation and extending the day by about 2 milliseconds every 100 years. The Earth and the Moon photographed from Mars orbit.jpg
A picture of Earth and the Moon from Mars. The presence of the Moon (which has about 1/81 the mass of Earth), is slowing Earth's rotation and extending the day by about 2 milliseconds every 100 years.

Tidal acceleration is an effect of the tidal forces between an orbiting natural satellite (e.g. the Moon), and the primary planet that it orbits (e.g. Earth). The acceleration causes a gradual recession of a satellite in a prograde orbit away from the primary, and a corresponding slowdown of the primary's rotation. The process eventually leads to tidal locking, usually of the smaller first, and later the larger body. The Earth–Moon system is the best-studied case.


The similar process of tidal deceleration occurs for satellites that have an orbital period that is shorter than the primary's rotational period, or that orbit in a retrograde direction.

The naming is somewhat confusing, because the speed of the satellite relative to the body it orbits is decreased as a result of tidal acceleration, and increased as a result of tidal deceleration.

Earth–Moon system

Discovery history of the secular acceleration

Edmond Halley was the first to suggest, in 1695, [1] that the mean motion of the Moon was apparently getting faster, by comparison with ancient eclipse observations, but he gave no data. (It was not yet known in Halley's time that what is actually occurring includes a slowing-down of Earth's rate of rotation: see also Ephemeris time – History. When measured as a function of mean solar time rather than uniform time, the effect appears as a positive acceleration.) In 1749 Richard Dunthorne confirmed Halley's suspicion after re-examining ancient records, and produced the first quantitative estimate for the size of this apparent effect: [2] a centurial rate of +10″ (arcseconds) in lunar longitude, which is a surprisingly accurate result for its time, not differing greatly from values assessed later, e.g. in 1786 by de Lalande, [3] and to compare with values from about 10″ to nearly 13″ being derived about a century later. [4] [5]

Pierre-Simon Laplace produced in 1786 a theoretical analysis giving a basis on which the Moon's mean motion should accelerate in response to perturbational changes in the eccentricity of the orbit of Earth around the Sun. Laplace's initial computation accounted for the whole effect, thus seeming to tie up the theory neatly with both modern and ancient observations. [6]

However, in 1854, John Couch Adams caused the question to be re-opened by finding an error in Laplace's computations: it turned out that only about half of the Moon's apparent acceleration could be accounted for on Laplace's basis by the change in Earth's orbital eccentricity. [7] Adams' finding provoked a sharp astronomical controversy that lasted some years, but the correctness of his result, agreed upon by other mathematical astronomers including C. E. Delaunay, was eventually accepted. [8] The question depended on correct analysis of the lunar motions, and received a further complication with another discovery, around the same time, that another significant long-term perturbation that had been calculated for the Moon (supposedly due to the action of Venus) was also in error, was found on re-examination to be almost negligible, and practically had to disappear from the theory. A part of the answer was suggested independently in the 1860s by Delaunay and by William Ferrel: tidal retardation of Earth's rotation rate was lengthening the unit of time and causing a lunar acceleration that was only apparent. [9]

It took some time for the astronomical community to accept the reality and the scale of tidal effects. But eventually it became clear that three effects are involved, when measured in terms of mean solar time. Beside the effects of perturbational changes in Earth's orbital eccentricity, as found by Laplace and corrected by Adams, there are two tidal effects (a combination first suggested by Emmanuel Liais). First there is a real retardation of the Moon's angular rate of orbital motion, due to tidal exchange of angular momentum between Earth and Moon. This increases the Moon's angular momentum around Earth (and moves the Moon to a higher orbit with a lower orbital speed). Secondly, there is an apparent increase in the Moon's angular rate of orbital motion (when measured in terms of mean solar time). This arises from Earth's loss of angular momentum and the consequent increase in length of day. [10]

A diagram of the Earth-Moon system showing how the tidal bulge is pushed ahead by Earth's rotation. This offset bulge exerts a net torque on the Moon, boosting it while slowing Earth's rotation. Tidal braking.svg
A diagram of the Earth–Moon system showing how the tidal bulge is pushed ahead by Earth's rotation. This offset bulge exerts a net torque on the Moon, boosting it while slowing Earth's rotation.

Effects of Moon's gravity

Because the Moon's mass is a considerable fraction of that of Earth (about 1:81), the two bodies can be regarded as a double planet system, rather than as a planet with a satellite. The plane of the Moon's orbit around Earth lies close to the plane of Earth's orbit around the Sun (the ecliptic), rather than in the plane of the earth's rotation (the equator) as is usually the case with planetary satellites. The mass of the Moon is sufficiently large, and it is sufficiently close, to raise tides in the matter of Earth. Foremost among such matter, the water of the oceans bulges out towards and away from the Moon. The average near-side bulge peak is moments after passing the Moon overhead, and Earth rotates over and under these bulges (in all but the far south and north split up by the north-south land masses) in just over a day. However, this rotation sees set in motion the bulges ahead of the position directly under the Moon as every tidal bulge has an inertia meaning it does not recede instantly. At every point in time this means the far more operative bulge, the moon-facing bulge, is offset from the line through the centers of Earth and the Moon, namely ahead of the Moon's orbit. Because of this, the gravitational pull of this nearer (so gravitationally important) tidal bulge and the Moon is not quite parallel to the Earth–Moon line, i.e. exerting torque between Earth and the Moon. This torque boosts the Moon in its orbit and slows the rotation of Earth.

As a result of this process, the mean solar day, which has to be 86,400 equal seconds, is actually getting longer when measured in SI seconds with stable atomic clocks. (The SI second, when adopted, was already a little shorter than the current value of the second of mean solar time. [11] ) The small difference accumulates over time, which leads to an increasing difference between our clock time (Universal Time) on the one hand, and Atomic Time and Ephemeris Time on the other hand: see ΔT. This led to the introduction of the leap second in 1972 [12] to compensate for differences in the bases for time standardization.

In addition to the effect of the ocean tides, there is also a tidal acceleration due to flexing of Earth's crust, but this accounts for only about 4% of the total effect when expressed in terms of heat dissipation. [13]

If other effects were ignored, tidal acceleration would continue until the rotational period of Earth matched the orbital period of the Moon. At that time, the Moon would always be overhead of a single fixed place on Earth. Such a situation already exists in the PlutoCharon system. However, the slowdown of Earth's rotation is not occurring fast enough for the rotation to lengthen to a month before other effects make this irrelevant: about 1 to 1.5 billion years from now, the continual increase of the Sun's radiation will likely cause Earth's oceans to vaporize, [14] removing the bulk of the tidal friction and acceleration. Even without this, the slowdown to a month-long day would still not have been completed by 4.5 billion years from now when the Sun will probably evolve into a red giant and likely destroy both Earth and the Moon. [15] [16]

Tidal acceleration is one of the few examples in the dynamics of the Solar System of a so-called secular perturbation of an orbit, i.e. a perturbation that continuously increases with time and is not periodic. Up to a high order of approximation, mutual gravitational perturbations between major or minor planets only cause periodic variations in their orbits, that is, parameters oscillate between maximum and minimum values. The tidal effect gives rise to a quadratic term in the equations, which leads to unbounded growth. In the mathematical theories of the planetary orbits that form the basis of ephemerides, quadratic and higher order secular terms do occur, but these are mostly Taylor expansions of very long time periodic terms. The reason that tidal effects are different is that unlike distant gravitational perturbations, friction is an essential part of tidal acceleration, and leads to permanent loss of energy from the dynamic system in the form of heat. In other words, we do not have a Hamiltonian system here.[ citation needed ]

Angular momentum and energy

The gravitational torque between the Moon and the tidal bulge of Earth causes the Moon to be constantly promoted to a slightly higher orbit and Earth to be decelerated in its rotation. As in any physical process within an isolated system, total energy and angular momentum are conserved. Effectively, energy and angular momentum are transferred from the rotation of Earth to the orbital motion of the Moon (however, most of the energy lost by Earth (−3.321 TW)[ citation needed ] is converted to heat by frictional losses in the oceans and their interaction with the solid Earth, and only about 1/30th (+0.121 TW) is transferred to the Moon). The Moon moves farther away from Earth (+38.247±0.004 mm/y), so its potential energy, which is still negative (in Earth's gravity well), increases, i. e. becomes less negative. It stays in orbit, and from Kepler's 3rd law it follows that its angular velocity actually decreases, so the tidal action on the Moon actually causes an angular deceleration, i.e. a negative acceleration (−25.858±0.003"/century2) of its rotation around Earth. The actual speed of the Moon also decreases. Although its kinetic energy decreases, its potential energy increases by a larger amount, i. e. Ep = -2Ec (Virial Theorem).

The rotational angular momentum of Earth decreases and consequently the length of the day increases. The net tide raised on Earth by the Moon is dragged ahead of the Moon by Earth's much faster rotation. Tidal friction is required to drag and maintain the bulge ahead of the Moon, and it dissipates the excess energy of the exchange of rotational and orbital energy between Earth and the Moon as heat. If the friction and heat dissipation were not present, the Moon's gravitational force on the tidal bulge would rapidly (within two days) bring the tide back into synchronization with the Moon, and the Moon would no longer recede. Most of the dissipation occurs in a turbulent bottom boundary layer in shallow seas such as the European Shelf around the British Isles, the Patagonian Shelf off Argentina, and the Bering Sea. [17]

The dissipation of energy by tidal friction averages about 3.75 terawatts, of which 2.5 terawatts are from the principal M2 lunar component and the remainder from other components, both lunar and solar. [18]

An equilibrium tidal bulge does not really exist on Earth because the continents do not allow this mathematical solution to take place. Oceanic tides actually rotate around the ocean basins as vast gyres around several amphidromic points where no tide exists. The Moon pulls on each individual undulation as Earth rotates—some undulations are ahead of the Moon, others are behind it, whereas still others are on either side. The "bulges" that actually do exist for the Moon to pull on (and which pull on the Moon) are the net result of integrating the actual undulations over all the world's oceans. Earth's net (or equivalent) equilibrium tide has an amplitude of only 3.23 cm, which is totally swamped by oceanic tides that can exceed one metre.

Historical evidence

This mechanism has been working for 4.5 billion years, since oceans first formed on Earth, but less so at times when much or most of the water was ice. There is geological and paleontological evidence that Earth rotated faster and that the Moon was closer to Earth in the remote past. Tidal rhythmites are alternating layers of sand and silt laid down offshore from estuaries having great tidal flows. Daily, monthly and seasonal cycles can be found in the deposits. This geological record is consistent with these conditions 620 million years ago: the day was 21.9±0.4 hours, and there were 13.1±0.1 synodic months/year and 400±7 solar days/year. The average recession rate of the Moon between then and now has been 2.17±0.31 cm/year, which is about half the present rate. The present high rate may be due to near resonance between natural ocean frequencies and tidal frequencies. [19]

Analysis of layering in fossil mollusc shells from 70 million years ago, in the Late Cretaceous period, shows that there were 372 days a year, and thus that the day was about 23.5 hours long then. [20] [21]

Quantitative description of the Earth–Moon case

The motion of the Moon can be followed with an accuracy of a few centimeters by lunar laser ranging (LLR). Laser pulses are bounced off mirrors on the surface of the Moon, emplaced during the Apollo missions of 1969 to 1972 and by Lunokhod 2 in 1973. [22] [23] Measuring the return time of the pulse yields a very accurate measure of the distance. These measurements are fitted to the equations of motion. This yields numerical values for the Moon's secular deceleration, i.e. negative acceleration, in longitude and the rate of change of the semimajor axis of the Earth–Moon ellipse. From the period 1970–2012, the results are:

−25.82 ± 0.03 arcsecond/century2 in ecliptic longitude [24]
+38.08 ± 0.04 mm/yr in the mean Earth–Moon distance [24]

This is consistent with results from satellite laser ranging (SLR), a similar technique applied to artificial satellites orbiting Earth, which yields a model for the gravitational field of Earth, including that of the tides. The model accurately predicts the changes in the motion of the Moon.

Finally, ancient observations of solar eclipses give fairly accurate positions for the Moon at those moments. Studies of these observations give results consistent with the value quoted above. [25]

The other consequence of tidal acceleration is the deceleration of the rotation of Earth. The rotation of Earth is somewhat erratic on all time scales (from hours to centuries) due to various causes. [26] The small tidal effect cannot be observed in a short period, but the cumulative effect on Earth's rotation as measured with a stable clock (ephemeris time, atomic time) of a shortfall of even a few milliseconds every day becomes readily noticeable in a few centuries. Since some event in the remote past, more days and hours have passed (as measured in full rotations of Earth) (Universal Time) than would be measured by stable clocks calibrated to the present, longer length of the day (ephemeris time). This is known as ΔT. Recent values can be obtained from the International Earth Rotation and Reference Systems Service (IERS). [27] A table of the actual length of the day in the past few centuries is also available. [28]

From the observed change in the Moon's orbit, the corresponding change in the length of the day can be computed:

+2.3 ms/d/century or +84 s/cy2 or +63 ns/d2.

However, from historical records over the past 2700 years the following average value is found:

+1.70 ± 0.05 ms/d/century [29] [30] or +62 s/cy2 or +46.5 ns/d2. (i.e. an accelerating cause is responsible of -0.6 ms/d/cy)

By twice integrating over the time, the corresponding cumulative value is a parabola having a coefficient of T2 (time in centuries squared) of (1/2) 62 s/cy2 :

ΔT = (1/2) 62 s/cy2 T2 = +31 s/cy2 T2.

Opposing the tidal deceleration of Earth is a mechanism that is in fact accelerating the rotation. Earth is not a sphere, but rather an ellipsoid that is flattened at the poles. SLR has shown that this flattening is decreasing. The explanation is that during the ice age large masses of ice collected at the poles, and depressed the underlying rocks. The ice mass started disappearing over 10000 years ago, but Earth's crust is still not in hydrostatic equilibrium and is still rebounding (the relaxation time is estimated to be about 4000 years). As a consequence, the polar diameter of Earth increases, and the equatorial diameter decreases (Earth's volume must remain the same). This means that mass moves closer to the rotation axis of Earth, and that Earth's moment of inertia decreases. This process alone leads to an increase of the rotation rate (phenomenon of a spinning figure skater who spins ever faster as they retract their arms). From the observed change in the moment of inertia the acceleration of rotation can be computed: the average value over the historical period must have been about −0.6 ms/century. This largely explains the historical observations.

Other cases of tidal acceleration

Most natural satellites of the planets undergo tidal acceleration to some degree (usually small), except for the two classes of tidally decelerated bodies. In most cases, however, the effect is small enough that even after billions of years most satellites will not actually be lost. The effect is probably most pronounced for Mars's second moon Deimos, which may become an Earth-crossing asteroid after it leaks out of Mars's grip.[ citation needed ] The effect also arises between different components in a binary star. [31]

Tidal deceleration

In tidal acceleration (1), a satellite orbits in the same direction as (but slower than) its parent body's rotation. The nearer tidal bulge (red) attracts the satellite more than the farther bulge (blue), imparting a net positive force (dotted arrows showing forces resolved into their components) in the direction of orbit, lifting it into a higher orbit.
In tidal deceleration (2) with the rotation reversed, the net force opposes the direction of orbit, lowering it. Tidal acceleration principle.svg
In tidal acceleration (1), a satellite orbits in the same direction as (but slower than) its parent body's rotation. The nearer tidal bulge (red) attracts the satellite more than the farther bulge (blue), imparting a net positive force (dotted arrows showing forces resolved into their components) in the direction of orbit, lifting it into a higher orbit.
In tidal deceleration (2) with the rotation reversed, the net force opposes the direction of orbit, lowering it.

This comes in two varieties:

  1. Fast satellites: Some inner moons of the giant planets and Phobos orbit within the synchronous orbit radius so that their orbital period is shorter than their planet's rotation. In other words, they orbit their planet faster than the planet rotates. In this case the tidal bulges raised by the moon on their planet lag behind the moon, and act to decelerate it in its orbit. The net effect is a decay of that moon's orbit as it gradually spirals towards the planet. The planet's rotation also speeds up slightly in the process. In the distant future these moons will strike the planet or cross within their Roche limit and be tidally disrupted into fragments. However, all such moons in the Solar System are very small bodies and the tidal bulges raised by them on the planet are also small, so the effect is usually weak and the orbit decays slowly. The moons affected are: Some hypothesize that after the Sun becomes a red giant, its surface rotation will be much slower and it will cause tidal deceleration of any remaining planets. [32]
  2. Retrograde satellites: All retrograde satellites experience tidal deceleration to some degree because their orbital motion and their planet's rotation are in opposite directions, causing restoring forces from their tidal bulges. A difference to the previous "fast satellite" case here is that the planet's rotation is also slowed down rather than sped up (angular momentum is still conserved because in such a case the values for the planet's rotation and the moon's revolution have opposite signs). The only satellite in the Solar System for which this effect is non-negligible is Neptune's moon Triton. All the other retrograde satellites are on distant orbits and tidal forces between them and the planet are negligible.

Mercury and Venus are believed to have no satellites chiefly because any hypothetical satellite would have suffered deceleration long ago and crashed into the planets due to the very slow rotation speeds of both planets; in addition, Venus also has retrograde rotation.


Size of the tidal bulge

Neglecting axial tilt, the tidal force a satellite (such as the Moon) exerts on a planet (such as Earth) can be described by the variation of its gravitational force over the distance from it, when this force is considered as applied to a unit mass :

where G is the universal gravitational constant, m is the satellite mass and r is the distance between the satellite and the planet.

Thus the satellite creates a disturbing potential on the planet, whose difference between the planet center and the closest (or farthest) point to the satellite is:

where A is the planet radius.

The size of the tidal bulge created on the planet can be estimated as roughly the ratio between this disturbing potential and the planet surface gravity:

A more exact calculation gives: [33]

assuming we neglect a second order effect due to rigidity of the planet material.

For the Moon-Earth system (m = 7.3 x 1022 kg, M = 6×1024 kg, A = 6.4 × 106 m, r = 3.8 × 108), this gives 0.7 meters, close to the true value for ocean tides height (roughly one meter).

Note that two bulges are formed, one centered roughly around the point nearest to the satellite and the other centered roughly around the point farthest from it.


Due to planet rotation, the bulges lags somewhat behind (?, is ahead of) the planet-satellite axis, which creates an angle between the two. The size of this lag angle depends on inertia and (much more importantly) on dissipation forces (e.g. friction) exerted on the bulge.

The satellite applies different forces on the close and far bulges. The difference is roughly times the planet diameter, where we replace the unit mass in the calculation above with the approximate mass of each bulge, (where ρ is the mass density of the bulge):

where we took into consideration the effect of the lag angle .

In order to get a rough estimation for the torque exerted by the satellite on the planet, we need to multiply this difference with the lever length (which is the planet diameter), and with the sine of the lag angle, giving:

A more exact calculation adds a 2/5 factor due to the planet spherical form and gives: [33]

Inserting the value of H found above this is:

This can be written as:

Where k is a factor related that can be expressed by Love numbers, taking into considerations non-uniformity in the planet mass density; corrections due to planet rigidity, neglected above, also enter here. For Earth, most of the bulge is made of sea water and has no correction for rigidity, but its mass density is 0.18 the average Earth mass density (1 g/cm3 vs. 5.5 g/cm3), so . The literature uses a close value of 0.2 ( [34] )

A similar calculation can be done for the tides created on the planet by the Sun. Here, m should be replaced by the mass of the Sun, and r by the distance to the Sun. Since α depends on the dissipation properties of Earth, it is expected to be the same for both. The resulting torque is 20% that exerted by the Moon.

Relation of the lag angle to energy dissipation

The work exerted by the satellite over the planet is created by a force F acting along the path of movement of a mass units moving in velocity u in the planet (in fact, in the bulge).

Forces and locations depend on the relative angle to the planet-satellite axis θ, that changes periodically with the angular momentum Ω. Since the force in the planet spherical coordinate system is symmetrical in the direction towards the satellite and in the opposite direction (it is outwards in both), the dependence is approximated as sinusoidal in 2θ. Thus the force exerted on a unit mass is of the form:

and the translation projected on the same direction is of the form:

due to the lag angle. The velocity component in the direction of the force is therefore:

And so the total work exerted over a unit mass during one cycle is: [34]

In fact, almost all of this is dissipated (e.g. as friction), as explained below.

Looking now at the total energy from the satellite potential in one of the bulges, this is equal to the total work performed on this in a quarter of the total angular range, i.e. from zero to maximal displacement:

where we have defined , and approximated for small α in the last equality, thus neglecting it.

The fraction of energy dissipated in each cycle is represented by the effective specific dissipation function, denoted by and defined as the total dissipation in one cycle divided by . This gives: [34]

The value of this is estimated as 1/13 for Earth, where the bulge is mainly liquid, 10−1-10−2 for the other inner planets and the Moon, where the bulge is mainly solid, and as 10−3–10−5 for the outer, mostly gaseous planets. [33] [34]

With this value for Earth at hand, the torque can be calculated to be 4.4×1016 N m, only 13% above the measured value of 3.9×1016 N m. [34]

Note that in the distant past, the value of for the Earth–Moon system was probably somewhat smaller. [34]

Retardation of the planet's rotation

Again neglecting axial tilt, The change over time in the planet angular momentum L is equal to the torque. L in turn is the product of the angular velocity Ω with the moment of inertia I.

For a spherical planet of approximately uniform mass density, , where f is a factor depending on the planet structure; a spherical planet of uniform density has f = 2/5 = 0.4. Since the angular momentum This gives:

Since the Earth density is larger at depth, its moment of inertia is somewhat smaller, with f = 0.33. [35]

For the Earth-Moon system, taking of 1/13 and k = 0.2, we get the deceleration of the Earth's rotation dΩ/dt = -4.5×10−22 radian sec−2 = -924.37 " cy−2 which corresponds to the acceleration of the length of the day (LOD) of 61 s/cy2 or 1.7 ms/d/cy or 46 ns/d2. For a 24-hour day, this is equivalent to an increase of 17 seconds in 1 million years for the LOD, or 1 hour (i.e. lengthening of the day by 1 hour) in 210 million years. Due to the additional 20% effect of the Sun, the day lengthens by 1 hour in approximately 180 million years. This calculation is pure theory, assumes no dissipation nor storage of forces through frictional heat, which is unrealistic given air masses, oceans & tectonics. Objects in earth-moon system orbit similarly can drain inertia, for example: 2020 CD3

A similar calculation shows that the Earth had exerted angular momentum through tidal friction on the Moon self-rotation, before this became tidally locked. At that period, one calculates the change in the Moon angular momentum ω in the same manner as for Ω above, except that m and M should be is switched, and A should be replaced by the Moon radius a = 1.7×106 meter. Taking of 10−1 — 10−2 as for the solid planets and k = 1, this gives the deceleration of the Moon's rotation dω/dt = -3×10−17 — −3×10−18 radian sec−2. For a 29.5-day long rotation period, this is equivalent to 1.5 – 15 minutes in 1 year, or 1 day in 102 — 103 years. Thus in astronomical timescales, the Moon became tidally locked very fast.

Effect on the satellite motion around the planet

Due to conservation of angular momentum, a torque of the same size as the one exerted by the satellite and of opposite direction is exerted by the planet on the satellite motion around the planet. Another effect, which will not be dealt with here, is the changes in the eccentricity and inclination of the orbit.

The moment of inertia of this motion is mr2. However now r itself depends on the angular velocity which we denote here n: according to Newtonian analysis of orbital motion:

Thus the satellite orbit angular momentum, , satisfies (neglecting eccentricity):

Additionally, since , we have:

Note that assuming all rotations are on the same direction and Ω > ω, as time passes, the angular momentum of the planet decreases and hence that of the satellite orbit increases. Due to its relation with the planet-satellite distance, the latter increases, so the angular velocity of the satellite orbit decreases.

For the Earth-Moon system, dr/dt gives 1.212×10−9 meter per second (or nm/s), or 3.8247 cm per year (or also m/cy)[ 24]. This is a 1% increase in the Earth-Moon distance in 100 million years. The deceleration of the Moon dn/dt is -1.2588×10−23 radian sec−2 or -25.858 "/cy2, and for a period of 29.5 days (a synodic month) is equivalent to an increase of 38 ms/cy, or 7 minutes in 1 million years, or 1 day (i.e. lengthening of the lunar period in 1 day) in 210 million years.

Effect of the Sun

The Sun-planet system has two tidal friction effects. One effect is that the Sun creates a tidal friction in the planet, which decreases its spinning angular momentum and hence also increases its orbital angular momentum around the Sun, hence increasing its distance and reducing its angular velocity (assuming the orbital angular velocity of the Sun is smaller than that of the planet spinning; otherwise directions of change are opposite).

If MS is the Sun mass and D is the distance to it, then the rate of change of D is given, similar to the above calculation, by:

The planet orbital angular velocity, ΩS, then changes as:

For the Earth-Sun system, this gives 1×10−13 meters per second, or 3 meters in 1 million years. This is a 1% increase in the Earth-Sun distance in half a billion years. The deceleration of the Earth's orbital angular velocity is -2×10−31 radian sec2 or -410×10−9 "/cy2, or equivalently for a 1-year period, 1 second in 1 billion years.

Another, relatively negligible, effect is that the planet creates tidal friction in the Sun. This creates a change in the distance to the Sun and the orbital angular velocity around it, as it does for the satellite in the satellite-planet system. Using the same equations but now for the planet-Sun system, with AS standing for the Sun radius (7×108 meters), we have:

where kS is a factor, presumably very small, due to the non-uniformity of mass densities of the Sun. Assuming this factor times sin(2αS) to be not larger than what is found in the outer planets, i.e. 10−3 — 10−5, [33] we have a negligible contribution from this effect.

A detailed calculation for the Earth–Moon system

Potential perturbation created by the Moon on Earth

The potential per mass unit that the Moon creates on Earth, whose center is located at distance r0 from the Moon along the z-axis, in the Earth–Moon rotating frame of reference, and in coordinates centered at the Earth center, is:

where is the distance from the Moon to the center of mass of the Earth–Moon system, ω is the angular velocity of the Earth around this point (the same as the lunar orbital angular velocity). The second term is the effective potential due to the centrifugal force of the Earth.

We expand the potential in Taylor series around the point. The linear term must vanish (at least on average in time) since otherwise the force on the Earth center would be non vanishing. Thus:

Moving to spherical coordinates this gives:

where are the Legendre polynomials.

The constant term has no mechanical importance, while the causes a fixed dilation, and is not directly involved in creating a torque.

Thus we focus on the other terms, whose sum we denote , and mainly on the term which is the largest, as is at most the ratio of the Earth radius to its distance from the Moon, which is less than 2%.

Form of the bulge I: response to a perturbative potential

We treat the potential created by the Moon as a perturbation to the Earth's gravitational potential. Thus the height on Earth at angles , is:

where , and the amplitude of δ is proportional to the perturbation. We expand δ in Legendre polynomials, where the constant term (which stands for dilation) will be ignored as we are not interested in it. Thus:

where δn are unknown constants we would like to find.

We assume for the moment total equilibrium, as well as no rigidity on Earth (e.g. as in a liquid Earth). Therefore, its surface is equipotential, and so is constant, where is the Earth potential per unit mass. Since δ is proportional to , which is much smaller than VE, This can be expanded in δ. Dropping non-linear terms we have:

Note that is the force per unit mass from Earth's gravity, i.e. is just the gravitational acceleration g.

Since the Legendre polynomials are orthogonal, we may equate their coefficients n both sides of the equation, giving:

Thus the height is the ratio between the perturbation potential and the force from the perturbated potential.

Form of the bulge II: the deformation creating a perturbative potential

So far we have neglected the fact that the deformation itself creates a perturbative potential. In order to account for this, we may calculate this perturbative potential, re-calculate the deformation and continue so iteratively.

Let us assume the mass density is uniform. Since δ is much smaller than A, the deformation can be treated as a thin shell added to the mass of the Earth, where the shell has a surface mass density ρδ (and can also be negative), with ρ being the mass density (if mass density is not uniform, then the change of shape of the planet creates differences in mass distribution in all depth, and this has to be taken into account as well). Since the gravitation potential has the same form as the electric potential, this is a simple problem in electrostatics. For the analogous electrostatic problem, the potential created by the shell has the form:

where the surface charge density is proportional to the discontinuity in the gradient of the potential:

is the vacuum permittivity, a constant relevant to electrostatics, related to the equation . The analogous equation in gravity is , so if charge density is replaced with mass density, should be replaced with .

Thus in the gravitational problem we have:

So that, again due to the orthogonality of Legendre polynomials:

Thus the perturbative potential per mass unit for is:

Note that since Earth's mass density is in fact not uniform, this result must be multiplied by a factor that is roughly the ratio of the bulge mass density and the average Earth mass, approximately 0.18. The actual factor is somewhat larger, since there is some deformation in the deeper solid layers of Earth as well. Let us denote this factor by x. Rigidity also lowers x, though this is less relevant for most of the bulge, made of sea water.

The deformation was created by the a perturbative potential of size . Thus for each coefficient of , the ratio of the original perturbative potential to that secondarily created by the deformation is:

with x = 1 for perfectly a non-rigid uniform planet.

This secondary perturbative potential creates another deformation which again creates a perturbative potential and so on ad infinitum, so that the total deformation is of the size:

For each mode, the ratio to δn, the naive estimation of the deformation, is and is denoted as Love number . For a perfectly a non-rigid uniform planet (e.g. a liquid Earth of non-compressible liquid), this is equal to , and for the main mode of n = 2, it is 5/2.

Similarly, n-th mode of the tidal perturbative potential per unit mass created by Earth at r = A is the Love number kn times the corresponding term in the original lunar tidal perturbative potential, where for a uniform mass density, zero rigidity planet kn is:

For a perfectly a non-rigid uniform planet (e.g. a liquid Earth of non-compressible liquid), this is equal to 3/2. In fact, for the main mode of n  2, the real value for Earth is a fifth of it, namely k2 = 0.3 [34] (which fits c2 = 0.23 or x = 0.38, roughly twice the density ratios of 0.18).

Calculation of the torque

Instead of calculating the torque exerted by the Moon on the Earth deformation, we calculate the reciprocal torque exerted by the Earth deformation on the Moon; both must be equal.

The potential created by the Earth bulge is the perturbative potential we have discussed above. Per unit mass, for r = A, this is the same as the lunar perturbative potential creating the bulge, with each mode multiplied by kn, with the n = 2 mode far dominating the potential. Thus at r = A the bulge perturbative potential per unit mass is: [34]

since the n-the mode it drops off as r−(n+1) for r > A, we have outside Earth:

However, the bulge actually lags at an angle α with respect to the direction to the Moon due to Earth's rotation. Thus we have:

The Moon is at r = r0, θ = 0. Thus the potential per unit mass at the Moon is:

Neglecting eccentricity and axial tilt, We get the torque exerted by the bulge on the Moon by multiplying : with the Moon's mass m, and differentiating with respect to θ at the Moon location. This is equivalent to differentiating with respect to α, [34] and gives:

This is the same formula used above, with r = r0 and k there defined as 2k2/3.

See also

Related Research Articles

A centripetal force is a force that makes a body follow a curved path. Its direction is always orthogonal to the motion of the body and towards the fixed point of the instantaneous center of curvature of the path. Isaac Newton described it as "a force by which bodies are drawn or impelled, or in any way tend, towards a point as to a centre". In Newtonian mechanics, gravity provides the centripetal force causing astronomical orbits.

Keplers laws of planetary motion Scientific laws describing motion of planets around the Sun

In astronomy, Kepler's laws of planetary motion, published by Johannes Kepler between 1609 and 1619, describe the orbits of planets around the Sun. The laws modified the heliocentric theory of Nicolaus Copernicus, replacing its circular orbits and epicycles with elliptical trajectories, and explaining how planetary velocities vary. The three laws state that:

  1. The orbit of a planet is an ellipse with the Sun at one of the two foci.
  2. A line segment joining a planet and the Sun sweeps out equal areas during equal intervals of time.
  3. The square of a planet's orbital period is proportional to the cube of the length of the semi-major axis of its orbit.

Nutation is a rocking, swaying, or nodding motion in the axis of rotation of a largely axially symmetric object, such as a gyroscope, planet, or bullet in flight, or as an intended behaviour of a mechanism. In an appropriate reference frame it can be defined as a change in the second Euler angle. If it is not caused by forces external to the body, it is called free nutation or Euler nutation. A pure nutation is a movement of a rotational axis such that the first Euler angle is constant. In spacecraft dynamics, precession is sometimes referred to as nutation.

Orbit Gravitationally curved path of an object around a point in outer space

In physics, an orbit is the gravitationally curved trajectory of an object, such as the trajectory of a planet around a star or a natural satellite around a planet. Normally, orbit refers to a regularly repeating trajectory, although it may also refer to a non-repeating trajectory. To a close approximation, planets and satellites follow elliptic orbits, with the center of mass being orbited at a focal point of the ellipse, as described by Kepler's laws of planetary motion.

In mechanics and physics, simple harmonic motion is a special type of periodic motion where the restoring force on the moving object is directly proportional to the object's displacement magnitude and acts towards the object's equilibrium position. It results in an oscillation which, if uninhibited by friction or any other dissipation of energy, continues indefinitely.

Kinematics is a subfield of physics, developed in classical mechanics, that describes the motion of points, bodies (objects), and systems of bodies without considering the forces that cause them to move. Kinematics, as a field of study, is often referred to as the "geometry of motion" and is occasionally seen as a branch of mathematics. A kinematics problem begins by describing the geometry of the system and declaring the initial conditions of any known values of position, velocity and/or acceleration of points within the system. Then, using arguments from geometry, the position, velocity and acceleration of any unknown parts of the system can be determined. The study of how forces act on bodies falls within kinetics, not kinematics. For further details, see analytical dynamics.

Angular displacement

Angular displacement of a body is the angle in radians through which a point revolves around a centre or line has been rotated in a specified sense about a specified axis. When a body rotates about its axis, the motion cannot simply be analyzed as a particle, as in circular motion it undergoes a changing velocity and acceleration at any time (t). When dealing with the rotation of a body, it becomes simpler to consider the body itself rigid. A body is generally considered rigid when the separations between all the particles remains constant throughout the body's motion, so for example parts of its mass are not flying off. In a realistic sense, all things can be deformable, however this impact is minimal and negligible. Thus the rotation of a rigid body over a fixed axis is referred to as rotational motion.

In geometry, a solid angle is a measure of the amount of the field of view from some particular point that a given object covers. That is, it is a measure of how large the object appears to an observer looking from that point. The point from which the object is viewed is called the apex of the solid angle, and the object is said to subtend its solid angle from that point.

In mechanics and geometry, the 3D rotation group, often denoted SO(3), is the group of all rotations about the origin of three-dimensional Euclidean space under the operation of composition. By definition, a rotation about the origin is a transformation that preserves the origin, Euclidean distance, and orientation. Every non-trivial rotation is determined by its axis of rotation and its angle of rotation. Composing two rotations results in another rotation; every rotation has a unique inverse rotation; and the identity map satisfies the definition of a rotation. Owing to the above properties, the set of all rotations is a group under composition. Rotations are not commutative, making it a nonabelian group. Moreover, the rotation group has a natural structure as a manifold for which the group operations are smoothly differentiable; so it is in fact a Lie group. It is compact and has dimension 3.

Orbital mechanics Field of classical mechanics concerned with the motion of spacecraft

Orbital mechanics or astrodynamics is the application of ballistics and celestial mechanics to the practical problems concerning the motion of rockets and other spacecraft. The motion of these objects is usually calculated from Newton's laws of motion and law of universal gravitation. Orbital mechanics is a core discipline within space-mission design and control.

In physics, circular motion is a movement of an object along the circumference of a circle or rotation along a circular path. It can be uniform, with constant angular rate of rotation and constant speed, or non-uniform with a changing rate of rotation. The rotation around a fixed axis of a three-dimensional body involves circular motion of its parts. The equations of motion describe the movement of the center of mass of a body.

Rotating reference frame

A rotating frame of reference is a special case of a non-inertial reference frame that is rotating relative to an inertial reference frame. An everyday example of a rotating reference frame is the surface of the Earth.

The rigid rotor is a mechanical model of rotating systems. An arbitrary rigid rotor is a 3-dimensional rigid object, such as a top. To orient such an object in space requires three angles, known as Euler angles. A special rigid rotor is the linear rotor requiring only two angles to describe, for example of a diatomic molecule. More general molecules are 3-dimensional, such as water, ammonia, or methane. The rigid-rotor Schroedinger equation is discussed in Section 11.2 on pages 240-253 of the textbook by Bunker and Jensen.

Flight dynamics (spacecraft) Application of mechanical dynamics to model the flight of space vehicles

Spacecraft flight dynamics is the application of mechanical dynamics to model how the external forces acting on a space vehicle or spacecraft determine its flight path. These forces are primarily of three types: propulsive force provided by the vehicle's engines; gravitational force exerted by the Earth and other celestial bodies; and aerodynamic lift and drag.

Angular distance is the angle between the two sightlines, or between two point objects as viewed from an observer.

A pendulum is a body suspended from a fixed support so that it swings freely back and forth under the influence of gravity. When a pendulum is displaced sideways from its resting, equilibrium position, it is subject to a restoring force due to gravity that will accelerate it back toward the equilibrium position. When released, the restoring force acting on the pendulum's mass causes it to oscillate about the equilibrium position, swinging it back and forth. The mathematics of pendulums are in general quite complicated. Simplifying assumptions can be made, which in the case of a simple pendulum allow the equations of motion to be solved analytically for small-angle oscillations.

In geometry, various formalisms exist to express a rotation in three dimensions as a mathematical transformation. In physics, this concept is applied to classical mechanics where rotational kinematics is the science of quantitative description of a purely rotational motion. The orientation of an object at a given instant is described with the same tools, as it is defined as an imaginary rotation from a reference placement in space, rather than an actually observed rotation from a previous placement in space.

Geodetic effect

The geodetic effect represents the effect of the curvature of spacetime, predicted by general relativity, on a vector carried along with an orbiting body. For example, the vector could be the angular momentum of a gyroscope orbiting the Earth, as carried out by the Gravity Probe B experiment. The geodetic effect was first predicted by Willem de Sitter in 1916, who provided relativistic corrections to the Earth–Moon system's motion. De Sitter's work was extended in 1918 by Jan Schouten and in 1920 by Adriaan Fokker. It can also be applied to a particular secular precession of astronomical orbits, equivalent to the rotation of the Laplace–Runge–Lenz vector.

Newtons theorem of revolving orbits

In classical mechanics, Newton's theorem of revolving orbits identifies the type of central force needed to multiply the angular speed of a particle by a factor k without affecting its radial motion. Newton applied his theorem to understanding the overall rotation of orbits that is observed for the Moon and planets. The term "radial motion" signifies the motion towards or away from the center of force, whereas the angular motion is perpendicular to the radial motion.

Kepler orbit

In celestial mechanics, a Kepler orbit is the motion of one body relative to another, as an ellipse, parabola, or hyperbola, which forms a two-dimensional orbital plane in three-dimensional space. A Kepler orbit can also form a straight line. It considers only the point-like gravitational attraction of two bodies, neglecting perturbations due to gravitational interactions with other objects, atmospheric drag, solar radiation pressure, a non-spherical central body, and so on. It is thus said to be a solution of a special case of the two-body problem, known as the Kepler problem. As a theory in classical mechanics, it also does not take into account the effects of general relativity. Keplerian orbits can be parametrized into six orbital elements in various ways.


  1. E Halley (1695), "Some Account of the Ancient State of the City of Palmyra, with Short Remarks upon the Inscriptions Found there", Phil. Trans., vol.19 (1695–1697), pages 160–175; esp. at pages 174–175. (see also transcription using a modern font here)
  2. Richard Dunthorne (1749), "A Letter from the Rev. Mr. Richard Dunthorne to the Reverend Mr. Richard Mason F. R. S. and Keeper of the Wood-Wardian Museum at Cambridge, concerning the Acceleration of the Moon", Philosophical Transactions, Vol. 46 (1749–1750) #492, pp.162–172; also given in Philosophical Transactions (abridgements) (1809), vol.9 (for 1744–49), p669–675 as "On the Acceleration of the Moon, by the Rev. Richard Dunthorne".
  3. J de Lalande (1786): "Sur les equations seculaires du soleil et de la lune", Memoires de l'Academie Royale des Sciences, pp.390–397, at page 395.
  4. J D North (2008), "Cosmos: an illustrated history of astronomy and cosmology", (University of Chicago Press, 2008), chapter 14, at page 454.
  5. See also P Puiseux (1879), "Sur l'acceleration seculaire du mouvement de la Lune", Annales Scientifiques de l'Ecole Normale Superieure, 2nd series vol.8 (1879), pp.361–444, at pages 361–365.
  6. Britton, John (1992). Models and Precision: The Quality of Ptolemy's Observations and Parameters. Garland Publishing Inc. p.  157. ISBN   978-0815302155.
  7. Adams, J C (1853). "On the Secular Variation of the Moon's Mean Motion". Phil. Trans. R. Soc. Lond. 143: 397–406. doi: 10.1098/rstl.1853.0017 .
  8. D. E. Cartwright (2001), "Tides: a scientific history", (Cambridge University Press 2001), chapter 10, section: "Lunar acceleration, Earth retardation and tidal friction" at pages 144–146.
  9. Khalid, M.; Sultana, M.; Zaidi, F. (2014). "Delta: Polynomial Approximation of Time Period 1620–2013". Journal of Astrophysics. 2014: 1–4. doi: 10.1155/2014/480964 .
  10. F R Stephenson (2002), "Harold Jeffreys Lecture 2002: Historical eclipses and Earth's rotation", in Astronomy & Geophysics , vol.44 (2002), pp. 2.22–2.27.
  11. (1) In McCarthy, D D; Hackman, C; Nelson, R A (2008). "The Physical Basis of the Leap Second". Astronomical Journal. 136 (5): 1906–1908. Bibcode:2008AJ....136.1906M. doi: 10.1088/0004-6256/136/5/1906 . it is stated (page 1908), that "the SI second is equivalent to an older measure of the second of UT1, which was too small to start with and further, as the duration of the UT1 second increases, the discrepancy widens." :(2) In the late 1950s, the cesium standard was used to measure both the current mean length of the second of mean solar time (UT2) (result: 9192631830 cycles) and also the second of ephemeris time (ET) (result:9192631770±20 cycles), see "Time Scales", by L. Essen, in Metrologia, vol.4 (1968), pp.161–165, on p.162. As is well known, the 9192631770 figure was chosen for the SI second. L Essen in the same 1968 article (p.162) stated that this "seemed reasonable in view of the variations in UT2".
  12. "What's a Leap Second". Timeanddate.com.
  13. Munk (1997). "Once again: once again—tidal friction". Progress in Oceanography. 40 (1–4): 7–35. Bibcode:1997PrOce..40....7M. doi:10.1016/S0079-6611(97)00021-9.
  14. Puneet Kollipara (22 January 2014), "Earth Won't Die as Soon as Thought", Science.
  15. Murray, C.D.; Dermott, Stanley F. (1999). Solar System Dynamics. Cambridge University Press. p. 184. ISBN   978-0-521-57295-8.
  16. Dickinson, Terence (1993). From the Big Bang to Planet X. Camden East, Ontario: Camden House. pp. 79–81. ISBN   978-0-921820-71-0.
  17. Munk, Walter (1997). "Once again: once again—tidal friction". Progress in Oceanography. 40 (1–4): 7–35. Bibcode:1997PrOce..40....7M. doi:10.1016/S0079-6611(97)00021-9.
  18. Munk, W.; Wunsch, C (1998). "Abyssal recipes II: energetics of tidal and wind mixing". Deep-Sea Research Part I. 45 (12): 1977–2010. Bibcode:1998DSRI...45.1977M. doi:10.1016/S0967-0637(98)00070-3.
  19. Williams, George E. (2000). "Geological constraints on the Precambrian history of Earth's rotation and the Moon's orbit". Reviews of Geophysics. 38 (1): 37–60. Bibcode:2000RvGeo..38...37W. CiteSeerX . doi:10.1029/1999RG900016.
  20. "Ancient shell shows days were half-hour shorter 70 million years ago: Beer stein-shaped distant relative of modern clams captured snapshots of hot days in the late Cretaceous". ScienceDaily. Retrieved 2020-03-14.
  21. Winter, Niels J. de; Goderis, Steven; Malderen, Stijn J. M. Van; Sinnesael, Matthias; Vansteenberge, Stef; Snoeck, Christophe; Belza, Joke; Vanhaecke, Frank; Claeys, Philippe (2020). "Subdaily-Scale Chemical Variability in a Torreites Sanchezi Rudist Shell: Implications for Rudist Paleobiology and the Cretaceous Day-Night Cycle". Paleoceanography and Paleoclimatology. 35 (2): e2019PA003723. doi: 10.1029/2019PA003723 . ISSN   2572-4525.
  22. Most laser pulses, 78%, are to the Apollo 15 site. See Williams, et al. (2008), p. 5.
  23. Another reflector emplaced by Lunokhod 1 in 1970 is no longer functioning. See Lunar Lost & Found: The Search for Old Spacecraft by Leonard David
  24. 1 2 J.G. Williams, D.H. Boggs and W. M.Folkner (2013). DE430 Lunar Orbit, Physical Librations, and Surface Coordinates p.10. "These derived values depend on a theory which is not accurate to the number of digits given." See also : Chapront, Chapront-Touzé, Francou (2002). A new determination of lunar orbital parameters, precession constant and tidal acceleration from LLR measurements
  25. Stephenson, F.R.; Morrison, L.V. (1995). "Long-term fluctuations in the Earth's rotation: 700 BC to AD 1990" (PDF). Philosophical Transactions of the Royal Society of London Series A. 351 (1695): 165–202. Bibcode:1995RSPTA.351..165S. doi:10.1098/rsta.1995.0028. S2CID   120718607.
  26. Jean O. Dickey (1995): "Earth Rotation Variations from Hours to Centuries". In: I. Appenzeller (ed.): Highlights of Astronomy. Vol. 10 pp.17..44.
  27. "IERS – Observed values of UT1-TAI, 1962-1999". www.iers.org.
  28. "LOD". Archived from the original on September 8, 2001.
  29. Dickey, Jean O.; Bender, PL; Faller, JE; Newhall, XX; Ricklefs, RL; Ries, JG; Shelus, PJ; Veillet, C; et al. (1994). "Lunar Laser ranging: a continuing legacy of the Apollo program" (PDF). Science. 265 (5171): 482–90. Bibcode:1994Sci...265..482D. doi:10.1126/science.265.5171.482. PMID   17781305. S2CID   10157934.
  30. F. R. Stephenson (1997). Historical Eclipses and Earth's Rotation. Cambridge University Press. ISBN   978-0-521-46194-8.
  31. Zahn, J.-P. (1977). "Tidal Friction in Close Binary Stars". Astron. Astrophys. 57: 383–394. Bibcode:1977A&A....57..383Z.
  32. Schröder, K.-P.; Smith, R.C. (2008). "Distant future of the Sun and Earth revisited". Monthly Notices of the Royal Astronomical Society . 386 (1): 155–163. arXiv: 0801.4031 . Bibcode:2008MNRAS.386..155S. doi:10.1111/j.1365-2966.2008.13022.x. S2CID   10073988. See also Palmer, J. (2008). "Hope dims that Earth will survive Sun's death". New Scientist . Retrieved 2008-03-24.
  33. 1 2 3 4 Goldreich, Peter; Soter, Steven (1966). "Q in the solar system" (PDF). Icarus. 5 (1–6): 375–389. doi:10.1016/0019-1035(66)90051-0.
  34. 1 2 3 4 5 6 7 8 9 MacDonald, G. J. (1964). "Tidal friction". Reviews of Geophysics. 2 (3): 467–541. doi:10.1029/rg002i003p00467.
  35. Introduction to Earth Sciences I, Columbia university