The VIRTIS (Visible, Infrared and Thermal Imaging Spectrometer) instrument on board the Rosetta spacecraft has provided evidence of carbon-bearing compounds on the nucleus of the comet 67P/Churyumov-Gerasimenko. The very low reflectance of the nucleus (normal albedo of 0.060 ± 0.003 at 0.55 micrometers), the spectral slopes in visible and infrared ranges (5 to 25 and 1.5 to 5% kÅ(-1)), and the broad absorption feature in the 2.9-to-3.6-micrometer range present across the entire illuminated surface are compatible with opaque minerals associated with nonvolatile organic macromolecular materials: a complex mixture of various types of carbon-hydrogen and/or oxygen-hydrogen chemical groups, with little contribution of nitrogen-hydrogen groups. In active areas, the changes in spectral slope and absorption feature width may suggest small amounts of water-ice. However, no ice-rich patches are observed, indicating a generally dehydrated nature for the surface currently illuminated by the Sun.

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 (H 2 O, CO 2 , 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 H 2 O as it was before equinox, CO 2 and CO, which had a poor correlation with respect to H 2 O pre-equinox, also became well correlated with H 2 O 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 H 2 O and CO 2 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 H 2 O and CO 2 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 H 2 O and CO 2 coma of comet 67P. Results. We obtain activity distribution of H 2 O and CO 2 showing a dominant source of H 2 O in the Hapi region, while more CO 2 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 H 2 O and CO 2 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 H 2 O and CO 2 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.

Although water vapour is the main species observed in the coma of comet 67P/Churyumov-Gerasimenko and water is the major constituent of cometary nuclei, limited evidence for exposed water-ice regions on the surface of the nucleus has been found so far. The absence of large regions of exposed water ice seems a common finding on the surfaces of many of the comets observed so far. The nucleus of 67P/Churyumov-Gerasimenko appears to be fairly uniformly coated with dark, dehydrated, refractory and organic-rich material. Here we report the identification at infrared wavelengths of water ice on two debris falls in the Imhotep region of the nucleus. The ice has been exposed on the walls of elevated structures and at the base of the walls. A quantitative derivation of the abundance of ice in these regions indicates the presence of millimetre-sized pure water-ice grains, considerably larger than in all previous observations. Although micrometre-sized water-ice grains are the usual result of vapour recondensation in ice-free layers, the occurrence of millimetre-sized grains of pure ice as observed in the Imhotep debris falls is best explained by grain growth by vapour diffusion in ice-rich layers, or by sintering. As a consequence of these processes, the nucleus can develop an extended and complex coating in which the outer dehydrated crust is superimposed on layers enriched in water ice. The stratigraphy observed on 67P/Churyumov-Gerasimenko is therefore the result of evolutionary processes affecting the uppermost metres of the nucleus and does not necessarily require a global layering to have occurred at the time of the comet's formation.