Uranus and Neptune, and their diverse satellite and ring systems, represent the least explored environments of our Solar System, and yet may provide the archetype for the most common outcome of planetary formation throughout our galaxy. Ice Giants will be the last remaining class of Solar System planet to have a dedicated orbital explorer, and international efforts are under way to realise such an ambitious mission in the coming decades. In 2019, the European Space Agency released a call for scientific themes for its strategic science planning process for the 2030s and 2040s, known as Voyage 2050. We used this opportunity to review our present-day knowledge of the Uranus and Neptune systems, producing a revised and updated set of scientific questions and motivations for their exploration. This review article describes how such a mission could explore their origins, ice-rich interiors, dynamic atmospheres, unique magnetospheres, and myriad icy satellites, to address questions at the heart of modern planetary science. These two worlds are superb examples of how planets with shared origins can exhibit remarkably different evolutionary paths: Neptune as the archetype for Ice Giants, whereas Uranus may be atypical. Exploring Uranus' natural satellites and Neptune's captured moon Triton could reveal how Ocean Worlds form and remain active, redefining the extent of the habitable zone in our Solar System. For these reasons and more, we advocate that an Ice Giant System explorer should become a strategic cornerstone mission within ESA's Voyage 2050 programme, in partnership with international collaborators, and targeting launch opportunities in the early 2030s.

The goal of Project GAUSS is to return samples from the dwarf planet Ceres. Ceres is the most accessible ocean world candidate and the largest reservoir of water in the inner solar system. It shows active cryovolcanism and hydrothermal activities in recent history that resulted in minerals not found in any other planets to date except for Earth's upper crust. The possible occurrence of recent subsurface ocean on Ceres and the complex geochemistry suggest possible past habitability and even the potential for ongoing habitability. Aiming to answer a broad spectrum of questions about the origin and evolution of Ceres and its potential habitability, GAUSS will return samples from this possible ocean world for the first time. The project will address the following top-level scientific questions: 1) What is the origin of Ceres and the origin and transfer of water and other volatiles in the inner solar system? 2) What are the physical properties and internal structure of Ceres? What do they tell us about the evolutionary and aqueous alteration history of icy dwarf planets? 3) What are the astrobiological implications of Ceres? Was it habitable in the past and is it still today? 4) What are the mineralogical connections between Ceres and our current collections of primitive meteorites? GAUSS will first perform a high-resolution global remote sensing investigation, characterizing the geophysical and geochemical properties of Ceres. Candidate sampling sites will then be identified, and observation campaigns will be run for an in-depth assessment of the candidate sites. Once the sampling site is selected, a lander will be deployed on the surface to collect samples and return them to Earth in cryogenic conditions that preserves the volatile and organic composition as well as the original physical status as much as possible.

The measured nitrogen-to-carbon ratio in comets is lower than for the Sun, a discrepancy which could be alleviated if there is an unknown reservoir of nitrogen in comets. The nucleus of comet 67P/Churyumov-Gerasimenko exhibits an unidentified broad spectral reflectance feature around 3.2 micrometers, which is ubiquitous across its surface. On the basis of laboratory experiments, we attribute this absorption band to ammonium salts mixed with dust on the surface. The depth of the band indicates that semivolatile ammonium salts are a substantial reservoir of nitrogen in the comet, potentially dominating over refractory organic matter and more volatile species. Similar absorption features appear in the spectra of some asteroids, implying a compositional link between asteroids, comets, and the parent interstellar cloud.

Data acquired at Ceres by the visible channel of the Visible and InfraRed mapping spectrometer (VIR) on board the NASA Dawn spacecraft are affected by the temperatures of both the visible (VIS) and the infrared (IR) sensors, which are respectively a CCD and a HgCdTe array. The variations of the visible channel temperatures measured during the sessions of acquisitions are correlated with variations in the spectral slope and shape for all the mission phases. The infrared channel (IR) temperature is more stable during the acquisitions, nonetheless it is characterized by a bi-modal distribution whether the cryocooler (and therefore the IR channel) is used or not during the visible channel operations. When the infrared channel temperature is high (175K, i.e. not in use and with crycooler off), an additional negative slope and a distortion are observed in the spectra of the visible channel. We developed an empirical correction based on a reference spectrum for the whole data set; it is designed to correct the two issues related to the sensor temperatures that we have identified. The reference spectrum is calculated to be representative of the global Ceres' surface. It is also made of data acquired when the visible and infrared channel temperatures are equal to the ones measured during an observation of the Arcturus star by VIR, which is consistent with several ground-based observations. The developed correction allows reliable analysis and mapping to be performed by minimizing the artifacts induced by fluctuations of the VIS temperature. Thanks to this correction, a direct comparison between different mission phases during which VIR experienced different visible and infrared channel temperatures is now possible.

The physical properties of cometary nuclei observed today relate to their complex history and help to constrain their formation and evolution. In this article, we review some of the main physical properties of cometary nuclei and focus in particular on the thermal, mechanical, structural and dielectric properties, emphasizing the progress made during the Rosetta mission. Comets have a low density of 480 $\pm$ 220 kg m-3 and a low permittivity of 1.9 - 2.0, consistent with a high porosity of 70 - 80 %, are weak with a very low global tensile strength $<$100 Pa, and have a low bulk thermal inertia of 0 - 60 J K-1 m-2 s-1/2 that allowed them to preserve highly volatiles species (e.g. CO, CO2, CH4, N2) into their interior since their formation. As revealed by 67P/Churyumov-Gerasimenko, the above physical properties vary across the nucleus, spatially at its surface but also with depth. The broad picture is that the bulk of the nucleus consists of a weakly bonded, rather homogeneous material that preserved primordial properties under a thin shell of processed material, and possibly covered by a granular material; this cover might in places reach a thickness of several meters. The properties of the top layer (the first meter) are not representative of that of the bulk nucleus. More globally, strong nucleus heterogeneities at a scale of a few meters are ruled out on 67P small lobe.

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.

Laboratory spectral measurements of relevant analogue materials were performed in the framework of the Rosetta mission in order to explain the surface spectral properties of comet 67P. Fine powders of coal, iron sulphides, silicates and their mixtures were prepared and their spectra measured in the Vis-IR range. These spectra are compared to a reference spectrum of 67P nucleus obtained with the VIRTIS/Rosetta instrument up to 2.7 {\mu}m, excluding the organics band centred at 3.2 {\mu}m. The species used are known to be chemical analogues for cometary materials which could be present at the surface of 67P. Grain sizes of the powders range from tens of nanometres to hundreds of micrometres. Some of the mixtures studied here actually reach the very low reflectance level observed by VIRTIS on 67P. The best match is provided by a mixture of sub-micron coal, pyrrhotite, and silicates. Grain sizes are in agreement with the sizes of the dust particles detected by the GIADA, MIDAS and COSIMA instruments on board Rosetta. The coal used in the experiment is responsible for the spectral slope in the visible and infrared ranges. Pyrrhotite, which is strongly absorbing, is responsible for the low albedo observed in the NIR. The darkest components dominate the spectra, especially within intimate mixtures. Depending on sample preparation, pyrrhotite can coat the coal and silicate aggregates. Such coating effects can affect the spectra as much as particle size. In contrast, silicates seem to play a minor role.

We study the spectrophotometric properties of dwarf planet Ceres in the VIS-IR spectral range by means of hyper-spectral images acquired by the VIR imaging spectrometer on board the NASA Dawn mission. Disk-resolved observations with a phase angle within the $7^{\circ}<\alpha<132^{\circ}$ interval were used to characterize Ceres' phase curve in the 0.465-4.05 $\mu$m spectral range. Hapke's model was applied to perform the photometric correction of the dataset, allowing us to produce albedo and color maps of the surface. The $V$-band magnitude phase function of Ceres was fitted with both the classical linear model and H-G formalism. The single-scattering albedo and the asymmetry parameter at 0.55$\mu$m are $w=0.14\pm0.02$ and $\xi=-0.11\pm0.08$, respectively (two-lobe Henyey-Greenstein phase function); the modeled geometric albedo is $0.094\pm0.007$; the roughness parameter is $\bar{\theta}=29^{\circ}\pm6^{\circ}$. Albedo maps indicate small variability on a global scale with an average reflectance of $0.034 \pm 0.003$. Isolated areas such as the Occator bright spots, Haulani, and Oxo show an albedo much higher than average. We measure a significant spectral phase reddening, and the average spectral slope of Ceres' surface after photometric correction is $1.1\%k\AA^{-1}$ and $0.85\%k\AA^{-1}$ at VIS and IR wavelengths, respectively. Broadband color indices are $V-R=0.38\pm0.01$ and $R-I=0.33\pm0.02$. H-G modeling of the $V$-band magnitude phase curve for $\alpha<30^{\circ}$ gives $H=3.14\pm0.04$ and $G=0.10\pm0.04$, while the classical linear model provides $V(1,1,0^{\circ})=3.48\pm0.03$ and $\beta=0.036\pm0.002$. The comparison with spectrophotometric properties of other minor bodies indicates that Ceres has a less back-scattering phase function and a slightly higher albedo than comets and C-type objects. However, the latter represents the closest match in the usual asteroid taxonomy.

Water ice-rich patches have been detected on the surface of comet 67P/Churyumov-Gerasimenko by the VIRTIS hyperspectral imager on-board the Rosetta spacecraft, since the orbital insertion in late August 2014. Among those, three icy patches have been selected, and VIRTIS data are used to analyse their properties and their temporal evolution while the comet was moving towards the Sun. We performed an extensive analysis of the spectral parameters, and we applied the Hapke radiative transfer model to retrieve the abundance and grain size of water ice, as well as the mixing modalities of water ice and dark terrains on the three selected water ice rich areas. Study of the spatial distribution of the spectral parameters within the ice-rich patches has revealed that water ice follows different patterns associated to a bimodal distribution of the grains: ~50 {\mu}m sized and ~2000 {\mu}m sized. In all three cases, after the first detections at about 3.5 AU heliocentric distance, the spatial extension and intensity of the water ice spectral features increased, it reached a maximum after 60-100 days at about 3.0 AU, and was followed by an approximately equally timed decrease and disappearanceat about ~2.2 AU, before perihelion. The behaviour of the analysed patches can be assimilated to a seasonal cycle. In addition we found evidence of short-term variability associated to a diurnal water cycle. The similar lifecycle of the three icy regions indicates that water ice is uniformly distributed in the subsurface layers, and no large water ice reservoirs are present.

Since the orbital insertion of the Rosetta spacecraft, comet 67P/Churyumov-Gerasimenko (67P/C-G) has been mapped by OSIRIS camera and VIRTIS spectro-imager, producing a huge quantity of images and spectra of the comet's nucleus. The aim of this work is to search for the presence of H$_2$O on the nucleus which, in general, appears very dark and rich in dehydrated organic material. After selecting images of the bright spots which could be good candidates to search for H$_2$O ice, taken at high resolution by OSIRIS, we check for spectral cubes of the selected coordinates to identify these spots observed by VIRTIS. The selected OSIRIS images were processed with the OSIRIS standard pipeline and corrected for the illumination conditions for each pixel using the Lommel-Seeliger disk law. The spots with higher I/F were selected and then analysed spectrophotometrically and compared with the surrounding area. We selected 13 spots as good targets to be analysed by VIRTIS to search for the 2 micron absorption band of water ice in the VIRTIS spectral cubes. Out of the 13 selected bright spots, eight of them present positive H$_2$O ice detection on the VIRTIS data. A spectral analysis was performed and the approximate temperature of each spot was computed. The H$_2$O ice content was confirmed by modeling the spectra with mixing (areal and intimate) of H$_2$O ice and dark terrain, using Hapke's radiative transfer modeling. We also present a detailed analysis of the detected spots.

From August to November 2014 the Rosetta orbiter has performed an extensive observation campaign aimed at the characterization of 67P/CG nucleus properties and to the selection of the Philae landing site. The campaign led to the production of a global map of the illuminated portion of 67P/CG nucleus. During this prelanding phase the comet's heliocentric distance decreased from 3.62 to 2.93 AU while Rosetta was orbiting around the nucleus at distances between 100 to 10 km. VIRTIS-M, the Visible and InfraRed Thermal Imaging Spectrometer - Mapping channel (Coradini et al. 2007) onboard the orbiter, has acquired 0.25-5.1 micron hyperspectral data of the entire illuminated surface, e.g. the north hemisphere and the equatorial regions, with spatial resolution between 2.5 and 25 m/pixel. I/F spectra have been corrected for thermal emission removal in the 3.5-5.1 micron range and for surface's photometric response. The resulting reflectance spectra have been used to compute several Cometary Spectral Indicators (CSI): single scattering albedo at 0.55 micron, 0.5-0.8 micron and 1.0-2.5 micron spectral slopes, 3.2 micron organic material and 2.0 micron water ice band parameters (center, depth) with the aim to map their spatial distribution on the surface and to study their temporal variability as the nucleus moved towards the Sun. Indeed, throughout the investigated period, the nucleus surface shows a significant increase of the single scattering albedo along with a decrease of the 0.5-0.8 and 1.0-2.5 micron spectral slopes, indicating a flattening of the reflectance. We attribute the origin of this effect to the partial removal of the dust layer caused by the increased contribution of water sublimation to the gaseous activity as comet crossed the frost-line.

Bright features have been recently discovered by Dawn on Ceres, which extend previous photometric and Space Telescope observations. These features should produce distortions of the line profiles of the reflected solar spectrum and therefore an apparent radial velocity variation modulated by the rotation of the dwarf planet. Here we report on two sequences of observations of Ceres performed in the nights of 31 July, 26-27 August 2015 by means of the high-precision HARPS spectrograph at the 3.6-m La Silla ESO telescope. The observations revealed a quite complex behaviour which likely combines a radial velocity modulation due to the rotation with an amplitude of approx +/- 6 m/s and an unexpected diurnal effect. The latter changes imply changes in the albedo of Occator's bright features due to the blaze produced by the exposure to solar radiation. The short-term variability of Ceres' albedo is on timescales ranging from hours to months and can both be confirmed and followed by means of dedicated radial velocity observations.

Earth and orbital-based radar observations of asteroids provide a unique opportunity to characterize surface roughness and the dielectric properties of their surfaces, as well as potentially explore some of their shallow subsurface physical properties. If the dielectric and topographic properties of asteroid's surfaces are defined, one can constrain their surface textural characteristics as well as potential subsurface volatile enrichment using the observed radar backscatter. To achieve this objective, we establish the first dielectric model of asteroid Vesta for the case of a dry, volatile-poor regolith--employing an analogy to the dielectric properties of lunar soil, and adjusted for the surface densities and temperatures deduced from Dawn's Visible and InfraRed mapping spectrometer (VIR). Our model suggests that the real part of the dielectric constant at the surface of Vesta is relatively constant, ranging from 2.3 to 2.5 from the night- to day-side of Vesta, while the loss tangent shows slight variation as a function of diurnal temperature, ranging from 0.006 to 0.008. We estimate the surface porosity to be ~55% in the upper meter of the regolith, as derived from VIR observations. This is ~12% higher than previous estimation of porosity derived from previous Earth-based X- and S-band radar observation. We suggest that the radar backscattering properties of asteroid Vesta will be mainly driven by the changes in surface roughness rather than potential dielectric variations in the upper regolith in the X- and S-band.

Over the last eight years, the Visual and Infrared Mapping Spectrometer (VIMS) aboard the Cassini orbiter has returned hyperspectral images in the 0.35-5.1 micron range of the icy satellites and rings of Saturn. These very different objects show significant variations in surface composition, roughness and regolith grain size as a result of their evolutionary histories, endogenic processes and interactions with exogenic particles. The distributions of surface water ice and chromophores, i.e. organic and non-icy materials, across the saturnian system, are traced using specific spectral indicators (spectral slopes and absorption band depths) obtained from rings mosaics and disk-integrated satellites observations by VIMS.

In the last few years Cassini-VIMS, the Visible and Infared Mapping Spectrometer, returned to us a comprehensive view of the Saturn's icy satellites and rings. After having analyzed the satellites' spectral properties (Filacchione et al. (2007a)) and their distribution across the satellites' hemispheres (Filacchione et al. (2010)), we proceed in this paper to investigate the radial variability of icy satellites (principal and minor) and main rings average spectral properties. This analysis is done by using 2,264 disk-integrated observations of the satellites and a 12x700 pixels-wide rings radial mosaic acquired with a spatial resolution of about 125 km/pixel. The comparative analysis of these data allows us to retrieve the amount of both water ice and red contaminant materials distributed across Saturn's system and the typical surface regolith grain sizes. These measurements highlight very striking differences in the population here analyzed, which vary from the almost uncontaminated and water ice-rich surfaces of Enceladus and Calypso to the metal/organic-rich and red surfaces of Iapetus' leading hemisphere and Phoebe. Rings spectra appear more red than the icy satellites in the visible range but show more intense 1.5-2.0 micron band depths. The correlations among spectral slopes, band depths, visual albedo and phase permit us to cluster the saturnian population in different spectral classes which are detected not only among the principal satellites and rings but among co-orbital minor moons as well. Finally, we have applied Hapke's theory to retrieve the best spectral fits to Saturn's inner regular satellites using the same methodology applied previously for Rhea data discussed in Ciarniello et al. (2011).

The origin of the irregular satellites of the giant planets has been long debated since their discovery. Their dynamical features argue against an in-situ formation suggesting they are captured bodies, yet there is no global consensus on the physical process at the basis of their capture. In this paper we explore the collisional capture scenario, where the actual satellites originated from impacts occurred within Saturn's influence sphere. By modeling the inverse capture problem, we estimated the families of orbits of the possible parent bodies and the specific impulse needed for their capture. The orbits of these putative parent bodies are compared to those of the minor bodies of the outer Solar System to outline their possible region of formation. Finally, we tested the collisional capture hypothesis on Phoebe by taking advantage of the data supplied by Cassini on its major crater, Jason. Our results presented a realistic range of solutions matching the observational and dynamical data.

Among the icy satellites of Saturn, Iapetus shows a striking dichotomy between its leading and trailing hemispheres, the former being significantly darker than the latter. Thanks to the VIMS imaging spectrometer on-board Cassini, it is now possible to investigate the spectral features of the satellites in Saturn system within a wider spectral range and with an enhanced accuracy than with previously available data. In this work, we present an application of the G-mode method to the high resolution, visible and near infrared data of Phoebe, Iapetus and Hyperion collected by Cassini/VIMS, to search for compositional correlations. We also present the results of a dynamical study on the efficiency of Iapetus in capturing dust grains travelling inward in Saturn system to evaluate the viability of Poynting-Robertson drag as the physical mechanism transferring the dark material to the satellite. The results of spectroscopic classification are used jointly with the ones of the dynamical study to describe a plausible physical scenario for the origin of Iapetus' dichotomy. Our work shows that mass transfer from the outer Saturnian system is an efficient mechanism, particularly for the range of sizes hypothesised for the particles composing the newly discovered outer ring around Saturn. Both spectral and dynamical data indicate Phoebe as the main source of the dark material. However, we suggest a multi-source scenario where now extinct prograde satellites and the disruptive impacts that generated the putative collisional families played a significant role in supplying the original amount of dark material.

In this paper we will summarize some of the most important results of the Cassini mission concerning the satellites of Saturn. Given the long duration of the mission, the complexity of the payload onboard the Cassini Orbiter and the amount of data gathered on the satellites of Saturn, it would be impossible to describe all the new discoveries made, therefore we will describe only some selected, paramount examples showing how Cassini's data confirmed and extended ground-based observations. In particular we will describe the achievements obtained for the satellites Phoebe, Enceladus and Titan. We will also put these examples in the perspective of the overall evolution of the system, stressing out why the selected satellites are representative of the overall evolution of the Saturn system.