Pole Tide

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

The instantaneous earth rotation axis shifts inside a square of about [math]\displaystyle{ 20 }[/math] meters in relation to a point with fixed coordinates on the earth (i.e., Chandler wobble with a period of [math]\displaystyle{ 14 }[/math] months). This entails a varying elastic response of the earth's crust. This has an effect smaller than [math]\displaystyle{ 2.5 }[/math] centimetres in vertical and [math]\displaystyle{ 0.7 }[/math] centimetres in horizontal, but must be taken into account if the observations are carried out over periods longer than two months.

From the IERS Conventions [Denis et al., 2004] , pages 83-84, the following expression can be derived for the displacement at a point of geographic latitude [math]\displaystyle{ \varphi }[/math] and longitude [math]\displaystyle{ \lambda }[/math]:

[math]\displaystyle{ \begin{array}{ll} \delta_{\hat{\mathbf r}}=&-\frac{\omega^2_E \, R^2}{2\,g}\,h \, \sin 2\varphi \left ( m_1 \cos \lambda + m_2 \sin \lambda \right ) \,\hat {\mathbf r}\qquad \\[0.3cm] \delta_{\hat{\mathbf \lambda}}= & -\frac{\omega^2_E \, R^2}{g}\, l \, \sin \varphi \left (-m_1 \sin \lambda + m_2 \cos \lambda \right ) \, \hat{\mathbf \lambda} \qquad \\[0.3cm] \delta_{\hat{\mathbf \varphi}}= & -\frac{\omega^2_E \, R^2}{g}\, l\, \cos 2\varphi \left (\, m_1 \cos \lambda + m_2 \sin \lambda \right )\, \hat{\mathbf \varphi}\qquad \end{array} \qquad\mbox{(1)} }[/math]

where [math]\displaystyle{ (m_1,m_2) }[/math] are the displacements from the 1903.0 CIO, pole position, and [math]\displaystyle{ h= 0.6027 }[/math], [math]\displaystyle{ l=0.0836 }[/math] are the Love numbers.

Taking the earth's angular rotation [math]\displaystyle{ \omega_E= 7.29\cdot 10^{-5} }[/math], the earth's equatorial radius [math]\displaystyle{ R=6378\, km }[/math] and the gravitational acceleration [math]\displaystyle{ g=9.8\, m/sec^2 }[/math], it follows:

[math]\displaystyle{ \begin{array}{ll} \delta_{\hat{\mathbf r}}=&-32 \, \sin 2\varphi \left ( m_1 \cos \lambda + m_2 \sin \lambda \right ) \, \hat {\mathbf r}\qquad \mbox{(mm)} \\[0.3cm] \delta_{\hat{\mathbf \lambda}}= & -9 \, \sin \varphi \left (-m_1 \sin \lambda + m_2 \cos \lambda \right ) \, \hat{\mathbf \lambda} \qquad \mbox{(mm)} \\[0.3cm] \delta_{\hat{\mathbf \varphi}}= & -9 \, \cos 2\varphi \left (\, m_1 \cos \lambda + m_2 \sin \lambda \right )\, \hat{\mathbf \varphi}\qquad \mbox{(mm)} \end{array} \qquad\mbox{(2)} }[/math]

where [math]\displaystyle{ (m_1,m_2) }[/math] are the displacements (given in seconds of arc). Pole displacements can be found at ftp://hpiers.obspm.fr/iers/eop/eop.others.

The displacement [math]\displaystyle{ \delta }[/math] is given in the radial, longitude and latitude [math]\displaystyle{ (\hat{\mathbf r}, \hat{\boldsymbol \lambda}, \hat{\boldsymbol \varphi}) }[/math] vectors, i.e., the UEN local system, and, thence, can be transformed to the TRS reference system (see Reference Frames in GNSS) applying the following latitude and longitude rotations indicated in figure 1:

[math]\displaystyle{ \left[ \begin{array}{c} \Delta x \\ \Delta y \\ \Delta z \end{array} \right]= \underbrace{\left( \begin{array}{ccc} \cos \lambda \cos \varphi & -\sin \lambda & -\cos \lambda \sin \varphi \\ \sin \lambda \cos \varphi & \cos \lambda & -\sin \lambda \sin \varphi \\ \sin \varphi & 0 & \cos \Phi \end{array} \right)}_{{\mathbf R}_3(-\lambda )\cdot {\mathbf R}_2(\varphi )} \left[ \begin{array}{l} \delta_{\hat{\mathbf r}}\\ \delta_{\hat{\mathbf \lambda}}\\ \delta_{\hat{\mathbf \varphi}} \end{array} \right] \qquad\mbox{(2)} }[/math]

Figure 1: Transformation from UEN [math]\displaystyle{ (\hat{\mathbf r}, \hat{\boldsymbol \lambda}, \hat{\boldsymbol \varphi}) }[/math] to TRS [math]\displaystyle{ (x,y,x) }[/math] coordinates.