SMARTnet: First Experience of Setting Up a Telescope System to Survey the Geostationary Ring

Space debris in geostationary orbits may be detected with optical telescopes when the objects are illuminated by the Sun. The advantage compared to Radar can be found in the illumination: radar illuminates the objects and thus the detection sensitivity depletest proportional to the fourth power of the d istance. The German Space Operation Center, GSOC, together with the Astronomical Institute of the University of Bern, AIUB, are setting up a telescope system called SMARTnet to demonstrate the capability of performing geostationary surveillance. Such a telescope system will consist of two telescopes on one mount: a smaller telescope with an aperture of 20cm will serve for fast survey while the larger one, a telescope with an aperture of 50cm, will be used for follow-up observations. The telescopes will be operated by GSOC from Oberpfaffenhofen by the internal monitoring and control system called SMARTnetMAC. The observation plan will be generated by MARTnetPlanning seven days in advance by applying an optimized planning scheduler, taking into account fault time like cloudy nights, priority of objects etc. From each picture taken, stars will be identified and everything not being a star is treated as a possible object. If the same object can be identified on multiple pictures within a short time span, the trace is called a tracklet. In the next step, several tracklets will be correlated to identify individual objects, ephemeris data for these objects are generated and catalogued . This will allow for services like collision avoidance to ensure safe operations for GSOC’s satellites. The complete data processing chain is handled by BACARDI, the backbone catalogue of relational debris information and is presented as a poster.

A covariance analysis to optimize the optical follow-up strategies

The Astronomical Institute of the University of Bern (AIUB) is conducting several search campaigns for space debris using optical sensors. The debris objects are discovered during systematic survey observations. In general, the result of a discovery consists in only a short observation arc, or tracklet, which is used to perform a first orbit determination in order to be able to observe t he object again in subsequent follow-up observations. The additional observations are used in the orbit improvement process to obtain accurate orbits to be included in a catalogue. In order to obtain the most accurate orbit within the time available it is necessary to optimize the follow-up observations strategy. In this paper an in‐depth study, using simulations and covariance analysis, is performed to identify the optimal sequence of follow-up observations to obtain the most accurate orbit propagation to be used for the space debris catalogue maintenance. The main factors that determine the accuracy of the results of an orbit determination/improvement process are: tracklet length, number of observations, type of orbit, astrometric error of the measurements, time interval between tracklets, and the relative position of the object along its orbit with respect to the observing station. The main aim of the covariance analysis is to optimize the follow-up strategy as a function of the object-observer geometry, the interval between follow-up observations and the shape of the orbit. This an alysis can be applied to every orbital regime but particular attention was dedicated to geostationary, Molniya, and geostationary transfer orbits. Finally the case with more than two follow-up observations and the influence of a second observing station are also analyzed.

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.

Relative estimate of accuracy in orbit determination from a short tracklet

The Astronomical Institute of the University of Bern (AIUB) is conducting several search campaigns for orbital debris. The debris objects are discovered during systematic survey observations. In general only a short observation arc, or tracklet, is available for most of these objects. From this discovery tracklet a first orbit determination is computed in order to be able to find the object again in subsequent follow-up observations. The additional observations are used in the orbit improvement process to obtain accurate orbits to be included in a catalogue. In this paper, the accuracy of the initial orbit determination is analyzed. This depends on a number of factors: tracklet length, number of observations, type of orbit, astrometric error, and observation geometry. The latter is characterized by both the position of the object along its orbit and the location of the observing station. Different positions involve different distances from the target object and a different observing angle with respect to its orbital plane and trajectory. The present analysis aims at optimizing the geometry of the discovery observation is depending on the considered orbit.

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.

The M-theory origin of global properties of gauge theories

We show that global properties of gauge groups can be understood as geometric properties in M-theory. Different wrappings of a system of N M5-branes on a torus reduce to four-dimensional theories with AN−1 gauge algebra and different unitary groups. The classical properties of the wrappings determine the global properties of the gauge theories without the need to impose any quantum conditions. We count the inequivalent wrappings as they fall into orbits of the modular group of the torus, which correspond to the S-duality orbits of the gauge theories.

Narrow dust jets in a diffuse gas coma: a natural product of small active regions on comets

Comets often display narrow dust jets but more diffuse gas comae when their eccentric orbits bring them into the inner solar system and sunlight sublimates the ice on the nucleus. Comets are also understood to have one or more active areas covering only a fraction of the total surface active with sublimating volatile ices. Calculations of the gas and dust distribution from a small active area on a comet’s nucleus show that as the gas moves out radially into the vacuum of space it expands tangentially, filling much of the hemisphere centered on the active region. The dust dragged by the gas remains more concentrated over the active area. This explains some puzzling appearances of comets having collimated dust jets but more diffuse gaseous atmospheres. Our test case is 67P/Churyumov–Gerasimenko, the Rosetta mission target comet, whose activity is dominated by a single area covering only 4% of its surface.

Modeled Interaction of Comet 67P/Churyumov-Gerasimenko with the Solar Wind Inside 2 AU

Periodic comets move around the Sun on elliptical orbits. As such comet 67P/Churyumov-Gerasimenko (hereafter 67P) spends a portion of time in the inner solar system where it is exposed to increased solar insolation. Therefore given the change in heliocentric distance, in case of 67P from aphelion at 5.68 AU to perihelion at ~1.24 AU, the comet’s activity—the production of neutral gas and dust—undergoes significant variations. As a consequence, during the inbound portion, the mass loading of the solar wind increases and extends to larger spatial scales. This paper investigates how this interaction changes the character of the plasma environment of the comet by means of multifluid MHD simulations. The multifluid MHD model is capable of separating the dynamics of the solar wind ions and the pick-up ions created through photoionization and electron impact ionization in the coma of the comet. We show how two of the major boundaries, the bow shock and the diamagnetic cavity, form and develop as the comet moves through the inner solar system. Likewise for 67P, although most likely shifted back in time with respect to perihelion passage, this process is reversed on the outbound portion of the orbit. The presented model herein is able to reproduce some of the key features previously only accessible to particle-based models that take full account of the ions’ gyration. The results shown herein are in decent agreement to these hybrid-type kinetic simulations.

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.

CODE rapid product series for the IGS

CODE, the Center for Orbit Determination in Europe, is a joint venture of the following four institutions: Astronomical Institute, University of Bern (AIUB), Bern, Switzerland;Federal Office of Topography swisstopo, Wabern, Switzerland; Federal Agency of Cartography and Geodesy (BKG), Frankfurt a. M., Germany; Institut für Astronomische und Physikalische Geodäsie, Technische Universität München (IAPG, TUM), Munich, Germany. It acts as a global analysis center of the International GNSS Service (IGS). The operational computations are performed at AIUB using the latest development version of the Bernese GNSS Software (Dach et al., 2015). In this context a rapid solution series is generated considering all active GPS and GLONASS satellites. It contains 24 hours of observed orbits and published at the day after the observations.

The CODE MGEX orbit and clock solution

The Center for Orbit Determination in Europe (CODE) is contributing as a global Analysis center to the International GNSS Service (IGS) since many years. The processing of GPS and GLONASS data is well established in CODE’s ultra-rapid, rapid, and final product lines. With the introduction of new signals for the established and new GNSS, new challenges and opportunities are arising for the GNSS data management and processing. The IGS started the Multi-GNSS-EXperiment (MGEX) in 2012 in order to gain first experience with the new data formats and to develop new strategies for making optimal use of these additional measurements. CODE has started to contribute to IGS MGEX with a consistent, rigorously combined triple-system orbit solution (GPS, GLONASS, and Galileo). SLR residuals for the computed Galileo satellite orbits are of the order of 10 cm. Furthermore CODE established a GPS and Galileo clock solution. A quality assessment shows that these experimental orbit and clock products allow even a Galileo-only precise point positioning (PPP) with accuracies on the decimeter- (static PPP) to meter-level (kinematic PPP) for selected stations.

Impact of the arc length on GNSS analysis results

Homogeneously reprocessed combined GPS/GLONASS 1- and 3-day solutions from 1994 to 2013, generated by the Center for Orbit Determination in Europe (CODE) in the frame of the second reprocessing campaign REPRO-2 of the International GNSS Service, as well as GPS- and GLONASS-only 1- and 3-day solutions for the years 2009 to 2011 are analyzed to assess the impact of the arc length on the estimated Earth Orientation Parameters (EOP, namely polar motion and length of day), on the geocenter, and on the orbits. The conventional CODE 3-day solutions assume continuity of orbits, polar motion components, and of other parameters at the day boundaries. An experimental 3-day solution, which assumes continuity of the orbits, but independence from day to day for all other parameters, as well as a non-overlapping 3-day solution, is included into our analysis. The time series of EOPs, geocenter coordinates, and orbit misclosures, are analyzed. The long-arc solutions were found to be superior to the 1-day solutions: the RMS values of EOP and geocenter series are typically reduced between 10 and 40 %, except for the polar motion rates, where RMS reductions by factors of 2–3 with respect to the 1-day solutions are achieved for the overlapping and the non-overlapping 3-day solutions. In the low-frequency part of the spectrum, the reduction is even more important. The better performance of the orbits of 3-day solutions with respect to 1-day solutions is also confirmed by the validation with satellite laser ranging.

Search for satellites near comet 67P/Churyumov-Gerasimenko using Rosetta/OSIRIS images

Context. The European Space Agency Rosetta mission reached and started escorting its main target, the Jupiter-family comet 67P/Churyumov-Gerasimenko, at the beginning of August 2014. Within the context of solar system small bodies, satellite searches from approaching spacecraft were extensively used in the past to study the nature of the visited bodies and their collisional environment. Aims. During the approaching phase to the comet in July 2014, the OSIRIS instrument onboard Rosetta performed a campaign aimed at detecting objects in the vicinity of the comet nucleus and at measuring these objects' possible bound orbits. In addition to the scientific purpose, the search also focused on spacecraft security to avoid hazardous material in the comet's environment. Methods. Images in the red spectral domain were acquired with the OSIRIS Narrow Angle Camera, when the spacecraft was at a distance between 5785 km and 5463 km to the comet, following an observational strategy tailored to maximize the scientific outcome. From the acquired images, sources were extracted and displayed to search for plausible displacements of all sources from image to image. After stars were identified, the remaining sources were thoroughly analyzed. To place constraints on the expected displacements of a potential satellite, we performed Monte Carlo simulations on the apparent motion of potential satellites within the Hill sphere. Results. We found no unambiguous detections of objects larger than similar to 6 m within similar to 20 km and larger than similar to 1 m between similar to 20 km and similar to 110 km from the nucleus, using images with an exposure time of 0.14 s and 1.36 s, respectively. Our conclusions are consistent with independent works on dust grains in the comet coma and on boulders counting on the nucleus surface. Moreover, our analysis shows that the comet outburst detected at the end of April 2014 was not strong enough to eject large objects and to place them into a stable orbit around the nucleus. Our findings underline that it is highly unlikely that large objects survive for a long time around cometary nuclei.

Orbital elements of the material surrounding comet 67P/Churyumov-Gerasimenko

Context. We investigate the dust coma within the Hill sphere of comet 67P/Churyumov-Gerasimenko. Aims. We aim to determine osculating orbital elements for individual distinguishable but unresolved slow-moving grains in the vicinity of the nucleus. In addition, we perform photometry and constrain grain sizes. Methods. We performed astrometry and photometry using images acquired by the OSIRIS Wide Angle Camera on the European Space Agency spacecraft Rosetta. Based on these measurements, we employed standard orbit determination and orbit improvement techniques. Results. Orbital elements and effective diameters of four grains were constrained, but we were unable to uniquely determine them. Two of the grains have light curves that indicate grain rotation. Conclusions. The four grains have diameters nominally in the range 0.14-0.50 m. For three of the grains, we found elliptic orbits, which is consistent with a cloud of bound particles around the nucleus. However, hyperbolic escape trajectories cannot be excluded for any of the grains, and for one grain this is the only known option. One grain may have originated from the surface shortly before observation. These results have possible implications for the understanding of the dispersal of the cloud of bound debris around comet nuclei, as well as for understanding the ejection of large grains far from the Sun.

Regional surface morphology of comet 67P/Churyumov-Gerasimenko from Rosetta/OSIRIS images

Aims. The OSIRIS camera onboard the Rosetta spacecraft has been acquiring images of the comet 67P/Churyumov-Gerasimenko (67P)'s nucleus at spatial resolutions down to similar to 0.17 m/px ever since Aug. 2014. These images have yielded unprecedented insight into the morphological diversity of the comet's surface. This paper presents an overview of the regional morphology of comet 67P. Methods. We used the images that were acquired at orbits similar to 20-30 km from the center of the comet to distinguish different regions on the surface and introduce the basic regional nomenclature adopted by all papers in this Rosetta special feature that address the comet's morphology and surface processes. We used anaglyphs to detect subtle regional and topographical boundaries and images from close orbit (similar to 10 km from the comet's center) to investigate the fine texture of the surface. Results. Nineteen regions have currently been defined on the nucleus based on morphological and/or structural boundaries, and they can be grouped into distinctive region types. Consolidated, fractured regions are the most common region type. Some of these regions enclose smooth units that appear to settle in gravitational sinks or topographically low areas. Both comet lobes have a significant portion of their surface covered by a dusty coating that appears to be recently placed and shows signs of mobilization by aeolian-like processes. The dusty coatings cover most of the regions on the surface but are notably absent from a couple of irregular large depressions that show sharp contacts with their surroundings and talus-like deposits in their interiors, which suggests that short-term explosive activity may play a significant role in shaping the comet's surface in addition to long-term sublimation loss. Finally, the presence of layered brittle units showing signs of mechanical failure predominantly in one of the comet's lobes can indicate a compositional heterogeneity between the two lobes.