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.

Narrow dust jets in a diffuse gas coma: a natural product of small active regions on comets

Comets often display narrow dust jets but more diffuse gas comae when their eccentric orbits bring them into the inner solar system and sunlight sublimates the ice on the nucleus. Comets are also understood to have one or more active areas covering only a fraction of the total surface active with sublimating volatile ices. Calculations of the gas and dust distribution from a small active area on a comet’s nucleus show that as the gas moves out radially into the vacuum of space it expands tangentially, filling much of the hemisphere centered on the active region. The dust dragged by the gas remains more concentrated over the active area. This explains some puzzling appearances of comets having collimated dust jets but more diffuse gaseous atmospheres. Our test case is 67P/Churyumov–Gerasimenko, the Rosetta mission target comet, whose activity is dominated by a single area covering only 4% of its surface.

Plasma environment of a weak comet – Predictions for Comet 67P/Churyumov–Gerasimenko from multifluid-MHD and Hybrid models

The interaction of a comet with the solar wind undergoes various stages as the comet’s activity varies along its orbit. For a comet like 67P/Churyumov–Gerasimenko, the target comet of ESA’s Rosetta mission, the various features include the formation of a Mach cone, the bow shock, and close to perihelion even a diamagnetic cavity. There are different approaches to simulate this complex interplay between the solar wind and the comet’s extended neutral gas coma which include magnetohydrodynamics (MHD) and hybrid-type models. The first treats the plasma as fluids (one fluid in basic single fluid MHD) and the latter treats the ions as individual particles under the influence of the local electric and magnetic fields. The electrons are treated as a charge-neutralizing fluid in both cases. Given the different approaches both models yield different results, in particular for a low production rate comet. In this paper we will show that these differences can be reduced when using a multifluid instead of a single-fluid MHD model and increase the resolution of the Hybrid model. We will show that some major features obtained with a hybrid type approach like the gyration of the cometary heavy ions and the formation of the Mach cone can be partially reproduced with the multifluid-type model.

The calculation of Afρ and mass loss rate for comets

Ab initio calculations of Afρ are presented using Mie scattering theory and a Direct Simulation Monte Carlo (DSMC) dust outflow model in support of the Rosetta mission and its target 67P/Churyumov-Gerasimenko (CG). These calculations are performed for particle sizes ranging from 0.010 μm to 1.0 cm. The present status of our knowledge of various differential particle size distributions is reviewed and a variety of particle size distributions is used to explore their effect on Afρ , and the dust mass production View the MathML sourcem˙. A new simple two parameter particle size distribution that curtails the effect of particles below 1 μm is developed. The contributions of all particle sizes are summed to get a resulting overall Afρ. The resultant Afρ could not easily be predicted a priori and turned out to be considerably more constraining regarding the mass loss rate than expected. It is found that a proper calculation of Afρ combined with a good Afρ measurement can constrain the dust/gas ratio in the coma of comets as well as other methods presently available. Phase curves of Afρ versus scattering angle are calculated and produce good agreement with observational data. The major conclusions of our calculations are: – The original definition of A in Afρ is problematical and Afρ should be: qsca(n,λ)×p(g)×f×ρqsca(n,λ)×p(g)×f×ρ. Nevertheless, we keep the present nomenclature of Afρ as a measured quantity for an ensemble of coma particles.– The ratio between Afρ and the dust mass loss rate View the MathML sourcem˙ is dominated by the particle size distribution. – For most particle size distributions presently in use, small particles in the range from 0.10 to 1.0 μm contribute a large fraction to Afρ. – Simplifying the calculation of Afρ by considering only large particles and approximating qsca does not represent a realistic model. Mie scattering theory or if necessary, more complex scattering calculations must be used. – For the commonly used particle size distribution, dn/da ∼ a−3.5 to a−4, there is a natural cut off in Afρ contribution for both small and large particles. – The scattering phase function must be taken into account for each particle size; otherwise the contribution of large particles can be over-estimated by a factor of 10. – Using an imaginary index of refraction of i = 0.10 does not produce sufficient backscattering to match observational data. – A mixture of dark particles with i ⩾ 0.10 and brighter silicate particles with i ⩽ 0.04 matches the observed phase curves quite well. – Using current observational constraints, we find the dust/gas mass-production ratio of CG at 1.3 AU is confined to a range of 0.03–0.5 with a reasonably likely value around 0.1.

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

ROSINA/DFMS and IES observations of 67P: Ion-neutral chemistry in the coma of a weakly outgassing comet

Context. The Rosetta encounter with comet 67P/Churyumov-Gerasimenko provides a unique opportunity for an in situ, up-close investigation of ion-neutral chemistry in the coma of a weakly outgassing comet far from the Sun. Aims. Observations of primary and secondary ions and modeling are used to investigate the role of ion-neutral chemistry within the thin coma. Methods. Observations from late October through mid-December 2014 show the continuous presence of the solar wind 30 km from the comet nucleus. These and other observations indicate that there is no contact surface and the solar wind has direct access to the nucleus. On several occasions during this time period, the Rosetta/ROSINA/Double Focusing Mass Spectrometer measured the low-energy ion composition in the coma. Organic volatiles and water group ions and their breakup products (masses 14 through 19), CO2+ (masses 28 and 44) and other mass peaks (at masses 26, 27, and possibly 30) were observed. Secondary ions include H3O+ and HCO+ (masses 19 and 29). These secondary ions indicate ion-neutral chemistry in the thin coma of the comet. A relatively simple model is constructed to account for the low H3O+/H2O+ and HCO+/CO+ ratios observed in a water dominated coma. Results from this simple model are compared with results from models that include a more detailed chemical reaction network. Results. At low outgassing rates, predictions from the simple model agree with observations and with results from more complex models that include much more chemistry. At higher outgassing rates, the ion-neutral chemistry is still limited and high HCO+/CO+ ratios are predicted and observed. However, at higher outgassing rates, the model predicts high H3O+/H2O+ ratios and the observed ratios are often low. These low ratios may be the result of the highly heterogeneous nature of the coma, where CO and CO2 number densities can exceed that of water.

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.

Composition-dependent outgassing of comet 67P/Churyumov-Gerasimenko from ROSINA/DFMS

Context. Early measurements of Rosetta’s target comet, 67P/Churyumov-Gerasimenko (67P), showed a strongly heterogeneous coma in H2O, CO, and CO2. Aims. The purpose of this work is to further investigate the coma heterogeneity of 67P, and to provide predictions for the near-perihelion outgassing profile based on the proposed explanations. Methods. Measurements of various minor volatile species by ROSINA/DFMS on board Rosetta are examined. The analysis focuses on the currently poorly illuminated winter (southern) hemisphere of 67P. Results. Coma heterogeneity is not limited to the major outgassing species. Minor species show better correlation with either H2O or CO2. The molecule CH4 shows a different diurnal pattern from all other analyzed species. Such features have implications for nucleus heterogeneity and thermal processing. Conclusions. Future analysis of additional volatiles and modeling the heterogeneity are required to better understand the observed coma profile.

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.

Sublimation of ice–tholins mixtures: A morphological and spectro-photometric study

Sublimation, the direct transition from solid to gas phase, is a process responsible for shaping and changing the reflectance properties of many Solar System surfaces. In this study, we have characterized the evolution of the structure/texture and of the visible and near-infrared (VIS–NIR) spectral reflectance of surfaces made of water ice mixed with analogues of complex extraterrestrial organic matter, named tholins, under low temperature (<-70° C) and pressure (10-⁵mbar) conditions. The experiments were carried out in the SCITEAS simulation setup recently built as part of the Laboratory for Outflow Studies of Sublimating Materials (LOSSy) at the University of Bern (Pommerol, A. et al. [2015a]. Planet. Space Sci. 109–110, 106–122). As the water ice sublimated, we observed in situ the formation of a sublimation lag deposit made of a water-free porous (>90% porosity) network of organic filaments on top of the ice. The temporal evolution of the tholins and water ice spectral features (reflectance at the absorption bands wavelengths, red slope, from 0.40 to 1.90lm) are analyzed throughout the sublimation of the samples. We studied how different mixtures of tholins with water (0.1 wt.% tholins as coating or inclusions within the water particles), and different ice particle sizes (4.5 ± 2.5 or 67 ± 31lm) influence the morphological and spectral evolutions of the samples. The sublimation of the ice below the mantle produces a gas flow responsible for the ejection of mm to cm-sized fragments of the deposit in outbursts-like events. The results show remarkable differences between these samples in term of mantle structure, speed of mantle building, rates and surface area of mantle ejections. These data provide useful references for interpreting remote-sensing observations of icy Solar System surfaces, in particular the activity of comet nuclei where sublimation of organic-rich ices and deposition of organic-dust particles likely play a major role. Consequently, the data presented here could be of high interest for the interpretation of Rosetta, and also New Horizons, observations.

Sensitivity and Fragmentation Calibration of the ROSINA Reflectron-type Time-Of-Flight Mass Spectrometer

The goal of this work was to increase the performance and to calibrate one of the ROSINA sensors, the Reflectron-type Time-Of-Flight mass spectrometer, currently flying aboard the ESA Rosetta spacecraft. Different optimization techniques were applied to both the lab and space models, and a static calibration was performed using different gas species expected to be detected in the vicinity of comet 67P/Churyumov-Gerasimenko. The database thus created was successfully applied to space data, giving consistent results with the other ROSINA sensors.

Sensitivity and fragmentation calibration of the ROSINA double focusing mass spectrometer

The goal of this work has been to calibrate sensitivities and fragmentation pattern of various molecules as well as further characterize the lab model of the ROSINA Double Focusing Mass Spectrometer (DFMS) on board ESA’s Rosetta spacecraft bound to comet 67P/Churyumov-Gerasimenko. The detailed calibration and characterization of the instrument is key to understand and interpret the results in the coma of the comet. A static calibration was performed for the following species: Ne, Ar, Kr, Xe, H2O, N2, CO2, CH4, C2H6, C3H8, C4H10, and C2H4. The purpose of the calibration was to obtain sensitivities for all detectors and emissions, the fragmentation behavior of the ion source and to show the capabilities to measure isotopic ratios at the comet. The calibration included the recording of different correction factors to evaluate the data, including a detailed investigation of the detector gain. The quality of the calibration that could be tested for different gas mixtures including the calibration of the density inside the ion source when calibration gas from the gas calibration unit is introduced. In conclusion the calibration shows that DFMS meets the design requirements and that DFMS will be able to measure the D/H at the comet and help shed more light on the puzzle about the origin of water on Earth.

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.