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.

We report the detection of a spectral signature observed at 682 cm$^{-1}$ by the Cassini Composite Infrared Spectrometer (CIRS) in nadir and limb geometry observations of Titan's southern stratospheric polar region in the middle of southern fall, while stratospheric temperatures are the coldest since the beginning of the Cassini mission. The 682 cm$^{-1}$ signature, which is only observed below an altitude of 300-km, is at least partly attributed to the benzene (C$_6$H$_6$) ice $\nu_{4}$ C-H bending mode. While we first observed it in CIRS nadir spectra of the southern polar region in early 2013, we focus here on the study of nadir data acquired in May 2013, which have a more favorable observation geometry. We derived the C$_6$H$_6$ ice mass mixing ratio in 5{\deg}S latitude bins from the south pole to 65{\deg}S and infer the C$_6$H$_6$ cloud top altitude to be located deeper with increasing distance from the pole. We additionally analyzed limb data acquired in March 2015, which were the first limb dataset available after the May 2013 nadir observation, in order to infer a vertical profile of its mass mixing ratio in the 0.1 - 1 mbar region (250 - 170 km). We derive an upper limit of $\sim$1.5 $\mu$m for the equivalent radius of pure C$_6$H$_6$ ice particles from the shape of the observed emission band. Several other unidentified signatures are observed near 687 and 702 cm$^{-1}$ and possibly 695 cm$^{-1}$, which could also be due to ice spectral signatures as they are observed in the deep stratosphere at pressure levels similar to the C$_6$H$_6$ ice ones. We could not reproduce these signatures with pure nitrile ice (HCN, HC$_3$N,CH$_3$CN, C$_2$H$_5$CN and C$_2$N$_2$) spectra available in the literature except the 695 cm$^{-1}$ feature that could possibly be due to C$_2$H$_3$CN ice.

We have analyzed the complete archive of the Visual and Infrared Mapping Spectrometer (VIMS) data in order to monitor and analyze the evolution of the clouds and haze coverage at both poles of Titan during the entire Cassini mission. Our objective is to give a cartographic synopsis from a VIMS perspective, to provide a global view of the seasonal evolution of Titan's atmosphere over the poles. We leave the detailed comparison with the Imaging Science Subsystem (ISS) and the Composite Infrared Spectrometer (CIRS) data sets to further studies. We have computed global hyperspectral mosaics for each of the 127 targeted flybys of Titan to produce synthetic color maps emphasizing the main atmospheric features. The north pole appears fully covered by a huge cloud as soon as the first observations in 2004 and up to the equinox in 2009 (Le Mou\'elic et al. 2012). The northern skies then became progressively clearer, after the circulation turnover in 2009, revealing the underlying lakes and seas to the optical instruments up to 2017. The reverse situation is observed over the south pole, which was mostly clear of such a high obscuring cloud during the first years of the mission, but started to develop a polar cloud in 2012. This feature grew up month after month until the end of the mission in 2017, with a poleward latitudinal extent of 75$^\circ$S in 2013 up to 58$^\circ$S in April 2017. Thanks to the spectral capabilities of VIMS, we have detected HCN spectral signatures over the north pole in almost all flybys between 2004 and 2008. These HCN signatures started then to show up over the south pole in almost all flybys between 2012 and 2017, so perfectly matching the timing and spatial extent of the northern and southern polar atmospheric features.

Titan's 'detached' haze, seen in Voyager images in 1980 and 1981 and monitored by the Cassini Imaging Science Subsystem (ISS), during the period 2004-2017, provides a measure of seasonal activity in Titan's mesosphere with observations over almost one half of Saturn's seasonal cycle. Here we report on retrieved haze extinction profiles that reveal a depleted layer that visually manifests as a thin layer detached from the main haze below. Our new measurements show the disappearance of the feature in 2012 and its reappearance in 2016, as well as details after the reappearance. These observations highlight the dynamical nature of the detached haze. The reappearance appears congruent but more complex than previously described by climate models. It occurs in two steps, first as haze reappearing at $450\pm20$ km and one year later at $510\pm20$ km. These new observations provide additional tight and valuable constraints about the underlying mechanisms, especially for Titan's mesosphere, that control Titan's haze cycle.

Titan's atmosphere is optically thick and hides the surface and the lower layers from the view at almost all wavelengths. However, because gaseous absorptions are spectrally selective, some narrow spectral intervals are relatively transparent and allow to probe the surface. To use these intervals (called windows) a good knowledge of atmospheric absorption is necessary. Once gas spectroscopic linelists are well established, the absorption inside windows depends on the way the far wings of the methane absorption lines are cut-off. We know that the intensity in all the windows can be explained with the same cut-off parameters, except for the window at 2 {\mu}m. This discrepancy is generally treated with a workaround which consists in using a different cut-off description for this specific window. This window is relatively transparent and surface may have specific spectral signatures that could be detected. Thus, a good knowledge of atmosphere opacities is essential and our scope is to better understand what causes the difference between the 2 {\mu}m window and the other windows. In this work, we used scattered light at the limb and transmissions in occultation observed with VIMS (Visual Infrared Mapping Spectrometer) onboard Cassini, around the 2 {\mu}m window. Data shows an absorption feature that participates to the shape of this window. Our atmospheric model fits well the VIMS data at 2 {\mu}m with the same cut-off than for the other windows, provided an additional absorption is introduced in the middle of the window around ~2.065 {\mu}m. It explains well the discrepancy between the cut-off used at 2 {\mu}m, and we show that a gas with a fairly constant mixing ratio, possibly ethane, may be the cause of this absorption. Finally, we studied the impact of this absorption on the retrieval of the surface reflectivity and found that it is significant.

UV observations with Cassini ISS Narrow Angle Camera of Titan's detached haze is an excellent tool to probe its aerosols content without being affected by the gas or the multiple scattering. Unfortunately, its low extent in altitude requires a high resolution calibration and limits the number of images available in the Cassini dataset. However, we show that it is possible to extract on each profile the local maximum of intensity of this layer and confirm its stability at $500 \pm 8$ km during the 2005-2007 period for all latitudes lower than 45$^\circ$N. Using the fractal aggregate scattering model of Tomasko et al. (2008) and a single scattering radiative transfer model, it is possible to derive the optical properties required to explain the observations made at different phase angles. Our results indicates that the aerosols have at least ten monomers of 60 nm radius, while the typical tangential column number density is about $2\cdot 10^{10}$ agg.m$^{-2}$. Moreover, we demonstrate that these properties are constant within the error bars in the southern hemisphere of Titan over the observed time period. In the northern hemisphere, the size of the aerosols tend to decrease relatively to the southern hemisphere and are associated with a higher tangential opacity. However, the lower number of observations available in this region due to the orbital constraints is a limiting factor in the accuracy of these results. Assuming a fixed homogeneous content we notice that the tangential opacity can fluctuate up to a factor 3 among the observations at the equator. These variations could be linked with short scale temporal and/or longitudinal events changing the local density of the layer.

Numerous geological features that could be evaporitic in origin have been identified on the surface of Titan. Although they seem to be water-ice poor, their main properties -chemical composition, thickness, stratification- are essentially unknown. In this paper, which follows on a previous one focusing on the surface composition (Cordier et al., 2013), we provide some answers to these questions derived from a new model. This model, based on the up-to-date thermodynamic theory known as "PC-SAFT", has been validated with available laboratory measurements and specifically developed for our purpose. 1-D models confirm the possibility of an acetylene and/or butane enriched central layer of evaporitic deposit. The estimated thickness of this acetylene-butane layer could explain the strong RADAR brightness of the evaporites. The 2-D computations indicate an accumulation of poorly soluble species at the deposit's margin. Among these species, HCN or aerosols similar to tholins could play a dominant role. Our model predicts the existence of chemically trimodal "bathtub rings" which is consistent with what it is observed at the south polar lake Ontario Lacus. This work also provides plausible explanations to the lack of evaporites in the south polar region and to the high radar reflectivity of dry lakebeds.

Between 2004 and 2007 the instruments of the CASSINI spacecraft discovered hydrocarbon lakes in the polar regions of Titan. We have developed a lake-atmosphere equilibrium model allowing the determination of the chemical composition of these liquid areas. The model is based on uncertain thermodynamic data and precipitation rates of organic species predicted to be present in the lakes and seas that are subject to spatial and temporal variations. Here we explore and discuss the influence of these uncertainties and variations. The errors and uncertainties relevant to thermodynamic data are simulated via Monte-Carlo simulations. Global Circulation Models (GCM) are also employed in order to investigate the possibility of chemical asymmetry between the south and the north poles, due to differences in precipitation rates. We find that mole fractions of compounds in the liquid phase have a high sensitivity to thermodynamic data used as inputs, in particular molar volumes and enthalpies of vaporization. When we combine all considered uncertainties, the ranges of obtained mole fractions are rather large (up to ~8500%) but the distributions of values are narrow. The relative standard deviations remain between 10% and ~300% depending on the compound considered. Compared to other sources of uncertainties and variability, deviation caused by surface pressure variations are clearly negligible, remaining of the order of a few percent up to ~20%. Moreover no significant difference is found between the composition of lakes located in north and south poles. Because the theory of regular solutions employed here is sensitive to thermodynamic data and is not suitable for polar molecules such as HCN and CH3CN, our work strongly underlines the need for experimental simulations and the improvement of Titan's atmospheric models.

Clouds on Titan result from the condensation of methane and ethane and, as on other planets, are primarily structured by circulation of the atmosphere. At present, cloud activity mainly occurs in the southern (summer) hemisphere, arising near the pole and at mid-latitudes from cumulus updrafts triggered by surface heating and/or local methane sources, and at the north (winter) pole, resulting from the subsidence and condensation of ethane-rich air into the colder troposphere. General circulation models predict that this distribution should change with the seasons on a 15-year timescale, and that clouds should develop under certain circumstances at temperate latitudes (~40\degree) in the winter hemisphere. The models, however, have hitherto been poorly constrained and their long-term predictions have not yet been observationally verified. Here we report that the global spatial cloud coverage on Titan is in general agreement with the models, confirming that cloud activity is mainly controlled by the global circulation. The non-detection of clouds at latitude ~40\degree N and the persistence of the southern clouds while the southern summer is ending are, however, both contrary to predictions. This suggests that Titan's equator-to-pole thermal contrast is overestimated in the models and that its atmosphere responds to the seasonal forcing with a greater inertia than expected.