Determination of frame-dragging using Earth gravity models from CHAMP and GRACE

Using the new generation Earth’s gravity field models EIGEN-2S, GGM01S and EIGEN-GRACE02S generated by the space missions CHAMP and GRACE, we have obtained an accurate measurement of the Lense–Thirring effect with the LAGEOS and LAGEOS II satellites analyzing about 10 years of data with the EIGEN-2S and GGM01S models and about 11 years of data with EIGEN-GRACE02S. This new analysis is in agreement with our previous measurements of the Lense–Thirring effect using the LAGEOS satellites and obtained with the JGM-3 and EGM96 Earth’s models. However, the new determinations are more accurate and, especially, more robust than our previous measurements. In the present analysis we are only using the nodal rates of the two satellites, making no use of the perigee rate, as in our previous analyses. The perigee is affected by a number of non-gravitational perturbations difficult to be modelled and whose impact in the total error budget is not easy to assess. Using the EIGEN-2S model, we obtain a total error budget between 18% and 36% of the Lense–Thirring effect due to all the error sources. Specifically, by using EIGEN-2S, we obtain: μ = 0.85, with a total error between ±0.18 and ±0.36, with GGM01S we get μ = 1.06 with a total error between ±0.19 and ±0.24 and with EIGEN-GRACE02S we obtain μ = 0.99, with a total error between ±0.05 and ±0.1, i.e., between 5% and 10% of the general relativistic predicted value of the Lense–Thirring effect. In addition to the analyses using EIGEN-2S, GGM01S and EIGEN-GRACE02S without the use of the perigee, we have also performed an analysis using the older model EGM96 with our previous method of combining the nodes of the LAGEOS satellites with the perigee of LAGEOS II. However, this analysis was performed over a period of about 10 years, i.e. about 2.5 times longer than any our previous analysis. The result using EGM96 over this longer period of observation agrees with our previous results over much shorter periods and with the EIGEN-2S, GGM01S and EIGEN-GRACE02S measurements of μ. In addition to the accurate determination of frame-dragging and in agreement with our previous analyses of the orbits of the LAGEOS satellites, we have observed, since 1998, an anomalous change in the Earth quadrupole coefficient, J₂ which agrees with recent findings of other authors. This anomalous variation of J₂ is accurately observed both on the node of LAGEOS and LAGEOS II and it is independent of the model used, i.e., it is observed by using the model EGM96 or by using EIGEN-2S, GGM01S or EIGEN-GRACE02S. However, this anomalous variation of the Earth quadrupole coefficient does not affect at all our determination of the Lense–Thirring effect thanks to the total elimination of the J₂-induced errors with our especially devised estimation technique.     

LAGEOS II perigee rate and eccentricity vector excitations residuals and the Yarkovsky-Schach effect

We have analysed LAGEOS II perigee rate and eccentricity vector excitation residuals over a period of about 7.8 years, adjusting and computing the satellite orbit with the full set of dynamical models included in the GEODYN II software code. The long-term behaviour of these orbital residuals appears to be characterised by several distinct frequencies which are a clear signature of the Yarkovsky-Schach perturbing effect. This non-gravitational perturbation is not included in the GEODYN II models for the orbit determination and analysis. Through an independent numerical analysis, and using the new LOSSAM model to represent the spin-axis behaviour of the satellite, we propagated the Yarkovsky-Schach effect on LAGEOS II perigee rate and compared the results obtained with the orbital residuals. We have thus been able to satisfactorily fit the amplitude of the Yarkovsky-Schach effect to the observed residuals. Our approach here has proven very successful with very positive results. We have been able to obtain a fractional reduction of about 40% of the post-fit rms with respect to the pre-fit value. When analysing the eccentricity vector residuals, we have been able to obtain a better result in the case of the real component, with a fractional reduction of the post-fit rms of about 49% of the initial value. The analysis of the effect's imaginary component in the eccentricity vector rate is more complicated and deserves additional scrutiny. In this case we need a deeper study which includes the analysis of other unmodelled and mismodelled effects acting on the imaginary component. The study performed in this paper will be of significant relevance not only for the geophysical applications involving LAGEOS II orbit analysis, but also for a refined re-analysis of the general relativistic precession produced by the Earth angular momentum, i.e., the Lense-Thirring effect.     

Possible Effect of the Earth’s Inertial Induction on the Orbital Decay of LAGEOS

The theory of velocity dependent inertial induction, based upon extended Mach’s principle, has been able to generate many interesting results related to celestial mechanics and cosmological problems. Because of the extremely minute magnitude of the effect its presence can be detected through the motion of accurately observed bodies like Earth satellites. LAGEOS I and II are medium altitude satellites with nearly circular orbits. The motions of these satellites are accurately recorded and the past data of a few decades help to test many theories including the general theory of relativity. Therefore, it is hoped that the effect of the Earth’s inertial induction can have any detectable effect on the motion of these satellites. It is established that the semi-major axis of LAGEOS I is decreasing at the rate of 1.3 mm/d. As the atmospheric drag is negligible at that altitude, a proper explanation of the secular change has been wanting, and, therefore, this paper examines the effect of the Earth’s inertial induction effect on LAGEOS I. Past researches have established that Yarkovsky thermal drag, charged and neutral particle drag might be the possible mechanisms for this orbital decay. Inertial induction is found to generate a perturbing force that results in 0.33 mm/d decay of the semi major axis. Some other changes are also predicted and the phenomenon also helps to explain the observed changes in the orbits of a few other satellites. The results indicate the feasibility of the theory of inertial induction i.e. the dynamic gravitation phenomenon of the Earth on its satellites as a possible partial cause for orbital decay.     

Tidal effects on the LAGEOS–LARES satellites and the LARASE program

The effects of solid and ocean tides have been computed on the right ascension of the ascending node of the two LAGEOS and LARES satellites and on the argument of pericenter of LAGEOS II. Their effects—together with the possible mis-modeling related to systematic errors in the estimate of the tidal coefficients, especially in the case of ocean tides—are quite important to be well established for the key role of the LAGEOS satellites, as well as of the newly LARES, in space geodesy and geophysics as well as in fundamental physics measurements. In the case of the measurement of the Lense–Thirring effect, the mis-modeling of long-period tides may mimic a secular effect on the cited orbital elements, thus producing a degradation in the measurement of the relativistic precession. A suitable combination of the orbital elements of the three satellites can help in avoiding the effects of the long-period tides of degree \(\ell =2\) (as for the Lunar solid tides with periods of 18.6 and 9.3 years) and \(\ell =4\), but other long-period tides, as the ocean \(K_1\) tide, which has the same periodicities of the right ascension of the ascending node \(\varOmega \) of the satellites, may strongly influence the measurement, especially if it is performed over a relatively short time span. These results are particularly important in the case of LARES, since they are new and because of the role that the orbit of LARES, and especially of its ascending node right ascension, will have in a new measurement of the Lense–Thirring effect by the joint analysis of its orbit with that of the two LAGEOS.     

The LAGEOS satellites orbital residuals determination and the Lense-Thirring effect measurement

The method applied since 1996 for the analysis of the orbital residuals of the LAGEOS satellites in order to measure the Lense-Thirring effect has been the subject of the present work. This method, based on the difference between the orbital elements of consecutive arcs, is explained and analysed also from the analytical point of view. It is proved that this “difference method” works well for the determination of the secular effects, as in the case of the relativistic precession induced by the Earth's angular momentum, but also very useful for the determination and study of the long-term periodic effects. Indeed, the only limitation in the determination of the periodic effects is the possibility of the reduction of their amplitude by a factor which depends from the periodicity of the given perturbation and from the orbital arc length chosen for the satellite during the data analysis. In the case of the Yarkovsky-Schach effect, the main non-gravitational perturbation seen in the LAGEOS satellites orbital residuals, in particular in its perigee rate and eccentricity vector excitation residuals, we show that the “difference method” is quite good also for the determination of the long-period perturbations induced by this subtle non-conservative force.     

On the reliability of the so-far performed tests for measuring the Lense–Thirring effect with the LAGEOS satellites

In this paper, we critically discuss the so-far performed attempts aimed at the detection of the general relativistic gravitomagnetic Lense–Thirring effect in the gravitational field of the Earth with the existing LAGEOS satellites. In the latest reported measurement of the gravitomagnetic shift with the nodes of the LAGEOS satellites and the second generation GRACE-only EIGEN-GRACE02S Earth gravity model over an observational time span of 11 years a 5–10% total accuracy is claimed at 1–3σ, respectively. We will show that, instead, it might be 15–45% (1–3σ) if the impact of the secular variations of the even zonal harmonics is considered. Possible strategies in order both to make more robust and reliable the tests with the node-only LAGEOS–LAGEOS II combination used and to overcome the problems affecting it with other alternative combinations are presented.     

The geomagnetic effects on the motion of an electrically charged artificial satellite

The orbital effects of the Lorentz force on the motion of an electrically charged artificial satellite moving in the Earth's magnetic field are determined. The geomagnetic field is considered as a multipole potential field and the satellite electrical charge is supposed to be constant. The relativistic perturbations of the main geomagnetic field are discussed briefly. The results are concentrated on the determination of the secular changes, and numerical values are computed for the case of the LAGEOS satellite. The results are discussed in the context of a possible detection of the Lense-Thirring effect analyzing the orbital perturbations of the LAGEOS and LAGEOS X satellites.     

On the use of Ajisai and Jason-1 satellites for tests of general relativity

In this paper we analyze in detail some aspects of the proposed use of Ajisai and Jason-1, together with the LAGEOS satellites, to measure the general relativistic Lense–Thirring effect in the gravitational field of the Earth. A linear combination of the nodes of such satellites is the proposed observable. The systematic error due to the mismodelling in the uncancelled even zonal harmonics would be ∼1% according to the latest present-day CHAMP/GRACE-based Earth gravity models. In regard to the non-gravitational perturbations especially affecting Jason-1, only relatively high-frequency harmonic perturbations should occur: neither semisecular nor secular bias of non-gravitational origin should affect the proposed combination: their maximum impact is evaluated to ∼4% over 2 years. Our estimation of the root-sum-square total error is about 4–5% over at least 3 years of data analysis required to average out the uncancelled tidal perturbations.     

利用双行根数计算太阳相位角的精度评估

光度特性测量是获取空间目标的物理特性的重要技术手段之一,无论是光变曲线的事后分析还是建立光度变化的仿真模型,都离不开一个重要的参数——太阳相位角（太阳-空间目标-测站的空间夹角）.目前空间目标的位置通常是通过双行根数（TLE）外推获得,存在一定误差,且随外推时间的延长而变大,因而有必要对其计算所得的太阳相位角的精度进行评估.以典型的不同高度的激光测距卫星LAGEOS1、AJISAI、STELLA为研究对象,以全球激光测距资料解算所得的高精度轨道作为参考轨道,对2012年全年利用双行根数计算所得的太阳相位角数据进行了比对分析,结果表明对于LAGEOS1、AJISAI这样的中高轨卫星,由于轨道较高,表征阻力的B＊恒定,计算所得的太阳相位角偏差较小,角分量级,且随外推时间的延长不会导致偏差明显增大;而对于STELLA这样的低轨卫星,因轨道较低、受变化的大气的影响显著,计算所得的太阳相位角偏差较大,尤其是当B＊比较大、变化较快时,偏差显著变大,且随外推时间的延长显著增大,在最差情况下：外推1d约为13＇,外推3d约为50＇,外推7d约为251＇,已超出目前的精度要求.因此,在事后分析中应尽可能使用1d之内的TLE计算太阳相位角,对于B＊较大且变化较快情况尤其需要注意.另外,针对UTC闰秒的情况,提出了一种处理方法,即在双行根数外推时判断外推时段是否跨越了闰秒时刻,若跨越了则进行修正：增加或减少1s,相应地需要修改结果对应的时间戳计算方法.     

The Parabolic Three-Body Problem

In this communication we present an analytical model for the restricted three-body problem, in the case where the perturber is in a parabolic orbit with respect to the central mass. The equations of motion are derived explicitly using the so-called Global Expansion of the disturbing function, and are valid for any eccentricity of the massless body, as well as in the case where both secondary masses have crossing orbits. Integrating the equations of motion over the complete passage of the perturber through the system, we are then able to construct a first-order algebraic mapping for the change in semimajor axis, eccentricity and inclination of the perturbed body.Comparisons with numerical solutions of the exact equations show that the map yields precise results, as long as the minimum distance between both bodies is not too small. Finally, we discuss several possible applications of this model, including the evolution of asteroidal satellites due to background bodies, and simulations of passing stars on extra-solar planets.     

Inversion of crater morphometric data to gain insight on the cratering process

Abstract— In recent years, morphometric data for Venus and several outer planet satellites have been collected, so we now have observational data of complex craters formed in a large range of target properties. We present general inversion techniques that can utilize the morphometric data to quantitatively test various models of complex crater formation. The morphometric data we use in this paper are depth of a complex crater, the diameter at which the depth-diameter ratio changes, and onset diameters for central peaks, terraces, and peak rings. We tested the roles of impactor velocities and hydrostatic pressure vs. crustal strength, and we tested the specific models of acoustic fluidization (Melosh, 1982) and nonproportional growth (Schultz, 1988). Neither the acoustic fluidization model nor the nonproportional growth in their published formulations are able to successfully reproduce the data. No dependence on impactor velocity is evident from our inversions. Most of the morphometric data is consistent with a linear dependence on the ratio of crustal strength to hydrostatic pressure on a planet, or the factor c/p g.     

Dynamical objectives of observation of mutual events

Recently the motion of the main satellites of Jupiter, Saturn and Uranus have been modelled in order to get accurate ephemerides. These models have been fitted over a large amount of observations. Among these ones, the positions issued from the observations of mutual events are the most accurate. We can then expect to obtain a new kind of dynamical informations directly linked to planetological questions. We have to determine what information is used in these observations to get the still unknown dynamical parameters. We look after these questions especially in the Jovian and Saturnian systems.     

Estimation of low degree geopotential coefficients using SLR data

Geodetic satellites have been providing the low frequency part of the geopotential models used for precise orbit determination purposes (e.g. JGM3, EGM96, …). Nevertheless they can be used to estimate the temporal variation of selected coefficients, helping to clarify the complex interrelations in the earth-ocean-atmosphere system. In this paper we present the two years long analysis of SLR data from the seven available geodetic satellites (Lageos I–II, Stella, Starlette, Ajisai, Etalon I–II) to recover monthly estimates of low degree geopotential coefficients; the results are obtained analysing the satellites separately and in proper combination. An accurate modelling of the satellite orbits is required in order to separate the geopotential coefficients: we assume as a priori geopotential the JGM3 model together with its associated tides and we take care of non-gravitational effects on the satellites by means of proper empirical estimated accelerations. The time series of the estimated coefficients (J₂, J₃, J₄, J₅) are inspected to detect the sub-annual perturbations related to seasonal variation of mass distribution. Huge residual seasonal signals in the orbit of Stella indicate a strong model deficiency related to the Sun's influence on the environment. The remaining six satellites are homogeneously modelled and build up a three cycles per year oscillation on J₂ and a seasonal oscillation (1 year and six month periods) revealed on the J₄. The origin and possible causes of these signals are further discussed in the text. We also present a preliminary estimate, using twelve years of Lageos-I and Lageos-II observations, that is compared with previous obtained values.     

SHORDE I program system and applications

As one of the Analysis Centers of SLR global data, a program system named SHORDE I has been developed at Shanghai Observatory. This software package is used for processing observational data of satellites and estimating satellite orbit and dynamical and geodetic parameters. The reference systems and the force model required for constructing the package were first realized at the cm-level. In order to save running time and the storage on the computer, a multistage-multiarc procedure has been especially designed.The technique has been applied to the SLR data from LAGEOS and earth rotation parameters and precise station coordinates were derived using the data from Nov. 1984 to Jan. 1988.     

A numerical solution for lageos thermal thrust: The rapid-spin case

We report the results of detailed numerical calculations of the thermal thrust on the rapidly-spinning LAGEOS spacecraft. This thrust results from anisotropic emission of thermal radiation from its surface. LAGEOS is a good test case for such calculations because of its relatively simple structure and because precise orbit determinations based on laser ranging give observed thrust effects for comparison.The numerical integration includes the varying heating over spacecraft-surface latitude from earth infrared radiation (for the earth-Yarkovsky force) and the varying solar heating as the spacecraft moves in and out of the earth's shadow (for the solar-Yarkovsky force). The computation allows for the poor thermal coupling between the spacecraft structure and individual surface elements (the fused-silica cube-corner reflectors and their aluminum retainer rings), and the poor conduction between structural hemispheres.A Fourier analysis of the computed force with respect to orbital longitude gives the important frequency components for the computation of long-term orbit perturbations. Empirical formulas fit to the numerical results accurately express the component amplitudes as simple functions of spin axis orbital latitude, the sun aspect angle from the spin axis, and the fraction of the orbit period spent in the earth's shadow. These results. based on first principles, are similar to those from simplified theories of the thermal thrust. but add the following new feature: The decrease in orbit-averaged satellite temperature when the orbit intersects the earth's shadow decreases the earth-Yarkovsky drag by 0.14 pm/s² from the no-eclipse value.The development of spacecraft-element thermal parameters is the most difficult part of the analysis; the paper tabulates the parameters that should be directly measured before the launch of future geodynamic satellites.     

Thermal force effects on slowly rotating, spherical artificial satellites-I. Solar heating

The heating of a spinning artificial satellite by natural radiation sources such as the Sun and the Earth results in temperature gradients arising across the satellite's surface. The corresponding anisotropic emission of thermal radiation leads to a recoil force, commonly referred to as “thermal force”. A quantitative theory of this effect is developed, based on more general assumptions than used so far, to model such radiation forces on spherically symmetric LAGEOS-like satellites. In particular, the theory holds for any ratio of the three basic timescales of the problem: the rotation period of the satellite, the orbital period around the Earth, and the relaxation time for the thermal processes. Thus, the simplifying assumption of a comparatively fast rotational motion is avoided, which will fail for LAGEOS within the next decade, owing to magnetic dissipation effects. A number of predictions about the future behaviour of non-gravitational long-term orbital perturbations of LAGEOS become possible with the new theory. In particular the Yarkovsky-Schach thermal force effects are studied arising as a consequence of the solar radiation flux onto the satellite, periodically interrupted by eclipses. Starting on about year 2005, the orbital perturbation effects predicted by the new theory are substantially different from those inferred in the fast-rotation case. This holds not only for the long-term semimajor axis effects, but also for eccentricity and inclination perturbations.     

Symmetric and Asymmetric Librations in Planetary and Satellite Systems at the 2/1 Resonance

Families of asymmetric periodic orbits at the 2/1 resonance are computed for different mass ratios. The existence of the asymmetric families depends on the ratio of the planetary (or satellite) masses. As models we used the Io-Europa system of the satellites of Jupiter for the case m₁>m₂, the system HD82943 for the new masses, for the case m₁=m₂ and the same system HD82943 for the values of the masses m₁<m₂ given in previous work. In the case m₁≥ m₂ there is a family of asymmetric orbits that bifurcates from a family of symmetric periodic orbits, but there exist also an asymmetric family that is independent of the symmetric families. In the case m₁<m₂ all the asymmetric families are independent from the symmetric families. In many cases the asymmetry, as measured by and by the mean anomaly M of the outer planet when the inner planet is at perihelion, is very large. The stability of these asymmetric families has been studied and it is found that there exist large regions in phase space where we have stable asymmetric librations. It is also shown that the asymmetry is a stabilizing factor. A shift from asymmetry to symmetry, other elements being the same, may destabilize the system.     

A semi-analytical model of the Galilean satellites’ dynamics

The Galilean satellites’ dynamics has been studied extensively during the last century. In the past it was common to use analytical expansions in order to get simple models to integrate, but with the new generation of computers it became prevalent the numerical integration of very sophisticated and almost complete equations of motion. In this article we aim to describe the resonant and secular motion of the Galilean satellites through a Hamiltonian, depending on the slow angles only, obtained with an analytical expansion of the perturbing functions and an averaging operation. In order to have a model as near as possible to the actual dynamics, we added perturbations and we considered terms that in similar studies of the past were neglected, such as the terms involving the inclinations and the Sun’s perturbation. Moreover, we added the tidal dissipation into the equations, in order to investigate how well the model captures the evolution of the system.