First laser measurements to space debris in Poland

The Borowiec Satellite Laser Ranging station (BORL 7811, Borowiec) being a part of the Space Research Centre of the Polish Academy of Sciences (SRC PAS) went through modernization in 2014–2015. One of the main tasks of the modernization was the installation of a high-energy laser module dedicated to space debris tracking. Surelite III by Continuum is a Nd:YAG pulse laser with 10?Hz repetition rate, a pulse width of 3–5?ns and a pulse energy of 450?mJ for green (532?nm). This new laser unit was integrated with the SLR system at Borowiec performing standard satellite tracking. In 2016 BORL 7811 participated actively to the observational campaigns related to the space debris targets from LEO region managed by the Space Debris Study Group (SDSG) of the International Laser Ranging Service (ILRS).Currently, Borowiec station regularly tracks 36 space debris from the LEO regime, including typical rocket bodies (Russian/Chinese) and cooperative targets like the inactive TOPEX/Poseidon, ENVISAT, OICETS and others. In this paper the first results of space debris laser measurements obtained by the Borowiec station in period August 2016 – January 2017 are presented. The results gained by the SRC PAS Borowiec station confirm the rotation of the defunct TOPEX/Poseidon satellite which spins with a period of approximately 10?s. The novelty of this work is the presentation of the sample results of the Chinese CZ-2C R/B target (NORAD catalogue number 31114) which is equipped (probably) with retroreflectors. Laser measurements to space debris is a very desirable topic for the next years, especially in the context of the Space Surveillance and Tracking (SST) activity. Some targets are very easy to track like defunct ENVISAT or TOPEX/Poseidon. On the other hand, there is a big population of different LEO targets with different orbital and physical parameters, which are challenging for laser ranging like small irregular debris and rocket boosters.     

A robotic telescope network for space debris identification and tracking

The two TAROT (Télescopes à Action Rapide pour les Objets Transitoires; Rapid Action Telescopes for Transient Objects) installations are fully robotic optical observatories with optimized observation scheduling, data processing and archiving. Zadko is a 1 m telescope in Western Australia. The fully robotisation of the Zadko telescope has just been completed; it is now included in the TAROT network. In this paper we provide an overview of this international network of robotic optical telescopes. We discuss the advantages of using the network to participate in a satellite and space debris tracking program. This network will access almost all geostationary belt objects, and provide the first real-time satellite positioning capability. The inclusion of the 1 m Zadko telescope into the network significantly extends the efficiency and sensitivity of the existing two telescope configuration.     

Prospects in the orbital and rotational dynamics of the Moon with the advent of sub-centimeter lunar laser ranging

Lunar laser ranging (LLR) measurements are crucial for advanced exploration of the laws of fundamental gravitational physics and geophysics as well as for future human and robotic missions to the Moon. The corner-cube reflectors (CCR) currently on the Moon require no power and still work perfectly since their installation during the project Apollo era. Current LLR technology allows us to measure distances to the Moon with a precision approaching 1 mm. As NASA pursues the vision of taking humans back to the Moon, new, more precise laser ranging applications will be demanded, including continuous tracking from more sites on Earth, placing new CCR arrays on the Moon, and possibly installing other devices such as transponders, etc. for multiple scientific and technical purposes. Since this effort involves humans in space, then in all situations the accuracy, fidelity, and robustness of the measurements, their adequate interpretation, and any products based on them, are of utmost importance. Successful achievement of this goal strongly demands further significant improvement of the theoretical model of the orbital and rotational dynamics of the Earth–Moon system. This model should inevitably be based on the theory of general relativity, fully incorporate the relevant geophysical processes, lunar librations, tides, and should rely upon the most recent standards and recommendations of the IAU for data analysis. This paper discusses methods and problems in developing such a mathematical model. The model will take into account all the classical and relativistic effects in the orbital and rotational motion of the Moon and Earth at the sub-centimeter level. The model is supposed to be implemented as a part of the computer code underlying NASA Goddard’s orbital analysis and geophysical parameter estimation package GEODYN and the ephemeris package PMOE 2003 of the Purple Mountain Observatory. The new model will allow us to navigate a spacecraft precisely to a location on the Moon. It will also greatly improve our understanding of the structure of the lunar interior and the nature of the physical interaction at the core–mantle interface layer. The new theory and upcoming millimeter LLR will give us the means to perform one of the most precise fundamental tests of general relativity in the solar system.     

Design of tracking mount and controller for mobile satellite laser ranging system

In this study, we have proposed and implemented a design for the tracking mount and controller of the ARGO-M (Accurate Ranging system for Geodetic Observation - Mobile) which is a mobile satellite laser ranging (SLR) system developed by the Korea Astronomy and Space Science Institute (KASI) and Korea Institute of Machinery and Materials (KIMM). The tracking mount comprises a few core components such as bearings, driving motors and encoders. These components were selected as per the technical specifications for the tracking mount of the ARGO-M. A three-dimensional model of the tracking mount was designed. The frequency analysis of the model predicted that the first natural frequency of the designed tracking mount was high enough. The tracking controller is simulated using MATLAB/xPC Target to achieve the required pointing and tracking accuracy. In order to evaluate the system repeatability and tracking accuracy of the tracking mount, a prototype of the ARGO-M was fabricated, and repeatability tests were carried out using a laser interferometer. Tracking tests were conducted using the trajectories of low earth orbit (LEO) and high earth orbit (HEO) satellites. Based on the test results, it was confirmed that the prototype of the tracking mount and controller of the ARGO-M could achieve the required repeatability along with a tracking accuracy of less than 1 arcsec.     

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.     

Uncertainty propagation for statistical impact prediction of space debris

Predictions of the impact time and location of space debris in a decaying trajectory are highly influenced by uncertainties. The traditional Monte Carlo (MC) method can be used to perform accurate statistical impact predictions, but requires a large computational effort. A method is investigated that directly propagates a Probability Density Function (PDF) in time, which has the potential to obtain more accurate results with less computational effort. The decaying trajectory of Delta-K rocket stages was used to test the methods using a six degrees-of-freedom state model. The PDF of the state of the body was propagated in time to obtain impact-time distributions. This Direct PDF Propagation (DPP) method results in a multi-dimensional scattered dataset of the PDF of the state, which is highly challenging to process. No accurate results could be obtained, because of the structure of the DPP data and the high dimensionality. Therefore, the DPP method is less suitable for practical uncontrolled entry problems and the traditional MC method remains superior. Additionally, the MC method was used with two improved uncertainty models to obtain impact-time distributions, which were validated using observations of true impacts. For one of the two uncertainty models, statistically more valid impact-time distributions were obtained than in previous research.     

An adaptive threshold method for improving astrometry of space debris CCD images

Optical survey is a main technique for observing space debris, and precisely measuring the positions of space debris is of great importance. Due to several factors, e.g. the angle object normal to the observer, the shape as well as the attitude of the object, the variations of observed characteristics for low earth orbital space debris are distinct. When we look at optical CCD images of observed objects, the size and brightness are varying, hence it’s difficult to decide the threshold during centroid measurement and precise astrometry. Traditionally the threshold is given empirically and constantly in data reduction, and obviously it’s not suitable for data reduction of space debris. Here we offer a solution to provide the threshold. Our method assumes that the PSF (point spread function) is Gaussian and estimates the signal flux by a directly two-dimensional Gaussian fit, then a cubic spline interpolation is performed to divide each initial pixel into several sub-pixels, at last the threshold is determined by the estimation of signal flux and the sub-pixels above threshold are separated to estimate the centroid. A trail observation of the fast spinning satellite Ajisai is made and the CCD frames are obtained to test our algorithm. The calibration precision of various threshold is obtained through the comparison between the observed equatorial position and the reference one, the latter are obtained from the precise ephemeris of the satellite. The results indicate that our method reduces the total errors of measurements, it works effectively in improving the centering precision of space debris images.     

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.     

Using parallel computing for the display and simulation of the space debris environment

Parallelism is becoming the leading paradigm in today’s computer architectures. In order to take full advantage of this development, new algorithms have to be specifically designed for parallel execution while many old ones have to be upgraded accordingly. One field in which parallel computing has been firmly established for many years is computer graphics. Calculating and displaying three-dimensional computer generated imagery in real time requires complex numerical operations to be performed at high speed on a large number of objects. Since most of these objects can be processed independently, parallel computing is applicable in this field. Modern graphics processing units (GPUs) have become capable of performing millions of matrix and vector operations per second on multiple objects simultaneously.     

A new approach to the telescope operation method for satellite tracking using a time synchronization technique

For the development of a telescope that is capable of precisely tracking satellites and high-speed operation such as satellite laser ranging, a special method of telescope operation is required. This study aims to propose a new telescope operation method and system configuration for the independent development of a mount and an operation system which includes the host computer. Considering that the tracking of a satellite is performed in real time, communication and synchronization between the two independent subsystems are important. Therefore, this study applied the concept of time synchronization, which is used in various fields of industry, to the communication between the command computer and the mount. In this case, communication delays do not need to be considered in general, and it is possible to cope with data loss. Above all, when the mount is replaced in the future, only the general communication interface needs to be modified, and thus, it is not limited by replacement in terms of the overall system management. The performance of the telescope operation method developed in this study was verified by applying the method to the first mobile SLR system in Korea. This study is significant in that it proposed a new operation method and system configuration, to which the concept of time synchronization was applied, for the observation system that requires an optical telescope.     

A new algorithm for optical observations of space debris with the TAROT telescopes

Since 2004, we observe satellites in the geostationary orbit with a network of robotic ground based fully automated telescopes called TAROT. One of them is located in France and the second at ESO, La Silla, Chile. The system processes the data in real time. Its wide field of view is useful for the discovery, the systematic survey and for the tracking of both catalogued and un-catalogued objects. We present a new source extraction algorithm based on morphological mathematic, which has been tested and is currently under implementation in the standard pipeline. Using this method, the observation strategy will correlate the measurements of the same object on successive images and give better detection rate and false alarm rate than the previous one. The overall efficiency and quality of the survey of the geostationary orbit has drastically improved and we can now detect satellites and debris in different orbits like Geostationary Transfer Orbit (GTO). Results obtained in real conditions with TAROT are presented.     

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.     

Accurate orbit predictions for debris orbit manoeuvre using ground-based lasers

Orbit manoeuvre of low Earth orbiting (LEO) debris using ground-based lasers has been proposed as a cost-effective means to avoid debris collisions. This requires the orbit of the debris object to be determined and predicted accurately so that the laser beam can be locked on the debris without the loss of valuable laser operation time. This paper presents the method and results of a short-term accurate LEO (<900 km in altitude) debris orbit prediction study using sparse laser ranging data collected by the EOS Space Debris Tracking System (SDTS). A main development is the estimation of the ballistic coefficients of the LEO objects from their archived long-term two line elements (TLE). When an object is laser tracked for two passes over about 24 h, orbit prediction (OP) accuracy of 10–20 arc seconds for the next 24–48 h can be achieved – the accuracy required for laser debris manoeuvre. The improvements in debris OP accuracy are significant in other applications such as debris conjunction analyses and the realisation of daytime debris laser tracking.     

Impact effects from small size meteoroids and space debris

When the impact risk from meteoroids and orbital debris is assessed the main concern is usually structural damage. With their high impact velocities of typically 10–20 km/s millimeter or centimeter sized objects can puncture pressure vessels and other walls or lead to destruction of complete subsystems or even whole spacecraft. Fortunately chances of collisions with such larger objects are small (at least at present). However, particles in the size range 1–100 μm are far more abundant than larger objects and every orbiting spacecraft will encounter them with certainty. Every solar cell (8 cm² area) of the Hubble Space Telescope encountered on average 12 impacts during its 8.25 years of space exposure. Most were from micron sized particles.     

半导体激光器干涉测距进展

全面评述了近年来半导体激光器应用于干涉测距的进展,对其中的代表性文献进行了详细的分析,指出了该技术今后的主要发展方向和关键问题。     

Fuzzy multi-criteria decision-making approach to prioritization of space debris for removal

As private companies and government space agencies begin to seriously consider the task of active space debris removal, it is becoming increasingly more important to determine the highest priority objects to deorbit. This work sets forth an approach for prioritization of space debris through the utilization of Multi-Criteria Decision-Making methodologies and fuzzy logic, as well as both quantitative and qualitative criteria. The proposed debris prioritization approach considers various criteria including the orbit, size, mass, pairwise and total collision probabilities, and decay timeframe of each debris object. The means of assigning attributes to each assessment criterion is discussed in detail. To determine the weighting scheme for the criteria, a questionnaire was prepared and shared with experts in the field of space situational awareness. The work examines over two thousand critical debris objects selected from the existing debris catalog with respect to these criteria. The quantified attributes for each debris object are then aggregated through the fuzzy versions of the Analytic Hierarchy Process and the Technique for Order Preference by Similarity to Ideal Solution. The results of the analysis identify high-priority debris objects for removal from Earth-bound orbits.     

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.     

Semi-analytical theory of mean orbital motion for geosynchronous space debris under gravitational influence

This paper provides a hamiltonian formulation of the equations of motion of an artificial satellite or space debris orbiting the geostationary ring. This theory of order 1 has been 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. The theory accounts for the influence of the Earth gravity field expanded in spherical harmonics, paying a particular attention to the resonance occurring for geosynchronous objects. The luni-solar perturbations are also taken into account. We present the resonant motion and its main characteristics: equilibria, stability, fundamental frequencies and width of the resonant area by comparison with a basic analytical model. Finally, we show some results concerning the long term dynamics of a typical space debris under the influence of the gravitational field of the Earth and the luni-solar interactions.     

Comparison of self-healing ionomer to aluminium-alloy bumpers for protecting spacecraft equipment from space debris impacts

This paper discusses the impact behavior of a self-healing ionomeric polymer and compares its protection capability against space debris impacts to that of simple aluminium-alloy bumpers. To this end, 14 impact experiments on both ionomer and Al-7075-T6 thin plates with similar surface density were made with 1.5 mm aluminium spheres at velocity between 1 and 4 km/s.