Waves at the electron plasma frequency are found throughout the heliosphere. They provide indicators of unstable electron distributions, are routinely used to estimate the local electron number density, and can lead to radio wave emission at the plasma frequency and its harmonics. Although they have been studied extensively in various solar and heliospheric plasma regions, there is a lack of statistical studies of plasma frequency waves in Earth’s magnetotail. Here, the occurrence and properties of plasma frequency waves, namely Langmuir and upper hybrid waves, are investigated in Earth’s magnetotail using the four Magnetospheric Multiscale spacecraft. In Earth’s magnetotail plasma frequency waves are observed about $1$~\% of the time. About $80$~\% of the waves are identified as Langmuir waves, while about $20$~\% are identified as upper hybrid waves. The waves are primarily found in the plasma sheet boundary layer. By comparing with the local electron distributions it is shown that the Langmuir waves are generated by the bump-on-tail instability, while upper hybrid waves are typically associated with broad electron beams or loss-cone-like distributions. The majority of the waves are found in close proximity to ion outflow regions associated with magnetic reconnection in the magnetotail. The waves are likely generated by plasma sheet electrons escaping along newly reconnected magnetic field lines or electron beams propagating toward the distant magnetotail.

Low Mach number collisionless shocks are routinely observed in the solar wind and upstream of planetary bodies. However, most in situ observations have lacked the necessary temporal resolution to directly study the kinetic behavior of ions across these shocks. We investigate a series of five low Mach number bow shock crossings observed by the Magnetospheric Multiscale (MMS) mission. The five shocks had comparable Mach numbers, but varying shock-normal angles ($66^{\circ} \lesssim \theta_{Bn} \lesssim 89^{\circ}$) and ramp widths ($5 \mathrm{km} \lesssim l \lesssim 100 \mathrm{km}$). The shock width is shown to be crucial in determining the fraction of protons reflected and energized by the shock, with proton reflection increasing with decreasing shock width. As the shock width increases proton reflection is arrested entirely. For nearly perpendicular shocks, reflected protons exhibit quasi-periodic structures, which persist far downstream of the shock. As the shock-normal angle becomes more oblique these periodic proton structures broaden to form an energetic halo population. Periodic fluctuations in the magnetic field downstream of the shocks are generated by fluctuations in dynamic pressure of alpha particles, which are decelerated by the cross-shock potential and subsequently undergo gyrophase bunching. These results demonstrate that complex kinetic-scale ion dynamics occur in low Mach number shocks, which depend significantly on the shock profile.

When the velocity shear between the two plasmas separated by Earth’s magnetopause is locally super-Alfvénic, the Kelvin-Helmholtz (KH) instability can develop. A crucial role is played by the interplanetary magnetic field (IMF) orientation, which can stabilize the velocity shear. Although, in a linear regime, the instability threshold is equally satisfied during both northward and southward IMF orientations, in-situ measurements show that KH instability is preferentially excited during the northward IMF orientation. We investigate this different behavior by means of a mixing parameter which we apply to two KH events to identify both boundaries and the center of waves/vortices. During the northward orientation, the waves/vortex boundaries have stronger electrons than ions mixing, while the opposite is observed at their center. During the southward orientation, instead, particle mixing is observed predominantly at the boundaries. In addition, stronger local ion and electron non-thermal features are observed during the northward than the southward IMF orientation. Specifically, ion distribution functions are more distorted, due to field-aligned beams, and electrons have a larger temperature anisotropy during the northward than the southward IMF orientation. The observed kinetic features provide an insight into both local and remote processes that affect the evolution of KH structures.

We investigate electrostatic waves in a magnetopause reconnection diffusion region event. In the electron diffusion region on the magnetospheric side, an oblique electrostatic wave is observed. The local distribution exhibits fast non-gyrotropic electron beams with drifts comparable to the electron thermal speed, but the wave has a much lower phase speed. Response of ions and possible cold electrons may contribute to wave excitation. Near the current sheet mid-plane, parallel electron beam-mode waves are modulated by whistler waves. In the separatrix region, parallel waves associated with field-aligned electron beams and perpendicular electron cyclotron waves with loss cone distributions exhibit modulation frequencies in the lower-hybrid wave frequency range. We infer that lower-hybrid waves scatter electrons to produce beams and alter loss cones to modulate electrostatic waves. The results advance our understanding about the regimes and mechanisms of electrostatic waves in reconnection and their coupling with lower-frequency waves.

Electron temperature anisotropy-driven instabilities such as the electron firehose instability (EFI) are especially significant in space collisionless plasmas, where collisions are so scarce that wave-particle interactions are the leading mechanisms in the isotropization of the distribution function and energy transfer. Observational statistical studies provided convincing evidence in favor of the EFI constraining the electron distribution function and limiting the electron temperature anisotropy. Magnetic reconnection is characterized by regions of enhanced temperature anisotropy that could drive instabilities – including the electron firehose instability – affecting the particle dynamics and the energy conversion. However, in situ observations of the fluctuations generated by the EFI are still lacking and the interplay between magnetic reconnection and EFI is still largely unknown. In this study, we use high-resolution in situ measurements by the Magnetospheric Multiscale (MMS) spacecraft to identify and investigate EFI fluctuations in the magnetic reconnection exhaust in the Earth’s magnetotail. We find that the wave properties of the observed fluctuations largely agree with theoretical predictions of the non-propagating EF mode. These findings are further supported by comparison with the linear kinetic dispersion relation. Our results demonstrate that the magnetic reconnection outflow can be the seedbed of EFI and provide the first direct in situ observations of EFI-generated fluctuations.

We report Magnetospheric Multiscale observations of large amplitude, parallel, electrostatic, proton plasma frequency waves on the magnetospheric side of the reconnecting magnetopause. The waves are often found in the magnetospheric separatrix region and in the outflow near the magnetospheric ion edge. Statistical results from five months of data show that these waves are closely tied to the presence of cold (typically tens of eV) ions, found for 88% of waves near the separatrix region, and that plasma properties are consistent with ion acoustic wavegrowth. We analyze one wave event in detail, concluding that the wave is ion acoustic. We provide a simple explanation for the mechanisms leading to the development of the ion acoustic instability. These waves can be important for separatrix dynamics by heating the cold ion component and providing a mechanism to damp the kinetic Alfvén waves propagating away from the reconnection site.

Lower hybrid waves are investigated at the magnetosheath separatrix region in asymmetric guide-field reconnection by using the Magnetospheric Multiscale (MMS) mission. Three of the four MMS spacecraft observe clear wave activities around the lower hybrid frequency across the magnetosheath separatrix, where a density gradient is present. The observed waves are consistent with generation by the lower hybrid drift instability. The characteristic properties of these waves include: (1) the waves propagate toward the x-line in the spacecraft frame due to the large out-of-plane magnetic field, which is in the same direction of the diamagnetic drift of the x-line; (2) the wave potential is about 20\% of the electron temperature. These drift waves effectively produce cross-field particle diffusion, enabling the transport of magnetosheath electrons into the exhaust region. The observations presented in this study indicate unique features of asymmetric guide-field reconnection.

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.

Whistler waves are often observed in magnetopause reconnection associated with electron beams. We analyze seven MMS crossings surrounding the electron diffusion region (EDR) to study the role of electron beams in whistler excitation. Waves have two major types: (1) Narrow-band waves with high ellipticities and (2) broad-band waves that are more electrostatic with significant variations in ellipticities and wave normal angles. While both types of waves are associated with electron beams, the key difference is the anisotropy of the background population, with perpendicular and parallel anisotropies, respectively. The linear instability analysis suggests that the first type of wave is mainly due to the background anisotropy, with the beam contributing additional cyclotron resonance to enhance the wave growth. The second type of distribution excites broadband waves via Landau resonance, and as seen in one event, the beam anisotropy induces an additional cyclotron mode. The results are supported by particle-in-cell simulations. We infer that the first type occurs downstream of the central EDR, where background electrons experience Betatron acceleration to form the perpendicular anisotropy; the second type occurs in the central EDR of guide field reconnection. A parametric study is conducted with linear instability analysis. A beam anisotropy alone of above ~3 likely excites the cyclotron mode waves. Large beam drifts cause Doppler shifts and may lead to left-hand polarizations in the ion frame. Future studies are needed to determine whether the observation covers a broader parameter regime and to understand the competition between whistler and other instabilities.

We report observations of the ion dynamics inside an Alfven branch wave that propagates near the reconnecting dayside magnetopause. The measured frequency, wave normal angle and polarization are within 1% with the predictions of a dispersion solver, and indicate that the wave is an electromagnetic ion cyclotron wave with very oblique wave vector. The magnetospheric plasma contains hot protons (keV), cold protons (eV), plus some heavy ions. The cold protons follow the magnetic field fluctuations and remain frozen-in, while the hot protons are at the limit of magnetization. The cold proton velocity fluctuations contribute to balance the Hall term in Ohm’s law, allowing the wave polarization to be highly-elliptical and right-handed, a necessary condition for propagation at oblique wave normal angles. The dispersion solver indicates that increasing the cold proton density facilitates generation and propagation of these waves at oblique angles, as it occurs for the observed wave.

Plasmas in Earth’s outer magnetosphere, magnetosheath, and solar wind are essentially collisionless. This means particle distributions are not typically in thermodynamic equilibrium and deviate significantly from Maxwellian distributions. The deviations of these distributions can be further enhanced by plasma processes, such as shocks, turbulence, and magnetic reconnection. Such distributions can be unstable to a wide variety of kinetic plasma instabilities, which in turn modify the electron distributions. In this paper the deviations of the observed electron distributions from a bi-Maxwellian distribution function is calculated and quantified using data from the Magnetospheric Multiscale (MMS) spacecraft. A statistical study from tens of millions of electron distributions shows that the primary source of the observed non-Maxwellianity are electron distributions consisting of distinct hot and cold components in Earth’s low-density magnetosphere. This results in large non-Maxwellianities in at low densities. However, after performing a stastical study we find regions where large non-Maxwellianities are observed for a given density. Highly non-Maxwellian distributions are routinely found are Earth’s bowshock, in Earth’s outer magnetosphere, and in the electron diffusion regions of magnetic reconnection. Enhanced non-Maxwellianities are observed in the turbulent magnetosheath, but are intermittent and are not correlated with local processes. The causes of enhanced non-Maxwellianities are investigated.

Mirror modes are ubiquitous in space plasma and grow from pressure anisotropy. Together with other instabilities, they play a fundamental role in constraining the free energy contained in the plasma. This study focuses on mirror modes observed in the solar wind by Solar Orbiter for heliocentric distances between 0.5 and 1 AU. Typically, mirror modes have timescales from several to tens of seconds and are considered quasi-MHD structures. In the solar wind, they also generally appear as isolated structures. However, in certain conditions, prolonged and bursty trains of higher frequency mirror modes are measured, which have been labeled previously as mirror mode storms. At present, only a handful of existing studies have focused on mirror mode storms, meaning that many open questions remain. In this study, Solar Orbiter has been used to investigate several key aspects of mirror mode storms: their dependence on heliocentric distance, association with local plasma properties, temporal/spatial scale, amplitude, and connections with larger-scale solar wind transients. The main results are that mirror mode storms often approach local ion scales and can no longer be treated as quasi-MHD, thus breaking the commonly used long-wavelength assumption. They are typically observed close to current sheets and downstream of interplanetary shocks. The events were observed during slow solar wind speeds and there was a tendency for higher occurrence closer to the Sun. The occurrence is low, so they do not play a fundamental role in regulating ambient solar wind but may play a larger role inside transients.