On 27 Apr 2015, when 67P/C-G was at 1.76 au from the Sun and moving towards perihelion, the OSIRIS and VIRTIS-M instruments on Rosetta observed the evolving dust and gas coma during a complete rotation of the comet. We aim to characterize the dust, H2O and CO2 gas spatial distribution in the inner coma. To do this we performed a quantitative analysis of the release of dust and gas and compared the observed H2O production rate with the one calculated using a thermo-physical model. For this study we selected OSIRIS WAC images at 612 nm (dust) and VIRTIS-M image cubes at 612 nm, 2700 nm (H2O) and 4200 nm (CO2). We measured the average signal in a circular annulus, to study spatial variation around the comet, and in a sector of the annulus, to study temporal variation in the sunward direction with comet rotation, both at a fixed distance of 3.1 km from the comet centre. The spatial correlation between dust and water, both coming from the sun-lit side of the comet, shows that water is the main driver of dust activity in this time period. The spatial distribution of CO2 is not correlated with water and dust. There is no strong temporal correlation between the dust brightness and water production rate as the comet rotates. The dust brightness shows a peak at 0deg sub-solar longitude, which is not pronounced in the water production. At the same epoch, there is also a maximum in CO2 production. An excess of measured water production, with respect to the value calculated using a simple thermo-physical model, is observed when the head lobe and regions of the Southern hemisphere with strong seasonal variations are illuminated. A drastic decrease in dust production, when the water production (both measured and from the model) displays a maximum, happens when typical Northern consolidated regions are illuminated and the Southern hemisphere regions with strong seasonal variations are instead in shadow.

We analyze 2-5 micrometre spectroscopic observations of the dust coma of comet 67P/Churyumov-Gerasimenko obtained with the VIRTIS-H instrument onboard Rosetta from 3 June to 29 October 2015 at heliocentric distances r_h = 1.24-1.55 AU. The 2-2.5 micrometre color, bolometric albedo, and color temperature are measured using spectral fitting. Data obtained at alpha = 90{\deg} solar phase angle show an increase of the bolometric albedo (0.05 to 0.14) with increasing altitude (0.5 to 8 km), accompanied by a possible marginal decrease of the color and color temperature. Possible explanations include the presence in the inner coma of dark particles on ballistic trajectories, and radial changes in particle composition. In the phase angle range 50-120{\deg}, phase reddening is significant (0.031 %/100 nm/{\deg}), for a mean color of 2 %/100 nm at alpha = 90{\deg}, that can be related to the roughness of the dust particles. Moreover, a decrease of the color temperature with decreasing phase angle is also observed at a rate of ~ 0.3 K/{\deg}, consistent with the presence of large porous particles, with low thermal inertia, and showing a significant day-to-night temperature contrast. Comparing data acquired at fixed phase angle (alpha = 90{\deg}), a 20% increase of the bolometric albedo is observed near perihelion. Heliocentric variations of the dust color are not significant in the analyzed time period. Measured color temperatures are varying from 260 to 320 K, and follow a r^0.6 variation in the r_h = 1.24-1.5 AU range, close to the expected r_h^0.5 value.

The latest developments of planetary exploration missions devoted to minor bodies required new solutions to correctly visualize and analyse data acquired over irregularly shaped bodies. ASI Space Science Data Center (SSDC-ASI, formerly ASDC-ASI Science Data Center) worked on this task since early 2013, when started developing the web tool MATISSE (Multi-purpose Advanced Tool for the Instruments of the Solar System Exploration) mainly focused on the Rosetta/ESA space mission data. In order to visualize very high-resolution shape models, MATISSE uses a Python module (vtpMaker), which can also be launched as a stand-alone command-line software. MATISSE and vtpMaker are part of the SSDC contribution to the new challenges imposed by the "orbital exploration" of minor bodies: 1) MATISSE allows to search for specific observations inside datasets and then analyse them in parallel, providing high-level outputs; 2) the 3D capabilities of both tools are critical in inferring information otherwise difficult to retrieve for non-spherical targets and, as in the case for the GIADA instrument onboard Rosetta, to visualize data related to the coma. New tasks and features adding valuable capabilities to the minor bodies SSDC tools are planned for the near future thanks to new collaborations.

We describe Castalia, a proposed mission to rendezvous with a Main Belt Comet (MBC), 133P/Elst-Pizarro. MBCs are a recently discovered population of apparently icy bodies within the main asteroid belt between Mars and Jupiter, which may represent the remnants of the population which supplied the early Earth with water. Castalia will perform the first exploration of this population by characterising 133P in detail, solving the puzzle of the MBC's activity, and making the first in situ measurements of water in the asteroid belt. In many ways a successor to ESA's highly successful Rosetta mission, Castalia will allow direct comparison between very different classes of comet, including measuring critical isotope ratios, plasma and dust properties. It will also feature the first radar system to visit a minor body, mapping the ice in the interior. Castalia was proposed, in slightly different versions, to the ESA M4 and M5 calls within the Cosmic Vision programme. We describe the science motivation for the mission, the measurements required to achieve the scientific goals, and the proposed instrument payload and spacecraft to achieve these.

We present 2-5 $\mu$m spectroscopic observations of the dust coma of 67P/Churyumov-Gerasimenko obtained with the VIRTIS-H instrument onboard Rosetta during two outbursts that occurred on 2015, 13 September 13.6 h UT and 14 September 18.8 h UT at 1.3 AU from the Sun. Scattering and thermal properties measured before the outburst are in the mean of values measured for moderately active comets. The colour temperature excess (or superheat factor) can be attributed to submicrometre-sized particles composed of absorbing material or to porous fractal-like aggregates such as those collected by the Rosetta in situ dust instruments. The power law index of the dust size distribution is in the range 2-3. The sudden increase of infrared emission associated to the outbursts is correlated with a large increase of the colour temperature (from 300 K to up to 630 K) and a change of the dust colour at 2-2.5 $\mu$m from red to blue colours, revealing the presence of very small grains ($\leq$ 100 nm) in the outburst material. In addition, the measured large bolometric albedos ($\sim$ 0.7) indicate bright grains in the ejecta, which could either be silicatic grains, implying the thermal degradation of the carbonaceous material, or icy grains. The 3-$\mu$m absorption band from water ice is not detected in the spectra acquired during the outbursts, whereas signatures of organic compounds near 3.4 $\mu$m are observed in emission. The H$_2$O 2.7-$\mu$m and CO$_2$ 4.3-$\mu$m vibrational bands do not show any enhancement during the outbursts.

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.

Infrared observations of the coma of 67P/Churyumov-Gerasimenko were carried out from July to September 2015, i.e., around perihelion (13 August 2015), with the high-resolution channel of the VIRTIS instrument onboard Rosetta. We present the analysis of fluorescence emission lines of H$_2$O, CO$_2$, $^{13}$CO$_2$, OCS, and CH$_4$ detected in limb sounding with the field of view at 2.7-5 km from the comet centre. Measurements are sampling outgassing from the illuminated southern hemisphere, as revealed by H$_2$O and CO$_2$ raster maps, which show anisotropic distributions, aligned along the projected rotation axis. An abrupt increase of water production is observed six days after perihelion. In the mean time, CO$_2$, CH$_4$, and OCS abundances relative to water increased by a factor of 2 to reach mean values of 32%, 0.47%, and 0.18%, respectively, averaging post-perihelion data. We interpret these changes as resulting from the erosion of volatile-poor surface layers. Sustained dust ablation due to the sublimation of water ice maintained volatile-rich layers near the surface until at least the end of the considered period, as expected for low thermal inertia surface layers. The large abundance measured for CO$_2$ should be representative of the 67P nucleus original composition, and indicates that 67P is a CO$_2$-rich comet. Comparison with abundance ratios measured in the northern hemisphere shows that seasons play an important role in comet outgassing. The low CO$_2$/H$_2$O values measured above the illuminated northern hemisphere are not original, but the result of the devolatilization of the uppermost layers.

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.

The increasing number and complexity of planetary exploration space missions require new tools to access, visualize and analyse data to improve their scientific return. ASI Science Data Center (ASDC) addresses this request with the web-tool MATISSE (Multi-purpose Advanced Tool for the Instruments of the Solar System Exploration), allowing the visualization of single observation or real-time computed high-order products, directly projected on the three-dimensional model of the selected target body. Using MATISSE it will be no longer needed to download huge quantity of data or to write down a specific code for every instrument analysed, greatly encouraging studies based on joint analysis of different datasets. In addition the extremely high-resolution output, to be used offline with a Python-based free software, together with the files to be read with specific GIS software, makes it a valuable tool to further process the data at the best spatial accuracy available. MATISSE modular structure permits addition of new missions or tasks and, thanks to dedicated future developments, it would be possible to make it compliant to the Planetary Virtual Observatory standards currently under definition. In this context the recent development of an interface to the NASA ODE REST API by which it is possible to access to public repositories is set.

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.

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.

The field of asteroid thermophysical modeling has experienced an extraordinary growth in the last ten years, as new thermal infrared data became available for hundreds of thousands of asteroids. The infrared emission of asteroids depends on the body's size, shape, albedo, thermal inertia, roughness and rotational properties. These parameters can therefore be derived by thermophysical modeling of infrared data. Thermophysical modeling led to asteroid size estimates that were confirmed at the few-percent level by later spacecraft visits. We discuss how instrumentation advances now allow mid-infrared interferometric observations as well as high-accuracy spectro-photometry, posing their own set of thermal-modeling challenges.We present major breakthroughs achieved in studies of the thermal inertia, a sensitive indicator for the nature of asteroids soils, allowing us, for instance, to determine the grain size of asteroidal regoliths. Thermal inertia also governs non-gravitational effects on asteroid orbits, requiring thermophysical modeling for precise asteroid dynamical studies. The radiative heating of asteroids, meteoroids, and comets from the Sun also governs the thermal stress in surface material; only recently has it been recognized as a significant weathering process. Asteroid space missions with thermal infrared instruments are currently undergoing study at all major space agencies. This will require a high level of sophistication of thermophysical models in order to analyze high-quality spacecraft data.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

In the framework of the Europlanet-RI program, a prototype of Virtual Observatory dedicated to Planetary Science was defined. Most of the activity was dedicated to the elaboration of standards to retrieve and visualize data in this field, and to provide light procedures to teams who wish to contribute with on-line data services. The architecture of this VO system and selected solutions are presented here, together with existing demonstrators.

Main Belt comets represent a recently discovered class of objects. They are quite intriguing because, while having a Tisserand invariant value higher than 3, are showing cometary activity. We study the activity of the Main Belt comets making the assumption that they are icy-bodies and that the activity has been triggered by an impact. We determine the characteristics of this activity and if the nowadays impact rate in the Main Asteroid Belt is compatible with the hypothesis of an activity triggered by a recent impact. Due to the fact that the Main Belt comets can be considered as a kind of comets, we apply a thermal evolution model developed for icy bodies in order to simulate their activity. We also apply a model to derive the impact rate, with respect to the size of the impactor, in the Main Belt. We demonstrate that a stable activity can result from a recent impact, able to expose ice-rich layers, and that the impact rate in the Main Belt is compatible with this explanation.