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.

MASER (Measurements, Analysis, and Simulation of Emission in the Radio range) is a comprehensive infrastructure dedicated to time-dependent low frequency radio astronomy (up to about 50 MHz). The main radio sources observed in this spectral range are the Sun, the magnetized planets (Earth, Jupiter, Saturn), and our Galaxy, which are observed either from ground or space. Ground observatories can capture high resolution data streams with a high sensitivity. Conversely, space-borne instruments can observe below the ionospheric cut-off (at about 10 MHz) and can be placed closer to the studied object. Several tools have been developed in the last decade for sharing space physics data. Data visualization tools developed by various institutes are available to share, display and analyse space physics time series and spectrograms. The MASER team has selected a sub-set of those tools and applied them to low frequency radio astronomy. MASER also includes a Python software library for reading raw data from agency archives.

The Europlanet-2020 programme, which ended on Aug 31st, 2019, included an activity called VESPA (Virtual European Solar and Planetary Access), which focused on adapting Virtual Observatory (VO) techniques to handle Planetary Science data. This paper describes some aspects of VESPA at the end of this 4-years development phase and at the onset of the newly selected Europlanet-2024 programme starting in 2020. The main objectives of VESPA are to facilitate searches both in big archives and in small databases, to enable data analysis by providing simple data access and online visualization functions, and to allow research teams to publish derived data in an interoperable environment as easily as possible. VESPA encompasses a wide scope, including surfaces, atmospheres, magnetospheres and planetary plasmas, small bodies, helio-physics, exoplanets, and spectroscopy in solid phase. This system relies in particular on standards and tools developed for the Astronomy VO (IVOA) and extends them where required to handle specificities of Solar System studies. It also aims at making the VO compatible with tools and protocols developed in different contexts, for instance GIS for planetary surfaces, or time series tools for plasma-related measurements. An essential part of the activity is to publish a significant amount of high-quality data in this system, with a focus on derived products resulting from data analysis or simulations.

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.

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.

The VESPA data access system focuses on applying Virtual Observatory (VO) standards and tools to Planetary Science. Building on a previous EC-funded Europlanet program, it has reached maturity during the first year of a new Europlanet 2020 program (started in 2015 for 4 years). The infrastructure has been upgraded to handle many fields of Solar System studies, with a focus both on users and data providers. This paper describes the broad lines of the current VESPA infrastructure as seen by a potential user, and provides examples of real use cases in several thematic areas. These use cases are also intended to identify hints for future developments and adaptations of VO tools to Planetary Science.

Since 2003 the Planetary Virtual Observatory and Laboratory (PVOL) has been storing and serving publicly through its web site a large database of amateur observations of the Giant Planets (Hueso et al., 2010a). These images are used for scientific research of the atmospheric dynamics and cloud structure on these planets and constitute a powerful resource to address time changing phenomena in their atmospheres. Advances over the last decade in observation techniques, and a wider conscience by professional astronomers of the quality of amateur observations, have resulted in the necessity to upgrade this database. We here present major advances in the PVOL database that has evolved into a full virtual planetary observatory encompassing also observations of Mercury, Venus, Mars, the Moon and the Galilean satellites. Besides the new objects, the images can be tagged and the database allows simple and complex searches over the data. The new web service: PVOL2 is available online in http://pvol2.ehu.eus/ , contains a fully functional search engine and constitutes one of the many services included in VESPA (Virtual Europan Solar and Planetary Access). Data from PVOL2 can be served from the VESPA portal using the EPN-TAP protocol. PVOL2 also provides long-term storage to amateur observations containing about 30,000 amateur observations starting in the year 2000. Current and past observations from the amateur community provide a global view of the Solar System planets over the years with several possibilities for scientific analysis and amateur astronomers involvement in planetary science.

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.

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.

A Data Access Protocol has been set up to search and retrieve Planetary Science data in general. This protocol will allow the user to select a subset of data from an archive in a standard way, based on the IVOA Table Access Protocol (TAP). The TAP mechanism is completed by an underlying Data Model and reference dictionaries. This paper describes the principle of the EPN- TAP protocol and interfaces, underlines the choices that have been made, and discusses possible evolutions.

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.

We acquired and analyzed adaptive-optics imaging observations of asteroid (2) Pallas from Keck II and the Very Large Telescope taken during four Pallas oppositions between 2003 and 2007, with spatial resolution spanning 32-88 km (image scales 13-20 km/pix). We improve our determination of the size, shape, and pole by a novel method that combines our AO data with 51 visual light-curves spanning 34 years of observations as well as occultation data. The shape model of Pallas derived here reproduces well both the projected shape of Pallas on the sky and light-curve behavior at all the epochs considered. We resolved the pole ambiguity and found the spin-vector coordinates to be within 5 deg. of [long, lat] = [30 deg., -16 deg.] in the ECJ2000.0 reference frame, indicating a high obliquity of ~84 deg., leading to high seasonal contrast. The best triaxial-ellipsoid fit returns radii of a=275 km, b= 258 km, and c= 238 km. From the mass of Pallas determined by gravitational perturbation on other minor bodies [(1.2 +/- 0.3) x 10-10 Solar Masses], we derive a density of 3.4 +/- 0.9 g.cm-3 significantly different from the density of C-type (1) Ceres of 2.2 +/- 0.1 g.cm-3. Considering the spectral similarities of Pallas and Ceres at visible and near-infrared wavelengths, this may point to fundamental differences in the interior composition or structure of these two bodies. We define a planetocentric longitude system for Pallas, following IAU guidelines. We also present the first albedo maps of Pallas covering ~80% of the surface in K-band. These maps reveal features with diameters in the 70-180 km range and an albedo contrast of about 6% wrt the mean surface albedo.