Earth’s precession–nutation motion: the error analysis of the theories IAU 2000 and IAU 2006 applying the VLBI data of the years 1984–2006

The long-term systematic errors of the analytical theories IAU 2000 and IAU 2006 of the Earth’s precession–nutational motion are studied making use of the VLBI data of 1984–2007. Several independent methods give indubitable evidence of the significant quadratic error in the IAU 2000 residuals of the precessional angle while the adopted value of the secular decrease /cy of the Earth’s ellipticity e (derived from Satellite Laser Ranging data) should manifest itself in the residuals of as the negative quadratic trend . The problem with the precession of the IAU 2006 theory adopted as a new international standard and based on the precession model P03 (Capitaine et al., Astron Astrophys 432:355–367, 2005) appears to be even more serious because the above mentioned quadratic term has already been incorporated into the P03 precession. Our analysis of the VLBI data demonstrates that the quadratic trend of the IAU 2006 residuals does amount to the expected value (30.0 ± 3) mas/cy². It means, first, that the theoretical precession rate of IAU 2006 should be augmented by the large secular correction and, second, that the available VLBI data have potentiality of estimating the rate . And indeed, processing these data by the numerical theory ERA of the Earth’s rotation (Krasinsky, Celest Mech Dyn Astron 96:169–217, 2006, Krasinsky and Vasilyev, Celest Mech Dyn Astron 96:219–237, 2006) yields the estimate /cy statistically in accordance with the satellite-based . On the other hand, applying IAU 2000/2006 models, the positive value /cy is found which is incompatible with the SLR estimate and, evidently, has no physical meaning. The large and steadily increasing error of the precession motion of the IAU 2006 theory makes the task of replacing IAU 2006 by a more accurate model be most pressing.     

Tidal effects in the earth–moon system and the earth's rotation

Analytical expressions for tidal torques induced by a tide‐arising planet which perturbs rotation of a nonrigid body are derived. Corresponding expressions both for secular and periodic perturbations of the Euler's angles are given for the case of the earth's rotation. Centennial secular rates of the nutation angle θ and of the earth's angular velocity ω, as well as the centennial logarithmic decrement ν of the Chandler wobble are evaluated:  mas, . In the Universal Time (UT) a large out‐of‐phase (sine) dissipative term with the period 18.6 years and the amplitude 2.3 ms is found. Corrections to nutation coefficients, which presumably have not been taken into account in IAU theory, are given. This revised version was published online in July 2006 with corrections to the Cover Date.     

Completely Integrable Systems Connected with Lie Algebras

In a recent paper Ballersteros and Ragnisco (1998) have proposed a new method of constructing integrable Hamiltonian systems. A new class of integrable systems may be devised using the following sequence: , where A is a Lie algebra is a Lie–Poisson structure on R ₃, C is a Casimir for is a reduced Poisson bracket and (A, ▵) is a bialgebra. We study the relation between a Lie-Poisson stucture Λ and a reduced Poisson bracket , which is a key element in using the Lie algebra A to constructing this sequence. New examples of Lie algebras and their related integrable Hamiltonian systems are given. This revised version was published online in July 2006 with corrections to the Cover Date.     

K emission-line widths and the solar chromosphere

Closely spaced microphotometer tracings parallel to the dispersion of one excellent frame of a K-line time sequence have been utilized for a study of the nature of the K_2v , K_2R intensities in the case of the solar chromosphere. The frequency of occurrence of the categories of intensity ratio are as follows: per cent; per cent; per cent; per cent; per cent. Two types of absorbing components are postulated to explain the pattern of observed K_2v , k_2R intensity ratios. One component with minor Doppler displacements acting on the normal K₂₃₂ profile, where K_2V >K_2R , produces the cases K_2v K_2R , K_2v = K_2R , K_2v <K_2R . The other component arises from dark condensations which are of size 3500 kms as seen in K_2R . They have principally large down flowing velocities in the range 5–8 km/sec and are seen on K₃ spectroheliograms with sizes of about 5000 kms, within the coarse network of emission. These dark condensations give rise to the situation K_2R = 0.K₂-line widths are measured for all tracings where K_2v , K_2R are measurable simultaneously. The distribution curve of these widths is extremely sharp. The K₂ emission source is identified with the bright fine mottles visible on the surface. Evidence for this interpretation comes from the study of auto-correlation functions of K₂ intensity variations and the spacing between the bright fine mottles from both spectrograms and spectroheliograms. The life time of the fine mottling is 200 sec.The supergranular boundaries which constitute the coarse network come in two intensity classes. A low intensity network has the fine mottles as its principal contributor to the K emission. When the network is bright, the enhancement is caused by increased K emission due to the accumulation of magnetic fields at the supergranule boundary. The K₂ widths of the low intensity supergranular boundary agree with the value found for the bright mottles. Those for the brighter network are lower than this value, similar to the K₂ widths as seen in the active regions.It is concluded that bright fine mottling is responsible for the relation, found by Wilson and Bappu, between K emission line widths and absolute magnitudes of the stars.The paper discusses the solar cycle equivalents that stellar chromospheres can demonstrate and indicates a possible line of approach for successful detection of cyclic activity in stellar chromospheres.     

Motion On the Sphere: Integrability and Families of Orbits

We study the problem of the motion of a unit mass on the unit sphere and examine the relation between integrability and certain monoparametric families of orbits. In particular we show that if the potential is compatible with a family of meridians, it is integrable with an integral linear in the velocities, while a family of parallels guarantees integrability with an integral quadratic in the velocities.     

Periodic solutions for resonance 4/3: Application to the construction of an intermediary orbit for Hyperion’s motion

The motion of Hyperion is an almost perfect application of second kind and second genius orbit, according to Poincaré’s classification. In order to construct such an orbit, we suppose that Titan’s motion is an elliptical one and that the observed frequencies are such that 4n _H−3n _T+3n _ω=0, where n _H, n _T are the mean motions of Hyperion and Titan, n _ω is the rate of rotation of Hyperion’s pericenter. We admit that the observed motion of Hyperion is a periodic motion such as . Then, .N _H, N _T, k∈ N ⁺. With that hypothesis we show that Hyperion’s orbit tends to a particular periodic solution among the periodic solutions of the Keplerian problem, when Titan’s mass tends to zero. The condition of periodicity allows us to construct this orbit which represents the real motion with a very good approximation. This revised version was published online in July 2006 with corrections to the Cover Date.     

Numerical theory of rotation of the deformable Earth with the two-layer fluid core. Part 2: Fitting to VLBI data

VLBI-based offsets of the Celestial Pole positions, as well as the variations of UT (series of Goddard Space Flight Center, 1984–2005) are processed applying the Earth’s rotation theory (ERA) 2005 constructed by the numerical integration of the differential equations of rotation of the deformable Earth. The equations were published earlier (Krasinsky 2006) as the first part of the work. The resulting weighted root mean square (WRMS) errors of the residuals , for the angles of nutation and precession are 0.136 and 0.129 mas, respectively. They are significantly less than the corresponding values 0.172 and 0.165 mas for the IAU 2000 model adopted as the international standard. In ERA 2005, the angles , are related to the inertial ecliptical frame J2000, the angle including the precessional secular motion. As the published observational data are theory-dependent being related to IAU 2000, a procedure to confront the numerical theory to the observed Celestial Pole offsets and UT variations is developed. Processing the VLBI data has shown that beside the well known 435-day FCN mode of the free core nutation, there exits a second mode, FICN, caused by the inner part of the fluid core, with the period of 420 day close to that of the FCN mode. Beatings between the two modes are responsible for the apparent damping and excitation of the free oscillations, and are implicitly modeled by ERA 2005. The nutational and precessional motions in ERA 2005 are proved to be mutually consistent but only in case the relativistic correction for the geodetic precession is applied. Otherwise, the overall WRMS error of the residuals would increase by 35%. Thus, the effect of the geodetic precession in the Earth rotation is confirmed experimentally. The other finding is the reliable estimation δ_c = 3.844 ± 0.028° of the phase lag δ_c of the tides in the fluid core. When processing the UT variations, a simple model of the elastic interaction between the mantle and fluid core at their common boundary made it possible to satisfactory describe the largest observed oscillations of UT with the period of 18.6 year, reducing the WRMS error of the UT residuals to the value 0.18 ms (after removing the secular, annual and semi-annual terms).     

Contribution of Vacuum Field to Angular Deviation of Light Path and Radar Echo Delay

The discovery of ‘twin quasistellar objects’ arose interests among astronomers and astrophysicists to study gravitational lensing problems. The deviation of light from its straight line path is caused by two sources according to the general theory of relativity: (i) the presence of massive objects, i.e. the presence of gravitational field and (ii) the presence of a ‘vacuum field’ which arises because there is a non-zero cosmological vacuum energy. Recently, the research on the relationship between cosmological constant and gravitational lensing process is rather active (see reference [1, 2, 3]. According to the Kottler space time metric, we have deduced an explicit representation of the angular deviation of light path. The deviation term is found to be simply , where M is the mass of the ‘astronomical lens’, r_min is the distance between the point of nearest approach and the centre of M, other symbols have their usual meaning. The presence of this term may be meaningful to the study of cosmological constant using the concept of gravitational lensing; however more sophisticated analysis awaits. Consider a signal radar to be sent from one planet to another. We have found that the radar echo delay contributed by the existence of the cosmological constant Λ is expressible as This revised version was published online in July 2006 with corrections to the Cover Date.     

[CII] Emission and Star Formation in Late-Type Galaxies

A linear correlation between the ratio of the[CII( $^{\text{2}}$ P A linear correlation between the ratio of the[CII( P → P )] line intensity to the [ CO(J:1 →0)] line emission, I /I and the equivalent width (EW) is found, over the range 2–71 ? in EW, for a sample of 21late-Type= galaxies. The latter is comprised of an optically selected sample of 12 normal Virgo Cluster spiral galaxies with [CII] detections obtained by us with ISOLWS, plus nine late-Type= galaxies with higher star formation rates (SFRs), for which [CII] data and, especially, EW data are available in the literature. As a result we infer I /I to be a reliable tracer of the current mass-normalized global SFR for non-starburst spiral galaxies. Moreover, the ratio of the [CII] line to the total far-infrared (FIR) continuum intensity, I /I , is found to decrease from ∼0.5% to ∼0.1% with decreasing SFR which we propose is due to a `[CII]-quiet' component of I from dust heated by the general interstellar radiation field (ISRF). The more `quiescent' galaxies in the sample have values of I /I different from those observed in `compact' Galactic interstellar regions. Their [CII]-emission is interpreted to be dominated by diffuse regions of the interstellar medium (ISM). For normal `star-forming' galaxies the diffuse component of the [CII] emission is estimated to account for at least 50% of the total. This revised version was published online in July 2006 with corrections to the Cover Date.     

3He enrichments by magnetosonic waves in current-driven instabilities

A theoretical model for³He enrichments in solar energetic particles is developed. First, current-driven, electrostatic instabilities that have frequencies ( is the cyclotron frequency of³He) are investigated for a plasma consisting of H,⁴He,³He, and electrons with the density of³He much lower than those of H and⁴He. It is found that in many cases the oblique ion-acoustic waves can have positive growth rates at frequencies and, at the same time, negative growth rates at and at _H. This can occur near the marginal state of the instability. The wave damping at these frequencies is caused by the cyclotron resonances of⁴He and H. The cyclotron damping at is negligible, however, because the abundance of³He is very small. The H cyclotron waves can be unstable at for a wide region of plasma parameters; the electron-to-ion temperature ratio must beT _e /T _H 1.5. To destabilize the₄He cyclotron waves with , high⁴He density and high electron temperature are both required. Then,³He enrichments are studied on the basis of the theory of nonlinear magnetosonic waves, which can promptly accelerate ions. The current-driven electrostatic waves with can enhance fluctuation velocities of³He. Thus, in the presence of these waves, magnetosonic waves can selectively accelerate³He particles to high energies. Finally, cyclotron resonances of heavy ions with the waves or are briefly discussed.     

Rotational evolution of exoplanets under the action of gravitational and magnetic perturbations

We investigate the evolution of the rotational axes of exoplanets under the action of gravitational and magnetic perturbations. The planet is assumed to be dynamically symmetrical and to be magnetised along its dynamical-symmetry axis. By qualitative methods of the bifurcation theory of multiparametric PDEs, we have derived a gallery of 69 phase portraits. The portraits illustrate evolutionary trajectories of the angular momentum of a planet for a variety of the initial conditions, for different values of the ratio between parameters describing gravitational and magnetic perturbations, and for different rates of the orbital evolution. We provide examples of the phase portraits, that reveal the differences in topology and the evolutionary track of in the vicinity of an equilibrium state. We determine the bifurcation properties, i.e., the way of reorganisation of phase trajectories in the vicinities of equilibria; and we point out the combinations of parameters’ values that permit ip-overs from a prograde to a retrograde spin mode.     

Secular frequencies of 3-D exoplanetary systems

Using a 12th order expansion of the perturbative potential in powers of the eccentricities and the inclinations, we study the secular effects of two non-coplanar planets which are not in mean–motion resonance. By means of Lie transformations (which introduce an action–angle formulation of the Hamiltonian), we find the four fundamental frequencies of the 3-D secular three-body problem and compute the long-term time evolutions of the Keplerian elements. To find the relations between these elements, the main combinations of the fundamental frequencies common to these evolutions are identified by frequency analysis. This study is performed for two different reference frames: a general one and the Laplace plane. We underline the known limitations of the linear Laplace–Lagrange theory and point out the great sensitivity of the 3-D secular three-body problem to its initial values. This analytical approach is applied to the exoplanetary system Andromedae in order to search whether the eccentricities evolutions and the apsidal configuration (libration of ) observed in the coplanar case are maintained for increasing initial values of the mutual inclination of the two orbital planes. Anne-Sophie Libert is FNRS Research Fellow.     

Deceleration of the Earth's rotation from old solar tables

The solar tables of ibn Yunis and of King Alfonso, those of Kepler, of G. D. and J.-J. Cassini, and of Lalande are compared with Newcomb's theory of the Sun to determine the deceleration of the Earth's rotation. Comparisons of mean motion and of longitude lead to separate determinations. A value for the deceleration, , assumed to be independent of time is obtained for the period — 146 through 1892. This result is based on the assumption that Newcomb's theory converts the Earth's orbital motion into a perfect clock, and the result is independent of lunar and tidal theory. The deceleration seems greater than that obtained from the literature for the tidal deceleration caused by the Moon.     

Meridional flow in thick accretion disks

We have studied the effect of the flow in the accretion disk. The specific angular momentum of the disk is assumed to be constant and the polytropic relation is used. We have solved the structure of the disk and the flow patterns of the irrotational perfect fluid.As far as the obtained results are concerned, the flow does not affect the shape of the configuration in the bulk of the disk, although the flow velocity reaches even a half of the sound velocity at the inner edge of the disk. Therefore, in order to study accretion disk models with the moderate mass accretion rate—i.e.,

Effect of oblateness on the non-linear stability of L 4 in the restricted three-body problem

Non-linear stability of the libration point L ₄ of the restricted three-body problem is studied when the more massive primary is an oblate spheroid with its equatorial plane coincident with the plane of motion, Moser's conditions are utilised in this study by employing the iterative scheme of Henrard for transforming the Hamiltonian to the Birkhoff's normal form with the help of double D'Alembert's series. It is found that L ₄ is stable for all mass ratios in the range of linear stability except for the three mass ratios: $$\begin{gathered} \mu _{c1} = 0.0242{\text{ }}...{\text{ }}{}^{{\text{\_\_}}}0.1790{\text{ }}...{\text{ }}A_1 , \hfill \\ \mu _{c2} = 0.0135{\text{ }}...{\text{ }}{}^{{\text{\_\_}}}0.0993{\text{ }}...{\text{ }}A_1 , \hfill \\ \mu _{c3} = 0.0109{\text{ }}...{\text{ }}{}^{{\text{\_\_}}}0.0294{\text{ }}...{\text{ }}A_1 . \hfill \\ \end{gathered} $$      

Impact of the Quadrupole Moment of the Sun on the Dynamics of the Earth-Moon System

Range of values of the Sun's mass quadrupole moment of coefficient J₂ arising both from experimental and theoretical determinations enlarge across literature on two orders of magnitude, from around 10^-7 until to 10^-5. The accurate knowledge of the Moon's physical librations, for which the Lunar Laser Ranging data reach an outstanding precision level, prove to be appropriate to reduce the interval of J₂ values by giving an upper bound of J₂. A solar quadrupole moment as high as 1.1 10^-5 given either from the upper bounds of the error bars of the observations, or from the Roche's theory, is not compatible with the knowledge of the lunar librations accurately modeled and observed with the LLR experiment. The suitable values of J₂ have to be smaller than 3.0 10^-6. As a consequence, this upper bound of 3.0 10^-6 is accepted to study the impact of the Sun's quadrupole moment of mass on the dynamics of the Earth-Moon system. Such as effect (with J₂ = 5.5±1.3 × 10^-6) has been already tested in 1983 by Campbell & Moffat using analytical approximate equations, and thus for the orbits of Mercury, Venus, the Earth and Icarus. The approximate equations are no longer sufficient compared with present observational data and exact equations are required. As if to compute the effect on the lunar librations, we have used our BJV relativistic model of solar system integration including the spin-orbit coupled motion of the Moon. The model is solved by numerical integration. The BJV model stems from general relativity by using the DSX formalism for purposes of celestial mechanics when it is about to deal with a system of n extended, weakly self-gravitating, rotating and deformable bodies in mutual interactions. The resulting effects on the orbital elements of the Earth have been computed and plotted over 160 and 1600 years. The impact of the quadrupole moment of the Sun on the Earth's orbital motion is mainly characterized by variations of , , and . As a consequence, the Sun's quadrupole moment of mass could play a sensible role over long time periods of integration of solar system models. This revised version was published online in July 2006 with corrections to the Cover Date.     

The restricted two-body problem in constant curvature spaces

We perform the bifurcation analysis of the Kepler problem on and . An analog of the Delaunay variables is introduced. We investigate the motion of a point mass in the field of a Newtonian center moving along a geodesic on and (the restricted two-body problem). For the case of a small curvature, the pericenter shift is computed using the perturbation theory. We also present the results of numerical analysis based on an analogy with the motion of a rigid body.