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.

Test Facility to Study Surface-Interaction Processes for Particle Detection in Space

The Imager for Low Energetic Neutral Atoms test facility at the University of Bern was developed to investigate, characterize, and quantify physical processes on surfaces that are used to ionize neutral atoms before their analysis in neutral particle-sensing instruments designed for space research. The facility has contributed valuable knowledge of the interaction of ions with surfaces (e.g., fraction of ions scattered from surfaces and angular scattering distribution) and employs a novel measurement principle for the determination of secondary electron emission yields as a function of energy, angle of incidence, particle species, and sample surface for low particle energies. Only because of this test facility it was possible to successfully apply surface-science processes for the new detection technique for low-energetic neutral particles with energies below about 1 keV used in space applications. All successfully flown spectrometers for the detection of low-energetic neutrals based on the particle–surface interaction process use surfaces evaluated, tested, and calibrated in this facility. Many instruments placed on different spacecraft (e.g., Imager for Magnetopause-to-Aurora Global Exploration, Chandrayaan-1, Interstellar Boundary Explorer, etc.) have successfully used this technique.

Lutetia׳s lineaments

The European Space Agency׳s Rosetta spacecraft flew by asteroid (21) Lutetia on July 10, 2010. Observations through the OSIRIS camera have revealed many geological features. Lineaments are identified on the entire observed surface of the asteroid. Many of these features are concentric around the North Pole Crater Cluster (NPCC). As observed on (433) Eros and (4) Vesta, this analysis of Lutetia assesses whether or not some of the lineaments could be created orthogonally to observed impact craters. The results indicate that the orientation of lineaments on Lutetia׳s surface could be explained by three impact craters: the Massilia and the NPCC craters observed in the northern hemisphere, and candidate crater Suspicio inferred to be in the southern hemisphere. The latter has not been observed during the Rosetta flyby. Of note, is that the inferred location of the Suspicio impact crater derived from lineaments matches locations where hydrated minerals have been detected from Earth-based observations in the southern hemisphere of Lutetia. Although the presence of these minerals has to be confirmed, this analysis shows that the topography may also have a significant contribution in the modification of the spectral shape and its interpretation. The cross-cutting relationships of craters with lineaments, or between lineaments themselves show that Massilia is the oldest of the three impact feature, the NPCC the youngest, and that the Suspicio impact crater is of intermediate age that is likely occurred closer in time to the Massilia event.

The spatial distribution of water in the inner coma of Comet 9P/Tempel 1: Comparison between models and observations

The near nucleus coma of Comet 9P/Tempel 1 has been simulated with the 3D Direct Simulation Monte Carlo (DSMC) code PDSC++ (Su, C.-C. [2013]. Parallel Direct Simulation Monte Carlo (DSMC) Methods for Modeling Rarefied Gas Dynamics. PhD Thesis, National Chiao Tung University, Taiwan) and the derived column densities have been compared to observations of the water vapour distribution found by using infrared imaging spectrometer on the Deep Impact spacecraft (Feaga, L.M., A’Hearn, M.F., Sunshine, J.M., Groussin, O., Farnham, T.L. [2007]. Icarus 191(2), 134–145. http://dx.doi.org/10.1016/j.icarus.2007.04.038). Modelled total production rates are also compared to various observations made at the time of the Deep Impact encounter. Three different models were tested. For all models, the shape model constructed from the Deep Impact observations by Thomas et al. (Thomas, P.C., Veverka, J., Belton, M.J.S., Hidy, A., A’Hearn, M.F., Farnham, T.L., et al. [2007]. Icarus, 187(1), 4–15. http://dx.doi.org/10.1016/j.icarus.2006.12.013) was used. Outgassing depending only on the cosine of the solar insolation angle on each shape model facet is shown to provide an unsatisfactory model. Models constructed on the basis of active areas suggested by Kossacki and Szutowicz (Kossacki, K., Szutowicz, S. [2008]. Icarus, 195(2), 705–724. http://dx.doi.org/10.1016/j.icarus.2007.12.014) are shown to be superior. The Kossacki and Szutowicz model, however, also shows deficits which we have sought to improve upon. For the best model we investigate the properties of the outflow.

Relating in situ gas measurements to the surface outgassing properties of cometary nuclei

The sensitivity of the gas flow field to changes in different initial conditions has been studied for the case of a highly simplified cometary nucleus model. The nucleus model simulated a homogeneously outgassing sphere with a more active ring around an axis of symmetry. The varied initial conditions were the number density of the homogeneous region, the surface temperature, and the composition of the flow (varying amounts of H2O and CO2) from the active ring. The sensitivity analysis was performed using the Polynomial Chaos Expansion (PCE) method. Direct Simulation Monte Carlo (DSMC) was used for the flow, thereby allowing strong deviations from local thermal equilibrium. The PCE approach can be used to produce a sensitivity analysis with only four runs per modified input parameter and allows one to study and quantify non-linear responses of measurable parameters to linear changes in the input over a wide range. Hence the PCE allows one to obtain a functional relationship between the flow field properties at every point in the inner coma and the input conditions. It is for example shown that the velocity and the temperature of the background gas are not simply linear functions of the initial number density at the source. As probably expected, the main influence on the resulting flow field parameter is the corresponding initial parameter (i.e. the initial number density determines the background number density, the temperature of the surface determines the flow field temperature, etc.). However, the velocity of the flow field is also influenced by the surface temperature while the number density is not sensitive to the surface temperature at all in our model set-up. Another example is the change in the composition of the flow over the active area. Such changes can be seen in the velocity but again not in the number density. Although this study uses only a simple test case, we suggest that the approach, when applied to a real case in 3D, should assist in identifying the sensitivity of gas parameters measured in situ by, for example, the Rosetta spacecraft to the surface boundary conditions and vice versa.

Imaging the Heliosphere using neutral Atoms from Solar Wind Energy down to 15 eV

We study the spatial and temporal distribution of hydrogen energetic neutral atoms (ENAs) from the heliosheath observed with the IBEX-Lo sensor of the Interstellar Boundary EXplorer (IBEX) from solar wind energies down to the lowest available energy (15 eV). All available IBEX-Lo data from 2009 January until 2013 June were included. The sky regions imaged when the spacecraft was outside of Earth's magnetosphere and when the Earth was moving toward the direction of observation offer a sufficient signal-to-noise ratio even at very low energies. We find that the ENA ribbon—a 20° wide region of high ENA intensities—is most prominent at solar wind energies whereas it fades at lower energies. The maximum emission in the ribbon is located near the poles for 2 keV and closer to the ecliptic plane for energies below 1 keV. This shift is an evidence that the ENA ribbon originates from the solar wind. Below 0.1 keV, the ribbon can no longer be identified against the globally distributed ENA signal. The ENA measurements in the downwind direction are affected by magnetospheric contamination below 0.5 keV, but a region of very low ENA intensities can be identified from 0.1 keV to 2 keV. The energy spectra of heliospheric ENAs follow a uniform power law down to 0.1 keV. Below this energy, they seem to become flatter, which is consistent with predictions. Due to the subtraction of local background, the ENA intensities measured with IBEX agree with the upper limit derived from Lyα observations.

Signal Processing for the Measurement of the Deuterium/Hydrogen Ratio in the Local Interstellar Medium

We report on a comprehensive signal processing procedure for very low signal levels for the measurement of neutral deuterium in the local interstellar medium from a spacecraft in Earth orbit. The deuterium measurements were performed with the IBEX-Lo camera on NASA’s Interstellar Boundary Explorer (IBEX) satellite. Our analysis technique for these data consists of creating a mass relation in three-dimensional time of flight space to accurately determine the position of the predicted D events, to precisely model the tail of the H events in the region where the H tail events are near the expected D events, and then to separate the H tail from the observations to extract the very faint D signal. This interstellar D signal, which is expected to be a few counts per year, 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. During the three years of the mission time when the deuterium signal was visible to IBEX, the observation geometry and orbit allowed for a total observation time of 115.3 days. Because of the spinning of the spacecraft and the stepping through eight energy channels the actual observing time of the interstellar wind was only 1.44 days. With the optimised data analysis we found three counts that could be attributed to interstellar deuterium. These results update our earlier work.

Doppler Orbit Determination of Deep Space Probes by the Bernese GNSS Software: First Results of the Combined Orbit Determination from DSN and Inter-Satellite Ka-Band Data from the Grail Mission

Navigation of deep space probes is most commonly operated using the spacecraft Doppler tracking technique. Orbital parameters are determined from a series of repeated measurements of the frequency shift of a microwave carrier over a given integration time. Currently, both ESA and NASA operate antennas at several sites around the world to ensure the tracking of deep space probes. Just a small number of software packages are nowadays used to process Doppler observations. The Astronomical Institute of the University of Bern (AIUB) has recently started the development of Doppler data processing capabilities within the Bernese GNSS Software. This software has been extensively used for Precise Orbit Determination of Earth orbiting satellites using GPS data collected by on-board receivers and for subsequent determination of the Earth gravity field. In this paper, we present the currently achieved status of the Doppler data modeling and orbit determination capabilities in the Bernese GNSS Software using GRAIL data. In particular we will focus on the implemented orbit determination procedure used for the combined analysis of Doppler and intersatellite Ka-band data. We show that even at this earlier stage of the development we can achieve an accuracy of few mHz on two-way S-band Doppler observation and of 2 µm/s on KBRR data from the GRAIL primary mission phase.

Optical Light Curve Observations to Determine Attitude States of Space Debris

The currently proposed space debris remediation measures include the active removal of large objects and “just in time” collision avoidance by deviating the objects using, e.g., ground-based lasers. Both techniques require precise knowledge of the attitude state and state changes of the target objects. In the former case, to devise methods to grapple the target by a tug spacecraft, in the latter, to precisely propagate the orbits of potential collision partners as disturbing forces like air drag and solar radiation pressure depend on the attitude of the objects. Non-resolving optical observations of the magnitude variations, so-called light curves, are a promising technique to determine rotation or tumbling rates and the orientations of the actual rotation axis of objects, as well as their temporal changes. The 1-meter telescope ZIMLAT of the Astronomical Institute of the University of Bern has been used to collect light curves of MEO and GEO objects for a considerable period of time. Recently, light curves of Low Earth Orbit (LEO) targets were acquired as well. We present different observation methods, including active tracking using a CCD subframe readout technique, and the use of a high-speed scientific CMOS camera. Technical challenges when tracking objects with poor orbit redictions, as well as different data reduction methods are addressed. Results from a survey of abandoned rocket upper stages in LEO, examples of abandoned payloads and observations of high area-to-mass ratio debris will be resented. Eventually, first results of the analysis of these light curves are provided.

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.

Development of a low energy ion source for ROSINA ion mode calibration

The European Rosetta mission on its way to comet 67P/Churyumov-Gerasimenko will remain for more than a year in the close vicinity (1 km) of the comet. The two ROSINA mass spectrometers on board Rosetta are designed to analyze the neutral and ionized volatile components of the cometary coma. However, the relative velocity between the comet and the spacecraft will be minimal and also the velocity of the outgassing particles is below 1km∕s. This combination leads to very low ion energies in the surrounding plasma of the comet, typically below 20eV. Additionally, the spacecraft may charge up to a few volts in this environment. In order to simulate such plasma and to calibrate the mass spectrometers, a source for ions with very low energies had to be developed for the use in the laboratory together with the different gases expected at the comet. In this paper we present the design of this ion source and we discuss the physical parameters of the ion beam like sensitivity, energy distribution, and beam shape. Finally, we show the first ion measurements that have been performed together with one of the two mass spectrometers.

Kinetic simulation of neutral∕ionized gas and electrically charged dust in the coma of comet 67P∕Churyumov-Gerasimenko

The cometary coma is a unique phenomenon in the solar system being a planetary atmosphere influenced by little or no gravity. As a comet approaches the sun, the water vapor with some fraction of other gases sublimate, generating a cloud of gas, ice and other refractory materials (rocky and organic dust) ejected from the surface of the nucleus. Sublimating gas molecules undergo frequent collisions and photochemical processes in the near‐nucleus region. Owing to its negligible gravity, comets produce a large and highly variable extensive dusty coma with a size much larger than the characteristic size of the cometary nucleus. The Rosetta spacecraft is en route to comet 67P/Churyumov‐Gerasimenko for a rendezvous, landing, and extensive orbital phase beginning in 2014. Both, interpretation of measurements and safety consideration of the spacecraft require modeling of the comet’s dusty gas environment. In this work we present results of a numerical study of multispecies gaseous and electrically charged dust environment of comet Chyuryumov‐Gerasimenko. Both, gas and dust phases of the coma are simulated kinetically. Photolytic reactions are taken into account. Parameters of the ambient plasma as well as the distribution of electric/magnetic fields are obtained from an MHD simulation [1] of the coma connected to the solar wind. Trajectories of ions and electrically charged dust grains are simulated by accounting for the Lorentz force and the nucleus gravity.

Comet 1P/Halley multifluid MHD model for the Giotto fly-by

The interaction of comets with the solar wind has been the focus of many studies including numerical modeling. We compare the results of our multifluid MHD simulation of comet 1P/Halley to data obtained during the flyby of the European Space Agency's Giotto spacecraft in 1986. The model solves the full set of MHD equations for the individual fluids representing the solar wind protons, the cometary light and heavy ions, and the electrons. The mass loading, charge-exchange, dissociative ion-electron recombination, and collisional interactions between the fluids are taken into account. The computational domain spans over several million kilometers, and the close vicinity of the comet is resolved to the details of the magnetic cavity. The model is validated by comparison to the corresponding Giotto observations obtained by the Ion Mass Spectrometer, the Neutral Mass Spectrometer, the Giotto magnetometer experiment, and the Johnstone Plasma Analyzer instrument. The model shows the formation of the bow shock, the ion pile-up, and the diamagnetic cavity and is able to reproduce the observed temperature differences between the pick-up ion populations and the solar wind protons. We give an overview of the global interaction of the comet with the solar wind and then show the effects of the Lorentz force interaction between the different plasma populations.

Molecular nitrogen in comet 67P/Churyumov-Gerasimenko indicates a low formation temperature

Molecular nitrogen (N2) is thought to have been the most abundant form of nitrogen in the protosolar nebula. It is the main N-bearing molecule in the atmospheres of Pluto and Triton and probably the main nitrogen reservoir from which the giant planets formed. Yet in comets, often considered the most primitive bodies in the solar system, N2 has not been detected. Here we report the direct in situ measurement of N2 in the Jupiter family comet 67P/Churyumov-Gerasimenko, made by the Rosetta Orbiter Spectrometer for Ion and Neutral Analysis mass spectrometer aboard the Rosetta spacecraft. A N2/CO ratio of Embedded Image (2σ standard deviation of the sampled mean) corresponds to depletion by a factor of ~25.4 ± 8.9 as compared to the protosolar value. This depletion suggests that cometary grains formed at low-temperature conditions below ~30 kelvin.

Time variability and heterogeneity in the coma of 67P/Churyumov-Gerasimenko

Comets contain the best-preserved material from the beginning of our planetary system. Their nuclei and comae composition reveal clues about physical and chemical conditions during the early solar system when comets formed. ROSINA (Rosetta Orbiter Spectrometer for Ion and Neutral Analysis) onboard the Rosetta spacecraft has measured the coma composition of comet 67P/Churyumov-Gerasimenko with well-sampled time resolution per rotation. Measurements were made over many comet rotation periods and a wide range of latitudes. These measurements show large fluctuations in composition in a heterogeneous coma that has diurnal and possibly seasonal variations in the major outgassing species: water, carbon monoxide, and carbon dioxide. These results indicate a complex coma-nucleus relationship where seasonal variations may be driven by temperature differences just below the comet surface.