Before Rosetta, the space missions Giotto and Stardust shaped our view on cometary dust, supported by plentiful data from Earth based observations and interplanetary dust particles collected in the Earth's atmosphere. The Rosetta mission at comet 67P/Churyumov-Gerasimenko was equipped with a multitude of instruments designed to study cometary dust. While an abundant amount of data was presented in several individual papers, many focused on a dedicated measurement or topic. Different instruments, methods, and data sources provide different measurement parameters and potentially introduce different biases. This can be an advantage if the complementary aspect of such a complex data set can be exploited. However, it also poses a challenge in the comparison of results in the first place. The aim of this work therefore is to summarise dust results from Rosetta and before. We establish a simple classification as a common framework for inter-comparison. This classification is based on a dust particle's structure, porosity, and strength as well as its size. Depending on the instrumentation, these are not direct measurement parameters but we chose them as they were the most reliable to derive our model. The proposed classification already proved helpful in the Rosetta dust community and we propose to take it into consideration also beyond. In this manner we hope to better identify synergies between different instruments and methods in the future.

The phase function of the dust coma of comet 67P has been determined from Rosetta/OSIRIS images \citep{Bertini17}. This function show a deep minimum at phase angles near 100$^\circ$, and a strong backscattering enhancement. These two properties cannot be reproduced by regular models of cometary dust, most of them based on wavelength-sized and randomly-oriented aggregate particles. We show, however, that an ensamble of oriented elongated particles of a wide variety of aspect ratios, with radii $r \gtrsim$10 $\mu$m, and whose long axes are perpendicular to the direction of the solar radiation, are capable of reproducing the observed phase function. These particles must be absorbing, with an imaginary part of the refractive index of about 0.1 to match the expected geometric albedo, and with porosity in the 60-70\% range.

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.

On 3 July 2016, several instruments on board ESA's Rosetta spacecraft detected signs of an outburst event on comet 67P, at a heliocentric distance of 3.32 AU from the sun, outbound from perihelion. We here report on the inferred properties of the ejected dust and the surface change at the site of the outburst. The activity coincided with the local sunrise and continued over a time interval of 14 - 68 minutes. It left a 10m-sized icy patch on the surface. The ejected material comprised refractory grains of several hundred microns in size, and sub-micron-sized water ice grains. The high dust mass production rate is incompatible with the free sublimation of crystalline water ice under solar illumination as the only acceleration process. Additional energy stored near the surface must have increased the gas density. We suggest a pressurized sub-surface gas reservoir, or the crystallization of amorphous water ice as possible causes.

The processes that led to the formation of the planetary bodies in the Solar System are still not fully understood. Using the results obtained with the comprehensive suite of instruments on-board ESA's Rosetta mission, we present evidence that comet 67P/Churyumov-Gerasimenko likely formed through the gentle gravitational collapse of a bound clump of mm-sized dust aggregates ("pebbles"), intermixed with microscopic ice particles. This formation scenario leads to a cometary make-up that is simultaneously compatible with the global porosity, homogeneity, tensile strength, thermal inertia, vertical temperature profiles, sizes and porosities of emitted dust, and the steep increase in water-vapour production rate with decreasing heliocentric distance, measured by the instruments on-board the Rosetta spacecraft and the Philae lander. Our findings suggest that the pebbles observed to be abundant in protoplanetary discs around young stars provide the building material for comets and other minor bodies.

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.

The Philae lander of the Rosetta mission, aimed at the in situ investigation of comet 67P/C-G, was deployed to the surface of the comet nucleus on 12 Nov 2014 at 2.99 AU heliocentric distance. The Dust Impact Monitor (DIM) as part of the Surface Electric Sounding and Acoustic Monitoring Experiment (SESAME) on the lander employed piezoelectric detectors to detect the submillimetre- and millimetre-sized dust and ice particles emitted from the nucleus. We determine the upper limit of the ambient flux of particles in the measurement range of DIM based on the measurements performed with the instrument during Philae's descent to its nominal landing site Agilkia at distances of about 22 km, 18 km, and 5 km from the nucleus barycentre and at the final landing site Abydos. The geometric factor of the DIM sensor is calculated assuming an isotropic ambient flux of the submillimetre- and millimetre-sized particles. For the measurement intervals when no particles were detected the maximum true impact rate was calculated by assuming Poisson distribution of the impacts, and it was given as the detection limit at a 95% confidence level. Based on measurements performed with DIM, the upper limit of the flux of particles in the measurement range of the instrument was of the order of $10^{-8}-10^{-7}m^{-2}s^{-1}sr^{-1}$ during descent and $1.6\cdot10^{-9}m^{-2}s^{-1}sr^{-1}$ at Abydos on 13 and 14 Nov 2014. Considering particle speeds below escape velocity, the upper limit for the volume density of particles in the measurement range of DIM was constrained to $10^{-11}m^{-3}-10^{-12}m^{-3}$. Results of the calculations performed with the GIPSI tool on the expected particle fluxes during the descent of Philae were compatible with the non-detection of compact particles by the DIM instrument.

We present a simple model for gas and dust flow from 67P/Churyumov-Gerasimenko that can be used to understand the grain impact observed by the DIM instrument on Philae (Krueger et al. 2015) We show how model results when applied to the GIADA measurements (Rotundi et al. 2015; Della Corte et al. 2015) can be used, in conjunction with the results found by the MIRO (Schloerb et al. 2015) and VIRTIS (De Sanctis et al. 2015) instruments to infer surface properties such as surface temperature and surface ice fraction.

The ESA Rosetta spacecraft, currently orbiting around comet 67P, has already provided in situ measurements of the dust grain properties from several instruments, particularly OSIRIS and GIADA. We propose adding value to those measurements by combining them with ground-based observations of the dust tail to monitor the overall, time-dependent dust-production rate and size distribution. To constrain the dust grain properties, we take Rosetta OSIRIS and GIADA results into account, and combine OSIRIS data during the approach phase (from late April to early June 2014) with a large data set of ground-based images that were acquired with the ESO Very Large Telescope (VLT) from February to November 2014. A Monte Carlo dust tail code has been applied to retrieve the dust parameters. Key properties of the grains (density, velocity, and size distribution) were obtained from Rosetta observations: these parameters were used as input of the code to considerably reduce the number of free parameters. In this way, the overall dust mass-loss rate and its dependence on the heliocentric distance could be obtained accurately. The dust parameters derived from the inner coma measurements by OSIRIS and GIADA and from distant imaging using VLT data are consistent, except for the power index of the size-distribution function, which is $\alpha$=--3, instead of $\alpha$=--2, for grains smaller than 1 mm. This is possibly linked to the presence of fluffy aggregates in the coma. The onset of cometary activity occurs at approximately 4.3 au, with a dust production rate of 0.5 kg/s, increasing up to 15 kg/s at 2.9 au. This implies a dust-to-gas mass ratio varying between 3.8 and 6.5 for the best-fit model when combined with water-production rates from the MIRO experiment.