空间碎片的空间密度算法研究

由于空间碎片的影响, 空间环境日益恶劣,有必要建立空间碎片环境工程模型对空间碎片撞击航天器进行风险评估.本文研究了空间碎片环境模型中的重要环节, 即碎片的空间密度问题,在统计理论与椭圆轨道理论基础上,分析及推导了空间碎片在空间中的分布状况及空间密度,并得出了碎片空间密度的空间坐标函数解析表达式.同时对引入假设条件的合理性进行了讨论,并利用双行元数据对结论进行了验证, 利用本文方法得出的结果与双行元数据吻合.      

基于期望最大算法的空间事件及异常值探测

美国ELSET数据库提供的TLE数据是目前使用最广泛的数据,在热层大气密度反演、弹道系数估计、碰撞预警等领域具有重要作用。受空间环境扰动、空间事件以及TLE产生过程等共同影响,ELSET数据库包含大量亟待清理的异常值和识别的空间事件,例如发布错误的TLE、轨道根数异常和Bstar异常。现有方法在清理异常轨道根数时缺乏统一性,需要使用不同的技术,清理流程较为繁杂,并且仅适用于特定轨道区域的少数目标。为克服现有方法的弊端,提出了一种基于期望最大算法的滑动窗口–多项式拟合预报方法,对含有轨道机动的碎片以及受空间环境影响的碎片进行异常值与空间事件探测。研究表明,该方法能够灵活处理不同空间环境下的异常值与空间事件探测,具有普适性,适用于所有轨道碎片。      

Debris flux comparisons from the Goldstone Radar,Haystack Radar,and Hax Radar prior,during, and after the last solar maximum

The continual monitoring of the low Earth orbit (LEO) debris environment using highly sensitive radars is essential for an accurate characterization of these dynamic populations. Debris populations are continually evolving since there are new debris sources, previously unrecognized debris sources, and debris loss mechanisms that are dependent on the dynamic space environment. Such radar data are used to supplement, update, and validate existing orbital debris models. NASA has been utilizing radar observations of the debris environment for over a decade from three complementary radars: the NASA JPL Goldstone radar, the MIT Lincoln Laboratory (MIT/LL) Long Range Imaging Radar (known as the Haystack radar), and the MIT/LL Haystack Auxiliary radar (HAX). All of these systems are highly sensitive radars that operate in a fixed staring mode to statistically sample orbital debris in the LEO environment. Each of these radars is ideally suited to measure debris within a specific size region. The Goldstone radar generally observes objects with sizes from 2 mm to 1 cm. The Haystack radar generally measures from 5 mm to several meters. The HAX radar generally measures from 2 cm to several meters. These overlapping size regions allow a continuous measurement of cumulative debris flux versus diameter from 2 mm to several meters for a given altitude window. This is demonstrated for all three radars by comparing the debris flux versus diameter over 200 km altitude windows for 3 nonconsecutive years from 1998 to 2003. These years correspond to periods before, during, and after the peak of the last solar cycle. Comparing the year to year flux from Haystack for each of these altitude regions indicate statistically significant changes in subsets of the debris populations. Potential causes of these changes are discussed. These analysis results include error bars that represent statistical sampling errors.     

Space debris cumulative flux considering the Interval Distance-based method

The availability of engineering models to estimate the risk from space debris is essential for space missions. According to current research, cumulative flux calculation is mostly carried out based on the equal-width interval discretization. The method discretizes the volume around the Earth into cells defined in earth centered inertial coordinates. The resulting debris flux onto a target object is shown to depend on the chosen size of the cells. To avoid a discretization error, this must be accounted for. In order to present reliable flux predictions for space mission, the algorithm improvement is an ongoing topic for the related research field. The aim of this study was to examine the discretization error during the cumulative flux determination process. Both the effect of interval step length and the orbital boundary are under investigation. Several typical orbits are selected as examples here and the 2018/01/03 TLE (Two Line Element) data published by the US Space Surveillance Network is used as the debris background in this paper. Furthermore, the Interval Distance-Based method for Discretization (IDD) is adopted in this paper. A position-centered flux determination method is introduced based on the IDD method. According to the example analysis, the IDD used in the flux calculation process provides results which are less affected by the interval step-size setup; and the orbital boundary has no effect on the calculation process. In other words, the discretization error is significantly reduced. The position-centered method provided a possible suggestion for the improvement of space debris environment models.     

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.     

Atmospheric drag modelling for orbital micro-debris at LEO altitudes

Data from satellite impact experiments and the scanning of recovered spacecraft offers an extended timebase to examine, using a consistent methodology, the microparticle fluxes. New penetration data from the TiCCE experiment on Eureca /1, 2/ adds to this database and shows that - despite an expected growth in the micro-debris flux - the observed flux is not greater than either LDEF or SMM. The question arises: “is this consistent with the micro particle flux being dominated by space debris or by meteoroids”.

To assist this assessment, numerical modelling using the Gear method /3/ of explicit time integration of the atmospheric drag lifetime of micron dimensioned orbital debris in both circular (LEO) and eccentric (GTO) orbits has been performed for the relevant space exposures. Results are applied to the data to examine whether the recent variations in flux can be attributed to varying levels of, orbital micro-debris caused by atmospheric drag and its changes during the solar cycle.     

The influence of solid rocket motor retro-burns on the space debris environment

The ESA space debris population model MASTER (Meteoroid and Space Debris Terrestrial Environment Reference) considers firings of solid rocket motors (SRM) as a debris source with the associated generation of slag and dust particles. The resulting slag and dust population is a major contribution to the sub-millimetre size debris environment in Earth orbit. The current model version, MASTER-2005, is based on the simulation of 1076 orbital SRM firings which contributed to the long-term debris environment. A comparison of the modelled flux with impact data from returned surfaces shows that the shape and quantity of the modelled SRM dust distribution matches that of recent Hubble Space Telescope (HST) solar array measurements very well. However, the absolute flux level for dust is under-predicted for some of the analysed Long Duration Exposure Facility (LDEF) surfaces. This points into the direction of some past SRM firings not included in the current event database. The most suitable candidates for these firings are the large number of SRM retro-burns of return capsules. Objects released by those firings have highly eccentric orbits with perigees in the lower regions of the atmosphere. Thus, they produce no long-term effect on the debris environment. However, a large number of those firings during the on-orbit time frame of LDEF might lead to an increase of the dust population for some of the LDEF surfaces. In this paper, the influence of SRM retro-burns on the short- and long-term debris environment is analysed. The existing firing database is updated with gathered information of some 800 Russian retro-firings. Each firing is simulated with the MASTER population generation module. The resulting population is compared against the existing background population of SRM slag and dust particles in terms of spatial density and flux predictions.     

Out-of-ecliptic distribution of interplanetary dust derived from near earth flux

The out-of-ecliptic distribution of interplanetary dust, i.e. its number density, mainly has been subject to optical or infrared remote sensing techniques. As the population in interplanetary space is made up of orbiting particles which will cross the ecliptic plane, determination of their orbital properties there gives a possibility also to derive their out-of-ecliptic distribution. Determination of orbital elements is provided by advanced detectors capable of measuring the vector of impact velocity. In a simple model, which applies for advanced detectors in near earth orbit, the feasibility of the method to determine the out-of-ecliptic spatial distribution of dust has been tested.     

The new ESA meteoroid model

The orbital distributions of meteoroids in interplanetary space are revised in the ESA meteoroid model to account for recently obtained observational data and to comply with the constraints due to the orbital evolution under planetary gravity and Poynting–Robertson effects. Infrared observations of the zodiacal cloud by the COBE DIRBE instrument, in situ flux measurements by the dust detectors on board Galileo and Ulysses spacecraft, and the crater size distributions on lunar rock samples retrieved by the Apollo missions are synthesized into a single model. Within the model, the orbital distributions are expanded into a sum of contributions due to a number of known sources, including the asteroid belt with the emphasis on the prominent families Themis, Koronis, Eos and Veritas, as well as comets on Jupiter-encountering orbits. An attempt to incorporate the meteor orbit database acquired by the AMOR radar is also discussed.     

Radiation measurements inside a human phantom aboard the International Space Station using Liulin-5 charged particle telescope

Liulin-5 is a particle telescope developed for the investigation of the radiation environment within the Russian spherical tissue-equivalent phantom on the International Space Station (ISS). Liulin-5 experiment is conducted aboard the Russian segment of ISS since 28 June 2007 as an adherent part of the international project MATROSHKA-R. The main objective of Liulin-5 experiment is to study the depth-dose distribution of the different components of the orbital radiation field in a human phantom. Additional objectives are mapping of the radiation environment in the phantom and its variations with time and orbital parameters (such as solar cycle, solar flare events, inclination and altitude). Liulin-5 is an active instrument, capable to provide real-time radiation data for the particle flux and dose rates, energy deposition and LET spectra. Data are recorded automatically on memory cards, periodically transported to ground by returning vehicles. In this report we present some first results from data analysis including energy deposition spectra, absorbed dose, dose rate and flux distribution measured simultaneously at 3 different depths of phantom’s radial channel and linear energy transfer (LET) spectrum. Data discussed are for the period July 2007–April 2008.     

WMAP extragalactic point sources as potential Space VLBI calibrators

The point source list of the Wilkinson Microwave Anisotropy Probe (WMAP) is a uniform, all-sky catalogue of bright sources with flux density measurements at high (up to 94 GHz) radio frequencies. We investigated the five-year WMAP list to compile a new catalogue of bright and compact extragalactic radio sources to be potentially studied with Very Long Baseline Interferometry at millimeter wavelengths (mm-VLBI) and Space VLBI (SVLBI). After comparing the WMAP data with the existing mm-VLBI catalogues, we sorted out the yet unexplored sources. Using the 41, 61 and 94 GHz WMAP flux densities, we calculated the spectral indices. By collecting optical identifications, lower-frequency radio flux densities and VLBI images from the literature, we created a list of objects which have not been investigated with VLBI at 86 GHz before. With total flux density at least 1 Jy and declination above −40°, we found 37 suitable new targets. It is a nearly 25% addition to the known mm-VLBI sources. Such objects are also potentially useful as phase-reference calibrators for the future Japanese SVLBI mission ASTRO-G at its highest observing frequency (43 GHz). The phase-referencing capability of ASTRO-G would allow long integrations and hence better sensitivity for observing faint target sources close to suitable phase calibrators in the sky.     

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

利用卫星两行轨道根数反演热层密度

两行轨道根数（TLEs）是基于一般摄动理论产生的用于预报地球轨道飞行器位置和速度的一组轨道参数,通过求解大气阻力微分方程,可反演出热层大气密度. 本文选取近圆轨道CHAMP卫星和椭圆轨道Explorer8卫星,以两行轨道根数数据为基础,计算反弹道系数,并根据不同轨道特征采用两种不同反演方法对热层大气密度进行研究. 结果表明,这两种方法反演得到的大气密度与实测值均符合较好,其中CHAMP卫星的反演结果和经验模式值相对于实测值的误差分别为7.94%和13.94%,Explorer8卫星的误差分别为9.04%和14.32%. 相比模式值,利用两行轨道根数数据反演的热层大气密度更接近于实测值,说明该方法可以作为获取大量可靠大气密度数据的一种有效途径.      

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.     

A material experiment for small satellites to characterise the behaviour of carbon nanotubes in space – development and ground validation

Over the last years, Carbon Nanotubes (CNT) drew interdisciplinary attention. Regarding space technologies a variety of potential applications were proposed and investigated. However, no complex data on the behaviour and degradation process of carbon nanotubes under space environment exist. Therefore, it is necessary to investigate the performance of these new materials in space environment and to revaluate the application potential of CNTs in space technologies.Hence, CiREX (Carbon Nanotubes – Resistance Experiment) was developed as a part of a student project. It is a small and compact experiment, which is designed for CubeSat class space satellites. These are a class of nanosatellites with a standardized size and shape. The CiREX design, electrical measurements and the satellites interfaces will be discussed in detail. CiREX is the first in-situ space material experiment for CNTs.To evaluate the data obtained from CiREX, ground validation tests are mandatory. As part of an extensive test series the behaviour of CNTs under solar ultra violet light (UV) and vacuum ultraviolet light (VUV) was examined. Single-walled carbon nanotubes (SWNT), multi-walled carbon nanotubes (MWNT) and MWNT/resin composite (ME) were exposed to different light sources. After the exposure, the defect density was investigated with Raman spectroscopy. There is a clear indication that UV and VUV light can increase the defect density of untreated CNTs and influence the electrical behaviour.     

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.     

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

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

Local orbital debris flux study in the geostationary ring

A local orbital debris flux analysis is performed in the geostationary (GEO) ring to investigate how frequently near-miss events occur for each longitude slot in the GEO ring. The current resident space object (RSO) environment at GEO is evaluated, and publicly-available two-line element (TLE) data are utilized in tandem with a geostationary torus configuration to simulate near-miss events incurred by the trackable RSO population at GEO. Methodology for determining near-miss events with this formulation is introduced, and the results of the analysis for a one-year time frame are provided to illustrate the need for active GEO remediation.     

Measuring atmospheric density using GPS–LEO tracking data

We present a method to estimate the total neutral atmospheric density from precise orbit determination of Low Earth Orbit (LEO) satellites. We derive the total atmospheric density by determining the drag force acting on the LEOs through centimeter-level reduced-dynamic precise orbit determination (POD) using onboard Global Positioning System (GPS) tracking data. The precision of the estimated drag accelerations is assessed using various metrics, including differences between estimated along-track accelerations from consecutive 30-h POD solutions which overlap by 6 h, comparison of the resulting accelerations with accelerometer measurements, and comparison against an existing atmospheric density model, DTM-2000. We apply the method to GPS tracking data from CHAMP, GRACE, SAC-C, Jason-2, TerraSAR-X and COSMIC satellites, spanning 12 years (2001–2012) and covering orbital heights from 400 km to 1300 km. Errors in the estimates, including those introduced by deficiencies in other modeled forces (such as solar radiation pressure and Earth radiation pressure), are evaluated and the signal and noise levels for each satellite are analyzed. The estimated density data from CHAMP, GRACE, SAC-C and TerraSAR-X are identified as having high signal and low noise levels. These data all have high correlations with anominal atmospheric density model and show common features in relative residuals with respect to the nominal model in related parameter space. On the contrary, the estimated density data from COSMIC and Jason-2 show errors larger than the actual signal at corresponding altitudes thus having little practical value for this study. The results demonstrate that this method is applicable to data from a variety of missions and can provide useful total neutral density measurements for atmospheric study up to altitude as high as 715 km, with precision and resolution between those derived from traditional special orbital perturbation analysis and those obtained from onboard accelerometers.