# Study of atmospheric and oceanic excitations in the motion of Earth's spin axis in space.

1. INTRODUCTIONThere are five Earth Orientation Parameters (EOP) that describe the full orientation of the Earth in space: two components of polar motion (fixing the position of the spin axis in the body), the angle of proper rotation around the spin axis, and two components of the celestial pole offsets (describing the misalignment of the spin axis with its position in space given by a precession-nutation model). In this study, we concentrate on the latter two. Thanks to a very accurate present IAU models of nutation (IAU2000--Mathews et al., 2002) and precession (IAU2006--Capitaine et al., 2003), these values are very small, typically around 0.1 milliarcsecond (mas), their dominant part being the retrograde Free Core Nutation (FCN) with the period of about 430 days. This motion is due to a strong resonance in a near-diurnal part of the spectrum (in terrestrial frame) in celestial motion of the Earth's spin axis, caused by the flattened fluid outer core of the Earth. This resonance, together with a partial viscosity of the mantle, leads to a significant modification of amplitudes and phases of the forced (by external torques exerted by the Moon, Sun and planets) nutation terms, with respect to the solution calculated for a rigid Earth model. The largest influence is observed for the retrograde annual term that is closest to the resonance.

Atmospheric and oceanic excitations have dominant effects in polar motion and proper rotation of the Earth, with significant seasonal periods in terrestrial frame. Their near-diurnal part, which is reflected into long periods in celestial frame, is rather weak. However, these small excitations are amplified by a resonance so that they become detectable with the recent observational techniques. Availability of atmospheric and oceanic excitation functions, defined by Barnes et al. (1983), with 6-hour resolution (that appeared more than ten years ago) enabled recently to study these effects in detail (see, e.g., Bizouard et al., 1998; Bizouard, 1999; Yseboodt et al., 2002; Brzezinski et al., 2002; Lambert, 2006; Vondrak and Ron, 2006 a, b, 2007, and 2008, or Vondrak, 2009). These studies show that the effect is most significant in annual and semi-annual terms of nutation, and is of the order of a hundred microarcseconds ([micro]as). The first part of the present study is devoted to this problem.

The value of the FCN period, used to derive the older model of nutation IAU1980, was equal to about 460 days (Wahr, 1981). It corresponded to the case where core is in hydrostatic equilibrium. The observations by Very Long-Baseline Interferometry (VLBI), so far the only observational technique capable of observing a celestial motion of the Earth spin axis with a sufficient accuracy, revealed that the period is somewhat smaller, only about 430 days (Roosbeek et al., 1999; Hinderer et al., 2000; Mathews et al., 2002; Vondrak and Ron, 2005; Vondrak et al., 2005). This corresponds to the core's flattening about 4 per cent larger than the hydrostatic equilibrium value (Herring et al., 1986; Gwinn et al., 1986). The question if the period P and quality factor Q of FCN are stable in time was recently addressed, e.g., by Lambert and Dehant (2007). They concluded that the resonant period is stable within less than half a day, but the differences in approach of different analysis VLBI centers have an impact of the same order of magnitude. They also found that the contribution of the atmosphere is negligible. We addressed the same problem recently (Vondrak and Ron, 2009) and found that all IVS analysis centers yield similar quasi-periodic changes that do not disappear when atmospheric and oceanic contributions are removed. In addition, different models of atmosphere/ocean forcing result in systematically different results. Here, we re-visit the problem of stability of P and Q and study the role of atmospheric/oceanic excitations in more detail.

2. THE DATA

In this work, we use the following data:

* Since different IVS analysis centers yield very similar results (Vondrak and Ron, 2009), we use the values of celestial pole offsets (CPO) only from the most recent IVS combined solution (Schluter and Behrend, 2007) ivs09q3X.eops, covering the interval 1984.1-2009.7. The data in the form dX, dY are given in unequally spaced intervals, 1-7 days long, sometimes with large outliers. Therefore, we first cleaned the data by removing CPO values exceeding in absolute value 1 mas, and then interpolated them to regular 3-day intervals, by using a weak smoothing and cubic spline function.

* Atmospheric and oceanic angular momentum functions are now available from different sources. Since the most important part of excitation is in near-diurnal part of the spectrum (in terrestrial frame), we need the data sampled with sub-diurnal interval. We use the following data, available in 6-hour intervals directly from the authors:

* Atmospheric angular momentum excitation functions (pressure + wind terms):

** NCEP/NCAR re-analysis, 1983.0-2009.5 (Salstein, 2005);

** ERA, 1979.0-2009.0 (Thomas et al., 2006; Dobslaw and Thomas, 2007; re-analysis model before 2001.0, operational model afterwards);

* Oceanic angular momentum excitation functions (matter + motion terms):

** ECCO model, 1993.0-2009.7 (Gross et al., 2005);

** OMCT model, 1979.0-2009.0 (Thomas et al., 2006; Dobslaw and Thomas, 2007, driven by re-analysis atmospheric model before 2001.0, by operational model afterwards).

These time series of [chi] (complex numbers) are given in terrestrial frame, so we have to transform them into the celestial (non-rotating) frame values [chi]. To do so, we first removed a constant part, which would lead to a big diurnal signal (with the period equal to a sidereal day) in celestial frame. Then a relatively simple formula in complex form [chi]' = -[chi][e.sup.i[PHI]] is applied, where [PHI] is the Greenwich sidereal time. The near-diurnal variations in terrestrial frame become long-periodic in celestial frame. Because we are interested in only long-periodic motion, that is comparable to nutation frequencies, we further apply the smoothing (Vondrak, 1977) to remove all periods shorter than 10 days. It is also necessary to say that the ECCO model is forced by the NCEP/NCAR data, so these two should be used together; a similar argument holds also for OMCT/ERA combination.

3. OBSERVED CELESTIAL POLE OFFSETS AND THEIR EXCITATION

The observed celestial pole offsets in mas, brought to regular 3-day intervals, are depicted in Figure 1. Both components display a clear quasi-periodic character, the accuracy of the results are evidently improving in time. The FFT spectral analysis, shown in Figure 2 reveals several dominant peaks. The largest one corresponds to the retrograde FCN; its double peak is the sign of the fact that both amplitude and phase of this motion are unstable. Significant are also long-periodic motions with periods of about 14 and 28 years. These are probably caused by mismodeled long-periodic nutation, but surely not by geophysical fluids (see below). Smaller peaks at prograde and retrograde annual (and to some extent also at semi-annual and terannual) periods reflect small differences from the adopted nutation model. These can be caused either by deficiency of the nutation model itself, or by additional geophysical excitations.

The spectra of atmospheric and oceanic angular momentum excitation functions are depicted in Figures 3 through 6. If we compare the first two, showing the atmospheric excitations from NCEP and ERA, there are practically no difference. The most prominent peaks are at prograde annual and semi-annual periods, less visible are retrograde excitations. Wind terms are much larger than the pressure ones. The situation is different for oceanic excitations, displayed in Figures 5 and 6. The peaks are smaller than for the atmosphere, and the differences between matter (equivalent to pressure) and motion (equivalent to wind) terms are almost negligible. OMCT peaks are however significantly higher than ECCO ones. It is caused by the fact that ECCO model is forced by the atmosphere only once per day while OMCT uses 6-hour forcing. Thus the ECCO model is somehow smoothed, with the subdiurnal signal (in terrestrial frame) suppressed. If compared to Figure 2, the atmospheric and oceanic excitations evidently lack a long-periodic part with decadal periods, and also the part with 430-day period. The first fact hints that celestial pole offsets at decadal periods are not caused by the atmospheric and/or oceanic excitations. The existence of FCN in Figure 2 is natural since it is a free motion that requires only a very weak excitation. To consider this, we recalculated the excitations by applying Brzezinski broad-band Liouville equations, transformed into the frequency domain (Brzezinski, 1994; Brzezinski et al., 2002), expressing the transfer function of the matter (with subscription p) and motion (with subscription w) terms, respectively. It has the complex form

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1)

where [sigma'.sub.c] = 6.32000 + 0.00237i, [[sigma]'.sub.f] = -0.0146011 + 0.0001533i (rad/day) are the complex Chandler and FCN frequencies in celestial frame, respectively, whose imaginary parts are closely related to the corresponding quality factors (these values are based on observed polar motion and celestial pole offsets), [[sigma].sub.C] = [a'.sub.C] - [OMEGA] is the Chandler frequency in terrestrial frame, where [OMEGA] = 6.30038 rad/day is the angular speed of Earth's rotation, [chi]'p,[chi]'w are excitations (matter and motion terms) in celestial frame, and [a.sub.p] = 9.509 x [10.sup.-2], [a.sub.w] = 5.489 x [10.sup.-4] are dimensionless numerical constants, expressing the response to a matter and motion excitation, respectively. Formula (1) has two resonant frequencies--near Chandler and FCN frequency. The first one is unimportant in our case, since it is far from nutation frequencies, but the FCN resonance can amplify atmospheric and oceanic excitations substantially. The resulting spectra for the atmospheric and oceanic excitations, added together, convoluted with the transfer function (1), and re-calculated into mas, are shown in Figures 7 (NCEP + ECCO) and 8 (ERA + OMCT). The matter and motion terms are shown separately to demonstrate the dominancy of the matter terms. Comparison of both graphs reveals that the amplitudes are significantly higher in the latter case, due to a larger oceanic OMCT contribution.

[FIGURE 1 OMITTED]

[FIGURE 2 OMITTED]

[FIGURE 3 OMITTED]

[FIGURE 4 OMITTED]

[FIGURE 5 OMITTED]

[FIGURE 6 OMITTED]

We can use Brzezinski transfer function (1) also to estimate the atmospheric and oceanic contribution to individual nutation terms, by estimating the most significant terms of the excitation by the least-squares method and then multiplying them by T([sigma]). We can only argue which arguments to take for the dominant annual and semi-annual terms. In principle, the atmospheric and oceanic excitation should have periods connected with changing seasons, i.e., with the tropical year. On the other hand, for practical reasons, they should be identified with the known dominant station terms, in order to facilitace their mutual combination. Here we take the nutation arguments l' (with the period 365.26d), and 2F - 2D + 2[OMEGA] (with the period 182.62d) that are both very close to the tropical year and its half. To convert the estimated cosine/sine terms in X, Y into prograde/retrograde terms, we use the formula

[C.sup.+] = 1 / 2 [[S.sub.x] - [C.sub.y] + i([C.sub.x] + [S.sub.y])],

[C.sup.-] = 1 / 2 [[-S.sub.x] - [C.sub.y] + i([C.sub.x] - [S.sub.y])] (2)

The results, based on data span 1993.0-2009.5 (NCEP + ECCO) and 1979.0-2009.0 (ERA + OMCT), are shown in Table 1.

Another possibility of how to compare the excitations with the observed celestial pole offsets is to numerically integrate the Brzezinski broad-band Liouville equations. In complex form they read as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3)

where P = dX + idY is the excited motion of Earth's spin axis in celestial frame.

[FIGURE 7 OMITTED]

The numerical integration of Eq. (3) is done by the fourth-order Runge-Kutta method with 6-hour step. We use the procedure rk4 from Numerical Recipes (Press et al., 1992) that we modified to our purpose by rewriting it into the complex form. To obtain two first-order equations, instead of a second-order one given by Eq. (3), we use the substitution [y.sub.1] = P, [y.sub.2] = [??] - i[[sigma]'.sub.c]P, which leads to differential equations for two complex functions [y.sub.1]; [y.sub.2]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4)

The solution generally yields two free damped circular motions: a rapid prograde Chandler wobble and a slow retrograde FCN with celestial frequencies [[sigma]'.sub.c] and [[sigma]'.sub.f,] respectively. To integrate the system (4), we need to choose the initial values, that is two complex constants, defining the amplitudes and phases of both free motions. We are not interested in rapid (nearly diurnal) Chandlerian motion, so we choose only one, pole position at initial epoch [P.sub.0]; its first derivative is constrained so that the Chandlerian amplitude disappears. This is easily assured by choosing the values [y.sub.1](0) = [P.sub.0] [y.sub.2](0) = i([[sigma]'.sub.f] - [a'.sub.c])[P.sub.0]. The final choice of [P.sub.0] is made by repeating the integration with different values [P.sub.0] until the fit of the integrated motion to VLBI observations reaches a minimum. The results of these integrations are displayed in Figure 9 (for NCEP + ECCO) and 10 (for ERA + OMCT) as full black lines, in both cases shown together with VLBI-based celestial pole offsets (gray points). Even if we have the data in a longer interval for ERA + OMCT, we use the interval of only 1993.0-2009.0 to facilitate the comparison of both integrations. Annual and semi-annual parts were removed from both series, and also long-periodic component from IVS, so that they represent practically only the FCN. Both figures show that the integration follows more or less the same pattern as the observations; the rms fit is slightly better for NCEP + ECCO. The amplitudes seem to be too high for ERA + OMCT. The integrated series can be used to estimate the annual and semi-annual sine/cosine terms. Using formula (2) we arrive again to prograde and retrograde nutations that are displayed in Table 2. It is necessary to say that the whole intervals common for IVS and geophysical fluids are used here, i.e., 1993.0-2009.5 for NCEP + ECCO and 1984.1-2009.0 for ERA+OMCT. Theoretically, the results should be identical with those shown in Table 1 (rows in bold). Table 1 also contains, in the last row, the estimated nutation corrections from IVS solution, in which we also considered that the IAU2000 model of nutation implicitly contains empirical, so called Sun-synchronous correction to prograde annual term, that is equal to (-0.0104 + 0.1082i) mas. This value was added to the term estimated from IVS, to be directly comparable with the atmospheric contribution.

[FIGURE 8 OMITTED]

The amplitudes are mostly similar, but there are larger differences in some cases, e.g., for retrograde annual term or prograde semi-annual term.

4. RESONANCES IN EARTH ROTATION AND ESTIMATION OF P, Q

The resonances in Earth orientation are most completely described by Mathews et al. (2002). They derived so called Mathews-Herring-Buffet (MHB) transfer function in the form

T([sigma]) = [e.sub.R] - [sigma] / [e.sub.R] + 1 [N.sub.0] [1 + (1 + [sigma]) ([Q.sub.0] + [[summation].sup.4.sub.j=1] [Q.sub.j] / [sigma] - [S.sub.j])], (5)

which expresses the ratio of the non-rigid amplitude of a forced nutation term with terrestrial frequency [sigma] (in cycles per sidereal day--cpsd) to its rigid Earth value. Here [e.sub.R] denotes the dynamical ellipticity of the rigid Earth used to compute the rigid solution, [N.sub.0], [Q.sub.j] are complex "strength" parameters, and [S.sub.j] are four complex resonance frequencies corresponding to Chandler Wobble (CW, with terrestrial period of about 435 days), retrograde FCN (with celestial period of about 430 days), Prograde Free Core Nutation (PFCN, with celestial period of about 1020 days) and Inner Core Wobble (ICW, with terrestrial period of about 2400 days), respectively. In our study, only [s.sub.2] [approximately equal to] -1.0023 cpsd (FCN frequency) is important since it is close to the frequencies of all nutation terms and, at the same time, the corresponding coefficient [Q.sub.2] (~4.89 x [10.sup.-2]) is two orders of magnitude larger than [Q.sub.3] (~2.96 x [10.sup.-4]), corresponding to PFCN. Equation (5) was used, together with a rigid-Earth solution by Souchay et al. (1999) to derive the presently adopted model of nutation IAU2000.

MHB transfer function (5) can be used to derive the parameters on its right-hand-side, provided the value T([sigma]) is known for several different frequencies [sigma] from the observations (for more details see, e.g., Vondrak et al., 2005). Here we use the amplitudes and phases of five dominant nutation terms (with periods 365.26, 182.62, 121.75, 27.55 and 13.66 days) as determined by IVS, both for positive and negative frequencies. We estimate them in six-year running intervals, using the weights computed from VLBI formal standard errors. To account for possible longer term variations, we insert also a bias and a linear trend to the least-squares fit, and, to avoid possible aliasing, we also estimated a term with FCN period. Dividing the estimated nutation amplitudes by their rigid-Earth values by Souchay et al. (1999) we get the values T. It is necessary to note that MHB solution also contains a relatively simple model of atmospheric excitation in the form of so called Sun-synchronous correction (see also above). It is a prograde annual term with the amplitude of about 100 [micro]as, obtained as an empirical correction removing obvious residuals at this frequency. We removed this correction from VLBI-based celestial pole offsets before determining the values T. Thus, we obtain ten different complex values of T (i.e., twenty different "observables"). These, in turn, were used to estimate the complex frequency [s.sub.2] in a weighted least-squares solution, in six-year moving intervals; all remaining parameters of Eq. (5) were fixed to MHB values. Once we have the complex values of resonance frequency [s.sub.2], we can compute the period and quality factor from the simple equations

P = 0.99727 / [Re([s.sub.2]) + 1], Q = -Re([s.sub.2]) / 2Im([s.sub.2]) (6)

where P is given in solar days.

[FIGURE 9 OMITTED]

[FIGURE 10 OMITTED]

The results for P and Q values are graphically depicted in Figures 11 and 12, together with their error bars. The uncertainties are computed from the dispersion among different values of T for different nutation frequencies. Full lines in both figures represent the variations based on only IVS observations and an assumption that the geophysical excitation is modeled by MHB Sun-synchronous correction. One question automatically arises: Can the observed variations of P and Q be caused by the influence of the atmosphere and oceans that can significantly perturb estimates of the forced nutation amplitudes? To verify this possibility, we use the time series of that part of nutation that is due to geophysical excitations, obtained in preceding section. We remove these integrated values from VLBI observations, instead of MHB Sun-synchronous correction, and repeated the same estimation of P, Q in six-year running intervals. The plots of these results are added to Figures 11 and 12, from which we see that using the real atmospheric and oceanic excitations instead of MHB model does not improve the temporal stability of the results. Different models of atmosphere/ocean yield results that differ significantly, and also error bars are larger. From this we also see how much the atmospheric and oceanic excitations at near-diurnal band from different sources diverge. It is in agreement with the findings by Yseboodt at al. (2002) that different agencies providing atmospheric excitations differ significantly at near-diurnal band. In addition to this, the average values of P and Q seem to be systematically different if real geophysical excitations are used instead of MHB Sun-synchronous correction.

5. CONCLUSIONS

Forced nutations due to excitation by the atmosphere and ocean are significant, especially at annual and semi-annual periods, but different models give different results. Slightly better agreement with VLBI-based celestial pole offsets is obtained for atmospheric and oceanic excitation NCEP + ECCO, if compared with ERA + OMCT (see Figs. 9 and 10). The quasi-periodic variations of estimated period P and quality factor Q of Free Core Nutation keep within the limit of 429.8-430.5 solar days in case in the period, and 17000-22000 in case of the quality factor. The observed variations exceed the formal uncertainty of the results, so they seem to be real. If the geophysically excited nutations are used instead of MHB Sun-synchronous correction in P and Q estimation, these variations are even larger. A systematically longer period of FCN is obtained for both models used, and systematic growth of quality factor is observed for ERA + OMCT excitations. The temporal variations of P, Q are thus heavily dependent on the removed AAM/OAM effect from different sources. We conclude that the current AAM and OAM data, at the present state of their modeling, are not capable of explaining the observed variations of P and Q. Other processes (as. e.g., fluctuations of the core flattening caused by mantle convection at core-mantle boundary) are probably responsible for this effect.

[FIGURE 11 OMITTED]

[FIGURE 12 OMITTED]

ACKNOWLEDGMENTS

This work was supported by the project LC506 "Recent dynamics of the Earth," financed by the Ministry of Education, Youth and Sports of the Czech Republic, and also by the grant No. 205/08/0908 "Analysis of the geophysical excitations in the Earth orientation parameters," awarded by the Grant Agency of the Czech Republic. The help of the referees, Z. Martinec and J. Hefty, is also appreciated.

REFERENCES

Barnes, R.T., Hide, R., White, A.A., and Wilson, C.A.: 1983, Atmospheric angular momentum fluctuations, length-of-day changes and polar motion. Proc. R. Soc. London A387, 31-73.

Bizouard, C.: 1999, The atmospheric contribution to precession-nutation. In: N. Capitaine (ed.) Journees 1998 Systemes de Reference Spatio-temporels, Obs. de Paris, 92-98.

Bizouard, C., Brzezinski, A., and Petrov, S.: 1998, Diurnal atmospheric forcing and temporal variations of the nutation amplitudes. J. Geod. 72, 561-577.

Brzezinski, A.: 1994, Polar motion excitation by variations of the effective angular momentum function: II. Extended model. Manuscripta Geodaetica. 19, 157-171.

Brzezinski, A., Bizouard, C., and Petrov, S.: 2002, Influence of the atmosphere on Earth rotation: What can be learned from the recent atmospheric angular momentum estimates? Surveys in Geophysics 23, 33-69.

Capitaine, N., Wallace, P., and Chapront, J.: 2003, Expressions for IAU 2000 precession quantities. Astron. Astrophys. 412, 567-586.

Dobslaw, H. and Thomas, M.: 2007, Simulation and observation of global ocean mass anomalies. J. Geophys. Res. 112, C05040, doi: 10.1029/2006JC004035.

Gross, R. S., Fukumori, I., and Menemenlis, D.: 2005, Atmospheric and oceanic excitation of decadal-scale Earth orientation variations, J. Geophys. Res., 110, B09405, doi:10.1029/2004JB003565.

Gwinn, C.R., Herring, T.A., and Shapiro, I.I: 1986, Geodesy by radio interferometry: Studies of the forced nutations of the Earth 2. Interpretation. J. Geophys. Res. 91, B5, 4755-4764.

Herring, T.A., Gwinn, C.R., and Shapiro, I.I: 1986, Geodesy by radio interferometry: Studies of the forced nutations of the Earth 1. Data analysis. J. Geophys. Res. 91, B5, 4745-4754.

Hinderer, J., Boy, J., Gegout, P., et al.: 2000, Are the free core nutation parameters variable in time? Physics of the Earth and Planetary Interiors 117, 37-49.

Lambert, S.B.: 2006, Atmospheric excitation of the Earth's free core nutation. Astron. Astrophys. 457, 717-720.

Lambert, S.B. and Dehant, V.: 2007, The Earth's core parameters as seen by VLBI. Astron. Astrophys. 469, 777-781.

Mathews, P.M., Herring, T.A., and Buffet, B.A.: 2002, Modeling of nutation and precession: New nutation series for nonrigid Earth and insights into the Earth's interior. J. Geophys. Res. 107, doi: 10.1029/2001JB000390.

Press, W.H., Teukolsky, S.A., Vetterling, W.T., and Flannery, B.P.: 1992, Numerical Recipes in Fortran, The art of scientific computing, 2nd edition, Cambridge University Press.

Roosbeek, F., Defraigne, P., Feissel, M., and Dehant, V.: 1999, The free core nutation period stays between 431 and 434 sidereal days. Geophys. Res. Lett. 26, 131-134.

Salstein, D.: 2005, Computing atmospheric excitation functions for Earth rotation/polar motion. Cahiers du Centre Europeen de Geodynamique et de Seismologie 24, Luxembourg, 83-88.

Schluter, W. and Behrend, D.: 2007, The International VLBI Service for Geodesy and Astrometry (IVS): Current capabilities and future prospects. J. Geod. 81, 379-387.

Souchay, J., Loysel, B., Kinoshita, H., and Folgueira, M.: 1999, Corrections and new development in rigid Earth nutation theory III. Final tables "REN-2000" including crossed-nutation and spin-orbit coupling effects. Astron. Astrophys. Suppl. Ser. 135, 111-131.

Thomas, M., Dobslaw, H., and Soffel, M.: 2006, The ocean's response to solar thermal and gravitational tides and impacts on EOP. In: A. Brzezinski, N. Capitaine, B. Kolaczek (eds.), Proc. Journees 2005 Systemes de reference spatio-temporels, Space Research Centre Warsaw and Observatoire de Paris, 203-206.

Vondrak, J.: 1977, Problem of smoothing observational data II, Bull. Astron. Inst. Czechosl. 28, 84-89.

Vondrak, J.: 2009, Geophysical contributions in precession-nutation. In: M.S. Dimitrijevic, M. Tsvetkov, L.C. Popovic, V. Golev (eds.) Proc. VI Serbian-Bulgarian Astronomical Conference, Publ. Astr. Soc. Rudjer Boskovic No. 9, 143-153.

Vondrak, J. and Ron, C.: 2005, Combining GPS and VLBI measurements of celestial motion of the Earth's spin axis and Universal time. Acta Geodyn. Geomater., Vol.2, No. 3, 87-94.

Vondrak, J. and Ron, C.: 2006a, Resonant period of free core nutation--its observed changes and excitations. Acta Geodyn. Geomater. 143, 53-60.

Vondrak, J. and Ron, C.: 2006b, Resonance effects and possible excitation of Free Core Nutation. In: A. Brzezinski, N. Capitaine, B. Kolaczek (eds.) Journees 2005 Systemes de Reference Spatio-temporels, Space Research Centre Warsaw and Obs. de Paris, 215-218.

Vondrak, J. and Ron, C.: 2007, Quasi-diurnal atmospheric and oceanic excitation of nutation. Acta Geodyn. Geomater. 148, 121-128.

Vondrak, J. and Ron, C.: 2008, VLBI observations of nutation, its geophysical excitations and determination of some Earth model parameters. In: N. Capitaine (ed.) Journees 2007 Systemes de reference spatio-temporels, Observatoire de Paris, 95-98.

Vondrak, J. and Ron, C.: 2009, Stability of period and quality factor of free core nutation. Acta Geodyn. Geomater. 155, 217-224.

Vondrak, J., Weber, R,. and Ron, C.: 2005, Free Core Nutation: direct observations and resonance effects. Astron. Astrophys. 444, 297-303.

Wahr, J.: 1981, The forced nutations of an elliptical, rotating, elastic and oceanless Earth. Geophys. J. R. Astron. Soc. 64, 705-727.

Yseboodt, M., de Viron, O., Cin, T.M., and Dehant, V.: 2002, Atmospheric excitation of the Earth's nutation: Comparison of different atmospheric models. J. Geophys. Res. 107, doi: 10.1029/2000JB000042.

Jan VONDRAK * and Cyril RON

Department of Galaxies and Planetary Systems, Astronomical Institute, Academy of Science of the Czech Republic, Bocni II, 141 31 Prague, Czech Republic, phone +420 267103043, fax +420 2 72769023

* Corresponding author's e-mail: vondrak@ig.cas.cz

(Received December 2009, accepted February 2010)

Table 1 Atmospheric and oceanic contribution in nutation [[micro]as], calculated by convolution with Brzezinski transfer function (1). Annual (/') Semi-annual (2F-2 D+2[OMEGA]) Excitation prograde retrograde prograde Re Im Re Im Re Im NCEP+ECCO p -62.0 58.1 -69.9 -2.7 -6.5 21.4 w 3.5 29.8 -8.5 -17.6 2.5 51.0 sum -58.5 87.9 -78.4 -20.3 -4.0 72.4 ERA+OMCT p -141.2 98.3 -46.8 63.1 -47.2 17.4 w 10.9 30.9 -6.1 -6.6 7.6 54.5 sum -130.3 129.2 -52.9 56.5 -39.6 71.9 Semi-annual (2F-2 D+2[OMEGA]) Excitation retrograde Re Im NCEP+ECCO p 11.0 13.1 w 2.1 -1.3 sum 13.1 11.8 ERA+OMCT p 3.4 -5.8 w 1.6 -4.1 sum 5.0 -9.9 Table 2 Atmospheric and oceanic contribution in nutation [[micro]as], computed from numerically integrated broadband Liouville equations (3), and comparison with IVS-based values. Annual (l') Semi-annual (2F-2D+2[OMEGA]) excitation prograde retrograde prograde retrograde Re Im Re Im Re Im Re Im NCEP+ECCO -58.6 87.8 -84.6 -21.6 -4.3 71.2 10.3 12.7 ERA+OMCT -126.7 97.6 -39.0 43.7 -43.3 12.7 -4.0 -6.7 IVS -22.5 104.7 -17.6 2.9 -0.6 19.8 -1.3 -16.9

Printer friendly Cite/link Email Feedback | |

Author: | Vondrak, Jan; Ron, Cyril |
---|---|

Publication: | Acta Geodynamica et Geromaterialia |

Article Type: | Report |

Geographic Code: | 4EXCZ |

Date: | Jan 1, 2010 |

Words: | 4738 |

Previous Article: | Precise near real-time GNSS analyses at Geodetic Observatory Pecny--precise orbit determination and water vapour monitoring. |

Next Article: | Earth orientation parameters and station coordinates from space geodesy techniques. |

Topics: |