Transformation between Celestial and Terrestrial Frames

Fundamentals
Title	Transformation between Celestial and Terrestrial Frames
Author(s)	J. Sanz Subirana, JM. Juan Zornoza and M. Hernandez-Pajares, University of Catalunia, Spain.
Level	Basic
Year of Publication	2011

Coordinate transformations between CRF and TRF frames are performed by mean of rotations corresponding to Precession, Nutation and Pole movement, briefly described as follows.

• Precession and nutation [forced rotation]: Earth rotation axis (and its equatorial plane) is not kept fixed in space, i.e., in relation to so called ”fixed-stars”, but it rotates about the pole of the ecliptic, as it is shown in figure 1. This movement is due to the effect of the gravitational attraction of the moon and the sun and major planets over the terrestrial ellipsoid. The total movement can be split into a secular component (precession, with a period of 26 000 years) and a periodic component (nutation, with a period of 18.6 years).
• Pole movement [free rotation]: Due to the structure of the earth’s mass distribution and its variation, the earth’s rotation axis is not fixed in relation to the earth’s crust. It moves around on the surface of the earth within a square of about 20 meters in relation to a point with fixed coordinates on earth. This movement has a period of about 430 sidereal days (Chandler period). On the other hand, earth rotation velocity is not constant, but it changes in time (although in very small quantities), as it was mentioned in the previous section.

The detailed expressions for transformation between the CRF and TRF frames are provided in Transforming celestial to terrestrial. The next equations briefly summarise such transformation: For a given epoch t, the coordinates transformation can be decomposed in a rotation matrices (i.e., orthogonal matrices) product as:

[math]\displaystyle{ \mathbf{[TRF]} = \mathbf{R_M}(t) \mathbf{R_S}(t)\mathbf{N}(t)\mathbf{P}(t)\mathbf{[CRF]} \qquad \mbox{(1)} }[/math]

being the inverse transformation:

[math]\displaystyle{ \mathbf{[CRF]} = \mathbf{P^T}(t) \mathbf{N^T} (t) \mathbf{R_S}^T (t) \mathbf{R_M}^T (t) \mathbf{[TRF]} \qquad \mbox{(2)} }[/math]

where: [CRF] Coordinates Vector in the Celestial Reference Frame. [TRF] Coordinates Vector in the Terrestrial Reference Frame. P Transformation matrix associated to the precession between the reference epoch and the epoch t. N Transformation matrix associated to the nutation at epoch t. R_S Transformation matrix associated to the earth rotation around the CEP axis. R_M Transformation matrix associated to the polar motion.

The matrices P and N are associated to the rotations needed to transform the coordinates from [CRF] to the [CEP]. They are provided by analytical expressions without requiring external parameters (see ICRF to CEP (link to article “ICRF to CEP”)), except for time conversions.

The matrices RS and RM are associated to the rotations needed to transform the coordinates from [CEP] to the [TRF]. Their computation requires the Earth Rotation Parameters files that are actualised periodically (see the website). More details can be found in CEP to ITRF (link to article “CEP to ITRF”)).

The transformation matrix for the polar motion is (see equation (link to equation eq:polar in the article CEP to ITRF)):

[math]\displaystyle{ RM(t) = R2(xp) R1(yp) \qquad \mbox{(3)} }[/math]

where xp and yp are the polar coordinates of the CEP in the TRF, and R1, R2 are the rotation matrices defined by (link to eq:TRANSF in the article Transformation between Terrestrial Frames). The transformation matrix associated to the earth rotation around the CEP axis is given by (see equation (link to equation eq:diurn in the article CEP to ITRF )):

[math]\displaystyle{ RS(t) = R3(ΘG ) \qquad \mbox{(4)} }[/math]

where ΘG is the Greenwich true sidereal time at epoch t and R3 is from (link to eq:TRANSF in the article Transformation between Terrestrial Frames).