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.     

Comparison of length of day with oceanic and atmospheric angular momentum series

This is a companion paper to earlier comparisons and study of operational polar motion series, published recently in the same journal. In this contribution, four operational, publicly available, length-of-day (LOD) time series have been compared to the atmospheric angular momentum (AAM) augmented with recent oceanic angular momentum (OAM) data during September 1997–July 2000, using several intervals ranging from 3 days to almost 3 years. Additionally, the LOD of the International GNSS Service (IGS) historical series and a new LOD combination (CMB) were also analyzed. All the six LOD series showed an overall correlation exceeding 0.99 for the complete interval of almost 3 years. Even for the shortest interval of only 3 days, the correlation was still higher than 0.60. The combined AAM + OAM series with inverted barometer corrections always gave the best correlation. The Rapid Service LOD of the International Earth Rotation and Reference Systems Service (IERS) compared the best at all intervals but the shortest one, where the CMB LOD was the best with a correlation of 0.73, followed by both IGS series with a correlation of about 0.71. Prior to all the correlation analyses, in addition to the removal of all the known (conventional) LOD tidal variations with periods ranging from 5.6 days to 18.6 years and lunar fortnightly and monthly oceanic tides, small corrections of lunar fortnightly and monthly tides, semi-annual, annual periodical signals, drift and scale had to be estimated with respect to the combined AAM + OAM series.     

Diurnal atmospheric forcing and temporal variations of the nutation amplitudes

A 29-year time-series of four-times-daily atmospheric effective angular momentum (EAM) estimates is used to study the atmospheric influence on nutation. The most important atmospheric contributions are found for the prograde annual (77 μas), retrograde annual (53 as), prograde semiannual (45 as), and for the constant offset of the pole (δψsinɛ₀=−86 as, δɛ=77 as). Among them only the prograde semiannual component is driven mostly by the wind term of the EAM function, while in all other cases the pressure term is dominant. These are nonnegligible quantities which should be taken into account in the new theory of nutation. Comparison with the VLBI corrections to the IAU 1980 nutation model taking into account the ocean tide contribution yields good agreement for the prograde annual and semiannual nutations. We also investigated time variability of the atmospheric contribution to the nutation amplitudes by performing the sliding-window least-squares analysis of both the atmospheric excitation and VLBI nutation data. Almost all detected variations of atmospheric origin can be attributed to the pressure term, the biggest being the in-phase annual prograde component (about 30 as) and the retrograde one (as much as 100200 as). These variations, if physical, limit the precision of classical modeling of nutation to the level of 0.1 mas. Comparison with the VLBI data shows significant correlation for the retrograde annual nutation after 1989, while for the prograde annual term there is a high correlation in shape but the size of the atmospherically driven variations is about three times less than deduced from the VLBI data. This discrepancy in size can be attributed either to inaccuracy of the theoretical transfer function or the frequency-dependent ocean response to the pressure variations. Our comparison also yields a considerably better agreement with the VLBI nutation data when using the EAM function without the IB correction for ocean response, which indicates that this correction is not adequate for nearly diurnal variations. Received: 10 September 1997 / Accepted: 5 March 1998     

Angular momentum of Arctic Ocean tides

Oceanic tidal angular momentum (OTAM) is calculated for the four major tides of the Arctic Ocean, based on the tidal elevations and current velocities from a recent two-dimensional numerical hydrodynamic model. The presented OTAM tables are meant to be complementary to other modeling studies that use satellite altimetry (which cannot observe Arctic Ocean tides because of ice cover and limited satellite inclinations). Although the Arctic Ocean's influence on earth rotation is, as may be expected, relatively small, the rapid advancement of the subject now calls for such small contributions to be explicitly accounted for. Received: 22 January 1996; Accepted: 5 December 1996     

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.     

Atmospheric Angular Momentum Time-Series: Characterization of their Internal Noise and Creation of a Combined Series

The atmosphere induces variations in Earth rotation. These effects are classically computed using the “angular momentum approach”. In this method, the variations in Earth rotation are estimated from the variations in the atmospheric angular momentum (AAM). Several AAM time-series are available from different meteorological centers. However, the estimation of atmospheric effects on Earth rotation differs when using one atmospheric model or the other. The purpose of this work is to build an objective criterion that justifies the use of one series in particular. Because the atmosphere is not the only cause of Earth rotation variations, this criterion cannot rely only on a comparison of AAM series with geodetic data. Instead, we determine the quality of each series by making an estimation of their noise level, using a generalized formulation of the “three-cornered hat method”. We show the existence of a link between the noise of the AAM series and their correlation with geodetic data: a noisy series is usually less correlated with Earth orientation data. As the quality of the series varies in time, we construct a combined AAM series, using time-dependent weights chosen so that the noise level of the combined series is minimal. To determine the influence of a minimal noise level on the correlation with geodetic data, we compute the correlation between the combined series and Earth orientation data. We note that the combined series is always amongst the best correlated series, which confirms the link established before. The quality criterion, while totally independent of Earth orientation observations, appears to be physically convincing when atmospheric and geodetic data are compared     

Atmospheric, oceanic and hydrological contributions to seasonal variations in length of day

 The annual and semiannual residuals derived in the axial angular momentum budget of the solid Earth–atmosphere system reflect significant signals. They must be caused by further excitation sources. Since, in particular, the contribution for the wind term from the atmospheric layer between the 10 and 0.3 hPa levels to the seasonal variations in length of day (LOD) is still missing, it is necessary to extend the top level into the upper stratosphere up to 0.3 hPa. Under the conservation of the total angular momentum of the entire Earth, variations in the oceanic angular momentum (OAM) and the hydrological angular momentum (HAM) are further significant excitation sources at seasonal time scales. Focusing on other contributions to the Earth's axial angular momentum budget, the following data are used in this study: axial atmospheric angular momentum (AAM) data derived for the 10–0.3 hPa layer from 1991 to 1997 for computing the missing wind effects; axial OAM functions as generated by oceanic general circulation models (GCMs), namely for the ECHAM3 and the MICOM models, available from 1975 to 1994 and from 1992 to 1994, respectively, for computing the oceanic contributions to LOD changes, and, concerning the HAM variations, the seasonal estimates of the hydrological contribution as derived by Chao and O'Connor [(1988) Geophys J 94: 263–270]. Using vector representation, it is shown that the vectors achieve a close balance in the global axial angular momentum budget within the estimated uncertainties of the momentum quantities on seasonal time scales. Received: 6 April 2000 / Accepted: 13 December 2000     

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.     

Temporal variations of the polar moment of inertia and the second-degree geopotential

In order to study the geodynamic behaviour of the Earth over short (elastic Earth) and long (almost perfectly liquid Earth) geological periodic variations, the changes of the moment of inertia are decomposed into two parts: the first, described by a volume integral, explains the effect of the density variations, while the second gives the impact of the surface variations using a surface integral. It is shown that both components have physical significance, but their contribution is different in case of short (lunisolar) and long (connected to secular despinning) periods.     

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.     

Wavelet analysis of interannual LOD, AAM, and ENSO: 1997–98 El Niño and 1998–99 La Niña signals

 On the basis of the data series of the length of day (LOD), the atmospheric angular momentum (AAM) and the Southern Oscillation Index (SOI) for January 1970–June 1999, the relationship among Interannual LOD, AAM, and the EL Ni?o/Southern Oscillation (ENSO) is analyzed by the wavelet transform method. The results suggest that they have similar time-varying spectral structures. The signals of 1997–98 El Ni?o and 1998–99 La Ni?a events can be detected from the LOD or AAM data. Received: 25 January 2000 / Accepted: 9 January 2001     

The impact of errors in polar motion and nutation on UT1 determinations from VLBI Intensive observations

The earth’s phase of rotation, expressed as Universal Time UT1, is the most variable component of the earth’s rotation. Continuous monitoring of this quantity is realised through daily single-baseline VLBI observations which are interleaved with VLBI network observations. The accuracy of these single-baseline observations is established mainly through statistically determined standard deviations of the adjustment process although the results of these measurements are prone to systematic errors. The two major effects are caused by inaccuracies in the polar motion and nutation angles introduced as a priori values which propagate into the UT1 results. In this paper, we analyse the transfer of these components into UT1 depending on the two VLBI baselines being used for short duration UT1 monitoring. We develop transfer functions of the errors in polar motion and nutation into the UT1 estimates. Maximum values reach 30 [μs per milliarcsecond] which is quite large considering that observations of nutation offsets w.r.t. the state-of-the-art nutation model show deviations of as much as one milliarcsecond.     

基于新参考系的极移改正

本文利用IAU2000决议中关于CIP的定义阐述了极移坐标系和极移坐标的定义,在总结微分旋转矩阵的性质的基础上给出了极移旋转矩阵的详细推导,结合极移产生的主要原因,详细介绍了目前国际上关于极移模型化工作的最新进展,最后给出了获取极移坐标和TIO位置的方法和途径。     

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.     

Analysis of the first year of Earth rotation parameters with a sub-daily resolution gained at the CODE processing center of the IGS

 The solutions of the CODE Analysis Center submitted to the IGS, the International Global Position System (GPS) Service for Geodynamics, are based on three days of observation of about 80–100 stations of the IGS network. The Earth rotation parameters (ERPs) are assumed to vary linearly over the three days with respect to an a priori model. Continuity at the day boundaries as well as the continuity of the first derivatives are enforced by constraints. Since early April 1995 CODE has calculated a new ERP series with an increased time resolution of 2 hours. Again continuity is enforced at the 2-hours-interval boundaries. The analysis method is described, particularly how to deal with retrograde diurnal terms in the ERP series which may not be estimated with satellite geodetic methods. The results obtained from the first year of data covered by the time series (time interval from 4 April 1995 to 30 June 1996) are also discussed. The series is relatively homogeneous in the sense of the used orbit model and the a priori model for the ERPs. The largest source of excitation at daily and sub-daily periods is likely to be the effect of the ocean tides. There is good agreement between the present results and Topex/Poseidon ocean tide models, as well as with models based on Very Long Baseline Interferometry (VLBI) and Satellite Laser Ranging (SLR) data. Non-oceanic periodic variations are also observed in the series. Their origin is most probably a consequence of the GPS solution strategy; other possible sources are the atmospheric tides. Received: 13 July 1999 / Accepted: 21 March 2000     

An historical empirical line method for the retrieval of surface reflectance factor from multi-temporal SPOT HRV,HRVIR and HRG multispectral satellite imagery

SPOT satellites have been imaging Earth's surface since SPOT 1 was launched in 1986. It is argued that absolute atmospheric correction is a prerequisite for quantitative remote sensing. Areas where land cover changes are occurring rapidly are also often areas most lacking in situ data which would allow full use of radiative transfer models for reflectance factor retrieval (RFR). Consequently, this study details the proposed historical empirical line method (HELM) for RFR from multi-temporal SPOT imagery. HELM is designed for use in landscape level studies in circumstances where no detailed overpass concurrent atmospheric or meteorological data are available, but where there is field access to the research site(s) and a goniometer or spectrometer is available. SPOT data are complicated by the ±27° off-nadir cross track viewing. Calibration to nadir only surface reflectance factor (ρ_s) is denoted as HELM-1, whilst calibration to ρ_s modelling imagery illumination and view geometries is termed HELM-2. Comparisons of field measured ρ_s with those derived from HELM corrected SPOT imagery, covering Helsinki, Finland, and Taita Hills, Kenya, indicated HELM-1 RFR absolute accuracy was ±0.02ρ_s in the visible and near infrared (VIS/NIR) bands and ±0.03ρ_s in the shortwave infrared (SWIR), whilst HELM-2 performance was ±0.03ρ_s in the VIS/NIR and ±0.04ρ_s in the SWIR. This represented band specific relative errors of 10–15%. HELM-1 and HELM-2 RFR were significantly better than at-satellite reflectance (ρ_SAT), indicating HELM was effective in reducing atmospheric effects. However, neither HELM approach reduced variability in mean ρ_s between multi-temporal images, compared to ρ_SAT. HELM-1 calibration error is dependent on surface characteristics and scene illumination and view geometry. Based on multiangular ρ_s measurements of vegetation-free ground targets, calibration error was negligible in the forward scattering direction, even at maximum off-nadir view. However, error exceeds 0.02ρ_s where off-nadir viewing was ≥20° in the backscattering direction within ±55° azimuth of the principal plane. Overall, HELM-1 results were commensurate with an identified VIS/NIR 0.02ρ_s accuracy benchmark. HELM thus increases applicability of SPOT data to quantitative remote sensing studies.