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.

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.

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

Charge exchange in cometary coma: Discovery of H- ions in the solar wind close to comet 67P/Churyumov-Gerasimenko

As Rosetta was orbiting comet 67P/Churyumov-Gerasimenko, the Ion and Electron Sensor detected negative particles with angular distributions like those of the concurrently measured solar wind protons but with fluxes of only about 10% of the proton fluxes and energies of about 90% of the proton energies. Using well-known cross sections and energy-loss data, it is determined that the fluxes and energies of the negative particles are consistent with the production of H- ions in the solar wind by double charge exchange with molecules in the coma.

Fractures on comet 67P/Churyumov-Gerasimenko observed by Rosetta/OSIRIS

The Optical, Spectroscopic, and Infrared Remote Imaging System (OSIRIS) experiment onboard the Rosetta spacecraft currently orbiting comet 67P/Churyumov-Gerasimenko has yielded unprecedented views of a comet's nucleus. We present here the first ever observations of meter-scale fractures on the surface of a comet. Some of these fractures form polygonal networks. We present an initial assessment of their morphology, topology, and regional distribution. Fractures are ubiquitous on the surface of the comet's nucleus. Furthermore, they occur in various settings and show different topologies suggesting numerous formation mechanisms, which include thermal insulation weathering, orbital-induced stresses, and possibly seasonal thermal contraction. However, we conclude that thermal insolation weathering is responsible for creating most of the observed fractures based on their morphology and setting in addition to thermal models that indicate diurnal temperature ranges exceeding 200K and thermal gradients of similar to 15K/min at perihelion are possible. Finally, we suggest that fractures could be a facilitator in surface evolution and long-term erosion.

Spectrophotometric properties of the nucleus of comet 67P/Churyumov-Gerasimenko from the OSIRIS instrument onboard the ROSETTA spacecraft

Context. The Rosetta mission of the European Space Agency has been orbiting the comet 67P/Churyumov-Gerasimenko (67P) since August 2014 and is now in its escort phase. A large complement of scientific experiments designed to complete the most detailed study of a comet ever attempted are onboard Rosetta. Aims. We present results for the photometric and spectrophotometric properties of the nucleus of 67P derived from the OSIRIS imaging system, which consists of a Wide Angle Camera (WAC) and a Narrow Angle Camera (NAC). The observations presented here were performed during July and the beginning of August 2014, during the approach phase, when OSIRIS was mapping the surface of the comet with several filters at different phase angles (1.3 degrees-54 degrees). The resolution reached up to 2.1 m/px. Methods. The OSIRIS images were processed with the OSIRIS standard pipeline, then converted into I/F. radiance factors and corrected for the illumination conditions at each pixel using the Lommel-Seeliger disk law. Color cubes of the surface were produced by stacking registered and illumination-corrected images. Furthermore, photometric analysis was performed both on disk-averaged photometry in several filters and on disk-resolved images acquired with the NAC orange filter, centered at 649 ran, using Hapke modeling. Results. The disk-averaged phase function of the nucleus of 67P shows a strong opposition surge with a G parameter value of -0.13 +/- 0.01 in the HG system formalism and an absolute magnitude H-v(1, 1, 0) = 15.74 +/- 0.02 mag. The integrated spectrophotometry in 20 filters covering the 250-1000 nm wavelength range shows a red spectral behavior, without clear absorption bands except for a potential absorption centered at similar to 290 rim that is possibly due to SO2 ice. The nucleus shows strong phase reddening, with disk-averaged spectral slopes increasing from 11%/( 100 nm) to 16%/(100 nm) in the 1.3 degrees-54 degrees phase angle range. The geometric albedo of the comet is 6.5 +/- 0.2% at 649 nm, with local variations of up to similar to 16% in the Hapi region. From the disk-resolved images we computed the spectral slope together with local spectrophotometry and identified three distinct groups of regions (blue, moderately red, and red). The Hapi region is the brightest, the bluest in term of spectral slope, and the most active surface on the comet. Local spectrophotometry shows an enhancement of the flux in the 700-750 nm that is associated with coma emissions.

Interstellar Flow and Temperature Determination with IBEX: Robustness and Sensitivity to Systematic Effects

The Interstellar Boundary Explorer (IBEX) samples the interstellar neutral (ISN) gas flow of several species every year from December through late March when the Earth moves into the incoming flow. The first quantitative analyses of these data resulted in a narrow tube in four-dimensional interstellar parameter space, which couples speed, flow latitude, flow longitude, and temperature, and center values with approximately 3° larger longitude and 3 km s⁻¹ lower speed, but with temperatures similar to those obtained from observations by the Ulysses spacecraft. IBEX has now recorded six years of ISN flow observations, providing a large database over increasing solar activity and using varying viewing strategies. In this paper, we evaluate systematic effects that are important for the ISN flow vector and temperature determination. We find that all models in use return ISN parameters well within the observational uncertainties and that the derived ISN flow direction is resilient against uncertainties in the ionization rate. We establish observationally an effective IBEX-Lo pointing uncertainty of ±0°18 in spin angle and confirm an uncertainty of ±0°1 in longitude. We also show that the IBEX viewing strategy with different spin-axis orientations minimizes the impact of several systematic uncertainties, and thus improves the robustness of the measurement. The Helium Warm Breeze has likely contributed substantially to the somewhat different center values of the ISN flow vector. By separating the flow vector and temperature determination, we can mitigate these effects on the analysis, which returns an ISN flow vector very close to the Ulysses results, but with a substantially higher temperature. Due to coupling with the ISN flow speed along the ISN parameter tube, we provide the temperature Tvisn∞=8710+440/-680 K for Visn∞=26 km s⁻¹ for comparison, where most of the uncertainty is systematic and likely due to the presence of the Warm Breeze.

Constraining planet structure from stellar chemistry: the cases of CoRoT-7, Kepler-10, and Kepler-93

Aims. We explore the possibility that the stellar relative abundances of different species can be used to constrain the bulk abundances of known transiting rocky planets. Methods. We use high resolution spectra to derive stellar parameters and chemical abundances for Fe, Si, Mg, O, and C in three stars hosting low mass, rocky planets: CoRoT-7, Kepler-10, and Kepler-93. These planets follow the same line along the mass-radius diagram, pointing toward a similar composition. The derived abundance ratios are compared with the solar values. With a simple stoichiometric model, we estimate the iron mass fraction in each planet, assuming stellar composition. Results. We show that in all cases, the iron mass fraction inferred from the mass-radius relationship seems to be in good agreement with the iron abundance derived from the host star's photospheric composition. Conclusions. The results suggest that stellar abundances can be used to add constraints on the composition of orbiting rocky planets.

On the nucleus structure and activity of comet 67P/Churyumov-Gerasimenko

Images from the OSIRIS scientific imaging system onboard Rosetta show that the nucleus of 67P/Churyumov-Gerasimenko consists of two lobes connected by a short neck. The nucleus has a bulk density less than half that of water. Activity at a distance from the Sun of >3 astronomical units is predominantly from the neck, where jets have been seen consistently. The nucleus rotates about the principal axis of momentum. The surface morphology suggests that the removal of larger volumes of material, possibly via explosive release of subsurface pressure or via creation of overhangs by sublimation, may be a major mass loss process. The shape raises the question of whether the two lobes represent a contact binary formed 4.5 billion years ago, or a single body where a gap has evolved via mass loss.