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.     

Accurate orbit predictions for debris orbit manoeuvre using ground-based lasers

Orbit manoeuvre of low Earth orbiting (LEO) debris using ground-based lasers has been proposed as a cost-effective means to avoid debris collisions. This requires the orbit of the debris object to be determined and predicted accurately so that the laser beam can be locked on the debris without the loss of valuable laser operation time. This paper presents the method and results of a short-term accurate LEO (<900 km in altitude) debris orbit prediction study using sparse laser ranging data collected by the EOS Space Debris Tracking System (SDTS). A main development is the estimation of the ballistic coefficients of the LEO objects from their archived long-term two line elements (TLE). When an object is laser tracked for two passes over about 24 h, orbit prediction (OP) accuracy of 10–20 arc seconds for the next 24–48 h can be achieved – the accuracy required for laser debris manoeuvre. The improvements in debris OP accuracy are significant in other applications such as debris conjunction analyses and the realisation of daytime debris laser tracking.     

MCMC-Particle-based group tracking of space objects within Bayesian framework

With the intense increase in space objects, especially space debris, it is necessary to efficiently track and catalog the extensive dense clusters of space objects. As the main instrument for low earth orbit (LEO) space surveillance, ground-based radar system is usually limited by its resolution while tracking small space debris with high density. Thus, the obtained measurement information could have been seriously missed, which makes the traditional tracking method inefficient. To address this issue, we conceived the concept of group tracking. For group tracking, the overall tendency of the group objects is expected to be revealed, and the trajectories of individual objects are simultaneously reconstructed explicitly. According to model the interaction between the group center and individual trajectories using the Markov random field (MRF) within Bayesian framework, the objects’ number and individual trajectory can be estimated more accurately in the condition of high miss alarm probability. The Markov chain Monte Carlo (MCMC)-Particle algorithm was utilized for solving the Bayesian integral problem. Furthermore, we introduced the mechanism for describing the behaviors of groups merging and splitting, which can expand the single group tracking algorithm to track variable multiple groups. Finally, simulation of the group tracking of space objects was carried out to validate the efficiency of the proposed method.     

Uncertainty-aware Cube algorithm for medium-term collision risk assessment

The number of Earth orbiting objects is constantly growing, and some orbital regions are becoming risky environments for space assets of interest, which are increasingly threatened by accidental collisions with other objects, especially in Low-Earth Orbit (LEO). Collision risk assessment is performed by various methods, both covariance and non-covariance based. The Cube algorithm is a non-covariance-based method used to estimate the collision rates between space objects, whose concept consists in dividing the space in cubes of fixed dimension and, at each time instant, checking if two or more objects share the same cube. Up to now its application has been limited to the long-term scenarios of orbital debris evolutionary models, where considering the uncertainties is not necessary and impractical. Within operative contexts, instead, medium-term collision risk analysis may be an important task, in which the propagation-related uncertainties play a prominent role, but the timescale poses challenges for the application of standard covariance-based conjunction analysis techniques. In this framework, this paper presents an approach for the evaluation of the medium-term collision frequency for objects in LEO, called Uncertainty-aware Cube method. It is a modified version of the Cube, able to take the possible errors in the space objects’ position into account for the detection of the conjunctions. As an object’s orbit is propagated, the along-track position error grows more and more, and each object could potentially be in a different position with respect to the one determined by numerical propagation and, thus, in a different cube. Considering the uncertainties, at each time instant the algorithm associates more than one cube to each object and checks if they share at least one cube. If so, a conjunction is detected and a degree of confidence is evaluated. The performance of the method is assessed in different LEO scenarios and compared to the original Cube method.     

Analysis of the orbit lifetime of CubeSats in low Earth orbits including periodic variation in drag due to attitude motion

It is estimated that more than 22,300 human-made objects are in orbit around the Earth, with a total mass above 8,400,000 kg. Around 89% of these objects are non-operational and without control, which makes them to be considered orbital debris. These numbers consider only objects with dimensions larger than 10 cm. Besides those numbers, there are also about 2000 operational satellites in orbit nowadays. The space debris represents a hazard to operational satellites and to the space operations. A major concern is that this number is growing, due to new launches and particles generated by collisions. Another important point is that the development of CubeSats has increased exponentially in the last years, increasing the number of objects in space, mainly in the Low Earth Orbits (LEO). Due to the short operational time, CubeSats boost the debris population. One of the requirements for space debris mitigation in LEO is the limitation of the orbital lifetime of the satellites, which needs to be lower than 25 years. However, there are space debris with longer estimated decay time. In LEÓs, the influence of the atmospheric drag is the main orbital perturbation, and is used in maneuvers to increment the losses in the satellite orbital energy, to locate satellites in constellations and to accelerate the decay.The goal of the present research is to study the influence of aerodynamic rotational maneuver in the CubeSat?s orbital lifetime. The rotational axis is orthogonal to the orbital plane of the CubeSat, which generates variations in the ballistic coefficient along the trajectory. The maneuver is proposed to accelerate the decay and to mitigate orbital debris generated by non-operational CubeSats. The panel method is selected to determine the drag coefficient as a function of the flow incident angle and the spinning rate. The pressure distribution is integrated from the satellite faces at hypersonic rarefied flow to calculate the drag coefficient. The mathematical model considers the gravitational potential of the Earth and the deceleration due to drag. To analyze the effects of the rotation during the decay, multiple trajectories were propagated, comparing the results obtained assuming a constant drag coefficient with trajectories where the drag coefficient changes periodically. The initial perigees selected were lower than 400 km of altitude with eccentricities ranging from 0.00 to 0.02. Six values for the angular velocity were applied in the maneuver. The technique of rotating the spacecraft is an interesting solution to increase the orbit decay of a CubeSat without implementing additional de-orbit devices. Significant changes in the decay time are presented due to the increase of the mean drag coefficient calculated by the panel method, when the maneuver is applied, reducing the orbital lifetime, however the results are independent of the angular velocity of the satellite.     

Characterizing space debris longitude-dependent distribution based on RAAN perturbation rate

The space debris environment is one of the major threats against payloads. Space debris orbital distribution is of great importance for space debris environment modeling. Due to perturbation factors, the Right Ascension of Ascending Node (RAAN) of space objects changes consistently, causing regular rotation of the orbit plane around Earth’s axis. Based on the investigation of the RAAN perturbation rate of concerned objects, this paper proposes a RAAN discretization method in order to present the space debris longitude-dependent distribution. Combined with two line element (TLE) data provided by the US Space Surveillance Network, the estimated value from RAAN discretization method is compared with the real case. The results suggest that using only the initial orbital data at the beginning of the time interval of interest, the RAAN discretization method is able to provide reliable longitude distribution of concerned targets in the next following period. Furthermore, spacecraft cumulative flux against space debris is calculated in this paper. The results suggest that the relevance between spacecraft RAAN setup and flux output is much smaller for LEO targets than MEO targets, which corresponds with the theory analysis. Since the nonspherical perturbation is the major factor for RAAN variation, the RAAN perturbation rate has little connection with the size of orbital objects. In other words, the RAAN discretization method introduced in this paper also applies to space debris of different size range, proposing a possible suggestion for the improvement of space debris environment engineering models.     

Instability of the present LEO satellite populations

Several studies conducted during 1991–2001 demonstrated, with some assumed launch rates, the future unintended growth potential of the Earth satellite population, resulting from random, accidental collisions among resident space objects. In some low Earth orbit (LEO) altitude regimes where the number density of satellites is above a critical spatial density, the production rate of new breakup debris due to collisions would exceed the loss of objects due to orbital decay.     

Removing orbital debris with lasers

Orbital debris in low Earth orbit (LEO) are now sufficiently dense that the use of LEO space is threatened by runaway collision cascading. A problem predicted more than thirty years ago, the threat from debris larger than about 1 cm demands serious attention. A promising proposed solution uses a high power pulsed laser system on the Earth to make plasma jets on the objects, slowing them slightly, and causing them to re-enter and burn up in the atmosphere. In this paper, we reassess this approach in light of recent advances in low-cost, light-weight modular design for large mirrors, calculations of laser-induced orbit changes and in design of repetitive, multi-kilojoules lasers, that build on inertial fusion research. These advances now suggest that laser orbital debris removal (LODR) is the most cost-effective way to mitigate the debris problem. No other solutions have been proposed that address the whole problem of large and small debris. A LODR system will have multiple uses beyond debris removal. International cooperation will be essential for building and operating such a system.     

Track detection of high-velocity resident space objects in Low Earth Orbit

The world’s economy has become heavily dependent on the services provided by satellites. With the exponential increase in satellite launches, the population of defunct or inactive hardware in space has grown substantially. This is especially true in sensitive orbits such as the Low Earth Orbit (LEO) and Geostationary Earth Orbit (GEO) regimes. These objects, collectively known as orbital debris, can reach speeds of up to 28 000km h^?1 in LEO. At these orbital speeds, even the smallest of objects can pose a considerable threat to operational satellites or astronauts. This makes the monitoring, and detection, of these objects of the utmost importance. This work describes the latest detection strategy used in one of Europe’s largest Space Situational Awareness (SSA) installation; the BIstatic RAdar for LEo Survey (BIRALES) space debris radar. We present a novel bottom-up approach that makes use of single-linkage clustering to identify faint radar streaks in spectrogram data. Tests on synthetic data have shown that the detection strategy presented in this study obtains a higher detection rate when it is compared against existing methods. Unlike other approaches, this detection strategy, using the Multi-beam streak detection strategy (MSDS) algorithm, was still able to recall 90% of the track information at an Signal-to-Noise Ratio (SNR) of 2dB.     

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.     

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.     

雷达波束驻留模式下空间碎片数量的置信区间估计

碎片数量估计是空间碎片环境统计特征描述的重要内容之一,对于空间碎片环境模型验证、航天器碰撞风险分析以及碎片数量增长趋势预测有重要意义.针对波束指向正东、正南任意仰角的雷达波束驻留(Beam-park)模式(天顶指向是波束指向仰角为90°时的特例),给出了一种估计碎片数量置信区间的方法.对于给定轨道高度范围内一个具有穿越雷达波束可能性(即雷达散射截面足够大,且轨道倾角相对测站纬度足够大)的碎片,将其是否真正穿越波束这一随机事件用(0-1)分布来建模,根据所采集的轨道高度和倾角数据,计算出该轨道高度范围内碎片穿越波束的平均概率,进而采用中心极限定理来估计碎片数量的置信区间.仿真结果表明了方法的有效性.      

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.     

Feasibility analysis of LEO and GEO large space debris de/re-orbiting taking into account launch mass of spacecraft-collector and its configuration layout

Analysis of the efficiency of two basic strategies for de/re-orbiting large space debris objects to disposal orbits (DO) is given. Large objects in LEO are classified into groups with similar orbital inclinations and comprise primarily last stages of launch vehicles, in GEO vicinity the paper studies upper stages. Under the first de/re-orbiting variant, it is assumed a spacecraft-collector is equipped with several thruster de/re-orbiting kits (TDKs); one of them can be fixed on an object and is capable of de/re-orbiting an object to a DO independently of the collector. In the second variant, a collector operates as a space tug: transfers objects to a DO and then returns to the next objects in line. The authors study possible configuration layouts of collectors in LEO and near GEO. The available analogous projects are analyzed. The efficiency of both de/re-orbiting variants can be properly compared using the estimations of collector's dry mass and having at one's disposal the parameters of the maneuvers required for transfers between all objects in the group. As reasonable criteria of effectiveness, one can consider (separately or jointly) the launch mass of an equipped collector, its ΔV budget, and the required number of such active spacecraft. Two de/re-orbiting variants are compared in terms of these criteria via mass-energy diagrams constructed for each group of objects in both altitude regions. Analysis of these diagrams shows that low Earth orbits can be more efficiently cleaned under the first de-orbiting variant by using a two-stage space system consisting of an active spacecraft carrying TDKs. For GEO, it is expedient to choose the second re-orbiting variant using a single-stage spacecraft. Our analysis shows that LEO cleaning is an order of magnitude more expensive than that for GEO, hence the problem of LEO population should be given increased attention.     

对一种利用人造粉尘清除空间碎片新方法的理论分析

近地轨道的空间碎片污染日益严重,目前碎片数量已达到历史最高值并可能引发一系列连锁反应.进行主动碎片清除十分必要.利用粉尘主动清除近地轨道空间碎片是一种主动碎片清除的新方法.本文基于此方法的基本原理进行分析研究,建立了单颗碎片与人造粉尘作用的基本假设和机理模型,并对其作用进行定量计算分析;结合碎片的空间密度分布,对该方法的作用效果进行了定量估算,得出一些基本分析结论,有助于对新方法的客观认识.      

A novel data association scheme for LEO space debris surveillance based on a double fence radar system

The increasing amount of space debris threatens to seriously deteriorate and damage space-based instruments in Low Earth Orbit (LEO) environments. Therefore, LEO space debris surveillance systems must be developed to provide situational awareness in space and issue warnings of collisions with LEO space debris. In this paper, a double fence radar system is proposed as an emerging paradigm for LEO space debris surveillance. This system exhibits several unique and promising characteristics compared with existing surveillance systems. In this paper, we also investigate the data association scheme for LEO space debris surveillance based on a double fence radar system. We also perform a theoretical analysis of the performance of our proposed scheme. The superiority and the effectiveness of our novel data association scheme is demonstrated by experimental results. The data used in our experiments is the LEO space debris catalog produced by the North American Air Defense Command (NORAD) up to 2009, especially for scenarios with high densities of LEO space debris, which were primarily produced by the collisions between Iridium 33 and Cosmos 2251. We hope that our work will stimulate and benefit future work on LEO space debris surveillance approaches and enable construction of the double fence radar system.     

Development of Operational Space Weather Prediction Models

In this report, we summarize the needs of space weather models, and recommend that developing operational prediction models, rather than transitioning from research to operation, is a more feasible and critical way for space weather services in the near future. Operational models for solar wind speed, geomagnetic indices, magnetopause, plasma sheet energetic electrons, inner boundary of ion plasma sheet, energetic electrons in outer radiation belt, and thermospheric density at low Earth orbit, have been developed and will be introduced briefly here. Their applications made a big progress in space weather services during the past two years in China.      

High accuracy satellite drag model (HASDM)

The dominant error source in force models used to predict low-perigee satellite trajectories is atmospheric drag. Errors in operational thermospheric density models cause significant errors in predicted satellite positions, since these models do not account for dynamic changes in atmospheric drag for orbit predictions. The Air Force Space Battlelab’s High Accuracy Satellite Drag Model (HASDM) estimates and predicts (out three days) a dynamically varying global density field. HASDM includes the Dynamic Calibration Atmosphere (DCA) algorithm that solves for the phases and amplitudes of the diurnal and semidiurnal variations of thermospheric density near real-time from the observed drag effects on a set of Low Earth Orbit (LEO) calibration satellites. The density correction is expressed as a function of latitude, local solar time and altitude. In HASDM, a time series prediction filter relates the extreme ultraviolet (EUV) energy index E_10.7 and the geomagnetic storm index a_p, to the DCA density correction parameters. The E_10.7 index is generated by the SOLAR2000 model, the first full spectrum model of solar irradiance. The estimated and predicted density fields will be used operationally to significantly improve the accuracy of predicted trajectories for all low-perigee satellites.     

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.