Equation of time

The equation of time is the difference over the course of a year between time as read from a sundial and time as read from a clock, measured in an ideal situation (ie. in a location at the centre of a time zone, and which does not use daylight saving time). The sundial can be ahead (fast) by as much as 16 min 33 s (around November 3) or fall behind by as much as 14 min 6 s (around February 12). It is caused by irregularity in the path of the Sun across the sky, due to a combination of the obliquity of the Earth's rotation axis and the eccentricity of its orbit. 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 was used historically to set clocks. Between the invention of accurate clocks in 1656 and the advent of commercial time distribution services around 1900, the common way to set clocks was by observing the passage of the sun across the local meridian at noon. The moment the sun passed overhead, the clock was set to noon, offset by the number of minutes given by the equation of time for that date. [cite book
last = Olmstead
first = Dennison
authorlink =
coauthors =
title = A Compendium of Astronomy
publisher = Collins & Brother
date = 1866
location = New York
pages =
url = http://books.google.com/books?id=QUwAAAAAYAAJ&pg=PA57
doi =
id =
isbn = p.57-58] The equation of time values for each day of the year, compiled by astronomical observatories, were widely listed in almanacs and ephemerides. [cite book
last = Milham
first = Willis I.
authorlink =
coauthors =
title = Time and Timekeepers
publisher = MacMillan
date = 1945
location = New York
pages =
url =
doi =
id =
isbn = 0780800087 p.11-15] [See for example, cite book
last = British Commission on Longitude
first =
authorlink =
coauthors =
title = Nautical Almanac and Astronomical Ephemeris for the year 1803
publisher = C. Bucton
date = 1794
location = London, UK
pages =
url = http://books.google.com/books?id=CPgNAAAAQAAJ&pg=PT2
doi =
id =
isbn = p.14]

Naturally, other planets will have an equation of time too. 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.

"Apparent" time versus "mean" time

The irregular daily movement of the Sun was known by the Babylonians, and Ptolemy has a whole chapter in the "Almagest" devoted to its calculation (Book III, chapter 9). However he did not consider the effect relevant for most calculations as the correction was negligible for the slow-moving luminaries. He only applied it for the fastest-moving luminary, the moon.

Until the invention of the pendulum and the development of reliable clocks towards the end of the 17th century, the equation of time as defined by Ptolemy remained a curiosity, not important to normal people except astronomers. Only when mechanical clocks started to take over timekeeping from sundials, which had served humanity for centuries, did the difference between clock time and solar time become an issue. Apparent solar time (or true or real solar time) is the time indicated by the Sun on a sundial, while "mean solar time" is the average as indicated by clocks. The first accurate tables for the equation of time were published in 1665 by Christiaan Huygens. [cite book
last = Huygens
first = Christiaan
authorlink =
coauthors =
title = Kort Onderwys aengaende het gebruyck der Horologien tot het vinden der Lenghten van Oost en West
publisher = [publisher unknown]
date = 1665
location = The Hague
pages =
url = http://www.xs4all.nl/~adcs/Huygens/17/kort.html
doi =
id =
isbn = ] Following the practice of earlier astronomers Huygens increased his values for the equaton of time by a constant to make all values positive throughout the year. Another set of tables published in 1672 by John Flamsteed, the first head of the new Greenwich Observatory, was the first to adopt the modern convention using positive and negative values for the equation of time. [cite book
last = Flamsteed
first = John
authorlink =
coauthors =
title = De Inaequalitate Dierum Solarium
publisher = William Godbid
date = 1672
location = London
pages =
url =
doi =
id =
isbn = ]

Until 1833, the equation of time was mean minus apparent solar time in the British "Nautical Almanac and Astronomical Ephemeris". Earlier, all times in the almanac were in apparent solar time because time aboard ship was determined by observing the Sun. In the unusual case that the mean solar time of an observation was needed, the extra step of adding the equation of time to "apparent" solar time was needed. Since 1834, all times have been in mean solar time because by then the time aboard most ships was determined by marine chronometers. In the unusual case that the apparent solar time of an observation was needed, the extra step of adding the equation of time to "mean" solar time was needed, requiring all differences in the equation of time to have the opposite sign.

As the daily movement of the Sun is one revolution per day, that is 360° every 24 hours or 1° every 4 minutes, 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 30 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 itself 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 statistically insignificant at the accuracy level of sundials.

Eccentricity of the Earth's orbit

The Earth revolves around the Sun. As such it appears that the Sun moves in one year around the Earth. If the Earth orbited the Sun with a constant speed in a plane perpendicular to its axis, then the apparent Sun would culminate every day at exactly 12 o'clock, and be a perfect time keeper (except for its slowing rotation). But the orbit of the Earth is an ellipse and its speed varies between 30.287 and 29.291 km/s, according to Kepler's laws of planetary motion, and as such the Sun seems to move faster at perihelion (currently around 3 January) and slower at aphelion a half year later. At these extreme instances this effect increases (respectively, decreases) the real solar day by 7.9 seconds. This accumulates every day. The final result is that the eccentricity of the Earth's orbit contributes a sine wave variation with an amplitude of 7.66 minutes and a period of one year to the equation of time. The zero points are reached at perihelion (at the beginning of January) and aphelion (beginning of July) while the maximum values are at the beginnings of April (negative) and October (positive).

Obliquity of the ecliptic

However, even if the Earth's orbit were circular, the motion of the Sun along the celestial equator would still not be uniform. This is a consequence of the tilt of the Earth's rotation with respect to its orbit, or equivalently, the tilt of the ecliptic (the path of the sun against the celestial sphere) with respect to the celestial equator. The projection of this motion onto the 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 and appears as a change in right ascension, and is a minimum at the equinoxes, when the Sun moves in a sloping direction and appears mainly as a change in declination, leaving less for the component in right ascension, which is the only component that affects the duration of the solar day. As a consequence of that, the daily shift of the shadow cast by the Sun in a sundial, due to obliquity, is smaller close to the equinoxes and greater close to the solstices. At the equinoxes, the Sun is seen slowing down by up to 20.3 seconds every day and at the solstices speeding up by the same amount.

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 another sine wave variation with an amplitude of 9.87 minutes and a period of a half year to the equation of time. The zero points of this sine wave are reached at the equinoxes and solstices, while the maxima are at the beginning of February and August (negative) and the beginning of May and November (positive).

ecular 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: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 changes over the centuries due to secular variations in both eccentricity and obliquity. At this moment both are slowly decreasing, but in reality they vary up and down over a timescale of hundreds of thousands of years. When the eccentricity, now 0.0167, reaches 0.047, the eccentricity effect may in some circumstances overshadow the obliquity effect, leaving the equation of time curve with only one maximum and minimum per year, as it is on Mars [http://www.giss.nasa.gov/research/briefs/allison_02/] .

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 21 March (at least at sufficient accuracy for our aim here). The shift of the perihelion is forwards, about 1.7 days every century. For example in 1246 the perihelion occurred on 22 December, the day of the solstice. At that time the two contributing waves had common zero points, and the resulting equation of time curve was symmetrical. Before that time the February minimum was larger than the November maximum, and the May maximum larger than the July minimum. The secular change is evident when one compares a current graph of the equation of time (see below) with one from about 2000 years ago, for example, one constructed from the data of Ptolemy.

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 (usually) be the conic section of the hyperbola, since the circle of the Sun's motion together with the gnomon point define a cone. At the spring and fall equinoxes, the cone degenerates into a plane and the hyperbola into a line.dubious 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.

More details

In general, the equation of time is equal to

:$alpha\; -\; M\; +\; psi!,$,

where $alpha!,$ is the Sun's right ascension, $M!,$ is the mean anomaly, and $psi!,$ is the angle from the periapsis to the vernal equinox. Using spherical trigonometry, the right ascension is given by

:$cos(alpha)\; =\; cos(\; u\; -\; psi)\; /\; cos(delta)!,$,

where $u!,$ is the true anomaly and $delta!,$ is the Sun's declination. The declination in turn is given by

:$sin(delta)\; =\; sin(\; u\; -\; psi)\; sin(epsilon)!,$,

where $epsilon!,$ is the obliquity.

In practice, it may be easier and faster to use an approximation for the curve rather than the exact formula. For Earth, the equation of time resembles the sum of two offset sine curves, with periods of one year and six months respectively. It can be approximated by

:$E\; =\; 9.87sin(2B)\; -\; 7.53cos(B)\; -\; 1.5sin(B)!,$

where $E!,$ is in minutes and

:$B\; =\; 360^circ(N\; -\; 81)/364!,$ if sin and cos have arguments in degrees,

or

:$B\; =\; 2pi(N\; -\; 81)/364!,$ if sin and cos have arguments in radians.

Here, $N!,$ is the so-called day number; i.e.,

:$N=1$ for January 1,

:$N=2$ for January 2,

and so on.

The following is a graph of the current equation of time.

Notice that the appearance of this graph can be directly deduced from the time evolution of the projection into the celestial equator of the Earth's Analemma loop trajectory.

From one year to the next, the equation of time can vary by as much as 20 seconds, mainly due to leap years. [http://www.sunlit-design.com/infosearch/sundialaccuracy.php] .

References

* cite book
last = Meeus
first = Jean
authorlink =
coauthors =
title = Mathematical Astronomy Morsels
publisher = Willmann-Bell
date = 1997
location = Richmond
pages = 337-346
url =
doi =
id =
isbn = 0-943396-51-4

External links

* [http://www.sunlit-design.com/gallery/visualisations/eot1.php Graphical Visualisation of Equation of Time] - Constantly updated
* [http://freepages.pavilion.net/users/aghelyar/sundat.htm Table] giving the Equation of Time and the declination of the sun for every day of the year
* [http://www.nmm.ac.uk/server/show/conWebDoc.351 The equation of time] described on the Royal Greenwich Observatory website
* [http://www.analemma.com/ An analemma site with many illustrations]
* [http://myweb.tiscali.co.uk/moonkmft/Articles/EquationOfTime.html The Equation of Time and the Analemma] , by Kieron Taylor
* [http://astro.isi.edu/games/analemma.html An article by Brian Tung] containing a link to a C program using a more accurate formula than most (particularly at high inclinations and eccentricities). The program can calculate solar declination, Equation of Time, or Analemma.
* [http://www.phys.uu.nl/~vgent/astro/almagestephemeris.htm Doing calculations using Ptolemy's geocentric planetary models with a discussion of his E.T. graph]
* [http://www.sunlit-design.com/products/thesunapi/documentation/sdxEOT.php Equation of Time function] for Excel, CAD or other programs. The Sun API is free and extremely accurate. For Windows computers.
* [http://theorderoftime.org/science/equationoftime.html The equation of time correction-table] A page describing how to correct a clock to a sundial.
*An [http://www.audemarspiguet.com/index-equation.html example] of an Audemars Piguet mechanical wristwatch containing this concept as a complication, including a description of the implementation in horology and several videos/animations.
*Two more examples of a mechanical wristwatch containing this complication, manufactured by Blancpain: [http://www.blancpain.ch/e/best_articles/EM_part_I.html Part 1] [http://www.blancpain.ch/e/best_articles/EM_part_II.html Part 2] .
* [http://www.pendulumofmayfair.co.uk/view.asp?pid=272&cat=Longcase%20Clocks Equation of Time Longcase Clock by John Topping C.1720]
* [http://theorderoftime.org/truetime/solartime.html Solar tempometer] - Calculate your solar time including the equation of time.

Footnotes