N. Fougere

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

67P/Churyumov-Gerasimenko (hereafter 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 × 1026 molecules s−1.

Three-dimensional direct simulation Monte-Carlo modeling of the coma of comet 67P/Churyumov-Gerasimenko observed by the VIRTIS and ROSINA instruments on board Rosetta

Since its rendezvous with comet 67P/Churyumov-Gerasimenko (67P), the Rosetta spacecraft has provided invaluable information contributing to our understanding of the cometary environment. On board, the VIRTIS and ROSINA instruments can both measure gas parameters in the rarefied cometary atmosphere, the so-called coma, and provide complementary results with remote sensing and in situ measurement techniques, respectively. The data from both ROSINA and VIRTIS instruments suggest that the source regions of H2O and CO2 are not uniformly distributed over the surface of the nucleus even after accounting for the changing solar illumination of the irregularly shaped rotating nucleus. The source regions of H2O and CO2 are also relatively different from one another. Aims. The use of a combination of a formal numerical data inversion method with a fully kinetic coma model is a way to correlate and interpret the information provided by these two instruments to fully understand the volatile environment and activity of comet 67P. Methods. In this work, the nonuniformity of the outgassing activity at the surface of the nucleus is described by spherical harmonics and constrained by ROSINA-DFMS data. This activity distribution is coupled with the local illumination to describe the inner boundary conditions of a 3D direct simulation Monte-Carlo (DSMC) approach using the Adaptive Mesh Particle Simulator (AMPS) code applied to the H2O and CO2 coma of comet 67P. Results. We obtain activity distribution of H2O and CO2 showing a dominant source of H2O in the Hapi region, while more CO2 is produced in the southern hemisphere. The resulting model outputs are analyzed and compared with VIRTIS-M/-H and ROSINA-DFMS measurements, showing much better agreement between model and data than a simpler model assuming a uniform surface activity. The evolution of the H2O and CO2 production rates with heliocentric distance are derived accurately from the coma model showing agreement between the observations from the different instruments and ground-based observations. Conclusions. We derive the activity distributions for H2O and CO2 at the surface of the nucleus described in spherical harmonics, which we couple to the local solar illumination to constitute the boundary conditions of our coma model. The model presented reproduces the coma observations made by the ROSINA and VIRTIS instruments on board the Rosetta spacecraft showing our understanding of the physics of 67P’s coma. This model can be used for further data analyses, such as dust modeling, in a future work.

First observations of H2O and CO2 vapor in comet 67P/Churyumov-Gerasimenko made by VIRTIS onboard Rosetta

Context. Outgassing from cometary nuclei involves complex surface and subsurface processes that need to be understood to investigate the composition of cometary ices from coma observations. Aims. We investigate the production of water, carbon dioxide, and carbon monoxide from the nucleus of comet 67P/Churyumov-Gerasimenko (67P). These species have different volatility and are key species of cometary ices. Methods. Using the high spectral-resolution channel of the Visible InfraRed Thermal Imaging Spectrometer (VIRTIS-H), we observed the ν3 vibrational bands of H2O and CO2 at 2.67 and 4.27 μm, respectively, from 24 November 2014 to 24 January 2015, when comet 67P was between 2.91 and 2.47 AU from the Sun. Observations were undertaken in limb-viewing geometry at distances from the surface of 0 to 1.5 km and with various line-of-sight (LOS) orientations in the body-fixed frame. A geometry tool was used to characterize the position of the LOS with respect to geomorphologic regions and the illumination properties of these regions. Results. The water production of 67P did not increase much from 2.9 to 2.5 AU. High water column densities are observed for LOS above the neck regions, suggesting they are the most productive in water vapor. While water production is weak in regions with low solar illumination, CO2 is outgassing from both illuminated and non-illuminated regions, which indicates that CO2 sublimates at a depth that is below the diurnal skin depth. The CO2/H2O column density ratio varies from 2 to 60%. For regions that are in sunlight, mean values between 2 and 7% are measured. The lower bound value is likely representative of the CO2/H2O production rate ratio from the neck regions. For carbon monoxide, we derive column density ratios CO/H2O < 1.9% and CO/CO2< 80%. An illumination-driven model, with a uniformly active surface releasing water at a mean rate of 8 × 1025 s-1, provides an overall agreement to VIRTIS-H data, although some mismatches show local surface inhomogeneities in water production. Rotational temperatures of 90–100 K are derived from H2O and CO2 averaged spectra.

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

UNUSUAL WATER PRODUCTION ACTIVITY OF COMET C/2012 S1 (ISON): OUTBURSTS AND CONTINUOUS FRAGMENTATION

The Solar Wind ANisotropies (SWAN) all-sky hydrogen Ly? camera on the SOlar and Heliospheric Observer (SOHO) satellite observed the hydrogen coma of comet C/2012 S1 (ISON) for most of the last month of its activity from 2013 October 24 to November 24, ending just 4?days before perihelion and its final disruption. The water production rate of the comet was determined from these observations. SOHO has been operating in a halo orbit around the Earth-Sun L1 Lagrange point since its launch in late 1995. Most water vapor produced by comets is ultimately photodissociated into two H atoms and one O atom producing a huge hydrogen coma that is routinely observed in the daily SWAN images in comets of sufficient brightness. Water production rates were calculated from 22 images over most of the last month of the pre-perihelion apparition. The water production rate increased very slowly on average from October 24.9 until November 12.9, staying between 1.8 and 3.4?? 1028?s?1, after which it increased dramatically, reaching 1.6 to 2?? 1030?s?1 from November 21.6 to 23.6. It was not detected after perihelion on December 3.7 when it should have been visible. We examine the active surface area necessary to explain the water production rate and its variation and are able to place constraints on the physical size of the original nucleus necessary to account for the large amount of activity from November 12.9 and until just before perihelion.

Evolution of CO2, CH4, and OCS abundances relative to H2O in the coma of comet 67P around perihelion from Rosetta/VIRTIS-H observations

Infrared observations of the coma of 67P/Churyumov–Gerasimenko were carried out from 2015 July to September, i.e. around perihelion (2015 August 13), with the high-resolution channel of the Visible and Infrared Thermal Imaging Spectrometer instrument onboard Rosetta. We present the analysis of fluorescence emission lines of H2O, CO2, 13CO2, OCS, and CH4 detected in limb sounding with the field of view at 2.7–5 km from the comet centre. Measurements are sampling outgassing from the illuminated Southern hemisphere, as revealed by H2O and CO2 raster maps, which show anisotropic distributions, aligned along the projected rotation axis. An abrupt increase of water production is observed 6 d after perihelion. In the meantime, CO2, CH4, and OCS abundances relative to water increased by a factor of 2 to reach mean values of 32, 0.47, and 0.18 per cent, respectively, averaging post-perihelion data. We interpret these changes as resulting from the erosion of volatile-poor surface layers. Sustained dust ablation due to the sublimation of water ice maintained volatile-rich layers near the surface until at least the end of the considered period, as expected for low thermal inertia surface layers. The large abundance measured for CO2 should be representative of the 67P nucleus original composition, and indicates that 67P is a CO2-rich comet. Comparison with abundance ratios measured in the Northern hemisphere shows that seasons play an important role in comet outgassing. The low CO2/H2O values measured above the illuminated Northern hemisphere are not original, but the result of the devolatilization of the uppermost layers.

Four-fluid MHD simulations of the plasma and neutral gas environment of comet 67P/Churyumov-Gerasimenko near perihelion

The neutral and plasma environment is critical in understanding the interaction of the solar wind and comet 67P/Churyumov-Gerasimenko (CG), the target of the European Space Agencys Rosetta mission. To serve this need and support the Rosetta mission, we have developed a 3-D four-fluid model, which is based on BATS-R-US (Block-Adaptive Tree Solarwind Roe-type Upwind Scheme) within SWMF (Space Weather Modeling Framework) that solves the governing multifluid MHD equations and the Euler equations for the neutral gas fluid. These equations describe the behavior and interactions of the cometary heavy ions, the solar wind protons, the electrons, and the neutrals. This model incorporates different mass loading processes, including photoionization and electron impact ionization, charge exchange, dissociative ion-electron recombination, and collisional interactions between different fluids. We simulated the plasma and neutral gas environment near perihelion in three different cases: an idealized comet with a spherical body and uniform neutral gas outflow, an idealized comet with a spherical body and illumination-driven neutral gas outflow, and comet CG with a realistic shape model and illumination-driven neutral gas outflow. We compared the results of the three cases and showed that the simulations with illumination-driven neutral gas outflow have magnetic reconnection, a magnetic pileup region and nucleus directed plasma flow inside the nightside reconnection region, which have not been reported in the literature.

IMAGING OBSERVATIONS of the HYDROGEN COMA of COMET 67P/CHURYUMOV-GERASIMENKO in 2015 SEPTEMBER by the PROCYON/LAICA

The water production rate of a comet is one of the fundamental parameters necessary to understand cometary activity when a comet approaches the Sun within 2.5 au, because water is the most abundant icy material in the cometary nucleus. Wide-field imaging observations of the hydrogen Lyα emission in comet 67P/Churyumov–Gerasimenko were performed by the Lyman Alpha Imaging Camera (LAICA) on board the 50 kg class micro spacecraft, the Proximate Object Close Flyby with Optical Navigation (PROCYON), on UT 2015 September 7.40, 12.37, and 13.17 (corresponding to 25.31, 30.28, and 31.08 days after the perihelion passage of the comet, respectively). We derive the water production rates, , of the comet from Lyα images of the comet by using a 2D axi-symmetric Direct Simulation Monte-Carlo model of the atomic hydrogen coma; (1.46 ± 0.47) × 1028, (1.24 ± 0.40) × 1028, and (1.30 ± 0.42) × 1028 molecules s−1 on 7.40, 12.37, and 13.17 September, respectively. These values are comparable to the values from in situ measurements by the Rosetta instruments in the 2015 apparition and the ground-based and space observations during the past apparitions. The comet did not show significant secular change in average water production rates just after the perihelion passage for the apparitions from 1982 to 2015. We emphasize that the measurements of absolute based on the wide field of view (e.g., by the LAICA/PROCYON) are so important to judge the soundness of the coma models used to infer based on in situ measurements by spacecraft, like the Rosetta.

Ultraviolet observations of the hydrogen coma of comet C/2013 A1 (Siding Spring) by MAVEN/IUVS

We used the Imaging Ultraviolet Spectrograph (IUVS) aboard the Mars Atmosphere and Volatile EvolutioN (MAVEN) orbiting spacecraft to construct images of the hydrogen coma of comet C/2013 A1 (Siding Spring) days before its close encounter with Mars. We obtain a water production rate of 1.1 ± 0.5 × 1028 molecules/s and determine the total impacting fluence of atoms and molecules corresponding to the photodissociation of water and its daughter species to be 2.4 ± 1.2 × 104 kg. We use these observations to confirm predictions that the mass of delivered hydrogen is comparable to the existing reservoir above 150 km. Furthermore, we reconcile disparity between observations and predictions about the detectability of the hydrogen perturbation and thermospheric response.

The heterogeneous coma of comet 67P/Churyumov-Gerasimenko as seen by ROSINA: H 2 O, CO 2, and CO from September 2014 to February 2016

Context. The ESA Rosetta mission has been investigating the environment of comet 67P/Churyumov-Gerasimenko (67P) since August 2014. Among the experiments on board the spacecraft, the ROSINA experiment (Rosetta Orbiter Spectrometer for Ion and Neutral Analysis) includes two mass spectrometers to analyse the composition of neutrals and ions and a COmet Pressure Sensor (COPS) to monitor the density and velocity of neutrals in the coma. Aims. We study heterogeneities in the coma during three periods starting in October 2014 (summer in the northern hemisphere) and ending in February 2016 (end of winter in the northern hemisphere). We provide a detailed description of the main volatiles dynamics (H 2 O, CO 2 , CO) and their abundance ratios. Methods. We analysed and compared the data of the Reflectron-type Time-Of-Flight (RTOF) mass spectrometer with data from both the Double Focusing Mass Spectrometer (DFMS) and COPS during the comet escort phase. This comparison has demonstrated that the observations performed with each ROSINA sensor are indeed consistent. Furthermore, we used a Direct Simulation Monte Carlo (DSMC) model to compare modelled densitites with in situ detections. Results. Our analysis shows how the active regions of the main volatiles evolve with the seasons with a variability mostly driven by the illumination conditions; this is the case except for an unexpected dichotomy suggesting the presence of a dust layer containing water deposited in the northern hemisphere during previous perihelions hiding the presence of CO 2 . The influence of various parameters is investigated in detail: distance to the comet, heliocentric distance, longitude and latitude of sub-satellite point, local time, and phase angle.