同时顾及章动和极移的地球自转方程

通过引进章动坐标系相对惯性参照系的转动角速度随时间的变化 ,导出了一个可同时解出章动和极移的地球自转方程 ,用这个方程可同时研究地球的强迫和自由转动。与现行研究地球自转的惯用方法相比 ,该方法综合性强 ,易于理解。     

液核地球极移和章动理论模型的改进

以经典的液核弹性地球自转动力学理论为基础,通过引入章动坐标系相对惯性空间的运动,建立液核地球极移和章动的联合动力学方程。由此定义液核地球的CIP轴,并对其进行解算。研究表明新方法克服了传统方法不能解决的一些问题(例如CEP和CIP的定义),简化了传统方法求解的复杂性。值得注意的是新方法可同时求解极移和章动,特别是在Smith M.L.(1977)理论中出现的倾斜模(TOM),在此只是作为一个特解而存在。     

Polar motion for an elastic earth model using a canonical theory

The purpose of this paper is to develop a canonical formulation of the rotational motion for an elastic Earth model. We have obtained the canonical equations for the precession and nutation motion in an inertial frame, and from this we have deduced the equations in an Earth-fixed frame. The linearized equations deduced for polar motion are equivalent to those obtained using Liouville's equations.     

Transformation of amplitudes and frequencies of precession and nutation of the earth’s rotation vector to amplitudes and frequencies of diurnal polar motion

The temporal change of the rotation vector of a rotating body is, in the first order, identical in a space-fixed system and in a body-fixed system. Therefore, if the motion of the rotation axis of the earth relative to a space-fixed system is given as a function of time, it should be possible to compute its motion relative to an earth-fixed system, and vice versa. This paper presents such a transformation. Two models of motion of the rotation axis in the space-fixed system are considered: one consisting only of a regular (i.e., strictly conical) precession and one extended by circular nutation components, which are superimposed upon the regular precession. The Euler angles describing the orientation of the earth-fixed system with respect to the space-fixed system are derived by an analytical solution of the kinematical Eulerian differential equations. In the first case (precession only), this is directly possible, and in the second case (precession and nutation), a solution is achieved by a perturbation approach, where the result of the first case serves as an approximation and nutation is regarded as a small perturbation, which is treated in a linearized form. The transformation by means of these Euler angles shows that the rotation axis performs in the earth-fixed system retrograde conical revolutions with small amplitudes, namely one revolution with a period of one sidereal day corresponding to precession and one revolution with a period which is slightly smaller or larger than one sidereal day corresponding to each (prograde or retrograde) circular nutation component. The peculiar feature of the derivation presented here is the analytical solution of the Eulerian differential equations.     

Influence of subdaily polar motion model on nutation offsets estimated by very long baseline interferometry

This paper studies the connection between the subdaily model for polar motion used in the processing of very long baseline interferometry (VLBI) observations and the estimated nutation offsets. By convention accepted by the International Earth Rotation Service, the subdaily model for polar motion recommended for routine processing of geodetic observations does not contain any daily retrograde terms due to their one-to-one correlation with the nutation. Nevertheless, for a 24-h VLBI solution a part of the signal contained in the polar motion given by the used subdaily model is numerically mistaken for a retrograde daily sidereal signal. This fictitious retrograde daily signal contributes to the estimated nutation, leading to systematic differences between the nutation offsets from VLBI solutions computed with different subdaily polar motion models. We demonstrate this effect using solutions for all suitable 24-h VLBI sessions over a time span of 11 years (2000–2011). By changing the amplitudes of one tidal term in the underlying subdaily model for polar motion and comparing the estimated parameters to the solutions computed with the unchanged subdaily model, the paper shows and explains theoretically the effects produced by the individual subdaily terms on the VLBI nutation estimates.     

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.     

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.     

Rotation of the Earth as a Triaxial Rigid Body

The Earth is taken as a triaxial rigid body, which rotates freely in the Euclidian space. The starting equations are the Euler dynamic equations, with A smaller than B and B smaller than C. The Euler equations are solved, and the numerical results are provided. In the calculations, the following parameters are used: (C-B)/A=0.003 273 53; (B-A)/C=0.000 021 96; (C-A)/B=0.003 295 49, and the mean angular velocity of the Earth's rotation, ω =0.000 072 921 15 rad/s. Calculations show that, besides the self-rotation of the Earth and the free Euler procession of its rotation, there exists the free nutation: the nutation angle, or the angle between the Earth's momentary rotation axis and the mean axis that periodically change with time. The free nutation is investigated.     

Effects of adopting new precession,nutation and equinox corrections on the terrestrial reference frame

First, the paper is devoted to the effects of adopting new definitive precession and equinox corrections on the terrestrial reference frame: The effect on polar motion is a diurnal periodic term with an amplitude increasing linearly in time; on UT1 it is a linear term. Second, general principles are given the use of which can determine the effects of small rotations (such as precession, nutation or equinox corrections) of the frame of a Conventional Inertial Reference System (CIS) on the frame of the Conventional Terrestrial Reference System (CTS). Next, seven CTS options are presented, one of which is necessary to accommodate such rotations (corrections). The last of these options requiring no changes in the origin of terrestrial longitudes and in UT1 is advocated; this option would be maintained by eventually referencing the Greenwich Mean Sidereal Time to a fixed point on the equator, instead of to the mean equinox of date, the current practice. Accommodating possible future changes in the astronomical nutation is discussed in the last section. The Appendix deals with the effects of differences which may exist between the various CTS's and CIS's (inherent in the various observational techniques) on earth rotation parameters (ERP) and how these differences can be determined. It is shown that the CTS differences can be determined from observations made at the same site, while the CIS differences by comparing the ERP's determined by the different techniques during the same time period. Presented at XVIII General Assembly of the International Astronomical Union Patras, Greece, August 17–26, 1982.     

三轴刚性地球自由旋转欧拉角变化数值模拟研究

假定地球是一个三轴刚性体,在Euclid空间中做自由旋转.在设定主惯性矩A小于B小于C的情况下,求解欧拉运动方程,得到数值解.计算结果表明:地球在除了自转和自由进动之外,同时还存在着自由章动.章动角会随着时间做周期性变化.重点讨论的是章动角的变化.     

Harmonic analysis of the Earth's gravitational field by means of semi-continuous ephemerides of a low Earth orbiting GPS-tracked satellite. Case study: CHAMP

An algorithm for the determination of the spherical harmonic coefficients of the terrestrial gravitational field representation from the analysis of a kinematic orbit solution of a low earth orbiting GPS-tracked satellite is presented and examined. A gain in accuracy is expected since the kinematic orbit of a LEO satellite can nowadays be determined with very high precision, in the range of a few centimeters. In particular, advantage is taken of Newton's Law of Motion, which balances the acceleration vector with respect to an inertial frame of reference (IRF) and the gradient of the gravitational potential. By means of triple differences, and in particular higher-order differences (seven-point scheme, nine-point scheme), based upon Newton's interpolation formula, the local acceleration vector is estimated from relative GPS position time series. The gradient of the gravitational potential is conventionally given in a body-fixed frame of reference (BRF) where it is nearly time independent or stationary. Accordingly, the gradient of the gravitational potential has to be transformed from spherical BRF to Cartesian IRF. Such a transformation is possible by differentiating the gravitational potential, given as a spherical harmonics series expansion, with respect to Cartesian coordinates by means of the chain rule, and expressing zero- and first-order Ferrer's associated Legendre functions in terms of Cartesian coordinates. Subsequently, the BRF Cartesian coordinates are transformed into IRF Cartesian coordinates by means of the polar motion matrix, the precession–nutation matrices and the Greenwich sidereal time angle (GAST). In such a way a spherical harmonic representation of the terrestrial gravitational field intensity with respect to an IRF is achieved. Numerical tests of a resulting Gauss–Markov model document not only the quality and the high resolution of such a space gravity spectroscopy, but also the problems resulting from noise amplification in the acceleration determination process.     

Results of three year observation with a superconducting gravimeter at the GeoForschungsZentrum Potsdam

The superconducting gravimeter (SG) TT70 has been continuously recording gravity data at the GeoForschungsZentrum (GFZ) Potsdam absolute gravity site since July 1992. The recorded data are edited and preprocessed by spike and step detection and elimination and gap filling. An atmospheric pressure correction is carried out on gravity data in the time domain with a complex admittance before tidal fitting. The atmospheric pressure admittance is calculated from tide free output of SG data and local atmospheric pressure using the cross spectral method. The ground water level admittance is determined by a single coefficient. Improvements with these corrections are shown in analysis results. New tidal parameters for Potsdam site are determined and compared with recordings of an Askania spring gravimeter at a nearby site. Deviations against the Wahr-Dehant-model are shown. Polar motion data of the IERS (International Earth Rotation Service, Paris) are used to calculate variations of centrifugal acceleration caused by polar motion (pole tide). These are compared with the corrected tide free output of SG series. For drift determination the polar motion correction is applied on SG data. The nutation period equivalent to the Earth's Nearly Diurnal Free Wobble is calculated from the SG data with a value of T_FCN = (437.4 ± 1.5) sidereal days. This result is compared with those obtained from other SG stations. Received 19 December 1995; Accepted 13 September 1996     

天球有效大气角动量中的周年和半年项估计

为了讨论大气对非刚体地球章动周年和半年项的贡献,本文将大气有效角动量(1980.0～1998.33)从地面参考系转换至空间参考系得到大气的天球有效角动量CEAM,低通滤波后用最小二乘方法拟合得到顺、逆向周年和半年项的振幅估计。     

Analytical versus semi-analytical determinations of the Oppolzer terms for a non-rigid Earth

Estimations of the Oppolzer terms for the angular momentum and rotation axes of a non-rigid Earth are obtained from two different approaches and compared. The first approach is an analytical method which relies on the solutions of the Liouville equations for a two-layer Earth model. The Oppolzer terms are evaluated from analytical expressions. The results are then compared to those calculated from Wahr's theory of nutation for a non-rigid Earth, which is the second approach used. Results are obtained for the main nutation frequencies and for the precession case. The differences between the two solutions are generally quite small (the relative error is most of the time under 8%) and are, for a large part, due to successive approximations and truncation effects during their determination. Departures of the results from the two methods are significantly larger for frequencies near the Free Core Nutation (FCN) resonance. This is particularly true for the Oppolzer terms of the angular momentum axis. The Earth model adopted is a little bit different in each case: for the Liouville system solution, we have limited the model to a homogeneous elastic mantle and a homogeneous liquid core. Another source of some of the small differences in the results is the presence of a solid inner core in Wahr's theory. We confirm through the analytical calculation the strong effect of the core on the Oppolzer terms of the angular momentum axis for a non-rigid Earth at the precession frequency. Finally, an application is given in the determination of the axes' position at J2000 for a non-rigid Earth. Received: 23 February 1998 /Accepted: 18 November 1998     

Hydrological excitation of polar motion by different variables from the GLDAS models

Continental hydrological loading by land water, snow and ice is a process that is important for the full understanding of the excitation of polar motion. In this study, we compute different estimations of hydrological excitation functions of polar motion (as hydrological angular momentum, HAM) using various variables from the Global Land Data Assimilation System (GLDAS) models of the land-based hydrosphere. The main aim of this study is to show the influence of variables from different hydrological processes including evapotranspiration, runoff, snowmelt and soil moisture, on polar motion excitations at annual and short-term timescales. Hydrological excitation functions of polar motion are determined using selected variables of these GLDAS realizations. Furthermore, we use time-variable gravity field solutions from the Gravity Recovery and Climate Experiment (GRACE) to determine the hydrological mass effects on polar motion excitation. We first conduct an intercomparison of the maps of variations of regional hydrological excitation functions, timing and phase diagrams of different regional and global HAMs. Next, we estimate the hydrological signal in geodetically observed polar motion excitation as a residual by subtracting the contributions of atmospheric angular momentum and oceanic angular momentum. Finally, the hydrological excitations are compared with those hydrological signals determined from residuals of the observed polar motion excitation series. The results will help us understand the relative importance of polar motion excitation within the individual hydrological processes, based on hydrological modeling. This method will allow us to estimate how well the polar motion excitation budget in the seasonal and inter-annual spectral ranges can be closed.     

On the accuracy assessment of celestial reference frame realizations

The scatter of the celestial pole offset (CPO) time-series obtained from very long baseline interferometry (VLBI) observations is used as a measure of the accuracy of celestial reference frame (CRF) realizations. Several scatter indices (SI), including some proposed for the first time, are investigated. The first SI is based on residual analysis of CPO series with respect to a free core nutation (FCN) model. The second group of SIs includes Allan deviation and its extensions, which allow the treatment of unequally weighted and multidimensional observations. Application of these criteria to several radio source catalogues (RSCs) showed their ability to perform a preliminary assessment of the quality of each RSC. The 2D Allan deviation estimate seems to be the most sensitive measure. The proposed extensions of Allan deviation, weighted and multidimensional, can also be used for the statistical analysis of other time-series.     

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     

利用SLR与伪距资料综合定轨

以GPS伪距为观测量对GPS35卫星进行定轨 ,然后将SLR与GPS伪距资料综合起来进行定轨 ,并将计算的轨道与IGS精密轨道进行了比较     

Atmospheric excitation of polar motion

The polar motion excited by the fluctuation of global atmospheric angular momentum (AAM) is investigated. Based on the global AAM data, numerical results demonstrate that the fluctuation of AAM can excite the seasonal wobbles (e.g., the 18-month wobble) and the Chandler wobble, which agree well with previous studies. In addition, by filtering the dominant low frequency components, some distinct polar wobbles corresponding to some great diurnal and semi-diurnal atmospheric tides are found.     

Atmospheric excitation of polar motion

The polar motion excited by the fluctuation of global atmospheric angular momentum (AAM) is investigated. Based on the global AAM data, numerical results demonstrate that the fluctuation of AAM can excite the seasonal wobbles (e.g., the 18-month wobble) and the Chandler wobble, which agree well with previous studies. In addition, by filtering the dominant low frequency components, some distinct polar wobbles corresponding to some great diurnal and semi-diurnal atmospheric tides are found.