Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2009 / Article
Special Issue

Space Dynamics

View this Special Issue

Research Article | Open Access

Volume 2009 |Article ID 347835 |

Luiz Danilo Damasceno Ferreira, Rodolpho Vilhena de Moraes, "GPS Satellites Orbits: Resonance", Mathematical Problems in Engineering, vol. 2009, Article ID 347835, 12 pages, 2009.

GPS Satellites Orbits: Resonance

Academic Editor: Antonio Prado
Received29 Jul 2009
Accepted28 Sep 2009
Published22 Dec 2009


The effects of perturbations due to resonant geopotential harmonics on the semimajor axis of GPS satellites are analyzed. For some GPS satellites, secular perturbations of about 4?m/day can be obtained by numerical integration of the Lagrange planetary equations considering in the disturbing potential the main secular resonant coefficients. Amplitudes for long-period terms due to resonant coefficients are also exhibited for some hypothetical satellites orbiting in the neighborhood of the GPS satellites orbits. The results are important to perform orbital maneuvers of GPS satellites such that they stay in their nominal orbits. Also, for the GPS satellites that are not active, the long-period effects due to the resonance must be taken into account in the surveillance of the orbital evolutions of such debris.

1. Introduction

The period of the orbits of the GPS satellites is about 12 hours, and the main perturbations acting on their orbits are caused by the nonuniform distribution of the Earth’s mass, by the lunar and solar gravitational attractions and by the solar radiation pressure. In this paper, it is analyzed just some perturbations due to resonant terms of the geopotential coefficients. The resonance considered here is the 2 : 1 commensurability between the orbital period of the GPS satellites and the period of the Earth’s rotation. As it is pointed out by Hugentobler [1] the resonance leads a daily drift rate in semimajor axis of up to 7?m/day.

Lagrange planetary equations, describing the temporal variation of the orbital elements, are used here to analyze the orbital perturbations of the GPS satellites under the influence of resonant coefficients.

2. Disturbing Potential

The geopotential acting on Earth’s artificial satellite can be expressed as [2]

Here GM is the geogravitational constant, r is geocentric distance of the satellite, ae is the semimajor axis of the adopted Earth’s ellipsoid, a is the orbital semimajor axis, Cnm and Snm are the spherical harmonic coefficients. Pnm denotes the associated Legendre functions, are the satellite’s geocentric latitude and the longitude, n and m are, respectively, the degree and order of the harmonic coefficients.

The geocentric distance, the latitude, and the longitude of the satellite can be expressed in terms of the orbital elements and (2.1) becomes [2, 3] where and are, respectively, functions of the satellite’s orbital inclination i and eccentricity e, O represents the right ascension of the orbital ascending node, ? is the argument of perigee, M is the mean anomaly, and ?T is the Greenwich sidereal time. The function can be expressed as and the argument is given by

The functions and are presented as tables [2] and can be adapted to be used in computers.

3. Resonance in GPS Satellites

Resonance is associated with small divisors. For artificial Earth satellites whose orbital periods are in commensurability with the period of the Earth’s rotation resonance can occur when [1, 3, 4]

where is the Earth’s sidereal rotation.

For the GPS satellites where the commensurability is 2 : 1, there are bounds among the parameters

Taking into account that , (3.1) can be put as

since , we get

Considering harmonics of order and degree up to order 4, Table 1, presents some values for the parameters satisfying the resonance condition for the GPS satellites.

Degree n Order m pq


Taking into account in the summations the conditions and and putting [5]


equation (2.2) can be written in general form as [1]


In order to analyze the orbital perturbations of GPS satellites due to the resonant coefficients presented in (Table 1), the Lagrange planetary equations will be used. However, here this analysis will be concentrated initially in the secular effects of such perturbation on the orbital semimajor axis. Therefore, we have [2]

with being the mean motion of the satellite.

Special care must be taken to compute . In fact, the parameters given by Table 1 must be considered and each set of them gives a unique solution for (3.8).

A numerical integration of (3.8) was performed for a period of one day using IGS/POE (The International GNSS Service Precise Orbital Ephemeredes) given on June, 07, 2007. Transformations were performed to get the corresponding orbital elements.

Figures 1, 2, 3, and 4 show the daily variations of the semimajor axis for all GPS satellites.

Figure 1 gives the variation due to the coefficients , which are responsible by the greatest variation in the semimajor axis that, in this case, corresponds to the satellite PRN 06. Figures 2, 3, and 4 present, respectively, the daily variations as function of the coefficients . Figure 5 represents the total daily variation due to the above considered coefficients.

Table 2 shows the maxima and minima values for the daily variation of the semimajor axis according to the resonant coefficients.

Resonants coefficientsMaximum drift rates m/day


Figures 6 and 7 show the semimajor axis variation for the satellites PRN 02 and PRN 06 for a time interval of 200 days and taking into account the harmonic coefficients . During this period and for the considered initial condition, it can be observed variations of about 600?m and 680?m, respectively.

4. Long-Period Perturbations

The effects of the resonance are enhanced when long periods are considered. Table 3 and Table 4, Figures 8, 9, 10, 11, 12, and 13 show some simulations using hypothetical satellites and a particular method to study effects of resonance on the orbits of artificial satellites [69].

Orbital elementsAmplitudePeriod

(°) t(days)

+ 0.4300.0112 km0.0148.7 10-42500
+ 0.7180.0112.5 km0.01458.6 10-47000
+ 0.2900.014.5 km0.00653.6 10-42000
+ 0.5000.013 km0.00452.8 10-41000
− km0.00854 10-41700
−50.011.4 km0.00231.3 10-6500
0.0555°1.8 km0.0052 10-33200
+ 0.7180.00555°2.75 km0.0062.8 10-33900
−1.7700.00555°1.2 km0.00451.4 10-32000
+ 5.1290.00555°900 m0.003310-31500
+ 3.7290.00555°1.8 km0.00551.9 10-32500
0.0563.4°7.5 km0.00310.8 10-31000
−20.0563.4°4 km0.00155.7 10-31000
0.0587°7 km0.003713.9 10-32000
−20.0587°3.5 km0.00210.3 10-31900

Orbital elementsAmplitudePeriod

(°)t (days)

ao −3.770.01120 m0.00041.78 10-3500
ao + 0.7180.01245 m0.00093.9 10-31000
ao + 0.9290.01110 m0.000451.78 10-3500
ao + 50.01100 m0.00031.15 10-3250
ao −2.50.005100 m0.00094.6 10-5500
ao + 0.4290.005700 m0.0062.4 x 10-41500
ao + 0.9290.00555°1.9 km0.00856.8 10-416000
ao + 1.730.00555°250 m0.00171.13 10-316000
ao + 50.00555°30 m0.000552.3 10-41500
ao −3.770.00555°47 m0.000353.43 10-5500
ao −0.4290.00587°110 m0.000852.29 10-41800
ao −2.230.00587°45 m0.000358.59 10-5800
ao −3.770.00587°20 m0.000144 10-5500

Table 3 contains the amplitude and period of the variations of orbital elements for hypothetical satellites considering low eccentricity, small and high inclination, and the influence of the harmonics and

Figures 8 and 9 represent, respectively, the temporary variation of the semimajor axis and of the eccentricity of artificial satellites of the GPS type in the neighborhood of the 2 : 1 resonance region when the harmonics and are considered. By Figure 10, assuming several values for the semimajor axes in the neighborhood of the 2 : 1 resonance, it can be observed distinct behavior for their temporary variations.

Figure 11 represents the temporary variation for the semimajor axis in the neighborhood of the resonance 2 : 1 considering the influence of the harmonics and for a satellite of the GPS type with inclination of about .

Table 4 gives the influence of the resonance due to the harmonics and considering different inclinations, including that of the GPS type satellite. It can be observed that the amplitudes of the variations are smaller when compared with those when the harmonics and are taken into account. The more the satellites approach the region that was defined as a resonant region, the more the variations increase.

Figure 12 represents the temporary variation of the semimajor axis of a satellite of the GPS type considering resonance due to the harmonics , , and .

Figure 13 represents the time variation of the orbital semimajor axis of a satellite with eccentricity about 0.01 and inclination about orbiting in a region near the 2 : 1 resonance. It is remarkable the oscillation of the semimajor axis in the region between by a = 26560.0?km and a = 26562.48?km. For instance, a variation of about 10?m in the initial semimajor axis (a = 26562.48?km) brings up a variation of about more than 10?km in its amplitude in a period of about 2000 days.

5. Conclusions

From the obtained results it can be seen that for the GPS satellites no negligible perturbations are provoked by resonant tesseral coefficients. The daily variation of the semimajor axis due to the is less than 4?m/day and those values found for a time interval of 200 days (~680?m) are compatible with the results presented by specialized literature (~7?m/day and ~1200?m) about this subject [1].

Taking into account the total effects of daily resonant perturbations, it can be observed that for some GPS satellites these effects are enhanced and for another satellites are attenuated but all of these effects are smaller than the amplitudes mentioned above. It was shown also that the effects of the resonance are very important for the analysis of long-period behavior of the GPS satellites orbits.

An important aspect to be considered is the necessity to perform orbital maneuvers of GPS satellites in such way that they stay in their nominal orbits. Also, for the GPS satellites that are not active, the long-term effects due to the resonance must be taken into account in the surveillance of the orbital evolutions of such debris.


This work was partially supported by CNPq under the contracts No. 305147/2005-6 and 300952/2008-2


  1. U. Hugentobler, “Astrometry and satellite orbits: theoretical consideration and typical applications,” Geodätisch-geophysikalische Arbeiten in der Schweiz, vol. 57, 1998. View at: Google Scholar
  2. W. M. Kaula, Theory of Satellite Geodesy, Blaisdell, Waltham, Mass, USA, 1966.
  3. G. Seeber, Satellite Geodesy: Foundations, Methods and Applications, Walter de Gruyter, Berlin, Germany, 2003.
  4. D. Ineichen, G. Beutler, and U. Hugentobler, “Sensitivity of GPS and GLONASS orbits with respect to resonant geopotential parameters,” Journal of Geodesy, vol. 77, no. 7-8, pp. 478–486, 2003. View at: Publisher Site | Google Scholar
  5. D. Delikaraoglu, On Principles, Methods and Recent Advances in Studies Towards a GPS-Based Control System for Geodesy and Geodynamics, NASA Technical Memorandum 100716, Greenbelt, Md, USA, 1989.
  6. P. H. C. L. Lima Jr., Sistemas ressonantes a altas excentricidades no movimento de satélites artificiais, Doctoral thesis, ITA, São José dos Campos, 1998.
  7. P. H. C. L. Lima Jr., R. Vilhena de Moraes, and S. S. Fernandes, “Semi analytical method to study geopotential perturbations conssidering high eccentric resonant orbits,” in Dinamics of Natural and Artificial Celestial Bodies, H. Pretka-Ziomek, E. Wnuk, P. K. Seidelmann, and D. Richardson, Eds., pp. 407–413, Kluwer, Dordrecht, The Netherlands, 2001. View at: Google Scholar
  8. J. K. S. Formiga, Study of resonances in the orbital motion of artificial satellites, Master dissertation, Faculdade de Engenharia do Campus de Guaratinguetá, Universidade Estadual Paulista, Guaratinguetá, Brazil, 2005.
  9. R. Vilhena de Moraes, J. K. S. Formiga, P. H. C. L. Lima Jr., and H. K. Kuga, “Orbital perturbations: resonance effects,” in Proceedings of 6th International Symposium of the IAA on Small Satellites for Earth Observation, p. 8, Berlin, Germany, April 2007, paper IAA-B6-0717P. View at: Google Scholar

Copyright © 2009 Luiz Danilo Damasceno Ferreira and Rodolpho Vilhena de Moraes. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2021, as selected by our Chief Editors. Read the winning articles.