Collision Risk Mitigation in Geostationary Orbit

The short- and long-term effects of spacecraft explosions, as a function of the end-of-life re-orbit altitude above the geostationary orbit (GEO), were analyzed in terms of their additional contribution to the debris flux in the GEO ring. The simulated debris clouds were propagated for 72yrs, taking into account all the relevant orbital perturbations.The results obtained show that 6–7 additional explosions in GEO would be sufficient, in the long term, to double the current collision risk with sizable objects in GEO. Unfortunately, even if spacecraft were to re-orbit between 300 and 500km above GEO, this would not significantly improve the situation. In fact, an altitude increase of at least 2000km would have to be adopted to reduce by one order of magnitude the long-term risk of collision among geostationary satellites and explosion fragments. The optimal debris mitigation strategy should be a compromise between the reliability and effectiveness of spacecraft end-of-life passivation, the re-orbit altitude and the acceptable debris background in the GEO ring. However, for as long as the re-orbit altitudes currently used are less than 500km above GEO, new spacecraft explosions must be avoided in order to preserve the geostationary environment over the long term.     

Track-While-Scan Algorithm in a Clutter Environment

A track-while scan (TWS) algorithm is developed for targets in a clutter environment. The problem has been studied using only the position measurements [1, 5-8], but the simulation results have not been satisfactory. Modern processing techniques (FFT processor) ) in air traffic control and surveillance radar receivers provide both position and radial velocity. The radial velocity measurement may be conveniently used in the target-track correlation process, which will reduce the association ambiguity in the clutter environment. t. In the clear environment the algorithm using the position and radial velocity measurements has been treated in [3, 4]. A TWS algorithm, using both position and radial velocity measurements for targets in a clutter environment, is presented here. The algorithm obtained is nonlinear and adaptive. In order to evaluate the improvement due to radial velocity measurement a simulation has been performed on a digital computer. The algorithm was run with and without radial velocity measurements to compare its performances. An improvement was noted especially when the target path included an accelerated portion.     

Dynamical evolution of high area-to-mass ratio debris released into GPS orbits

A large set of simulations, including all the relevant perturbations, was carried out to investigate the long-term dynamical evolution of fictitious high area-to-mass ratio (A/M) objects released, with a negligible velocity variation, in each of the six orbital planes used by Global Positioning System (GPS) satellites. As with similar objects discovered in near synchronous trajectories, long lifetime orbits, with mean motions of about 2 revolutions per day, were found possible for debris characterized by extremely high area-to-mass ratios. Often the lifetime exceeds 100 years up to A/M ∼ 45 m²/kg, decreasing rapidly to a few months above such a threshold. However, the details of the evolution, which are conditioned by the complex interplay of solar radiation pressure and geopotential plus luni-solar resonances, depend on the initial conditions. Different behaviors are thus possible. In any case, objects like those discovered in synchronous orbits, with A/M as high as 20–40 m²/kg, could also survive in this orbital regime, with semi-major axes close to the semi-synchronous values, with maximum eccentricities between 0.3 and 0.7, and with significant orbit pole precessions (faster and wider for increasing values of A/M), leading to inclinations between 30° and more than 90°.     

Benefits and risks of using electrodynamic tethers to de-orbit spacecraft

By using electrodynamic drag to greatly increase the orbital decay rate, an electrodynamic space tether can remove spent or dysfunctional spacecraft from low Earth orbit (LEO) rapidly and safely. Moreover, the low mass requirements of such tether devices make them highly advantageous compared to conventional rocket-based de-orbit systems. However, a tether system is much more vulnerable to space debris impacts than a typical spacecraft and its design must be proved to be safe up to a certain confidence level before being adopted for potential applications. To assess space debris related concerns, in March 2001 a new task (Action Item 19.1) on the “Potential Benefits and Risks of Using Electrodynamic Tethers for End-of-life De-orbit of LEO Spacecraft” was defined by the Inter-Agency Space Debris Coordination Committee (IADC). Two tests were proposed to compute the fatal impact rate of meteoroids and orbital debris on space tethers in circular orbits, at different altitudes and inclinations, as a function of the tether diameter to assess the survival probability of an electrodynamic tether system during typical de-orbiting missions. IADC members from three agencies, the Italian Space Agency (ASI), the Japan Aerospace Exploration Agency (JAXA) and the US National Aeronautics and Space Administration (NASA), participated in the study and different computational approaches were specifically developed within the framework of the IADC task. This paper summarizes the content of the IADC AI 19.1 Final Report. In particular, it introduces the potential benefits and risks of using tethers in space, it describes the assumptions made in the study plan, it compares and discusses the results obtained by ASI, JAXA and NASA for the two tests proposed. Some general conclusions and recommendations are finally extrapolated from this massive and intensive piece of research.     

Impact of the time span selected to calibrate the ballistic parameter on spacecraft re-entry predictions

The trajectory modeling of satellites that are re-entering the Earth’s atmosphere, as a result of natural orbital decay, has always been a challenging task. Residual lifetime estimations and re-entry predictions are affected by substantial uncertainties, associated with atmospheric density models, with the forecasts of the relevant solar and geomagnetic activity indices and with tracking data, which for uncontrolled re-entries are usually sparse and not particularly accurate. Furthermore, modeling the aerodynamic forces that act on low altitude satellites is a formidable task, especially for objects of a complex shape and unknown attitude evolution.     

Orbital evolution of the first upper stages used for the new European and Chinese navigation satellite systems

The solutions adopted for the disposal of the upper stages used to put in orbit the first satellites of the new European (Galileo) and Chinese (Beidou) navigation constellations were analyzed. The orbit evolution of the rocket bodies was modeled for 200 years, taking into account all relevant perturbations, and the chosen disposal options were evaluated in terms of their long-term consequences for the debris environment. The results obtained, when applicable, were also discussed in the context of the eccentricity instability problem, pointed out in previous studies. In addition, the long-term evolution of the fragments resulting from a Beidou rocket body breakup, and of simulated high area-to-mass ratio objects released in the disposal orbits of the first two Galileo upper stages, was investigated.Eight out of ten Beidou upper stages were found to have an orbital lifetime <25 years and the other two resulted in a dwell time of approximately 6 years below 2000 km. It was also found that the perigee heights of the two upper stages used to deploy the first Galileo test spacecraft will remain more than 169 km above the constellation nominal altitude, never crossing the existing or planned navigation systems. In spite of an inclination resonance possibly leading to the exponential growth of the eccentricity over several decades, the optimal choice of the disposal orbital elements was able to prevent such an outcome, by maintaining the orbit nearly circular. Therefore, the upper stage disposal strategies used so far for Beidou and Galileo have generally been quite successful in averting the long-term interference of such rocket bodies with the navigation constellations, provided that accidental breakups are prevented.     

Assessing the Risk of Orbital Debris Impact

A Space Debris Impact Risk Analysis Tool (SDIRAT) was developed and implemented to assess the orbital debris impact risk on a specified target in Earth orbit, in terms of flux, relative velocity, impact velocity, direction of the incoming particles, debris mass and diameter. Based on a deterministic approach, SDIRAT uses a realistic orbital debris population where each representative particle is identified by its rectangular coordinates (position and velocity) at a reference epoch. Using this information, some geometrical algorithms were developed and implemented to evaluate the contribution of each particle to the incoming flux. The position of the particle with respect to a specified target drives the selection criteria to reject, or select, it as a possible projectile. On the other hand, the relative velocity vector can be used to estimate the impact direction of the incoming flux. SDIRAT was conceived as a general tool for a variety of scenarios, such as low circular and elliptical orbits, up to the geosynchronous ring. This paper presents some examples of possible applications, including the computation of the incoming debris flux on SAX (low Earth orbit), SIRIO (geosynchronous orbit) and the IRIS upper stage (elliptical orbit). Other applications assess the impact risk for the Soviet Radar Ocean Reconnaissance Satellites Cosmos 1900 and Cosmos 1932.     

Drag and energy accommodation coefficients during sunspot maximum

Conditions appropriate to gas-surface interactions on satellite surfaces in orbit have not been successfully duplicated in the laboratory. However, measurements by pressure gauges and mass spectrometers in orbit have revealed enough of the basic physical chemistry that realistic theoretical models of the gas-surface interaction can now be used to calculate physical drag coefficients. The dependence of these drag coefficients on conditions in space can be inferred by comparing the physical drag coefficient of a satellite with a drag coefficient fitted to its observed orbital decay. This study takes advantage of recent data on spheres and attitude stabilized satellites to compare physical drag coefficients with the histories of the orbital decay of several satellites during the recent sunspot maximum. The orbital decay was obtained by fitting, in a least squares sense, the semi-major axis decay inferred from the historical two-line elements acquired by the US Space Surveillance Network. All the principal orbital perturbations were included, namely geopotential harmonics up to the 16th degree and order, third body attraction of the Moon and the Sun, direct solar radiation pressure (with eclipses), and aerodynamic drag, using the Jacchia-Bowman 2006 (JB2006) model to describe the atmospheric density. After adjusting for density model bias, a comparison of the fitted drag coefficient with the physical drag coefficient has yielded values for the energy accommodation coefficient as well as for the physical drag coefficient as a function of altitude during solar maximum conditions. The results are consistent with the altitude and solar cycle variation of atomic oxygen, which is known to be adsorbed on satellite surfaces, affecting both the energy accommodation and angular distribution of the reemitted molecules.     

USA-193 decay predictions using public domain trajectory data and assessment of the post-intercept orbital debris cloud

In early 2008, the need arose to predict the orbital decay of the American spacecraft USA-193, whose characteristics, function and orbit were classified information. With no orbit data and independent Italian tracking capability available, we turned our attention on the orbits determined by a worldwide network of about 20 visual satellite observers. The orbits of USA-193 obtained from their visual observations were therefore used as the sole source of orbit information. Contrary to our expectations, this exercise was extremely successful and we learned a lot in the process. The orbits provided by the visual observers were very accurate for such a low satellite (although the minimum and very stable level of solar activity helped considerably); however, data gaps of a few days were sometimes possible, due to unfavorable pass geometry or weather and light conditions. In any case, the orbital period and the semimajor axis were so accurate that it was possible for us to obtain very good decay fits using special perturbation software, including various atmospheric density models together with all the other relevant perturbing accelerations. We were therefore able to estimate accurate values of the ballistic parameter and the resulting decay and reentry predictions were extremely stable. Amateur optical observations and images of USA-193 had also led to a rough estimation of the shape and sizes of the satellite, revealing that the solar arrays had never been deployed. With this information, and taking into account our estimates of the ballistic parameter, we obtained reasonable and consistent values of the spacecraft mass. Based on previous reentry fragmentation analyses, we were then able to guess the expected USA-193 casualty area, casualty expectancy, debris ground footprint and probability of impact in Italy. Lastly, after the decision by the US Government to destroy the satellite, we independently predicted the interception time windows and the post-event ground tracks. Following the successful spacecraft breakup, we analyzed the evolution of the resulting debris cloud and assessed its (very limited) adverse impact on the circumterrestrial environment.