Achievable debris orbit prediction accuracy using laser ranging data from a single station

Earlier studies have shown that an orbit prediction accuracy of 20 arc sec ground station pointing error for 1–2 day predictions was achievable for low Earth orbit (LEO) debris using two passes of debris laser ranging (DLR) data from a single station, separated by about 24 h. The accuracy was determined by comparing the predicted orbits with subsequent tracking data from the same station. This accuracy statement might be over-optimistic for other parts of orbit far away from the station. This paper presents the achievable orbit prediction accuracy using satellite laser ranging (SLR) data of Starlette and Larets under a similar data scenario as that of DLR. The SLR data is corrupted with random errors of 1 m standard deviation so that its accuracy is similar to that of DLR data. The accurate ILRS Consolidated Prediction Format orbits are used as reference to compute the orbit prediction errors. The study demonstrates that accuracy of 20 arc sec for 1–2 day predictions is achievable.     

Improved kHz-SLR tracking techniques and orbit quality analysis for LEO missions

Satellite gravity field missions such as CHAMP, GRACE and GOCE are designed as low Earth orbiting spacecraft (LEO) with orbit heights of about 250–500 km. The challenging mission objectives require a very precise knowledge of the satellite orbit position in space. For these missions precise orbit information is typically provided by GPS satellite-to-satellite tracking (SST) observations supported by satellite laser ranging (SLR).     

A simple and valid analysis method for orbit anomaly detection

A simple analysis method for orbit anomaly detection, called semi-major axis change method (SACM) was presented by using a relationship between the change of orbit parameters and velocity increments. In this method, the mean value and standard deviation of the semi-major axis change in different time intervals were first calculated according to historical data. Then, these two parameters, the mean value and standard deviation of the semi-major axis change, are chosen as basis variables and combined as an anomalous criterion. For orbit objects with different characteristic, anomalous thresholds were given in different time intervals for identifying the anomalies of the orbital objects. Finally, this method is used for low earth orbit (LEO) satellites and American–Russian breakup debris. By adopting this method, the characteristics of the orbit change were given. The accuracy rate of anomaly analysis for LEO satellites and American–Russian breakup debris can reach to 100%, which demonstrates that the method was rapid and valid.     

Orbital debris–debris collision avoidance

We focus on preventing collisions between debris and debris, for which there is no current, effective mitigation strategy. We investigate the feasibility of using a medium-powered (5 kW) ground-based laser combined with a ground-based telescope to prevent collisions between debris objects in low-Earth orbit (LEO). The scheme utilizes photon pressure alone as a means to perturb the orbit of a debris object. Applied over multiple engagements, this alters the debris orbit sufficiently to reduce the risk of an upcoming conjunction. We employ standard assumptions for atmospheric conditions and the resulting beam propagation. Using case studies designed to represent the properties (e.g. area and mass) of the current debris population, we show that one could significantly reduce the risk of nearly half of all catastrophic collisions involving debris using only one such laser/telescope facility. We speculate on whether this could mitigate the debris fragmentation rate such that it falls below the natural debris re-entry rate due to atmospheric drag, and thus whether continuous long-term operation could entirely mitigate the Kessler syndrome in LEO, without need for relatively expensive active debris removal.     

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.     

Onboard orbit determination using GPS observations based on the unscented Kalman filter

Spaceborne GPS receivers are used for real-time navigation by most low Earth orbit (LEO) satellites. In general, the position and velocity accuracy of GPS navigation solutions without a dynamic filter are 25 m (1σ) and 0.5 m/s (1σ), respectively. However, GPS navigation solutions, which consist of position, velocity, and GPS receiver clock bias, have many abnormal excursions from the normal error range for space operation. These excursions lessen the accuracy of attitude control and onboard time synchronization. In this research, a new onboard orbit determination algorithm designed with the unscented Kalman filter (UKF) was developed to improve the performance. Because the UKF is able to obtain the posterior mean and covariance accurately by using the second-order Taylor series expansion through the sampled sigma points that are propagated by using the true nonlinear system, its performance can be better than that of the extended Kalman filter (EKF), which uses the linearized state transition matrix to predict the covariance. The dynamic models for orbit propagation applied perturbations due to the 40 × 40 geo-potential, the gravity of the Sun and Moon, solar radiation pressure, and atmospheric drag. The 7(8)th-order Runge–Kutta numerical integration was applied for orbit propagation. Two types of observations, navigation solutions and C/A code pseudorange, can be used at the user’s discretion. The performances of the onboard orbit determination were verified using real GPS data of the CHAMP and KOMPSAT-2 satellites. The results of the orbit determination were compared with the precision orbit ephemeris (POE) of the CHAMP and KOMPSAT-2 satellites.     

Laser measurements to space debris from Graz SLR station

In order to test laser ranging possibilities to space debris objects, the Satellite Laser Ranging (SLR) Station Graz installed a frequency doubled Nd:YAG pulse laser with a 1 kHz repetition rate, a pulse width of 10 ns, and a pulse energy of 25 mJ at 532 nm (on loan from German Aerospace Center Stuttgart – DLR). We developed and built low-noise single-photon detection units to enable laser ranging to targets with inaccurate orbit predictions, and adapted our standard SLR software to include a few hundred space debris targets. With this configuration, we successfully tracked – within 13 early-evening sessions of each about 1.5 h – 85 passes of 43 different space debris targets, in distances between 600 km and up to more than 2500 km, with radar cross sections from >15 m² down to <0.3 m², and measured their distances with an average precision of about 0.7 m RMS.     

Removing orbital debris with lasers

Orbital debris in low Earth orbit (LEO) are now sufficiently dense that the use of LEO space is threatened by runaway collision cascading. A problem predicted more than thirty years ago, the threat from debris larger than about 1 cm demands serious attention. A promising proposed solution uses a high power pulsed laser system on the Earth to make plasma jets on the objects, slowing them slightly, and causing them to re-enter and burn up in the atmosphere. In this paper, we reassess this approach in light of recent advances in low-cost, light-weight modular design for large mirrors, calculations of laser-induced orbit changes and in design of repetitive, multi-kilojoules lasers, that build on inertial fusion research. These advances now suggest that laser orbital debris removal (LODR) is the most cost-effective way to mitigate the debris problem. No other solutions have been proposed that address the whole problem of large and small debris. A LODR system will have multiple uses beyond debris removal. International cooperation will be essential for building and operating such a system.     

Modeling of LEO orbital debris populations for ORDEM2008

The NASA Orbital Debris Engineering Model, ORDEM2000, is in the process of being updated to a new version: ORDEM2008. The data-driven ORDEM covers a spectrum of object size from 10 μm to greater than 1 m, and ranging from LEO (low Earth orbit) to GEO (geosynchronous orbit) altitude regimes. ORDEM2008 centimeter-sized populations are statistically derived from Haystack and HAX (the Haystack Auxiliary) radar data, while micron-sized populations are estimated from shuttle impact records. Each of the model populations consists of a large number of orbits with specified orbital elements, the number of objects on each orbit (with corresponding uncertainty), and the size, type, and material assignment for each object. This paper describes the general methodology and procedure commonly used in the statistical inference of the ORDEM2008 LEO debris populations. Major steps in the population derivations include data analysis, reference-population construction, definition of model parameters in terms of reference populations, linking model parameters with data, seeking best estimates for the model parameters, uncertainty analysis, and assessment of the outcomes. To demonstrate the population-derivation process and to validate the Bayesian statistical model applied in the population derivations throughout, this paper uses illustrative examples for the special cases of large-size (>1 m, >32 cm, and >10 cm) populations that are tracked by SSN (the Space Surveillance Network) and also monitored by Haystack and HAX radars operating in a staring mode.     

Geostationary orbit determination using SATRE

A new strategy of precise orbit determination (POD) for GEO (Geostationary Earth Orbit) satellite using SATRE (SAtellite Time and Ranging Equipment) is presented. Two observation modes are proposed and different channels of the same instruments are used to construct different observation modes, one mode receiving time signals from their own station and the other mode receiving time signals from each other for two stations called pairs of combined observations. Using data from such a tracking network in China, the results for both modes are compared. The precise orbit determination for the Sino-1 satellite using the data from 6 June 2005 to 13 June 2005 has been carried out in this work. The RMS (Root-Mean-Square) of observing residuals for 3-day solutions with the former mode is better than 9.1 cm. The RMS of observing residuals for 3-day solutions with the latter mode is better than 4.8 cm, much better than the former mode. Orbital overlapping (3-day orbit solution with 1-day orbit overlap) tests show that the RMS of the orbit difference for the former mode is 0.16 m in the radial direction, 0.53 m in the along-track direction, 0.97 m in the cross-track direction and 1.12 m in the 3-dimension position and the RMS of the orbit difference for the latter mode is 0.36 m in the radial direction, 0.89 m in the along-track direction, 1.18 m in the cross-track direction and 1.52 m in the 3-dimension position, almost the same as the former mode. All the experiments indicate that a meter-level accuracy of orbit determination for geostationary satellite is achievable.     

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.     

Analysis of laser ranges and angular measurements data fusion for space debris orbit determination

In the framework of space debris, the orbit determination process is a fundamental step, both, for researchers and for satellite operators. The accurate knowledge of the orbit of space debris objects is needed to allow space debris characterization studies and to avoid unnecessary collision avoidance maneuvers.The accuracy of the results of an orbit determination process depends on several factors as the number, the accuracy, the kind of processed measurements, their distribution along the orbit, and the object-observer relative geometry. When the observation coverage of the target orbit is not homogeneous, the accuracy of the orbit determination can be improved processing different kind of observables. Recent studies showed that the satellite laser ranging technique can be successfully applied to space debris.In this paper, we will investigate the benefits of using laser ranges and angular measurements for the orbit determination process. We will analyze the influence of the number of used observations, of the covered arc of orbit, of each observable, and of the observation geometry on the estimated parameters. Finally, using data acquired on short observation arcs, we analyze the achievable accuracies for the orbital regimes with the highest space debris density, and to the consequences of the data fusion on catalog maintenance operations. The results shown are obtained using only real data (both angular and laser measurements) provided by sensors of the Swiss Optical Ground Station and Geodynamics Observatory Zimmerwald owned by the Astronomical Institute of the University of Bern (AIUB) and for some studies also using ranges provided from other stations of the International Laser Ranging Service (ILRS).     

The LAGEOS satellites orbit and Yukawa-like interactions

LAGEOS II general relativity pericenter precession has been analysed in terms of the errors produced by the mismodelling of both the gravitational and non-gravitational perturbations acting on the satellite orbit. The accuracy in the pericenter determination may be considered as an upper-bound value for the estimate of the strength α of a possible new-long-range-interaction described by a Yukawa-like potential. In the present work we have focused on the constraints in α that can be obtained with the current best multi-satellites gravity field model EGM96 (α < 2.6 × 10⁻¹⁰) and also with the first promising models from the CHAMP (α < 1.8 × 10⁻¹⁰) and GRACE (α < 1.2 × 10⁻¹⁰) gravimetric missions. These results represent, potentially, an improvement of two or three orders-of-magnitude with respect to the best constraints obtained in the past with Earth–LAGEOS and Lunar–LAGEOS data (|α| < 10⁻⁵–10⁻⁸). The impact of the non-gravitational perturbations mismodelling in the final error budget has been determined together with the improvements obtainable in the constraint of the strength α with the proposed LARES satellite.     

Assessing attitude error of FORMOSAT-3/COSMIC satellites and its impact on orbit determination

An attitude determination and control system (ADCS) is critical to satellite attitude maneuvers and to the coordinate transformation from the inertial frame to the spacecraft frame. This paper shows specific sensors in the ADCS of the satellite mission FORMOSAT-3/COSMIC (F3/C) and the impact of the ADCS quality on orbit accuracy. The selection of main POD antenna depends on the beta angles of the different F3/C satellites (for FM2 and FM4) during the inflight phase. In particular, under the eclipse, alternative attitude sensors are activated to replace the Sun sensors, and such a sensor change leads to anomalous GPS phase residuals and a degraded orbit accuracy. Since the nominal attitude serves as a reference for ADCS, the 3-dimensional attitude-induced errors in reduced dynamic orbits over selected days in 2010 show 9.35, 10.78, 4.97, 5.48, 7.18, and 6.89 cm for FM1–FM6. Besides, the 3-dimensional velocity errors induced by the attitude effect are 0.10, 0.10, 0.07, 0.08, 0.09, and 0.10 for FM1–FM6. We analyze the quality of the observed attitude transformation matrix of F3/C and its impact on kinematic orbit determination. With 249 days of GPS in 2008, the analysis leads to the following averaged 3-dimensional attitude-induced orbit errors: 2.72, 2.62, 2.37, 1.90, 1.70, and 1.99 cm for satellites FM1–FM6. Critical suggestions of geodetic payloads for the follow-on mission of F3/C are presented based on the current result.     

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.     

The fast debris evolution model

The ‘particles-in-a-box’ (PIB) model introduced by Talent [Talent, D.L. Analytic model for orbital debris environmental management. J. Spacecraft Rocket, 29 (4), 508–513, 1992.] removed the need for computer-intensive Monte Carlo simulation to predict the gross characteristics of an evolving debris environment. The PIB model was described using a differential equation that allows the stability of the low Earth orbit (LEO) environment to be tested by a straightforward analysis of the equation’s coefficients. As part of an ongoing research effort to investigate more efficient approaches to evolutionary modelling and to develop a suite of educational tools, a new PIB model has been developed. The model, entitled Fast Debris Evolution (FADE), employs a first-order differential equation to describe the rate at which new objects ?10 cm are added and removed from the environment. Whilst Talent [Talent, D.L. Analytic model for orbital debris environmental management. J. Spacecraft Rocket, 29 (4), 508–513, 1992.] based the collision theory for the PIB approach on collisions between gas particles and adopted specific values for the parameters of the model from a number of references, the form and coefficients of the FADE model equations can be inferred from the outputs of future projections produced by high-fidelity models, such as the DAMAGE model.     

Highly Elliptical Orbits for Arctic observations: Assessment of ionizing radiation

The ionizing radiation environment was analyzed for a variety of potential Highly Elliptical Orbits (HEOs) with orbital periods ranging from 6 h to 24 h suitable to continuously monitor the Arctic region. Several models available from the ESA Space Environment Information System (SPENVIS) online tool were employed, including the new-generation AE9/AP9 model for trapped radiation. Results showed that the Total Ionizing Dose (TID) has a well-pronounced local minimum for the 14-h orbit, which is nearly identical to the overall minimum observed for the longest orbital period (24 h). The thickness of slab aluminum shielding required to keep the annual TID below 10, 5 and 3.33 krad (i.e. 150, 75 and 50 krad for 15 years of mission duration) for a 14-h orbit is 2.1, 2.7 and 3.1 mm respectively. The 16-h orbit requires an additional 0.5 mm of aluminum to achieve the same results, while the 24-h orbit requires less shielding in the order of 0.2–0.3 mm. Comparison between the AE8/AP8 and AE9/AP9 models was conducted for all selected orbits. Results demonstrated that differences ranged from −70% to +170% depending on orbit geometry.     

Precise orbit determination using the batch filter based on particle filtering with genetic resampling approach

In this study, genetic resampling (GRS) approach is utilized for precise orbit determination (POD) using the batch filter based on particle filtering (PF). Two genetic operations, which are arithmetic crossover and residual mutation, are used for GRS of the batch filter based on PF (PF batch filter). For POD, Laser-ranging Precise Orbit Determination System (LPODS) and satellite laser ranging (SLR) observations of the CHAMP satellite are used. Monte Carlo trials for POD are performed by one hundred times. The characteristics of the POD results by PF batch filter with GRS are compared with those of a PF batch filter with minimum residual resampling (MRRS). The post-fit residual, 3D error by external orbit comparison, and POD repeatability are analyzed for orbit quality assessments. The POD results are externally checked by NASA JPL’s orbits using totally different software, measurements, and techniques. For post-fit residuals and 3D errors, both MRRS and GRS give accurate estimation results whose mean root mean square (RMS) values are at a level of 5 cm and 10–13 cm, respectively. The mean radial orbit errors of both methods are at a level of 5 cm. For POD repeatability represented as the standard deviations of post-fit residuals and 3D errors by repetitive PODs, however, GRS yields 25% and 13% more robust estimation results than MRRS for post-fit residual and 3D error, respectively. This study shows that PF batch filter with GRS approach using genetic operations is superior to PF batch filter with MRRS in terms of robustness in POD with SLR observations.     

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