Pickup Ions Research Articles

Statistical distribution of mirror-mode-like structures in the magnetosheaths of unmagnetised planets – Part 1: Mars as observed by the MAVEN spacecraft

Abstract. In this series of papers, we present statistical maps of mirror-mode-like (MM) structures in the magnetosheaths of Mars and Venus and calculate the probability of detecting them in spacecraft data. We aim to study and compare them with the same tools and a similar payload at both planets. We consider their dependence on extreme ultraviolet (EUV) solar flux levels (high and low) and, specific to Mars, on Mars Year (MY) as well as atmospheric seasons (four solar longitudes Ls). We first use magnetic-field-only criteria to detect these structures and present ways to mitigate ambiguities in their nature. In line with many previous studies at Earth, this technique has the advantage of using one instrument (a magnetometer) with good time resolution, facilitating comparisons between planetary and cometary environments. Applied to the magnetometer data of the Mars Atmosphere and Volatile EvolutioN (MAVEN) spacecraft from November 2014 to February 2021 (MY32–MY35), we detect events closely resembling MMs lasting in total more than 170 000 s, corresponding to about 0.1 % of MAVEN's total time spent in the Martian plasma environment. We calculate MM-like occurrences normalised to the spacecraft's residence time during the course of the mission. Detection probabilities are about 1 % at most for any given controlling parameter. In general, MM-like structures appear in two main regions: one behind the shock and the other close to the induced magnetospheric boundary, as expected from theory. Detection probabilities are higher on average in low-solar-EUV conditions, whereas high-solar-EUV conditions see an increase in detections within the magnetospheric tail. We tentatively link the former tendency to two combining effects: the favouring of ion cyclotron waves the closer to perihelion due to plasma beta effects and, possibly, the non-gyrotropy of pickup ion distributions. This study is the first of two on the magnetosheaths of Mars and Venus.

Pickup Ion–Mediated Magnetic Reconnection in the Outer Heliosphere

Pickup ions (PUIs) play a crucial role in the heliosphere, contributing to the mediation of large-scale structures such as the distant solar wind, the heliospheric termination shock, and the heliopause. While magnetic reconnection is thought to be a common process in the heliosphere due to the presence of heliospheric current sheets, it is poorly understood how PUIs might affect the evolution of magnetic reconnection. Although it is reasonable to suppose that PUIs decrease the reconnection rate since the plasma beta becomes much larger than 1 when PUIs are included, we show for the first time that such a supposition is invalid and that PUI-induced turbulence, heat conduction, and viscosity can preferentially boost magnetic reconnection in heliospheric current sheets in the distant solar wind. This suggests that it is critical to include the effect of the turbulence, heat conduction, and viscosity caused by PUIs to understand the dynamics of magnetic reconnection in the outer heliosphere.

How catching the interstellar wind in the inner solar system led the way on a road to interdisciplinary research between heliophysics and astrophysics

Combined in situ observations of the interstellar wind through the solar system and of its pickup ions (PUIs), implanted after ionization in the solar wind, explain, in comparison with interstellar absorption lines of nearby stars, that the Sun is in an interaction region of the two nearest interstellar clouds. This new finding disrupts the long-held understanding that we are inside the local interstellar cloud (LIC). We discuss how space physics evolved toward such interdisciplinary studies between heliophysics and astrophysics. In 1984, the discovery of interstellar He+ PUIs exposed the very local interstellar medium to in situ diagnostics at 1AU. These PUIs provide the interstellar gas composition and form a stepping stone for the acceleration of ions, especially into anomalous cosmic rays. Using the Sun as a gravitational spectrograph, direct imaging of the neutral interstellar wind, first for He and then for H, O, and Ne, provides the interstellar gas velocity vector and temperature at the heliopause. Combining the interstellar gas flow vectors, those of secondary neutral He and O, and the interstellar magnetic field direction deduced from the interstellar H deflection and termination shock anisotropy seen by the Voyagers provides synergistically the heliosphere’s shape, its interaction with the interstellar medium, and constrains our radiation environment. This ISMF organizes the bright Ribbon seen in all-sky images of energetic neutral atoms with the potential to provide its precision determination. The elemental and isotopic composition from PUI and neutral gas observations constrains the galactic evolution and Big Bang cosmology, opening additional interdisciplinary opportunities.

The Role of Electrons and Helium Atoms in Global Modeling of the Heliosphere

We present a new three-dimensional, MHD-plasma/kinetic-neutrals model of the solar wind (SW) interaction with the local interstellar medium (LISM), which self-consistently includes neutral hydrogen and helium atoms. This new model also treats electrons as a separate fluid and includes the effect of Coulomb collisions. While the properties of electrons in the distant SW and in the LISM are mostly unknown due to the lack of in situ observations, a common assumption for any global, single-ion model is to assume that electrons have the temperature of the ion mixture, which includes pickup ions. In the new model, electrons in the SW are colder, which results in a better agreement with New Horizons observations in the supersonic SW. In the LISM, however, ions and electrons are almost in thermal equilibrium. As for the plasma mixture, the major differences between the models are in the inner heliosheath, where the new model predicts a charge-exchange-driven cooling and a decrease of the heliosheath thickness. The filtration of interstellar neutral atoms at the heliospheric interface is discussed. The new model predicts an increase in the H density by ∼2% at 1 au. However, the fraction of pristine H atoms decreases by ∼12%, while the density of atoms born in the outer and inner heliosheath increases by 5% and ∼35%, respectively. While at 1 au the density of He atoms remains unchanged, the contribution from the “warm breeze” increases by ∼3%.

Probing the Length of the Heliospheric Tail with Energetic Neutral Atoms (ENAs) from 0.52 to 80 keV

The shape of the heliosphere is currently under active debate. Energetic neutral atoms (ENAs) offer the best method for investigating the global structure of the heliosphere. To date, the Interstellar Boundary Explorer (IBEX) and the Ion and Neutral Camera (INCA) that was on board Cassini provide the only global ENA observations of the heliosphere. While extensive modeling has been done at IBEX-Hi energies (0.52–6 keV), no global ENA modeling has been conducted for INCA energies (5.2–55 keV). Here, we use an ENA model of the heliosphere based on hybrid results that capture the heating and acceleration of pickup ions (PUIs) at the termination shock to compare modeled global ENA results with IBEX-Hi and INCA observations using both a long- and short-tail model of the heliosphere. We find that the modeled ENA results for the two heliotail configurations produce similar results from the IBEX-Hi through the INCA energies. We conclude from our modeled ENAs, which only include PUI acceleration at the termination shock, that ENA observations in currently available energy ranges are insufficient for probing the shape and length of the heliotail. However, as a prediction for the future IMAP-Ultra mission (3–300 keV) we present modeled ENA maps at 80 keV, where the cooling length (∼600 au) is greater than the distance where the long- and short-heliotail models differ (∼400 au), and find that IMAP-Ultra should be able to identify the shape of the heliotail, predicting differences in the north lobe to downwind flux ratio between the models at 48%.

Low-frequency Waves Due to Newborn Interstellar Pickup Ions Observed from 43 to 47 au by the Voyager 1 Spacecraft

Interstellar neutral atoms enter the heliosphere at a relatively slow speed corresponding to the motion of the Sun through the local interstellar medium, which is approximately 25 km s−1. Neutral hydrogen atoms enter from the approximate location of the Voyager spacecraft and are eventually ionized primarily by collision with thermal solar wind ions. An earlier analysis by Hollick et al. examined low-frequency magnetic waves observed by the Voyager spacecraft from launch through 1990 that are thought to arise from the scattering of newborn interstellar pickup H+ and He+. We report an analysis of Voyager 1 observations in 1991, which is the last year of high-resolution magnetic field data that are publicly available, and find 70 examples of low-frequency waves with the characteristics that suggest excitation by pickup H+ and 10 examples of waves consistent with excitation by pickup He+. We find a particularly dense cluster of observations at the tail end of what is thought to be a Merged Interaction Region (MIR) that was previously studied by Burlaga & Ness using Voyager 2 observations. This is not unexpected if the MIR is followed by a large rarefaction region, as they tend to be regions of reduced turbulence levels that permit the growth of the waves over the long time periods that are generally required of this instability.

Solar Control of the Pickup Ion Plume in the Dayside Magnetosheath of Venus

AbstractUsing the 8.5‐year Venus Express measurements, we demonstrate the asymmetric plasma distributions in the Venusian magnetosheath. An escaping plume is formed by pickup oxygen ions in the hemisphere where the motional electric field points outward from Venus, while the velocity of solar wind protons is faster in the opposite hemisphere. The pickup O+ escape rate is estimated to be (3.6 ± 1.4) × 1024 s−1 at solar maximum, which is comparable to the ion loss rate through the magnetotail, and (1.3 ± 0.4) × 1024 s−1 at solar minimum. The increase of O+ fluxes with extreme ultraviolet (EUV) intensity is significant upstream of the bow shock, partially attributed to the increase of exospheric neutral oxygen density. However, the solar wind velocity just has a slight effect on the pickup O+ escape rate in the magnetosheath, while the effect of solar wind density is not observed. Our results suggest the pickup O+ escape rate is mainly controlled by EUV radiation.

Hybrid Simulation of Proton Cyclotron Waves Upstream of Mars Generated by Pickup Ion Beam Distribution Data

Hybrid Simulation of Proton Cyclotron Waves Upstream of Mars Generated by Pickup Ion Beam Distribution Data

Relating Energetic Ion Spectra to Energetic Neutral Atoms

Heliospheric energetic neutral atoms (ENAs) originate from energetic ions that are neutralized by charge exchange with neutral atoms in the heliosheath and very local interstellar medium (VLISM). Since neutral atoms are unaffected by electromagnetic fields, they propagate ballistically with the same speeds as parent particles. Consequently, measurements of ENA distributions allow one to remotely image the energetic ion distributions in the heliosheath and VLISM. The origin of the energetic ions that spawn ENAs is still debated, particularly at energies higher than ∼keV. In this work, we summarize five possible sources of energetic ions in the heliosheath that cover the ENA energy from a few keV to hundreds of keV. Three sources of the energetic ions are related to pickup ions (PUIs): those PUIs transmitted across the heliospheric termination shock (HTS), those reflected once or multiple times at the HTS, i.e., reflected PUIs, and those PUIs multiply reflected and further accelerated by the HTS. Two other kinds of ions that can be considered are ions transmitted from the suprathermal tail of the PUI distribution and other particles accelerated at the HTS. By way of illustration, we use these energetic particle distributions, taking account of their evolution in the heliosheath, to calculate the ENA intensities and to analyze the characteristics of ENA spectra observed at 1 au.

An Anomalous Cosmic-Ray Mediated Termination Shock: Implications for Energetic Neutral Atoms

The Voyager 2 crossing of the termination shock indicated that most of the upstream energy from the thermal solar wind ions was transferred to pickup ions (PUIs) and other energetic particles downstream of the shock. We use hybrid simulations at the termination shock for the Voyager 2, flank, and tail directions to evaluate the distributions of different ion species downstream of the shock over the energy range of 0.52–55 keV. Here, we extend the work of Gkioulidou et al., which showed an energy-dependent discrepancy between modeled and energetic neutral atom (ENA) observations, and fit distributions to a hybrid model to show that a population of PUIs accelerated via diffusive shock acceleration (DSA) to become low-energy anomalous cosmic rays (ACRs) can bridge the gap between modeled and observed ENA fluxes. Our results with the inclusion of DSA via hybrid fitting give entirely new and novel evidence that DSA at the termination shock is likely to be an important physical process. These ACRs carry a significant fraction of the energy density at the termination shock (22%, 13%, and 19% in the Voyager 2, flank, and tail directions, respectively). Using these ACRs in global ENA modeling of the heliosphere from 0.52 to 55 keV, we find that scaling factors as large as 1.8–2.5 are no longer required to match ENA observations at energies of ∼1–4 keV. Large discrepancies between modeled and observed ENAs only remain over energies of 4–20 keV, indicating that there may be a further acceleration mechanism in the heliosheath at these energies.

Turbulence Driving by Interstellar Pickup Ions in the Outer Solar Wind

We revisit the question of how the unstable scattering of interstellar pickup ions (PUIs) may drive turbulence in the outer solar wind and why the energy released into fluctuations by this scattering appears to be significantly less than the standard bispherical prediction. We suggest that energization of the newly picked-up ions by the ambient turbulence during the scattering process can result in a more spherical distribution of PUIs and reduce the generated fluctuation energy to a level consistent with the observations of turbulent intensities and core solar wind heating. This scenario implies the operation of a self-regulation mechanism that maintains the observed conditions of turbulence and heating in the PUI-dominated solar wind.

Solar Wind Protons Forming Partial Ring Distributions at Comet 67P

AbstractWe 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 19 April 2016, and the partial rings occur at a spatial scale comparable to the ion gyroradius. We use a new visualization 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.

Statistical Properties of Plateau‐Like Turbulence Spectra in the Martian Magnetosheath: MAVEN Observations

AbstractThe Martian magnetosheath provides us with a natural laboratory to study plasma turbulence in the presence of pickup ions and locally generated instabilities. Unlike the typical magnetic‐field spectra with a single spectral scaling at magnetohydrodynamics (MHD) scales in Earth's magnetosheath, the magnetic‐field spectra in the Martian magnetosheath during 4 years of Mars Atmosphere and Volatile EvolutioN observations frequently present an additional spectral break‐point with a shallow slope at MHD scales which we define as a plateau‐like spectral feature. The average occurrence rate of plateau‐like magnetic‐field spectra is 56.6% of our measurement intervals. At moderate pick‐up angles, the occurrence rate increases to a maximum of ∼70.0%. Furthermore, we present a positive correlation with the local ion density and anti‐correlations with the local βi and the solar Extreme Ultra Violet irradiance. A similar occurrence rate in the quasi‐perpendicular and the quasi‐parallel magnetosheath (60.1% vs. 52.9%) indicates that the plateau‐like spectra are more likely formed locally than in the upstream solar wind. Our results suggest that energy injection from pickup ion driven micro‐instabilities, for example, in the form of proton cyclotron waves, has insufficient time to evolve into a fully developed cascade in such a confined space like the Martian magnetosheath.

A MAVEN Case Study of Radial IMF at Mars: Impacts on the Dayside Ionosphere

AbstractThe solar wind interaction with Mars controls the transfer of energy and momentum from the solar wind into the magnetosphere, ionosphere and atmosphere, driving structure, and dynamics within each. This interaction is highly dependent on the upstream Interplanetary Magnetic Field (IMF) orientation. We use in‐situ plasma measurements made by the Mars Atmosphere and Volatile EvolutioN (MAVEN) mission to identify several prominent features that arise when the IMF is aligned approximately parallel or antiparallel to solar wind flow (conditions known as “radial IMF”). In particular, solar wind protons and alphas are observed to directly penetrate down to periapsis altitudes, while the magnetic barrier forms deep within the dayside ionosphere. The MAVEN observations are consistent with either an ionopause‐like boundary or diamagnetic cavity forming beneath the barrier, as a consequence of the dense cold ionosphere and the absence of significant crustal magnetic fields at this periapsis location. The planetary ions above the magnetic barrier are exposed to solar wind flow and subsequent mass‐loading. The (convective electric field or “ion pickup”) force is weak and highly variable during radial IMF. While wave particle interactions and subsequent wave heating contribute to incorporating the heavy planetary ions into the solar wind flow, the solar wind momentum is not fully deflected around the obstacle and is delivered into the collisional atmosphere. Significant ion heating is observed deep within the dayside ionosphere, and observed ionospheric density and temperature profiles demonstrate that these ion energization mechanisms drive significant erosion and likely escape to space.

MAVEN Proton Observations Near the Martian Moon Phobos: Does Phobos Backscatter Solar Wind Protons?

AbstractESA's Mars Express (MEX) may have observed twice solar wind protons backscattered by Phobos. However, these detections remain uncertain and call for an independent confirmation. Here, we analyze the proton measurements collected by the Suprathermal And Thermal Ion Composition experiment onboard NASA's Mars Atmosphere and Volatile EvolutioN mission (MAVEN) during the 15 closest encounters of MAVEN with Phobos (closest approach between 80 and 370 km). We use a model of magnetic connectivity to the Martian bow shock and an exospheric pickup ion model to identify the origin of the protons observed close to the Martian moon. We find no evidence of protons backscattered by Phobos. Instead, foreshock‐shock foot protons and Martian exospheric pickup protons are often observed near Phobos. These results put into new context not only the past MEX detections, but also the future observation attempts of the upcoming JAXA Martian Moons Exploration mission.

Observations of a Solar Energetic Particle Event From Inside and Outside the Coma of Comet 67P.

We analyze observations of a solar energetic particle (SEP) event at Rosetta's target comet 67P/Churyumov-Gerasimenko during 6-10 March 2015. The comet was 2.15AU from the Sun, with the Rosetta spacecraft approximately 70km from the nucleus placing it deep inside the comet's coma and allowing us to study its response. The Eastern flank of an interplanetary coronal mass ejection (ICME) also encountered Rosetta on 6 and 7 March. Rosetta Plasma Consortium data indicate increases in ionization rates, and cometary water group pickup ions exceeding 1keV. Increased charge exchange reactions between solar wind ions and cometary neutrals also indicate increased upstream neutral populations consistent with enhanced SEP induced surface activity. In addition, the most intense parts of the event coincide with observations interpreted as an infant cometary bow shock, indicating that the SEPs may have enhanced the formation and/or intensified the observations. These solar transient events may also have pushed the cometopause closer to the nucleus. We track and discuss characteristics of the SEP event using remote observations by SOHO, WIND, and GOES at the Sun, in situ measurements at Solar Terrestrial Relations Observatory Ahead, Mars and Rosetta, and ENLIL modeling. Based on its relatively prolonged duration, gradual and anisotropic nature, and broad angular spread in the heliosphere, we determine the main particle acceleration source to be a distant ICME which emerged from the Sun on 6 March 2015 and was detected locally in the Martian ionosphere but was never encountered by 67P directly. The ICME's shock produced SEPs for several days which traveled to the in situ observation sites via magnetic field line connections.

Exploring turbulence from the Sun to the local interstellar medium: Current challenges and perspectives for future space missions

Turbulence is ubiquitous in space plasmas. It is one of the most important subjects in heliospheric physics, as it plays a fundamental role in the solar wind—local interstellar medium interaction and in controlling energetic particle transport and acceleration processes. Understanding the properties of turbulence in various regions of the heliosphere with vastly different conditions can lead to answers to many unsolved questions opened up by observations of the magnetic field, plasma, pickup ions, energetic particles, radio and UV emissions, and so on. Several space missions have helped us gain preliminary knowledge on turbulence in the outer heliosphere and the very local interstellar medium. Among the past few missions, theVoyagershave paved the way for such investigations. This paper summarizes the open challenges and voices our support for the development of future missions dedicated to the study of turbulence throughout the heliosphere and beyond.

Oblique Electromagnetic Instabilities Driven by Pickup Ion Ring-beam Distributions in the Outer Heliosheath. II. Hybrid Simulations

Hybrid simulations are carried out to study the instabilities driven by the outer heliosheath pickup ions of ring-beam velocity distributions. Following the linear instability analysis in the companion paper, the electromagnetic instabilities at various propagation directions (0° ≤ θ ≤ 180°) are studied for pickup ion ring-beam distributions of different pickup angles (0° ≤ α ≤ 90°). Previous studies of the ring-beam pickup ion distributions have been limited to the parallel- and antiparallel-propagating modes (θ = 0° and 180°). In this paper, several one-dimensional hybrid simulations are performed to investigate the most unstable modes of different propagation angles predicted by the linear instability analysis. Moreover, two-dimensional hybrid simulations are conducted to examine the simultaneous excitation of the parallel, antiparallel, and oblique unstable modes and their role in the scattering of the outer heliosheath pickup ions. Our recent one-dimensional hybrid and particle-in-cell simulations with simulation domains along the background magnetic field demonstrated that at small pickup angles, the outer heliosheath pickup ions experience an incomplete scattering limited to the hemisphere of positive parallel velocities with respect to the background magnetic field. In contrast, the present two-dimensional hybrid simulation results show that when both parallel/antiparallel and oblique unstable modes are acting together, the pickup ions can be scattered into the hemisphere of negative parallel velocities and approach a near-isotropic velocity distribution. This is in conflict with the spatial retention scenario of the energetic neutral atom ribbon observed by the Interstellar Boundary Explorer, which requires pickup ions at small pickup angles to stream away from the ribbon directions.

Oblique Electromagnetic Instabilities Driven by Pickup Ion Ring-beam Distributions in the Outer Heliosheath. I. Linear Instability Analysis

The energetic neutral atom ribbon observed by the Interstellar Boundary Explorer spacecraft is believed to originate from the pickup ions in the outer heliosheath. The outer heliosheath pickup ions generally have a ring-beam velocity distribution at a certain pickup angle, α, the angle at which these ions are picked up by the interstellar magnetic field. The pickup ion ring-beam distributions can drive unstable waves of different propagation angles with respect to the background interstellar magnetic field, θ. Previous studies of the outer heliosheath pickup ion dynamics were mainly focused on ring-like pickup ion distributions with α ≈ 90° and/or the parallel- and antiparallel-propagating unstable waves (θ = 0° and 180°). The present study carries out linear kinetic instability analysis to investigate both the parallel and oblique unstable modes (0° ≤ θ ≤ 180°) driven by ring-beam pickup ion distributions of different pickup angles between 0° and 90°. Our linear instability analysis reveals that ring-beam pickup ion distributions can excite oblique mirror waves as well as parallel/quasi-parallel and oblique right- and left-helicity waves. The maximum growth rate among all the instabilities belongs to the parallel-propagating left-helicity waves at most pickup angles. Furthermore, the evolution of the unstable mirror waves by varying pickup angle indicates that as the pickup angle increases, the maximum growth rate of the mirror modes increases, while its propagation angle decreases.

Shock-wave heating mechanism of the distant solar wind: Explanation of Voyager-2 data

Aims. One of the important discoveries made by Voyager-2 is the nonadiabatic radial profile of the solar wind proton temperature. This phenomenon has been studied for several decades. The dissipation of turbulence energy has been proposed as the main physical process responsible for the temperature profile. The turbulence is both convected with the solar wind and originated in the solar wind by the compressions and shears in the flows and by pick-up ions. The compression source of the solar wind heating in the outer heliosphere appears due to shock waves, which originated either in the solar corona or in the solar wind itself. The goal of this work is to demonstrate that the shock-wave heating itself is enough to explain the temperature profile obtained by Voyager-2. Methods. The effect of shock-wave heating is demonstrated in the frame of a very simple spherically symmetric high-resolution (in both space and time) gas-dynamical data-driven solar wind model. This data-driven model employs the solar-wind parameters at 1 AU with minute resolution. The data are taken from the NASA OMNIWeb database. It is important to underline that (1) the model captures the shocks traveling and/or originating in the solar wind, and (2) other sources of heating are not taken into account in the model. We extended this simple model to the magnetohydrodynamic (MHD) and two-component models and found very similar results. Results. The results of the numerical modeling with the one-minute OMNI data as the boundary condition show very good agreement with the solar-wind temperature profiles obtained by Voyager-2. It is also noteworthy that the numerical results with daily averaged OMNI data show a very similar temperature profile, while the numerical runs with 27-day-averaged OMNI data demonstrate the adiabatic behavior of the temperature.