Equation of time

Last updated

The equation of time: above the axis a sundial will appear fast relative to a clock showing local mean time, and below the axis a sundial will appear slow Equation of time.svg
The equation of time: above the axis a sundial will appear fast relative to a clock showing local mean time, and below the axis a sundial will appear slow
This graph shows how many minutes the clock is ahead (+) or behind (-) the apparent sun. See the section "Notation" below. Tijdvereffening-equation of time-en.jpg
This graph shows how many minutes the clock is ahead (+) or behind (−) the apparent sun. See the section "Notation" below.

The equation of time describes the discrepancy between two kinds of solar time. The word equation is used in the medieval sense of "reconciliation of a difference". The two times that differ are the apparent solar time, which directly tracks the diurnal motion of the Sun, and mean solar time, which tracks a theoretical mean Sun with uniform motion along the celestial equator. Apparent solar time can be obtained by measurement of the current position (hour angle) of the Sun, as indicated (with limited accuracy) by a sundial. Mean solar time, for the same place, would be the time indicated by a steady clock set so that over the year its differences from apparent solar time would have a mean of zero. [1]

Contents

The equation of time is the east or west component of the analemma, a curve representing the angular offset of the Sun from its mean position on the celestial sphere as viewed from Earth. The equation of time values for each day of the year, compiled by astronomical observatories, were widely listed in almanacs and ephemerides. [2] [3] :14

The equation of time can be approximated by a sum of two sine waves (see explanation below):

[minutes]

where represents the number of days since January 1 of the current year, .

Concept

Clock with auxiliary dial displaying the equation of time. Piazza Dante, Naples (1853) L' equazione del tempo. Foro Carolino in piazza Dante, Napoli.jpg
Clock with auxiliary dial displaying the equation of time. Piazza Dante, Naples (1853)

During a year the equation of time varies as shown on the graph; its change from one year to the next is slight. Apparent time, and the sundial, can be ahead (fast) by as much as 16  min  33  s (around 3 November), or behind (slow) by as much as 14 min 6 s (around 11 February). The equation of time has zeros near 15 April, 13 June, 1 September, and 25 December. Ignoring very slow changes in the Earth's orbit and rotation, these events are repeated at the same times every tropical year. However, due to the non-integral number of days in a year, these dates can vary by a day or so from year to year. As an example of the inexactness of the dates, according to the U.S. Naval Observatory's Multiyear Interactive Computer Almanac the equation of time was zero at 02:00 UT1 on 16 April 2011. [4] :277

The graph of the equation of time is closely approximated by the sum of two sine curves, one with a period of a year and one with a period of half a year. The curves reflect two astronomical effects, each causing a different non-uniformity in the apparent daily motion of the Sun relative to the stars:

The equation of time vanishes only for a planet with zero axial tilt and zero orbital eccentricity. [5] Two examples of planets with large equations of time are Mars and Uranus. On Mars the difference between sundial time and clock time can be as much as 50 minutes, due to the considerably greater eccentricity of its orbit. The planet Uranus, which has an extremely large axial tilt, has an equation of time that makes its days start and finish several hours earlier or later depending on where it is in its orbit.

Notation

Animation showing equation of time and analemma path over one year. EquationofTimeandAnalemma.gif
Animation showing equation of time and analemma path over one year.

The United States Naval Observatory states "the Equation of Time is the difference apparent solar time minus mean solar time", i.e. if the sun is ahead of the clock the sign is positive, and if the clock is ahead of the sun the sign is negative. [6] [7] The equation of time is shown in the upper graph above for a period of slightly more than a year. The lower graph (which covers exactly one calendar year) has the same absolute values but the sign is reversed as it shows how far the clock is ahead of the sun. Publications may use either format: in the English-speaking world, the former usage is the more common, but is not always followed. Anyone who makes use of a published table or graph should first check its sign usage. Often, there is a note or caption which explains it. Otherwise, the usage can be determined by knowing that, during the first three months of each year, the clock is ahead of the sundial. The mnemonic "NYSS" (pronounced "nice"), for "new year, sundial slow", can be useful. Some published tables avoid the ambiguity by not using signs, but by showing phrases such as "sundial fast" or "sundial slow" instead. [8]

History

The phrase "equation of time" is derived from the medieval Latin aequātiō diērum, meaning "equation of days" or "difference of days". The word aequātiō (and Middle English equation ) was used in medieval astronomy to tabulate the difference between an observed value and the expected value (as in the equation of the centre, the equation of the equinoxes, the equation of the epicycle). Gerald J. Toomer uses the medieval term "equation", from the Latin aequātiō (equalization or adjustment), for Ptolemy's difference between the mean solar time and the apparent solar time. Johannes Kepler's definition of the equation is "the difference between the number of degrees and minutes of the mean anomaly and the degrees and minutes of the corrected anomaly." [9] :155

The difference between apparent solar time and mean time was recognized by astronomers since antiquity, but prior to the invention of accurate mechanical clocks in the mid-17th century, sundials were the only reliable timepieces, and apparent solar time was the generally accepted standard. Mean time did not supplant apparent time in national almanacs and ephemerides until the early 19th century. [10]

Early astronomy

The irregular daily movement of the Sun was known to the Babylonians.[ citation needed ]

Book III of Ptolemy's Almagest (2nd century) is primarily concerned with the Sun's anomaly, and he tabulated the equation of time in his Handy Tables. [11] Ptolemy discusses the correction needed to convert the meridian crossing of the Sun to mean solar time and takes into consideration the nonuniform motion of the Sun along the ecliptic and the meridian correction for the Sun's ecliptic longitude. He states the maximum correction is 8+13 time-degrees or 59 of an hour (Book III, chapter 9). [12] However he did not consider the effect to be relevant for most calculations since it was negligible for the slow-moving luminaries and only applied it for the fastest-moving luminary, the Moon.

Based on Ptolemy's discussion in the Almagest , values for the equation of time (Arabic taʿdīl al-ayyām bi layālayhā) were standard for the tables (zij) in the works of medieval Islamic astronomy. [13]

Early modern period

A description of apparent and mean time was given by Nevil Maskelyne in the Nautical Almanac for 1767: "Apparent Time is that deduced immediately from the Sun, whether from the Observation of his passing the Meridian, or from his observed Rising or Setting. This Time is different from that shewn by Clocks and Watches well regulated at Land, which is called equated or mean Time." He went on to say that, at sea, the apparent time found from observation of the Sun must be corrected by the equation of time, if the observer requires the mean time. [1]

The right time was originally considered to be that which was shown by a sundial. When good mechanical clocks were introduced, they agreed with sundials only near four dates each year, so the equation of time was used to "correct" their readings to obtain sundial time. Some clocks, called equation clocks, included an internal mechanism to perform this "correction". Later, as clocks became the dominant good timepieces, uncorrected clock time, i.e., "mean time", became the accepted standard. The readings of sundials, when they were used, were then, and often still are, corrected with the equation of time, used in the reverse direction from previously, to obtain clock time. Many sundials, therefore, have tables or graphs of the equation of time engraved on them to allow the user to make this correction. [8] :123

The equation of time was used historically to set clocks. Between the invention of accurate clocks in 1656 and the advent of commercial time distribution services around 1900, there were several common land-based ways to set clocks. A sundial was read and corrected with the table or graph of the equation of time. If a transit instrument was available, the sun's transit across the meridian (the moment the sun appears to be due south or north of the observer) was noted; the clock was then set to noon and offset by the number of minutes given by the equation of time for that date. A third method did not use the equation of time; instead, it used stellar observations to give sidereal time, exploiting the relationship between sidereal time and mean solar time. [14] :57–58

The first tables to give the equation of time in an essentially correct way were published in 1665 by Christiaan Huygens. [15] Huygens, following the tradition of Ptolemy and medieval astronomers in general, set his values for the equation of time so as to make all values positive throughout the year. [15] This meant that any clock being set to mean time by Huygens's tables was consistently about 15 minutes slow compared to today's mean time.

Another set of tables was published in 1672–73 by John Flamsteed, who later became the first Astronomer Royal of the new Royal Greenwich Observatory. These appear to have been the first essentially correct tables that gave today's meaning of Mean Time (previously, as noted above, the sign of the equation was always positive and it was set at zero when the apparent time of sunrise was earliest relative to the clock time of sunrise). Flamsteed adopted the convention of tabulating and naming the correction in the sense that it was to be applied to the apparent time to give mean time. [16]

The equation of time, correctly based on the two major components of the Sun's irregularity of apparent motion, was not generally adopted until after Flamsteed's tables of 1672–73, published with the posthumous edition of the works of Jeremiah Horrocks. [17] :49

Robert Hooke (1635–1703), who mathematically analyzed the universal joint, was the first to note that the geometry and mathematical description of the (non-secular) equation of time and the universal joint were identical, and proposed the use of a universal joint in the construction of a "mechanical sundial". [18] :219

18th and early 19th centuries

The corrections in Flamsteed's tables of 1672–1673 and 1680 gave mean time computed essentially correctly and without need for further offset. But the numerical values in tables of the equation of time have somewhat changed since then, owing to three factors:

A sundial made in 1812 by Whitehurst & Son, with a circular scale showing the equation of time correction. This is now on display in Derby Museum and Art Gallery. Derby Sundial C 5810.JPG
A sundial made in 1812 by Whitehurst & Son, with a circular scale showing the equation of time correction. This is now on display in Derby Museum and Art Gallery.

From 1767 to 1833, the British Nautical Almanac and Astronomical Ephemeris tabulated the equation of time in the sense 'add or subtract (as directed) the number of minutes and seconds stated to or from the apparent time to obtain the mean time'. Times in the Almanac were in apparent solar time, because time aboard ship was most often determined by observing the Sun. This operation would be performed in the unusual case that the mean solar time of an observation was needed. In the issues since 1834, all times have been in mean solar time, because by then the time aboard ship was increasingly often determined by marine chronometers. The instructions were consequently to add or subtract (as directed) the number of minutes stated to or from the mean time to obtain the apparent time. So now addition corresponded to the equation being positive and subtraction corresponded to it being negative.

As the apparent daily movement of the Sun is one revolution per day, that is 360° every 24 hours, and the Sun itself appears as a disc of about 0.5° in the sky, simple sundials can be read to a maximum accuracy of about one minute. Since the equation of time has a range of about 33 minutes, the difference between sundial time and clock time cannot be ignored. In addition to the equation of time, one also has to apply corrections due to one's distance from the local time zone meridian and summer time, if any.

The tiny increase of the mean solar day due to the slowing down of the Earth's rotation, by about 2 ms per day per century, which currently accumulates up to about 1 second every year, is not taken into account in traditional definitions of the equation of time, as it is imperceptible at the accuracy level of sundials.

Major components

Eccentricity of the Earth's orbit

Equation of time (red solid line) and its two main components plotted separately, the part due to the obliquity of the ecliptic (mauve dashed line) and the part due to the Sun's varying apparent speed along the ecliptic due to eccentricity of the Earth's orbit (dark blue dash & dot line) Zeitgleichung.png
Equation of time (red solid line) and its two main components plotted separately, the part due to the obliquity of the ecliptic (mauve dashed line) and the part due to the Sun's varying apparent speed along the ecliptic due to eccentricity of the Earth's orbit (dark blue dash & dot line)

The Earth revolves around the Sun. As seen from Earth, the Sun appears to revolve once around the Earth through the background stars in one year. If the Earth orbited the Sun with a constant speed, in a circular orbit in a plane perpendicular to the Earth's axis, then the Sun would culminate every day at exactly the same time, and be a perfect time keeper (except for the very small effect of the slowing rotation of the Earth). But the orbit of the Earth is an ellipse not centered on the Sun, and its speed varies between 30.287 and 29.291 km/s, according to Kepler's laws of planetary motion, and its angular speed also varies, and thus the Sun appears to move faster (relative to the background stars) at perihelion (currently around 3 January) and slower at aphelion a half year later. [20] [21] [22]

At these extreme points, this effect varies the apparent solar day by 7.9 s/day from its mean. Consequently, the smaller daily differences on other days in speed are cumulative until these points, reflecting how the planet accelerates and decelerates compared to the mean.

As a result, the eccentricity of the Earth's orbit contributes a periodic variation which is (in the first-order approximation) a sine wave with:

This component of the EoT is represented by aforementioned factor a:

Obliquity of the ecliptic

Sun and planets at local apparent noon (Ecliptic in red, Sun and Mercury in yellow, Venus in white, Mars in red, Jupiter in yellow with red spot, Saturn in white with rings). Middaysun.gif
Sun and planets at local apparent noon (Ecliptic in red, Sun and Mercury in yellow, Venus in white, Mars in red, Jupiter in yellow with red spot, Saturn in white with rings).

Even if the Earth's orbit were circular, the perceived motion of the Sun along our celestial equator would still not be uniform. [5] This is a consequence of the tilt of the Earth's rotational axis with respect to the plane of its orbit, or equivalently, the tilt of the ecliptic (the path the Sun appears to take in the celestial sphere) with respect to the celestial equator. The projection of this motion onto our celestial equator, along which "clock time" is measured, is a maximum at the solstices, when the yearly movement of the Sun is parallel to the equator (causing amplification of perceived speed) and yields mainly a change in right ascension. It is a minimum at the equinoxes, when the Sun's apparent motion is more sloped and yields more change in declination, leaving less for the component in right ascension, which is the only component that affects the duration of the solar day. A practical illustration of obliquity is that the daily shift of the shadow cast by the Sun in a sundial even on the equator is smaller close to the solstices and greater close to the equinoxes. If this effect operated alone, then days would be up to 24 hours and 20.3 seconds long (measured solar noon to solar noon) near the solstices, and as much as 20.3 seconds shorter than 24 hours near the equinoxes. [20] [23] [22]

In the figure on the right, we can see the monthly variation of the apparent slope of the plane of the ecliptic at solar midday as seen from Earth. This variation is due to the apparent precession of the rotating Earth through the year, as seen from the Sun at solar midday.

In terms of the equation of time, the inclination of the ecliptic results in the contribution of a sine wave variation with:

This component of the EoT is represented by the aforementioned factor "b":

Secular effects

The two above mentioned factors have different wavelengths, amplitudes and phases, so their combined contribution is an irregular wave. At epoch 2000 these are the values (in minutes and seconds with UT dates):

PointValueDate
minimum−14 min 15 s11 February
zero0 min 0 s15 April
maximum+3 min 41 s14 May
zero0 min 0 s13 June
minimum−6 min 30 s26 July
zero0 min 0 s1 September
maximum+16 min 25 s3 November
zero0 min 0 s25 December

[ citation needed ]

E.T. = apparent − mean. Positive means: Sun runs fast and culminates earlier, or the sundial is ahead of mean time. A slight yearly variation occurs due to presence of leap years, resetting itself every 4 years. The exact shape of the equation of time curve and the associated analemma slowly change over the centuries, due to secular variations in both eccentricity and obliquity. At this moment both are slowly decreasing, but they increase and decrease over a timescale of hundreds of thousands of years. [24]

On shorter timescales (thousands of years) the shifts in the dates of equinox and perihelion will be more important. The former is caused by precession, and shifts the equinox backwards compared to the stars. But it can be ignored in the current discussion as our Gregorian calendar is constructed in such a way as to keep the vernal equinox date at 20 March (at least at sufficient accuracy for our aim here). The shift of the perihelion is forwards, about 1.7 days every century. In 1246 the perihelion occurred on 22 December, the day of the solstice, so the two contributing waves had common zero points and the equation of time curve was symmetrical: in Astronomical Algorithms Meeus gives February and November extrema of 15 m 39 s and May and July ones of 4 m 58 s. Before then the February minimum was larger than the November maximum, and the May maximum larger than the July minimum. In fact, in years before −1900 (1901 BCE) the May maximum was larger than the November maximum. In the year −2000 (2001 BCE) the May maximum was +12 minutes and a couple seconds while the November maximum was just less than 10 minutes. The secular change is evident when one compares a current graph of the equation of time (see below) with one from 2000 years ago, e.g., one constructed from the data of Ptolemy. [25]

Practical use

If the gnomon (the shadow-casting object) is not an edge but a point (e.g., a hole in a plate), the shadow (or spot of light) will trace out a curve during the course of a day. If the shadow is cast on a plane surface, this curve will be a conic section (usually a hyperbola), since the circle of the Sun's motion together with the gnomon point define a cone. At the spring and autumnal equinoxes, the cone degenerates into a plane and the hyperbola into a line. With a different hyperbola for each day, hour marks can be put on each hyperbola which include any necessary corrections. Unfortunately, each hyperbola corresponds to two different days, one in each half of the year, and these two days will require different corrections. A convenient compromise is to draw the line for the "mean time" and add a curve showing the exact position of the shadow points at noon during the course of the year. This curve will take the form of a figure eight and is known as an analemma. By comparing the analemma to the mean noon line, the amount of correction to be applied generally on that day can be determined.

The equation of time is used not only in connection with sundials and similar devices, but also for many applications of solar energy. Machines such as solar trackers and heliostats have to move in ways that are influenced by the equation of time.

Civil time is the local mean time for a meridian that often passes near the center of the time zone, and may possibly be further altered by daylight saving time. When the apparent solar time that corresponds to a given civil time is to be found, the difference in longitude between the site of interest and the time zone meridian, daylight saving time, and the equation of time must all be considered. [26]

Calculation

The equation of time is obtained from a published table, or a graph. For dates in the past such tables are produced from historical measurements, or by calculation; for future dates, of course, tables can only be calculated. In devices such as computer-controlled heliostats the computer is often programmed to calculate the equation of time. The calculation can be numerical or analytical. The former are based on numerical integration of the differential equations of motion, including all significant gravitational and relativistic effects. The results are accurate to better than 1 second and are the basis for modern almanac data. The latter are based on a solution that includes only the gravitational interaction between the Sun and Earth, simpler than but not as accurate as the former. Its accuracy can be improved by including small corrections.

The following discussion describes a reasonably accurate (agreeing with almanac data to within 3 seconds over a wide range of years) algorithm for the equation of time that is well known to astronomers. [27] :89 It also shows how to obtain a simple approximate formula (accurate to within 1 minute over a large time interval), that can be easily evaluated with a calculator and provides the simple explanation of the phenomenon that was used previously in this article.

Mathematical description

The precise definition of the equation of time is: [28] :1529

The quantities occurring in this equation are:

Here time and angle are quantities that are related by factors such as: 2π radians = 360° = 1 day = 24 hours. The difference, EOT, is measurable since GHA is an angle that can be measured and Universal Time, UT, is a scale for the measurement of time. The offset by π = 180° = 12 hours from UT is needed because UT is zero at mean midnight while GMHA = 0 at mean noon. Universal Time is discontinuous at mean midnight so another quantity day number N, an integer, is required in order to form the continuous quantity time t: t = N + UT/24 hr days. Both GHA and GMHA, like all physical angles, have a mathematical, but not a physical discontinuity at their respective (apparent and mean) noon. Despite the mathematical discontinuities of its components, EOT is defined as a continuous function by adding (or subtracting) 24 hours in the small time interval between the discontinuities in GHA and GMHA.

According to the definitions of the angles on the celestial sphere GHA = GAST − α (see hour angle)
where:

On substituting into the equation of time, it is

Like the formula for GHA above, one can write GMHA = GAST − αM, where the last term is the right ascension of the mean Sun. The equation is often written in these terms as [4] :275 [30] :45

where αM = GAST − UT + offset. In this formulation a measurement or calculation of EOT at a certain value of time depends on a measurement or calculation of α at that time. Both α and αM vary from 0 to 24 hours during the course of a year. The former has a discontinuity at a time that depends on the value of UT, while the latter has its at a slightly later time. As a consequence, when calculated this way EOT has two, artificial, discontinuities. They can both be removed by subtracting 24 hours from the value of EOT in the small time interval after the discontinuity in α and before the one in αM. The resulting EOT is a continuous function of time.

Another definition, denoted E to distinguish it from EOT, is

Here GMST = GAST − eqeq, is the Greenwich mean sidereal time (the angle between the mean vernal equinox and the mean Sun in the plane of the equator). Therefore, GMST is an approximation to GAST (and E is an approximation to EOT); eqeq is called the equation of the equinoxes and is due to the wobbling, or nutation of the Earth's axis of rotation about its precessional motion. Since the amplitude of the nutational motion is only about 1.2 s (18″ of longitude) the difference between EOT and E can be ignored unless one is interested in subsecond accuracy.

A third definition, denoted Δt to distinguish it from EOT and E, and now called the Equation of Ephemeris Time [28] :1532 (prior to the distinction that is now made between EOT, E, and Δt the latter was known as the equation of time) is

here Λ is the ecliptic longitude of the mean Sun (the angle from the mean vernal equinox to the mean Sun in the plane of the ecliptic).

The difference Λ − (GMST − UT + offset) is 1.3 s from 1960 to 2040. Therefore, over this restricted range of years Δt is an approximation to EOT whose error is in the range 0.1 to 2.5 s depending on the longitude correction in the equation of the equinoxes; for many purposes, for example correcting a sundial, this accuracy is more than good enough.

Right ascension calculation

The right ascension, and hence the equation of time, can be calculated from Newton's two-body theory of celestial motion, in which the bodies (Earth and Sun) describe elliptical orbits about their common mass center. Using this theory, the equation of time becomes:

where the new angles that appear are:

To complete the calculation three additional angles are required:

The celestial sphere and the Sun's elliptical orbit as seen by a geocentric observer looking normal to the ecliptic showing the 6 angles (M, lp, a, n, l, E) needed for the calculation of the equation of time. For the sake of clarity the drawings are not to scale. EquationOfTimeGeom.svg
The celestial sphere and the Sun's elliptical orbit as seen by a geocentric observer looking normal to the ecliptic showing the 6 angles (M, λp, α, ν, λ, E) needed for the calculation of the equation of time. For the sake of clarity the drawings are not to scale.

All these angles are shown in the figure on the right, which shows the celestial sphere and the Sun's elliptical orbit seen from the Earth (the same as the Earth's orbit seen from the Sun). In this figure ε is the obliquity, while e = 1 − (b/a)2 is the eccentricity of the ellipse.

Now given a value of 0 ≤ M ≤ 2π, one can calculate α(M) by means of the following well-known procedure: [27] :89

First, given M, calculate E from Kepler's equation: [31] :159

Although this equation cannot be solved exactly in closed form, values of E(M) can be obtained from infinite (power or trigonometric) series, graphical, or numerical methods. Alternatively, note that for e = 0, E = M, and by iteration: [32] :2

This approximation can be improved, for small e, by iterating again:

,

and continued iteration produces successively higher order terms of the power series expansion in e. For small values of e (much less than 1) two or three terms of the series give a good approximation for E; the smaller e, the better the approximation.

Next, knowing E, calculate the true anomaly ν from an elliptical orbit relation [31] :165

The correct branch of the multiple valued function arctan x to use is the one that makes ν a continuous function of E(M) starting from νE=0 = 0. Thus for 0 ≤ E < π use arctan x = arctan x, and for π < E ≤ 2π use arctan x = arctan x + π. At the specific value E = π for which the argument of tan is infinite, use ν = E. Here arctan x is the principal branch, |arctan x| < π/2; the function that is returned by calculators and computer applications. Alternatively, this function can be expressed in terms of its Taylor series in e, the first three terms of which are:

.

For small e this approximation (or even just the first two terms) is a good one. Combining the approximation for E(M) with this one for ν(E) produces:

.

The relation ν(M) is called the equation of the center; the expression written here is a second-order approximation in e. For the small value of e that characterises the Earth's orbit this gives a very good approximation for ν(M).

Next, knowing ν, calculate λ from its definition:

The value of λ varies non-linearly with M because the orbit is elliptical and not circular. From the approximation for ν:

.

Finally, knowing λ calculate α from a relation for the right triangle on the celestial sphere shown above [33] :22

Note that the quadrant of α is the same as that of λ, therefore reduce λ to the range 0 to 2π and write

,

where k is 0 if λ is in quadrant 1, it is 1 if λ is in quadrants 2 or 3 and it is 2 if λ is in quadrant 4. For the values at which tan is infinite, α = λ.

Although approximate values for α can be obtained from truncated Taylor series like those for ν, [34] :32 it is more efficacious to use the equation [35] :374

where y = tan2(ε/2). Note that for ε = y = 0, α = λ and iterating twice:

.

Final calculation

The equation of time is obtained by substituting the result of the right ascension calculation into an equation of time formula. Here Δt(M) = M + λpα[λ(M)] is used; in part because small corrections (of the order of 1 second), that would justify using E, are not included, and in part because the goal is to obtain a simple analytical expression. Using two-term approximations for λ(M) and α(λ) allows Δt to be written as an explicit expression of two terms, which is designated Δtey because it is a first order approximation in e and in y.

1) minutes

This equation was first derived by Milne, [35] :375 who wrote it in terms of λ = M + λp. The numerical values written here result from using the orbital parameter values, e = 0.016709, ε = 23.4393° = 0.409093 radians, and λp = 282.9381° = 4.938201 radians that correspond to the epoch 1 January 2000 at 12 noon UT1. When evaluating the numerical expression for Δtey as given above, a calculator must be in radian mode to obtain correct values because the value of 2λp − 2π in the argument of the second term is written there in radians. Higher order approximations can also be written, [36] :Eqs (45) and (46) but they necessarily have more terms. For example, the second order approximation in both e and y consists of five terms [28] :1535

2)

This approximation has the potential for high accuracy, however, in order to achieve it over a wide range of years, the parameters e, ε, and λp must be allowed to vary with time. [27] :86 [28] :1531,1535 This creates additional calculational complications. Other approximations have been proposed, for example, Δte [27] :86 [37] which uses the first order equation of the center but no other approximation to determine α, and Δte2 [38] which uses the second order equation of the center.

The time variable, M, can be written either in terms of n, the number of days past perihelion, or D, the number of days past a specific date and time (epoch):

3) days days
4)
Curves of Dt and Dtey along with symbols locating the daily values at noon (at 10-day intervals) obtained from the Multiyear Interactive Computer Almanac vs d (day) for the year 2000 EquationOfTime612.png
Curves of Δt and Δtey along with symbols locating the daily values at noon (at 10-day intervals) obtained from the Multiyear Interactive Computer Almanac vs d (day) for the year 2000
Derivative of -Dt. The axis on the right shows the length of the solar day. Length of solar day.png
Derivative of −Δt. The axis on the right shows the length of the solar day.

Here MD is the value of M at the chosen date and time. For the values given here, in radians, MD is that measured for the actual Sun at the epoch, 1 January 2000 at 12 noon UT1, and D is the number of days past that epoch. At periapsis M = 2π, so solving gives D = Dp = 2.508109. This puts the periapsis on 4 January 2000 at 00:11:41 while the actual periapsis is, according to results from the Multiyear Interactive Computer Almanac [39] (abbreviated as MICA), on 3 January 2000 at 05:17:30. This large discrepancy happens because the difference between the orbital radius at the two locations is only 1 part in a million; in other words, radius is a very weak function of time near periapsis. As a practical matter this means that one cannot get a highly accurate result for the equation of time by using n and adding the actual periapsis date for a given year. However, high accuracy can be achieved by using the formulation in terms of D.

When D > Dp, M is greater than 2π and one must subtract a multiple of 2π (that depends on the year) from it to bring it into the range 0 to 2π. Likewise for years prior to 2000 one must add multiples of 2π. For example, for the year 2010, D varies from 3653 on 1 January at noon to 4017 on 31 December at noon; the corresponding M values are 69.0789468 and 75.3404748 and are reduced to the range 0 to 2π by subtracting 10 and 11 times 2π respectively.

One can always write:

5) D = nY + d

where:

The resulting equation for years after 2000, written as a sum of two terms, given 1), 4) and 5), is:

6) [minutes]

In plain text format:

7) EoT =  -7.659sin(6.24004077 + 0.01720197(365*(y-2000) + d)) + 9.863sin( 2 (6.24004077 + 0.01720197 (365*(y-2000) + d)) + 3.5932 ) [minutes]

Term "a" represents the contribution of eccentricty, term "b" represents contribution of obliquity.

The result of the computations is usually given as either a set of tabular values, or a graph of the equation of time as a function of d. A comparison of plots of Δt, Δtey, and results from MICA all for the year 2000 is shown in the figure. The plot of Δtey is seen to be close to the results produced by MICA, the absolute error, Err = |Δtey − MICA2000|, is less than 1 minute throughout the year; its largest value is 43.2 seconds and occurs on day 276 (3 October). The plot of Δt is indistinguishable from the results of MICA, the largest absolute error between the two is 2.46 s on day 324 (20 November).

Continuity

For the choice of the appropriate branch of the arctan relation with respect to function continuity a modified version of the arctangent function is helpful. It brings in previous knowledge about the expected value by a parameter. The modified arctangent function is defined as:

.

It produces a value that is as close to η as possible. The function round rounds to the nearest integer.

Applying this yields:

.

The parameter M + λp arranges here to set Δt to the zero nearest value which is the desired one.

Secular change

The difference between the MICA and Δt results was checked every 5 years over the range from 1960 to 2040. In every instance the maximum absolute error was less than 3 s; the largest difference, 2.91 s, occurred on 22 May 1965 (day 141). However, in order to achieve this level of accuracy over this range of years it is necessary to account for the secular change in the orbital parameters with time. The equations that describe this variation are: [27] :86 [28] :1531,1535

According to these relations, in 100 years (D = 36525), λp increases by about 0.5% (1.7°), e decreases by about 0.25%, and ε decreases by about 0.05%.

As a result, the number of calculations required for any of the higher-order approximations of the equation of time requires a computer to complete them, if one wants to achieve their inherent accuracy over a wide range of time. In this event it is no more difficult to evaluate Δt using a computer than any of its approximations.

In all this note that Δtey as written above is easy to evaluate, even with a calculator, is accurate enough (better than 1 minute over the 80-year range) for correcting sundials, and has the nice physical explanation as the sum of two terms, one due to obliquity and the other to eccentricity that was used previously in the article. This is not true either for Δt considered as a function of M or for any of its higher-order approximations.

Alternative calculation

Another procedure for calculating the equation of time can be done as follows. [37] Angles are in degrees; the conventional order of operations applies.

n = 360°/365.24 days,

where n is the Earth's mean angular orbital velocity in degrees per day, a.k.a. "the mean daily motion".

where D is the date, counted in days starting at 1 on 1 January (i.e. the days part of the ordinal date in the year). 9 is the approximate number of days from the December solstice to 31 December. A is the angle the Earth would move on its orbit at its average speed from the December solstice to date D.

B is the angle the Earth moves from the solstice to date D, including a first-order correction for the Earth's orbital eccentricity, 0.0167 . The number 3 is the approximate number of days from 31 December to the current date of the Earth's perihelion. This expression for B can be simplified by combining constants to:

.

Here, C is the difference between the angle moved at mean speed, and at the angle at the corrected speed projected onto the equatorial plane, and divided by 180° to get the difference in "half-turns". The value 23.44° is the tilt of the Earth's axis ("obliquity"). The subtraction gives the conventional sign to the equation of time. For any given value of x, arctan x (sometimes written as tan−1x) has multiple values, differing from each other by integer numbers of half turns. The value generated by a calculator or computer may not be the appropriate one for this calculation. This may cause C to be wrong by an integer number of half-turns. The excess half-turns are removed in the next step of the calculation to give the equation of time:

 minutes

The expression nint(C) means the nearest integer to C. On a computer, it can be programmed, for example, as INT(C + 0.5). Its value is 0, 1, or 2 at different times of the year. Subtracting it leaves a small positive or negative fractional number of half turns, which is multiplied by 720, the number of minutes (12 hours) that the Earth takes to rotate one half turn relative to the Sun, to get the equation of time.

Compared with published values, [8] this calculation has a root mean square error of only 3.7 s. The greatest error is 6.0 s. This is much more accurate than the approximation described above, but not as accurate as the elaborate calculation.

Solar declination

The value of B in the above calculation is an accurate value for the Sun's ecliptic longitude (shifted by 90°), so the solar declination δ becomes readily available:

which is accurate to within a fraction of a degree.

See also

Notes

  1. 1 2 Maskelyne, Nevil (1767). The Nautical Almanac and Astronomical Ephemeris. London: Commissioners of Longitude.
  2. Milham, Willis I. (1945). Time and Timekeepers. New York: Macmillan. pp. 11–15. ISBN   978-0780800083.
  3. British Commission on Longitude (1794). Nautical Almanac and Astronomical Ephemeris for the year 1803. London, UK: C. Bucton.
  4. 1 2 Heilbron, J. L. (1999). The Sun in the Church: Cathedrals as Solar Observatories. Cambridge, MA: Harvard University Press. ISBN   9780674005365.
  5. 1 2 Jenkins, Alejandro (2013). "The Sun's position in the sky". European Journal of Physics. 34 (3): 633–652. arXiv: 1208.1043 . Bibcode:2013EJPh...34..633J. doi:10.1088/0143-0807/34/3/633. S2CID   119282288.
  6. Astronomical Applications Department. "The Equation of Time". United States Naval Observatory. United States Navy. Retrieved 1 August 2022.
  7. Astronomical Applications Department. "Astronomical Almanac Glossary". United States Naval Observatory. United States Navy. Retrieved 1 August 2022.
  8. 1 2 3 Waugh, Albert E. (1973). Sundials, Their Theory and Construction . New York: Dover Publications. p.  205. ISBN   978-0-486-22947-8.
  9. Kepler, Johannes (1995). Epitome of Copernican Astronomy & Harmonies of the World. Prometheus Books. ISBN   978-1-57392-036-0.
  10. McCarthy & Seidelmann 2009, p. 9.
  11. Neugebauer, Otto (1975), A History of Ancient Mathematical Astronomy, Studies in the History of Mathematics and Physical Sciences, vol. 1, New York / Heidelberg / Berlin: Springer-Verlag, pp. 984–986, doi:10.1007/978-3-642-61910-6, ISBN   978-0-387-06995-1
  12. Toomer, G.J. (1998). Ptolemy's Almagest. Princeton University Press. p. 171. ISBN   978-0-691-00260-6.
  13. Kennedy, E. S. (1956). "A Survey of Islamic Astronomical Tables". Transactions of the American Philosophical Society. 46 (2): 141. doi:10.2307/1005726. hdl: 2027/mdp.39076006359272 . JSTOR   1005726.
    Reprinted in: Kennedy, E. S. (1989). A survey of Islamic astronomical tables (2nd ed.). Philadelphia, PA: American Philosophical Society. p. 19. ISBN   9780871694621.
  14. Olmstead, Dennison (1866). A Compendium of Astronomy. New York: Collins & Brother.
  15. 1 2 Huygens, Christiaan (1665). Kort Onderwys aengaende het gebruyck der Horologien tot het vinden der Lenghten van Oost en West. The Hague: [publisher unknown].
  16. Flamsteed, John (1673) [1672 for the imprint, and bound with other sections printed 1673]. De Inaequalitate Dierum Solarium. London: William Godbid.
  17. Vince, S. "A Complete System of Astronomy". 2nd edition, volume 1, 1814
  18. Mills, Allan (2007). "Robert Hooke's 'universal joint' and its application to sundials and the sundial-clock". Notes Rec. R. Soc. 61 (2). Royal Society Publishing: 219–236. doi: 10.1098/rsnr.2006.0172 .
  19. Maskelyne, Nevil (1764). "Some Remarks upon the Equation of Time, and the True Manner of Computing It". Philosophical Transactions. 54. Royal Society: 336–347. JSTOR   105569 .
  20. 1 2 "The Equation of Time". Royal Museums Greenwich . Archived from the original on 10 September 2015. Retrieved 29 January 2021.
  21. "Eccentricity". Analemma. Retrieved 29 January 2021.
  22. 1 2 Taylor, Kieran (4 November 2018). "The Equation of Time: Why Sundial time Differs From Clock Time Depending On Time Of Year". moonkmft. Retrieved 29 January 2021.
  23. "Obliquity". Analemma. Retrieved 29 January 2021.
  24. Karney, Kevin (December 2005). "Variation in the Equation of Time" (PDF).
  25. Meeus 1997.
  26. "How to find the exact time of solar noon, wherever you are in the world". Spot-On Sundials. London. Retrieved 23 July 2013.
  27. 1 2 3 4 5 Duffett-Smith, Peter; Zwart, Jonathan (2017). Practical astronomy with your calculator or spreadsheet (4th ed.). Cambridge, UK: Cambridge University Press. ISBN   9781108436076.
  28. 1 2 3 4 5 6 Hughes, David W.; Yallop, B. D.; Hohenkerk, C. Y. (June 1989). "The Equation of Time". Monthly Notices of the Royal Astronomical Society . 238 (4): 1529–1535. doi: 10.1093/mnras/238.4.1529 . ISSN   0035-8711.
  29. 1 2 Astronomical Applications Department. "Computing Approximate Sidereal Time". United States Naval Observatory. United States Navy. Retrieved 1 August 2022.
  30. Roy, A. E. (1988). Orbital motion (3rd ed.). Bristol, England: A. Hilger. ISBN   0-85274-228-2.
  31. 1 2 Moulton, Forest Ray (1914). An introduction to celestial mechanics (2nd ed.). Macmillan.
  32. Hinch, E. J. (2002). Perturbation methods. Cambridge, UK: Cambridge University Press. doi:10.1017/CBO9781139172189. ISBN   9781139172189.
  33. Burington, Richard S. (1965) [1933]. Handbook of Mathematical Tables and Formulas (4th ed.). McGraw-Hill. LCCN   63-23531.
  34. Whitman, Alan M (2007). "A Simple Expression for the Equation of Time". Compendium. 14. North American Sundial Society: 29–33. CiteSeerX   10.1.1.558.1314 . ISSN   1074-8059.
  35. 1 2 Milne, R. M. (December 1921). "Note on the Equation of Time". The Mathematical Gazette. 10 (155). Mathematical Association: 372–375. doi:10.1017/S0025557200232944. S2CID   126276982.
  36. Müller, M. (1995). "Equation of Time: Problem in Astronomy" (PDF). Acta Physica Polonica A. 88 (Supplement). Institute of Physics, Polish Academy of Sciences: S49–S67.
  37. 1 2 Williams, David O. (2009). "The Latitude and Longitude of the Sun". Archived from the original on 23 March 2012.
  38. "Approximate Solar Coordinates Archived 24 July 2011 at the Wayback Machine ", "Naval Oceanography Portal".
  39. United States Naval Observatory April 2010, Multiyear Interactive Computer Almanac (version 2.2.1), Richmond VA: Willmann-Bell.

Related Research Articles

<span class="mw-page-title-main">Gravitational redshift</span> Shift of wavelength of a photon to longer wavelength

In physics and general relativity, gravitational redshift is the phenomenon that electromagnetic waves or photons travelling out of a gravitational well lose energy. This loss of energy corresponds to a decrease in the wave frequency and increase in the wavelength, known more generally as a redshift. The opposite effect, in which photons gain energy when travelling into a gravitational well, is known as a gravitational blueshift. The effect was first described by Einstein in 1907, eight years before his publication of the full theory of relativity.

<span class="mw-page-title-main">Kepler's laws of planetary motion</span> Laws describing the motion of planets

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.
<span class="mw-page-title-main">Astronomical coordinate systems</span> System for specifying positions of celestial objects

In astronomy, coordinate systems are used for specifying positions of celestial objects relative to a given reference frame, based on physical reference points available to a situated observer. Coordinate systems in astronomy can specify an object's relative position in three-dimensional space or plot merely by its direction on a celestial sphere, if the object's distance is unknown or trivial.

<span class="mw-page-title-main">Great-circle distance</span> Shortest distance between two points on the surface of a sphere

The great-circle distance, orthodromic distance, or spherical distance is the distance between two points on a sphere, measured along the great-circle arc between them. This arc is the shortest path between the two points on the surface of the sphere.

In mathematics and its applications, a Sturm–Liouville problem is a second-order linear ordinary differential equation of the form for given functions , and , together with some boundary conditions at extreme values of . The goals of a given Sturm–Liouville problem are:

<span class="mw-page-title-main">Orthographic map projection</span> Azimuthal perspective map projection

Orthographic projection in cartography has been used since antiquity. Like the stereographic projection and gnomonic projection, orthographic projection is a perspective projection in which the sphere is projected onto a tangent plane or secant plane. The point of perspective for the orthographic projection is at infinite distance. It depicts a hemisphere of the globe as it appears from outer space, where the horizon is a great circle. The shapes and areas are distorted, particularly near the edges.

The solar zenith angle is the zenith angle of the sun, i.e., the angle between the sun’s rays and the vertical direction. It is the complement to the solar altitude or solar elevation, which is the altitude angle or elevation angle between the sun’s rays and a horizontal plane. At solar noon, the zenith angle is at a minimum and is equal to latitude minus solar declination angle. This is the basis by which ancient mariners navigated the oceans.

The solar azimuth angle is the azimuth of the Sun's position. This horizontal coordinate defines the Sun's relative direction along the local horizon, whereas the solar zenith angle defines the Sun's apparent altitude.

<span class="mw-page-title-main">Tissot's indicatrix</span> Characterization of distortion in map projections

In cartography, a Tissot's indicatrix is a mathematical contrivance presented by French mathematician Nicolas Auguste Tissot in 1859 and 1871 in order to characterize local distortions due to map projection. It is the geometry that results from projecting a circle of infinitesimal radius from a curved geometric model, such as a globe, onto a map. Tissot proved that the resulting diagram is an ellipse whose axes indicate the two principal directions along which scale is maximal and minimal at that point on the map.

<span class="mw-page-title-main">Sunrise equation</span> Equation to derive time of sunset and sunrise

The sunrise equation or sunset equation can be used to derive the time of sunrise or sunset for any solar declination and latitude in terms of local solar time when sunrise and sunset actually occur.

<span class="mw-page-title-main">Cassini projection</span> Cylindrical equidistant map projection

The Cassini projection is a map projection first described in an approximate form by César-François Cassini de Thury in 1745. Its precise formulas were found through later analysis by Johann Georg von Soldner around 1810. It is the transverse aspect of the equirectangular projection, in that the globe is first rotated so the central meridian becomes the "equator", and then the normal equirectangular projection is applied. Considering the earth as a sphere, the projection is composed of the operations:

Atmospheric tides are global-scale periodic oscillations of the atmosphere. In many ways they are analogous to ocean tides. They can be excited by:

<span class="mw-page-title-main">Vibrations of a circular membrane</span> Equations of waves in a drumhead-like disc

A two-dimensional elastic membrane under tension can support transverse vibrations. The properties of an idealized drumhead can be modeled by the vibrations of a circular membrane of uniform thickness, attached to a rigid frame. Due to the phenomenon of resonance, at certain vibration frequencies, its resonant frequencies, the membrane can store vibrational energy, the surface moving in a characteristic pattern of standing waves. This is called a normal mode. A membrane has an infinite number of these normal modes, starting with a lowest frequency one called the fundamental mode.

<span class="mw-page-title-main">Kepler orbit</span> Celestial orbit whose trajectory is a conic section in the orbital plane

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.

<span class="mw-page-title-main">Geographical distance</span> Distance measured along the surface of the Earth

Geographical distance or geodetic distance is the distance measured along the surface of the Earth, or the shortest arch length.

<span class="mw-page-title-main">Poisson distribution</span> Discrete probability distribution

In probability theory and statistics, the Poisson distribution is a discrete probability distribution that expresses the probability of a given number of events occurring in a fixed interval of time if these events occur with a known constant mean rate and independently of the time since the last event. It can also be used for the number of events in other types of intervals than time, and in dimension greater than 1.

<span class="mw-page-title-main">Envelope (waves)</span> Smooth curve outlining the extremes of an oscillating signal

In physics and engineering, the envelope of an oscillating signal is a smooth curve outlining its extremes. The envelope thus generalizes the concept of a constant amplitude into an instantaneous amplitude. The figure illustrates a modulated sine wave varying between an upper envelope and a lower envelope. The envelope function may be a function of time, space, angle, or indeed of any variable.

In orbital mechanics, a frozen orbit is an orbit for an artificial satellite in which perturbations have been minimized by careful selection of the orbital parameters. Perturbations can result from natural drifting due to the central body's shape, or other factors. Typically, the altitude of a satellite in a frozen orbit remains constant at the same point in each revolution over a long period of time. Variations in the inclination, position of the apsis of the orbit, and eccentricity have been minimized by choosing initial values so that their perturbations cancel out. This results in a long-term stable orbit that minimizes the use of station-keeping propellant.

<span class="mw-page-title-main">Position of the Sun</span> Calculating the Suns location in the sky at a given time and place

The position of the Sun in the sky is a function of both the time and the geographic location of observation on Earth's surface. As Earth orbits the Sun over the course of a year, the Sun appears to move with respect to the fixed stars on the celestial sphere, along a circular path called the ecliptic.

In mathematics, a linear recurrence with constant coefficients sets equal to 0 a polynomial that is linear in the various iterates of a variable—that is, in the values of the elements of a sequence. The polynomial's linearity means that each of its terms has degree 0 or 1. A linear recurrence denotes the evolution of some variable over time, with the current time period or discrete moment in time denoted as t, one period earlier denoted as t − 1, one period later as t + 1, etc.

References