The difference between the Uniform Dynamical Time and Universal Time is referred to as (delta ). Delta is used in numerous astronomical calculations, that is, eclipses,and length of day. It is additionally required to reduce quantified positions of minor planets to a uniform timescale for the purpose of orbital determination. Since Universal Time is established on the basis of the variable rotation of planet Earth, the quantity mirrors the unevenness of that rotation, and so it changes slowly, but rather irregularly, as time passes. We have worked on empirical formulae for estimating and have discovered a set of polynomials of the 4th order with nine intervals which is accurate within the range of ±0.6 seconds for the duration of years 1620–2013.

1. Introduction

The expression “timescale” is quite frequently used in astronomical contexts. To define it in astronomical terms, it may be put as a way of measuring time based on a particular periodic natural phenomenon. Two main distinct groups of timescales are used in astronomy. The first group of timescales is based on second which are known as International Atomic Time (IAT). It is the standard for the SI (System International) second. The -based timescales are comparatively new in the history of timekeeping, since they depend on atomic clocks that were first put to regular use in the 1950s era. Prior to that, all timescales were associated somewhat with the rotation of the Earth.

Timescales that rely on the rotation of the Earth are used for astronomical purposes as well. A relevant example would be a telescope pointing that relies on the geographic orientation of the observer. Universal Time mostly refers to the specific timescale . Historically, Universal Time (earlier known as the Greenwich Mean Time) has been achieved from Greenwich sidereal time using a general expression. However, is not fit for the computation of positions of the Moon, Sun, and planets using gravitational theories of their respective movements. Such theories prohibit variations in the rate of the rotation of Earth on its axis. Modern astrodynamical theories of the motions of the Sun, the Moon, and the planets are based on an evenly increasing and uniform timescale referred to as Terrestrial Time . runs a little ahead of (a refined measure of mean solar time at Greenwich) by an amount known as delta . As Earth’s rotation does not decelerate at a uniform rate, nontidal effects make it inconceivable to predict the precise values of in the distant past or remote future. Unfortunately, estimating the standard error in before 1600 is a tough task. It depends on several factors including the accuracy of determining from previous eclipse records and designing the models of physical processes creating changes in Earth’s rotation.

In the recent past, various polynomial representations for the values for the last few centuries have been suggested by Meeus [1], Islam et al. [2], and Meeus and Simons [3]. In this paper, we aim to present set of the 4th degree polynomials for delta with the least possible absolute error within seconds for the duration of 1620–2013.

2. Literature Review

Proving that the Earth rotates was no easy task for the scientists of former times. This problem goes back to the 17th century, when Halley [4] found that quadratic terms had to be added to the Moon’s mean longitude to match the times recorded for ancient eclipses. Laplace [5] announced that the acceleration term was due to perturbations from the Earth’s orbital eccentricity. Adams [6] then determined that Laplace had not included many higher order terms, which reduced Laplace’s final result to about half of Halley’s empirical value. To explain the source for the remaining observed effect, Ferrell [7] and Delaunay [8] independently assigned this discrepancy to tidal interactions between the Earth and the Moon. Newcomb [9] analyzed variations of the Earth’s rotation to explain some of the lunar residuals, but he could not obtain verification of these variations from inner planet data (Newcomb [10]).

In the latter part of the 19th century, Chandler discovered minute variations in Earth’s axis of rotation or polar wobble (Chandler [11, 12]). Newcomb then announced that these, and possibly other irregular variations of the Earth’s rotation, might be an explanation of the residuals in the lunar mean longitude (Newcomb [10]).

In all these cases, calculated positions were ahead of the observed positions by amounts proportional to the frequencies. Astronomers came to believe that discrepancies in ephemerides were not due to errors in the expressions for the mean longitude but were due to unmodeled irregularities and a deceleration of the Earth’s rotation, on which depends. Spencer Jones [13] examined residuals in the mean longitudes of the Sun, Moon, and two other planets and then concluded that the error was due to a slow deceleration in Earth’s rotation (Spencer Jones [14]). Since Newton’s theory of gravitation requires a nonaccelerating uniform timescale for the computation of orbital motions and because the Earth’s rotation was assumed to be decelerating, astronomers thought a timescale determined by the orbits themselves would be the uniform scale they needed. Thus, their proposition was a timescale called Ephemeris Time (), based on orbital motions, to be used for all dynamical calculations. In 1984, was replaced by Terrestrial Dynamical Time () as an independent argument for apparent geocentric ephemerides. The unit of is a day of 86 400 SI seconds at mean sea level. For practical purposes, is seconds. In 1991, was renamed Terrestrial Time (). is considered to be a uniform timescale and used as the time argument for the predictions of the astronomical events in dynamical theories (Seidelmann [15]; Islam et al. [2]).

Eclipse predictions are computed in . To convert predictions to , the difference between and must be in our knowledge. This parameter is known as or delta.

The mathematical modeling of the Earth’s delta has been ongoing for some time and there are multiple models to choose from. Before the advancement of AD 1600, values of had to rely on historical records of naked-eye observations of eclipses and occultations. Stephenson and a few other researchers have identified hundreds of these types of observations in early European, Middle Eastern, and Chinese annals, manuscripts, canons, and records. Moreover, some of these scientists have fit a lot of records with simple polynomials to achieve best fits for describing the value of for the centuries occurring before AD 1600 (Morrison and Ward [16]; Stephenson and Fatoohi [17]; Stephenson and Fatoohi [18]; Stephenson and Morrison [19]; Stephenson and Morrison [20]; Stephenson [21]; Morrison and Stephenson [22]; Steele and Stephenson [23]). Despite their relatively low precision and inaccuracy, these data represent our only record of the value of during the past several millennia.

Close-to-accurate values of only exist sometime after the invention of the telescope (1610). A careful analysis of telescopic timings of stellar occultation by the Moon permits the direct measurement of during this time period. In fact, all values of before 1955 depend on observations of the Moon, either via solar eclipses or via lunar occultations (Seidelmann [15]). Nowadays, is representative of the observed rotation angle of the Earth relative to an inertial reference frame formed by extragalactic radio sources. Its measurement, with the help of several observatories, is coordinated by the International Earth Rotation and Reference Systems Service ().

Numerous polynomial representations for the values of the last few centuries have been suggested in the recent past to evade the need of incorporating lengthy tables in a computer programs. Meeus [1] presents a 12th order polynomial liable for the time span of 1800–1997 with the maximum absolute error of seconds and two 9th or 10th order polynomials covering a similar time span with a maximum absolute error of seconds (Meeus [1]). Montenbruck and Pfleger [24] provide seven 3rd order polynomials to accommodate the period between 1825 and 2000 with maximum absolute error of 2.13 seconds. Meeus and Simons [3] give fourth-order polynomials with eight segments to cover the period between 1620 and 2000 with maximum absolute error of 3.2 seconds (Meeus and Simons [3]; Islam et al. [2]). Islam et al. [2] presented eight 4th degree polynomials to the entire curve of with precisely the same intervals as were used by Meeus and Simons [3] compromising maximum absolute error of seconds on time span from 1620 to 2000 .

3. New Proposed Polynomial Approximation

The data used in this research was excerpted from the pages and of the Astronomical Almanac published by the Nautical Almanac Offices of the US Naval Observatory and the Bulletins issued by the International Earth Rotation Service, , Paris. For a new set of polynomials, we apply method of the least square on the curve of delta to get appropriate polynomials for epoch 1620–2013 (Figure 1). The quantity , as defined by Meeus and Simons [3], is

In this equation, the purpose of the quantities is basically to make the independent variable as small as possible during a certain interval of time, so we simply call it the “scaling factor.”

The empirical formulae that we found here can calculate the value of at any instant of the years 1620–2013 . The 4th degree polynomial approximation of is where , and are coefficients of respective polynomials taken in seconds, given in Table 1.

4. Discussion of Results and Conclusion

The value of scaling parameter and coefficients of the 4th degree polynomial with nine intervals that we have calculated are displayed in Table 1. It is clear from the table that the values of and coefficients are only valid for adjacent intervals. We choose the 4th degrees polynomials in order to compare our results with previous approximations (Meeus and Simons [3] and Islam et al. [2]).

Figure 2 shows error analysis of the 4th degree polynomial with respect to its nine segments. Error analysis revealed that, in the 4th degree polynomial, minimum error of 0.000168 occurs in the year 1712, while maximum error of 0.598961 appears in the year 1692.

The maximum error found in this study is so minute that it turns out to be less than all sets of polynomials presented in Montenbruck and Pfleger [24], Meeus and Simons [3], and Islam et al. [2]. Therefore, it is concluded that a set of nine fourth degree polynomials has proven to be the most appropriate set covering a time span of 1620–2013 with seconds of accuracy. The time accuracy of astronomical events can predict the timescale in a much better fashion with this polynomial approximation. The uses include the solving of predicaments ranging from astronomical software to eliminating the steps in the search of the long table of the observed values of during the aforementioned period.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


The authors would like to thank Ms. Wishaal Khalid and other, anonymous, reviewers for their valuable comments and suggestions for the improvement of the quality of the paper.