A high order method for orbital conjunctions analysis: Sensitivity to initial uncertainties

A high order method to quickly assess the effect that uncertainties produce on orbital conjunctions through a numerical high-fidelity propagator is presented. In particular, the dependency of time and distance of closest approach to initial uncertainties on position and velocity of both objects involved in a conjunction is studied. The approach relies on a numerical integration based on differential algebraic techniques and a high-order algorithm that expands the time and distance of closest approach in Taylor series with respect to relevant uncertainties. The modeled perturbations are atmospheric drag, using NRLMSISE-00 air density model, solar radiation pressure with shadow, third body perturbation using JPL’s DE405 ephemeris, and EGM2008 gravity model. The polynomial approximation of the final position is used as an input to compute analytically the expansion of time and distance of closest approach. As a result, the analysis of a close encounter can be performed through fast, multiple evaluations of Taylor polynomials. Test cases with objects ranging from LEO to GEO regimes are considered to assess the performances and the accuracy of the proposed method.     

A Partially Orthogonal EnKF approach to atmospheric density estimation using orbital debris

A key requirement for accurate trajectory prediction and space situational awareness is knowledge of how non-conservative forces affect space object motion. These forces vary temporally and spatially, and are driven by the underlying behavior of space weather particularly in Low Earth Orbit (LEO). Existing trajectory prediction algorithms adjust space weather models based on calibration satellite observations. However, lack of sufficient data and mismodeling of non-conservative forces cause inaccuracies in space object motion prediction, especially for uncontrolled debris objects. The uncontrolled nature of debris objects makes them particularly sensitive to the variations in space weather. Our research takes advantage of this behavior by utilizing observations of debris objects to infer the space environment parameters influencing their motion.The hypothesis of this research is that it is possible to utilize debris objects as passive, indirect sensors of the space environment. We focus on estimating atmospheric density and its spatial variability to allow for more precise prediction of LEO object motion. The estimated density is parameterized as a grid of values, distributed by latitude and local sidereal time over a spherical shell encompassing Earth at a fixed altitude of 400 km. The position and velocity of each debris object are also estimated. A Partially Orthogonal Ensemble Kalman Filter (POEnKF) is used for assimilation of space object measurements to estimate density.For performance comparison, the scenario characteristics (number of objects, measurement cadence, etc.) are based on a sensor tasking campaign executed for the High Accuracy Satellite Drag Model project. The POEnKF analysis details spatial comparisons between the true and estimated density fields, and quantifies the improved accuracy in debris object motion predictions due to more accurate drag force models from density estimates. It is shown that there is an advantage to utilizing multiple debris objects instead of just one object. Although the work presented here explores the POEnKF performance when using information from only 16 debris objects, the research vision is to utilize information from all routinely observed debris objects. Overall, the filter demonstrates the ability to estimate density to within a threshold of accuracy dependent on measurement/sensor error. In the case of a geomagnetic storm, the filter is able to track the storm and provide more accurate density estimates than would be achieved using a simple exponential atmospheric density model or MSIS Atmospheric Model (when calm conditions are assumed).     

一种基于TLE数据的轨道异常分析方法

空间在轨物体的轨道异常是航天工程及预警领域普遍关注的问题,及时发现轨道异常意义重大,通过分析空间物体的轨道异常,可以及时发现和识别规避事件或碰撞事件,还可以了解监测网的能力.本文提出一种基于TLE数据的简单的轨道异常分析方法——长半轴变化法.该方法快速有效,应用到低轨在用卫星和美俄解体碎片的异常分析中,异常物体正确识别率可达到100%;对美俄解体碎片进行轨道异常分析后得出,美国空间监视网可以稳定探测90%以上的解体碎片.      

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.     

Improved orbit predictions using two-line elements

The density of orbital space debris constitutes an increasing environmental challenge. There are two ways to alleviate the problem: debris mitigation and debris removal. This paper addresses collision avoidance, a key aspect of debris mitigation. We describe a method that contributes to achieving a requisite increase in orbit prediction accuracy for objects in the publicly available two-line element (TLE) catalog. Batch least-squares differential correction is applied to the TLEs. Using a high-precision numerical propagator, we fit an orbit to state vectors derived from successive TLEs. We then propagate the fitted orbit further forward in time. These predictions are validated against precision ephemeris data derived from the international laser ranging service (ILRS) for several satellites, including objects in the congested sun-synchronous orbital region. The method leads to a predicted range error that increases at a typical rate of 100 m per day, approximately a 10-fold improvement over individual TLE’s propagated with their associated analytic propagator (SGP4). Corresponding improvements for debris trajectories could potentially provide conjunction analysis sufficiently accurate for an operationally viable collision avoidance system based on TLEs only.     

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.     

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.     

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.     

利用TLE数据分析LEO卫星轨道异常的新方法——————综合判据法

及时准确地发现在轨卫星的轨道异常意义重大. 通过有效的异常算法, 能够找出发生轨道异常的碎片或航天器, 为空间碎片碰撞预警系统分析和验证碰撞事件提供数据支持. 通过对利用TLE (Two Line Elements)数据分析LEO在轨卫星轨道异常的方法研究, 提出了一个利用单个卫星相邻根数时间差控制加综合判据的判别方法. 分析表明, 相对于取单一因素阈值的判别方法, 综合判据法能够最大限度地减少漏判, 并且保持相对较高的判断准确率.      

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.     

Assessment of the close approach frequency and collision probability for satellites in different configurations of large constellations

The Earth orbital environment is drastically changing due to an intensification of the space activities. In particular, several projects of large constellations, proposed for the next years for communications purpose like global internet access, Internet of Things, or for Earth observations, will lead to the deployment of several thousands of new satellites at an unprecedented rate. It is a crucial challenge for space traffic management, which will deal with a great number of satellite conjunctions, potentially causing a collision with damaging consequences for the constellation itself and the space environment sustainability.In this paper, we investigate the close approach frequency and the cumulative collision probability for each referenced constellation. For this purpose, we compute the orbital evolution of satellites in different constellations during the lifecycle, from the deployment to the decommissioning, and we apply the CUBE algorithm and the Foster method to assess the collision probability with the background space debris population assuming a constant uncertainty in position. We show the variation of risk defined by the close approach frequency and the cumulative collision probability as a function of the proposed configuration. In particular, satellites of the Iridium and Kuiper constellation, but also satellite of the Telesat constellation on polar orbits are the most exposed at a collision. Moreover, the decommissioning phase contribute for a major part to the final cumulative collision probability.     

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.     

Assessing and minimizing collisions in satellite mega-constellations

We aim to provide satellite operators and researchers with an efficient means for evaluating and mitigating collision risk during the design process of mega-constellations. We first introduce a novel algorithm for conjunction prediction that relies on large-scale numerical simulations and uses a sequence of filters to greatly reduce its computational expense. We then use this brute-force algorithm to establish baselines of endogenous (intra-constellation), or self-induced, conjunction events for the FCC-reported designs of the OneWeb LEO and SpaceX Starlink mega-constellations. We demonstrate how these deterministic results can be used to validate more computationally efficient, stochastic techniques for close-encounter prediction by adopting a new probabilistic approach from Solar-System dynamics as a simple test case. Finally, we show how our methodology can be applied during the design phase of large constellations by investigating Minimum Space Occupancy (MiSO) orbits, a generalization of classical frozen orbits that holistically account for the perturbed-Keplerian dynamics of the Earth-satellite-Moon-Sun system. The results indicate that the adoption of MiSO orbital configurations of the proposed mega-constellations can significantly reduce the risk of endogenous collisions with nearly indistinguishable adjustments to the nominal orbital elements of the constellation satellites.     

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.     

空间碎片碰撞预警研究

介绍了空间碎片的分布和危害，阐明了航天器进行碎片预警规避的必要性，讨论了碎片碰撞预警技术中的关键问题，用已有的碰撞和规避实例对自行研制的软件进行了计算检验，结果证明了软件的正确性，也说明了国内进行碎片预警工作的可行性．     

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.     

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.     

基于随机点模拟方法分析空间目标危险交会

基于SGP4模型在空间目标轨道预报中的应用, 在预报的位置速度信息和误差 信息基础上, 提出一种空间两目标碰撞预警的分析方法, 即随机点模拟方法. 与传统的交会平面积分方法相比, 其主要有两点不同: 一是在误差信息中考虑 了误差均值的影响, 即误差椭球不再以预报位置为中心分布; 二是在分析方 法上侧重于真实模拟可能的交会情形, 而不忽略任一方向上的误差. 通过算例分 析验证了该方法的可行性, 结果表明误差均值的非零性使得最大碰撞概率不 一定出现在预报的最近交会距离时刻. 同时仿真结果还表明, 两目标在相对速 度方向上的相对位置仍然存在误差, 这可能造成随机点模拟的碰撞概率计算 值较交会平面积分方法偏小. 不同的碰撞预警分析方法对应不同的预警门限, 根据文中实例, 初步确定10-6为随机点模拟方法的红色预警值.      

地气光辐射对空间目标成像特性影响分析

为分析地气光辐射对空间目标成像特性的影响,以地球同步轨道(GEO)卫星搭载的可见光成像器为探测平台,利用卫星工具包(STK)设计高椭圆轨道(HEO)及近地轨道(LEO)目标运动场景,根据空间目标、地球、太阳、探测平台之间的位置关系,采用微元法建立空间目标与地气光背景等效星等模型,推导出空间目标信噪比(SNR)计算公式。分析了距离、角度参数变化对不同轨道空间目标、地气光背景等效星等及空间目标信噪比的影响。仿真结果表明:当探测平台距离空间目标较远时,地气光背景等效星等低于空间目标等效星等,地气光辐射比空间目标信号强。当地气光辐射进入和离开空间目标探测视场时,空间目标信噪比最大,该时间段是进行空间目标探测的最佳“观测窗口”。仿真得出的空间目标信噪比值为空间目标探测识别提供了理论计算依据。