形变地球对不同驱动力的响应

地球形变位移场和重力场的时空变化无论在基础理论研究,还是在地理空间信息建设中都具有重要的意义.地球在各种力学机制的作用下产生了形变,形变又导致地球引力位的变化,即形变附加位或Euler引力位增量.基于矢量球函数的基本理论,讨论了引潮力、负载力和地表应力对地球形变和引力位增量的影响,给出了均匀不可压缩地球模型的Euler引力位增量的具体表达式和Love数的理论关系.可为地球形变的理论研究提供参考和依据.     

Mutual gravitational potential ofN solid bodies

The mutual gravitational potential ofN solid bodies is expanded without approximation in terms of harmonic coefficients of each body. As an application the Euler dynamical equations for the motion of the axis of figure of the rigid Earth are integrated analytically by the method of variation of parameters.     

Nonsingular expansions of the gravity potential and its derivatives at satellite altitudes in the ellipsoidal coordinate system

The series in ellipsoidal harmonics for derivatives of the Earth’s gravity potential are used only on the reference ellipsoid enveloping the Earth due to their very complex mathematical structure. In the current study, the series in ellipsoidal harmonics are constructed for first- and second-order derivatives of the potential at satellite altitudes; their structure is similar to the series on the reference ellipsoid. The point P is chosen at a random satellite altitude; then, the ellipsoid of revolution is described, which passes through this point and is confocal to the reference ellipsoid. An object-centered coordinate system with the origin at the point P is considered. Using a sequence of transformations, the nonsingular series in ellipsoidal harmonics is constructed for first and second derivatives of the potential in the object-centered coordinate system. These series can be applied to develop a model of the Earth’s potential, based on combined use of surface gravitational force measurements, data on the satellite orbital position, its acceleration, or measurements of the gravitational force gradients of the first and second order. The technique is applicable to any other planet of the Solar System.     

地球的弹性形变能

采用较完善的地球模型PREM，通过求解弹性体的基本运动方程，得到形变位移矢量，由此分别计算出因日月引潮力势和地球自转离心力势引起的地球的弹性形变能，它将为弹性地球运动的Hamilton表达式提供摄动项的量级估计，并将是理论上探讨地月系演化的一个重要方面。     

A Hamiltonian theory for an elastic earth: Elastic energy of deformation

In this paper we study only the perturbation due to the deformation of the elastic mantle by a tidal body force. In a previous publication (Getino and Ferrándiz, 1989a) we defined two canonical systems of variables - we gave them the names ofelastic variables of Euler and Andoyer respectively. Next, using them, we obtained the canonical expression of rotational kinetic energy. In the present paper, using the same variables, we build up the elastic energy which is produced by the deformation of the elastic mantle. We show that the three termsm = 0, 1, 2 corresponding to the second order of the development in spherical harmonics of the perturbing potential, a tidal potential, are of the same order of magnitude. In addition, the numerical integration for a particular Earth Model (Takeuchi's Model 2) is performed, with the aim of obtaining a numerical estimate of the coefficients which intervene in both this energy and the previously mentioned kinetic energy.     

A Hamiltonian theory for an elastic earth: Canonical variables and kinetic energy

This paper is the first of a set of four, in which we shall develop the first part of a project dedicated to elaborating a Hamiltonian theory for the rotational motion of a deformable Earth. Here we study only the perturbation due to the deformation of the elastic mantle by tidal body force. In the present paper, we define two canonical systems of variables—we give these variables the names of elastic variables of Euler and Andoyer respectively. Next, using them, we obtain the canonical expression of rotational kinetic energy, which is valid for any Earth model satisfying hypotheses as general as those established in Section 2.     

The Earth-Moon tidal force function

The concept “the tidal force function of the Earth-Moon system” is introduced and its exact determination based on the Stokes constants (harmonic coefficients) in the external gravitational potential of both bodies is outlined. The exact determination of the torque due to the Moon exerted on the Earth may be performed in terms of the Stokes constants of both bodies and the mutual position of both ellipsoids of inertia.     

月球卫星轨道变化的分析解

由于月球自转缓慢及其引力位的特点,使得讨论月球卫星与人造地球卫星轨道变化的方法有所不同。     

引力常数变化对地球自转长期变化的影响

探讨和估计了各种引力常数变化理论对地球角速度和日长变化的影响。各种引力常数变化理论包括了引力常数G随时间、空间以及速度变化等几个方面的影响。另外也估计了对地球自转角速度和日长变化产生的效应。其中有些研究对探讨地球自转变化也有启发意义。     

Motion of a Meteoroid Released from an Asteroid

Evidence of asteroid surface features as regolith grains and larger boulders implies resurfacing possibility due to external forces such as gravitational tidal force during close planet encounters. Motion of a meteoroid released from an asteroid in the gravitational fields of the asteroid and the Earth is modeled. We are interested mainly in a distance between the meteoroid and the asteroid as a function of the time. Applications to Itokawa and some close approaching NEAs are presented.     

Solution to the rotation of the elastic earth by method of rigid dynamics

Hamiltonian mechanics is applied to the problem of the rotation of the elastic Earth. We first show the process for the formulation of the Hamiltonian for rotation of a deformable body and the derivation of the equations of motion from it. Then, based on a simple model of deformation, the solution is given for the period of Euler motion, UT1 and the nutation of the elastic Earth. In particular it is shown that the elasticity of the Earth acts on the nutation so as to decrease the Oppolzer terms of the nutation of the rigid Earth by about 30 per cent. The solution is in good agreement with results which have been obtained by other, different approaches.     

Effect of the Earth’s compression on the physical libration of the Moon

The effect of the Earth??s compression on the physical libration of the Moon is studied using a new vector method. The moment of gravitational forces exerted on the Moon by the oblate Earth is derived considering second order harmonics. The terms in the expression for this moment are arranged according to their order of magnitude. The contribution due to a spherically symmetric Earth proves to be greater by a factor of 1.34 × 10⁶ than a typical term allowing for the oblateness. A linearized Euler system of equations to describe the Moon??s rotation with allowance for external gravitational forces is given. A full solution of the differential equation describing the Moon??s libration in longitude is derived. This solution includes both arbitrary and forced oscillation harmonics that we studied earlier (perturbations due to a spherically symmetric Earth and the Sun) and new harmonics due to the Earth??s compression. We posed and solved the problem of spinorbital motion considering the orientation of the Earth??s rotation axis with regard to the axes of inertia of the Moon when it is at a random point in its orbit. The rotation axes of the Earth and the Moon are shown to become coplanar with each other when the orbiting Moon has an ecliptic longitude of L _? = 90° or L _? = 270°. The famous Cassini??s laws describing the motion of the Moon are supplemented by the rule for coplanarity when proper rotations in the Earth-Moon system are taken into account. When we consider the effect of the Earth??s compression on the Moon??s libration in longitude, a harmonic with an amplitude of 0.03?? and period of T ₈ = 9.300 Julian years appears. This amplitude exceeds the most noticeable harmonic due to the Sun by a factor of nearly 2.7. The effect of the Earth??s compression on the variation in spin angular velocity of the Moon proves to be negligible.     

Physics of the Earth and otons

Short-time variations of gravitational potential derivatives (otonic gravity-impulses) are described which are produced by fast-moving otons (objects of general relativity) in the Earth. Expressions for oton mass are obtained through measurable physical quantities. The question of otonic gravity-impulses registration is analysed.     

Motion in a modified Chermnykh’s restricted three-body problem with oblateness

In this paper, the restricted problem of three bodies is generalized to include a case when the passively gravitating test particle is an oblate spheroid under effect of small perturbations in the Coriolis and centrifugal forces when the first primary is a source of radiation and the second one an oblate spheroid, coupled with the influence of the gravitational potential from the belt. The equilibrium points are found and it is seen that, in addition to the usual three collinear equilibrium points, there appear two new ones due to the potential from the belt and the mass ratio. Two triangular equilibrium points exist. These equilibria are affected by radiation of the first primary, small perturbation in the centrifugal force, oblateness of both the test particle and second primary and the effect arising from the mass of the belt. The linear stability of the equilibrium points is explored and the stability outcome of the collinear equilibrium points remains unstable. In the case of the triangular points, motion is stable with respect to some conditions which depend on the critical mass parameter; influenced by the small perturbations, radiating effect of the first primary, oblateness of the test body and second primary and the gravitational potential from the belt. The effects of each of the imposed free parameters are analyzed. The potential from the belt and small perturbation in the Coriolis force are stabilizing parameters while radiation, small perturbation in the centrifugal force and oblateness reduce the stable regions. The overall effect is that the region of stable motion increases under the combine action of these parameters. We have also found the frequencies of the long and short periodic motion around stable triangular points. Illustrative numerical exploration is rendered in the Sun–Jupiter and Sun–Earth systems where we show that in reality, for some values of the system parameters, the additional equilibrium points do not in general exist even when there is a belt to interact with.     

Yukawa-type potential effects in the anomalistic period of celestial bodies

Several contemporary modified models of gravity predict the existence of a non-Newtonian Yukawa-type correction to the classical gravitational potential. We study the motion of a secondary celestial body under the influence of the corrected gravitational force of a primary. We derive two equations to approximate the periastron time rate of change and its total variation over one revolution (i.e., the difference between the anomalistic period and the Keplerian period) under the influence of the non-Newtonian radial acceleration. Kinematically, this influence produces apsidal motion. We performed numerical estimations for Mercury, for the companion star of the pulsar PSR 1913+16, and for the extrasolar Planet b of the star HD 80606. We also considered the case of the artificial Earth satellite GRACE-A, but the results present a low degree of reliability from a practical standpoint.     

The lunar capture hypothesis revisited

Recent work on planetary formation processes have suggested that ancient planetary bodies could have been warmer and, therefore, more easily deformable soon after formation than at present. By use of the estimates for the elastic parameters believed to be appropriate for a warm ancient Moon and Earth, it is shown that the energy of deformation of the planetary bodies during a close gravitational encounter was sufficient to effect capture.     

A new model of angular momentum transfer from the rotating central body of a two‐body system to the orbital motion of the system

In a previous paper we treated within the framework of our Projective Unified Field Theory (Schmutzer 2004, 2005a) the 2‐body system (e.g. Earth‐Moon system) with a rotating central body in a rather abstract manner. Here a concrete model of the transfer of angular momentum from the rotating central body to the orbital motion of the whole 2‐body system is presented, where particularly the transfer is caused by the inhomogeneous gravitational force of the Moon acting on the oceanic waters of the Earth, being modeled by a spherical shell around the solid Earth. The theory is numerically tested. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Approximate analytic method for high-apogee twelve-hour orbits of artificial Earth’s satellites

We propose an approach to the study of the evolution of high-apogee twelve-hour orbits of artificial Earth’s satellites. We describe parameters of the motion model used for the artificial Earth’s satellite such that the principal gravitational perturbations of the Moon and Sun, nonsphericity of the Earth, and perturbations from the light pressure force are approximately taken into account. To solve the system of averaged equations describing the evolution of the orbit parameters of an artificial satellite, we use both numeric and analytic methods. To select initial parameters of the twelve-hour orbit, we assume that the path of the satellite along the surface of the Earth is stable. Results obtained by the analytic method and by the numerical integration of the evolving system are compared. For intervals of several years, we obtain estimates of oscillation periods and amplitudes for orbital elements. To verify the results and estimate the precision of the method, we use the numerical integration of rigorous (not averaged) equations of motion of the artificial satellite: they take into account forces acting on the satellite substantially more completely and precisely. The described method can be applied not only to the investigation of orbit evolutions of artificial satellites of the Earth; it can be applied to the investigation of the orbit evolution for other planets of the Solar system provided that the corresponding research problem will arise in the future and the considered special class of resonance orbits of satellites will be used for that purpose.     

Analytical Solution of Coupled Perturbation of Tesseral Harmonic Terms of Mars's Non-Spherical Gravitational Potential

The non-spherical gravitational potential of the planet Mars is sig- nificantly different from that of the Earth. The magnitudes of Mars’ tesseral harmonic coefficients are basically ten times larger than the corresponding val- ues of the Earth. Especially, the magnitude of its second degree and order tesseral harmonic coefficient J_2,2 is nearly 40 times that of the Earth, and approaches to the one tenth of its second zonal harmonic coefficient J₂. For a low-orbit Mars probe, if the required accuracy of orbit prediction of 1-day arc length is within 500 m (equivalent to the order of magnitude of 10−4 standard unit), then the coupled terms of J2 with the tesseral harmonics, and even those of the tesseral harmonics themselves, which are negligible for the Earth satellites, should be considered when the analytical perturbation solution of its orbit is built. In this paper, the analytical solutions of the coupled terms are presented. The anal- ysis and numerical verification indicate that the effect of the above-mentioned coupled perturbation on the orbit may exceed 10⁻⁴ in the along-track direc- tion. The conclusion is that the solutions of Earth satellites cannot be simply used without any modification when dealing with the analytical perturbation solutions of Mars-orbiting satellites, and that the effect of the coupled terms of Mars's non-spherical gravitational potential discussed in this paper should be taken into consideration.     

Optimizing expansions of the Earth’s gravity potential and its derivatives over ellipsoidal harmonics

A set of spherical harmonics is the most widely used representation of the Earth’s gravity potential. This series converges outside and on the surface of a reference sphere enveloping the Earth. However, the Earth’s surface is better approximated by the reference ellipsoid—a compressed ellipsoid of revolution that covers the entire Earth. The gravity potential can be expanded in a series over ellipsoidal harmonics on the surface of the reference ellipsoid and on the surface of other external confocal ellipsoids of revolution. In contrast to spherical harmonics, depending on the associated Legendre functions of the first kind, ellipsoidal harmonics depend also on the associated Legendre functions of the second kind. The latter contain the very slowly converging hypergeometric Gauss series. The number of series increases with increasing the order of their derivatives. In this work, we derived new series for the gravitational potential of the Earth and its derivatives over ellipsoidal harmonics. Starting from the first order derivative, all the series corresponding to higher order derivatives depend on the same two hypergeometric Gauss series. The latter converges considerably faster than that for the hypergeometric series previously used when computing the gravity potential and its derivatives.