Shape and motion estimate of LEO debris using light curves

We have derived a tri-axial ellipsoidal model of an LEO object, a Cosmos 2082 rocket body, including its rotational axis direction, rotation period, precession, and a compositional parameter, using only light curve data from an optical telescope. The brightness of the object was monitored for two days and least-squares fitting was used to determine these values. The derived axial ratios are 100:18:18, the coordinates of the rotational axis direction on the celestial sphere are R.A. = 305.8° and Dec. = 2.6°, and its observed average rotation period is 41 s. When precession is considered, its amplitude and precession period are 30.5° and 29.4 min. These results show that optical light curve data are sufficient to determine the shapes and the motions of LEO objects.     

First laser measurements to space debris in Poland

The Borowiec Satellite Laser Ranging station (BORL 7811, Borowiec) being a part of the Space Research Centre of the Polish Academy of Sciences (SRC PAS) went through modernization in 2014–2015. One of the main tasks of the modernization was the installation of a high-energy laser module dedicated to space debris tracking. Surelite III by Continuum is a Nd:YAG pulse laser with 10?Hz repetition rate, a pulse width of 3–5?ns and a pulse energy of 450?mJ for green (532?nm). This new laser unit was integrated with the SLR system at Borowiec performing standard satellite tracking. In 2016 BORL 7811 participated actively to the observational campaigns related to the space debris targets from LEO region managed by the Space Debris Study Group (SDSG) of the International Laser Ranging Service (ILRS).Currently, Borowiec station regularly tracks 36 space debris from the LEO regime, including typical rocket bodies (Russian/Chinese) and cooperative targets like the inactive TOPEX/Poseidon, ENVISAT, OICETS and others. In this paper the first results of space debris laser measurements obtained by the Borowiec station in period August 2016 – January 2017 are presented. The results gained by the SRC PAS Borowiec station confirm the rotation of the defunct TOPEX/Poseidon satellite which spins with a period of approximately 10?s. The novelty of this work is the presentation of the sample results of the Chinese CZ-2C R/B target (NORAD catalogue number 31114) which is equipped (probably) with retroreflectors. Laser measurements to space debris is a very desirable topic for the next years, especially in the context of the Space Surveillance and Tracking (SST) activity. Some targets are very easy to track like defunct ENVISAT or TOPEX/Poseidon. On the other hand, there is a big population of different LEO targets with different orbital and physical parameters, which are challenging for laser ranging like small irregular debris and rocket boosters.     

Estimation of orbital parameters of broken-up objects from in-situ debris measurements

Even sub-millimeter-size debris could cause a fatal damage on a spacecraft. Such tiny debris cannot be followed up or tracked from the ground. Therefore, Kyushu University has initiated IDEA the project for In-situ Debris Environmental Awareness, which conducts in-situ measurements of sub-millimeter-size debris. One of the objectives is to estimate the location of on-orbit satellite fragmentations from in-situ measurements. The previous studies revealed that it is important to find out the right nodal precession rate to estimate the orbital parameters of a broken-up object properly. Therefore, this study derives a constraint equation that applies to the nodal precession rate of the broken-up object. This study also establishes an effective procedure to estimate properly the orbital parameters of a broken-up object with the constraint equation.     

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.     

Impact effects from small size meteoroids and space debris

When the impact risk from meteoroids and orbital debris is assessed the main concern is usually structural damage. With their high impact velocities of typically 10–20 km/s millimeter or centimeter sized objects can puncture pressure vessels and other walls or lead to destruction of complete subsystems or even whole spacecraft. Fortunately chances of collisions with such larger objects are small (at least at present). However, particles in the size range 1–100 μm are far more abundant than larger objects and every orbiting spacecraft will encounter them with certainty. Every solar cell (8 cm² area) of the Hubble Space Telescope encountered on average 12 impacts during its 8.25 years of space exposure. Most were from micron sized particles.     

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.     

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.     

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.     

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.     

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 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.     

Comparison of self-healing ionomer to aluminium-alloy bumpers for protecting spacecraft equipment from space debris impacts

This paper discusses the impact behavior of a self-healing ionomeric polymer and compares its protection capability against space debris impacts to that of simple aluminium-alloy bumpers. To this end, 14 impact experiments on both ionomer and Al-7075-T6 thin plates with similar surface density were made with 1.5 mm aluminium spheres at velocity between 1 and 4 km/s.     

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.     

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.     

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.     

The bistatic radar capabilities of the Medicina radiotelescopes in space debris detection and tracking

An accurate measurement of the position and trajectory of the space debris fragments is of primary importance for the characterization of the orbital debris environment. The Medicina Radioastronomical Station is a radio observation facility that is here proposed as receiving part of a ground-based space surveillance system for detecting and tracking space debris at different orbital regions (from Low Earth Orbits up to Geostationary Earth Orbits). The proposed system consists of two bistatic radars formed by the existing Medicina receiving antennas coupled with appropriate transmitters. This paper focuses on the current features and future technical development of the receiving part of the observational setup. Outlines of possible transmitting systems will also be given together with the evaluation of the observation strategies achievable with the proposed facilities.     

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.     

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.     

Semi-analytical theory of mean orbital motion for geosynchronous space debris under gravitational influence

This paper provides a hamiltonian formulation of the equations of motion of an artificial satellite or space debris orbiting the geostationary ring. This theory of order 1 has been formulated using canonical and non-singular elements for eccentricity and inclination. The analysis is based on an expansion in powers of the eccentricity and of the inclination. The theory accounts for the influence of the Earth gravity field expanded in spherical harmonics, paying a particular attention to the resonance occurring for geosynchronous objects. The luni-solar perturbations are also taken into account. We present the resonant motion and its main characteristics: equilibria, stability, fundamental frequencies and width of the resonant area by comparison with a basic analytical model. Finally, we show some results concerning the long term dynamics of a typical space debris under the influence of the gravitational field of the Earth and the luni-solar interactions.