Comparison of GOCE-GPS gravity fields derived by different approaches

Several techniques have been proposed to exploit GNSS-derived kinematic orbit information for the determination of long-wavelength gravity field features. These methods include the (i) celestial mechanics approach, (ii) short-arc approach, (iii) point-wise acceleration approach, (iv) averaged acceleration approach, and (v) energy balance approach. Although there is a general consensus that—except for energy balance—these methods theoretically provide equivalent results, real data gravity field solutions from kinematic orbit analysis have never been evaluated against each other within a consistent data processing environment. This contribution strives to close this gap. Target consistency criteria for our study are the input data sets, period of investigation, spherical harmonic resolution, a priori gravity field information, etc. We compare GOCE gravity field estimates based on the aforementioned approaches as computed at the Graz University of Technology, the University of Bern, the University of Stuttgart/Austrian Academy of Sciences, and by RHEA Systems for the European Space Agency. The involved research groups complied with most of the consistency criterions. Deviations only occur where technical unfeasibility exists. Performance measures include formal errors, differences with respect to a state-of-the-art GRACE gravity field, (cumulative) geoid height differences, and SLR residuals from precise orbit determination of geodetic satellites. We found that for the approaches (i) to (iv), the cumulative geoid height differences at spherical harmonic degree 100 differ by only ≈10 % ; in the absence of the polar data gap, SLR residuals agree by ≈96 % . From our investigations, we conclude that real data analysis results are in agreement with the theoretical considerations concerning the (relative) performance of the different approaches.

Optimizing the aquatic toxicity assessment under REACH through an integrated testing strategy (ITS).

To satisfy REACH requirements a high number of data on chemical of interest should be supplied to the European Chemicals Agency. To organize the various kinds of information and help the registrants to choose the best strategy to obtain the needed information limiting at the minimum the use of animal testing, integrated testing strategies (ITSs) schemes can be used. The present work deals with regulatory data requirements for assessing the hazards of chemicals to the aquatic pelagic environment. We present an ITS scheme for organizing and using the complex existing data available for aquatic toxicity assessment. An ITS to optimize the choice of the correct prediction strategy for aquatic pelagic toxicity is described. All existing information (like physico-chemical information), and all the alternative methods (like in silico, in vitro or the acute-to-chronic ratio) are considered. Moreover the weight of evidence approach to combine the available data is included.

GOCE: precise orbit determination for the entire mission

The Gravity field and steady-state Ocean Circulation Explorer (GOCE) was the first Earth explorer core mission of the European Space Agency. It was launched on March 17, 2009 into a Sun-synchronous dusk-dawn orbit and re-entered into the Earth’s atmosphere on November 11, 2013. The satellite altitude was between 255 and 225 km for the measurement phases. The European GOCE Gravity consortium is responsible for the Level 1b to Level 2 data processing in the frame of the GOCE High-level processing facility (HPF). The Precise Science Orbit (PSO) is one Level 2 product, which was produced under the responsibility of the Astronomical Institute of the University of Bern within the HPF. This PSO product has been continuously delivered during the entire mission. Regular checks guaranteed a high consistency and quality of the orbits. A correlation between solar activity, GPS data availability and quality of the orbits was found. The accuracy of the kinematic orbit primarily suffers from this. Improvements in modeling the range corrections at the retro-reflector array for the SLR measurements were made and implemented in the independent SLR validation for the GOCE PSO products. The satellite laser ranging (SLR) validation finally states an orbit accuracy of 2.42 cm for the kinematic and 1.84 cm for the reduced-dynamic orbits over the entire mission. The common-mode accelerations from the GOCE gradiometer were not used for the official PSO product, but in addition to the operational HPF work a study was performed to investigate to which extent common-mode accelerations improve the reduced-dynamic orbit determination results. The accelerometer data may be used to derive realistic constraints for the empirical accelerations estimated for the reduced-dynamic orbit determination, which already improves the orbit quality. On top of that the accelerometer data may further improve the orbit quality if realistic constraints and state-of-the-art background models such as gravity field and ocean tide models are used for the reduced-dynamic orbit determination.

Kelvin-Helmholtz instabilities at the magnetic cavity boundary of comet 67P/Churyumov-Gerasimenko

We investigate the plasma environment of comet 67P/Churyumov-Gerasimenko, the target of the European Space Agency's Rosetta mission. Rosetta will rendezvous with the comet in 2014 at almost 3.5 AU and follow it all the way to and past perihelion at 1.3 AU. During its journey towards the inner solar system the comet's environment will significantly change. The interaction of the solar wind with a well developed neutral coma leads to the formation of an upstream bow shock and, closer to the comet, the inner shock separating the solar wind, with cometary pick-up ions mass-loaded, from the inner cometary ions which are dragged outward through abundant collisions and charge exchange with the expanding neutral gas. As a consequence the interplanetary magnetic field is prevented from penetrating the innermost region of the comet, the so-called magnetic cavity. We use our magnetohydrodynamics model BATSRUS (Block-Adaptive-Tree-Solarwind-Roe-Upwind-Scheme) to simulate the solar wind - comet interaction. The model includes photoionization, ion-electron recombination, and charge exchange. Under certain conditions our model predicts an unstable plasma flow at the inner shock. We show that the plasma shear flow around the magnetic cavity can lead to Kelvin-Helmholtz instabilities. We investigate the onset of this phenomenon with change of heliocentric distance and furthermore show that a previously stable magnetic cavity boundary can become unstable when the neutral gas is predominately released from the dayside of the comet.

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.

Comparison of 3D kinetic and hydrodynamic models to ROSINA-COPS measurements of the neutral coma of 67P/Churyumov-Gerasimenko

67P/Churyumov-Gerasimenko (67P) is a Jupiter-family comet and the object of investigation of the European Space Agency mission Rosetta. This report presents the first full 3D simulation results of 67P’s neutral gas coma. In this study we include results from a direct simulation Monte Carlo method, a hydrodynamic code, and a purely geometric calculation which computes the total illuminated surface area on the nucleus. All models include the triangulated 3D shape model of 67P as well as realistic illumination and shadowing conditions. The basic concept is the assumption that these illumination conditions on the nucleus are the main driver for the gas activity of the comet. As a consequence, the total production rate of 67P varies as a function of solar insolation. The best agreement between the model and the data is achieved when gas fluxes on the night side are in the range of 7% to 10% of the maximum flux, accounting for contributions from the most volatile components. To validate the output of our numerical simulations we compare the results of all three models to in situ gas number density measurements from the ROSINA COPS instrument. We are able to reproduce the overall features of these local neutral number density measurements of ROSINA COPS for the time period between early August 2014 and January 1 2015 with all three models. Some details in the measurements are not reproduced and warrant further investigation and refinement of the models. However, the overall assumption that illumination conditions on the nucleus are at least an important driver of the gas activity is validated by the models. According to our simulation results we find the total production rate of 67P to be constant between August and November 2014 with a value of about 1 × 10²⁶ molecules s⁻¹.

Molecular oxygen in oort cloud comet 1P/Halley

Recently, the ROSINA mass spectrometer suite on board the European Space Agency's Rosetta spacecraft discovered an abundant amount of molecular oxygen, O2, in the coma of Jupiter family comet 67P/Churyumov-Gerasimenko of O2/H2O = 3.80 ± 0.85%. It could be shown that O2 is indeed a parent species and that the derived abundances point to a primordial origin. Crucial questions are whether the O2 abundance is peculiar to comet 67P/Churyumov-Gerasimenko or Jupiter family comets in general, and also whether Oort cloud comets such as comet 1P/Halley contain similar amounts of molecular oxygen. We investigated mass spectra obtained by the Neutral Mass Spectrometer instrument during the flyby by the European Space Agency's Giotto probe of comet 1P/Halley. Our investigation indicates that a production rate of O2 of 3.7 ± 1.7% with respect to water is indeed compatible with the obtained Halley data and therefore that O2 might be a rather common and abundant parent species.

The capabilities of ROSINA/DFMS to measure argon isotopes at comet 67P/Churyumov–Gerasimenko

Little is known about the noble gas abundances in comets. These highly volatile atoms are possible tracers of the history of cometary matter including the thermal evolution. They can help quantify the contribution of cometary impacts to terrestrial oceans and help elucidate on the formation history of comets and their role in the formation and evolution of planetary atmospheres. This paper focuses on argon and the capabilities to measure this noble gas with in situ mass spectrometry at comet 67P/Churyumov–Gerasimenko, the target of the European Space Agency׳s spacecraft Rosetta. Argon may have been detected by remote sensing in a single Oort cloud comet but to date nothing is known about the isotopic abundances of argon in comets. Furthermore, no detection of argon in a Jupiter-family comet has been reported. Comet 67P/Churyumov–Gerasimenko belongs to the group of Jupiter-family comets and originates most likely in the Kuiper belt. Onboard Rosetta is ROSINA/DFMS (Rosetta Orbiter Spectrometer for Ion and Neutral Analysis/Double Focusing Mass Spectrometer). DFMS has unprecedented mass resolution and high sensitivity and is designed to measure isotopic ratios including argon (Balsiger et al., 2007). Argon measurements using the DFMS lab model (identical to the flight model) demonstrate this capability. At very least, this mass spectrometer has the resolution and sensitivity to reduce the upper limit on the argon outgassing rate relative to water by more than three orders of magnitude (for 38Ar). Most likely, ROSINA/DFMS will provide the first detection of argon in a Jupiter-family comet together with the first determination of the ³⁶Ar/³⁸Ar ratio at a comet.

Sulfur Isotopic Ratios at 67P/Churyumov-Gerasimenko and Characterization of ROSINA-DFMS FM & FS

Comets are thought to be the most pristine bodies present in the Solar System. In consequence of spending the majority of their existence beyond 30 AU, their composition can give insights on the physical and chemical conditions during their formation. Since August 2014 the European Space Agency spacecraft Rosetta accompanies the Jupiter family comet 67P/Churyumov-Gerasimenko on its way to perihelion and beyond. In this study the isotope fractionation of 34S are reported in H2S, OCS, SO2, S2, and CS2 at 67P. In addition for the first time the isotope fractionation for 33S is presented for cometary volatiles. The ratio 32S/33S is given for H2S, SO2 and a tentative value is given for CS2. With a mean value of -50 ± 22‰ and -306 ± 31‰ for δ34S and δ33S respectively, H2S shows a significant depletion in both 34S and 33S. For SO2 the depletion is less distinct with δ34S and δ33S being -67 ± 40‰ and -130 ± 53‰, respectively. The strongest depletion is present for CS2 with -114 ± 21‰and -276 ± 55‰, respectively. For OCS and S2 only δ34S could be determined which is -252 ± 77‰ and -357 ± 145‰, respectively. A comparison with sulfur isotopic ratios measured in SiC grains revealed that both SiC grains and the five volatile species have similar sulfur isotopic ratios. However, it is beyond the scope of this work to investigate the possibility of a link between SiC grains and cometary ices. Nevertheless, mass-dependent or mass-independent fractionation due to photo dissociation can be ruled out as sole cause of the seen depletion of 33S and 34S. Furthermore, an upper limit of (9.64 ± 0.19)·10.4 for D/H in HDS has been determined. This value is about a factor two higher than D/H in H2O for the same comet reported by (Altwegg et al., 2015). Besides the investigation concerning isotopic ratios of sulfur bearing species in this work the calibration and characterization of ROSINA/DFMS has been continued. Here it is reported about the deviation of the mass scale for MCP/LEDA low resolution spectra and the calibration measurements performed in the laboratory. Furthermore the outcome of the attempt to describe the sensitivity of DFMS with an empirical function will be discussed. The last part of the characterization of DFMS is dedicated to determine the so-called individual pixel gain for the laboratory and the flight model. Moreover, correlation between the depletion’s manifestation of the MCP with respect to the applied voltages has been investigated for both models. It has been found that further measurements are needed to understand the manifestation of depletion at the laboratory model. For the model on board of Rosetta it could be shown that most of the present feature are due to the usage of the MCP and suggestions have been made in order to answer the remaining question considering the depletion of the MCP.

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.

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.

Dust measurements in the coma of comet 67P/Churyumov-Gerasimenko inbound to the Sun

Critical measurements for understanding accretion and the dust/gas ratio in the solar nebula, where planets were forming 4.5 billion years ago, are being obtained by the GIADA (Grain Impact Analyser and Dust Accumulator) experiment on the European Space Agency's Rosetta spacecraft orbiting comet 67P/Churyumov-Gerasimenko. Between 3.6 and 3.4 astronomical units inbound, GIADA and OSIRIS (Optical, Spectroscopic, and Infrared Remote Imaging System) detected 35 outflowing grains of mass 10(-10) to 10(-7) kilograms, and 48 grains of mass 10(-5) to 10(-2) kilograms, respectively. Combined with gas data from the MIRO (Microwave Instrument for the Rosetta Orbiter) and ROSINA (Rosetta Orbiter Spectrometer for Ion and Neutral Analysis) instruments, we find a dust/gas mass ratio of 4 +/- 2 averaged over the sunlit nucleus surface. A cloud of larger grains also encircles the nucleus in bound orbits from the previous perihelion. The largest orbiting clumps are meter-sized, confirming the dust/gas ratio of 3 inferred at perihelion from models of dust comae and trails.

The morphological diversity of comet 67P/Churyumov-Gerasimenko

Images of comet 67P/Churyumov-Gerasimenko acquired by the OSIRIS (Optical, Spectroscopic and Infrared Remote Imaging System) imaging system onboard the European Space Agency's Rosetta spacecraft at scales of better than 0.8 meter per pixel show a wide variety of different structures and textures. The data show the importance of airfall, surface dust transport, mass wasting, and insolation weathering for cometary surface evolution, and they offer some support for subsurface fluidization models and mass loss through the ejection of large chunks of material.

Detection efficiency of microchannel plates for e(-) and pi(-) in the momentum range from 17.5 to 345 MeV/c

High-energy e(-) and pi(-) were measured by the multichannel plate (MCP) detector at the PiM1 beam line of the High Intensity Proton Accelerator Facilities located at the Paul Scherrer Institute, Villigen, Switzerland. The measurements provide the absolute detection efficiencies for these particles: 5.8% +/- 0.5% for electrons in the beam momenta range 17.5-300 MeV/c and 6.0% +/- 1.3% for pions in the beam momenta range 172-345 MeV/c. The pulse height distribution determined from the measurements is close to an exponential function with negative exponent, indicating that the particles penetrated the MCP material before producing the signal somewhere inside the channel. Low charge extraction and nominal gains of the MCP detector observed in this study are consistent with the proposed mechanism of the signal formation by penetrating radiation. A very similar MCP ion detector will be used in the Neutral Ion Mass (NIM) spectrometer designed for the JUICE mission of European Space Agency (ESA) to the Jupiter system, to perform measurements of the chemical composition of the Galilean moon exospheres. The detection efficiency for penetrating radiation determined in the present studies is important for the optimisation of the radiation shielding of the NIM detector against the high-rate and high-energy electrons trapped in Jupiter's magnetic field. Furthermore, the current studies indicate that MCP detectors can be useful to measure high-energy particle beams at high temporal resolution. (C) 2015 AIP Publishing LLC.