Evaluation of disturbances from solar radiation in orbital elements of geosynchronous satellites based on harmonics

Using specialized codes for the search of periodic and linear components we show that direct solar radiation leads to short-period variations of all the orbital elements of geosynchronous satellites. The variation period of the semimajor axis a, orbit inclination i and the longitude of the ascending node Ω is 1 day. Eccentricity e, the argument of perigee ω and the mean anomaly M vary with a period of 0.5 days. Direct solar radiation also leads to long-period variations in e, ω and M with a period of 1 year. The elements a, i and Ω undergo variations only in the amplitude of diurnal variations with a period of 1 or 0.5 years. Secular variability (linear components) are not detected. To obtain the initial value array of the orbital elements we used the Lagrange equations of perturbed motion in the form of a Gaussian with their subsequent integration via a special method of harmonics: the values of the derived orbital elements, obtained from the Lagrange equations, were presented through the periodic functions that are easy to integrate.