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.     

On the tidal effects in the motion of artificial satellites

The general perturbations in the elliptic and vectorial elements of a satellite as caused by the tidal deformations of the non-spherical Earth are developed into trigonometric series in the standard ecliptical arguments of Hill-Brown lunar theory and in the equatorinal elements ω and Ω of the satellite. The integration of the differential equations for variation of elements of the satellite in this theory is easy because all arguments are linear or nearly linear in time. The trigonometrical expansion permits a judgment about the relative significance of the amplitudes and periods of different tidal ‘waves’ over a long period of time. Graphs are presented of the tidal perturbations in the elliptic elements of the BE-C satellite which illustrate long term periodic behavior. The tidal effects are clearly noticeable in the observations and their comparison with the theory permits improvement of the ‘global’ Love numbers for the Earth.     

Sun-synchronous orbits near critical inclination

The long period dynamics of Sun-synchronous orbits near the critical inclination 116.6° are investigated. It is known that, at the critical inclination, the average perigee location is unchanged by Earth oblateness. For certain values of semimajor axis and eccentricity, orbit plane precession caused by Earth oblateness is synchronous with the mean orbital motion of the apparent Sun (a Sun-synchronism). Sun-synchronous orbits have been used extensively in meteorological and remote sensing satellite missions. Gravitational perturbations arising from an aspherical Earth, the Moon, and the Sun cause long period fluctuations in the mean argument of perigee, eccentricity, inclination, and ascending node. Double resonance occurs because slow oscillations in the perigee and Sun-referenced ascending node are coupled through the solar gravity gradient. It is shown that the total number and infinitesimal stability of equilibrium solutions can change abruptly over the Sun-synchronous range of semimajor axis values (1.54 to 1.70 Earth radii). The effect of direct solar radiation pressure upon certain stable equilibria is investigated.     

Perturbations of the orbit of Explorer 19 due to solar radiation

The orbit of the balloon satellite, Explorer 19, is analysed to determine the effects of direct solar radiation pressure over one revolution of the satellite (111 min) for MJD 42822 and MJD 42966. At the earlier date, the satellite entered the Earth's shadow, presenting an opportunity to examine the effectiveness of two different shadow models. The reflectivity of the surface of the satellite was estimated from analysis of the variations in orbital eccentricity over a period of 236 days.Although many of the parameters associated with the shape and orientation of the satellite are unknown, the theory for a non-spherical satellite is applied using trial and error methods to determine the parameters of best fit. The paper concludes with an examination of the perturbations in orbital eccentricity and inclination due to incident, specularly reflected, and diffusely reflected radiation.     

Exospheric zonal winds between 540 and 620 km from the orbit of Explorer 24

Data on the variation of the orbital inclination of the balloon satellite Explorer 24 (1964-76A) from 1964 to 1968 have been used to determine zonal winds between 540 and 620 km. In this height region the effect of zonal winds on the orbital inclination may become very small compared to other perturbations like accelerations due to the geopotential, lunisolar gravitation and the solar radiation pressure. It is demonstrated especially that the solar radiation pressure may become the most significant force changing the orbital inclination. The diurnal mean zonal winds derived from Explorer 24 point to an exospheric rotation rate which is about 6–10% less than the rotation rate of the Earth in the analyzed height region. Since the possible errors of the data analysis are of a similar order of magnitude, it can not be excluded that the exosphere corotates with the Earth. Furthermore, a local time dependence of the zonal winds could be detected. The diurnal varitation of the zonal wind is shown to be in good agreement with the theoretical model of Blum and Harris. Our results are discussed and compared with all previous investigations of orbital inclination changes of satellites above 350 km.     

Variations of the thermal accommodation coefficient from the analysis of the atmospheric lift effects on the satellite 1974-70 A

The changes of the orbital inclination of the satellite 1974-70 A show some peculiarities which cannot be explained by the usual disturbing effects (odd zonal harmonics, lunisolar perturbations, rotation of the atmosphere). The effect of a lift force normal to the orbital plane explains the residual inclination variations if a changing value of the thermal accommodation coefficient is introduced. The theory agrees well with the observations if the coefficient rises from a low value (～0.1) corresponding to a quasi-specular reflection of the molecules incident to the satellite surface to higher values (>0.9) characterizing a diffuse reflection.     

Long-periodic and secular perturbations to the orbit of a spherical satellite due to direct solar radiation pressure

A theory is developed for the perturbations to the orbit of a spherically symmetric satellite which accounts for the changes in the perigee and nodal positions and the variations of the Sun-Earth distance and direction over an orbital revolution. The theory is semi-analytical, the equations of motion being integrated with respect to time over the sunlit period of each orbital revolution. Long-periodic and short-periodic perturbations may be treated separately, and this is important for long-term analyses in terms of mean elements where short-period terms are averaged or omitted.     

Albedo perturbation models: General formalism and applications to LAGEOS

The force due to radiation pressure on a satellite of arbitrary shape is written in a general form within a formalism similar to that used in the theory of radiative transfer in atmospheres. Then the corresponding integrals are evaluated for the simple case of a spherically symmetric satellite, and applied to model the perturbation due to the Earth-reflected radiation flux on LAGEOS. For this purpose, the optical behaviour of the Earth's surface and atmosphere is described as a combination of Lambertian diffusion (continents), partial specular reflection consistent with Fresnel law (oceans) and anisotropic diffusion according to Chandrasekhar's radiative transfer theory (clouds). The in-plane Gauss componentsT andS vs. mean anomaly are computed for a simple orbital geometry and for different models of the Earth's optical properties. A sensitive dependence is found on the assumed cloud distribution, with significant perturbations possibly arising from oceanic specular reflection when the satellite is close to the Earth's shadow boundaries.On leave from Astronomical Institute, Charles University, védská 8, 15000 Prague 5, Czechoslovakia     

Perturbations of a spheroidal satellite due to direct solar radiation pressure

The orbital accelerations of certain balloon satellites exhibit marked oscillations caused by solar radiation impinging on the surface of the satellites, which, once spherical, have assumed a spheroidal shape producing a component of force at right-angles to the Sun-satellite direction. Given the characteristics and orientation of the satellite, the equations of force are determined by the formulae of Lucas. Otherwise the phase-angle and magnitude of the right-angle force are determined by trial and error, or best-fit techniques. Using a variation of the approach developed by Aksnes, a semi-analytical algorithm is presented for evaluating the perturbations of the Keplerian elements by direct solar radiation pressure on a spheroidal satellite. The perturbations are obtained by summing over the sunlit part of each orbit and allow for a linear variation in the phase-angle. The algorithm is used to determine the orbital accelerations of 1963-30D due to direct solar radiation pressure, and these results are compared to the observed values over two separate periods of the satellite's lifetime.     

Perturbations of a close-Earth satellite due to sunlight diffusely reflected from the Earth

A semianalytic method has been developed to calculate the radiation-pressure perturbations of a close-Earth satellite due to sunlight reflected from the Earth. The assumptions made are that the satellite is spherically symmetric and that the solar radiation is reflected from the Earth according to Lambert's Law with uniform albedo. By using expressions for the components of the radiation-pressure force due to Lochry, the expressions for the perturbations of the elements were developed into series in the true anomalyv. The perturbations within a given revolution can be obtained analytically by integrating with respect tov while holding all slowly varying quantities constant. The long-range perturbations are then obtained by accumulating the net perturbations at the end of each revolution.     

Orbital perturbations due to radiation pressure for a spacecraft of complex shape

We analyze the perturbations due to solar radiation pressure on the orbit of a high artificial satellite. The latter is modelled in a simplified way (axisymmetric body plus despun antenna emitting a radio beam), which seems suitable to describe the main effects for existing telecommunication satellites. We use the regularized general perturbation equations, by expressing the force in the moving Gauss' reference frame and by expanding the results in terms of some small parameters, referring both to the orbit (small eccentricity and inclination) and to the spacecraft's attitude. Some interesting results are derived, which assess the relative importance of different physical effects and of different parts of the spacecraft in determining the long-term evolution of the orbital elements.     

Perturbations in satellite orbits due to variations in rotation of the central body

Variations in satellite orbital elements are derived due to perturbations in the external gravitational field of the central body caused by mass deformations of the body occurring from variations in its rotation; the central body is assumed to be perfectly elastic. General theory derived is applied to the actual Earth, as an example; possible resonances are discussed.     

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.     

Geostationary secular dynamics revisited: application to high area-to-mass ratio objects

The long-term dynamics of the geostationary Earth orbits (GEO) is revisited through the application of canonical perturbation theory. We consider a Hamiltonian model accounting for all major perturbations: geopotential at order and degree two, lunisolar perturbations with a realistic model for the Sun and Moon orbits, and solar radiation pressure. The long-term dynamics of the GEO region has been studied both numerically and analytically, in view of the relevance of such studies to the issue of space debris or to the disposal of GEO satellites. Past studies focused on the orbital evolution of objects around a nominal solution, hereafter called the forced equilibrium solution, which shows a particularly strong dependence on the area-to-mass ratio. Here, we (i) give theoretical estimates for the long-term behavior of such orbits, and (ii) we examine the nature of the forced equilibrium itself. In the lowest approximation, the forced equilibrium implies motion with a constant non-zero average ‘forced eccentricity’, as well as a constant non-zero average inclination, otherwise known in satellite dynamics as the inclination of the invariant ‘Laplace plane’. Using a higher order normal form, we demonstrate that this equilibrium actually represents not a point in phase space, but a trajectory taking place on a lower-dimensional torus. We give analytical expressions for this special trajectory, and we compare our results to those found by numerical orbit propagation. We finally discuss the use of proper elements, i.e., approximate integrals of motion for the GEO orbits.     

Analytical solutions to the four post-Newtonian effects in a near-Earth satellite orbit

On the basis of the results by Huang et al. (1990), this paper further discusses and analyses the four post-Newtonian effects in a near-Earth satellite orbit: the Schwarzschild solution, the post-Newtonian effects of the geodesic precession, the Lense-Thirring precession and the oblateness of the Earth. A full analytical solution to the effects including their direct perturbations and mixed perturbations due to the Newtonian oblateness (J ₂) perturbation and the Schwarzschild solution is obtained using the quasi-mean orbital element method analogous to the Kozai's mean orbital element one. Some perturbation properties of the post-Newtonian effects are revealed. The results obtained not only can provide a sound scientific basis for the precise determination of a man-made satellite orbit but also is suitable for similar mechanics systems, such as the motions of planets, asteroids and natural satellites.     

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.     

Perturbations of a close-earth satellite due to sunlight diffusely reflected from the Earth

A semianalytical method has been developed to calculate the radiation-pressure perturbations of a close-Earth satellite due to sunlight reflected from the Earth. It is assumed that the satellite is spherically symmetric and that the solar radiation is reflected from the Earth according to Lambert's Law. To account for the increasing reflectivity of the Earth toward the poles, its albedo is assumed to have a latitudinal dependence given bya=a ₀ +a ₂ sin². The effect of the terminator on the perturbations has been neglected. The perturbations within a particular revolution are given analytically, while the long-range perturbations are obtained by accumulation.     

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.     

The Artificial Earth Satellite Motion in an Oblate,Rotating Atmosphere with Symmetrical Diurnal Effect

The equations describing the disturbed motion of an artificial Earth satellite in the atmosphere are integrated by using a new variable instead of the true anomaly. The atmospheric flattening and rotation, the linear variation of the density scale height with the altitude, the symmetrical diurnal effect and the variation of its amplitude with the heigth are taken into account. Approximate formulae for the perturbations in all the orbital elements over a revolution of the satellite are given.     

Luni-solar perturbations of an Earth satellite

Luni-solar perturbations of the orbit of an artificial Earth satellite are given by modifying the analytical theory of an artificial lunar satellite derived by the author in recent papers. Expressions for the first-order changes, both secular and periodic, in the elements of the geocentric Keplerian orbit of the earth satellite are given, the moon's geocentric orbit, including solar perturbations in it, being found by using Brown's lunar theory.The effects of Sun and Moon on the satellite orbit are described to a high order of accuracy so that the theory may be used for distant earth satellites.