Evolution of the orbital elements for objects with high area-to-mass ratios in geostationary transfer orbits

A new population of uncatalogued objects in geosynchronous Earth orbits (GEO), with a mean motion of about 1 rev/day and eccentricities up to 0.6, has been identified recently. The first observations of this new type of objects were acquired in the framework of the European Space Agency’s (ESA) search for space debris in GEO and the geostationary transfer orbit (GTO) using the ESA 1-m telescope on Tenerife. Earlier studies have postulated that the perturbations due to the solar radiation pressure can lead to such large eccentricities for GEO objects with a high area-to-mass ratio (A/M). The simulations showed that the eccentricities of GEO objects with large A/M exhibit periodic variations with periods of about one year and amplitudes depending on the value of A/M. The findings of these studies could be confirmed by observations from the ESA 1-m telescope on Tenerife.     

Orbit improvement for GTO objects using follow-up observations

The Astronomical Institute of the University of Bern was and is conducting several search campaigns for space debris in the geostationary ring (GEO) and the geostationary transfer orbit (GTO). First tests to build up a catalogue of satellites and space debris were performed recently. As the discovery observations of an object usually cover a time interval of a few minutes only, the first orbits determined are assumed to be circular. For GEO objects a circular orbit is a rather good approximation. This is certainly not the case for GTO objects. Nevertheless, a circular orbit seems to be an acceptable approximation for a short time interval around the observations. The first orbits can be used to acquire follow-up observations, which allow the determination of elliptical orbits. The maximum allowed time span between the detection and the follow-up for a successful recovery is found using simulated GTO orbits and observations. Further follow-up observations are simulated in order to study the orbit improvement process and the required force model. The accuracy of orbits needed to build up a catalogue is studied. The simulations are compared to the results achieved with observations acquired by the 1 m telescope on Tenerife.     

The situation in the geostationary ring

For more than 25 years, the practice of reorbiting of a geostationary satellite at the end of its mission in order to protect the geostationary orbit (GEO) environment has been recommended and performed by a number of operators. In recent years, an internationally recognised re-orbiting altitude has been defined by the Inter-Agency Space Debris Coordination Committee (IADC). Based on orbital data contained in the DISCOS database, the situation on the geostationary ring is analysed. In January 2004, from 1036 catalogued objects passing through the geostationary region, 340 are controlled within their allocated longitude slots, 395 are drifting above, below or through GEO, and 140 are in a libration orbit. In the periods 1997–2003, 103 spacecraft reached their end of life; 34 were reorbited in compliance with the IADC recommendation, 35 were reorbited below the minimum recommended altitude, and 34 were abandoned without any end-of-life disposal manoeuvre. Apart from these catalogued objects, the ESA 1-m telescope has observed many smaller debris (down to 10–15 cm) in this orbital region representing a collision risk for GEO spacecraft which is difficult to quantify.     

Apparent rotation properties of space debris extracted from photometric measurements

Knowledge about the rotation properties of space debris objects is essential for the active debris removal missions, accurate re-entry predictions and to investigate the long-term effects of the space environment on the attitude motion change. Different orbital regions and object’s physical properties lead to different attitude states and their change over time.Since 2007 the Astronomical Institute of the University of Bern (AIUB) performs photometric measurements of space debris objects. To June 2016 almost 2000 light curves of more than 400 individual objects have been acquired and processed. These objects are situated in all orbital regions, from low Earth orbit (LEO), via global navigation systems orbits and high eccentricity orbit (HEO), to geosynchronous Earth orbit (GEO). All types of objects were observed including the non-functional spacecraft, rocket bodies, fragmentation debris and uncorrelated objects discovered during dedicated surveys. For data acquisition, we used the 1-meter Zimmerwald Laser and Astrometry Telescope (ZIMLAT) at the Swiss Optical Ground Station and Geodynamics Observatory Zimmerwald, Switzerland. We applied our own method of phase-diagram reconstruction to extract the apparent rotation period from the light curve. Presented is the AIUB’s light curve database and the obtained rotation properties of space debris as a function of object type and orbit.     

International scientific optical network for space debris research

A joint team of researchers under the auspices of the Center for Space Debris Information Collection, Processing and Analysis of the Russian Academy of Sciences collaborates with 15 observatories around the world to perform observations of space debris. For this purpose, 14 telescopes were equipped with charge-coupled device (CCD) cameras, Global Positioning System (GPS) receivers, CCD frame processing and ephemeris computation software, with the support of the European and Russian grants. Many of the observation campaigns were carried out in collaboration with the Astronomical Institute of the University of Bern (AIUB) team operating at the Zimmerwald observatory and conducting research for the European Space Agency (ESA), using the Tenerife/Teide telescope for searching and tracking of unknown objects in the geostationary region (GEO). More than 130,000 measurements of space objects along a GEO arc of 340.9°, collected and processed at Space Debris Data Base in the Ballistic Center of the Keldysh Institute of Applied Mathematics (KIAM) in 2005–2006, allowed us to find 288 GEO objects that are absent in the public orbital databases and to determine their orbital elements. Methods of discovering and tracking small space debris fragments at high orbits were developed and tested. About 40 of 150 detected unknown objects of magnitudes 15–20.5 were tracked during many months. A series of dedicated 22-cm telescopes with large field of view for GEO survey tasks is in process of construction. 7 60-cm telescopes will be modernized in 2007.     

Feasibility of performing space surveillance tasks with a proposed space-based optical architecture

Under ESA contract an industrial consortium including Aboa Space Research Oy (ASRO), the Astronomical Institute of the University of Bern (AIUB), and the Dutch National Aerospace Laboratory (NLR), proposed the observation concept, developed a suitable sensor architecture, and assessed the performance of a space-based optical (SBO) telescope in 2005. The goal of the SBO study was to analyse how the existing knowledge gap in the space debris population in the millimetre and centimetre regime may be closed by means of a passive optical instrument. The SBO instrument was requested to provide statistical information on the space debris population in terms of number of objects and size distribution. The SBO instrument was considered to be a cost-efficient with 20 cm aperture and 6° field-of-view and having flexible integration requirements. It should be possible to integrate the SBO instrument easily as a secondary payload on satellites launched into low-Earth orbits (LEO), or into geostationary orbit (GEO). Thus the selected mission concept only allowed for fix-mounted telescopes, and the pointing direction could be requested freely. Since 2007 ESA focuses space surveillance and tracking activities in the Space Situational Awareness (SSA) preparatory program. Ground-based radars and optical telescopes are studied for the build-up and maintenance of a catalogue of objects. In this paper we analyse how the proposed SBO architecture could contribute to the space surveillance tasks survey and tracking. We assume that the SBO instrumentation is placed into a circular sun-synchronous orbit at 800 km altitude. We discuss the observation conditions of objects at higher altitude, and select an orbit close to the terminator plane. A pointing of the sensor orthogonal to the orbital plane with optimal elevation slightly in positive direction (0° and +5°) is found optimal for accessing the entire GEO regime within one day, implying a very good coverage of controlled objects in GEO, too. Simulations using ESA’s Program for Radar and Optical Observation Forecasting (PROOF) in the version 2005 and a GEO reference population extracted from DISCOS revealed that the proposed pointing scenario provides low phase angles together with low angular velocities of the objects crossing the field-of-view. Radiometric simulations show that the optimal exposure time is 1–2 s, and that spherical objects in GEO with a diameter of below 1 m can be detected. The GEO population can be covered under proper illumination nearly completely, but seasonal drops of the coverage are possible. Subsequent observations of objects are on average at least every 1.5 days, not exceeding 3 days at maximum. A single observation arc spans 3° to 5° on average. Using a simulation environment that connects PROOF to AIUB’s program system CelMech we verify the consistency of the initial orbit determination for five selected test objects on subsequent days as a function of realistic astrometric noise levels. The initial orbit determination is possible. We define requirements for a correlator process essential for catalogue build-up and maintenance. Each single observation should provide an astrometric accuracy of at least 1”–1.5” so that the initially determined orbits are consistent within a few hundred kilometres for the semi-major axis, 0.01 for the eccentricity, and 0.1° for the inclination.     

Collision probability and spacecraft disposition in the geostationary orbit

Since 1963 approximately 300 satellites have been launched into the geostationary orbit, followed possibly by another additional 200 satellites up to the year 2000. Ground surveillance with radar and optical sensors able to detect objects of 1 m minimum size in the geostationary ring indicates a total population of several hundred which includes active and defunct satellites and spent upper stages. In addition, a population of untrackable objects is conjectured, whose size can only be estimated, possibly several thousand of smaller objects.

The purpose of this paper is to review the long-term evolution of orbits in the geostationary ring and at higher altitude, the collision probabilities and disposition options.

The major perturbations are considered including attitude-orbit cross-coupling effects which could cause secular orbit perturbations.

Collision probabilities for current and projected populations are reviewed considering different approaches, such as a deterministic treatment of the uncontrolled population and a stochastic modeling for the controlled satellites. Also, colocation, that is sharing of the same longitude slot by several operational satellites, is a potential source for collision, if no preventive measures are taken.

As regards spacecraft disposition options, the conclusion is that reorbiting is currently the only practical measure to safeguard the geostationary orbit. In this recommended procedure the defunct satellites are inserted into a so-called graveyard orbit, located suffieciently high above the geostationary orbit.     

Space debris mitigation strategies and practices in geosynchronous transfer orbits

Missions to geosynchronous orbits remain one of the most important elements of space launch traffic, accounting for 40% of all missions to Earth orbit and beyond during the four-year period 2000–2003. The vast majority of these missions leave one or more objects in geosynchronous transfer orbits (GTOs), contributing on a short-term or long-term basis to the space debris population. National and international space debris mitigation guidelines seek to curtail the accumulation of debris in orbits which penetrate the regions of low Earth orbit and of geosynchronous orbit. The orbital lifetime of objects in GTO can be greatly influenced by the initial values of perigee, inclination, and right ascension of the orbital plane, leading to orbital lifetimes of from less than one month to more than 100 years. An examination of the characteristic GTOs employed by launch vehicles from around the world has been conducted. The consequences of using perigees above 300 km and super-synchronous apogees, typically above 40,000 km, have been identified. In addition, the differences in orbital behavior of launch vehicle stages and mission-related debris in GTOs have been investigated. Greater coordination and cooperation between space launch service providers and spacecraft designers and owners could significantly improve overall compliance with guidelines to mitigate the accumulation of debris in Earth orbit.     

Active debris removal of multiple priority targets

Today’s space debris environment shows major concentrations of objects within distinct orbital regions for nearly all size regimes. The most critical region is found at orbital altitudes near 800 km with high declinations. Within this region many satellites are operated in so called sun-synchronous orbits (SSO). Among those, there are Earth observation, communication and weather satellites. Due to the orbital geometry in SSO, head-on encounters with relative velocities of about 15 km/s are most probable and would thus result in highly energetic collisions, which are often referred to as catastrophic collisions, leading to the complete fragmentation of the participating objects. So called feedback collisions can then be triggered by the newly generated fragments, thus leading to a further population increase in the affected orbital region. This effect is known as the Kessler syndrome.     

Dynamical evolution of high area-to-mass ratio objects in Molniya orbits

For special demands, some notable orbit types have been developed by human, including the Molniya orbits, which have a relatively high eccentricity up to about 0.7, and a period of 12 h. Considering that space debris with high area-to-mass ratio (A/M) has been discovered, such objects may also exist in Molniya orbits due to spacecraft and upper stages fragmentation events. However, there are not sufficient studies of the complex dynamical phenomena of such orbits. These studies can enrich the knowledge about the long-term evolution of these orbits, be helpful to propose uncatalogued objects observation and identification, and also set the protected region as well as active debris removal. In this paper, the characteristics of 2:1 resonance of Molniya satellite orbits are studied. A large set of numerical simulations, including all the relevant perturbations, is carried out to further investigate the main characteristics, and special attention is payed to the dynamical evolution of objects with high A/M, particularly affected by the direct solar radiation pressure. The long-term dynamical evolution of orbital elements, as well as the dependency of lifetime on the A/M value, is discussed.     

Influence of Earth’s shadowing effects on space debris stability

Solar radiation pressure affects the evolution of high area-to-mass geostationary space debris. In this work, we extend the stability study of Valk et al. (2009) by considering the influence of Earth’s shadows on the short- and long-term time evolutions of space debris. To assess the orbits stability, we use the Mean Exponential Growth factor of Nearby Orbits (MEGNO), which is an efficient numerical tool to distinguish between regular and chaotic behaviors. To reliably compute long-term space debris motion, we resort to the Global Symplectic Integrator (GSI) of Libert et al. (2011) which consists in the symplectic integration of both Hamiltonian equations of motion and variational equations. We show how to efficiently compute the MEGNO indicator in a complete symplectic framework, and we also discuss the choice of a symplectic integrator, since propagators adapted to the structure of the Hamiltonian equations of motion are not necessarily suited for the associated variational equations. The performances of our method are illustrated and validated through the study of the Arnold diffusion problem. We then analyze the effects of Earth’s shadows, using the adapted conical and cylindrical Earth’s shadowing models introduced by Hubaux et al. (2012) as the smooth shadow function deriving from these models can be easily included into the variational equations. Our stability study shows that Earth’s shadows greatly affect the global behaviour of space debris orbits by increasing the size of chaotic regions around the geostationary altitude. We also emphasize the differences in the results given by conical or cylindrical Earth’s shadowing models. Finally, such results are compared with a non-symplectic integration scheme.     

Bayesian linking of geosynchronous orbital debris tracks as seen by the Large Synoptic Survey Telescope

We describe a Bayesian sampling model for linking and constraining orbit models from angular observations of “streaks” in optical telescope images. Our algorithm is particularly suited to situations where the observation times are small fractions of the orbital periods of the observed objects or when there is significant confusion of objects in the observation field. We use Markov Chain Monte Carlo to sample from the joint posterior distribution of the parameters of multiple orbit models (up to the number of observed tracks) and parameters describing which tracks are linked with which orbit models. Using this algorithm, we forecast the constraints on geosynchronous (GEO) debris orbits achievable with the planned Large Synoptic Survey Telescope (LSST). Because of the short 15 s exposure times, preliminary orbit determinations of GEO objects from LSST will have large and degenerate errors on the orbital elements. Combined with the expected crowded fields of GEO debris it will be challenging to reliably link orbital tracks in LSST observations given the currently planned observing cadence.     

Space debris observations with the Slovak AGO70 telescope: Astrometry and light curves

The Faculty of Mathematics, Physics and Informatics of Comenius University in Bratislava, Slovakia (FMPI) operates its own 0.7-m Newtonian telescope (AGO70) dedicated to the space surveillance tracking and research, with an emphasis on space debris. The observation planning focuses on objects on geosynchronous (GEO), eccentric (GTO and Molniya) and global navigation satellite system (GNSS) orbits. To verify the system’s capabilities, we conducted an observation campaign in 2017, 2018 and 2019 focused on astrometric and photometric measurements. In last two years we have built up a light curve catalogue of space debris which is now freely available for the scientific community. We report periodic signals extracted from more than 285 light curves of 226 individual objects. We constructed phase diagrams for 153 light curves for which we obtained apparent amplitudes.     

A robotic telescope network for space debris identification and tracking

The two TAROT (Télescopes à Action Rapide pour les Objets Transitoires; Rapid Action Telescopes for Transient Objects) installations are fully robotic optical observatories with optimized observation scheduling, data processing and archiving. Zadko is a 1 m telescope in Western Australia. The fully robotisation of the Zadko telescope has just been completed; it is now included in the TAROT network. In this paper we provide an overview of this international network of robotic optical telescopes. We discuss the advantages of using the network to participate in a satellite and space debris tracking program. This network will access almost all geostationary belt objects, and provide the first real-time satellite positioning capability. The inclusion of the 1 m Zadko telescope into the network significantly extends the efficiency and sensitivity of the existing two telescope configuration.     

On prediction of re-entry time of an upper stage from GTO

The evolution of objects in geostationary transfer orbit (GTO) is determined by a complex interplay of atmospheric drag and luni-solar gravity. These orbits are highly eccentric (eccentricity >0.7) and have large variations in velocity and perturbations during a revolution. The periodic changes in the perigee altitudes of these orbits are mainly due to the gravitational perturbations of the Sun and the Moon. The re-entry time of the objects in such orbits is sensitive to the initial conditions. The aim of this paper is to study the re-entry time of the cryogenic stage of the Indian geo-synchronous launch vehicle, GSLV-F04/CS, which has been decaying since 2 September 2007 from initial orbit with eccentricity equal to 0.706. Two parameters, initial eccentricity and ballistic coefficient, are chosen for optimal estimation. It is known that the errors are more in eccentricity for the observations based on two line elements (TLEs). These two parameters are computed with response surface method using a genetic algorithm for the selected eight different zones, based on rough linear variation of the mean apogee altitude during 200 days orbit evolution. The study shows that the GSLV-F04/CS will re-enter between 5 December 2010 and 7 January 2011. The methodology is also applied to study the re-entry of six decayed objects (cryogenic stages of GSLV and Molniya satellites). Good agreement is noticed between the actual and the predicted re-entry times. The absolute percentage error in re-entry prediction time for all the six objects is found to be less than 7%. The present methodology is being adopted at Vikram Sarabhai Space Centre (VSSC) to predict the re-entry time of GSLV-F04/CS.     

Symplectic integration of space debris motion considering several Earth’s shadowing models

In this work, we present a symplectic integration scheme to numerically compute space debris motion. Such an integrator is particularly suitable to obtain reliable trajectories of objects lying on high orbits, especially geostationary ones. Indeed, it has already been demonstrated that such objects could stay there for hundreds of years. Our model takes into account the Earth’s gravitational potential, luni-solar and planetary gravitational perturbations and direct solar radiation pressure. Based on the analysis of the energy conservation and on a comparison with a high order non-symplectic integrator, we show that our algorithm allows us to use large time steps and keep accurate results. We also propose an innovative method to model Earth’s shadow crossings by means of a smooth shadow function. In the particular framework of symplectic integration, such a function needs to be included analytically in the equations of motion in order to prevent numerical drifts of the energy. For the sake of completeness, both cylindrical shadows and penumbra transitions models are considered. We show that both models are not equivalent and that big discrepancies actually appear between associated orbits, especially for high area-to-mass ratios.     

Charge exchange interactions on near-Earth proton radiation for orbit perturbation of high area-to-mass ratio objects

This novel concept expels neutral gas in the presence of geomagnetically-trapped protons in near-Earth orbit. The expelled neutral gas acts to induce charge exchange collisions with the geomagnetically-trapped protons and induce drag on objects which pass through it. The charge exchange collisions between the neutral gas and the geomagnetically-trapped protons create neutrals with similar kinetic energy that are not confined by the geomagnetic field. The charge exchange neutrals are able to collide with orbital objects and perturb their orbits. The delta-v applied by the charge exchange neutral flux is greatest on high area-to-mass objects. Numerical simulation shows charge exchange neutral impacts produce a delta-v on objects on the order of 3.8 x 10⁻¹¹ m/s at a distance of 1 km from the center of the expelled gas in a 1,000 km orbit. The impulse imparted by charge exchange neutral impacts is at least six orders of magnitude smaller than that provided by the induced drag caused by gas expulsion. The localized drag increase can force a majority of small objects into the orbit of the expelled gas cloud, even if that orbit is retrograde to the initial orbit of the objects. This new technique can be applied to the remediation of space debris.     

Orbital debris–debris collision avoidance

We focus on preventing collisions between debris and debris, for which there is no current, effective mitigation strategy. We investigate the feasibility of using a medium-powered (5 kW) ground-based laser combined with a ground-based telescope to prevent collisions between debris objects in low-Earth orbit (LEO). The scheme utilizes photon pressure alone as a means to perturb the orbit of a debris object. Applied over multiple engagements, this alters the debris orbit sufficiently to reduce the risk of an upcoming conjunction. We employ standard assumptions for atmospheric conditions and the resulting beam propagation. Using case studies designed to represent the properties (e.g. area and mass) of the current debris population, we show that one could significantly reduce the risk of nearly half of all catastrophic collisions involving debris using only one such laser/telescope facility. We speculate on whether this could mitigate the debris fragmentation rate such that it falls below the natural debris re-entry rate due to atmospheric drag, and thus whether continuous long-term operation could entirely mitigate the Kessler syndrome in LEO, without need for relatively expensive active debris removal.     

Analysis of laser ranges and angular measurements data fusion for space debris orbit determination

In the framework of space debris, the orbit determination process is a fundamental step, both, for researchers and for satellite operators. The accurate knowledge of the orbit of space debris objects is needed to allow space debris characterization studies and to avoid unnecessary collision avoidance maneuvers.The accuracy of the results of an orbit determination process depends on several factors as the number, the accuracy, the kind of processed measurements, their distribution along the orbit, and the object-observer relative geometry. When the observation coverage of the target orbit is not homogeneous, the accuracy of the orbit determination can be improved processing different kind of observables. Recent studies showed that the satellite laser ranging technique can be successfully applied to space debris.In this paper, we will investigate the benefits of using laser ranges and angular measurements for the orbit determination process. We will analyze the influence of the number of used observations, of the covered arc of orbit, of each observable, and of the observation geometry on the estimated parameters. Finally, using data acquired on short observation arcs, we analyze the achievable accuracies for the orbital regimes with the highest space debris density, and to the consequences of the data fusion on catalog maintenance operations. The results shown are obtained using only real data (both angular and laser measurements) provided by sensors of the Swiss Optical Ground Station and Geodynamics Observatory Zimmerwald owned by the Astronomical Institute of the University of Bern (AIUB) and for some studies also using ranges provided from other stations of the International Laser Ranging Service (ILRS).