Рассмотрена динамическая эволюция космического мусора в окрестности активно используемых орбит: геостационарной и глобальных навигационных спутниковых систем ГЛОНАСС, GPS, Бэйдоу и Галилео, с целью оценить интервалы времени, необходимые для проникновения объектов космического мусора в область движения геостационарных и навигационных спутников. Начальные значения больших полуосей орбит выбирались на 500~км выше и ниже номинальных орбит навигационных спутников в окрестности резонансов высоких порядков, также были рассмотрены супергеосинхронные орбиты в окрестности резонансов высоких порядков. Отношение миделева сечения к массе варьировалось от малых значений, соответствующих спутникам, до больших значений, соответствующих фрагментам космического мусора. Орбитальная эволюция объектов моделировалась на интервале времени 240 лет с помощью "Численной модели движения искусственных спутников Земли", разработанной в Томском государственном университете. Модель возмущающих сил учитывала: гравитационное поле Земли, притяжение Луны и Солнца, приливы в теле Земли, световое давление с учетом тени Земли, эффект Пойнтинга - Робертсона, сопротивление атмосферы. Уравнения движения интегрировались методом Эверхарта 19 порядка. Результаты существенно зависят от положения плоскости орбиты. Большие значения отношения миделева сечения к массе приводят к прохождению объектов через геостационарную область и область движения спутников глобальных навигационных систем вследствие светового давления и эффекта Лидова - Кодзаи, приводящих к долгопериодическим колебаниям эксцентриситета и наклона орбит.

Dynamical evolution of space debris was considered in the vicinity of actively used orbits: geostationary and global navigation satellite systems GLONASS, GPS, BeiDou, Galileo. The aim of the research is to estimate time intervals which are necessary for space debris objects to get into the region of navigation satellites motion and the geostationary orbit. Initial orbits of objects were chosen 500 km above and below with respect to nominal semi-major axes of navigation orbits and on the resonant super-geostationary orbits. The area-to-mass ratio was varied from small values corresponding to satellites to high values which correspond to the space debris fragments. The orbital evolution of objects was modeled for 240 years using “Numerical Model of Motion of Artificial Satellites” elaborated at the Tomsk State University. The model of perturbing forces takes into account the major perturbing factors: gravitational field of the Earth, gravitation of the Moon and the Sun, tides of the Earth, direct radiation pressure by taking into consideration the Earth’s shadow, the Poynting--Robertson effect, and the atmospheric drag. The equations of motion were integrated by Everhart's method of the 19 th order. Results significantly depend on the orbital plane position. High values of the area-to-mass ratio cause the objects to pass through the geostationary region and the region of navigation satellites motion due to the direct radiation pressure and the Lidov-Kozai effect, which produce long-period oscillations of eccentricity and orbit inclination.