An Instrument Design for Space-Based Optical Observation of Space Debris

Currently, observations of space debris are primarily performed with ground-based sensors. These sensors have a detection limit at some centimetres diameter for objects in Low Earth Orbit (LEO) and at about two decimetres diameter for objects in Geostationary Orbit (GEO). The few space-based debris observations stem mainly from in-situ measurements and from the analysis of returned spacecraft surfaces. Both provide information about mostly sub-millimetre-sized debris particles. As a consequence the population of centimetre- and millimetre-sized debris objects remains poorly understood. The development, validation and improvement of debris reference models drive the need for measurements covering the whole diameter range. In 2003 the European Space Agency (ESA) initiated a study entitled “Space-Based Optical Observation of Space Debris”. The first tasks of the study were to define user requirements and to develop an observation strategy for a space-based instrument capable of observing uncatalogued millimetre-sized debris objects. Only passive optical observations were considered, focussing on mission concepts for the LEO, and GEO regions respectively. Starting from the requirements and the observation strategy, an instrument system architecture and an associated operations concept have been elaborated. The instrument system architecture covers the telescope, camera and onboard processing electronics. The proposed telescope is a folded Schmidt design, characterised by a 20 cm aperture and a large field of view of 6°. The camera design is based on the use of either a frame-transfer charge coupled device (CCD), or on a cooled hybrid sensor with fast read-out. A four megapixel sensor is foreseen. For the onboard processing, a scalable architecture has been selected. Performance simulations have been executed for the system as designed, focussing on the orbit determination of observed debris particles, and on the analysis of the object detection algorithms. In this paper we present some of the main results of the study. A short overview of the user requirements and observation strategy is given. The architectural design of the instrument is discussed, and the main tradeoffs are outlined. An insight into the results of the performance simulations is provided.

Orbital manifestation of whipple's disease: an atypical case

BACKGROUND: Whipple's disease is a systemic disorder caused by an infection with a gram-positive bacillus, Tropheryma whipplei. Almost every organ system can be affected in Whipple's disease, resulting in varying clinical symptoms. CASE REPORT: As far as we are aware, this report of a 61-year-old male is the first presenting with a periorbital manifestation of the disease, with severe exophthalmos and optic nerve involvement, leading to rapid visual loss. This emergency case was successfully treated by a surgical orbital decompression combined with systemic use of antibiotics and steroids. CONCLUSION: Whipple's disease can affect the periorbital tissues and the optic nerve, causing massive exophthalmos and serious transient visual loss. In such a case surgical decompression of the affected orbit combined with antibiotics and steroids is a recommended valid treatment option.

Late-onset posttraumatic septal hematoma and abscess formation in a six-year-old Tamil girl--case report and literature review

Nasal septal hematoma with abscess (NSHA) is an uncommon complication of trauma and studies on children are especially rare. We discuss the case of a 6-year-old girl, who was initially evaluated independently by three doctors for minor nasal trauma but had to be re-hospitalized 6 days later with NSHA. Although septal hematoma had initially been excluded (5, 7 and 24 hours after trauma), a secondary accumulation of blood seems to have occured. Delayed hematoma formation has been described in the orbit as a result of possible venous injuries after endoscopic sinus surgery. However, such an observation is new for septal hematoma in children. Thus, we recommend re-evaluation for septal hematoma 48h to 72h after paediatric nasal trauma. Such a scheduled re-examination offers a chance to treat delayed subperichondral hematoma on time before almost inevitable superinfection leads to abscess formation and destruction of the nasal infrastructure. We suggest that parents should be vigilant for delayed nasal obstruction as possible herald of hematoma accumulation within the first week.

GPS-only gravity field determination from GOCE data

The Gravity field and steady-state Ocean Circulation Explorer (GOCE) is now in orbit for more than four years. This is longer than the originally planned lifetime of the satellite and after three years on the same altitude the satellite has been lowered to 235 km in several steps. In the frame of the GOCE High-level Processing Facility the Astronomical Institute of the University of Bern (AIUB) is responsible for the determination of the official Precise Science Orbit (PSO) product. Kinematic GOCE orbits are part of this product and are used by several institutions in- and outside the HPF for determining the low degrees of the Earth’s gravity field. AIUB GOCE GPS-only gravity field solutions using the Celestial Mechanics Approach and covering the Release 4 period as well as a more recent time interval at the lower orbit altitude are shown and discussed. Special attention is paid to the impact of systematic deficiencies in the kinematic orbits on the resulting gravity fields, e.g., related to the geomagnetic equator, and on possibilities to get rid of them.

Kinetic modeling of sodium in the lunar exosphere

Our knowledge about the lunar environment is based on a large volume of ground-based, remote, and in situ observations. These observations have been conducted at different times and sampled different pieces of such a complex system as the surface-bound exosphere of the Moon. Numerical modeling is the tool that can link results of these separate observations into a single picture. Being validated against previous measurements, models can be used for predictions and interpretation of future observations results. In this paper we present a kinetic model of the sodium exosphere of the Moon as well as results of its validation against a set of ground-based and remote observations. The unique characteristic of the model is that it takes the orbital motion of the Moon and the Earth into consideration and simulates both the exosphere as well as the sodium tail self-consistently. The extended computational domain covers the part of the Earth’s orbit at new Moon, which allows us to study the effect of Earth’s gravity on the lunar sodium tail. The model is fitted to a set of ground-based and remote observations by tuning sodium source rate as well as values of sticking, and accommodation coefficients. The best agreement of the model results with the observations is reached when all sodium atoms returning from the exosphere stick to the surface and the net sodium escape rate is about 5.3 × 1022 s−1.

Evidence of direct detection of interstellar deuterium in the local interstellar medium by IBEX

We report the first in situ measurements of neutral deuterium originating in the local interstellar medium (LISM) in Earth’s orbit. These measurements were performed with the IBEX-Lo camera on NASA’s interstellar boundary explorer (IBEX) satellite. All data from the spring observation periods of 2009 through 2011 have been analysed. In the three years of the IBEX mission time, the observation geometry and orbit allowed for a total observation time of 115.3 days for the LISM. However, the effects of the spinning spacecraft and the stepping through 8 energy channels mean that we are only observing the interstellar wind for a total time of 1.44 days, in which 2 counts for interstellar deuterium were collected. We report here a conservative number, because a possibility of systematic error or additional noise, though eliminated in our analysis to the best of our knowledge, only supports detection at a 1-sigma level. From these observations, we derive a ratio D/H = (5.8 ± 4.4) × 10-4 at 1 AU. After modelling the transport and loss of D and H from the termination shock to Earth’s orbit, we find that our result of D/HLISM = (1.6 ± 1.2) × 10-5 agrees with D/HLIC = (1.6 ± 0.4) × 10-5 for the local interstellar cloud. This weak interstellar signal is extracted from a strong terrestrial background signal consisting of sputter products from the sensor’s conversion surface. As reference, we accurately measure the terrestrial D/H ratio in these sputtered products and then discriminate this terrestrial background source. Because of the diminishing D and H signal at Earth’s orbit during the rising solar activity due to photoionisation losses and increased photon pressure, our result demonstrates that in situ measurements of interstellar deuterium in the inner heliosphere are only possible during solar minimum conditions.

The free escape continuum of diffuse ions upstream of the Earth's quasi-parallel bow shock

The Earth's bow shock is very efficient in accelerating ions out of the incident solar wind distribution to high energies (≈ 200 keV/e). Fluxes of energetic ions accelerated at the quasi-parallel bow shock, also known as diffuse ions, are best represented by exponential spectra in energy/charge, which require additional assumptions to be incorporated into these model spectra. One of these assumptions is a so-called "free escape boundary" along the interplanetary magnetic field into the upstream direction. Locations along the IBEX orbit are ideally suited for in situ measurements to investigate the existence of an upstream free escape boundary for bow shock accelerated ions. In this study we use 2 years of ion measurements from the background monitor on the IBEX spacecraft, supported by ACE solar wind observations. The IBEX Background Monitor is sensitive to protons > 14 keV, which includes the energy of the maximum flux for diffuse ions. With increasing distance from the bow shock along the interplanetary magnetic field, the count rates for diffuse ions stay constant for ions streaming away from the bow shock, while count rates for diffuse ions streaming toward the shock gradually decrease from a maximum value to ~1/e at distances of about 10 RE to 14 RE. These observations of a gradual decrease support the transition to a free escape continuum for ions of energy >14 keV at distances from 10 RE to 14 RE from the bow shock.

Solar wind reflection from the lunar surface: The view from far and near

The Moon appears bright in the sky as a source of energetic neutral atoms (ENAs). These ENAs have recently been imaged over a broad energy range both from near the lunar surface, by India's Chandrayaan-1 mission (CH-1), and from a much more distant Earth orbit by NASA's Interstellar Boundary Explorer (IBEX) satellite. Both sets of observations have indicated that a relatively large fraction of the solar wind is reflected from the Moon as energetic neutral hydrogen. CH-1's angular resolution over different viewing angles of the lunar surface has enabled measurement of the emission as a function of angle. IBEX in contrast views not just a swath but a whole quadrant of the Moon as effectively a single pixel, as it subtends even at the closest approach no more than a few degrees on the sky. Here we use the scattering function measured by CH-1 to model global lunar ENA emission and combine these with IBEX observations. The deduced global reflection is modestly larger (by a factor of 1.25) when the angular scattering function is included. This provides a slightly updated IBEX estimate of AH=0.11±0.06 for the global neutralized albedo, which is ˜25% larger than the previous values of 0.09±0.05, based on an assumed uniform scattering distribution.

An Optical Survey for Space Debris on Highly Eccentric MEO Orbits

Optical surveys for space debris in high-altitude orbits have been conducted since more than ten years. Originally these efforts concentrated mainly on the geostationary ring (GEO). Corresponding observation strategies, processing techniques and cataloguing approaches have been developed and successfully applied. The ESA GEO surveys, e.g., resulted in the detection of a significant population of small-size debris and later in the discovery of high area-to-mass ratio objects in GEO-like orbits. The observation scenarios were successively adapted to survey the geostationary transfer orbit (GTO) region; and recently surveys to search for debris in the medium Earth orbit (MEO) region of the global navigation satellite constellations were successfully conducted. Comparably less experience (both, in terms of practical observation and strategy definition) is available for eccentric orbits that (at least partly) are in the MEO region, in particular for the Molniya-type orbits. Several breakup events and deliberate fragmentations are known to have taken place in such orbits. Different survey and follow-up strategies for searching space debris objects in highly-eccentric MEO orbits, and to acquire orbits which are sufficiently accurate to catalogue such objects and to maintain their orbits over longer time spans were developed. Simulations were performed to compare the performance of different survey and cataloguing strategies. Eventually, optical observations were conducted in the framework of an ESA study using ESA’s Space Debris Telescope (ESASDT) the 1-m Zeiss telescope located at the Optical Ground Station (OGS) at the Teide Observatory at Tenerife, Spain. A first series of surveys of Molnjya-type orbits was performed between January and April 2013. During these four months survey observations were performed during nine nights. A basic survey consisted of observing a single geocentric field for 10 minutes. If a faint object was found, follow-up observations were performed during the same night to ensure a save rediscovery of the object during the next nights. Additional follow-up observations to maintain the orbits of these newly discovered faint objects were also acquired with AIUB ́s 1m ZIMLAT telescope in Zimmerwald, Switzerland. Eventually 195 basic surveys were performed during these nine nights corresponding to about 32.5 hours of observations. In total 24 uncorrelated faint objects were discovered and all known catalogue objects in the survey fields were detected. On average one uncorrelated object was found every 80 minutes. Some of these objects show a considerable brightness variation and have a high area-to-mass ratio as determined in the orbit estimation process.

Assessment of detectability of neutral interstellar deuterium by IBEX observations

Context. The abundance of deuterium in the interstellar gas in front of the Sun gives insight into the processes of filtration of neutral interstellar species through the heliospheric interface and potentially into the chemical evolution of the Galactic gas. Aims: We investigate the possibility of detection of neutral interstellar deuterium at 1 AU from the Sun by direct sampling by the Interstellar Boundary Explorer (IBEX). Methods: Using both previous and the most recent determinations of the flow parameters of neutral gas in the local interstellar cloud (LIC) and an observation-based model of solar radiation pressure and ionization in the heliosphere, we simulated the flux of neutral interstellar D at IBEX for the actual measurement conditions. We assessed the number of interstellar D atom counts expected during the first three years of IBEX operation. We also simulated the observations expected during an epoch of high solar activity. In addition, we calculated the expected counts of D atoms from the thin terrestrial water layer covering the IBEX-Lo conversion surface, sputtered by neutral interstellar He atoms. Results: Most D counts registered by IBEX-Lo are expected to come from the water layer, exceeding the interstellar signal by 2 orders of magnitude. However, the sputtering should stop once the Earth leaves the portion of orbit traversed by interstellar He atoms. We identify seasons during the year when mostly the genuine interstellar D atoms are expected in the signal. During the first 3 years of IBEX operations about 2 detectable interstellar D atoms are expected. This number is comparable to the expected number of sputtered D atoms registered during the same time intervals. Conclusions: The most favorable conditions for the detection occur during low solar activity, in an interval including March and April each year. The detection chances could be improved by extending the instrument duty cycle, say, by making observations in the special deuterium mode of IBEX-Lo.

Geocenter coordinates estimated from GNSS data as viewed by perturbation theory

Time series of geocenter coordinates were determined with data of two global navigation satellite systems (GNSSs), namely the U.S. GPS (Global Positioning System) and the Russian GLONASS (Global’naya Nawigatsionnaya Sputnikowaya Sistema). The data was recorded in the years 2008–2011 by a global network of 92 permanently observing GPS/GLONASS receivers. Two types of daily solutions were generated independently for each GNSS, one including the estimation of geocenter coordinates and one without these parameters. A fair agreement for GPS and GLONASS was found in the geocenter x- and y-coordinate series. Our tests, however, clearly reveal artifacts in the z-component determined with the GLONASS data. Large periodic excursions in the GLONASS geocenter z-coordinates of about 40 cm peak-to-peak are related to the maximum elevation angles of the Sun above/below the orbital planes of the satellite system and thus have a period of about 4 months (third of a year). A detailed analysis revealed that the artifacts are almost uniquely governed by the differences of the estimates of direct solar radiation pressure (SRP) in the two solution series (with and without geocenter estimation). A simple formula is derived, describing the relation between the geocenter z-coordinate and the corresponding parameter of the SRP. The effect can be explained by first-order perturbation theory of celestial mechanics. The theory also predicts a heavy impact on the GNSS-derived geocenter if once-per-revolution SRP parameters are estimated in the direction of the satellite’s solar panel axis. Specific experiments using GPS observations revealed that this is indeed the case. Although the main focus of this article is on GNSS, the theory developed is applicable to all satellite observing techniques. We applied the theory to satellite laser ranging (SLR) solutions using LAGEOS. It turns out that the correlation between geocenter and SRP parameters is not a critical issue for the SLR solutions. The reasons are threefold: The direct SRP is about a factor of 30–40 smaller for typical geodetic SLR satellites than for GNSS satellites, allowing it in most cases to not solve for SRP parameters (ruling out the correlation between these parameters and the geocenter coordinates); the orbital arc length of 7 days (which is typically used in SLR analysis) contains more than 50 revolutions of the LAGEOS satellites as compared to about two revolutions of GNSS satellites for the daily arcs used in GNSS analysis; the orbit geometry is not as critical for LAGEOS as for GNSS satellites, because the elevation angle of the Sun w.r.t. the orbital plane is usually significantly changing over 7 days.

Earth gravity field recovery using GPS, GLONASS, and SLR satellites

The time variable Earth’s gravity field provides the information about mass transport within the system Earth, i.e., the relationship of mass transport between atmosphere, oceans, and land hydrology. We recover the low-degree parameters of the time variable gravity field using microwave observations from GPS and GLONASS satellites and from SLR data to five geodetic satellites, namely LAGEOS-1/2, Starlette, Stella, and AJISAI. GPS satellites are particularly sensitive to specific coefficients of the Earth's gravity field, because of the deep 2:1 orbital resonance with Earth rotation (two revolutions of the GPS satellites per sidereal day). The resonant coefficients cause, among other, a “secular” drift (actually periodic variations of very long periods) of the semi-major axes of up to 5.3 m/day of GPS satellites. We processed 10 years of GPS and GLONASS data using the standard orbit models from the Center of Orbit Determination in Europe (CODE) with a simultaneous estimation of the Earth gravity field coefficients and other parameters, e.g., satellite orbit parameters, station coordinates, Earth rotation parameters, troposphere delays, etc. The weekly GNSS gravity solutions up to degree and order 4/4 are compared to the weekly SLR gravity field solutions. The SLR-derived geopotential coefficients are compared to monthly GRACE and CHAMP results.

SLR-derived terrestrial reference frame using the observations to LAGEOS-1/2, Starlette, Stella, and AJISAI

Currently, the contributions of Starlette, Stella, and AJISAI are not taken into account when defining the International Terrestrial Reference Frame (ITRF), despite the large amount of data collected in a long time-span. Consequently, the SLR-derived parameters and the SLR part of the ITRF are almost exclusively defined by LAGEOS-1 and LAGEOS-2. We investigate the potential of combining the observations to several SLR satellites with different orbital characteristics. Ten years of SLR data are homogeneously processed using the development version 5.3 of the Bernese GNSS Software. Special emphasis is put on orbit parameterization and the impact of LEO data on the estimation of the geocenter coordinates, Earth rotation parameters, Earth gravity field coefficients, and the station coordinates in one common adjustment procedure. We find that the parameters derived from the multi-satellite solutions are of better quality than those obtained in single satellite solutions or solutions based on the two LAGEOS satellites. A spectral analysis of the SLR network scale w.r.t. SLRF2008 shows that artifacts related to orbit perturbations in the LAGEOS-1/2 solutions, i.e., periods related to the draconitic years of the LAGEOS satellites, are greatly reduced in the combined solutions.

Observations of the northern seasonal polar cap on Mars: I. Spring sublimation activity and processes

Spring sublimation of the seasonal CO2 northern polar cap is a dynamic process in the current Mars climate. Phenomena include dark fans of dune material propelled out onto the seasonal ice layer, polygonal cracks in the seasonal ice, sand flow down slipfaces, and outbreaks of gas and sand around the dune margins. These phenomena are concentrated on the north polar erg that encircles the northern residual polar cap. The Mars Reconnaissance Orbiter has been in orbit for three Mars years, allowing us to observe three northern spring seasons. Activity is consistent with and well described by the Kieffer model of basal sublimation of the seasonal layer of ice applied originally in the southern hemisphere. Three typical weak spots have been identified on the dunes for escape of gas sublimed from the bottom of the seasonal ice layer: the crest of the dune, the interface of the dune with the interdune substrate, and through polygonal cracks in the ice. Pressurized gas flows through these vents and carries out material entrained from the dune. Furrows in the dunes channel gas to outbreak points and may be the northern equivalent of southern radially-organized channels ("araneiform" terrain), albeit not permanent. Properties of the seasonal CO2 ice layer are derived from timing of seasonal events such as when final sublimation occurs. Modification of dune morphology shows that landscape evolution is occurring on Mars today, driven by seasonal activity associated with sublimation of the seasonal CO2 polar cap.

Homogeneous reprocessing of GPS, GLONASS and SLR observations

The International GNSS Service (IGS) provides operational products for the GPS and GLONASS constellation. Homogeneously processed time series of parameters from the IGS are only available for GPS. Reprocessed GLONASS series are provided only by individual Analysis Centers (i. e. CODE and ESA), making it difficult to fully include the GLONASS system into a rigorous GNSS analysis. In view of the increasing number of active GLONASS satellites and a steadily growing number of GPS+GLONASS-tracking stations available over the past few years, Technische Universität Dresden, Technische Universität München, Universität Bern and Eidgenössische Technische Hochschule Zürich performed a combined reprocessing of GPS and GLONASS observations. Also, SLR observations to GPS and GLONASS are included in this reprocessing effort. Here, we show only SLR results from a GNSS orbit validation. In total, 18 years of data (1994–2011) have been processed from altogether 340 GNSS and 70 SLR stations. The use of GLONASS observations in addition to GPS has no impact on the estimated linear terrestrial reference frame parameters. However, daily station positions show an RMS reduction of 0.3 mm on average for the height component when additional GLONASS observations can be used for the time series determination. Analyzing satellite orbit overlaps, the rigorous combination of GPS and GLONASS neither improves nor degrades the GPS orbit precision. For GLONASS, however, the quality of the microwave-derived GLONASS orbits improves due to the combination. These findings are confirmed using independent SLR observations for a GNSS orbit validation. In comparison to previous studies, mean SLR biases for satellites GPS-35 and GPS-36 could be reduced in magnitude from −35 and −38 mm to −12 and −13 mm, respectively. Our results show that remaining SLR biases depend on the satellite type and the use of coated or uncoated retro-reflectors. For Earth rotation parameters, the increasing number of GLONASS satellites and tracking stations over the past few years leads to differences between GPS-only and GPS+GLONASS combined solutions which are most pronounced in the pole rate estimates with maximum 0.2 mas/day in magnitude. At the same time, the difference between GLONASS-only and combined solutions decreases. Derived GNSS orbits are used to estimate combined GPS+GLONASS satellite clocks, with first results presented in this paper. Phase observation residuals from a precise point positioning are at the level of 2 mm and particularly reveal poorly modeled yaw maneuver periods.