An active debris removal parametric study for LEO environment remediation

Recent analyses on the instability of the orbital debris population in the low Earth orbit (LEO) region and the collision between Iridium 33 and Cosmos 2251 have reignited interest in using active debris removal (ADR) to remediate the environment. There are, however, monumental technical, resource, operational, legal, and political challenges in making economically viable ADR a reality. Before a consensus on the need for ADR can be reached, a careful analysis of its effectiveness must be conducted. The goal is to demonstrate the need and feasibility of using ADR to better preserve the future environment and to explore different operational options to maximize the benefit-to-cost ratio. This paper describes a new sensitivity study on using ADR to stabilize the future LEO debris environment. The NASA long-term orbital debris evolutionary model, LEGEND, is used to quantify the effects of several key parameters, including target selection criteria/constraints and the starting epoch of ADR implementation. Additional analyses on potential ADR targets among the existing satellites and the benefits of collision avoidance maneuvers are also included.     

Ariane 5 GTO debris mitigation using natural perturbations

GTO objects can potentially collide with operative satellites in LEO and GEO protected regions. Internationally accepted debris mitigation guidelines require that these objects exit these protected regions within 25?years, e.g. by re-entering and burning up in Earth’s atmosphere. In this paper, an inventory of the GTO debris generated from Ariane 5 launches in the period 2012–2017 is provided, and it is expected that none of these objects will re-enter within 25?years. For future launches, natural perturbations can be exploited to increase compliance with mitigation guidelines without the use of extra propellant or complex de-orbiting systems, which is attractive from an economic point of view. The lifetime of GTO objects is very sensitive to initial conditions and some environmental and body-related parameters, mainly due to the effect of solar gravity on the perigee altitude. As a consequence, the lifetime of a specific GTO object cannot be predicted accurately, but its probability of re-entering in less than 25?years can be estimated with proper accuracy by following a statistical approach. By propagating the orbits of over 800,000 simulated Ariane 5 GTO objects, it was found that the launch time leading to the highest probability of compliance with debris mitigation guidelines for GEO launches from Kourou corresponds to about 2 PM local time, regardless of the date of launch, which leads to compliance rates ranging from 60 to 100%. Current practice is to launch at around 5–9?PM, so a change in procedures would be required in order to reach a higher degree of compliance with debris mitigation guidelines, which was predicted to be on average below 20% for the objects generated in the period 2012–2017.     

Stabilization of an injected conducting layer for artificially enhancing drag on orbital debris

The evolution of a magnetized conducting medium suspended in magnetic and gravitational fields is examined. In this paper some effects of the influence of velocity fields on the linear stability properties of such layers are investigated. A fully compressible, three-dimensional analysis of the layer is described. The relevant equations are derived and then solved by the MagnetoHydroDynamic SPEctral Compressible Linear Stability (MHDSPECLS) algorithm, a Chebyshev collocation code. The code allows for the computation of magnetic and thermal effects. A complete stabilization of the system is found above a critical velocity of approximately 2500 m/s.     

Identification of a debris cloud from the nuclear powered SNAPSHOT satellite with haystack radar measurements

Data from the Massachusetts Institute of Technology Lincoln Laboratory Long Range Imaging Radar (known as the Haystack radar) have been used in the past to examine families of objects from individual satellite breakups or families of orbiting objects that can be isolated in altitude and inclination. This is possible because, for some time after a breakup, the debris cloud of particles can remain grouped together in similar orbit planes. This cloud will be visible to the radar, in fixed staring mode, for a short time twice each day, as the orbit plane moves through the field of view. There should be a unique three-dimensional pattern in observation time, range, and range rate which can identify the cloud. Eventually, through slightly differing precession rates of the right ascension of ascending node of the debris cloud, the observation time becomes distributed so that event identification becomes much more difficult.     

Assessment of space debris collisions against spacecraft with deorbit devices

Deorbit methods have been employed to remove space debris from orbit. One of these methods is to utilize atmospheric drag. In this method, a membrane loaded into the spacecraft is expanded to increase atmospheric drag. Although this method works without requiring fuel, it has the disadvantage of a high risk of collision with other debris owing to its larger area. Area-time product and energy-to-mass ratio have been used as indices to evaluate the risk of collisions between spacecraft and debris. However, the evaluation criteria were uncertain because these two indices are independent. In this paper, we propose a new evaluation index, single-sheet collision factor (SSCF), that comprehensively evaluates the collision risk based on experiments simulating debris collisions. As a result of the hypervelocity collision experiment, we found that the penetration-area mass of the spacecraft affects the severity of debris collisions. In this paper, the product of the exterior-wall thickness, the exterior-wall density, and the space debris cross-sectional area defines the penetration-area mass of the spacecraft. Furthermore, we compare and evaluate various deorbit methods using SSCF. The comparison showed that the penetration-area mass of the SSCF could be quantitatively determined for the debris-collision severity due to difference in structural materials of spacecraft. SSCF will be used to create rules for space-environment conservation with the expansion of the space-development market.     

A new algorithm for optical observations of space debris with the TAROT telescopes

Since 2004, we observe satellites in the geostationary orbit with a network of robotic ground based fully automated telescopes called TAROT. One of them is located in France and the second at ESO, La Silla, Chile. The system processes the data in real time. Its wide field of view is useful for the discovery, the systematic survey and for the tracking of both catalogued and un-catalogued objects. We present a new source extraction algorithm based on morphological mathematic, which has been tested and is currently under implementation in the standard pipeline. Using this method, the observation strategy will correlate the measurements of the same object on successive images and give better detection rate and false alarm rate than the previous one. The overall efficiency and quality of the survey of the geostationary orbit has drastically improved and we can now detect satellites and debris in different orbits like Geostationary Transfer Orbit (GTO). Results obtained in real conditions with TAROT are presented.     

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.     

Faster algorithm of debris cloud orbital character from spacecraft collision breakup

Space debris is polluting the space environment. Collision fragment is its important source. NASA standard breakup model, including size distributions, area-to-mass distributions, and delta velocity distributions, is a statistic experimental model used widely. The general algorithm based on the model is introduced. But this algorithm is difficult when debris quantity is more than hundreds or thousands. So a new faster algorithm for calculating debris cloud orbital lifetime and character from spacecraft collision breakup is presented first. For validating the faster algorithm, USA 193 satellite breakup event is simulated and compared with general algorithm. Contrast result indicates that calculation speed and efficiency of faster algorithm is very good. When debris size is in 0.01–0.05 m, the faster algorithm is almost a hundred times faster than general algorithm. And at the same time, its calculation precision is held well. The difference between corresponding orbital debris ratios from two algorithms is less than 1% generally.     

Uncertainty propagation for statistical impact prediction of space debris

Predictions of the impact time and location of space debris in a decaying trajectory are highly influenced by uncertainties. The traditional Monte Carlo (MC) method can be used to perform accurate statistical impact predictions, but requires a large computational effort. A method is investigated that directly propagates a Probability Density Function (PDF) in time, which has the potential to obtain more accurate results with less computational effort. The decaying trajectory of Delta-K rocket stages was used to test the methods using a six degrees-of-freedom state model. The PDF of the state of the body was propagated in time to obtain impact-time distributions. This Direct PDF Propagation (DPP) method results in a multi-dimensional scattered dataset of the PDF of the state, which is highly challenging to process. No accurate results could be obtained, because of the structure of the DPP data and the high dimensionality. Therefore, the DPP method is less suitable for practical uncontrolled entry problems and the traditional MC method remains superior. Additionally, the MC method was used with two improved uncertainty models to obtain impact-time distributions, which were validated using observations of true impacts. For one of the two uncertainty models, statistically more valid impact-time distributions were obtained than in previous research.     

Properties of the high area-to-mass ratio space debris population at high altitudes

In the framework of its space debris research activities ESA established an optical survey program to study the space debris environment at high altitudes, in particular in the geostationary ring and in the geostationary transfer orbit region. The Astronomical Institute of the University of Bern (AIUB) performs these surveys on behalf of ESA using ESA’s 1-m telescope in Tenerife. Regular observations were started in 1999 and are continued during about 120–140 nights per year. Results from these surveys revealed a substantial amount of space debris at high altitudes in the size range from 0.1 to 1 m. Several space debris populations with different dynamical properties were identified in the geostationary ring. During the searches for debris in the geostationary transfer orbit region a new population of objects in unexpected orbits, where no potential progenitors exist, was found. The orbital periods of these objects are clustered around one revolution per day; the eccentricities, however, are scattered between 0 and 0.6. By following-up some of these objects using the ESA telescope and AIUB’s 1-m telescope in Zimmerwald, Switzerland, it was possible to study the properties of this new population. One spectacular finding from monitoring the orbits over time spans of days to months is the fact that these objects must have extreme area-to-mass ratios, which are by several orders of magnitudes higher than for ‘normal-type’ debris. This in turn supports the hypothesis that the new population actually is debris generated in or near the geostationary ring and which is in orbits with periodically varying eccentricity and inclination due to perturbations by solar radiation pressure. In order to further study the nature of these debris, multi-color and temporal photometry (light curves) were acquired with the Zimmerwald telescope. The light curves show strong variations over short time intervals, including signals typical for specular reflections. Some objects exhibit distinct periodic variations with periods ranging from 10 to several 100 s. All this is indicative for objects with complicated shapes and some highly reflective surfaces.     

A novel signal processing approach for LEO space debris based on a fence-type space surveillance radar system

The increase in space debris can seriously threaten regular activities in the Low Earth Orbit (LEO) environment. Therefore, it is necessary to develop robust, efficient and reliable techniques to understand the potential motions of the LEO debris. In this paper, we propose a novel signal processing approach to detect and estimate the motions of LEO space debris that is based on a fence-type space surveillance radar system. Because of the sparse distribution of the orbiting debris through the fence in our observations, we formulate the signal detection and the motion parameter estimation as a sparse signal reconstruction problem with respect to an over-complete dictionary. Moreover, we propose a new scheme to reduce the size of the original over-complete dictionary without the loss of the important information. This new scheme is based on a careful analysis of the relations between the acceleration and the directions of arrival for the corresponding LEO space debris. Our simulation results show that the proposed approach can achieve extremely good performance in terms of the accuracy for detection and estimation. Furthermore, our simulation results demonstrate the robustness of the approach in scenarios with a low Signal-to-Noise Ratio (SNR) and the super-resolution properties. We hope our signal processing approach can stimulate further work on monitoring LEO space debris.     

Using parallel computing for the display and simulation of the space debris environment

Parallelism is becoming the leading paradigm in today’s computer architectures. In order to take full advantage of this development, new algorithms have to be specifically designed for parallel execution while many old ones have to be upgraded accordingly. One field in which parallel computing has been firmly established for many years is computer graphics. Calculating and displaying three-dimensional computer generated imagery in real time requires complex numerical operations to be performed at high speed on a large number of objects. Since most of these objects can be processed independently, parallel computing is applicable in this field. Modern graphics processing units (GPUs) have become capable of performing millions of matrix and vector operations per second on multiple objects simultaneously.     

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.     

An adaptive threshold method for improving astrometry of space debris CCD images

Optical survey is a main technique for observing space debris, and precisely measuring the positions of space debris is of great importance. Due to several factors, e.g. the angle object normal to the observer, the shape as well as the attitude of the object, the variations of observed characteristics for low earth orbital space debris are distinct. When we look at optical CCD images of observed objects, the size and brightness are varying, hence it’s difficult to decide the threshold during centroid measurement and precise astrometry. Traditionally the threshold is given empirically and constantly in data reduction, and obviously it’s not suitable for data reduction of space debris. Here we offer a solution to provide the threshold. Our method assumes that the PSF (point spread function) is Gaussian and estimates the signal flux by a directly two-dimensional Gaussian fit, then a cubic spline interpolation is performed to divide each initial pixel into several sub-pixels, at last the threshold is determined by the estimation of signal flux and the sub-pixels above threshold are separated to estimate the centroid. A trail observation of the fast spinning satellite Ajisai is made and the CCD frames are obtained to test our algorithm. The calibration precision of various threshold is obtained through the comparison between the observed equatorial position and the reference one, the latter are obtained from the precise ephemeris of the satellite. The results indicate that our method reduces the total errors of measurements, it works effectively in improving the centering precision of space debris images.     

Fuzzy multi-criteria decision-making approach to prioritization of space debris for removal

As private companies and government space agencies begin to seriously consider the task of active space debris removal, it is becoming increasingly more important to determine the highest priority objects to deorbit. This work sets forth an approach for prioritization of space debris through the utilization of Multi-Criteria Decision-Making methodologies and fuzzy logic, as well as both quantitative and qualitative criteria. The proposed debris prioritization approach considers various criteria including the orbit, size, mass, pairwise and total collision probabilities, and decay timeframe of each debris object. The means of assigning attributes to each assessment criterion is discussed in detail. To determine the weighting scheme for the criteria, a questionnaire was prepared and shared with experts in the field of space situational awareness. The work examines over two thousand critical debris objects selected from the existing debris catalog with respect to these criteria. The quantified attributes for each debris object are then aggregated through the fuzzy versions of the Analytic Hierarchy Process and the Technique for Order Preference by Similarity to Ideal Solution. The results of the analysis identify high-priority debris objects for removal from Earth-bound orbits.     

Two-dimensional phase plane structure and the stability of the orbital motion for space debris in the geosynchronous ring

Based on the orbital resonance model, we study the two-dimensional phase plane structure of the motion of space debris orbiting the geosynchronous ring under the combined effects of the tesseral harmonics J₂₂, J₃₁ and J₃₃ of the Earth’s gravitational field. We present the main characteristic parameters of the two-dimensional phase plane structure. We also analyze the stability of the two-dimensional phase plane structure with numerical method. Our main findings indicate that the combined effects of the tesseral harmonics J₂₂, J₃₁ and J₃₃ fully determine the two-dimensional phase plane structure of the space debris, and it remains robust under the effect of the Earth’s actual gravitational field, the luni-solar perturbations and the solar radiation pressure with the normal area-to-mass ratios.     

First astrometric observations of space debris with the MéO telescope

The MéO (for Métrologie Optique) telescope is the Satellite and Lunar Laser Ranging (SLR) dedicated telescope of Observatoire de la Côte d’Azur (France) located at plateau de Calern. The telescope uses an altazimuth mount. The motorization of the mount has a capability of 6 deg/s allowing the follow up of Low Earth Orbits (LEO) satellites, as well as Medium Earth Orbits (MEO) and geostationary (GEO) satellites, and the Moon. The telescope has a primary mirror of 1.54 m. It uses a Nasmyth focus equipped with an EMCCD camera. The telescope field of view, defined by the equivalent focal length and the size of the camera, is currently 3.4 arcmin × 3.4 arcmin.     

On the practical exploitation of perturbative effects in low Earth orbit for space debris mitigation

This paper presents the results of a numerical evaluation of the natural lifetime reduction in low Earth orbit, due to dynamical perturbations. The study considers two values for the area-to-mass ratio, a nominal ratio which resembles a typical value of spacecraft in orbit today, and an enhanced ratio which covers the surface augmentation. The results were obtained with two orbit propagators, one of a semi-analytical nature and the second one using non-averaged equations of motion. The simulations for both propagators were set up similarly to allow comparison. They both use the solar radiation pressure and the secular terms of the geopotential ($J_2,J_4$ and $J_6$). The atmospheric drag was turned on and off in both propagators to alternatively study the eccentricity build up and the residual lifetime. The non-averaging case also covers a validation with the full 6?×?6 geopotential. The results confirm the findings in previous publications, that is, the possibility for de-orbiting from altitudes above the residual atmosphere if a solar sail is deployed at the end-of-life, due to the combined effect of solar radiation pressure and the oblateness of the Earth. At near polar inclinations, shadowing effects can be exploited to the same end. The results obtained with the full, non-averaging propagator revealed additional de-orbiting corridors associated with solar radiation pressure which were not found by previous work on space debris mitigation. The results of both tools are compared for specific initial conditions. For nominal values of area-to-mass ratio, instead, it is confirmed that this resonance effect is negligible.The paper then puts the findings in the perspective of the current satellite catalogue. It identifies space missions which are currently close to a resonance corridor and shows the orbit evolution within the resonances with a significantly shorter residual orbital lifetime. The paper finishes with a discussion on the exploitation of these effects with regards to the long-term simulation of the space debris environment and a flux and collision probability comparison.     

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.