国内天文地球动力学中的潮汐研究

简要说明了天文地球动力学范畴内所研究的潮汐现象，包括由日月引潮力引起的固体潮、海洋潮、大气潮和由于地球自转轴的极移引起的极潮，以及这些潮汐对地球自转和地球自转的测量产生的效应。重点阐述中国天文学界在这一领域里的研究成果。这些研究涉及潮汐影响地球自转的机制，也就是各种潮汐效应与极移、自转速率变化和章动的关系，包括构建这类关系的理论模型，分析潮汐对它们的影响，利用中国古代丰富的天象记录计算地球自转的长期减慢，计算弹性或滞弹地球的洛夫数，依据某一地球模型计算潮汐效应或章动序列等等。研究也涉及在测量地球自转参数的不同技术中各种潮汐效应对测量结果产生的影响及其改正，并涉及与潮汐有关的观测方法的优化和数据处理过程的改进。最后介绍了中国学者所发现的脉冲星的周期和周期变率测量中的潮汐效应，尽管它们的量级甚微，但不容忽视。     

Influence of the Atmospheric Tides on the Earth Rotation

Spectral analysis of the components of the relative atmospheric angular momentum vector is performed based on the series of these components for the 6 h intervals within the period of 1958–2000. These series have been computed in the Subbureau of the Atmospheric Angular Momentum of the International Earth Rotation Service using the NCEP/NCAR reanalysis of atmospheric observations. The basic harmonics of diurnal tides are determined. New results on the fortnight's and week's duration oscillations of the equatorial components of the atmospheric angular momentum are obtained. The zonal tides transformation mechanisms in the atmosphere are discussed. It is shown that the main mechanism of the zonal tides effect on the atmospheric variability is the amplitude modulation of daily oscillations of the relative atmospheric angular momentum. The effects of the atmospheric tides on the Earth rotation are discussed.     

由PREM模型参数计算地球自转的周期变化

弹性地球在日月引潮力势作用下的形变引起其转动惯量的改变，从而导致地球自转速率的变化．本文利用PREM地球模型所给的物质密度和弹性等参数分布．计算日月引潮力势产生的地球形变附加势，进而计算转动惯量的变化．最后得到一系列包含周期同引潮势带谐项、振幅大于1微秒的自转速率周期变化系数．     

Dynamical History of the Earth–Moon System

Differential equations describing the tidal evolution of the earth's rotation and of the lunar orbital motion are presented in a simple close form. The equations differ in form for orbits fixed to the terrestrial equator and for orbits with the nodes precessing along the ecliptic due to solar perturbations. Analytical considerations show that if the contemporary lunar orbit were equatorial the evolution would develop from an unstable geosynchronous orbit of the period about 4.42 h (in the past) to a stable geosynchronous orbit of the period about 44.8 days (in the future). It is also demonstrated that at the contemporary epoch the orbital plane of the fictitious equatorial moon would be unstable in the Liapunov's sense, being asymptotically stable at early stages of the evolution. Evolution of the currently near-ecliptical lunar orbit and of the terrestrial rotation is traced backward in time by numerical integration of the evolutional equations. It is confirmed that about 1.8 billion years ago a critical phase of the evolution took place when the equatorial inclination of the moon reached small values and the moon was in a near vicinity of the earth. Before the critical epoch t _cr two types of the evolution are possible, which at present cannot be unambiguously distinguished with the help of the purely dynamical considerations. In the scenario that seems to be the most realistic from the physical point of view, the evolution also has started from a geosynchronous equatorial lunar orbit of the period 4.19 h. At t < t _cr the lunar orbit has been fixed to the precessing terrestrial equator by strong perturbations from the earth's flattening and by tidal effects; at the critical epoch the solar perturbations begin to dominate and transfer the moon to its contemporary near-ecliptical orbit which evolves now to the stable geosynchronous state. Probably this scenario is in favour of the Darwin's hypothesis about originating the moon by its separation from the earth. Too much short time scale of the evolution in this model might be enlarged if the dissipative Q factor had somewhat larger values in the past than in the present epoch. Values of the length of day and the length of month, estimated from paleontological data, are confronted with the results of the developed model.     

A New Symbolic Processor for the Earth Rotation Theory

In this paper we present a new symbolic processor specially suited for the Earth rotation theory. This processor works with a more general kind of Poisson series called Kinoshita series, which has resulted to be very useful in the Earth rotation theory. Its structure is adapted for dealing with the more general analytical expressions that appear in the Earth rotation theory. This new algebraic processor has been successfully used for computing different contributions to the nutation series of the rigid Earth.This revised version was published online in October 2005 with corrections to the Cover Date.     

对月日潮汐摩擦耗散的估计

月日潮汐摩擦和地球惯量矩变化是日长长期变化的主要原因．在本文中，利用最新的地球物理和古生物钟数据，对过去15亿年以来的月日潮汐摩擦、地球惯量矩变化和日长长期变化等作了数值对比研究．由此得到二个重要结论：一是仅利用地球的自转形变不能解释J2的变化，这说明地球的重力分异现象至今仍存在着；其二是在几亿年前的潮汐摩擦比现在大得多，若取潮汐耗散与距离的立方成反比时，理论结果与由古生物钟得到的回归年日数和朔望月日数数据较为符合。     

地球自转动力学中的两个问题

介绍了地球自动力学中的两个目前正在研究的问题，1．Chandler摆动的随机激发，在分析Chandler摆动各种激发的可能性后，认为随机运动是最可能的激发源，在此基础,提出了一个Chandler摆动激发的动力学模型，并从理论和数值模拟两方面对此模型做了统计分析研究，描述了今后对此问题的研究思路，2．地球内部动力学是目前国际地球动力学界的一个热点研究课题，介绍了它的现状和最近的发展动态以及准备在这方面开展研究工作的打算。     

Secular changes of LOD associated with a growth of the inner core

From recent estimates of the age of the inner core based on the theory of thermal evolution of the core, we estimate that nowadays the growth of the inner core may perhaps contribute to the observed overall secular increase of LOD caused mainly by tidal friction (i.e., 1.72 ms per century) by a relative decrease of 2 to 7 µs per century. Another, albeit much less plausible, hypothesis is that crystallization of the inner core does not produce any change of LOD, but makes the inner core rotate differentially with respect to the outer core and mantle. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Improved Estimation of Earth Rotation Parameters Using the Adaptive Ridge Regression

The multicollinearity among regression variables is a common phenomenon in the reduction of astronomical data. The phenomenon of multicollinearity and the diagnostic factors are introduced first. As a remedy, a new method, called adaptive ridge regression (ARR), which is an improved method of choosing the departure constant θ in ridge regression, is suggested and applied in a case that the Earth orientation parameters (EOP) are determined by lunar laser ranging (LLR). It is pointed out, via a diagnosis, the variance inflation factors (VIFs), that there exists serious multicollinearity among the regression variables. It is shown that the ARR method is effective in reducing the multicollinearity and makes the regression coefficients more stable than that of using ordinary least squares estimation (LS), especially when there is serious multicollinearity. This revised version was published online in July 2006 with corrections to the Cover Date.     

A Check on the Variations of Earth‘s Rotation with an Ancient Solar Eclipse

We address the relation between an ancient total eclipse, which occurred on A.D.1542 August 11 and the variation of Earth‘‘s rotation. The total eclipse was recorded in some ancient Chinese books, especially in local chronicles. Some of the documents include useful information for determining the location of the totality zone. The parameters of the eclipse are calculated by using the DE406 Ephemeris.A high-precision value of AT which expresses the variation of the Earth‘‘s rotation,of about 300 ～ 380 s, is obtained.     

地球自转变化的重力效应

极移对重力测量的影响在前人的研究和工作中已得到了充分的考虑并有不同形式的改正公式。但是由于地轴转动速率变化的重力效应的量级很小,且限于以前的测量精度以及资料处理方法, 因而常常被忽略。随着重力测量精度的提高以及资料处理方法的改进, 有必要重新讨论该项影响。本文从推导地球瞬时自转极运动对地面台站重力测量值影响的公式入手, 将地球自转角速率变化的重力效应和极移的重力效应分离开来, 并对它们的量级分别作了估计。本文着重讨论自转角速率变化对重力测量的影响。用 Ｉ Ｅ Ｒ Ｓ 的日长变化序列和极移序列分别具体地计算了１９８２ ～１９９６ 年期间对武汉的重力测量值的影响, 结果表明, 自转角速率变化对武汉地区重力值的影响最大可达０ ．２５ ｍｉｃｒｏｇａｌ（ｐｅａｋ －ｔｏ － ｐｅａｋ） , 这在高精度的绝对重力测量和相对重力测量中不应再被忽略     

中国地球自转和地壳运动监测的研究工作

主要介绍了１９９５年至１９９８年期间有关中国地球自转和地壳运动监测的研究工作及取 得的进展。     

On the effect of the mantle elasticity on the earth's rotation

Starting from the Hamiltonian model for a solid Earth with an elastic mantle previously developped by the authors, analytical expressions are derived which give the nutation series corresponding to the plane perpendicular to the angular momentum vector, to the plane perpendicular to the rotational axis and to the equator of figure, as well as the series that give the polar motion. The effects of the different perturbations — solid Earth, centrifugal and tidal potentials — are calculated separately. The corrections due to the elasticity of the mantle, which mostly correspond to the Oppolzer terms, are calculated with an accuracy of 10^–6 arc sec., given that the intrinsic observational accuracy has reached 0.01 mas.     

On the “rotational angular momentum” of the oceans and the corresponding polar motion

Periodic polar motions caused by ocean tides are predicted. In the Liouville equations for rotational motion the complete excitation functions for the ocean tides have to be used. This does not depend on the fact that hydrodynamical ocean tide models do not consider the centrifugal acceleration. The observable polar motion of the Celestial Ephemeris Pole CEP (more exactly: the terrestrial location of the CEP) is tabulated for the ten ocean tides M₂, S₂, N₂, K₁, O₁, P₁, M f, M f′, M m, Ssa. Typical amplitudes for the largest ocean tides are 0.4 milliarcseconds. This is within the reach of geodetic VLBI and SLR observations.     

大气和带谐潮汐对日长高频变化的影响

用国际地球自转服务提供的高精度日长序列ΔLOD(IERS(EOP99C04),和由美国国家环境预报中心(NCEP)对过去的气象资料进行重新分析和计算,提供的轴向大气角动量序列(AAM)来研究大气和带谐潮汐对地球自转变化(ΔLOD)中高频振荡激发的贡献。由Multitaper功率谱分析表明,ΔLOD序列在9.13d(天)、13.63d、13.66d、27.55d周期时呈现出比较强的谱峰,这些谱都在95%的显著性水平以上,有比较强的信噪比。从ΔLOD和AAM之差值序列谱分析结果及地球带谐响应系数k的估计,表明大气不是ΔLOD的高频变化的主要激发源。根据理论潮汐周期可知,有可能存在其它的激发源,如海洋对地球自转变化的激发。这还有待于进一步深入探讨。     

Effects of the Triaxiality on the Rotation of Celestial Bodies: Application to the Earth,Mars and Eros

In this paper we discuss the influence of the triaxiality of a celestialbody on its free rotation, i.e. in absence of any external gravitationalperturbation. We compare the results obtained through two different analytical formalisms, one established from Andoyer variables by usingHamiltonian theory, the other one from Euler's variables by usingLagrangian equations. We also give a very accurate formulation of thepolar motion (polhody) in the case of a small amplitude of this motion.Then, we carry out a numerical integration of the problem, with aRunge–Kutta–Felberg algorithm, and for the two kinds of methods above, that we apply to three different celestial bodies considered as rigid : the Earth, Mars, and Eros. The reason of this choice is that each of this body corresponds to a more or less triaxial shape.In the case of the Earth and Mars we show the good agreement betweenanalytical and numerical determinations of the polar motion, and theamplitude of the effect related to the triaxial shape of the body, whichis far from being negligible, with some influence on the polhody of theorder of 10 cm for the Earth, and 1 m for Mars. In the case of Eros, weuse recent output data given by the NEAR probe, to determine in detailthe nature of its free rotational motion, characterized by the presence ofimportant oscillations for the Euler angles due to the particularly largetriaxial shape of the asteroid.     

在后牛顿引力理论(PPN形式)中引力常数变化对地球自转产生的效应

介绍和论述了在后牛顿引力理论(PPN形式)中在优越参考系和非优越参考系中经过参数化后引力常数变化对地球自转产生的效应,其中特别重点介绍了年周期变化的效应。此外也将理论结果同观测结果相对比。     

The theory of canonical perturbations applied to attitude dynamics and to the Earth rotation. Osculating and nonosculating Andoyer variables

In the method of variation of parameters we express the Cartesian coordinates or the Euler angles as functions of the time and six constants. If, under disturbance, we endow the “constants” with time dependence, the perturbed orbital or angular velocity will consist of a partial time derivative and a convective term that includes time derivatives of the “constants”. The Lagrange constraint, often imposed for convenience, nullifies the convective term and thereby guarantees that the functional dependence of the velocity on the time and “constants” stays unaltered under disturbance. “Constants” satisfying this constraint are called osculating elements. Otherwise, they are simply termed orbital or rotational elements. When the equations for the elements are required to be canonical, it is normally the Delaunay variables that are chosen to be the orbital elements, and it is the Andoyer variables that are typically chosen to play the role of rotational elements. (Since some of the Andoyer elements are time-dependent even in the unperturbed setting, the role of “constants” is actually played by their initial values.) The Delaunay and Andoyer sets of variables share a subtle peculiarity: under certain circumstances the standard equations render the elements nonosculating. In the theory of orbits, the planetary equations yield nonosculating elements when perturbations depend on velocities. To keep the elements osculating, the equations must be amended with extra terms that are not parts of the disturbing function [Efroimsky, M., Goldreich, P.: J. Math. Phys. 44, 5958–5977 (2003); Astron. Astrophys. 415, 1187–1199 (2004); Efroimsky, M.: Celest. Mech. Dyn. Astron. 91, 75–108 (2005); Ann. New York Acad. Sci. 1065, 346–374 (2006)]. It complicates both the Lagrange- and Delaunay-type planetary equations and makes the Delaunay equations noncanonical. In attitude dynamics, whenever a perturbation depends upon the angular velocity (like a switch to a noninertial frame), a mere amendment of the Hamiltonian makes the equations yield nonosculating Andoyer elements. To make them osculating, extra terms should be added to the equations (but then the equations will no longer be canonical). Calculations in nonosculating variables are mathematically valid, but their physical interpretation is not easy. Nonosculating orbital elements parameterise instantaneous conics not tangent to the orbit. (A nonosculating i may differ much from the real inclination of the orbit, given by the osculating i.) Nonosculating Andoyer elements correctly describe perturbed attitude, but their interconnection with the angular velocity is a nontrivial issue. The Kinoshita–Souchay theory tacitly employs nonosculating Andoyer elements. For this reason, even though the elements are introduced in a precessing frame, they nevertheless return the inertial velocity, not the velocity relative to the precessing frame. To amend the Kinoshita–Souchay theory, we derive the precessing-frame-related directional angles of the angular velocity relative to the precessing frame. The loss of osculation should not necessarily be considered a flaw of the Kinoshita–Souchay theory, because in some situations it is the inertial, not the relative, angular velocity that is measurable [Schreiber, K. U. et al.: J. Geophys. Res. 109, B06405 (2004); Petrov, L.: Astron. Astrophys. 467, 359–369 (2007)]. Under these circumstances, the Kinoshita–Souchay formulae for the angular velocity should be employed (as long as they are rightly identified as the formulae for the inertial angular velocity).     

液核地球自转速率变化的带谐潮效应Ⅰ--理论公式

根据液核地球动力学原理,重新推导了在日月子引潮力作用下的地球自转速率变化的理论公式,由此明确地引入了液核影响因子,它与地球和液核的极转动惯量的变化有关。不同于Yoder.C.F(1981)的理论,本文给出了基于引潮力位Doodson展开的地球自转速率变化、日长变化和世界时变化的理论公式,指出尺度因子的转动惯量应取为地幔的有效极转动惯量,Love数应取地幔的有效Love数的理由。     

液核地球自转速率变化的带谐潮效应Ⅱ--数值计算与比较

地球自转速率的潮汐变化可由无量纲参数k/cm)(k和cm分别为壳幔的有效Love数和有效极转动惯量)来表征。对于一个具有弹性地幔、平衡海潮和核幔不耦合的地球k/cm=0.944,且与潮波频率无关。海潮的非平衡扰动使k/cm为复数,且与频率有关。大气对自转速率有效勒夫数的贡献约为Δkat=0.0075。同时地幔滞弹性对勒夫数也产生扰动。利用本文得到的理论公式和最新的潮汐数据计算了地球自转速率的潮汐变化,及其有关地球物理机制的影响。