A&A 428, 255-260 (2004)

DOI: 10.1051/0004-6361:20035952

**S. Lambert ^{1} - N. Capitaine^{2}**

1 - Department of the Navy, US Naval Observatory, Earth Orientation Department, 3450 Massachusetts Avenue NW, Washington DC 20392, USA

*On leave from* 2

2 -
SYRTE - UMR 8630/CNRS - Observatoire de Paris, 61 avenue de l'Observatoire, 75014 Paris, France

Received 29 December 2003 / Accepted 9 August 2004

**Abstract**

Because of their importance in the accurate modeling of the Earth's orientation in space, some non-negligible predictable effects on precession-nutation are investigated. This paper considers the coupling effects between the axial and the equatorial components of the Earth's rotation vector in the dynamical equations, and the effects of the second order lunisolar torque due to the Earth's zonal deformations. Firstly, the coupling effects are shown to contribute for less than 0.1 as and are therefore negligible. Secondly, we demonstrate that the 0.7 mas contribution of the rotation rate variations due to zonal tides to the nutation in obliquity deduced by Bretagnon et al. (2000, Proc. IAU Coll., 180, 230; 2001, Celest. Mech. Dyn. Astr., 80, 177) is an artefact which comes from an incomplete way of taking into account the effect of the rotation rate variations. The net contribution is shown to be negligibly small. Thirdly, for an Earth model with an elastic mantle and decoupled liquid core, the contribution of the second-order lunisolar torque due to the Earth's zonal deformations is shown to be 207.9 as and -9.7 as on the 18.6-year nutation respectively in longitude and in obliquity, and a correction of -4925.9 as/century on the precession in longitude.

**Key words: **reference systems - astrometry - time

Very long baseline radio interferometry (VLBI) observations provide on a regular basis "celestial pole offsets'' that contain errors in the models for the position of the celestial pole in the celestial reference system. The most accurate available precession-nutation model is the IAU 2000A model recommended by the International Astronomical Union (IAU Resolution B1.8) adopted at the IAU 24th General Assembly (2000) and implemented by the IERS beginning on 1 January 2003.

Comparison of this new model against VLBI observations shows differences of the order of 200 as. These differences are the consequence of various influences as well as geophysical processes that are summarized below. The main oscillation appearing in these residuals comes from the Free Core Nutation (FCN). It gives rise to a 430-day periodic term, the amplitude of which varies strongly and is therefore not predictable but can be described using an empirical model (Herring et al. 2002). As the nutation modeling contains only permanent predictable effects, the FCN is not included in IAU 2000A and therefore appears in the residuals. Other unpredictable effects come from the variable atmospheric and oceanic forcing which induces time varying nutation amplitudes of about 80 as for the annual term (Bizouard et al. 1998). Some theoretical predictable effects are taken into account in the IAU 2000A precession-nutation model whereas they were not included in the previous models. These effects include:

- the couplings between axial and equatorial components appearing when developing the dynamical equations of the Earth rotation at the second order;
- the effects of variations of the dynamical ellipticity, giving rise to changes in the lunisolar torque exerted on the Earth and therefore on precession-nutation.

The deformations due to zonal tides have also direct effects on precession and nutation. The effect on nutation was first studied by Souchay & Folgueira (2000). It was concluded that, for a basic Earth model, this should modify the amplitude of the 18.6-year nutation in longitude by 168 as, perturbations on other nutation terms all being less than 10 as. Another work by Krasinsky (1999) showed larger corrections based on a computation of the gravitational torque exerted on the Earth. The largest correction provided is 5.11 mas on the 13.66-day term, whereas the semi-annual term and the 18.6-year term are changed respectively by 2.38 mas and 1.07 mas. Note that VLBI observations do not show such large discrepancies in the nutation model.

Considering such a large disparity in the amplitudes of the contributions to nutation due to zonal tides that are reported in various studies, a rigorous computation of these effects is necessary. In this paper, we compute the effects due to the torque produced by the action of the tesseral potential on the deformations produced by the zonal potential together with the effects on the rotation rate produced by these deformations. In Sect. 2 we recall the expressions of the Euler-Liouville equations developed up to the second order, emphasing the couplings between variations of the axial component and the equatorial component of the instantaneous rotation vector. In Sect. 3 we introduce the second order part of the lunisolar torque coming from Earth's zonal deformations, and in Sect. 4, we compute all the predictable effects due to zonal tides on the Earth's precession-nutation.

2 Dynamical equations of the Earth's rotation up to the second order

We consider a deformable Earth with fluid core. In a rotating frame with its axes oriented towards the Earth's mean axes of inertia, the angular momentum equation is:

where is the whole Earth's angular momentum vector, corresponding to the product of the inertia tensor :

wherein second-order quantities

(3) |

being the mean Earth's rotation rate.

is the tidal torque expressed in the terrestrial frame. The equatorial component
of the torque is computed from the expression
in the celestial frame by the relation:

and the Euler's kinematical relations:

relate the terrestrial motion of the Earth's axis of rotation and the relative variations of the Earth's rotation rate

Using the equations of Sasao et al. (1980), one obtains the following expressions:

in which , , and , quantities subscripted by f being relevant to the core.

The previous equations reflect two different couplings between the axial and the equatorial component of the rotation vector of the Earth. One coupling comes from the products
or
in the left hand side of (6). Another coupling appears through the sidereal rotation angle
depending upon *m*_{3}, in Eqs. (4) and (5) in the opposite sense to transform coordinates between the terrestrial to the celestial frame.

Variations of Euler's angles are related to the instantaneous vector of rotation by Euler's kinematical relations (5). One gets the nutation angles referred to the ecliptic of epoch:

= | (7) | ||

= |

Note that the periodic part of the nutation angles referred to the ecliptic of date ( and ) is the same as the one of and , and that the precession is expressed along the fixed ecliptic. The contribution to the precession is equivalent to the secular term of .

3 Torque including Earth's zonal deformations

The tidal gravitational torque exerted on the Earth is derived from the tesseral part of the lunisolar potential
in the celestial reference frame (Sasao et al. 1980):

where the superscripts (1) and (2) stand for the first and second order parts of the torque respectively. The second order part includes the deformations

The complex-valued quantity
depends on the right ascension
and declination
of the perturbing body:

(9) |

in which

**Table 1:**
Excess of LOD (s) due to zonal tides from IERS Conventions 2003, relative variations of the Earth's rotation rate *m*_{3}^{z} (10^{-12} rad) and *c*_{33}/*C* (as) from Eqs. (11) and (12).

(10) |

The variations of the dynamical ellipticity are related to corresponding changes in the excess of length-of-day (LOD) and to the axial component

where the superscript

and being undimensional parameters representing the deformabilities of the mantle and the core-mantle boundary under degree 2 tidal forcing (see e.g. Sasao et al. 1980; Mathews et al. 1991) and the axial moment of inertia of the core.

The second order torque due to the zonal deformations is finally expressed using quantity *c*_{33} or equivalently *m*_{3}^{z}:

Direct computation of the changes in the dynamical ellipticity can be based on expression for

Our computation of *c*_{33} uses this model for LOD variations and the mean moments of inertia of the whole Earth provided by the IERS Conventions 2003:

(14) | |||

Using the compliances from Mathews et al. (2002) estimates, one gets:

(15) |

Table 1 gives the excess of the length-of-day due to tidal deformations provided by the IERS Conventions 2003 together with the variations

4 Solution of the dynamical equations

All the terms showing products or in the left-hand side of Eq. (6) express coupling between the motion of the axis of rotation and the Earth's rotation rate (note that the use of a decoupled core involves the relation ). The largest contribution of these terms to nutation has been evaluated as being lower than 0.1 as for the 18.6-year term. This confirms that this coupling effect is negligible and it will not be considered further in the following.

The effect of taking into account the Earth's rotation irregularities in the sidereal rotation angle in computing the solution for precession-nutation has to be investigated.

Bretagnon et al. (2000, 2001) constructed non-rigid nutation series by applying the transfer function MHB 2000 (Mathews et al. 2002) to the SMART97 rigid nutation series (Bretagnon et al. 1998). Given that the transfer function of MHB 2000 does not take into account the variations in the Earth's rotation rate, these effects were added to the third component of the Earth's rotation vector and the Euler's angles were recomputed from Euler's kinematical relationships. This method gave rise in particular to a 18.6-year term in obliquity with an amplitude reaching 712 as.

We investigate this effect considering the equations for a rigid ellipsoidal Earth (the quantity *c* is null) together with variations in Earth's rotation rate. Equation (6) is reduced to:

and therefore:

where is a constant depending on initial values of the problem.

The first term of solution (17) is the free motion of the rotation axis of the Earth, known as the Euler motion. The second term of Eq. (17) corresponds to the forced motion of the rotation axis. Variations of the Euler's angles are therefore:

Using an iterative process, the integration of the system (18) converges rapidly after a few iterations. Note that the sidereal rotation angle is affected by the irregularities

By a numerical evaluation, one can see that omitting the variations in the Earth's rotation rate in the two rotations of angle plus or minus , that is to say writing instead of , leads to an error below 0.1 as on and . The effect of the former rotation is actually mostly cancelled out by the latter so that the net effect is negligibly small.

However, omitting to take into account the Earth's rotation rate irregularities in one of the rotations of plus or minus
leads to large terms which could be interpreted as being coupling effects. For example, if the computation of the torque in the terrestrial frame neglects these irregularities, one gets a substantial contribution in the in-phase 18.6-year term of the obliquity (-729 as) and on the 9.3-year in-phase term of the longitude (71 as). This 18.6-year oscillation of the angle
is due to a coupling between the precession in longitude
and the 18.6-year term in the variations of the Earth's rotation rate of amplitude *a*_{18.6}:

(19) |

The precession in longitude is also affected by a spurious effect of 4 mas/c. This explains the erroneous values provided in the above-mentioned studies that followed such an approach. These "artefacts" are displayed in Table 2 together with the terms provided by Bretagnon et al. (2000, 2001) and the net effect.

**Table 2:**
Effect of zonal variations in Earth's rotation: artefacts induced by omitting the variations when transforming the torque from the celestial frame to the terrestrial frame compared to the values of Bretagnon et al. (2001), and net contribution.

The computation of the effects of the second order torque induced by zonal tides on nutation is done using Eq. (6) in the frequency domain once the torque is computed according to Eq. (13). The calculation procedure of the lunisolar torque made in the present paper is based on the lunar theory ELP2000 (Chapront-Touzé & Chapront 1983) and the solar system semi-analytical solution VSOP87 (Bretagnon & Francou 1988). Poisson series manipulations have been processed using the GREGOIRE software package developed by Chapront (2003). We use standard values (IERS Conventions 2003) for other orbital parameters involved in the computation of the torque.

The computed variations of the nutations angles are the following:

**(i) Secular term.**- The constant term of the torque leads to a secular variation in longitude of -4926 as/c. Note that such a contribution would be automatically included in the observed value for the precession. The obliquity is not affected.
**(ii) Periodical terms.**- The main second-order contribution to the periodic terms comes from the 18.6-year term of the torque which is a coupling between the 18.6-year variation of the dynamical ellipticity and the constant term of the lunisolar potential. Other contributions are below 5 as. It is interesting to note that, according to Table 1, the 18.6-year term is in the opposite sense compared to the others. Namely, it is negative and decreases the excess of LOD (or equivalently, decreases the difference between axial and equatorial moment of inertia so that the Earth becomes less elliptic). This means that the effect of the 18.6-year variation is to concentrate the masses of the Earth around the rotation axis, so that the Earth rotation speed increases to satisfy the angular momentum conservation law, contrarily to the other contributions which are positive and make the Earth oblateness increase. A consequence of the oblateness decrease is that the torque exerted on the equatorial bulge is lowered and the 18.6-year nutational response to this excitation will consequently be lowered. The 18.6-year nutation coefficient should be decreased in absolute value. In fact, this contribution reaches 208 as in the nutation in longitude whereas it produces only a small effect of almost -10 as on the nutation in obliquity. Since the 18.6-year nutation amplitude in longitude is negative (-17 ), such a result is in the right direction.

Mathews et al. (2002) investigated a refined Earth model (elastic mantle and decoupled liquid core) to evaluate the influence of the so-called "non-linear'' terms in the dynamical equations in order to remove them from the observations before fitting the geophysical parameters used in the transfer function. This gave a contribution to the 18.6-year term of 94 as and -29 as respectively in longitude and in obliquity, including the second order torque due to both zonal and sectorial tides. Note that the differences between the values obtained in this study and those of Mathews et al. (2002) appearing in Table 3 are in fact considerably reduced with respect to revised values obtained by Mathews (2004). Excluding the sectorial part, the revised contributions in the 18.6-year nutation are 194 as and -10 as respectively in longitude and in obliquity (Mathews 2004). The remaining difference with respect to our values is due to anelasticity and ocean tide effects taken into account in the MHB work. It should moreover be noted that since this work has been carried out, Mathews (2003) pointed out that the above computed effects could be cancelled out by the sectorial and zonal parts of the potential acting on the deformations due to the tesseral potential. This point needs further investigation.

**Table 3:**
Effects of the second order contribution of the tidal potential (zonal tides contribution) on nutation angles (as) and precession (as/c) from different studies. MHB 2000Z: values of Mathews et al. (2002) including only the zonal tides contribution; SF99: Souchay & Folgueira (2000).

Previous studies (Bretagnon et al. 2000, 2001) pointed out that the Earth's rotation rate variations due to zonal tides could have a noticeable effect on precession-nutation, with an amplitude of the order of 700 as in the 18.6-year nutation in obliquity, although such an effect was not detected in the observations.

The present paper has shown that this effect is in fact cancelled out by taking into account the variations of the rotation rate in the expression of the torque as seen from the terrestrial frame used in the dynamical equations. This point resolves a serious dilemma of several years concerning the magnitude of these effects. Moreover, this term does not contain the effects induced by the variations of the dynamical ellipticity, although it may sometimes have been understood as being so.

We also recomputed the effects due to the coupling between axial and equatorial components of the rotation vector in the dynamical equations of the Earth rotation and we took into account the second-order contribution of the torque due to zonal deformations of the Earth's shape. These results are obtained by solving the dynamical Eqs. (6) using the model for variations in the rotation rate due to zonal tides of the IERS Conventions 2003. For an elastic Earth with a decoupled liquid core, we have concluded that:

- the coupling effects between axial and equatorial components of the rotation vector do not bring any contribution larger than 0.1 as and are therefore negligible;
- the second-order torque induced by Earth's zonal deformations produces larger effects. The effect on the 18.6-year nutation is of 207.9 as in longitude and -9.7 as in obliquity. Other contributions are about a few microarcseconds. The correction to the precession in longitude is -4925.9 as/c.

The authors are grateful to J. Souchay and P. M. Mathews for valuable discussions and to J. Chapront for having kindly provided the software package GREGOIRE. They thank the anonymous referee for very useful suggestions in improving the manuscript.

- Bizouard, C., Brzezinski, A., & Petrov, S. D. 1998, J. Geodesy, 72, 561 In the text
- Bretagnon, P., & Francou, G. 1988, A&A, 202, 309 In the text
- Bretagnon, P., Francou, G., Rocher, P., & Simon, J.-L. 1998, A&A, 329, 329 In the text
- Bretagnon, P., Rocher, P., & Simon, J.-L. 2000, Proc. IAU Coll. 180, ed. K. J. Johnston et al. (Washington DC: US Naval Obs.), 230 In the text
- Bretagnon, P., Rocher, P., & Simon, J.-L. 2001, Celest. Mech. Dyn. Astr., 80, 177 In the text
- Chapront-Touzé, M., & Chapront, J. 1983, A&A, 124, 50 In the text
- Chapront, J. 2003, Notice for GREGOIRE, Paris Observatory In the text
- Defraigne, P., & Smits, I. 1999, Geophys. J. Int., 139, 563 In the text
- Hefty, J., & Capitaine, N. 1990, Geophys. J. Int., 103, 219 In the text
- Herring, T., Mathews P. M., & Buffett, B. A. 2002, J. Geophys. Res., 107, B4, 10.1029/2001JB000165 In the text
- IERS Conventions 2003, IERS Technical Note 32, ed. D. D. McCarthy & G. Petit (Frankfurt am Main: Verlag des Bundesamts für Kartographie und Geodäsie, 2004, 127 pp.), in press In the text
- Krasinsky, G. A. 1999, Celest. Mech. Dyn. Astr., 75, 39 In the text
- Mathews, P. M., Buffett, B. A., Herring, T. A., & Shapiro, I. I. 1991, J. Geophys. Res., 96, B5, 8243 In the text
- Mathews, P. M., Herring, T. A., & Buffett, B. A. 2002, J. Geophys. Res., 107, B4, 10.1029/2001JB000390 In the text
- Mathews, P. M. 2003, AGU Fall Meeting In the text
- Mathews, P. M. 2004, private communication In the text
- McCarthy, D. D. 1996, IERS Conventions 1996, Observatoire de Paris In the text
- Melchior, P. 1978, The Tides of the Planet Earth (Pergamon Press) In the text
- Sasao, T., Okubo, S., & Saito, M. 1980, Proc. IAU Symp. 78, ed. E. P. Federov, M. L. Smith, & P. L. Bender (Hingham, Mass.: D. Reidel), 165 In the text
- Souchay, J., & Folgueira, M. 2000, Earth, Moon Planets, 81, 201 In the text
- Yoder, C. F., Williams, J. G., & Parke, M. E. 1981, J. Geophys. Res., 86, 881 In the text

Copyright ESO 2004