A&A 428, 255-260 (2004)
S. Lambert1 - N. Capitaine2
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
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 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.
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:
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:
Using the equations of Sasao et al. (1980), one obtains the following expressions:
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 m3, 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:
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):
The complex-valued quantity
depends on the right ascension
of the perturbing body:
Table 1: Excess of LOD (s) due to zonal tides from IERS Conventions 2003, relative variations of the Earth's rotation rate m3z (10-12 rad) and c33/C (as) from Eqs. (11) and (12).Assuming as generally (see e.g. Melchior 1978; Yoder et al. 1981) that the trace of the inertia tensor is constant, one has in the case of zonal deformations:
The second order torque due to the zonal deformations is finally expressed using quantity c33 or equivalently m3z:
Our computation of c33 uses this model for LOD variations and the mean moments of inertia of the whole Earth provided by the IERS Conventions 2003:
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:
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:
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 a18.6:
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:
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 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.