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.     

Short-period and long-period perturbations of a spherical satellite due to direct solar radiation

On the basis of expressions derived by Kozai, and new ones developed here, a detailed, semianalytic algorithm is presented for calculating radiation-pressure perturbations in the Keplerian elements. Through some simple modifications, the algorithm is also made to hold whene=0 and/ori=0. The perturbations are obtained by summing over the sunlit segment of the satellite's orbit during each revolution or part thereof. The end points of this segment are evaluated numerically once per revolution. The effect of the inherent uncertainties in the boundaries of the Earth's shadow is discussed. The algorithm is tested by means of numerical integration of the equations of motion and through comparisons with observations of the balloon satellite 1963 30D during a 200-day interval.     

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.     

Lunar perturbations of artificial satellites of the earth

Lunisolar perturbations of an artificial satellite for general terms of the disturbing function were derived by Kaula (1962). However, his formulas use equatorial elements for the Moon and do not give a definite algorithm for computational procedures. As Kozai (1966, 1973) noted, both inclination and node of the Moon's orbit with respect to the equator of the Earth are not simple functions of time, while the same elements with respect to the ecliptic are well approximated by a constant and a linear function of time, respectively. In the present work, we obtain the disturbing function for the Lunar perturbations using ecliptic elements for the Moon and equatorial elements for the satellite. Secular, long-period, and short-period perturbations are then computed, with the expressions kept in closed form in both inclination and eccentricity of the satellite. Alternative expressions for short-period perturbations of high satellites are also given, assuming small values of the eccentricity. The Moon's position is specified by the inclination, node, argument of perigee, true (or mean) longitude, and its radius vector from the center of the Earth. We can then apply the results to numerical integration by using coordinates of the Moon from ephemeris tapes or to analytical representation by using results from lunar theory, with the Moon's motion represented by a precessing and rotating elliptical orbit.     

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.     

The effects of diffuse solar radiation on the motion of an artificial satellite

Previous analyses of the orbits of spherical balloon satellites have attempted to satisfy residuals in observed perturbations in Keplerian elements, assumed to be caused by diffuse radiation pressure, by introducing small variations ins, the parameter representing the reflection characteristics of the satellite's surface. It is difficult to distinguish, however, between those perturbations caused by diffuse radiation and those caused by reflected radiation, as a result of the deformation of the assumed sphere. Following the derivation by Lucas of exact expressions for both incident and reflected radiation forces on a prolate spheroidal satellite, and the subsequent work of Aksnes pertaining to spherical satellites, the theory is extended to include the effects of diffuse radiation whilst at the same time qualifying the assumption that the radiation force acting along the Sun-satellite line can be taken as parallel to the Sun-Earth line.     

Analytical short-term orbit predictions withJ 3 andJ 4 in terms of KS elements

Analytical theory for short-term orbit motion of satellite orbits with Earth's zonal harmonicsJ ₃ andJ ₄ is developed in terms of KS elements. Due to symmetry in KS element equations, only two of the nine equations are integrated analytically. The series expansions include terms of third power in the eccentricity. Numerical studies with two test cases reveal that orbital elements obtained from the analytical expressions match quite well with numerically integrated values during a revolution. Typically for an orbit with perigee height, eccentricity and inclination of 421.9 km, 0.17524 and 30 degrees, respectively, maximum differences of 27 and 25 cm in semimajor axis computation are noted withJ ₃ andJ ₄ term during a revolution. For application purposes, the analytical solutions can be used for accurate onboard computation of state vector in navigation and guidance packages.     

The lunar orbit revisited,III

In this paper we present an investigation on the tidal evolution of a system of three bodies: the Earth, the Moon and the Sun. Equations are derived including dissipation in the planet caused by the tidal interaction between the planet and the satellite and between the planet and the sun. Dissipation within the Moon is included as well. The set of differential equations obtained is valid as long as the solar disturbances dominate the perturbations on the satellite's motion due to the oblateness of the planet, namelya/R _e greater than 15, and closer than that point equations derived in a preceding paper are used.The result shows the Moon was closer to the Earth in the past than now with an inclination to the ecliptic greater than today, whereas the obliquity was smaller. Toward the past, the inclination to the Earth's equator begins decreasing to 12° fora/R _e=12 and suddenly grows. During the first stage the results are weakly dependant on the magnitude of the dissipation within the satellite, whereas the distance of the closest approach and the prior history are strongly dependent on that dissipation. In particular, the crossing of the Roche limit can be avoided.     

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.     

Radial,transverse and normal satellite position perturbations due to the geopotential

Perturbations in the position of a satellite due to the Earth's gravitational effects are presented. The perturbations are given in the radial, transverse (or alongtrack) and normal (or cross-track) components. The solution is obtained by projecting the Kepler element perturbations obtained by Kaula [Kaula, 1966] into each of the three components. The resulting perturbations are presented in a form analogous to the form of Kaula's solution which facilitates implementation and interpretation.     

The motion of a satellite in resonance with the second-degree sectorial harmonic

The solution to the motion of a satellite in an eccentric orbit and in resonance with the second-degree sectorial harmonic of the potential field is developed. The method of solution used parallels the well known von Zeipel method of general perturbations. The solution consists of expressions for the variations of the Delaunay variables. These expressions are composed of the perturbations developed by Brouwer in 1959 for the motion of an artificial satellite plus first-order perturbations due to the second-degree sectorial harmonic (in terms of the Legendre normal elliptic integrals of the first and second kind).This paper presents the results of one phase of research carried out at the Jet Propulsion Laboratory, California Institute of Technology, under Contract No. NAS 7-100, sponsored by the National Aeronautics and Space Administration.     

Applying an Expansion of the Second Derivatives of the Earth’s Gravitational Potential in Modified Spherical Harmonics for Constructing and Estimating Its Models

Second-order derivatives of the Earth’s potential in a local north-oriented coordinate system are expanded in series of modified spherical harmonics. Linear relations are derived between the spectral coefficients of these series and the spectrum of the geopotential. Based on these relations, recurrent procedures are developed for estimating the geopotential coefficients from the spectrum of each derivative and, conversely, for simulating the spectrum from a known geopotential model. The very simple structure of the expressions for the derivatives is convenient for estimating the coefficients of the geopotential by the least squares method at a certain step of processing satellite gradiometry data. Since the new series are orthogonal, the method with a quadrature formula can be applied, which helps avoid aliasing errors caused by the truncation of the series. The spectral coefficients of the derivatives are estimated using the derived relations for different models on an average orbital sphere of the GOCE satellite and at other altitudes above the Earth’s surface.

     

Estimation of lunar and solar tides in the Earth’s thermosphere from density measurements by the CACTUS microaccelerometer

The densities measured by the CACTUS microaccelerometer at altitudes from 270 to 600 km are used to analyze the effect of tidal perturbations in the Earth’s thermosphere caused by the gravitational attraction of the Moon and the Sun. These tidal perturbations are considered a priori small and are not taken into account in modern atmospheric density models. The residuals between the densities measured by the CACTUS microaccelerometer and calculated by models are analyzed, and the density variations correlating with variations of the zenith angles from the Moon to the center of the Earth to the satellite and from the Sun to the center of the Earth to the satellite are found at altitudes from 270 to 600 km. The amplitude of the perturbations revealed in the study grows with height. The phase of the tidal perturbations also varies with height. The amplitude of the density variations is about 30% at 270–320 km and increases to 80% at 520–570 km. The results agree with a priori theoretical estimates obtained for tidal motion of gaseous matter with a variable density.     

Spectral decomposition of geopotential,earth and ocean tidal perturbations in linear satellite theory

The spectra of geopotential, Earth and ocean tidal perturbations on a satellite can be obtained using Kaula's linear theory, or an extension thereof, as summations of terms depending on four indices l, m, p, q. In this work algorithms are presented that generate the equivalence classes induced by the composition rule of frequency on the set of all (l, m, p, q) combinations up to a maximum degree L and maximum value Q of the last index. These algorithms eliminate the need to search the set of frequencies when the linear theory is programmed on a computer.     

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.     

Precision Analytical Calculation of Geodynamical Effects on Satellite Motion

A new analytical method for calculating satellite orbital perturbations due to different disturbing forces is developed. It is based on the Poincaré method of small parameter but takes advantages of modern high-performance computers and of the tools of computer algebra. All perturbations proportional up to and including the 5th-order of small parameters are obtained. The method can precisely calculate the effects of all geodynamical forces on satellite motion given by the most up-to-date IAU and IERS models, such as non-central Earth gravity potential, precession and nutation of the geoequator, polar motion and irregularities in the Earth's rotation, effect of ocean and solid Earth tides, pole tide, and secular variations of gravity coefficients.Numerical tests prove the method's accuracy to be equivalent to 1–2 cm when calculating positions of high altitude geodetic satellites (like ETALON), and/or of GLONASS navigational spacecraft. The accuracy is stable over 1 year at least and comparable to that of the best tracking measurements of satellites.Positions of low altitude geodynamical satellites (like STARLETTE) by the analytical method are calculated to an accuracy of about 70cm over a month's interval. The method is developed for future use in GLONASS/GPS on-board ephemeris computation where it can improve the current scheme of their flight control.This revised version was published online in October 2005 with corrections to the Cover Date.     

On the determination of the long period tidal perturbations in the elements of artificial earth satellites

In the present article we develop the theory of the long period tidal effects in the motion of artificial satellites assuming the variability of elastic parameters of the Earth (Love numbers) across the parallels. The dependence of Love numbers on the longitude produces perturbations of the period of one day or less and hence is neglected in the present theory. In this respect we follow in the footsteps of Kaula (1969). If the deviations ofk ₂ andk ₃ from pure constants are not taken into consideration, then the perturbations caused by the variability ofk ₂ andk ₃ across the parallels will be misinterpreted as the perturbations caused byk ₄...-terms, and the spurious values ofk ₄... will be deduced. It is extremely doubtful, however, that the real effects caused byk ₄,k ₅,..., are significant enough to be detected. The short period effects with the period of the revolution of the satellite, or less, were removed from the differential equations for the variation of elements of the satellite by the averaging over the orbit of the satellite. These differential equations are in the form convenient for numerical integration over a long interval of time and also suitable for developing the tidal effects into trigonometric series with the arguments ω, Ω of the satellite andl, l′, F, D, Γ of the Moon. The numerical integration can be performed using some simple quadrature formula, without resorting to a predictor-corrector system.     

Tesseral harmonic perturbations in radial transverse and binormal components

The formulae for the perturbations in radial, transverse and binormal components of the Earth artificial satellite motion have been derived. Perturbations due to the tesseral part of the geopotential are considered. The geopotential expressed in terms of the orbital elements has the form proposed by Wnuk (1988). The formulae for the perturbations have been obtained using the Hori (1966) method. They can be effectively applied in calculation of the perturbations in the components including the coefficients of the high order and degree tesseral harmonics. The derived formulae reveal no singularities at zero eccentricity.     

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.     

Analytical method for the effects of the asteroid belt on planetary orbits

Analytic expressions are derived for the perturbation of planetary orbits due to a thick constant density asteroid belt. The derivations include extensions and adaptations of Plakhov's analytic expressions for the perturbations in five of the orbital elements for closed orbits around Saturn's rings. The equations of Plakhov are modified to include the effect of ring thickness and additional equations are derived for the perturbations in the sixth orbital element, the mean anomaly. The gravitational potential and orbital perturbations are derived for the asteroid belt with and without thickness, and for a hoop approximation to the belt. The procedures are also applicable to Saturn's rings and the newly discovered rings of Uranus.The effects of the asteroid belt thickness on the gravitational potential coefficients and the orbital motions are demonstrated. Comparisons between the Mars orbital perturbations obtained using the analytic expressions and those obtained using numerical integration are discussed. The effects of the asteroid belt on the Earth based ranging to Mars are also demonstrated.