At a low activity comet the plasma is distributed in an asymmetric way. The hybrid simulation code Amitis is used to look at the spatial evolution of ion velocity distribution functions (VDFs), from the upstream solar wind to within the comet magnetosphere where the solar wind is heavily mass-loaded by the cometary plasma. We find that the spatial structures of the ions and fields form a highly asymmetric, half-open induced magnetosphere. The VDFs of solar wind and cometary ions vary drastically for different locations in the comet magnetosphere. The shape of the VDFs differ for different species. The solar wind protons show high anisotropies that occasionally resemble partial rings, in particular at small cometocentric distances. A second, decoupled, proton population is also found. Solar wind alpha particles show similar anisotropies, although less pronounced and at different spatial scales. The VDFs of cometary ions are mostly determined by the structure of the electric field. We perform supplementary dynamic particle backtracing to understand the flow patterns of solar wind ions that lead to these anisotropic distributions. This tracing is needed to understand the origin of cometary ions in a given part of the comet magnetosphere. The particle tracing also aids in interpreting observed VDFs and relating them to spatial features in the electric and magnetic fields of the comet environment.

We apply a recently proposed method to estimate heavy ion escape from Mars. The method combines in situ observations with a hybrid plasma model, which treats ions as particles and electrons as a fluid. With this method, we investigate how solar upstream conditions, including solar extreme ultraviolet (EUV) radiation, solar wind dynamic pressure, and interplanetary magnetic field (IMF) strength and cone angle, affect the heavy ion loss. The results indicate that the heavy ion escape rate is greater in high EUV conditions. The escape rate increases with increasing solar wind dynamic pressure, and decreases as the IMF strength increases. The ion escape rate is highest when the solar wind is parallel to the IMF and lowest when they are perpendicular. The plume escape rate decreases when the solar wind convective electric field increases.

No spacecraft visiting a comet has been equipped with instruments to directly measure the static electric field. However, the electric field can occasionally be estimated indirectly by observing its effects on the ion velocity distribution. We present such observations made by the Rosetta spacecraft on 19th of April 2016 when comet 67P was at a low outgassing rate and the plasma environment was relatively homogeneous. The ion velocity distributions show the cometary ions on the first half of their gyration. We estimate the bulk drift velocity and the gyration speed from the distributions. By using the local measured magnetic field and assuming an E x B drift of the gyrocentre, we get an estimate for the average electric field driving this ion motion. We analyse a period of 13h, during which the plasma environment does not change drastically. We find that the average strength of the electric field is 0.21mV/m. The direction of the electric field is mostly anti-sunward. This is in agreement with previous results based on different methods

The present atmosphere of Venus contains almost no water, but recent measurements indicate that in its early history Venus had an Earth-like ocean. Understanding how the Venusian atmosphere evolved is important not only for Venus itself, but also for understanding the evolution of other planetary atmospheres. In this study, we quantify the escape rates of oxygen ions from the present Venus to infer the past of the Venusian atmosphere. We show that an extrapolation of the current escape rates back in time leads to the total escape of 0.02-0.6 m of a global equivalent layer of water. This implies that the loss of ions to space, inferred from the present state, cannot account for the loss of an historical Earth-like ocean. We find that the O+ escape rate increases with solar wind energy flux, where more energy available leads to a higher escape rate. Oppositely, the escape rate decrease slightly with increased EUV flux, though the small variation of EUV flux over the measured solar cycle may explain the weak dependency. These results indicate that there isn’t enough energy transferred from the solar wind to Venus’ upper atmosphere that can lead to the escape of the atmosphere over the past 3.9 billion years. This means that the Venusian atmosphere didn’t have as much water in its atmosphere as previously assumed or the present-day escape rates don’t represent the historical escape rates at Venus. Otherwise, some other mechanisms have acted to more effectively remove the water from the Venusian atmosphere.

We report the detection of large-amplitude, quasi-harmonic density fluctuations with associated magnetic field oscillations in the region surrounding the diamagnetic cavity of comet 67P. Typical frequencies are ~0.1 Hz, corresponding to ~10 times the water and ≲0.5 times the proton gyro-frequencies, respectively. Magnetic field oscillations are not always clearly observed in association with these density fluctuations, but when they are, they consistently have wave vectors perpendicular to the background magnetic field, with the principal axis of polarization close to field-aligned and with a ~90° phase shift w.r.t. the density fluctuations. The fluctuations are observed in association with asymmetric plasma density and magnetic field enhancements previously found in the region surrounding the diamagnetic cavity, occurring predominantly on their descending slopes. This is a new type of waves not previously observed at comets. They are likely Ion Bernstein waves, and we propose that they are excited by unstable ring, ring-beam or spherical shell distributions of cometary ions just outside the cavity boundary. These waves may play an important role in redistributing energy between different particle populations and reshape the plasma environment of the comet.

We investigate sunward planetary ions in the Martian magnetotail that potentially reduce the amount of escaping ions. The global properties of sunward flows in the Martian magnetotail are characterized, based on over 13-years of ion data (May 2007–December 2020) collected by the ASPERA-3 instrument on Mars Express. We find that sunward flows mainly occur in the vicinity of the crustal fields, implying that crustal fields may play a key role in producing such flows. The occurrence rate and sunward flux are higher during solar maximum rather than solar minimum. However, we identify a relatively low occurrence rate of sunward flows and low sunward flux, suggesting that sunward flows have negligible influence on total ion escape at Mars. This is different from those at Venus, where sunward flows can significantly decrease the total escape rates of ions.

Using over 6 years of magnetic field data (2014.10-2020.12) collected by the Mars Atmosphere and Volatile EvolutioN (MAVEN), we conduct a statistical study on the three-dimensional average magnetic field structure around Mars. We find that this magnetic field structure conforms to the pattern typical of an induced magnetosphere, that is, the interplanetary magnetic field (IMF) which is carried by the solar wind and which drapes, piles up, slips around the planet, and eventually forms a tail in the wake. The draped field lines from both hemispheres along the direction of the solar wind electric field (E) are directed towards the nightside magnetic equatorial plane, which looks like they are “sinking” toward the wake. These “sinking” field lines from the +E-hemisphere (E pointing away from the plane) are more flared and dominant in the tail, while the field lines from the –E-hemisphere (E pointing towards) are more stretched and “pinched” towards the plasma sheet. Such highly “pinched” field lines even form a loop over the pole of the –E-hemisphere. The tail current sheet also shows an E-asymmetry: the sheet is thicker with a stronger tailward J×B force at +E-flank, but much thinner and with a weaker J×B (even turns sunward) at –E-flank. Additionally, we find that IMF Bx can induce a kink-like field structure at the boundary layer; the field strength is globally enhanced and the field lines flare less during high dynamic pressure; however, the rotation of the planet, against expectations, modulate the configuration of the tail current sheet insignificantly.

We present partial ring distributions of solar wind protons observed by the Rosetta spacecraft at comet 67P/Churyumov-Gerasimenko. The formation of ring distributions is usually associated with high activity comets, where the spatial scales are larger than multiple ion gyroradii. Our observations are made at a low-activity comet at a heliocentric distance of 2.8 AU on April 19th, 2016, and the partial rings occur at a spatial scale comparable to the ion gyroradius. We use a new visualisation method to simultaneously show the angular distribution of median energy and differential flux. A fitting procedure extracts the bulk speed of the solar wind protons, separated into components parallel and perpendicular to the gyration plane, as well as the gyration velocity. The results are compared with models and put into context of the global comet environment. We find that the formation mechanism of these partial rings of solar wind protons is entirely different from the well-known partial rings of cometary pickup ions at high-activity comets. A density enhancement layer of solar wind protons around the comet is a focal point for proton trajectories originating from different regions of the upstream solar wind. If the spacecraft location coincides with this density enhancement layer, the different trajectories are observed as an energy-angle dispersion and manifest as partial rings in velocity space.