The coma of comet 67P/Churyumov-Gerasimenko has been probed by the Rosetta spacecraft and shows a variety of different molecules. The ROSINA COmet Pressure Sensor and the Double Focusing Mass Spectrometer provide in-situ densities for many volatile compounds including the 14 gas species H2O, CO2, CO, H2S, O2, C2H6, CH3OH, H2CO, CH4, NH3, HCN, C2H5OH, OCS, and CS2. We fit the observed densities during the entire comet mission between August 2014 and September 2016 to an inverse coma model. We retrieve surface emissions on a cometary shape with 3996 triangular elements for 50 separated time intervals. For each gas we derive systematic error bounds and report the temporal evolution of the production, peak production, and the time-integrated total production. We discuss the production for the two lobes of the nucleus and for the northern and southern hemispheres. Moreover we provide a comparison of the gas production with the seasonal illumination.

In 2004, the Rosetta spacecraft was sent to comet 67P/Churyumov-Gerasimenko for the first ever long-term investigation of a comet. After its arrival in 2014, the spacecraft spent more than two years in immediate proximity to the comet. During these two years, the ROSINA Double Focusing Mass Spectrometer (DFMS) onboard Rosetta discovered a coma with an unexpectedly complex chemical composition that included many oxygenated molecules. Determining the exact cometary composition is an essential first step to understanding of the organic rich chemistry in star forming regions and protoplanetary disks that are ultimately conserved in cometary ices. In this study a joint approach of laboratory calibration and space data analysis was used to perform a detailed identification and quantification of CHO-compounds in the coma of 67P/Churyumov-Gerasimenko. The goal was to derive the CHO-compound abundances relative to water for masses up to 100 u. For this study, the May 2015 post-equinox period represent the best bulk abundances of comet 67P/Churyumov-Gerasimenko. A wide variety of CHO-compounds were discovered and their bulk abundances were derived. Finally, these results are compared to abundances of CHO-bearing molecules in other comets, obtained mostly from ground-based observations and modelling.

Cometary comae are generally depleted in nitrogen. The main carriers for volatile nitrogen in comets are NH3 and HCN. It is known that ammonia readily combines with many acids like e.g. HCN, HNCO, HCOOH, etc. encountered in the interstellar medium as well as in cometary ice to form ammonium salts (NH4+X-) at low temperatures. Ammonium salts, which can play a significant role in prebiotic chemistry, are hard to detect in space as they are unstable in the gas phase and their infrared signature is often hidden by thermal radiation or by e.g. OH in minerals. Here we report the presence of all possible sublimation products of five different ammonium salts at comet 67P/Churyumov-Gerasimenko measured by the ROSINA instrument on Rosetta. The relatively high sublimation temperatures of the salts leads to an apparent lack of volatile nitrogen in the coma. This then also explains the observed trend of higher NH3/H2O ratios with decreasing perihelion distances in comets.

To understand how Phosphorus-bearing molecules are formed in star-forming regions, we have analysed ALMA observations of PN and PO towards the massive star-forming region AFGL 5142, combined with a new analysis of the data of the comet 67P/Churyumov-Gerasimenko taken with the ROSINA instrument onboard Rosetta. The ALMA maps show that the emission of PN and PO arises from several spots associated with low-velocity gas with narrow linewidths in the cavity walls of a bipolar outflow. PO is more abundant than PN in most of the spots, with the PO/PN ratio increasing as a function of the distance to the protostar. Our data favor a formation scenario in which shocks sputter phosphorus from the surface of dust grains, and gas-phase photochemistry induced by UV photons from the protostar allows efficient formation of the two species in the cavity walls. Our analysis of the ROSINA data has revealed that PO is the main carrier of P in the comet, with PO/PN>10. Since comets may have delivered a significant amount of prebiotic material to the early Earth, this finding suggests that PO could contribute significantly to the phosphorus reservoir during the dawn of our planet. There is evidence that PO was already in the cometary ices prior to the birth of the Sun, so the chemical budget of the comet might be inherited from the natal environment of the Solar System, which is thought to be a stellar cluster including also massive stars.

The Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA) suite of instruments operated throughout the over two years of the Rosetta mission operations in the vicinity of comet 67P/Churyumov-Gerasimenko. It measured gas densities and composition throughout the comet's atmosphere, or coma. Here we present two-years' worth of measurements of the relative densities of the four major volatile species in the coma of the comet, H2O. CO2, CO and O2, by one of the ROSINA sub-systems called the Double Focusing Mass Spectrometer (DFMS). The absolute total gas densities were provided by the Comet Pressure Sensor (COPS), another ROSINA sub-system. DFMS is a very high mass resolution and high sensitivity mass spectrometer able to resolve at a tiny fraction of an atomic mass unit. We have analyzed the combined DFMS and COPS measurements using an inversion scheme based on spherical harmonics that solves for the distribution of potential surface activity of each species as the comet rotates, changing solar illumination, over short intervals and as the comet changes distance from the sun and orientation of its spin axis over long time intervals. We also use the surface boundary conditions derived from the inversion scheme to simulate the whole coma with our fully kinetic Direct Simulation Monte Carlo model and calculate the production rates of the four major species throughout the mission. We compare the derived production rates with revised remote sensing observations by the Visible and Infrared Thermal Imaging Spectrometer (VIRTIS) as well as with published observations from the Microwave Instrument for the Rosetta Orbiter (MIRO). Finally we use the variation of the surface production of the major species to calculate the total mass loss over the mission and, for different estimates of the dust/gas ratio, calculate the variation of surface loss over the nucleus.

Our modern day Solar System has $4.6\times10^9$ yrs of evolution behind it with just a few relics of its birth conditions remaining. Comets are thought to be some of the most pristine tracers of the initial ingredients that were combined to produce the Earth and the other planets. Other low-mass protostars may be analogous to our proto-Sun and hence, could be used to study the building blocks necessary to form Solar-like systems. This study tests this idea on the basis of new high sensitivity, high spatial resolution ALMA data on the protoplanetary disc-scales (~70 au) of IRAS 16293-2422 and the bulk composition of comet 67P/Churyumov-Gerasimenko, as determined for the first time with the unique in situ monitoring carried out by Rosetta. The comparative analysis of the observations from the Protostellar Interferometric Line Survey (PILS) and the measurements made with Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA) shows that the relative abundances of CHO-, N-, and S-bearing molecules correlate, with some scatter, between protostellar and cometary data. A tentative correlation is seen for the first time for P- and Cl-bearing compounds. The results imply that the volatile composition of cometesimals and planetesimals is partially inherited from the pre- and protostellar phases of evolution.

Comets are considered to be some of the most pristine and unprocessed solar system objects accessible to in-situ exploration. Investigating their molecular and elemental composition takes us on a journey back to the early period of our solar system and possibly even further. In this work, we deduce the bulk abundances of the major volatile species in comet 67P/Churyumov-Gerasimenko, the target of the European Space Agency's Rosetta mission. The basis are measurements obtained with the ROSINA instrument suite on board the Rosetta orbiter during a suitable period of high outgassing near perihelion. The results are combined with both gas and dust composition measurements published in the literature. This provides an integrated inventory of the major elements present in the nucleus of 67P/Churyumov-Gerasimenko. Similar to comet 1P/Halley, which was visited by ESA's Giotto spacecraft in 1986, comet 67P/Churyumov-Gerasimenko also shows near-solar abundances of oxygen and carbon, whereas hydrogen and nitrogen are depleted compared to solar. Still, the degree of devolatilization is lower than that of inner solar system objects, including meteorites and the Earth. This supports the idea that comets are among the most pristine objects in our solar system.

In situ research on cometary chemistry began when measurements from the Giotto mission at comet 1P/Halley revealed the presence of complex organics in the coma. New telescopes and space missions have provided detailed remote and in situ measurements of the composition of cometary volatiles. Recently, the Rosetta mission to comet 67P/Churyumov-Gerasimenko more than doubled the number of parent species and the number of isotopic ratios known in comets. 39 of the 66 parent species have also been detected in pre- and proto-stellar clouds, making the similarities of cometary ices with pre-stellar material very intriguing. Most isotopic ratios are non-solar. The variations in D/H in water between different comets indicate a large range in the protoplanetary disk where comets formed. All of these results point to a non-homogeneized protoplanetary disk where comets received their material. This diverse origin is in contrast to the Sun, who received its material from the bulk of the collapsing cloud. Detection of N2 and Ar in the coma of 67P and the presence of very volatile S2 point to low formation temperature of 20-30 K of the comet. The xenon isotopic ratios measured in 67P can explain the long standing question about the origin of the terrestrial atmospheric xenon. While we can exclude comets as being the source of the bulk terrestrial water due to their generally high D/H, the noble gases in the Earth atmosphere are compatible with a cometary delivery. This means that the amount of organics delivered by comets may be highly significant.

This white paper proposes that AMBITION, a Comet Nucleus Sample Return mission, be a cornerstone of ESA's Voyage 2050 programme. We summarise some of the most important questions still open in cometary science after the successes of the Rosetta mission, many of which require sample analysis using techniques that are only possible in laboratories on Earth. We then summarise measurements, instrumentation and mission scenarios that can address these questions, with a recommendation that ESA select an ambitious cryogenic sample return mission. Rendezvous missions to Main Belt comets and Centaurs are compelling cases for M-class missions, expanding our knowledge by exploring new classes of comets. AMBITION would engage a wide community, drawing expertise from a vast range of disciplines within planetary science and astrophysics. With AMBITION, Europe will continue its leadership in the exploration of the most primitive Solar System bodies.

The nucleus of the Jupiter-family comet 67P/Churyumov-Gerasimenko was discovered to be bi-lobate in shape when the European Space Agency spacecraft Rosetta first approached it in July 2014. The bi-lobate structure of the cometary nucleus has led to much discussion regarding the possible manner of its formation and on how the composition of each lobe might compare with that of the other. During its two-year-long mission from 2014 to 2016, Rosetta remained in close proximity to 67P/Churyumov-Gerasimenko, studying its coma and nucleus in situ. Based on lobe-specific measurements of HDO and H2O performed with the ROSINA DFMS mass spectrometer on board Rosetta, the Deuterium-to-Hydrogen ratios in water from the two lobes could be compared. No appreciable difference was observed, suggesting that both lobes formed in the same region and are homogeneous in their Deuterium-to-Hydrogen ratios.

The European Space Agency (ESA) spacecraft Rosetta accompanied the Jupiter-family comet (JFC) 67P / Churyumov-Gerasimenko for over two years along its trajectory through the inner solar system. Between 2014 and 2016, it performed almost continuous in-situ measurements of the comet's gaseous atmosphere in close proximity to its nucleus. In this study, the $^{16}$O/$^{18}$O ratio of H$_2$O in the coma of 67P / Churyumov-Gerasimenko, as measured by the ROSINA DFMS mass spectrometer on board Rosetta, was determined from the ratio of H${_2}^{16}$O / H${_2}^{18}$O and $^{16}$OH / $^{18}$OH. The value of 445 $\pm$ 35 represents an $\sim$ 11% enrichment of $^{18}$O compared with the terrestrial ratio of 498.7 $\pm$ 0.1. This cometary value is consistent with the comet containing primordial water, in accordance with leading self-shielding models. These models predict primordial water to be between 5% to 20% enriched in heavier oxygen isotopes compared to terrestrial water. Addendum: The $^{16}$O/$^{17}$O ratio of H$_2$O in the coma of 67P / Churyumov-Gerasimenko.

We reconstruct the temporal evolution of the source distribution for the four major gas species H2O, CO2, CO, and O2 on the surface of comet 67P/Churyumov-Gerasimenko during its 2015 apparition. The analysis applies an inverse coma model and fits to data between August 6th 2014 and September 5th 2016 measured with the Double Focusing Mass Spectrometer (DFMS) of the Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA) and the COmet Pressure Sensor (COPS). The spatial distribution of gas sources with their temporal variation allows one to construct surface maps for gas emissions and to evaluate integrated productions rates. For all species peak production rates and integrated productions rates per orbit are evaluated separately for the northern and the southern hemisphere. The nine most active emitting areas on the comet's surface are defined and their correlation to emissions for each of the species is discussed.

The ROSINA mass spectrometer DFMS on board ESA's Rosetta spacecraft detected the major isotopes of the noble gases argon, krypton, and xenon in the coma of comet 67P/Churyumov-Gerasimenko. Earlier, it has been shown that xenon exhibits an isotopic composition distinct from anywhere else in the solar system. However, argon isotopes, within error, were shown to be consistent with solar isotope abundances. This discrepancy suggested an additional exotic component of xenon in comet 67P/Churyumov-Gerasimenko. Here we show that also krypton exhibits an isotopic composition close to solar. Furthermore, we found a depletion compared to solar of argon with respect to krypton and of krypton with respect to xenon, which is a necessity to postulate an addition of exotic xenon in the comet.

Recent measurements carried out at comet 67P/C-G with the ${\it Rosetta}$ probe revealed that molecular oxygen, O$_2$, is the fourth most abundant molecule in comets. Models show that O$_2$ is likely of primordial nature, coming from the interstellar cloud from which our Solar System was formed. However, gaseous O$_2$ is an elusive molecule in the interstellar medium with only one detection towards quiescent molecular clouds, in the $\rho$ Oph A core. We perform a deep search for molecular oxygen, through the $2_1 - 0_1$ rotational transition at 234 GHz of its $^{16}$O$^{18}$O isotopologue, towards the warm compact gas surrounding the nearby Class 0 protostar IRAS 16293--2422 B with the ALMA interferometer. The targeted $^{16}$O$^{18}$O transition is surrounded by two brighter transitions at $\pm 1$ km s$^{-1}$ relative to the expected $^{16}$O$^{18}$O transition frequency. After subtraction of these two transitions, residual emission at a 3$\sigma$ level remains, but with a velocity offset of $0.3 - 0.5$ km s$^{-1}$ relative to the source velocity, rendering the detection "tentative". We derive the O$_2$ column density for two excitation temperatures $T_{\rm ex}$ of 125 and 300 K, as indicated by other molecules, in order to compare the O$_2$ abundance between IRAS16293 and comet 67P/C-G. Assuming that $^{16}$O$^{18}$O is not detected and using methanol CH$_3$OH as a reference species, we obtain a [O$_2$]/[CH$_3$OH] abundance ratio lower than $2-5$, depending on the assumed $T_{\rm ex}$, a three to four times lower abundance than the [O$_2$]/[CH$_3$OH] ratio of $5-15$ found in comet 67P/C-G. Such a low O$_2$ abundance could be explained by the lower temperature of the dense cloud precursor of IRAS16293 with respect to the one at the origin of our Solar System that prevented an efficient formation of O$_2$ in interstellar ices.

The evolutionary past of our Solar System can be pieced together by comparing analogous low-mass protostars with remnants of our Protosolar Nebula - comets. Sulphur-bearing molecules may be unique tracers of the joint evolution of the volatile and refractory components. ALMA Band 7 data from the large unbiased Protostellar Interferometric Line Survey (PILS) are used to search for S-bearing molecules in the outer disc-like structure, 60 au from IRAS 16293-2422 B, and are compared with data on 67P/C-G stemming from the ROSINA instrument aboard Rosetta. Species such as SO$_{2}$, SO, OCS, CS, H$_{2}$CS, H$_{2}$S and CH$_{3}$SH are detected via at least one of their isotopologues towards IRAS 16293-2422 B. The search reveals a first-time detection of OC$^{33}$S towards this source and a tentative first-time detection of C$^{36}$S towards a low-mass protostar. The data show that IRAS 16293-2422 B contains much more OCS than H$_{2}$S in comparison to 67P/C-G; meanwhile, the SO/SO$_{2}$ ratio is in close agreement between the two targets. IRAS 16293-2422 B has a CH$_{3}$SH/H$_{2}$CS ratio in range of that of our Solar System (differences by a factor of 0.7-5.3). It is suggested that the levels of UV radiation during the initial collapse of the systems may have varied and have potentially been higher for IRAS 16293-2422 B due to its binary nature; thereby, converting more H$_{2}$S into OCS. It remains to be conclusively tested if this also promotes the formation of S-bearing complex organics. Elevated UV levels of IRAS 16293-2422 B and a warmer birth cloud of our Solar System may jointly explain the variations between the two low-mass systems.

The chemical evolution of a star- and planet-forming system begins in the prestellar phase and proceeds across the subsequent evolutionary phases. The chemical trail from cores to protoplanetary disks to planetary embryos can be studied by comparing distant young protostars and comets in our Solar System. One particularly chemically rich system that is thought to be analogous to our own is the low-mass IRAS 16293-2422. ALMA-PILS observations have made the study of chemistry on the disk scales (< 100 AU) of this system possible. Under the assumption that comets are pristine tracers of the outer parts of the innate protosolar disk, it is possible to compare the composition of our infant Solar System to that of IRAS 16293-2422. The Rosetta mission has yielded a wealth of unique in situ measurements on comet 67P/C-G, making it the best probe to date. Herein, the initial comparisons in terms of the chemical composition and isotopic ratios are summarized. Much work is still to be carried out in the future as the analysis of both of these data sets is still ongoing.

Magnetohydrodynamics simulations have been carried out in studying the solar wind and cometary plasma interactions for decades. Various plasma boundaries have been simulated and compared well with observations for comet 1P/Halley. The Rosetta mission, which studies comet 67P/Churyumov-Gerasimenko, challenges our understanding of the solar wind and comet interactions. The Rosetta Plasma Consortium observed regions of very weak magnetic field outside the predicted diamagnetic cavity. In this paper, we simulate the inner coma with the Hall magnetohydrodynamics equations and show that the Hall effect is important in the inner coma environment. The magnetic field topology becomes complex and magnetic reconnection occurs on the dayside when the Hall effect is taken into account. The magnetic reconnection on the dayside can generate weak magnetic filed regions outside the global diamagnetic cavity, which may explain the Rosetta Plasma Consortium observations. We conclude that the substantial change in the inner coma environment is due to the fact that the ion inertial length (or gyro radius) is not much smaller than the size of the diamagnetic cavity.

Context. The Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA) was designed to measure the composition of the gas in the coma of comet 67P/Churyumov-Gerasimenko, the target of the European Space Agency's Rosetta mission. In addition to the volatiles, ROSINA measured refractories sputtered off the comet by the interaction of solar wind protons with the surface of the comet. Aims. The origin of different solar system materials is still heavily debated. Isotopic ratios can be used to distinguish between different reservoirs and investigate processes occurring during the formation of the solar system. Methods. ROSINA consisted of two mass spectrometers and a pressure sensor. In the ROSINA Double Focusing Mass Spectrometer (DFMS), the neutral gas of cometary origin was ionized and then deflected in an electric and a magnetic field that separated the ions based on their mass-to-charge ratio. The DFMS had a high mass resolution, dynamic range, and sensitivity that allowed detection of rare species and the known major volatiles. Results. We measured the relative abundance of all three stable silicon isotopes with the ROSINA instrument on board the Rosetta spacecraft. Furthermore, we measured $^{13}$C/$^{12}$C in C$_2$H$_4$, C$_2$H$_5$, and CO. The DFMS in situ measurements indicate that the average silicon isotopic composition shows depletion in the heavy isotopes $^{29}$Si and $^{30}$Si with respect to $^{28}$Si and solar abundances, while $^{13}$C to $^{12}$C is analytically indistinguishable from bulk planetary and meteorite compositions. Although the origin of the deficiency of the heavy silicon isotopes cannot be explained unambiguously, we discuss mechanisms that could have contributed to the measured depletion of the isotopes $^{29}$Si and $^{30}$Si.

Starting from several monthly data sets of Rosetta's COmetary Pressure Sensor we reconstruct the gas density in the coma around comet 67P/Churyumov-Gerasimenko. The underlying inverse gas model is constructed by fitting ten thousands of measurements to thousands of potential gas sources distributed across the entire nucleus surface. The ensuing self-consistent solution for the entire coma density and surface activity reproduces the temporal and spatial variations seen in the data for monthly periods with Pearson correlation coefficients of 0.93 and higher. For different seasonal illumination conditions before and after perihelion we observe a systematic shift of gas sources on the nucleus.

Because of the high fraction of refractory material present in comets, the heat produced by the radiogenic decay of elements such as aluminium and iron can be high enough to induce the loss of ultravolatile species such as nitrogen, argon or carbon monoxide during their accretion phase in the protosolar nebula. Here, we investigate how heat generated by the radioactive decay of 26Al and 60Fe influences the formation of comet 67P/Churyumov-Gerasimenko, as a function of its accretion time and size of parent body. We use an existing thermal evolution model that includes various phase transitions, heat transfer in the ice-dust matrix, and gas diffusion throughout the porous material, based on thermodynamic parameters derived from Rosetta observations. Two possibilities are considered: either, to account for its bilobate shape, 67P/Churyumov-Gerasimenko was assembled from two primordial ~2 kilometer-sized planetesimals, or it resulted from the disruption of a larger parent body with a size corresponding to that of comet Hale-Bopp (~70 km). To fully preserve its volatile content, we find that either 67P/Churyumov-Gerasimenko's formation was delayed between ~2.2 and 7.7 Myr after that of Ca-Al-rich Inclusions in the protosolar nebula or the comet's accretion phase took place over the entire time interval, depending on the primordial size of its parent body and the composition of the icy material considered. Our calculations suggest that the formation of 67P/Churyumov-Gerasimenko is consistent with both its accretion from primordial building blocks formed in the nebula or from debris issued from the disruption of a Hale-Bopp-like body.

The nature of the icy material accreted by comets during their formation in the outer regions of the protosolar nebula is a major open question in planetary science. Some scenarios of comet formation predict that these bodies agglomerated from crystalline ices condensed in the protosolar nebula. Concurrently, alternative scenarios suggest that comets accreted amorphous ice originating from the interstellar cloud or from the very distant regions of the protosolar nebula. On the basis of existing laboratory and modeling data, we find that the N$_2$/CO and Ar/CO ratios measured in the coma of the Jupiter family comet 67P/Churyumov-Gerasimenko by the ROSINA instrument aboard the European Space Agency's Rosetta spacecraft match those predicted for gases trapped in clathrates. If these measurements are representative of the bulk N$_2$/CO and Ar/CO ratios in 67P/Churyumov-Gerasimenko, it implies that the ices accreted by the comet formed in the nebula and do not originate from the interstellar medium, supporting the idea that the building blocks of outer solar system bodies have been formed from clathrates and possibly from pure crystalline ices. Moreover, because 67P/Churyumov-Gerasimenko is impoverished in Ar and N$_2$, the volatile enrichments observed in Jupiter's atmosphere cannot be explained solely via the accretion of building blocks with similar compositions and require an additional delivery source. A potential source may be the accretion of gas from the nebula that has been progressively enriched in heavy elements due to photoevaporation.

Molecular oxygen has been detected in the coma of comet 67P/Churyumov-Gerasimenko with abundances in the 1-10% range by the ROSINA-DFMS instrument on board the Rosetta spacecraft. Here we find that the radiolysis of icy grains in low-density environments such as the presolar cloud may induce the production of large amounts of molecular oxygen. We also show that molecular oxygen can be efficiently trapped in clathrates formed in the protosolar nebula, and that its incorporation as crystalline ice is highly implausible because this would imply much larger abundances of Ar and N2 than those observed in the coma. Assuming that radiolysis has been the only O2 production mechanism at work, we conclude that the formation of comet 67P/Churyumov-Gerasimenko is possible in a dense and early protosolar nebula in the framework of two extreme scenarios: (1) agglomeration from pristine amorphous icy grains/particles formed in ISM and (2) agglomeration from clathrates that formed during the disk's cooling. The former scenario is found consistent with the strong correlation between O2 and H2O observed in 67P/C-G's coma while the latter scenario requires that clathrates formed from ISM icy grains that crystallized when entering the protosolar nebula.

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. One crucial question is whether the O2 abundance is peculiar to comet 67P/Churyumov-Gerasimenko or Jupiter family comets in general or whether also 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 obtained during the flyby by the European Space Agency's Giotto probe at 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.

Grand Canonical Monte Carlo simulations are used to reproduce the N$_2$/CO ratio ranging between 1.7 $\times$ 10$^{-3}$ and 1.6 $\times$ 10$^{-2}$ observed {\it in situ} in the Jupiter family comet 67P/Churyumov-Gerasimenko by the ROSINA mass spectrometer aboard the Rosetta spacecraft, assuming that this body has been agglomerated from clathrates in the protosolar nebula. Simulations are done using an elaborated interatomic potentials for investigating the temperature dependence of the trapping within a multiple guest clathrate formed from a gas mixture of CO and N$_2$ in proportions corresponding to those expected for the protosolar nebula. By assuming that 67P/Churyumov-Gerasimenko agglomerated from clathrates, our calculations suggest the cometary grains must have been formed at temperatures ranging between $\sim$31.8 and 69.9 K in the protosolar nebula to match the N$_2$/CO ratio measured by the ROSINA mass spectrometer. The presence of clathrates in Jupiter family comets could then explain the potential N$_2$ depletion (factor up to $\sim$87 compared to the protosolar value) measured in 67P/Churyumov-Gerasimenko.

Remote sensing observations meet some limitations when used to study the bulk atmospheric composition of the giant planets of our solar system. A remarkable example of the superiority of in situ probe measurements is illustrated by the exploration of Jupiter, where key measurements such as the determination of the noble gases abundances and the precise measurement of the helium mixing ratio have only been made available through in situ measurements by the Galileo probe. This paper describes the main scientific goals to be addressed by the future in situ exploration of Saturn placing the Galileo probe exploration of Jupiter in a broader context and before the future probe exploration of the more remote ice giants. In situ exploration of Saturn's atmosphere addresses two broad themes that are discussed throughout this paper: first, the formation history of our solar system and second, the processes at play in planetary atmospheres. In this context, we detail the reasons why measurements of Saturn's bulk elemental and isotopic composition would place important constraints on the volatile reservoirs in the protosolar nebula. We also show that the in situ measurement of CO (or any other disequilibrium species that is depleted by reaction with water) in Saturn's upper troposphere would constrain its bulk O/H ratio. We highlight the key measurements required to distinguish competing theories to shed light on giant planet formation as a common process in planetary systems with potential applications to most extrasolar systems. In situ measurements of Saturn's stratospheric and tropospheric dynamics, chemistry and cloud-forming processes will provide access to phenomena unreachable to remote sensing studies. Different mission architectures are envisaged, which would benefit from strong international collaborations.