13.07.2013 Views

USNO Circular 179 - U.S. Naval Observatory

USNO Circular 179 - U.S. Naval Observatory

USNO Circular 179 - U.S. Naval Observatory

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

PRECESSION & NUTATION 41<br />

an intermediate step in the computations, or the final step where only very low accuracy is needed<br />

(10 arcseconds or worse) and nutation can be ignored.<br />

Thus, a precession transformation is applied to celestial coordinates to convert them from<br />

the mean equator and equinox of J2000.0 to the mean equator and equinox of another date, t.<br />

Nutation is applied to the resulting coordinates to transform them to the true equator and equinox<br />

of t. These transformations should be understood to be inherently geocentric rotations and they<br />

originate in dynamical theories. Generally we will be starting with celestial coordinates in the<br />

GCRS, which are obtained from basic ICRS data by applying the usual algorithms for proper<br />

place. 4 As discussed in Chapter 3, the ICRS is not based on a dynamically defined equator and<br />

equinox and so neither is the GCRS. Therefore, before we apply precession and nutation — and<br />

if we require a final accuracy of better than 0.02 arcsecond — we must first apply the frame bias<br />

correction (section 3.5) to transform the GCRS coordinates to the dynamical mean equator and<br />

equinox of J2000.0. Schematically,<br />

GCRS<br />

<br />

frame bias<br />

⇓<br />

mean equator & equinox of J2000.0<br />

<br />

precession<br />

⇓<br />

mean equator & equinox of t<br />

<br />

nutation<br />

⇓<br />

true equator & equinox of t<br />

Mathematically, this sequence can be represented as follows:<br />

r true(t) = N(t) P(t) B r GCRS<br />

where rGCRS is a direction vector with respect to the GCRS and rtrue(t) is the equivalent vector<br />

with respect to the true equator and equinox of t. N(t) and P(t) are the nutation and precession<br />

rotation matrices, respectively. The remainder of this chapter shows how to compute the elements<br />

of these matrices. B is the (constant) frame-bias matrix given in section 3.5.<br />

The transformation from the mean equator and equinox of J2000.0 to the mean equator and<br />

equinox of t is simply<br />

(5.4)<br />

r mean(t) = P(t) r mean(J2000.0)<br />

4 Computing proper place involves adjusting the catalog place of a star or other extra-solar system object for proper<br />

motion and parallax (where known), gravitational light deflection within the solar system, and aberration due to the<br />

Earth’s motions. For a solar system object there are comparable adjustments to its position vector taken from a<br />

barycentric gravitational ephemeris. See section 1.3. In conventional usage, an apparent place can be considered to<br />

be a proper place that has been transformed to the true equator and equinox of date. The details of the proper<br />

place computations are beyond the scope of this circular but are described in detail in many textbooks on positional<br />

astronomy, and in Hohenkerk et al. (1992), Kaplan, et al. (1989), and Klioner (2003).<br />

(5.3)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!