Feasibility of performing space surveillance tasks with a proposed space-based optical architecture

Under ESA contract an industrial consortium including Aboa Space Research Oy (ASRO), the Astronomical Institute of the University of Bern (AIUB), and the Dutch National Aerospace Laboratory (NLR), proposed the observation concept, developed a suitable sensor architecture, and assessed the performance of a space-based optical (SBO) telescope in 2005. The goal of the SBO study was to analyse how the existing knowledge gap in the space debris population in the millimetre and centimetre regime may be closed by means of a passive optical instrument. The SBO instrument was requested to provide statistical information on the space debris population in terms of number of objects and size distribution. The SBO instrument was considered to be a cost-efficient with 20 cm aperture and 6° field-of-view and having flexible integration requirements. It should be possible to integrate the SBO instrument easily as a secondary payload on satellites launched into low-Earth orbits (LEO), or into geostationary orbit (GEO). Thus the selected mission concept only allowed for fix-mounted telescopes, and the pointing direction could be requested freely. Since 2007 ESA focuses space surveillance and tracking activities in the Space Situational Awareness (SSA) preparatory program. Ground-based radars and optical telescopes are studied for the build-up and maintenance of a catalogue of objects. In this paper we analyse how the proposed SBO architecture could contribute to the space surveillance tasks survey and tracking. We assume that the SBO instrumentation is placed into a circular sun-synchronous orbit at 800 km altitude. We discuss the observation conditions of objects at higher altitude, and select an orbit close to the terminator plane. A pointing of the sensor orthogonal to the orbital plane with optimal elevation slightly in positive direction (0° and +5°) is found optimal for accessing the entire GEO regime within one day, implying a very good coverage of controlled objects in GEO, too. Simulations using ESA’s Program for Radar and Optical Observation Forecasting (PROOF) in the version 2005 and a GEO reference population extracted from DISCOS revealed that the proposed pointing scenario provides low phase angles together with low angular velocities of the objects crossing the field-of-view. Radiometric simulations show that the optimal exposure time is 1–2 s, and that spherical objects in GEO with a diameter of below 1 m can be detected. The GEO population can be covered under proper illumination nearly completely, but seasonal drops of the coverage are possible. Subsequent observations of objects are on average at least every 1.5 days, not exceeding 3 days at maximum. A single observation arc spans 3° to 5° on average. Using a simulation environment that connects PROOF to AIUB’s program system CelMech we verify the consistency of the initial orbit determination for five selected test objects on subsequent days as a function of realistic astrometric noise levels. The initial orbit determination is possible. We define requirements for a correlator process essential for catalogue build-up and maintenance. Each single observation should provide an astrometric accuracy of at least 1”–1.5” so that the initially determined orbits are consistent within a few hundred kilometres for the semi-major axis, 0.01 for the eccentricity, and 0.1° for the inclination.     

Mathematical description of the NaK model for MASTER-2005

The release of NaK droplets has been modeled for the new version of the European Meteoroid and Space Debris Terrestrial Environment Reference model MASTER-2005. Previously published versions of the model have been revised. The parameters of the model are introduced and discussed. NaK droplets consist of eutectic sodium–potassium alloy and have been released during RORSAT reactor core ejections. They contributed to the space debris environment in the centimeter and millimeter size regime. Sixteen nuclear powered RORSATs launched between 1980 and 1988 activated a reactor core ejection system in Sufficiently High Orbits (SHO), mostly between 900 and 950 km altitude. The core ejection caused an opening of the primary coolant circuit. The liquid coolant has been released into space during these core ejections. The outflow is considered as a discrete event for each of the sixteen core ejections in total. The NaK coolant has been forming droplets up to a diameter of 5.5 cm. NaK releases are restricted to a very narrow region near 65° inclination. This paper gives the parameters of the NaK release model as it is implemented in MASTER-2005. The quantitative values of all model parameters including characteristic diameter and uniformity parameter are presented. The ratio of the characteristic droplet size to the orifice diameter is discussed. It is estimated that altogether 128 kg of NaK-78 (8 kg per RORSAT) was released on orbit. Simulation runs show that there are still 45,000 droplets with a total mass of 97 kg in orbit at the reference epoch 1 May 2005, whereas the smallest droplet has a diameter of 5 mm. Results of orbit propagation simulation runs are presented in terms of spatial density.     

ESA’s process for the identification and assessment of high-risk conjunction events

ESA’s Space Debris Office provides an operational service for the assessment of collision risks of ESA satellites. Currently, the ENVISAT and ERS-2 missions in low Earth orbits are covered by this service. If an upcoming high-risk conjunction event is predicted based on analysis of Two-Line Element (TLE) data from the US Space Surveillance Network, then independent tracking data of the potential high-risk conjunction object are acquired to improve the knowledge of its orbit. This improved knowledge and the associated small error covariances derived from the orbit determination process scale down the position error ellipsoid at the conjunction epoch. Hence, for the same miss-distance, in most cases an avoidance manoeuvre can be suppressed with an acceptable residual risk.     

Orbit prediction and shuttle rendezvous strategy for eureca

The European Retrievable Carrier (EURECA) is a platform to be launched, deployed and retrieved in low Earth orbit by the Space Shuttle.A newly developed analytical orbit prediction method is described which meets the severe requirements for EURECA's orbit propagation. It is based on an averaging procedure including the Earth's zonal harmonics J₂, J₃ and J₄ and a refined treatment of the air drag perturbation where EURECA's large solar panels are taken into account. Some orbit prediction results are included.In order to offer more flexibility for the Shuttle retrieval of EURECA, it is proposed to execute a part of the rendezvous manoeuvres by EURECA. A corresponding strategy is described.     

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.     

NaK release model for MASTER-2009

Sodium–potassium droplets from the primary coolant loop of Russian orbital reactors have been released into space. These droplets are called NaK droplets. Sixteen nuclear powered satellites of the type RORSAT launched between 1980 and 1988 activated a reactor core ejection system, mostly between 900 and 950 km altitude. The core ejection causes an opening of the primary coolant loop. The liquid coolant consists of eutectic sodium–potassium alloy and has been released into space during these core ejections. The NaK coolant has been forming droplets up to a diameter of 5.5 cm. NaK droplets have been modeled before in ESA's MASTER Debris and Meteoroid Environment Model. The approach is currently revised for the MASTER-2009 upgrade. A mathematical improvement is introduced by substituting the current size distribution function by the modified Rosin–Rammler equation. A bimodal size distribution is derived which is based on the modified mass based Rosin–Rammler equation. The equation is modified by truncating the size range and normalizing over the finite range between the size limits of the smallest and the biggest droplet. The parameters of the model are introduced and discussed. For the validation of the NaK release model, sixteen release events are simulated. The resulting size distribution is compared with radar measurement data. The size distribution model fits well with revised published measurement data of radar observations. Results of orbit propagation simulation runs are presented in terms of spatial density.     

DISCOS - ESA''s database and information system characterising objects in space

This paper presents an overview of the main features of ESA's future space debris database DISCOS (Database and Information System Characterising Objects in Space). The DISCOS system has been developed around an ORACLE relational database management software by the University of Kent (UK) under an ESA contract. The DISCOS catalogue will be installed at ESOC, the European Space Operations Centre, and serve as a common ESA information system for the space debris environment.     

Design and development of correlation techniques to maintain a space surveillance system catalogue

A growing interest exists in a future, autonomous European Space Surveillance System (ESSS). Currently, most of the knowledge about Earth-orbiting space objects is based on information provided by the USASPACECOM. This paper presents the required initial orbit determination (IOD) and correlation techniques to process optical measurements. Former studies were focused on the handling of radar measurements, which are summarised with the aim of describing a global procedure for processing hybrid measurement types (combination of radar and optic data for catalogue maintenance). The introduction of manoeuvres are presented due to their importance in the space object catalogue maintenance.The detection of uncatalogued objects and the successful correlation of already catalogued objects involve two different tasks for telescopes: survey and tasking. Assumptions for both strategies are developed on the basis of the previous work developed at the University of Berne (see [T. Flohrer, T. Schildknecht, R. Musci, E. Stöveken, Performance estimation for GEO space surveillance, Advances in Space Research 35 (2005). [1]; T. Flohrer, T. Schildknecht, R. Musci, Proposed strategies for optical observations in a future European Space Surveillance Network, presented in the 36th COSPAR Scientific Assembly (2006). [2]; R. Musci, T. Schildknecht, M. Ploner, Orbit improvement for GEO objects using follow-up observations, Advances in Space Research 34 (2004). [3]; R. Musci, T. Schildknecht, M. Ploner, G. Beutler, Orbit improvement for GTO objects using follow-up observations, Advances in Space Research 35 (2005). [4]; R. Musci, T. Schildknecht, T. Flohrer, G. Beutler, Concept for a catalogue of space debris in GEO, Proceedings of the Fourth European Conference on Space Debris, (ESA SP-587, 2005). [5]]). When a new object appears in the field of view, initial orbit determination must be performed. When only one telescope per site is available, the initial measurements are separated by only a few seconds. Therefore, the initial orbit determination is quite inaccurate due to bad mathematical conditioning of the problem. In order to improve the accuracy of the initial orbit determination, several follow-up observations of the object are required. This implies that the telescope needs to track the detected objects over a long period, and therefore the time available for surveying is reduced. By processing the additional follow-up measurements, separated now by a few hours, the initial orbit determination gives more accurate results and the object can be recovered after an orbital revolution. When several telescopes per site are available, the optical strategies may be modified. The survey tasks can be distributed between the available telescopes. In this way the number of images corresponding to each object increases and to track the detected object over long periods is not always needed. Numerical results will be shown in order to evaluate the accuracy and features of the different telescope strategies. A key point for performing efficiently the cataloguing process is the calculation of the estimated state vector covariance matrix. The covariance matrix analysis allows an adaptive tasking-survey telescope scheduling. Moreover, the implementation of a proper batch orbit determination process by means of a square root information filter (SRIF) requires a realistic initial covariance matrix.Hybrid measurements are available from objects that can be observed through both radar and optical sensors (e.g. GTO objects). The batch orbit determination and correlation process of hybrid measurements is also based on SRIF using an extended measurement model. Both the initial orbit determination methods using radar and optical measurements have to be sufficiently accurate to initialise SRIF correctly. In order to avoid filter divergence, the estimated covariance must be correctly updated after processing both kinds of measurements. The implemented algorithms are explained and their performance is shown through realistic simulations.Techniques to detect and characterise object manoeuvres during the cataloguing process have been developed and implemented. Four main groups of manoeuvre objects have been established by means of their observed permitted orbital ranges (GEO, LEO, MEO–GPS, Molniya). The study is based on the historical TLEs files. When an object with an uncatalogued orbit appears, a comparison between the new orbit and the orbits contained in the permitted ranges of one of the manoeuvre groups is performed. If the required Δa and/or Δi to convert the lost orbit into the detected orbit seems to be feasible, a manoeuvre will be identified and the orbit will be updated in the catalogue. Otherwise, it will be decided that a new object was found. For this purpose, a procedure to estimate the manoeuvres and reset orbits have been developed.     

Use of two-line element data for thermosphere neutral density model calibration

Traditional empirical thermospheric density models are widely used in orbit determination and prediction of low-Earth satellites. Unfortunately, these models often exhibit large density errors of up to around 30% RMS. Density errors translate into orbit errors, adversely affecting applications such as re-entry operations, manoeuvre planning, collision avoidance and precise orbit determination for geodetic missions. The extensive database of two-line element (TLE) orbit data contains a wealth of information on satellite drag, at a sufficiently high spatial and temporal resolution to allow a calibration of existing neutral density models with a latency of one to two days. In our calibration software, new TLE data for selected objects is converted to satellite drag data on a daily basis. The resulting drag data is then used in a daily adjustment of density model calibration parameters, which modify the output of an existing empirical density model with the aim of increasing its accuracy. Two different calibration schemes have been tested using TLE data for about 50 objects during the year 2000. The schemes involve either height-dependent scale factors to the density or corrections to CIRA-72 model temperatures, which affect the density output based on a physical model. Both schemes have been applied with different spherical harmonic expansions of the parameters in latitude and local solar time. Five TLE objects, varying in perigee altitude between 280 and 530 km, were deliberately not used during calibration, in order to provide independent validation. Even with a single daily parameter, the RMS density model error along their tracks can already be reduced from the 30% to the 15% level. Adding additional parameters results in RMS errors lower than 12%.     

Additional historical solid rocket motor burns

The use of orbital solid rocket motors (SRM) is responsible for the release of a high number of slag and Al₂O₃ dust particles which contribute to the space debris environment. This contribution has been modeled for the ESA space debris model MASTER (Meteoroid and Space Debris Terrestrial Environment Reference). The current model version, MASTER-2005, is based on the simulation of 1076 orbital SRM firings which mainly contributed to the long-term debris environment. SRM firings on very low earth orbits which produce only short living particles are not considered. A comparison of the modeled flux with impact data from returned surfaces shows that the shape and quantity of the modeled 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 analyzed Long Duration Exposure Facility (LDEF) surfaces. This indicates that some past SRM firings are not included in the current event database. Thus it is necessary to investigate, if additional historical SRM burns, like the retro-burn of low orbiting re-entry capsules, may be responsible for these dust impacts. The most suitable candidates for these firings are the large number of SRM retro-burns of return capsules. This paper focuses on the SRM retro-burns of Russian photoreconnaissance satellites, which were used in high numbers during the time of the LDEF mission. It is discussed which types of satellites and motors may have been responsible for this historical contribution. Altogether, 870 additional SRM retro-burns have been identified. An important task is the identification of such missions to complete the current event data base. Different types of motors have been used to de-orbit both large satellites and small film return capsules. The results of simulation runs are presented.