Issue 
A&A
Volume 526, February 2011



Article Number  A106  
Number of page(s)  4  
Section  Celestial mechanics and astrometry  
DOI  https://doi.org/10.1051/00046361/201015894  
Published online  06 January 2011 
The Earth’s variable Chandler wobble
Observatoire de Paris, Département Systèmes de Référence Temps Espace
(SYRTE), CNRS/UMR 8630, 75014
Paris, France
email: christian.bizouard@obspm.fr
Received:
8
October
2010
Accepted:
17
November
2010
Aims. We investigated the causes of the Earth’s Chandler wobble variability over the past 60 years. Our approach is based on integrating of the atmospheric and oceanic angular momentum computed by global circulation models. We directly compared the result of the integration with the Earth’s pole coordinate observed by precise astrometric, space, and geodetic techniques. This approach differs from the traditional approach in which the observed polar motion is transformed into a socalled geodetic excitation function, and compared afterwards with the angular momentum of the external geophysical fluid layers.
Methods. In the time domain, we integrated the atmospheric angular momentum time series from the National Center for Environmental Prediction/National Center for Atmospheric Research Reanalysis project and the oceanic angular momentum data from the ECCO consortium. We extracted the Chandler wobble from this modeled polar motion by singular spectrum analysis, and compared it with the Chandler wobble extracted from the observed polar motion given by the International Earth Rotation and Reference Systems Service data.
Results. We showed that the combination of the atmosphere and the oceans explains most of the observed Chandler wobble variations, and is consistent with results reported in the literature and obtained with the traditional approach. Our approach allows one to appreciate the separate contributions of the atmosphere and the oceans to the various bumps and valleys observed in the Chandler wobble. Though the atmosphere explains the Chandler wobble amplitude variations between 1949 and 1970, the reexcitation of the Chandler wobble that begins in the 1980s, after a minimum around 1970, and that reaches its maximum in the late 1990s is due to the oceans, while the atmospheric contribution remains stable within the same period.
Key words: reference systems / Earth
© ESO, 2011
1. Introduction
The Chandler wobble (CW) is a free rotational mode of the Earth associated with its dynamical ellipticity. In the absence of forcing, the CW would have a period of 430.3 days (Mathews et al. 2002) fixed mainly by the whole planet’s dynamical ellipticity and, more marginally, by mantle and ocean deformabilities, and would lose most of its energy after a few tens of years because of dissipation in the mantle and in the oceans. Observation of the Earth’s polar motion, however, reveals a prograde oscillation whose pseudo period can be as far as 20 days from the above value (e.g., Vondrák 1988; Vicente & Wilson 1997). It gains energy at some times (e.g., Danjon & Guinot 1954) so that it never disappears.
Plag (1997) interpreted the CW as a “forced free” mode and proposed a climatic forcing. He associated the temporal characteristics of the fourteen to sixteen months oscillation over the northern hemisphere with frequency or phase shifts of the CW. According to Celaya et al. (1999), Gross (2000), and Brzeziński & Nastula (2002), the CW excitation is accounted for, on average, by the atmosphere and oceans. Following Aoyama et al. (2003), atmospheric and oceanbottom pressure and winds are comparable sources of excitation, intermittently playing a prominent role. Seitz & Schmidt (2005) conclude that the atmospheric and oceanic bottom pressure fluctuations are the prominent cause of the CW variability. They could account for some features in the band of 400–460 days by wavelet analysis: an amplitude increase between 1980 and 1990 and a decrease after 1993. However, they did not precisely investigate the amplitude and phase variations of the Chandler period and neglected the wind contribution. Liao et al. (2007) noticed the variability in the wind contribution around the Chandler period. The lack of clear conclusions is partly due to discrepancies up to 0.7 millisecond of arc (mas) between the various atmospheric angular momentum data sets in the Chandler frequency band associated with inaccurately determining the wind term.
Though most of the above studies implement comparisons between the fluid layer angular momentum and the socalled geodetic excitation function derived from the observed pole coordinates, very few studies used the inverse approach that integrates the former and directly compares it with the observed polar motion (see, e.g., Vondrák 1989, 1990). We revisited the CW variability in light of the latter approach, as described in Sect. 2. Data are presented in Sect. 3 and results discussed in Sect. 4.
2. The polar motion excitation by fluid layers
At low frequency and for a deformable Earth surrounded by fluid layers, the time evolution of the complex pole coordinates p = x − iy is given by (e.g., Barnes et al. 1983) (1)where σ_{c} = σ_{0}(1 + i/2Q) is the Chandler frequency, c = c_{13} + ic_{23} and h = h_{1} + ih_{2} are the offdiagonal moment of inertia and the equatorial relative angular momentum of the fluid layer, respectively, and A and C are the mean equatorial and axial moments of inertia of the Earth, respectively. The LHS of (1) is referred to as the geodetic excitation χ_{g}, while χ_{f} is referred to as the geophysical (or fluid) excitation. This equation does not account for the existence of fluid outer and inner cores. However, the resonances produced by these internal layers lie in the diurnal frequency band and do not substantially affect seasonal and interannual polar motion. Mathews et al. (2002) predicts a resonance due to the inner core with a theoretical period of ~2400 days but with an amplitude expected to be insignificant for our purpose.
The resonant frequency σ_{c}, or eigenfrequency, can be expressed in terms of the whole Earth’s dynamical ellipticity e = (C − A)/C and the Love number k, with additional contributions from anelasticity, and ocean loading and currents to the mantle deformability (see, e.g., Sasao et al. 1980; Mathews et al. 1991, 2002). The theoretical computation of σ_{c} leads to a period of 430.3 days and a quality factor Q = 88 (Mathews et al. 2002). The theoretical period is consistent with other estimates of the observed Chandler period based on analyses of polar motion time series (e.g., Vicente & Wilson 1997). However, there is no consensus on the value of Q: depending on the time interval and the analysis method, it ranges from 50 (Furuya & Chao 1996; Kuehne et al. 1996) to 180 (Vicente & Wilson 1997; Aoyama et al. 2003). Gross (2000) and Brzeziński & Nastula (2002) find intermediate values of ~140.
The most common approach to explaining the Chandler wobble excitation by external fluid layers consists of computing the LHS of (1) from observed pole coordinates using, e.g., the discrete formula proposed by Wilson (1985) (also Eq. (1) of Vicente & Wilson 1997), and comparing it against the RHS directly given by fluid layer angular momentum time series. Nevertheless, the noise behavior of χ_{f} and χ_{g} around the Chandler frequency and their weak amplitude compared with the seasonal excitation can make this approach difficult. Conversely, a modeled polar motion obtained by time integration of the RHS of (1) can be directly compared with the observed pole coordinates.
The solution of (1) reads (2)wherein t_{0} is the time at which one starts the integration. The first term represents a damped, free wobble of relaxation time 2Q/σ_{0} ~ 40 years for Q = 100. This progressive damping is not observed, but balanced by the integral expressing the forced polar motion. The latter term mainly contains the effects of the seasonal forcing and an oscillation emerging in the Chandler period.
3. The data set
We used the pole coordinates from the IERS EOP C 01 data set, computed at the International Earth Rotation and Reference Systems Service (IERS) Earth Orientation Parameter Product Center^{1}. These series are the result of the combination and smoothing of polar motion series obtained from optical astrometry, very long baseline radio interferometry, and various space and geodetic techniques (mainly global navigation and satellite system). Using other Earth orientation parameter series produced by other institutes (e.g., by the IERS Rapid Service/Prediction Center of the United States Naval Observatory, Washington, DC) would lead to the same conclusions since the agreement between the various data sets reaches a few mas, which is far below the CW amplitude variations investigated in this paper.
We considered the geophysical forcing made up of (i) atmospheric angular momentum (AAM) as given by the National Center for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis project (Kalnay et al. 1996), and (ii) oceanic angular momentum (OAM) output from the ECCO model that employs the MIT global circulation model in a nearglobal domain (Fukumori et al. 2002). The model is forced by NCEP/NCAR reanalysis products (12hourly wind stresses, daily adiabatic airsea fluxes). In some versions of the ECCO run, the temperature and salinity of the model sea surface can be relaxed toward observed values, and sea level anomalies can be assimilated into the model to correct errors associated with inaccuracies in timevarying wind forcing. However, the longest available series over 1949–2002 (kf049f) were produced with no data assimilation. To complete the series up to 2009.5, we therefore concatenated with the series kf079, which is also obtained without data assimilation.
Fig. 1 The observed polar motion (top) and the Chandler mode (bottom) between 1949.0 and 2009.5 from the IERS EOP C 01 data set. 

Open with DEXTER 
The data sets are made available through the IERS Special Bureaus for the Atmosphere^{2} (Salstein et al. 1993) and the Oceans^{3} (chaired by Gross at JPL, Pasadena). Since we investigated seasonal and interannual time scales, the AAM was used under the inverted barometer hypothesis. All time series spanned 1949.0–2009.5. They were made to be homogeneous by removing the high frequency part of the signal below 10 days with a Vondrák (1977) filter of admittance larger than 99% above 10 days and smaller than 1% below 4 days, and resampling them at 10day intervals.
4. Results and discussion
Using (2), we computed the modeled polar motion p_{f}(t) from the geophysical excitation by a trapezoidal integration. The value of the observed CW at t_{0} was taken as initial condition. Doing so, we assumed that the atmosphere and the oceans are the only sources of the observed CW, thus neglecting the contribution of continental waters. However, based on a predictive model, Celaya et al. (1999) concluded that the hydrological forcing of the CW does not compete with atmosphere and oceans.
The CW component of this modeled polar motion had to be compared with the CW component of the observed, or geodetic, polar motion. The CW was extracted from the observed and modeled polar motions by singular spectrum analysis (SSA, Vautard & Ghil 1989), a technique that allows one to isolate quasi periodic modes without assumptions on their periods. We chose a covariance lag close to 6.4 years in order to separate the CW from the annual oscillation whose amplitude is of the same order of magnitude (see Fig. 1). To figure out the error propagated by the SSA decomposition, we ran the SSA on a large number of replicated pole coordinate time series perturbed by a Gaussian noise of the same variance as the residuals (i.e., signal cleaned up from the first five reconstructed components including CW, annual wobble, and long term irregular variation), allowing one to make a reliable statistics on the stability of the principal components. We found that the relative error is ~1%.
We repeated the integration of the geophysical data sets and the comparison with observed polar motion using different values for Q and computed the regression coefficient between the observed and modeled Chandler modes (Fig. 2). One sees that Q ~ 84 produces an optimal regression.
Fig. 2 The regression coefficient between the observed and modeled Chandler modes as a function of Q. 

Open with DEXTER 
Fig. 3 The observed and modeled variations in the amplitude and the phase of the Chandler mode. 

Open with DEXTER 
Fig. 4 The regression coefficient between the observed and modeled Chandler modes as a function of the starting epoch t_{0}. 

Open with DEXTER 
The amplitude and phase of the observed and modeled CW obtained with Q = 84 are displayed in Fig. 3. The regression coefficient between the observed and modeled CW is 1.001 + 0.027 i. The atmosphere alone poorly accounts for the CW irregularities. Though it explains the large amplitude before 1960, the minimum around 1970 and the bump in 1980−1990 is mainly due to the oceans, which confirms the results of Gross (2000). The atmospheric contribution remains stable within the same period. Similar remarks can be made for the phase variations before 1995. Around 2000, the phase of the observed CW presents a bump around −70° that is likely due to the oceans. Nevertheless, the combination of atmosphere and oceans predicts a phase variation that is different by ~20°.
The phase shifts of the CW are equivalent to variations of the observed Chandler period that can reach 33 days. The observed Chandler period, which is the pseudo period of the “forced free” CW, is distinct from the eigenperiod defined by the Earth’s structure properties. It corresponds to the period for which the forcing exhibits maximal energy in the vicinity of the eigenperiod. The eigenperiod cannot change by such an amount on time scales of a few years.
We repeated the integration using different starting epochs t_{0}, and we computed the regression coefficient between the observed and modeled Chandler modes (Fig. 4). The combination of the atmosphere with the oceans generally drew the amplitude of the regression coefficient towards unity and the phase closer to zero.
As an alternative to the SSA technique, we extracted the Chandler mode from observed or modeled polar motion by leastsquare adjustments over a 6.4yr sliding window. This method demands, however, making an assumption about the period of the Chandler mode (fixed to 430.3 days). The results were very close to those presented above. Another alternative consisted of slicing the wavelet spectrum at the Chandler frequency. This method was used in Seitz & Schmidt (2005) over a shorter time period, and it assessed the pressure terms of the atmospheric and oceanic angular momenta only, yielding results that were consistent with ours over the common time period.
Fig. 5 The real and imaginary parts of the complex amplitude A(t) of the observed and modeled excitation functions at the Chandler frequency. 

Open with DEXTER 
To assess the consistency of the integration approach, we finally transformed the observed and modeled CW extracted by SSA into excitation functions using the Wilson (1985) formula. Let (3)be the modeled or observed CW, wherein the time variable complex amplitude a(t) is estimated by least squares over a 6.4yr sliding window. As previously, the Chandler period is fixed to 430.3 days. Possible, but small, variations in the period are considered in the phase of a(t). From Eq. (1), the corresponding excitation is given by (4)where (5)\newpageFigure 5 displays the real and imaginary parts of the complex amplitude A(t) deduced from observed and modeled CW. Since the formal error on a(t) is of the order of 1 mas, one expects the error on A(t) to be majored by 0.3 mas, which is far smaller than the overall variability of the signals, so the integration approach provides an original way of isolating the CW excitation.
Our results show that the integration of the NCEP/NCAR reanalysis project’s atmospheric data, combined with the ECCO oceanic data, explain most of the observed variations in the CW amplitude and phase since the middle of the twentieth century. The originality of this work resides in its using the integration approach, which complements the classical approach in most other works and is based on comparing excitation functions. The integration approach therefore constitutes an interesting alternative method, which has the advantage being more intuitive from the point of view of the astronomer who directly observes the position of the pole.
Acknowledgments
The authors are grateful to Drs. O. de Viron, L. Zotov, and Prof. V. Dehant for useful discussions and suggestions.
References
 Aoyama, Y., Naito, I., Iwabuchi, T., & Yamazaki, N. 2003, Earth Plan. Sp., 55, 25 [NASA ADS] [Google Scholar]
 Barnes, R. T. H., Hide, R., White, A. A., & Wilson, C. A. 1983, Proc. R. Soc. London A, 387, 31 [NASA ADS] [CrossRef] [Google Scholar]
 Brzeziński, A., & Nastula, J. 2002, Adv. Sp. Res., 30, 195 [NASA ADS] [CrossRef] [Google Scholar]
 Celaya, M., Wahr, J. M., & Bryan, F. O. 1999, J. Geophys. Res., 104, 12813 [NASA ADS] [CrossRef] [Google Scholar]
 Danjon, A., & Guinot, B. 1954, C. R. Acad. Sci., 238, 1081 [Google Scholar]
 Furuya, M., & Chao, B. F. 1996, Geophys. J. Int., 127, 693 [NASA ADS] [CrossRef] [Google Scholar]
 Fukumori, I., Lee, T., & Menemenlis, D. et al. 2002, International Symp. En Route to GODAE, Biarritz, France [Google Scholar]
 Gross, R. S. 2000, Geophys. Res. Lett., 27, 2329 [NASA ADS] [CrossRef] [Google Scholar]
 Kalnay, E., Kanamitsu, M., & Kistler, R. 1996, Bull. Amer. Meteorol. Soc., 77, 437 [NASA ADS] [CrossRef] [Google Scholar]
 Kuehne, J., Wilson, C. R., & Johnson, S. 1996, J. Geophys. Res., 101, 13573 [NASA ADS] [CrossRef] [Google Scholar]
 Liao, D., Zhou, Y., & Liao, X. 2007, Chinese Astron. Astrophys., 31, 57 [NASA ADS] [CrossRef] [Google Scholar]
 Mathews, P. M., Buffett, B. A., Herring, T. A., & Shapiro, I. I. 1991, J. Geophys. Res., 96, B5, 8243 [Google Scholar]
 Mathews, P. M., Herring, T. A., & Buffett, B. A. 2002, J. Geophys. Res., 107, B4, ETG 31 [NASA ADS] [CrossRef] [Google Scholar]
 Plag, H. P. 1997, In Tidal Phenomena, Lecture Notes in Earth Sciences, 66, (Berlin/Heidelberg: Springer), 183 [Google Scholar]
 Salstein, D. A., Kann, D. M., Miller, A. J., & Rosen, R. D. 1993, Bull. Amer. Meteor. Soc., 10, 67 [CrossRef] [Google Scholar]
 Sasao, T., Okubo, S., & Saito, M. 1980, In Proc. IAU Symp. (Hingham, Mass.: D. Reidel), ed. E. P. Federov, et al., 78, 165 [Google Scholar]
 Seitz, F., & Schmidt, M. 2005, J. Geophys. Res., 110, B11406 [NASA ADS] [CrossRef] [Google Scholar]
 Vautard, R., & Ghil, M. 1989, Phys. Rev. D, 35, 395 [Google Scholar]
 Vicente, R. O., & Wilson, C. R. 1997, J. Geophys. Res., 102, 20439 [NASA ADS] [CrossRef] [Google Scholar]
 Vondrák, J. 1977, Bull. Astron. Inst. Czech, 28, 84 [NASA ADS] [Google Scholar]
 Vondrák, J. 1988, In International Astronomical Union Symp. (Dordrecht: Kluwer Academic Publishers), ed. A. Kay, & G. A. Wilkins, 128, 359 [Google Scholar]
 Vondrák, J. 1989, Wiss. Z. Tech. Universit. Dresden, 38, 49 [Google Scholar]
 Vondrák, J. 1990, Bull. Astron. Inst. Czechosl., 49, 211 [NASA ADS] [Google Scholar]
 Wilson, C. R. 1985, Geophys. J. R. Astron. Soc., 80, 551 [Google Scholar]
 Wilson, C. R., & Chen, J. 1996, J. Geod., 70, 581 [NASA ADS] [Google Scholar]
All Figures
Fig. 1 The observed polar motion (top) and the Chandler mode (bottom) between 1949.0 and 2009.5 from the IERS EOP C 01 data set. 

Open with DEXTER  
In the text 
Fig. 2 The regression coefficient between the observed and modeled Chandler modes as a function of Q. 

Open with DEXTER  
In the text 
Fig. 3 The observed and modeled variations in the amplitude and the phase of the Chandler mode. 

Open with DEXTER  
In the text 
Fig. 4 The regression coefficient between the observed and modeled Chandler modes as a function of the starting epoch t_{0}. 

Open with DEXTER  
In the text 
Fig. 5 The real and imaginary parts of the complex amplitude A(t) of the observed and modeled excitation functions at the Chandler frequency. 

Open with DEXTER  
In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.