Optical Light Curve Observations to Determine Attitude States of Space Debris

The currently proposed space debris remediation measures include the active removal of large objects and “just in time” collision avoidance by deviating the objects using, e.g., ground-based lasers. Both techniques require precise knowledge of the attitude state and state changes of the target objects. In the former case, to devise methods to grapple the target by a tug spacecraft, in the latter, to precisely propagate the orbits of potential collision partners as disturbing forces like air drag and solar radiation pressure depend on the attitude of the objects. Non-resolving optical observations of the magnitude variations, so-called light curves, are a promising technique to determine rotation or tumbling rates and the orientations of the actual rotation axis of objects, as well as their temporal changes. The 1-meter telescope ZIMLAT of the Astronomical Institute of the University of Bern has been used to collect light curves of MEO and GEO objects for a considerable period of time. Recently, light curves of Low Earth Orbit (LEO) targets were acquired as well. We present different observation methods, including active tracking using a CCD subframe readout technique, and the use of a high-speed scientific CMOS camera. Technical challenges when tracking objects with poor orbit redictions, as well as different data reduction methods are addressed. Results from a survey of abandoned rocket upper stages in LEO, examples of abandoned payloads and observations of high area-to-mass ratio debris will be resented. Eventually, first results of the analysis of these light curves are provided.

Extraction of spin periods of space debris from optical light curves

The population of space debris increased drastically during the last years. These objects have become a great threat for active satellites. Because the relative velocities between space debris and satellites are high, space debris objects may destroy active satellites through collisions. Furthermore, collisions involving massive objects produce large number of fragments leading to significant growth of the space debris population. The long term evolution of the debris population is essentially driven by so-called catastrophic collisions. An effective remediation measure in order to stabilize the population in Low Earth Orbit (LEO) is therefore the removal of large, massive space debris. To remove these objects, not only precise orbits, but also more detailed information about their attitude states will be required. One important property of an object targeted for removal is its spin period, spin axis orientation and their change over time. Rotating objects will produce periodic brightness variations with frequencies which are related to the spin periods. Such a brightness variation over time is called a light curve. Collecting, but also processing light curves is challenging due to several reasons. Light curves may be undersampled, low frequency components due to phase angle and atmospheric extinction changes may be present, and beat frequencies may occur when the rotation period is close to a multiple of the sampling period. Depending on the method which is used to extract the frequencies, also method-specific properties have to be taken into account. The astronomical Institute of the University of Bern (AIUB) light curve database will be introduced, which contains more than 1,300 light curves acquired over more than seven years. We will discuss properties and reliability of different time series analysis methods tested and currently used by AIUB for the light curve processing. Extracted frequencies and reconstructed phases for some interesting targets, e.g. GLONASS satellites, for which also SLR data were available for the period confirmation, will be presented. Finally we will present the reconstructed phase and its evolution over time of a High-Area-to-Mass-Ratio (HAMR) object, which AIUB observed for several years.

A study of the time complexity of an Elitist Genetic Algorithm used for associating optical observations of MEO and GEO Space Debris

Currently several thousands of objects are being tracked in the MEO and GEO regions through optical means. The problem faced in this framework is that of Multiple Target Tracking (MTT). In this context both the correct associations among the observations, and the orbits of the objects have to be determined. The complexity of the MTT problem is defined by its dimension S. Where S stands for the number of ’fences’ used in the problem, each fence consists of a set of observations that all originate from dierent targets. For a dimension of S ˃ the MTT problem becomes NP-hard. As of now no algorithm exists that can solve an NP-hard problem in an optimal manner within a reasonable (polynomial) computation time. However, there are algorithms that can approximate the solution with a realistic computational e ort. To this end an Elitist Genetic Algorithm is implemented to approximately solve the S ˃ MTT problem in an e cient manner. Its complexity is studied and it is found that an approximate solution can be obtained in a polynomial time. With the advent of improved sensors and a heightened interest in the problem of space debris, it is expected that the number of tracked objects will grow by an order of magnitude in the near future. This research aims to provide a method that can treat the correlation and orbit determination problems simultaneously, and is able to e ciently process large data sets with minimal manual intervention.

Space debris attitude simulation - ιOTA (In-Orbit Tumbling Analysis)

Today, there is little knowledge on the attitude state of decommissioned intact objects in Earth orbit. Observational means have advanced in the past years, but are still limited with respect to an accurate estimate of motion vector orientations and magnitude. Especially for the preparation of Active Debris Removal (ADR) missions as planned by ESA’s Clean Space initiative or contingency scenarios for ESA spacecraft like ENVISAT, such knowledge is needed. ESA's “Debris Attitude Motion Measurements and Modelling” project (ESA Contract No. 40000112447), led by the Astronomical Institute of the University of Bern (AIUB), addresses this problem. The goal of the project is to achieve a good understanding of the attitude evolution and the considerable internal and external effects which occur. To characterize the attitude state of selected targets in LEO and GTO, multiple observation methods are combined. Optical observations are carried out by AIUB, Satellite Laser Ranging (SLR) is performed by the Space Research Institute of the Austrian Academy of Sciences (IWF) and radar measurements and signal level determination are provided by the Fraunhofer Institute for High Frequency Physics and Radar Techniques (FHR). The In-Orbit Tumbling Analysis tool (ιOTA) is a prototype software, currently in development by Hyperschall Technologie Göttingen GmbH (HTG) within the framework of the project. ιOTA will be a highly modular software tool to perform short-(days), medium-(months) and long-term (years) propagation of the orbit and attitude motion (six degrees-of-freedom) of spacecraft in Earth orbit. The simulation takes into account all relevant acting forces and torques, including aerodynamic drag, solar radiation pressure, gravitational influences of Earth, Sun and Moon, eddy current damping, impulse and momentum transfer from space debris or micro meteoroid impact, as well as the optional definition of particular spacecraft specific influences like tank sloshing, reaction wheel behaviour, magnetic torquer activity and thruster firing. The purpose of ιOTA is to provide high accuracy short-term simulations to support observers and potential ADR missions, as well as medium-and long-term simulations to study the significance of the particular internal and external influences on the attitude, especially damping factors and momentum transfer. The simulation will also enable the investigation of the altitude dependency of the particular external influences. ιOTA's post-processing modules will generate synthetic measurements for observers and for software validation. The validation of the software will be done by cross-calibration with observations and measurements acquired by the project partners.

Streak detection algorithm for space debris detection on optical images.

Any image processing object detection algorithm somehow tries to integrate the object light (Recognition Step) and applies statistical criteria to distinguish objects of interest from other objects or from pure background (Decision Step). There are various possibilities how these two basic steps can be realized, as can be seen in the different proposed detection methods in the literature. An ideal detection algorithm should provide high recognition sensitiv ity with high decision accuracy and require a reasonable computation effort . In reality, a gain in sensitivity is usually only possible with a loss in decision accuracy and with a higher computational effort. So, automatic detection of faint streaks is still a challenge. This paper presents a detection algorithm using spatial filters simulating the geometrical form of possible streaks on a CCD image. This is realized by image convolution. The goal of this method is to generate a more or less perfect match between a streak and a filter by varying the length and orientation of the filters. The convolution answers are accepted or rejected according to an overall threshold given by the ackground statistics. This approach yields as a first result a huge amount of accepted answers due to filters partially covering streaks or remaining stars. To avoid this, a set of additional acceptance criteria has been included in the detection method. All criteria parameters are justified by background and streak statistics and they affect the detection sensitivity only marginally. Tests on images containing simulated streaks and on real images containing satellite streaks show a very promising sensitivity, reliability and running speed for this detection method. Since all method parameters are based on statistics, the true alarm, as well as the false alarm probability, are well controllable. Moreover, the proposed method does not pose any extraordinary demands on the computer hardware and on the image acquisition process.

Photometric monitoring of non-resolved space debris and databases of optical light curves.

The population of space debris increased drastically during the last years. Collisions involving massive objects may produce large number of fragments leading to significantly growth of the space debris population. An effective remediation measure in order to stabilize the population in LEO, is therefore the removal of large, massive space debris. To remove these objects, not only precise orbits, but also more detailed information about their attitude states will be required. One important property of an object targeted for removal is its spin period and spin axis orientation. If we observe a rotating object, the observer sees different surface areas of the object which leads to changes in the measured intensity. Rotating objects will produce periodic brightness vari ations with frequencies which are related to the spin periods. Photometric monitoring is the real tool for remote diagnostics of the satellite rotation around its center of mass. This information is also useful, for example, in case of contingency. Moreover, it is also important to take into account the orientation of non-spherical body (e.g. space debris) in the numerical integration of its motion when a close approach with the another spacecr aft is predicted. We introduce the two databases of light curves: the AIUB data base, which contains about a thousand light curves of LEO, MEO and high-altitude debris objects (including a few functional objects) obtained over more than seven years, and the data base of the Astronomical Observatory of Odessa University (Ukraine), which contains the results of more than 10 years of photometric monitoring of functioning satellites and large space debris objects in low Earth orbit. AIUB used its 1m ZIMLAT telescope for all light curves. For tracking low-orbit satellites, the Astronomical Observatory of Odessa used the KT-50 telescope, which has an alt-azimuth mount and allows tracking objects moving at a high angular velocity. The diameter of the KT-50 main mirror is 0.5 m, and the focal length is 3 m. The Odessa's Atlas of light curves includes almost 5,5 thousand light curves for ~500 correlated objects from a time period of 2005-2014. The processing of light curves and the determination of the rotation period in the inertial frame is challenging. Extracted frequencies and reconstructed phases for some interesting targets, e.g. GLONASS satellites, for which also SLR data were available for confirmation, will be presented. The rotation of the Envisat satellite after its sudden failure will be analyzed. The deceleration of its rotation rate within 3 years is studied together with the attempt to determine the orientation of the rotation axis.

Bone-added periodontal plastic surgery: a new approach in esthetic dentistry

This article proposes a combined technique including bone grafting, connective tissue graft, and coronally advanced flap to create some space for simultaneous bone regrowth and root coverage. A 23 year-old female was referred to our private clinic with a severe class II Miller recession and lack of attached gingiva. The suggested treatment plan comprised of root coverage combined with xenograft bone particles. The grafted area healed well and full coverage was achieved at 12-month follow-up visit. Bone-added periodontal plastic surgery can be considered as a practical procedure for management of deep gingival recession without buccal bone plate.

Black Degradable Plastic Mulch Evaluation

Black polyethylene plastic mulch provides many benefits to fruit and vegetable producers. It increases earliness by increasing spring soil temperatures, conserves soil moisture, and reduces pesticide usage by decreasing weed and disease pressure. Furthermore, during seasons of high precipitation, it protects fertilizer from leaching below the root zone. Unfortunately, polyethylene mulches do not degrade and must be removed from the field and discarded each season. This is a labor-intensive process whether it is done mechanically or by hand. Several degradable plastic mulches have been developed that are designed to be incorporated into the soil profile, eliminating the need for removal, with no negative impact on soil quality or health. However, these degradable plastics often do not meet degradation expectations (either degrade too quickly or degrade incompletely and require manual removal). The objective of this project was to evaluate several degradable mulches for storage life, ease of use, and influence on tomato production.

Accumulation rates of volcanic glass and ice-rafted debris in sediments of the North Pacific Ocean

Mass accumulation rates (MAR) of different components of North Pacific deep-sea sediment provide detailed information about the timing of the onset of major Northern Hemisphere glaciation that occurred at 2.65 Ma. An increase in explosive volcanism in the Kamchatka-Kurile and Aleutian arcs occured at this same time, suggesting a link between volcanism and glaciation. Sediments recovered by piston-coring techniques during ODP Leg 145 provide a unique opportunity to undertake a detailed test of this possibility. Here we use volcanic glass as a proxy for explosive volcanism and ice-rafted debris (IRD) as a proxy for glaciation. The MAR of both glass and IRD increase markedly at 2.65 Ma. Further, the flux of the volcanic glass increased just prior the flix of ice-radted material, suggesting that the cooling resulting from explosive volcanic eruptions may have been the ultimate trigger for the mid-Pliocene glacial intensification.

Nd-Sr isotopes in fossil fish debris from IODP Exp302

Strontium and neodymium radiogenic isotope ratios in early to middle Eocene fossil fish debris (ichthyoliths) from Lomonosov Ridge (Integrated Ocean Drilling Program Expedition 302) help constrain water mass compositions in the Eocene Arctic Ocean between 55 and 45 Ma. The inferred paleodepositional setting was a shallow, offshore marine to marginal marine environment with limited connections to surrounding ocean basins. The new data demonstrate that sources of Nd and Sr in fish debris were distinct from each other, consistent with a salinity-stratified water column above Lomonosov Ridge in the Eocene. The 87Sr/86Sr values of ichthyoliths (0.7079 - 0.7087) are more radiogenic than Eocene seawater, requiring brackish to fresh water conditions in the environment where fish metabolized Sr. The 87Sr/86Sr variations probably record changes in the overall balance of river Sr flux to the Eocene Arctic Ocean between 55 and 45 Ma and are used here to reconstruct surface water salinity values. The eNd values of ichthyoliths vary between -5.7 and -7.8, compatible with periodic (or intermittent) supply of Nd to Eocene Arctic intermediate water (AIW) from adjacent seas. Although the Norwegian-Greenland Sea and North Atlantic Ocean were the most likely sources of Eocene AIW Nd, input from the Tethys Sea (via the Turgay Strait in early Eocene time) and the North Pacific Ocean (via a proto-Bering Strait) also contributed.

Ultimate and failure strengths of normal muds and serpentinite muds of ODP Leg 125 samples

Serpentinite seamounts in the Mariana forearc have been explained as diapirs rising from the Benioff zone. This hypothesis predicts that the serpentinites should have low strengths as well as low densities relative to the surrounding rocks. Drilling during Leg 125 showed that the materials forming Conical Seamount in the Mariana forearc and Torishima Forearc Seamount in the Izu-Bonin forearc are water-charged serpentinite muds of density <2 g/cm**3. Wykeham-Farrance torsion-vane tests showed that they are plastic solids with a rheology that bears many similarities to the idealized Cam clay soil model and is well described by critical-state soil mechanics. The serpentinite muds have ultimate strengths of 1.3 to 273.7 kPa and yield strengths of approximately 1.0 to 50 kPa. These muds thus are orders of magnitude weaker than salt and are, in fact, comparable in density and strength to common deep-sea clay muds. Such weak and low-density materials easily become diapiric. Serpentinite muds from the summit of Conical Seamount are weaker and more ductile than those on its flanks or on Torishima Forearc Seamount. Moreover, the summit muds do not contain the pronounced pure- and simple-shear fabrics that characterize those on the seamount flanks. The seamounts are morphologically similar to shield volcanoes, and anastomosing serpentinite debris flows descending from their summits are similar in map view to pahoehoe flows. These morphologic features, together with the physical properties of the muds and their similarities to other oceanic muds and the geochemistry of the entrained waters, suggest that many forearc serpentinite seamounts are gigantic (10-20 km wide, 1.5-2.0 km high) mud volcanoes that formed by the eruption of highly liquid serpentinite muds. Torishima Forearc Seamount, which is blanketed by more ìnormalî pelagic/volcaniclastic sediment, has probably been inactive since the Miocene. Conical Seamount, which seems to consist entirely of serpentinite mud and is venting fresh water of unusual chemistry from its summit, is presently active. Muds from the flanks of Conical Seamount are stronger and more brittle than those from the summit site, and muds from Torishima Forearc Seamount are stronger yet; this suggests that the serpentinite debris flows are compacted and dewatered as they mature. The shear fabrics probably result from downslope creep and flow, but may also be inherited.