Welcome to the NavList Message Boards.


A Community Devoted to the Preservation and Practice of Celestial Navigation and Other Methods of Traditional Wayfinding

Compose Your Message

Add Images & Files
    Delta T long term formula
    From: Paul Hirose
    Date: 2018 Feb 3, 09:38 -0800

    A 2016 paper (Stephenson FR, Morrison LV, Hohenkerk CY., "Measurement of
    the Earth’s rotation: 720 BC to AD 2015," Proc. R. Soc. A 472: 20160404,
    http://dx.doi.org/10.1098/rspa.2016.0404) provided a set of polynomial
    approximations for ∆T, the time scale difference TT-UT1.
    The polynomials are valid from -720.0 to 2016.0. Outside that span we
    can use the parabola that best fits the whole data set: if t = (year -
    1825.0) / 100, then ∆T = -320 + 32.5 * t^2.
    Although a constant deceleration of Earth's rotation implies a parabolic
    graph of ∆T, a parabola may be a poor fit in the short term. For
    example, figure 10 of the paper (a graph of ∆T from 1550 to the present)
    doesn't show the parabola because it's off the graph. At 2018.0 the
    parabolic formula gives -200 s, whereas the true value is +69.
    We can get more reasonable values by starting at the last polynomial,
    which ends at 2016.0, when ∆T = 68.041. Then extrapolate from that point
    by integrating the length of day expression in the paper. That's how the
    values on the UKHO page ("ΔT & lod from −2000 to +2500") were obtained.
    It doesn't give a formula for the integral, but I have derived one,
    where t = year - 1825.0. To extrapolate into the future from 2016.0:
    ∆T = -293.600 + .00325 * t^2 + 349 * cos(.00419 * t)
    Note t is years, not centuries. The input to the cosine function is
    radians, not degrees.
    To extrapolate into the past from -720.0, change the constant from
    -293.600 to -385.979. With those constants there is less than one
    millisecond discontinuity at the junctions with the polynomials, and the
    extrapolations back to -2000 and forward to +2500 are practically
    identical (within the estimated error) to the values on the UKHO page.

    Browse Files

    Drop Files


    What is NavList?

    Join NavList

    (please, no nicknames or handles)
    Do you want to receive all group messages by email?
    Yes No

    You can also join by posting. Your first on-topic post automatically makes you a member.

    Posting Code

    Enter the email address associated with your NavList messages. Your posting code will be emailed to you immediately.

    Email Settings

    Posting Code:

    Custom Index

    Start date: (yyyymm dd)
    End date: (yyyymm dd)

    Visit this site
    Visit this site
    Visit this site
    Visit this site
    Visit this site
    Visit this site