The outer solar system is home to a diverse range of objects, holding important clues about the formation and evolution of our solar system, the emergence and current distribution of life, and the physical processes controlling both our own and exoplanetary systems. This White Paper summarizes the Outer Planets Analysis Group's (OPAG's) priorities as they relate to the Decadal Survey. Taking into account the science to be achieved, the timing of solar system events, technological readiness, and programmatic factors, our mission recommendations are as follows. OPAG strongly endorses the completion and launch of the Europa Clipper mission, maintaining the science capabilities identified upon its selection, and a Juno extended mission at Jupiter. For the decade 2023-2032, OPAG endorses a new start for two directed missions: first, a mission to Neptune or Uranus (the ice giants) with atmospheric probe(s), and second, a life detection Ocean World mission, along with additional technological development for life detection. A Neptune mission is preferred because, while the Neptune and Uranus systems provide equally compelling opportunities to address the Origins and Processes Questions, Triton is a higher priority ocean world target than the Uranian satellites. The mission to Neptune or Uranus should fly first because a delay threatens key science objectives, and additional technological development is required for a directed life detection mission. Along with missions, we emphasize the necessity of maintaining a healthy Research and Analysis (R&A) program as well as a robust Earth-based observing program. OPAG's top two technology priorities are rapid development of a next-generation radioisotope power source for a mission to Neptune or Uranus, and development of key life detection technologies in support of an Ocean World mission.

We investigate the enrichment patterns of several delivery scenarios of the volatiles to the atmospheres of ice giants, having in mind that the only well constrained determination made remotely, i.e. the carbon abundance measurement, suggests that their envelopes possess highly supersolar metallicities, i.e. close to two orders of magnitude above that of the protosolar nebula. In the framework of the core accretion model, only the delivery of volatiles in solid forms (amorphous ice, clathrates, pure condensates) to these planets can account for the apparent supersolar metallicity of their envelopes. In contrast, because of the inward drift of icy particles through various snowlines, all mechanisms invoking the delivery of volatiles in vapor forms predict subsolar abundances in the envelopes of Uranus and Neptune. Alternatively, even if the disk instability mechanism remains questionable in our solar system, it might be consistent with the supersolar metallicities observed in Uranus and Neptune, assuming the two planets suffered subsequent erosion of their H-He envelopes. The enrichment patterns derived for each delivery scenario considered should be useful to interpret future in situ measurements by atmospheric entry probes.

Interstellar comets offer direct samples of volatiles from distant protoplanetary disks. 2I/Borisov is the first notably active interstellar comet discovered in our solar system[1]. Comets are condensed samples of the gas, ice, and dust that were in a star's protoplanetary disk during the formation of its planets and inform our understanding on how chemical compositions and abundances vary with distance from the central star. Their orbital migration moves volatiles[2], organic material, and prebiotic chemicals in their host system[3]. In our solar system, hundreds of comets have been observed remotely, and a few have been studied up close by space missions[4]. However, knowledge of extrasolar comets has been limited to what could be gleaned from distant, unresolved observations of cometary regions around other stars, with only one detection of carbon monoxide[5]. Here we report that the coma of 2I/Borisov contains significantly more CO than H2O gas, with abundances of at least 173%, more than three times higher than previously measured for any comet in the inner (<2.5 au) solar system[4]. Our ultraviolet observations of 2I/Borisov provide the first glimpse into the ice content and chemical composition of the protoplanetary disk of another star that is substantially different from our own.

This White Paper presents a mission concept called MIRECLE - the Mid-InfraRed Exoplanet CLimate Explorer. With a moderately sized aperture of 2 meters, broad wavelength coverage (4 - 25 um), and next generation instruments, MIRECLE will be capable of efficiently characterizing a statistically significant sample of terrestrial planets, many of which will be in their host stars's habitable zones. Spectroscopic characterization of terrestrial atmospheres will provide constraints for the distribution of planets with tenuous vs. substantial atmospheres, on the inner and outer edges of the habitable zone, and climate models to assess the potential for habitability. For the few brightest targets, the detection of specific combinations of molecules would provide evidence of biosignatures. For all other targets, this comprehensive survey would filter out the airless, desiccated, or lifeless worlds, thus providing a subset of potentially habitable worlds ready for in-depth atmospheric characterization using a larger aperture telescope.

In this white paper, we present a cross-section of important scientific questions that remain partially or completely unanswered, ranging from Titan exosphere to the deep interior, and we detail which instrumentation and mission scenarios should be used to answer them. Our intention is to formulate the science goals for the next generation of planetary missions to Titan in order to prepare the future exploration of the moon. The ESA L-class mission concept that we propose is composed of a Titan orbiter and at least an in situ element (lake lander and/or drone(s)).

Remote sensing observations suffer significant limitations when used to study the bulk atmospheric composition of the giant planets of our solar system. This impacts our knowledge of the formation of these planets and the physics of their atmospheres. A remarkable example of the superiority of in situ probe measurements was 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 were only made available through in situ measurements by the Galileo probe. Here we describe the main scientific goals to be addressed by the future in situ exploration of Saturn, Uranus, and Neptune, placing the Galileo probe exploration of Jupiter in a broader context. An atmospheric entry probe targeting the 10-bar level would yield insight into two broad themes: i) the formation history of the giant planets and that of the Solar System, and ii) the processes at play in planetary atmospheres. The probe would descend under parachute to measure composition, structure, and dynamics, with data returned to Earth using a Carrier Relay Spacecraft as a relay station. An atmospheric probe could represent a significant ESA contribution to a future NASA New Frontiers or flagship mission to be launched toward Saturn, Uranus, and/or Neptune.

Planned astronomical observatories of the 2020s will be capable of obtaining reflected light photometry and spectroscopy of cool extrasolar giant planets. Here we explain that such data are valuable both for understanding the origin and evolution of giant planets as a whole and for preparing for the interpretation of similar datasets from potentially habitable extrasolar terrestrial planets in the decades to follow.

Unexpected dynamic phenomena have surprised solar system observers in the past and have led to important discoveries about solar system workings. Observations at the initial stages of these events provide crucial information on the physical processes at work. We advocate for long-term/permanent programs on ground-based and space-based telescopes of all sizes - including Extremely Large Telescopes (ELTs) - to conduct observations of high-priority dynamic phenomena, based on a predefined set of triggering conditions. These programs will ensure that the best initial dataset of the triggering event are taken; separate additional observing programs will be required to study the temporal evolution of these phenomena. While not a comprehensive list, the following are notional examples of phenomena that are rare, that cannot be anticipated, and that provide high-impact advances to our understandings of planetary processes. Examples include: new cryovolcanic eruptions or plumes on ocean worlds; impacts on Jupiter, Saturn, Uranus, or Neptune; extreme eruptions on Io; convective superstorms on Saturn, Uranus, or Neptune; collisions within the asteroid belt or other small-body populations; discovery of an interstellar object passing through our solar system (e.g. 'Oumuamua); and responses of planetary atmospheres to major solar flares or coronal mass ejections.

The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events.

Future remote sensing of exoplanets will be enhanced by a thorough investigation of our solar system Ice Giants (Neptune-size planets). What can the configuration of the magnetic field tell us (remotely) about the interior, and what implications does that field have for the structure of the magnetosphere; energy input into the atmosphere, and surface geophysics (for example surface weathering of satellites that might harbour sub-surface oceans). How can monitoring of auroral emission help inform future remote observations of emission from exoplanets? Our Solar System provides the only laboratory in which we can perform in-situ experiments to understand exoplanet formation, dynamos, systems and magnetospheres.

The Origins Space Telescope (OST) is one of four mission concepts currently being studied by NASA in preparation for the Astrophysics 2020 Decadal Survey. With active cooling (~4 K), OST will be sensitive in mid- to far-IR wavelengths, using imaging and spectroscopy to probe the furthest reaches of our galaxies, trace the path of water through star and planet formation, and place thermochemical constraints on the atmospheres of exoplanets ranging in size from Jupiter to Earth. This contribution to the Exoplanet Science Strategy committee discusses the significant advancements that the OST Mid-Infrared Imager, Spectrometer, and Coronagraph (MISC) instrument can make in studying cool planetary atmospheres. We particularly focus on the atmospheres of transiting rocky planets in the habitable zones of mid-to-late M stars. We discuss how OST thermal infrared observations can significantly enhance our understanding of the temperature structure and molecular abundances of biologically interesting gases on these worlds, including O3, CH4, H2O, and CO2.

This whitepaper discusses the diversity of exoplanets that could be detected by future observations, so that comparative exoplanetology can be performed in the upcoming era of large space-based flagship missions. The primary focus will be on characterizing Earth-like worlds around Sun-like stars. However, we will also be able to characterize companion planets in the system simultaneously. This will not only provide a contextual picture with regards to our Solar system, but also presents a unique opportunity to observe size dependent planetary atmospheres at different orbital distances. We propose a preliminary scheme based on chemical behavior of gases and condensates in a planet's atmosphere that classifies them with respect to planetary radius and incident stellar flux.

It is currently unknown how common life is on exoplanets, or how long planets can remain viable for life. To date, we have a superficial notion of habitability, a necessary first step, but so far lacking an understanding of the detailed interaction between stars and planets over geological timescales, dynamical evolution of planetary systems, and atmospheric evolution on planets in other systems. A planet mass, net insolation, and atmospheric composition alone are insufficient to determine the probability that life on a planet could arise or be detected. The latter set of planetary considerations, among others, underpin the concept of the habitable zone (HZ), defined as the circumstellar region where standing bodies of liquid water could be supported on the surface of a rocky planet. However, stars within the same spectral class are often treated in the same way in HZ studies, without any regard for variations in activity among individual stars. Such formulations ignore differences in how nonthermal emission and magnetic energy of transient events in different stars affect the ability of an exoplanet to retain its atmosphere.In the last few years there has been a growing appreciation that the atmospheric chemistry, and even retention of an atmosphere in many cases, depends critically on the high-energy radiation and particle environments around these stars. Indeed, recent studies have shown stellar activity and the extreme space weather, such as that created by the frequent flares and coronal mass ejections (CMEs) from the active stars and young Sun, may have profoundly affected the chemistry and climate and thus habitability of the early Earth and terrestrial type exoplanets. The goal of this white paper is to identify and describe promising key research goals to aid the field of the exoplanetary habitability for the next 20 years.

Strong electron cooling on the neutral gas in cometary comae has been predicted for a long time, but actual measurements of low electron temperature are scarce. We present in situ measurements of plasma density, electron temperature and spacecraft potential by the Rosetta Langmuir probe instrument, LAP. Data acquired within a few hundred km from the nucleus are dominated by a warm component with electron temperature typically 5--10 eV at all heliocentric distances covered (1.25 to 3.83 AU). A cold component, with temperature no higher than about 0.1 eV, appears in the data as short (few to few tens of seconds) pulses of high probe current, indicating local enhancement of plasma density as well as a decrease in electron temperature. These pulses first appeared around 3 AU and were seen for longer periods close to perihelion. The general pattern of pulse appearance follows that of neutral gas and plasma density. We have not identified any periods with only cold electrons present. The electron flux to Rosetta was always dominated by higher energies, driving the spacecraft potential to order -10 V. The warm (5--10 eV) electron population is interpreted as electrons retaining the energy they obtained when released in the ionisation process. The sometimes observed cold populations with electron temperatures below 0.1 eV verify collisional cooling in the coma. The cold electrons were only observed together with the warm population. The general appearance of the cold population appears to be consistent with a Haser-like model, implicitly supporting also the coupling of ions to the neutral gas. The expanding cold plasma is unstable, forming filaments that we observe as pulses.

We use measurements from the Rosetta plasma consortium (RPC) Langmuir probe (LAP) and mutual impedance probe (MIP) to study the spatial distribution of low-energy plasma in the near-nucleus coma of comet 67P/Churyumov-Gerasimenko. The spatial distribution is highly structured with the highest density in the summer hemisphere and above the region connecting the two main lobes of the comet, i.e. the neck region. There is a clear correlation with the neutral density and the plasma to neutral density ratio is found to be about 1-2x10^-6, at a cometocentric distance of 10 km and at 3.1 AU from the sun. A clear 6.2 h modulation of the plasma is seen as the neck is exposed twice per rotation. The electron density of the collisonless plasma within 260 km from the nucleus falls of with radial distance as about 1/r. The spatial structure indicates that local ionization of neutral gas is the dominant source of low-energy plasma around the comet.

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.

On November 12, 2014, the ESA/Rosetta descent module Philae landed on the Abydos site of comet 67P/Churyumov-Gerasimenko. Aboard this module, the Ptolemy mass spectrometer measured a CO/CO2 ratio of 0.07 +/- 0.04 which differs substantially from the value obtained in the coma by the Rosetta/ROSINA instrument, suggesting a heterogeneity in the comet nucleus. To understand this difference, we investigated the physico-chemical properties of the Abydos subsurface leading to CO/CO2 ratios close to that observed by Ptolemy at the surface of this region. We used a comet nucleus model that takes into account different water ice phase changes (amorphous ice, crystalline ice and clathrates), as well as diffusion of molecules throughout the pores of the matrix. The input parameters of the model were optimized for the Abydos site and the ROSINA CO/CO2 measured ratio is assumed to correspond to the bulk value in the nucleus. We find that all considered structures of water ice are able to reproduce the Ptolemy observation with a time difference not exceeding ~50 days, i.e. lower than ~2% on 67P/Churyumov-Gerasimenko's orbital period. The suspected heterogeneity of 67P/Churyumov-Gerasimenko's nucleus is also found possible only if it is constituted of crystalline ices. If the icy phase is made of amorphous ice or clathrates, the difference between Ptolemy and ROSINA's measurements would rather originate from the spatial variations in illumination on the nucleus surface. An eventual new measurement of the CO/CO2 ratio at Abydos by Ptolemy could be decisive to distinguish between the three water ice structures.

Bimolecular diffusion coefficients are important parameters used by atmospheric models to calculate altitude profiles of minor constituents in an atmosphere. Unfortunately, laboratory measurements of these coefficients were never conducted at temperature conditions relevant to the atmosphere of Titan. Here we conduct a detailed uncertainty analysis of the bimolecular diffusion coefficient parameters as applied to Titan's upper atmosphere to provide a better understanding of the impact of uncertainty for this parameter on models. Because temperature and pressure conditions are much lower than the laboratory conditions in which bimolecular diffusion parameters were measured, we apply a Bayesian framework, a problem-agnostic framework, to determine parameter estimates and associated uncertainties. We solve the Bayesian calibration problem using the open-source QUESO library which also performs a propagation of uncertainties in the calibrated parameters to temperature and pressure conditions observed in Titan's upper atmosphere. Our results show that, after propagating uncertainty through the Massman model, the uncertainty in molecular diffusion is highly correlated to temperature and we observe no noticeable correlation with pressure. We propagate the calibrated molecular diffusion estimate and associated uncertainty to obtain an estimate with uncertainty due to bimolecular diffusion for the methane molar fraction as a function of altitude. Results show that the uncertainty in methane abundance due to molecular diffusion is in general small compared to eddy diffusion and the chemical kinetics description. However, methane abundance is most sensitive to uncertainty in molecular diffusion above 1200 km where the errors are nontrivial and could have important implications for scientific research based on diffusion models in this altitude range.

The recent derivation of a lower limit for the $^{14}$N/$^{15}$N ratio in Saturn's ammonia, which is found to be consistent with the Jovian value, prompted us to revise models of Saturn's formation using as constraints the supersolar abundances of heavy elements measured in its atmosphere. Here we find that it is possible to account for both Saturn's chemical and isotopic compositions if one assumes the formation of its building blocks at $\sim$45 K in the protosolar nebula, provided that the O abundance was $\sim$2.6 times protosolar in its feeding zone. To do so, we used a statistical thermodynamic model to investigate the composition of the clathrate phase that formed during the cooling of the protosolar nebula and from which the building blocks of Saturn were agglomerated. We find that Saturn's O/H is at least $\sim$34.9 times protosolar and that the corresponding mass of heavy elements ($\sim$43.1 \Mearth) is within the range predicted by semi-convective interior models.

We use a statistical-thermodynamic model to investigate the formation and composition of noble-gas-rich clathrates on Pluto's surface. By considering an atmospheric composition close to that of today's Pluto and a broad range of surface pressures, we find that Ar, Kr and Xe can be efficiently trapped in clathrates if they formed at the surface, in a way similar to what has been proposed for Titan. The formation on Pluto of clathrates rich in noble gases could then induce a strong decrease in their atmospheric abundances relative to their initial values. A clathrate thickness of order of a few centimeters globally averaged on the planet is enough to trap all Ar, Kr and Xe if these noble gases were in protosolar proportions in Pluto's early atmosphere. Because atmospheric escape over an extended period of time (millions of years) should lead to a noble gas abundance that either remains constant or increases with time, we find that a potential depletion of Ar, Kr and Xe in the atmosphere would best be explained by their trapping in clathrates. A key observational test is the measurement of Ar since the Alice UV spectrometer aboard the New Horizons spacecraft will be sensitive enough to detect its abundance $\sim$10 times smaller than in the case considered here.

A striking feature of the atmosphere of Titan is that no heavy noble gases other than argon were detected by the Gas Chromatograph Mass Spectrometer (GCMS) aboard the Huygens probe during its descent to Titan's surface in January 2005. Here we provide an explanation of the mysterious absence or rarity of these noble gases in Titan's atmosphere: the thermodynamic conditions prevailing at the surface-atmosphere interface of the satellite allow the formation of multiple guest clathrates that preferentially store some species, including all heavy noble gases, over others. The clean water ice needed for formation of these clathrates could be delivered by successive episodes of cryovolcanic lavas that have been hypothesized to regularly cover the surface of Titan. The formation of clathrates in the porous lavas and their propensity for trapping Ar, Kr and Xe would progressively remove these species from the atmosphere of Titan over its history. In some circumstances, a global clathrate crust with an average thickness not exceeding a few meters could be sufficient on Titan for a complete removal of the heavy noble gases from the atmosphere.

The origin of Titan's atmospheric methane is a key issue for understanding the origin of the Saturnian satellite system. It has been proposed that serpentinization reactions in Titan's interior could lead to the formation of the observed methane. Meanwhile, alternative scenarios suggest that methane was incorporated in Titan's planetesimals before its formation. Here, we point out that serpentinization reactions in Titan's interior are not able to reproduce the deuterium over hydrogen (D/H) ratio observed at present in methane in its atmosphere, and would require a maximum D/H ratio in Titan's water ice 30% lower than the value likely acquired by the satellite during its formation, based on Cassini observations at Enceladus. Alternatively, production of methane in Titan's interior via radiolytic reactions with water can be envisaged but the associated production rates remain uncertain. On the other hand, a mechanism that easily explains the presence of large amounts of methane trapped in Titan in a way consistent with its measured atmospheric D/H ratio is its direct capture in the satellite's planetesimals at the time of their formation in the solar nebula. In this case, the mass of methane trapped in Titan's interior can be up to 1,300 times the current mass of atmospheric methane.