Parameter variability of the observed periodic oscillations of polar motion with smaller amplitudes

Compared to the Chandler and annual wobbles, the higher-frequency components of polar motion (PM) have substantially smaller amplitudes. Therefore, their study has had to wait until higher-quality time series with high temporal resolution, as measured by space geodetic techniques, have become available. Based on the combined Earth orientation series SPACE99 computed by the Jet Propulsion Laboratory (JPL) from 1976 to 2000 at daily intervals, the periodic PM terms, in particular at the quasi-biennial, 300-day, semi-Chandler, semi-annual, 4-month, 90-day, 2-month and 1.5-month periods, have been separated by band-pass filtering and it has been found that the persistence of oscillations becomes less with increasing frequency. In order to quantify and better describe the parameter variability of these PM components over time, the radii, direction angles and period lengths were computed from the periodic terms filtered out from the time series. The results clearly show the characteristics and time evolution of the periodic PM components. The largest elliptic oscillation is the semi-annual wobble with a maximum semi-major axis of up to 13 mas (milliarc seconds). The other wobbles are smaller. They have maximum semi-major axes of between 3 and 8 mas. If the oscillations have period lengths of 4 months and less, then they are elapsed not only progradly, but also retrogradly. AcknowledgementsThis paper was presented at the 27th General Assembly of the European Geophysical Society in Nice, France, 22–26 April 2002. Thanks go to Kevin Fleming for his linguistic advice. The author would also like to thank Barbara Koaczek for suggesting some valuable improvements.     

Excitation of non-atmospheric polar motion by the migration of the Pacific Warm Pool

Changes in the annual variation of the Earths polar motion are found to be largely caused by the variation of the atmospheric angular momentum (AAM). Recent simulation results of oceanic general circulation models further suggest global oceanic effects on the annual polar motion in addition to the atmosphere. In comparison with previous model studies of global oceanic effects, this research particularly singles out a large-scale ocean anomaly and investigates its effect on the annual polar motion, determined from satellite observations of the movement of the Western Pacific Warm Pool (WPWP). Although the scale of the warm pool is much smaller than that of the solid Earth, analysis of the non-atmospheric polar motion excitation has shown that the WPWP contributes non-negligibly to the annual polar motion. The analysis consists of over 30 years of WPWP data (1970–2000) and shows values of polar motion excitation for the x-component of (2.5 mas, –79°) and for the y-component of (0.6 mas, 173°). Comparison of this result with the total geodetic non-atmospheric polar motion excitation of (10.3 mas, 59°) for the x-component and (10.6 mas, 62°) for the y-component shows the significance of the WPWP. Changes in the Earths polar motion have attracted significant attention, not only because it is an important geodetic issue, but also because it has significant value as a global measure of variations within the hydrosphere, atmosphere, cryosphere, and solid Earth, and hence global changes.Tel: 86–21–64386191 Fax: 86–21–64384618Acknowledgments. The authors are grateful to Dr. R. Gross (JPL) and two anonymous reviewers for providing invaluable comments. They also thank Dr. J.L. Chen (CSR) for helpful discussions. Y. Zhou, D. Zheng and X. Liao were supported by the National Natural Science Foundation of China (10273018, 10133010) and Key Project of Chinese Academy of Sciences (KJCX2-SW-T1). X-H. Yan was supported by the National Aeronautics and Space Administration (NASA) through Grant NAG5–12745, and by the National Science Foundation (NSF) through the Presidential Faculty Fellow award to X-H. Yan (OCE-9453499). W.T. Liu was supported by the NASA Physical Oceanography Program.     

Modeling and forecast of the polar motion excitation functions for short-term polar motion prediction

Short-term forecast of the polar motion is considered by introducing a prediction model for the excitation function that drives the polar motion dynamics. The excitation function model consists of a slowly varying trend, periodic modes with annual and several sub-annual frequencies (down to the 13.6-day fortnightly tidal period), and a transient decay function with a time constant of 1.5 days. Each periodic mode is stochastically specified using a second-order auto-regression process, allowing its frequency, phase, and amplitude to vary in time within a statistical tolerance. The model is used to time-extrapolate the excitation function series, which is then used to generate a polar motion forecast dynamically. The skills of this forecast method are evaluated by comparison to the C-04 polar motion series. Over the lead-time horizon of four months, the proposed method has performed equally well to some of the state-of-art polar motion prediction methods, none of which specifically features forecasting of the excitation function. The annual mode in the ₂ component is energetically the most dominant periodicity. The modes with longer periods, annual and semi-annual in particular, are found to contribute more significantly to forecast accuracy than those with shorter periods.     

7个月周期极移的非线性动力机制

从滞弹性阻尼形变摄动造成CW频率调制假设出发，对CW的共振激发模型加上参数的时变调制，变成了参数共振模型。经正演计算发现，参数共振模型完全符合CW的实际，表明滞弹性阻尼形变摄动造成频率的3％调制，进一步使得CW振幅调制可达70％以上。这一参数共振模型是一个非线性动力系统，在非线性情况下，运动将发生分岔，即多解。     

A study of gravity variations caused by polar motion using superconducting gravimeter data from the GGP network

Long-term continuous gravity observations, recorded at five superconducting gravimeter (SG) stations in the Global Geodynamic Project (GGP) network, as well as data on orientation variations in the Earths rotation axis (i.e. polar motion), have been used to investigate the characteristics of gravity variations on the Earths surface caused by polar motion. All the SG gravity data sets were pre-processed using identical techniques to remove the luni-solar gravity tides, the long-term trends of the instrumental drift, and the effects of atmospheric pressure. The analysis indicates that the spectral peaks, related to the Chandler and annual wobbles, were identified in both the power and product spectral density estimates. The magnitude of gravity variations, as well as the gravimetric amplitude factor associated with the Chandler wobble, changed significantly at different SG stations and during different observation periods. However, when all the SG observations at these five sites were combined, the gravimetric parameters of the Chandler wobble were retrieved accurately: 1.1613 ± 0.0737 for the amplitude factor and –1°.30 ± 1°.33 for the phase difference. The value of the estimated amplitude factor is in agreement with that predicted theoretically for the zonal tides of an elastic Earth model.     

New Solution for the Earth’s Free Wobble and Its Geophysical Implications

In this paper, the theory of the free wobble of the triaxial Earth is developed and new conclusions are drawn: the Euler period should be actually expressed by the first kind of complete elliptic integral; the trace of the free polar motion is elliptic and the orientations of its semi-minor and major axes are approximately parallel to the Earth’s principal axes A and B, respectively. In addition, the present theory shows that there is a mechanism of frequency-amplitude modulation in the Chandler wobble, whi...     

Analysis of decade-long time series of GPS-based polar motion estimates at 15-min temporal resolution

We present results from the generation of 10-year-long continuous time series of the Earth’s polar motion at 15-min temporal resolution using Global Positioning System ground data. From our results, we infer an overall noise level in our high-rate polar motion time series of 60 \(\upmu \hbox {as}\) (RMS). However, a spectral decomposition of our estimates indicates a noise floor of 4 \(\upmu \hbox {as}\) at periods shorter than 2 days, which enables recovery of diurnal and semidiurnal tidally induced polar motion. We deliberately place no constraints on retrograde diurnal polar motion despite its inherent ambiguity with long-period nutation. With this approach, we are able to resolve damped manifestations of the effects of the diurnal ocean tides on retrograde polar motion. As such, our approach is at least capable of discriminating between a historical background nutation model that excludes the effects of the diurnal ocean tides and modern models that include those effects. To assess the quality of our polar motion solution outside of the retrograde diurnal frequency band, we focus on its capability to recover tidally driven and non-tidal variations manifesting at the ultra-rapid (intra-daily) and rapid (characterized by periods ranging from 2 to 20 days) periods. We find that our best estimates of diurnal and semidiurnal tidally induced polar motion result from an approach that adopts, at the observation level, a reasonable background model of these effects. We also demonstrate that our high-rate polar motion estimates yield similar results to daily-resolved polar motion estimates, and therefore do not compromise the ability to resolve polar motion at periods of 2–20 days.     

Prediction of polar motion

Based on an analysis of polar motion behavior, we found the possibility of predicting polar motion up to one year in advance. Comparing these predicted polar coordinates with the observed ones (smoothed), the rms of the differences is about 0".02. The differences of the relative polar motion are much smaller. For any time interval of 20–30 days throughout the whole year, the rms of the relative polar motion differences is about 0".01. It appears that 80–90% of the polar motion is composed of the stable, predictable Chandler and annual terms.     

Linear drift and periodic variations observed in long time series of polar motion

 Two long time series were analysed: the C01 series of the International Earth Rotation Service and the pole series obtained by re-analysis of the classical astronomical observations using the HIPPARCOS reference frame. The linear drift of the pole was determined to be 3.31 ± 0.05 milliarcseconds/year towards 76.1 ± 0.80° west longitude. For the least-squares fit the a priori correlations between simultaneous pole coordinates x _p , y _p were taken into account, and the weighting function was calculated by estimating empirical variance components. The decadal variations of the pole path were investigated by Fourier and wavelet analysis. Using sliding windows, the periods and amplitudes of the Chandler wobble and annual wobble were determined. Typical periods in the variable Chandler wobble and annual wobble parameters were obtained from wavelet analyses. Received: 21 January 2000 / Accepted: 28 August 2000     

Possible improvement of Earth orientation forecast using autocovariance prediction procedures

 Autocovariance prediction has been applied to attempt to improve polar motion and UT1-UTC predictions. The predicted polar motion is the sum of the least-squares extrapolation model based on the Chandler circle, annual and semiannual ellipses, and a bias fit to the past 3 years of observations and the autocovariance prediction of these extrapolation residuals computed after subtraction of this model from pole coordinate data. This prediction method has been applied also to the UT1-UTC data, from which all known predictable effects were removed, but the prediction error has not been reduced with respect to the error of the current prediction model. However, the results show the possibility of decreasing polar motion prediction errors by about 50 for different prediction lengths from 50 to 200 days with respect to the errors of the current prediction model. Because of irregular variations in polar motion and UT1-UTC, the accuracy of the autocovariance prediction does depend on the epoch of the prediction. To explain irregular variations in x, y pole coordinate data, time-variable spectra of the equatorial components of the effective atmospheric angular momentum, determined by the National Center for Environmental Prediction, were computed. These time-variable spectra maxima for oscillations with periods of 100–140 days, which occurred in 1985, 1988, and 1990 could be responsible for excitation of the irregular short-period variations in pole coordinate data. Additionally, time-variable coherence between geodetic and atmospheric excitation function was computed, and the coherence maxima coincide also with the greatest irregular variations in polar motion extrapolation residuals. Received: 22 October 1996 / Accepted: 16 September 1997     

用极移数据直接测算钱德勒摆动Q值

运用已经发展起来的共振激发和参数共振模型,在频率调制的基础上,对IERS极移数据测算出钱德勒摆动的Q值,得到平均值为63,并以每百年0.8的速率在增大。这一结果与Q的滞弹性PREM模型理论值69十分接近,表明共振激发模型和参数共振模型完全与钱德勒摆动滞弹理论相容。     

Comparison of polar motion with oceanic and atmospheric angular momentum time series for 2-day to Chandler periods

During a 4-year period starting in July 1996 and using intervals ranging from 3 days to 4 years, four precise polar motion (PM) series have been compared to excitation by atmospheric angular momentum (AAM) augmented with oceanic angular momentum (OAM) data. The first three series (C03, C04 and Bulletin A) are multi-technique combinations generated by the International Earth Rotation and Reference Systems Service (IERS) and the fourth combined series (IGS00P02) is produced by the International GPS Service (IGS) using only GPS data. The IGS PM compared the best with the combined excitations of atmosphere and oceans (AAM+OAM) at all intervals, showing high overall correlation of 0.8–0.9. Even for the interval of only three days, the IGS PM gave a significant correlation of about 0.6. Moreover, during the interval of February 1999 – July 2000, which should be representative of the current precision of the IGS PM, a significant correlation (>0.4) extended to periods as short as 2.2 days and 2.5 days for the x_p and y_p PM components, respectively. When using the IERS Bulletin B (C04) PM and an interval of almost 6 years, starting in November 1994, the combined OAM+AAM accounted for practically all the annual, semi-annual and Chandler wobble (CW) PM signals. When only AAM was used, either the US National Centers for Environment Prediction reanalysis data, which were used throughout this study, or the Japanese Meteorological Agency data, two large and well-resolved amplitude peaks of about 0.1 mas/day, remained at the retrograde annual and CW periods.     

An improved empirical model for the effect of long-period ocean tides on polar motion

Because the tide-raising potential is symmetric about the Earth’s polar axis it can excite polar motion only by acting upon non-axisymmetric features of the Earth like the oceans. In fact, after removing atmospheric and non-tidal oceanic effects, polar motion excitation observations show a strong fortnightly tidal signal that is not completely explained by existing dynamical and empirical ocean tide models. So a new empirical model for the effect of the termensual (Mtm and mtm), fortnightly (Mf and mf), and monthly (Mm) tides on polar motion is derived here by fitting periodic terms at these tidal frequencies to polar motion excitation observations that span 2 January 1980 to 8 September 2006 and from which atmospheric and non-tidal oceanic effects have been removed. While this new empirical tide model can fully explain the observed fortnightly polar motion excitation signal during this time interval it would still be desirable to have a model for the effect of long-period ocean tides on polar motion that is determined from a dynamical ocean tide model and that is therefore independent of polar motion observations.     

Atmospheric tides and rotation of the Earth

The six-hourly values of the atmospheric angular momentum (AAM) functions computed by the U.S. National Meteorological Center (NMC) were used to estimate the effects of the atmospheric tides on the Earth's rotation. Variations of the equatorial components ₁ and ₂ of the AAM have periods close to gravitational tidesP ₁ andK ₁.The amplitudes of the detected variations in ₁ and ₂ functions have been found to be much larger than the theoretical ones, the reason of this amplification remains unexplained. According to theoretical formulations, these waves can be expressed only as retrograde motions. Because of frame effects, there is a correspondance between diurnal retrograde polar motion and precession-nutations and the atmospheric effect on polar motion cannot be detected from observations.The second part of this paper deals the effects of atmospheric tides in Earth rotation. High-frequency UT1 variations have been derived from VLBI and GPS techniques during the SEARCH'92 campaign (Study ofEarth-AtmosphereRapidCHanges) (Dickey et al. 1994). They have been compared to values derived by Ray et al. (1994) from global ocean tide model. The results obtained in the present paper show the existence of variations of thermal origin with an amplitude of about 1µs in Universal Time UT1. The agreement between observed and theoretical values is better when the determined thermal atmospheric tides are taken into account.Oceanic tidal signal explains a large part (60% of the signal variance) of the diurnal and sub-diurnal variations. Our results show that only a small part of the residuals (5%) accounts for the atmospheric tidal effects. The residual signal remains unexplained; it might be due to mismodelization of oceanic or atmospheric tides or effect of other geophysical phenomena.     

Predictability of the Earth’s polar motion

The numerical prediction of the Earth’s polar motion is of both theoretical and practical interest. The present paper is aimed at a comprehensive, experimental study of the predictability of polar motion using a homogeneous BIH (Bureau International de l’Heure) data set for the period 1967–1983. Based on our knowledge of the physics of the annual and the Chandler wobbles, we build the numerical model for the polar motion by allowing the wobble period to vary. Using an optimum base length of six years for prediction, this “floating-period” model, equipped with a nonlinear least-squares estimator, is found to yield polar motion predictions accurate to within 0″.012 to 0″.024 depending on the prediction length up to one year, corresponding to a predictability of 89–82%. This represents a considerable improvement over the conventional fixed-period predictor, which, by its nature, does not respond to variations in the apparent wobble periods (in particular, a dramatic decrease in the periods of both the annual and the Chandler wobbles after the year 1980). The superiority of the floating-period predictor to other predictors based on critically different numerical models is also demonstrated.     

Combinations of Earth-orientation measurements: SPACE97, COMB97, and POLE97

Independent Earth-orientation measurements taken by the space-geodetic techniques of lunar and satellite laser ranging, very-long-baseline interferometry, and the global positioning system have been combined using a Kalman filter. The resulting combined Earth-orientation series, SPACE97, consists of values and uncertainties for universal time, polar motion, and their rates spanning the period 28 September 1976 to 3 January 1998 at daily intervals. The space-geodetic measurements used to generate SPACE97 have then been combined with optical astrometric measurements to form two additional combined Earth-orientation series: (1) COMB97, consisting of values and uncertainties for universal time, polar motion, and their rates spanning the period 20 January 1962 to 1 January 1998 at 5-day intervals, and (2) POLE97, consisting of values and uncertainties for polar motion and its rate spanning the period 20 January 1900 to 21 December 1997 at monthly intervals. Received: 10 August 1998 / Accepted: 31 May 1999     

Polar motion modeling,analysis, and prediction with time dependent harmonic coefficients

A time dependent amplitude model was proposed for the analysis and prediction of polar motion time series. The formulation was implemented to analyze part of the new combined solution, EOP (IERS) C 04, daily polar motion time series of 14 years length using a statistical model with first order autoregressive disturbances. A new solution approach, where the serial correlations of the disturbances are eliminated by sequentially differencing the measurements, was used to estimate the model parameters using weighted least squares. The new model parsimoniously represents the 14-year time series with 0.5 mas rms fit, close to the reported 0.1 mas observed pole position precisions for the x and y components. The model can also predict 6 months into the future with less than 4 mas rms prediction error for both polar motion components, and down to sub mas for one-step ahead prediction as validated using a set of daily time series data that are not used in the estimation. This study is dedicated to the memory of Prof. Urho Uotila (1923–2006) whose teaching of “Adjustment Computations” over the years influenced so much, so many of us who had the privilege of being his students.     

EOP and scale from continuous VLBI observing: CONT campaigns to future VGOS networks

Continuous (CONT) VLBI campaigns have been carried out about every 3 years since 2002. The basic idea of these campaigns is to acquire state-of-the-art VLBI data over a continuous time period of about 2 weeks to demonstrate the highest accuracy of which the current VLBI system is capable. In addition, these campaigns support scientific studies such as investigations of high-resolution Earth rotation, reference frame stability, and daily to sub-daily site motions. The size of the CONT networks and the observing data rate have increased steadily since 1994. Performance of these networks based on reference frame scale precision and polar motion/LOD comparison with global navigation satellite system (GNSS) earth orientation parameters (EOP) has been substantially better than the weekly operational R1 and R4 series. The precisions of CONT EOP and scale have improved by more than a factor of two since 2002. Polar motion precision based on the WRMS difference between VLBI and GNSS for the most recent CONT campaigns is at the 30 \(\upmu \)as level, which is comparable to that of GNSS. The CONT campaigns are a natural precursor to the planned future VLBI observing networks, which are expected to observe continuously. We compare the performance of the most recent CONT campaigns in 2011 and 2014 with the expected performance of the future VLBI global observing system network using simulations. These simulations indicate that the expected future precision of scale and EOP will be at least 3 times better than the current CONT precision.