Earth orientation parameters explained

A rotation of a body in a 3D space in general is expressed by a three dimensional function of time. There exist more than one way to describe a rotation. The most commonly used way is to represent a rotation with three Euler angles.

Roughly speaking, the Earth's rotation can be considered as consisting of two components, the tidally driven component with precisely known frequencies and the component driven by an exchange of the angular momentum between the solid Earth and geophysical fluids. The latter component is not predictable in principle. The atmosphere contributes to the UT1 at a level of 10-6 rad, more than three orders of magnitude higher than the accuracy of observations. The first component is also affected by the atmosphere and can be predicted only at a level of 10-9 rad. Therefore, the Earth's rotation has to be continuously measured with modern space geodesy techniques.

The Earth rotation is mathematically expressed as a transformation of a vector in the rotating terrestrial coordinate system rt to the inertial celestial coordinate rc. This can be expressed as a product of the rotation matrix with a vector

rc =   M̂3(E3(t)) · M̂2(E2(t)) · M̂1(E1(t))   rt,

where, E1(t), E2(t), E3(t) are Euler angles with respect to axes 1,2,3 and M̂x(Ex) is a rotation matrix with respect to axis x:

                  1       0       0

       M1(E1) =   0  cos E1  sin E1

                  0 -sin E1  cos E1 
                                                 cos E2  0 -sin E2
                                      M2(E2) =   0       1       0

                                                 sin E2  0  cos E2 
                                                                               cos E3  sin E3  0

                                                                     M3(E3) = -sin E3  0  cos E3

                                                                               0       1       0
However, accordint to the adopted so-called Newcomb-Andoyer formalism, the rotation matrix is decomposed as a product of 12 elementary rotations:

M̂(t) = M̂3(ζt)) · M̂2(-θ(t)) · M̂3(z(t)) · M̂1(-ε0(t)) · M̂3(Δψ(t)) · M̂10(t) + Δε(t)) · M̂3(-S1(t) + E3(t)) · M̂2(E2(t)) · M̂1(E1(t)) · M̂3(H3(t)) · M̂2(H2(t)) · M̂1(H1(t))


The choice of Earth rotation parameterization is not logical, not economical, not optimal. This choice follows a historical tradition. Decomposition of a product of three matrix into a product of 12 matrix can be done by more than one way. There is an alternative decomposition Ginot-Capitaine. That decomposition is entirely equivalent to the Newcomb-Andoyer formalism.

NERS library keeps numerical coefficients of expansion ζ(t), θ(t), z(t); ε0(t), Δψ(t), Δε(t), and S1(t) and has the code that computes them on the specified moment of time. Empirical functions E(t) and H(t) are taken from the NERS server EOP message. Function E(t) comes as a table of values on specified, in general non-equidistant epochs. The tables are updated several times a day. Function H(t) comes in a form of a table of expansion coefficients determined from analysis of observations. It is updated 4–6 times a year. See NERS how for explanation how NERS server generates the EOP message.

NERS client automatically downloads the EOP message and extracts from there E(t) and H(t) functions relevant to the request, and computes the Earth's rotation matrix.

This page was last time edited by Leonid Petrov ()
Last update: [an error occurred while processing this directive]