Observations provided from RPC-ICA in combination with the data from RPC-MAG and ROSINA-COPS show that many aspects of the time variability of the detected ions is correlated with the magnetic field or -- to a smaller extent -- with neutral atmosphere density. We also show that not all changes in the cometary ion data reflect the nature of the plasma dynamics, but are a consequence of the instrumental limitations. The main outcome of the article in Appendix 1 is that the cometary ions can be divided into two populations with distinct characteristics. One population we termed the convecting population, is accelerated to higher energies through the interaction with the solar wind. The other population we termed the expanding population is moving radially away from the nucleus in the terminator plane. Both populations exhibit a significant anti-sunward component.In addition we present in this thesis a case with observations day-side of the terminator plane. There we show how the expanding population has a sunward component, consistent with initial radial expansion of the ions from the nucleus which gradually turn into an anti-sunward flow which is then observed in the terminator plane.

We analyse the Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA) - the Double Focusing Mass Spectrometer data between 2014 August and 2016 February to examine the effect of seasonal variations on the four major species within the coma of 67P/Churyumov-Gerasimenko (H2O, CO2, CO, and O-2), resulting from the tilt in the orientation of the comet's spin axis. Using a numerical data inversion, we derive the non-uniform activity distribution at the surface of the nucleus for these species, suggesting that the activity distribution at the surface of the nucleus has not significantly been changed and that the differences observed in the coma are solely due to the variations in illumination conditions. A three-dimensional Direct Simulation Monte Carlo model is applied where the boundary conditions are computed with a coupling of the surface activity distributions and the local illumination. The model is able to reproduce the evolution of the densities observed by ROSINA including the changes happening at equinox. While O-2 stays correlated with H2O as it was before equinox, CO2 and CO, which had a poor correlation with respect to H2O pre-equinox, also became well correlated with H2O post-equinox. The integration of the densities from the model along the line of sight results in column densities directly comparable to the VIRTIS-H observations. Also, the evolution of the volatiles' production rates is derived from the coma model showing a steepening in the production rate curves after equinox. The model/data comparison suggests that the seasonal effects result in the Northern hemisphere of 67P's nucleus being more processed with a layered structure while the Southern hemisphere constantly exposes new material.

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