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.

Liulin silicon semiconductor spectrometers as cosmic ray monitors at the high mountain observatories Jungfraujoch and Lomnický stit

Currently, most cosmic ray data are obtained by detectors on satellites, aircraft, high-altitude balloons and ground (neutron monitors). In our work, we examined whether Liulin semiconductor spectrometers (simple silicon planar diode detectors with spectrometric properties) located at high mountain observatories could contribute new information to the monitoring of cosmic rays by analyzing data from selected solar events between 2005 and 2013. The decision thresholds and detection limits of these detectors placed at Jungfraujoch (Switzerland; 3475 m a.s.l.; vertical cut-off rigidity 4.5 GV) and Lomnicky stıt (Slovakia; 2633 m a.s.l.; vertical cut-off rigidity 3.84 GV) highmountain observatories were determined. The data showed that only the strongest variations of the cosmic ray flux in this period were detectable. The main limitation in the performance of these detectors is their small sensitive volume and low sensitivity of the PIN photodiode to neutrons.

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.

Trajectory mapping of middle atmospheric water vapor by a mini network of NDACC instruments

The important task to observe the global coverage of middle atmospheric trace gases like water vapor or ozone usually is accomplished by satellites. Climate and atmospheric studies rely upon the knowledge of trace gas distributions throughout the stratosphere and mesosphere. Many of these gases are currently measured from satellites, but it is not clear whether this capability will be maintained in the future. This could lead to a significant knowledge gap of the state of the atmosphere. We explore the possibilities of mapping middle atmospheric water vapor in the Northern Hemisphere by using Lagrangian trajectory calculations and water vapor profile data from a small network of five ground-based microwave radiometers. Four of them are operated within the frame of NDACC (Network for the Detection of Atmospheric Composition Change). Keeping in mind that the instruments are based on different hardware and calibration setups, a height-dependent bias of the retrieved water vapor profiles has to be expected among the microwave radiometers. In order to correct and harmonize the different data sets, the Microwave Limb Sounder (MLS) on the Aura satellite is used to serve as a kind of traveling standard. A domain-averaging TM (trajectory mapping) method is applied which simplifies the subsequent validation of the quality of the trajectory-mapped water vapor distribution towards direct satellite observations. Trajectories are calculated forwards and backwards in time for up to 10 days using 6 hourly meteorological wind analysis fields. Overall, a total of four case studies of trajectory mapping in different meteorological regimes are discussed. One of the case studies takes place during a major sudden stratospheric warming (SSW) accompanied by the polar vortex breakdown; a second takes place after the reformation of stable circulation system. TM cases close to the fall equinox and June solstice event from the year 2012 complete the study, showing the high potential of a network of ground-based remote sensing instruments to synthesize hemispheric maps of water vapor.

Sublimation of water ice mixed with silicates and tholins: Evolution of surface texture and reflectance spectra, with implications for comets

The surfaces of many objects in the Solar System comprise substantial quantities of water ice sometimes mixed with minerals and/or organic molecules. The sublimation of the ice changes the structural and optical properties of these objects. We present laboratory data on the evolution of the structure and the visible and near-infrared spectral reflectance of icy surface analogues of cometary ices, made of water ice, complex organic matter (tholins) and silicates, as they undergo sublimation under low temperature (<-70°C) and pressure (10-⁵mbar) conditions inside the SCITEAS simulation chamber. As the water ice sublimated, we observed in situ the formation of a porous sublimation lag deposit, or sublimation mantle, at the top of the ice. This mantle is a network of filaments made of the non-volatile particles. Organics or phyllosilicates grains, able to interact via stronger inter-particulate forces than olivine grains, can form a foam-like structure having internal cohesiveness, holding olivine grains together. As this mantle builds-up, the band depths of the sub-surface water ice are attenuated until complete extinction under only few millimeters of mantle. Optically thick sublimation mantles are mainly featureless in the near infrared. The absorption bands of the minerals present in the mantle are weak, or even totally absent if minerals are mixed with organics which largely dominate the VIS–NIR reflectance spectrum. During sublimation, ejections of large fragments of mantle, triggered by the gas flow, expose ice particles to the surface. The contrast of brightness between mantled and ice-exposed areas depends on the wavelength range and the dust/ice ratio considered. We describe how the chemical nature of the non-volatiles, the size of their particles, the way they are mixed with the ice and the dust/ice mass ratio influence the texture, activity and spectro-photometric properties of the sublimation mantles. These data provide useful references for interpreting remote-sensing observations of comets and also icy satellites or trans-neptunian objects.

CODE ultra-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. In this context an ultra-rapid solution series is generated considering GPS and GLONASS satellites. It is updated several times per day and contains 24 hours of observed and 24 hours of predicted orbit interval. More details are available in: Lutz, S., G. Beutler, S. Schaer, R. Dach, A. Jäggi; 2014: CODE's new ultra-rapid orbit and ERP products for the IGS. GPS Solutions. DOI 10.1007/s10291-014-0432-2

CODE repro2 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, Dow et al, 2009). The operational computations are performed at AIUB using the latest development version of the Bernese GNSS Software (Dach et al., 2015). In this context the contribution to the IGS repro02 effort is generated considering only the GPS satellites between 1994 and 2001 as well as the GPS and GLONASS satellites from 2002 to the end of 2013.

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.

CODE final 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. In this context a final solution series is generated considering all active GPS and GLONASS satellites. It is published in daily files with a delay of about two weeks.

CODE product series for the IGS-MGEX project

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. In this context a multi-GNSS solution is generated considering all active GPS, GLONASS, Galileo, BeiDou (expect for GEOs), and QZSS satellites as a contribution to the IGS-MGEX project. The results are published with a delay of about two weeks.

Mapping spectroscopic uncertainties into prospective methane retrieval errors from Sentinel-5 and its precursor

Sentinel-5 (S5) and its precursor (S5P) are future European satellite missions aiming at global monitoring of methane (CH4) column-average dry air mole fractions (XCH4). The spectrometers to be deployed onboard the satellites record spectra of sunlight backscattered from the Earth's surface and atmosphere. In particular, they exploit CH4 absorption in the shortwave infrared spectral range around 1.65 mu m (S5 only) and 2.35 mu m (both S5 and S5P) wavelength. Given an accuracy goal of better than 2% for XCH4 to be delivered on regional scales, assessment and reduction of potential sources of systematic error such as spectroscopic uncertainties is crucial. Here, we investigate how spectroscopic errors propagate into retrieval errors on the global scale. To this end, absorption spectra of a ground-based Fourier transform spectrometer (FTS) operating at very high spectral resolution serve as estimate for the quality of the spectroscopic parameters. Feeding the FTS fitting residuals as a perturbation into a global ensemble of simulated S5- and S5P-like spectra at relatively low spectral resolution, XCH4 retrieval errors exceed 0.6% in large parts of the world and show systematic correlations on regional scales, calling for improved spectroscopic parameters.

Estimation of satellite antenna phase center offsets for Galileo

Satellite antenna phase center offsets for the GalileoInOrbitValidation(IOV) and FullOperationalCapability (FOC) satellites are estimated by two different analysiscenters based on tracking data of a global GNSS network. The mean x- and y-offsets could be determined with a precision of a few centimeters. However, daily estimates of thex-offsets of the IOV satellites show pronounced systematic effects with a peak-to-peak amplitude of up to 70 cm that depend on the orbit model and the elevation of the Sun above the orbital plane. For the IOV y-offsets, no dependence on the orbit model exists but the scatter strongly depends on the elevation of the Sun above the orbital plane. In general, these systematic effects are significantly smaller for the FOC satellites. The z-offsets of the two analysis centers agree within the 10–15 cm level, and the time series do not show systematic effects. The application of an averaged Galileo satellite antenna model obtained from the two solutions results in a reduction of orbit day boundary discontinuities by up to one third—even if an independent software package is used.

CODE final 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. In this context a final solution series is generated considering all active GPS and GLONASS satellites. It is published in daily files with a delay of about two weeks.