The texture, composition, and morphology of dunes observed in the equatorial regions of Titan may reflect present and/or past climatic conditions. Determining the physio-chemical properties and the morphodynamics of Titan's dunes is therefore essential to understanding of the climatic and geological history of the largest moon of Saturn. We quantitatively derived average surface properties of dune and interdune areas (texture, composition) from modeling of the microwave backscattered signal and Monte-Carlo inversion of the despeckled Cassini/SAR data over Titan's three largest sand seas: Belet, Shangri-La and Fensal. We present the first analysis of the backscatter functions extracted from despeckled SAR images that cover such a large range in incidence angles, including data from the beginning of the Cassini mission up to its Grand Finale. We show that dunes and interdunes have significantly different physical properties. Dunes are found to be more microwave absorbent than interdunes. Additionally, potential secondary bedforms, such as ripples and avalanches, may have been detected, providing potential evidence for currently active dunes and sediment transport. Our modelling shows that the interdunes have multi-scale roughnesses with higher dielectric constants than the dunes which have a low dielectric constant consistent with organic sand. The radar brightness of the interdunes can be explained by the presence of a shallow layer of significantly larger organic grains, possibly non-mobilized by the winds. {Together, our} findings suggest that Titan's sand seas evolve under the current multi-directional wind regimes with dunes that elongate with their crests aligned in the residual drift direction.

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

Unraveling the stratigraphic record is the key to understanding ancient climate and past climate changes on Mars. River deposits when placed in stratigraphic order could constrain the number, magnitudes, and durations of the wettest climates in Mars history. We establish the stratigraphic context of river deposits in Aeolis Dorsa sedimentary basin, 10E of Gale crater. Here, wind has exhumed a stratigraphic section of >=4 unconformity-bounded sedimentary rock packages, recording >=3 distinct episodes of surface runoff. Early deposits (>700m thick) are embayed by river deposits (>400m), which are in turn unconformably draped by fan-shaped deposits (<100m) which we interpret as alluvial fans. Yardang-forming deposits (>900 m) unconformably drape all previous deposits. River deposits embay a dissected sedimentary-rock landscape, and comprise >=2 distinguishable units. The total interval spanned by river deposits is >(1x10^6-2x10^7) yr; more if we include alluvial-fan deposits. Alluvial-fan deposits unconformably postdate thrust faults which crosscut river deposits. We infer a relatively dry interval of >4x10^7 yr after river deposits formed and before fan-shaped deposits formed. The time gap between the end of river deposition and the onset of yardang-forming deposits is constrained to >10^8 yr by the density of impact craters embedded at the unconformity. We correlate yardang-forming deposits to the upper layers of Gale crater's mound (Mt. Sharp/Aeolis Mons), and fan-shaped deposits to Peace Vallis fan. Alternations between periods of low vs. high mean obliquity may have modulated erosion-deposition cycling in Aeolis. This is consistent with results from an ensemble of simulations of Solar System orbital evolution and the resulting history of Mars obliquity. Almost all simulations yield intervals of continuously low mean Mars obliquity that are long enough to match our unconformity data.

The decay of the martian atmosphere - which is dominated by carbon dioxide - is a component of the long-term environmental change on Mars from a climate that once allowed rivers to flow to the cold and dry conditions of today. The minimum size of craters serves as a proxy for palaeopressure of planetary atmospheres, because thinner atmospheres permit smaller objects to reach the surface at high velocities and form craters. The Aeolis Dorsa region near Gale crater on Mars contains a high density of preserved ancient craters interbedded with river deposits and thus can provide constraints on atmospheric density around the time of fluvial activity. Here we use high-resolution orthophotos and digital terrain models from the Mars Reconnaissance Orbiter to identify ancient craters in Aeolis Dorsa that date to about 3.6 Gyr ago and compare their size distribution with models of atmospheric filtering of impactors. We obtain an upper limit of 0.9$\pm$0.1 bar, rising to 1.9$\pm$0.2 bar if rimmed circular mesas - interpreted to be erosionally-resistant fills of floors of impact craters - are excluded. We assume target properties appropriate for desert alluvium: if sediment had rock-like rock-mass strength similar to bedrock at the time of impact, the upper limit increases by a factor of up to two. If Mars did not have a stable multibar atmosphere at the time that the rivers were flowing - as suggested by our results - then the warm and wet CO2/H2O greenhouse is ruled out, and long-term average temperatures were most likely below freezing.

River deposits are one of the main lines of evidence that tell us that Mars once had a climate different from today, and so changes in river deposits with time tell us something about how Mars climate changed with time. In this study, we focus in on one sedimentary basin - Aeolis Dorsa - which contains an exceptionally high number of exceptionally well-preserved river deposits that appear to have formed over an interval of $>$0.5 Myr. We use changes in the river deposits' scale with stratigraphic elevation as a proxy for changes in river paleodischarge. Meander wavelengths tighten upwards and channel widths narrow upwards, and there is some evidence for a return to wide large-wavelength channels higher in the stratigraphy. Meander wavelength and channel width covary with stratigraphic elevation. The factor of 1.5-2 variations in paleochannel dimensions with stratigraphic elevation correspond to $\sim$2.6-fold variability in bank-forming discharge (using standard wavelength-discharge scalings and width-discharge scalings). Taken together with evidence from a marker bed for discharge variability at $\sim$10m stratigraphic distances, the variation in the scale of river deposits indicates that bank-forming discharge varied at both 10m stratigraphic (10$^2$ - 10$^6$ yr) and $\sim$100 m stratigraphic (10$^3$ - 10$^9$ yr) scales. Because these variations are correlated across the basin, they record a change in basin-scale forcing, rather than smaller-scale internal feedbacks. Changing sediment input leading to a change in characteristic slopes and/or drainage area could be responsible, and another possibility is changing climate ($\pm$50 W/m$^2$ in peak energy available for snow/ice melt).

The 163 comets observed during the WISE/NEOWISE prime mission represent the largest infrared survey to date of comets, providing constraints on dust, nucleus sizes, and CO+CO2 production. We present detailed analyses of the WISE/NEOWISE comet discoveries, and discuss observations of the active comets showing 4.6 $\mu$m band excess. We find a possible relation between dust and CO+CO2 production, as well as possible differences in the sizes of long and short period comet nuclei.

Titan's equatorial regions are covered by eastward propagating linear dunes. This direction is opposite to mean surface winds simulated by Global Climate Models (GCMs), which are oriented westward at these latitudes, similar to trade winds on Earth. Different hypotheses have been proposed to address this apparent contradiction, involving Saturn's gravitational tides, large scale topography or wind statistics, but none of them can explain a global eastward dune propagation in the equatorial band. Here we analyse the impact of equinoctial tropical methane storms developing in the superrotating atmosphere (i.e. the eastward winds at high altitude) on Titan's dune orientation. Using mesoscale simulations of convective methane clouds with a GCM wind profile featuring superrotation, we show that Titan's storms should produce fast eastward gust fronts above the surface. Such gusts dominate the aeolian transport, allowing dunes to extend eastward. This analysis therefore suggests a coupling between superrotation, tropical methane storms and dune formation on Titan. Furthermore, together with GCM predictions and analogies to some terrestrial dune fields, this work provides a general framework explaining several major features of Titan's dunes: linear shape, eastward propagation and poleward divergence, and implies an equatorial origin of Titan's dune sand.

The Wide-field Infrared Survey Explorer (WISE) observed 52 Centaurs and Scattered Disk Objects in the thermal infrared, including the 15 discoveries that were new. We present analyses of these observations to estimate sizes and mean optical albedos. We find mean albedos of 0.08 +/- 0.04 for the entire data set. Thermal fits yield average beaming parameters of 0.9 +/- 0.2 that are similar for both SDO and Centaur sub-classes. Biased cumulative size distributions yield size-frequency distribution power law indices ~ -1.7 +/- 0.3. The data also reveal a relation between albedo and color at the 3-sigma level. No significant relation between diameter and albedos is found.

The impactor flux early in Mars history was much higher than today, so sedimentary sequences include many buried craters. In combination with models for the impactor flux, observations of the number of buried craters can constrain sedimentation rates. Using the frequency of crater-river interactions, we find net sedimentation rate \lesssim 20-300 {\mu}m/yr at Aeolis Dorsa. This sets a lower bound of 1-15 Myr on the total interval spanned by fluvial activity around the Noachian-Hesperian transition. We predict that Gale Crater's mound (Aeolis Mons) took at least 10-100 Myr to accumulate, which is testable by the Mars Science Laboratory.