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.

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.

Ice-rafted debris record in the Prydz Bay

Ice-rafted debris mass accumulation rates (IRD MAR) at a drill site on the Antarctic continental margin are investigated to evaluate the linkages between East Antarctic Ice Sheet extent and Southern Ocean temperatures in the early to mid-Pliocene. ODP Site 1165 is within 400 km of the Antarctic coastline and in the direct pathway of icebergs released by the Amery Ice Shelf. The Amery Ice Shelf is the largest ice shelf in East Antarctica and it buttresses the Lambert Glacier drainage system, which accounts for 14% of the outflow from the East Antarctic Ice Sheet. IRD MAR were low during peak Southern Ocean warming in the early Pliocene. After a brief precursor, a tenfold increase in IRD MAR at 3.3 Ma marks the termination of the early Pliocene ice sheet minimum, coincident with the M2 glacial. For the mid-Pliocene, a strong correlation exists between the high-amplitude signal in the LR04 benthic stack and IRD MAR, suggesting linkages between East Antarctic ice extent, global ice volume and deep-water temperatures. The IRD record at Site 1165 provides evidence of greater sensitivity of the Lambert Glacier-Amery Ice Shelf system to Southern Ocean warming than is currently predicted by ice sheet models, which may relate to uncertainties in the understanding of ocean heat uptake, poleward heat transport and ice sheet-ocean interactions.

Neogene ice-rafted debris record of ODP Hole 120-751A

One of the primary objectives of Leg 120 was to obtain a high-resolution Neogene stratigraphic section from the Kerguelen Plateau. Site 751, located in the central part of the Raggatt Basin on the Southern Kerguelen Plateau in 1633.8 m of water (57°43.56'S; 79°48.89'E), was selected as the dedicated Neogene site for this objective. High-resolution sampling at Site 751 was used to delineate in detail the Neogene ice-rafted debris (IRD) occurrences on the Kerguelen Plateau. The oldest IRD found at Site 751 was approximately 9.9 Ma, and it was not until approximately 8.5 Ma that significant concentrations of IRD were detected. The first major IRD event at this site occurred in the uppermost Miocene between 6.0 and 5.5 Ma. During this time period, a general climatic cooling and glacial expansion occurred on Antarctica. The late Miocene IRD event was followed by a continuous episode of elevated IRD deposition in the lowermost Pliocene between 4.5 and 4.1 Ma. The 0.4-m.y. duration and the timing of the early Pliocene IRD event on the Kerguelen Plateau corresponds with IRD fluxes observed on the Falkland Plateau and in the Weddell Abyssal Plain. This correspondence of data indicates that a major global climatic event occurred during the early Pliocene. The East Antarctic Ice Sheet may have experienced deglaciation between 4.5 and 4.1 Ma and, as a result, released large volumes of sediment-laden ice into the Southern Ocean.

Coarse-sand ice-rafted debris components and mass accumulation rates of ODP Site 152-918

To understand the late Cenozoic glacial history of the Northern Hemisphere, continuous long-term proxy records from climatically sensitive regions must be examined. Ice-rafted debris (IRD) from Ocean Drilling Program (ODP) Site 918, located in the Irminger Basin, is one such record. IRD in marine sediments is a direct indicator of the presence of glacial ice extending to sea level on adjacent landmasses, and, therefore, is an important paleoclimatic signal from the mid- to high latitudes. The IRD record at Site 918 is the first long-term ice-rafting record available for southeast Greenland, a region that may have been a key nucleation area for widespread glaciation during the late Cenozoic (Larsen et al, 1994, doi:10.2973/odp.proc.ir.152.1994). This data report presents the results of coarse sand-size IRD mass accumulation rate (MAR) analyses for Site 918 from the late Miocene through the Pleistocene. In addition, a preliminary analysis of IRD compositions is included. Detailed discussions of the local, regional, and global paleoclimatic implications of this data, and of the companion Site 919 Pleistocene IRD MAR data (Krissek, 1999, doi:10.2973/odp.proc.sr.163.118.1999), are in preparation. Such future work will include comparisons of these IRD MAR data sets to the Site 919 oxygen isotope stratigraphy developed by Flower (1998, doi:10.2973/odp.proc.sr.152.219.1998).

Ice rafted debris in the Nordic Seas

A continuous 3.5 Myr IRD record was produced from Ocean Drilling Program (ODP) Site 907. A timescale based on magnetic polarity chrons, oxygen isotope stratigraphy (for the last 1Myr) and orbital tuning was developed. The record documents a stepwise inception of large-scale glacial cycles in the Nordic Seas region, the first being a marked expansion of the Greenland ice sheet at 3.3 Ma. A second step occurred at 2.74 Ma by an expansion of large scale ice sheets in the Northern Hemisphere. Ice sheet variability around the Nordic Seas was tightly coupled to global ice volume over the past 3.3 Myr. Between 3 and 1 Ma, most of the variance of the IRD signal is in the 41 kyr band, whereas the last 1 Myr is characterized by stronger 100 kyr variance. The Gamma Ray Porosity Evaluator (GRAPE) density record is closely linked with IRD variations and documents sub orbital variability resembling the late Quaternary Heinrich/Bond cycles.

Geochemistry and provenance of basaltic clasts within volcaniclastic debris flows at DSDP Site 89-585

Pebble-sized basaltic and glassy clasts were extracted from seamount-derived volcaniclastic debris flows and analyzed for various trace elements, including the rare earths, to determine their genetic relationships and provenance. All the clasts were originally derived from relatively shallow submarine lava flows prior to sedimentary reworking, and have undergone minor low-grade alteration. They are classified into three petrographic groups (A, B, and C) characterized by different phenocryst assemblages and variable abundances and ratios of incompatible elements. Group A (clast from Hole 585) is a hyaloclastite fragment which is olivine-normative and distinct from the other clasts, with incompatibleelement ratios characteristic of transitional or alkali basalts. Groups B and C (clasts from Hole 585A) are quartz-normative, variably plagioclase-clinopyroxene-olivine phyric tholeiites, all with essentially similar ratios of highly incompatible elements and patterns of enrichment in light rare earth elements (chrondrite-normalized). Variation within Groups B and C was governed by low-pressure fractionation of the observed phenocryst phases, whereas the most primitive compositions of each group may be related by variable partial melting of a common source. The clasts have intraplate chemical characteristics, although relative to oceanic hot-spot-related volcanics (e.g., Hawaiian tholeiites) they are marginally depleted in most incompatible elements. The source region was enriched in all incompatible elements, compared with a depleted mid-ocean-ridge basalt source.

Geochemical composition of surface snow, ice and debris from glaciers in Norway, Nepal and New Zealand

Alpine glacier samples were collected in four contrasting regions to measure supraglacial dust and debris geochemical composition. A total of 70 surface glacier ice, snow and debris samples were collected in 2009 and 2010 in Svalbard, Norway, Nepal and New Zealand. Trace elemental abundances in snow and ice samples were measured via inductively coupled plasma mass spectrometry (ICP-MS). Supraglacial debris mineral, bulk oxide and trace element composition were determined via X-ray diffraction (XRD) and X-ray fluorescence spectroscopy (XRF). A total of 45 elements and 10 oxide compound abundances are reported. The uniform data collection procedure, analytical measurement methods and geochemical comparison techniques are used to evaluate supraglacial dust and debris composition variability in the contrasting glacier study regions. Elemental abundances revealed sea salt aerosol and metal enrichment in Svalbard, low levels of crustal dust and marine influences to southern Norway, high crustal dust and anthropogenic enrichment in the Khumbu Himalayas, and sulfur and metals attributed to quiescent degassing and volcanic activity in northern New Zealand. Rare earth element and Al/Ti elemental ratios demonstrated distinct provenance of particulates in each study region. Ca/S elemental ratio data showed seasonal denudation in Svalbard and Norway. Ablation season atmospheric particulate transport trajectories were mapped in each of the study regions and suggest provenance pathways. The in situ data presented provides first order glacier surface geochemical variability as measured from four diverse alpine glacier regions. This geochemical surface glacier data is relevant to glaciologic ablation rate understanding as well as satellite atmospheric and land-surface mapping techniques currently in development.