Physical characterization of the deep-space debris WT1190F: A testbed for advanced SSA techniques

We report on extensive ${\mathit{BVR}}_c{I}_c$ photometry and low-resolution ($\lambda /\mathrm{\Delta }\lambda ～250$) spectroscopy of the deep-space debris WT1190F, which impacted Earth offshore from Sri Lanka, on 2015 November 13. In spite of its likely artificial origin (as a relic of some past lunar mission), the case offered important points of discussion for its suggestive connection with the envisaged scenario for a (potentially far more dangerous) natural impactor, like an asteroid or a comet.Our observations indicate for WT1190F an absolute magnitude $R_c={32.45}_{\pm 0.31}$, with a flat dependence of reflectance on the phase angle, such as ${\mathit{dR}}_c/d?～{0.007}_{\pm 2}$?mag?deg^?1. The detected short-timescale variability suggests that the body was likely spinning with a period twice the nominal figure of $P_{\mathrm{flash}}={1.4547}_{\pm 0.0005}\mathrm{s}$, as from the observed lightcurve. In the ${\mathit{BVR}}_c{I}_c$ color domain, WT1190F closely resembled the Planck deep-space probe. This match, together with a depressed reflectance around 4000 and 8500 Å may be suggestive of a “grey” (aluminized) surface texture.The spinning pattern remained in place also along the object fiery entry in the atmosphere, a feature that may have partly shielded the body along its fireball phase perhaps leading a large fraction of its mass to survive intact, now lying underwater along a tight ($～1\times 80$?km) strip of sea, at a depth of 1500?m or less.Under the assumption of Lambertian scatter, an inferred size of ${216}_{\pm 30}/\sqrt{\alpha /0.1}$?cm is obtained for WT1190F. By accounting for non-gravitational dynamical perturbations, the Area-to-Mass ratio of the body was in the range $(0.006?\mathit{AMR}?0.011)$?m²?kg^?1.Both these figures resulted compatible with the two prevailing candidates to WT1190F’s identity, namely the Athena II Trans-Lunar Injection Stage of the Lunar Prospector mission, and the ascent stage of the Apollo 10 lunar module, callsign “Snoopy”. Both candidates have been analyzed in some detail here through accurate 3D CAD design mockup modelling and BRDF reflectance rendering to derive the inherent photometric properties to be compared with the observations.     

On the practical exploitation of perturbative effects in low Earth orbit for space debris mitigation

This paper presents the results of a numerical evaluation of the natural lifetime reduction in low Earth orbit, due to dynamical perturbations. The study considers two values for the area-to-mass ratio, a nominal ratio which resembles a typical value of spacecraft in orbit today, and an enhanced ratio which covers the surface augmentation. The results were obtained with two orbit propagators, one of a semi-analytical nature and the second one using non-averaged equations of motion. The simulations for both propagators were set up similarly to allow comparison. They both use the solar radiation pressure and the secular terms of the geopotential ($J_2,J_4$ and $J_6$). The atmospheric drag was turned on and off in both propagators to alternatively study the eccentricity build up and the residual lifetime. The non-averaging case also covers a validation with the full 6?×?6 geopotential. The results confirm the findings in previous publications, that is, the possibility for de-orbiting from altitudes above the residual atmosphere if a solar sail is deployed at the end-of-life, due to the combined effect of solar radiation pressure and the oblateness of the Earth. At near polar inclinations, shadowing effects can be exploited to the same end. The results obtained with the full, non-averaging propagator revealed additional de-orbiting corridors associated with solar radiation pressure which were not found by previous work on space debris mitigation. The results of both tools are compared for specific initial conditions. For nominal values of area-to-mass ratio, instead, it is confirmed that this resonance effect is negligible.The paper then puts the findings in the perspective of the current satellite catalogue. It identifies space missions which are currently close to a resonance corridor and shows the orbit evolution within the resonances with a significantly shorter residual orbital lifetime. The paper finishes with a discussion on the exploitation of these effects with regards to the long-term simulation of the space debris environment and a flux and collision probability comparison.     

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.     

Satellite collision probability estimation using polynomial chaos expansions

This paper presents the application of polynomial chaos (PC) to estimating the probability of collision between two spacecraft. Common methods of quantifying this probability for conjunction assessment use either Monte Carlo analyses or include simplifying assumptions to improve tractability. A PC expansion, or PCE, provides a means for approximating the solution to a large set of stochastic ordinary differential equations, which includes orbit propagation. When compared to Monte Carlo methods, non-intrusive, i.e., sampling-based, PCE generation techniques may greatly reduce the number of orbit propagations required to approximate the possibly non-Gaussian, a posteriori probability density function. The presented PC-based method of computing collision probability requires no fundamental simplifying assumptions, and reduces the computation time compared to Monte Carlo. This paper considers two cases where the common conjunction assessment assumptions are no longer valid. The results indeed demonstrate a reduction in computation time when compared to Monte Carlo, and improved accuracy when compared to simplified techniques.     

Estimating the spin axis orientation of the Echostar-2 box-wing geosynchronous satellite

For the first time, the spin axis orientation of an inactive box-wing geosynchronous satellite has been estimated from ground-based optical photometric observations of Echostar-2’s specular reflections. Recent photometric light curves obtained of Echostar-2 over four years suggest that unusually bright and brief specular reflections were occurring twice within an observed spin period. These bright and brief specular reflections suggested two satellite surfaces with surface normals separated by approximately 180°. The geometry between the satellite, the Sun, and the observing location at the time of each of the brightest observed reflections, was used to estimate Echostar-2’s equatorial spin axis orientation coordinates. When considering prograde and retrograde rotation, Echostar-2’s spin axis orientation was estimated to have been located within 30° of either equatorial coordinate pole. Echostar-2’s spin axis was observed to have moved approximately 180° in right ascension, within a time span of six months, suggesting a roughly one year spin axis precession period about the satellite’s angular momentum vector.     

An optimization approach for observation association with systemic uncertainty applied to electro-optical systems

The observation to observation measurement association problem for dynamical systems can be addressed by determining if the uncertain admissible regions produced from each observation have one or more points of intersection in state space. An observation association method is developed which uses an optimization based approach to identify local Mahalanobis distance minima in state space between two uncertain admissible regions. A binary hypothesis test with a selected false alarm rate is used to assess the probability that an intersection exists at the point(s) of minimum distance. The systemic uncertainties, such as measurement uncertainties, timing errors, and other parameter errors, define a distribution about a state estimate located at the local Mahalanobis distance minima. If local minima do not exist, then the observations are not associated. The proposed method utilizes an optimization approach defined on a reduced dimension state space to reduce the computational load of the algorithm. The efficacy and efficiency of the proposed method is demonstrated on observation data collected from the Georgia Tech Space Object Research Telescope.     

空间碎片自然衰减和太阳活动

分析研究了空间碎片数随太阳辐射流量Ｆ１０．７的变化;给出预报Ｆ１０．７长期变化的计算方法和预测空间碎片数的数学模型。结果显示：①强太阳活动造成空间碎片年增长率下降;②空间碎片数与太阳活动１１年变化密切相关,相关数为０．９;③空间碎片增长率约为发射率的两倍;④若发射率保持不变,则到２０２０年,大于１０ｃｍ的碎片数将达到１４５００;⑤若小碎片的增长为大碎片增长的两倍,则到２０２０年,大于１ｃｍ的碎片数可达１２５０００。     

Comparison of the transport codes HZETRN,HETC and FLUKA for a solar particle event

The protection of astronauts and instrumentation from galactic cosmic rays and solar particle events is one of the primary constraints associated with mission planning in low earth orbit or deep space. To help satisfy this constraint, several computational tools have been developed to analyze the effectiveness of various shielding materials and structures exposed to space radiation. These tools are now being carefully scrutinized through a systematic effort of verification, validation, and uncertainty quantification. In this benchmark study, the deterministic transport code HZETRN is compared to the Monte Carlo transport codes HETC-HEDS and FLUKA for a 30 g/cm² water target protected by a 20 g/cm² aluminum shield exposed to a parameterization of the February 1956 solar particle event. Neutron and proton fluences as well as dose and dose equivalent are compared at various depths in the water target. The regions of agreement and disagreement between the three codes are quantified and discussed, and recommendations for future work are given.     

地球同步轨道长寿命卫星热控涂层太阳吸收率性能退化研究

文章介绍了15年地球同步轨道环境对卫星表面太阳吸收率性能影响的试验模拟研究,研究结果为长寿命卫星热设计及热控涂层选择和研制提供可靠依据;介绍自行研制的空间低能综合环境试验设备、太阳吸收率原位测试系统和空间低能综合环境模拟试验方法,并对航天器常用的S781白漆、SR107ZK白漆、F46镀银和OSR二次表面镜热控材料进行空间低能综合环境模拟试验,获得了这些热控涂层在地球同步轨道15年期间太阳吸收率性能退化模拟数据,与已有的飞行试验数据进行对比研究,取得满意的结果。     

Perspective on space radiation for space flights in 2020–2040

The Sun undergoes several well known periodicities in activity, such as the Schwabe 11 year cycle, the Gleissberg 80–90 year cycle, the Suess 200–210 year cycle and the Halstatt 2200–2300 year cycle. In addition, there is evidence that the 20th century levels of solar activity are unusually high. The years 2020–2040 are expected to coincide with increased activity in human space flight beyond low Earth orbit. The solar cycles and the present level of solar activity are reviewed and their activities during the years 2020–2040 are discussed with a perspective on space radiation and the future program of space flight. It is prudent to prepare for continuing levels of high solar activity as well as for the low levels of the current deep minimum, which has corresponded to high galactic cosmic ray flux.     

International scientific optical network for space debris research

A joint team of researchers under the auspices of the Center for Space Debris Information Collection, Processing and Analysis of the Russian Academy of Sciences collaborates with 15 observatories around the world to perform observations of space debris. For this purpose, 14 telescopes were equipped with charge-coupled device (CCD) cameras, Global Positioning System (GPS) receivers, CCD frame processing and ephemeris computation software, with the support of the European and Russian grants. Many of the observation campaigns were carried out in collaboration with the Astronomical Institute of the University of Bern (AIUB) team operating at the Zimmerwald observatory and conducting research for the European Space Agency (ESA), using the Tenerife/Teide telescope for searching and tracking of unknown objects in the geostationary region (GEO). More than 130,000 measurements of space objects along a GEO arc of 340.9°, collected and processed at Space Debris Data Base in the Ballistic Center of the Keldysh Institute of Applied Mathematics (KIAM) in 2005–2006, allowed us to find 288 GEO objects that are absent in the public orbital databases and to determine their orbital elements. Methods of discovering and tracking small space debris fragments at high orbits were developed and tested. About 40 of 150 detected unknown objects of magnitudes 15–20.5 were tracked during many months. A series of dedicated 22-cm telescopes with large field of view for GEO survey tasks is in process of construction. 7 60-cm telescopes will be modernized in 2007.     

Short-arc tracklet association for geostationary objects

Measurement association and initial orbit determination is a fundamental task when building up a database of space objects. This paper proposes an efficient and robust method to determine the orbit using the available information of two tracklets, i.e. their line-of-sights and their derivatives. The approach works with a boundary-value formulation to represent hypothesized orbital states and uses an optimization scheme to find the best fitting orbits. The method is assessed and compared to an initial-value formulation using a measurement set taken by the Zimmerwald Small Aperture Robotic Telescope of the Astronomical Institute at the University of Bern. False associations of closely spaced objects on similar orbits cannot be completely eliminated due to the short duration of the measurement arcs. However, the presented approach uses the available information optimally and the overall association performance and robustness is very promising. The boundary-value optimization takes only around 2% of computational time when compared to optimization approaches using an initial-value formulation. The full potential of the method in terms of run-time is additionally illustrated by comparing it to other published association methods.     

Microscopic approach to analyze solar-sail space-environment effects

Near-sun space-environment effects on metallic thin films solar sails as well as hollow-body sails with inflation fill gas are considered. Analysis of the interaction of the solar radiation with the solar-sail materials is presented. This analysis evaluates worst-case solar radiation effects during solar-radiation-pressure acceleration. The dependence of the thickness of solar sail on temperature and on wavelength of the electromagnetic spectrum of solar radiation is investigated. Physical processes of the interactions of photons, electrons, protons and α-particles with sail material atoms and nuclei, and inflation fill gas molecules are analyzed. Calculations utilized conservative assumptions with the highest values for the available cross sections for interactions of solar photons, electrons and protons with atoms, nuclei and hydrogen molecules. It is shown that for high-energy photons, electrons and protons the beryllium sail is mostly transparent. Sail material will be partially ionized by solar UV and low-energy solar electrons. For a hollow-body photon sail effects including hydrogen diffusion through the solar-sail walls, and electrostatic pressure is considered. Electrostatic pressure caused by the electrically charged sail’s electric field may require mitigation since sail material tensile strength decreases with elevated temperature. It also can substitute inflation-gas pressure loss due to gas diffusion and perforation by micrometeoroids impact to keep the sail inflated.     

数据同化在空间天气学中的应用

由太阳活动引起的耀斑和日冕物质抛射等短时间尺度变化的空间天气事件会影响并危害地球磁层、电离层、中高层大气、卫星运行安全以及人类健康,因此对这些空间天气事件的预测显得尤为重要。数据同化在稀疏观测和异步采集的情况下能够增加模型的预测能力,对模型变量进行自洽分析。在数值预报中引入数据同化方法,能够提高预测可信度。本文从数据同化方法的角度出发,主要分析了数据同化目前在大气、电离层、磁层、太阳及其他行星科学研究中的应用,并初步讨论了数据同化未来在空间天气方面的应用。