Comets hold the key to the understanding of our solar system, its formation and its evolution, and to the fundamental plasma processes at work both in it and beyond it. A comet nucleus emits gas as it is heated by the sunlight. The gas forms the coma, where it is ionised, becomes a plasma and eventually interacts with the solar wind. Besides these neutral and ionised gases, the coma also contains dust grains, released from the comet nucleus. As a cometary atmosphere develops when the comet travels through the solar system, large-scale structures, such as the plasma boundaries, develop and disappear, while at planets such large-scale structures are only accessible in their fully grown, quasi-steady state. In situ measurements at comets enable us to learn both how such large-scale structures are formed or reformed and how small-scale processes in the plasma affect the formation and properties of these large scale structures. Furthermore, a comet goes through a wide range of parameter regimes during its life cycle, where either collisional processes, involving neutrals and charged particles, or collisionless processes are at play, and might even compete in complicated transitional regimes. Thus a comet presents a unique opportunity to study this parameter space, from an asteroid-like to a Mars- and Venus-like interaction. Fast flybys of comets have made many new discoveries, setting the stage for a multi-spacecraft mission to accompany a comet on its journey through the solar system. This white paper reviews the present-day knowledge of cometary plasmas, discusses the many questions that remain unanswered, and outlines a multi-spacecraft ESA mission to accompany a comet that will answer these questions by combining both multi-spacecraft observations and a rendezvous mission, and at the same time advance our understanding of fundamental plasma physics and its role in planetary systems.

We use measurements from the Rosetta plasma consortium (RPC) Langmuir probe (LAP) and mutual impedance probe (MIP) to study the spatial distribution of low-energy plasma in the near-nucleus coma of comet 67P/Churyumov-Gerasimenko. The spatial distribution is highly structured with the highest density in the summer hemisphere and above the region connecting the two main lobes of the comet, i.e. the neck region. There is a clear correlation with the neutral density and the plasma to neutral density ratio is found to be about 1-2x10^-6, at a cometocentric distance of 10 km and at 3.1 AU from the sun. A clear 6.2 h modulation of the plasma is seen as the neck is exposed twice per rotation. The electron density of the collisonless plasma within 260 km from the nucleus falls of with radial distance as about 1/r. The spatial structure indicates that local ionization of neutral gas is the dominant source of low-energy plasma around the comet.

Magnetic reconnection is a fundamental process in solar system and astrophysical plasmas, through which stored magnetic energy associated with current sheets is converted into thermal, kinetic and wave energy. Magnetic reconnection is also thought to be a key process involved in shedding internally produced plasma from the giant magnetospheres at Jupiter and Saturn through topological reconfiguration of the magnetic field. The region where magnetic fields reconnect is known as the diffusion region and in this letter we report on the first encounter of the Cassini spacecraft with a diffusion region in Saturn's magnetotail. The data also show evidence of magnetic reconnection over a period of 19 h revealing that reconnection can, in fact, act for prolonged intervals in a rapidly rotating magnetosphere. We show that reconnection can be a significant pathway for internal plasma loss at Saturn. This counters the view of reconnection as a transient method of internal plasma loss at Saturn. These results, although directly relating to the magnetosphere of Saturn, have applications in the understanding of other rapidly rotating magnetospheres, including that of Jupiter and other astrophysical bodies.

Current disruption (CD) and the related kinetic instabilities in the near-Earth magnetosphere represent physical mechanisms which can trigger multi-scale substorm activity including global reorganizations of the magnetosphere. Lui et al. (2008) proposed a CD scenario in which the kinetic scale linear modes grow and reach the typical dipolarization scales through an inverse cascade. The experimental verification of the inverse nonlinear cascade is based on wavelet analysis. In this paper the Hilbert-Huang transform is used which is suitable for nonlinear systems and allows to reconstruct the time-frequency representation of empirical decomposed modes in an adaptive manner. It was found that, in the Lui et al. (2008) event, the modes evolve globally from high-frequencies to low-frequencies. However, there are also local frequency evolution trends oriented towards high-frequencies, indicating that the underlying processes involve multi-scale physics and non-stationary fluctuations for which the simple inverse cascade scenario is not correct.

Recent research has shown that distinct physical regions in the Venusian induced magnetosphere are recognizable from the variations of strength of the magnetic field and its wave/fluctuation activity. In this paper the statistical properties of magnetic fluctuations are investigated in the Venusian magnetosheath and wake regions. The main goal is to identify the characteristic scaling features of fluctuations along Venus Express (VEX) trajectory and to understand the specific circumstances of the occurrence of different types of scalings. For the latter task we also use the results of measurements from the previous missions to Venus. Our main result is that the changing character of physical interactions between the solar wind and the planetary obstacle is leading to different types of spectral scaling in the near-Venusian space. Noisy fluctuations are observed in the magnetosheath, wavy structures near the terminator and in the nightside near-planet wake. Multi-scale turbulence is observed at the magnetosheath boundary layer and near the quasi-parallel bow shock. Magnetosheath boundary layer turbulence is associated with an average magnetic field which is nearly aligned with the Sun-Venus line. Noisy magnetic fluctuations are well described with the Gaussian statistics. Both magnetosheath boundary layer and near shock turbulence statistics exhibit non-Gaussian features and intermittency over small spatio-temporal scales. The occurrence of turbulence near magnetosheath boundaries can be responsible for the local heating of plasma observed by previous missions.

Recent research has shown that distinct physical regions in the Venusian induced magnetosphere are recognizable from the variations of strength and of wave/fluctuation activity of the magnetic field. In this paper the statistical properties of magnetic fluctuations are investigated in the Venusian magnetosheath, terminator, and wake regions. The latter two regions were not visited by previous missions. We found 1/f fluctuations in the magnetosheath, large-scale structures near the terminator and more developed turbulence further downstream in the wake. Location independent short-tailed non-Gaussian statistics was observed.