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)).

We have analyzed the complete Visual and Infrared Mapping Spectrometer (VIMS) data archive of Titan. Our objective is to build global surface cartographic products, by combining all the data gathered during the 127 targeted flybys of Titan into synthetic global maps interpolated on a grid at 32 pixels per degree (~1.4 km/pixel at the equator), in seven infrared spectral atmospheric windows. Multispectral summary images have been computed for each single VIMS cube in order to rapidly identify their scientific content and assess their quality. These summary images are made available to the community on a public website (vims.univ-nantes.fr). The global mapping work faced several challenges due to the strong absorbing and scattering effects of the atmosphere coupled to the changing observing conditions linked to the orbital tour of the Cassini mission. We determined a surface photometric function which accounts for variations in incidence, emergence and phase angles, and which is able to mitigate brightness variations linked to the viewing geometry of the flybys. The atmospheric contribution has been reduced using the subtraction of the methane absorption band wings, considered as proxies for atmospheric haze scattering. We present a new global three color composite map of band ratios (red: 1.59/1.27 {\mu}m; green: 2.03/1.27 {\mu}m; blue: 1.27/1.08 {\mu}m), which has also been empirically corrected from an airmass (the solar photon path length through the atmosphere) dependence. This map provides a detailed global color view of Titan's surface partially corrected from the atmosphere and gives a global insight of the spectral variability, with the equatorial dunes fields appearing in brownish tones, and several occurrences of bluish tones localized in areas such as Sinlap, Menvra and Selk craters. (...)

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.

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.

The James Webb Space Telescope (JWST), scheduled for launch in 2018, is the successor to the Hubble Space Telescope (HST) but with a significantly larger aperture (6.5 m) and advanced instrumentation focusing on infrared science (0.6-28.0 $\mu$m ). In this paper we examine the potential for scientific investigation of Titan using JWST, primarily with three of the four instruments: NIRSpec, NIRCam and MIRI, noting that science with NIRISS will be complementary. Five core scientific themes are identified: (i) surface (ii) tropospheric clouds (iii) tropospheric gases (iv) stratospheric composition and (v) stratospheric hazes. We discuss each theme in depth, including the scientific purpose, capabilities and limitations of the instrument suite, and suggested observing schemes. We pay particular attention to saturation, which is a problem for all three instruments, but may be alleviated for NIRCam through use of selecting small sub-arrays of the detectors - sufficient to encompass Titan, but with significantly faster read-out times. We find that JWST has very significant potential for advancing Titan science, with a spectral resolution exceeding the Cassini instrument suite at near-infrared wavelengths, and a spatial resolution exceeding HST at the same wavelengths. In particular, JWST will be valuable for time-domain monitoring of Titan, given a five to ten year expected lifetime for the observatory, for example monitoring the seasonal appearance of clouds. JWST observations in the post-Cassini period will complement those of other large facilities such as HST, ALMA, SOFIA and next-generation ground-based telescopes (TMT, GMT, EELT).

Titan's polar surface is dotted with hundreds of lacustrine depressions. Based on the hypothesis that they are karstic in origin, we aim at determining the efficiency of surface dissolution as a landshaping process on Titan, in a comparative planetology perspective with the Earth as reference. Our approach is based on the calculation of solutional denudation rates and allow inference of formation timescales for topographic depressions developed by chemical erosion on both planetary bodies. The model depends on the solubility of solids in liquids, the density of solids and liquids, and the average annual net rainfall rates. We compute and compare the denudation rates of pure solid organics in liquid hydrocarbons and of minerals in liquid water over Titan and Earth timescales. We then investigate the denudation rates of a superficial organic layer in liquid methane over one Titan year. At this timescale, such a layer on Titan would behave like salts or carbonates on Earth depending on its composition, which means that dissolution processes would likely occur but would be 30 times slower on Titan compared to the Earth due to the seasonality of precipitation. Assuming an average depth of 100 m for Titan's lacustrine depressions, these could have developed in a few tens of millions of years at polar latitudes higher than 70{\deg} N and S, and a few hundreds of million years at lower polar latitudes. The ages determined are consistent with the youth of the surface (<1 Gyr) and the repartition of dissolution-related landforms on Titan.