Analytical and semi-analytical investigations of geosynchronous space debris with high area-to-mass ratios

This paper provides a Hamiltonian formulation of the averaged equations of motion with respect to short periods (1 day) of a space debris subjected to direct solar radiation pressure and orbiting near the geostationary ring. This theory is based on a semi-analytical theory of order 1 regarding the averaging process, 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, truncated at an arbitrary high order.     

大面质比空间碎片在太阳光压和引力作用下的轨道演化

分析了较高轨道(a > 10000km)大面质比空间碎片的轨道动力学演化问题. 重点讨论了位于地球同步轨道的空间碎片轨道演化问题, 并给出轨道偏心率 随时间演化的表达式. 通过进一步分析得出, 倾角大于63°26'的GTO轨 道空间碎片, 仅在J₂和第三体摄动影响下, 会出现轨道偏心率升高; 而对 于大面质比空间碎片, 在J₂项和太阳光压同时作用下, 当近地点指向的角 变率与太阳平黄经变化率接近时, 会出现长期共振现象, 导致轨道偏心率升 高, 近地点降低. 分析还得出, 轨道演化过程中, 偏心率的最大值与初始轨 道近地点的指向有关.      

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.     

The classical Laplace plane as a stable disposal orbit for geostationary satellites

The classical Laplace plane is a frozen orbit, or equilibrium solution for the averaged dynamics arising from Earth oblateness and lunisolar gravitational perturbations. The pole of the orbital plane of uncontrolled GEO satellites regress around the pole of the Laplace plane at nearly constant inclination and rate. In accordance with Friesen et al. (1993), we show how this stable plane can be used as a robust long-term disposal orbit. The current graveyard regions for end-of-life retirement of GEO payloads, which is several hundred kilometers above GEO depending on the spacecraft characteristics, cannot contain the newly discovered high area-to-mass ratio debris population. Such objects are highly susceptible to the effects of solar radiation pressure exhibiting dramatic variations in eccentricity and inclination over short periods of time. The Laplace plane graveyard, on the contrary, would trap this debris and would not allow these objects to rain down through GEO. Since placing a satellite in this inclined orbit can be expensive, we discuss some alternative disposal schemes that have acceptable cost-to-benefit ratios.     

Space debris mitigation strategies and practices in geosynchronous transfer orbits

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

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

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.     

基于相控阵雷达波束篱笆的空间碎片数量与分布估计方法

随着载人航天与空间站等航天活动的增多,不能有效防护、也无法定期跟踪和编目的小尺寸(尤其是1～10 cm)碎片的危害越来越受到关注,这些碎片信息的获取依赖于统计采样技术.针对简化的相控阵雷达波束篱笆空间碎片探测模式,提出了一种采用统计技术估计空间碎片总数量以及高度和倾角分布的方法.将碎片穿越波束篱笆的过程用Poisson分布来建模,根据观测时段内穿越波束篱笆目标的平均到达率及测量的轨道高度和倾角数据来估计给定轨道高度范围或倾角范围内碎片的数量,进而得到碎片的总数量以及碎片数量随轨道高度或倾角的分布.在获取雷达散射截面信息时,该方法还可用于估计碎片数量随尺寸的分布.通过仿真实验验证了该方法的有效性.      

Physical properties and long-term evolution of the debris clouds produced by two catastrophic collisions in Earth orbit

The population of cataloged orbital debris increased by approximately 40% in just a couple of years, from January 2007 to February 2009. This was due to two collisions in space, which involved the catastrophic destruction of three intact satellites (Fengyun 1C, Cosmos 2251 and Iridium 33) in high inclination orbits. Both events occurred in the altitude range already most affected by previous launch activity and breakup events, thus boosting the cataloged population in low Earth orbit by more than 60%.     

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.     

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.     

太阳中纬度探测共振借力轨道设计

针对太阳中纬度探测任务,设计了共振借力轨道方案.根据离开超越速度的解析表达式得到借力后日心轨道的可行集并将结论推广到任意借力行星.进而提出RRpInc图形化设计方法,通过RRpInc图分析得到共振借力的超越速度大小直接决定了最终轨道倾角.以共振借力行星作为目标星体设计借力转移轨道,以逃逸速度和到达超速作为性能指标,采用多目标优化算法优化.给出了两个设计方案:方案1选择金星作为共振借力行星,探测器通过"金星-地球-地球"借力飞行增加金星共振借力前的超越速度;方案2选择地球为共振借力行星,探测器在木星借力后返回地球共振借力.两个方案最终轨道倾角分别约为30°和35°.     

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.     

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

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

Impact of the orbital eccentricity on the attitude performance before and after the deorbiting phase for Alsat-1

Due to the presence of periodic forcing terms in the gravity gradient torque, orbit eccentricity may produce large response for the roll, yaw and pitch angles. This paper investigates the influence of the orbit eccentricity on the performance of the attitude determination and control subsystem (ADCS) pointing of passive Low Earth Orbit (LEO) satellites stabilized by a gravity gradient boom or having long appendages before and after the deorbiting operation. The contribution of this work is twofold. First, the satellite attitude dynamics and kinematics are modeled by introducing the orbit eccentricity in the equations of motion of a LEO satellite in order to provide the best scenario in which satellite operators can keep the nominal functionality of LEO satellites with a gravity gradient boom after the deorbiting operation. Second, a Quaternion-based Extended Kalman Filter (EKF) is analyzed when the orbit eccentricity is considered in order to determine the influence of this disturbance on the convergence and stability of the filter. The simulations in this work are based on the true parameters of Alsat-1 which is a typical LEO satellite stabilized by a gravity gradient boom. The results show that the orbit eccentricity has a big influence on the pointing system accuracy causing micro-vibrations that affect the geocentric pointing particularly after the deorbiting phase. In this case, satellites have no orbital correction option. The Quaternion-based Extended Kalman Filter analyzed in this paper, achieved satisfactory results for eccentricity values less than 0.4 with respect to pointing system accuracy. However, singularities were observed for eccentricity values greater than 0.4.     

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.     

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.     

Dynamical cartography of Earth satellite orbits

We have carried out a numerical investigation of the coupled gravitational and non-gravitational perturbations acting on Earth satellite orbits in an extensive grid, covering the whole circumterrestrial space, using an appropriately modified version of the SWIFT symplectic integrator, which is suitable for long-term ($～$120?years) integrations of the non-averaged equations of motion. Hence, we characterize the long-term dynamics and the phase-space structure of the Earth-orbiter environment, starting from low altitudes ($～$400?km) and going up to the GEO region and beyond. This investigation was done in the framework of the EC-funded “ReDSHIFT” project, with the purpose of enabling the definition of passive debris removal strategies, based on the use of physical mechanisms inherent in the complex dynamics of the problem (i.e., resonances). Accordingly, the complicated interactions among resonances, generated by different perturbing forces (i.e., lunisolar gravity, solar radiation pressure, tesseral harmonics in the geopotential) are accurately depicted in our results, where we can identify the regions of phase space where the motion is regular and long-term stable and regions for which eccentricity growth and even instability due to chaotic behavior can emerge. The results are presented in an “atlas” of dynamical stability maps for different orbital zones, with a particular focus on the (drag-free) range of semimajor axes, where the perturbing effects of the Earth’s oblateness and lunisolar gravity are of comparable order. In some regions, the overlapping of the predominant lunisolar secular and semi-secular resonances furnish a number of interesting disposal hatches at moderate to low eccentricity orbits. All computations were repeated for an increased area-to-mass ratio, simulating the case of a satellite equipped with an on-board, area-augmenting device. We find that this would generally promote the deorbiting process, particularly at the transition region between LEO and MEO. Although direct reentry from very low eccentricities is very unlikely in most cases of interest, we find that a modest “delta-v” ($\mathrm{\Delta }V$) budget would be enough for satellites to be steered into a relatively short-lived resonance and achieve reentry into the Earth’s atmosphere within reasonable timescales ($～$50?years).     

Influence of the 2:1 resonance in the orbits of GPS satellites

Due to the characteristics of their orbits the GPS satellites are submitted to the following main perturbations: terrestrial gravitational field, luni-solar gravitational attraction and solar radiation pressure (including the effects of the Earth's shadow). An additional perturbation arises due to the 2:1 commensurability of the orbital period of the satellite with the period of the Earth's rotation. An analytical theory is briefly presented to solve the equations of motion including the previously mentioned effects. The analytical solution, based on the Lie-Hori method, is compared with a numerical integration of the equations.