Ion Inertial Scale Research Articles

Nonlinear dynamics of large-amplitude, small-scale Alfvén waves

We study large-amplitude, very oblique Alfvén waves at low β, with small gradient length scales, comparable to the ion inertial scale di. Such waves have large density fluctuations and slight dispersion from finite-frequency and finite-ion sound radius effects. We derive a weakly nonlinear evolution equation governing the behavior of the waves in one dimension and categorize the different solitons appearing in different regimes: the regular solitons involve full rotations of the transverse magnetic field similar to modified Korteweg–de Vries (mKdV) solitons (our nonlinear equation reduces to the mKdV equation in the long-wavelength limit). However, for sufficiently small soliton widths, some become singular, small-amplitude solitons with density discontinuities and, are, thus expected to become strongly dissipative in a real plasma. These solutions may be useful in explaining some aspects of the sharp, ion-scale magnetic field rotations (switchbacks) observed in the near-Sun solar wind by Parker Solar Probe.

Multi‐Scale Observation of Magnetotail Reconnection Onset: 2. Microscopic Dynamics

AbstractWe analyze the local dynamics of magnetotail reconnection onset using Magnetospheric Multiscale (MMS) data. In conjunction with MMS, the macroscopic dynamics of this event were captured by a number of other ground and space‐based observatories, as is reported in a companion paper. We find that the local dynamics of the onset were characterized by the rapid thinning of the cross‐tail current sheet below the ion inertial scale, accompanied by the growth of flapping waves and the subsequent onset of electron tearing. Multiple kinetic‐scale magnetic islands were detected coincident with the growth of an initially sub‐Alfvénic, demagnetized tailward ion exhaust. The onset and rapid enhancement of parallel electron inflow at the exhaust boundary was a remote signature of the intensification of reconnection Earthward of the spacecraft. Two secondary reconnection sites are found embedded within the exhaust from a primary X‐line. The primary X‐line was designated as such on the basis that (a) while multiple jet reversals were observed in the current sheet, only one reversal of the electron inflow was observed at the high‐latitude exhaust boundary, (b) the reconnection electric field was roughly five times larger at the primary X‐line than the secondary X‐lines, and (c) energetic electron fluxes increased and transitioned from anti‐field‐aligned to isotropic during the primary X‐line crossing, indicating a change in magnetic topology. The results are consistent with the idea that a primary X‐line mediates the reconnection of lobe magnetic field lines and accelerates electrons more efficiently than its secondary X‐line counterparts.

Efficient kinetic Lattice Boltzmann simulation of three-dimensional Hall-MHD turbulence

Simulating plasmas in the Hall-magnetohydrodynamics (Hall-MHD) regime represents a valuable approach for the investigation of complex nonlinear dynamics developing in astrophysical frameworks and fusion machines. The Hall electric field is computationally very challenging as it involves the integration of an additional term, proportional to $\boldsymbol {\nabla } \times ((\boldsymbol {\nabla }\times \boldsymbol {B})\times \boldsymbol {B})$ , in Faraday's induction law. The latter feeds back on the magnetic field $B$ at small scales (between the ion and electron inertial scales), requiring very high resolutions in both space and time to properly describe its dynamics. The computational advantage provided by the kinetic lattice Boltzmann (LB) approach is exploited here to develop a new code, the fast lattice-Boltzmann algorithm for MHD experiments (flame). The flame code integrates the plasma dynamics in lattice units coupling two kinetic schemes, one for the fluid protons (including the Lorentz force), the other to solve the induction equation describing the evolution of the magnetic field. Here, the newly developed algorithm is tested against an analytical wave-solution of the dissipative Hall-MHD equations, pointing out its stability and second-order convergence, over a wide range of the control parameters. Spectral properties of the simulated plasma are finally compared with those obtained from numerical solutions from the well-established pseudo-spectral code ghost. Furthermore, the LB simulations we present, varying the Hall parameter, highlight the transition from the MHD to the Hall-MHD regime, in excellent agreement with the magnetic field spectra measured in the solar wind.

Extending the Spectral Difference Method with Divergence Cleaning (SDDC) to the Hall MHD Equations

The Hall Magnetohydrodynamic (MHD) equations are an extension of the standard MHD equations that include the “Hall” term from the general Ohm’s law. The Hall term decouples ion and electron motion physically on the ion inertial length scales. Implementing the Hall MHD equations in a numerical solver allows more physical simulations for plasma dynamics on length scales less than the ion inertial scale length but greater than the electron inertial length. The present effort is an important step towards producing physically correct results to important problems, such as the Geospace Environmental Modeling (GEM) Magnetic Reconnection problem. The solver that is being modified is currently capable of solving the resistive MHD equations on unstructured grids using the spectral difference scheme which is an arbitrarily high-order method that is relatively simple to parallelize. The GEM Magnetic Reconnection problem is used to evaluate whether the Hall MHD equations have been correctly implemented in the solver using the spectral difference method with divergence cleaning (SDDC) algorithm by comparing against the reconnection rates reported in the literature.

Magnetic Field Spectral Evolution in the Inner Heliosphere

Parker Solar Probe and Solar Orbiter data are used to investigate the radial evolution of magnetic turbulence between 0.06 ≲ R ≲ 1 au. The spectrum is studied as a function of scale, normalized to the ion inertial scale d i . In the vicinity of the Sun, the inertial range is limited to a narrow range of scales and exhibits a power-law exponent of, α B = −3/2, independent of plasma parameters. The inertial range grows with distance, progressively extending to larger spatial scales, while steepening toward a α B = −5/3 scaling. It is observed that spectra for intervals with large magnetic energy excesses and low Alfvénic content steepen significantly with distance, in contrast to highly Alfvénic intervals that retain their near-Sun scaling. The occurrence of steeper spectra in slower wind streams may be attributed to the observed positive correlation between solar wind speed and Alfvénicity.

Spatiotemporal dynamics of high-wavenumber turbulence in a basic laboratory plasma

High-spatial resolution observation of high-wavenumber broadband turbulence is achieved by controlling the magnetic field to be relatively low and measuring with a azimuthally arranged multi-channel Langmuir array in a basic laboratory plasma. The observed turbulence consists of narrowband low-frequency fluctuations and broadband high-frequency turbulent fluctuations. The low-frequency fluctuations have a frequency of about 0.7 times the ion cyclotron frequency and a spatial scale of 1/10 of the ion inertial scale. In comparison, high-frequency fluctuations have a higher frequency than the ion cyclotron frequency and spatial scales of 1/10–1/40 of the ion inertial scale. Two-dimensional correlation analysis evaluates the spatial and temporal correlation lengths and reveals that the high-wavenumber broadband fluctuations have turbulent characteristics. The measurements give us further understanding of small scale turbulence in space and fusion plasmas.

Statistical properties of three-dimensional Hall magnetohydrodynamics turbulence

The three-dimensional (3D) Hall magnetohydrodynamics (HMHD) equations are often used to study turbulence in the solar wind. Some earlier studies have investigated the statistical properties of 3D HMHD turbulence by using simple shell models or pseudospectral direct numerical simulations (DNSs) of the 3D HMHD equations; these DNSs have been restricted to modest spatial resolutions and have covered a limited parameter range. To explore the dependence of 3D HMHD turbulence on the Reynolds number Re and the ion-inertial scale di, we have carried out detailed pseudospectral DNSs of the 3D HMHD equations and their counterparts for 3D MHD (di = 0). We present several statistical properties of 3D HMHD turbulence, which we compare with 3D MHD turbulence by calculating (a) the temporal evolution of the energy-dissipation rates and the energy; (b) the wave-number dependence of fluid and magnetic spectra; (c) the probability distribution functions of the cosines of the angles between various pairs of vectors, such as the velocity and the magnetic field; and (d) various measures of the intermittency in 3D HMHD and 3D MHD turbulence.

Hybrid Simulations of a Tangential Discontinuity‐Driven Foreshock Bubble Formation in Comparison With a Hot Flow Anomaly Formation

AbstractHot flow anomalies (HFAs) and foreshock bubbles (FBs) are significant foreshock transients that can accelerate particles and disturb the magnetosphere‐ionosphere system. Yet, their early formation mechanisms are still not fully understood. To investigate the formation of tangential discontinuity (TD)‐driven FBs and HFAs, we use 2‐D local hybrid simulations where a reflected or an injected warm foreshock ion beam can interact with a TD whose half‐thickness is comparable to the ion inertial scale. We show that the foreshock ions perform a partial gyration within, or across, the TD. Bulk motion differences between partially‐gyrating foreshock ions and fluid‐electrons lead to the generation of currents. As the trigger, these foreshock‐driven currents change the magnetic field topology around the TD and force the frozen‐in solar wind plasma to redistribute along with the field lines, shaping the foreshock transient. This confirms a recently proposed kinetic formation model. The extent of the magnetic field direction change across the TD within the foreshock ion gyromotion determines the current profile and thus the type of foreshock transient that forms. For a thin TD, the foreshock ions generate a current that is much stronger on the upstream side than the downstream side, forming an FB with one upstream compressional boundary. For the same foreshock ion gyroradius and magnetic shear, a thick TD yields comparable foreshock‐driven currents on the upstream and downstream sides, forming an HFA with two compressional boundaries. Our study suggests that the TD thickness is one of the factors that determine the formation of FBs and HFAs.

Measurement of the Taylor scale in a magnetized turbulent laboratory plasma wind-tunnel

The fluid Taylor scale is measured in the Bryn Mawr Experiment (BMX) of the Bryn Mawr Plasma Laboratory and examined as a potential dissipation scale of magnetic turbulence within the plasma. We present the first laboratory measurements of the Taylor scale of a turbulent magnetized plasma through multi-point correlations of broadband magnetic fluctuations. From spatial and temporal correlations, respectively, the measured Taylor scales are 2±1 and 3±1 cm. These measurements are on the same order of magnitude as estimated ion dissipation scales within the BMX plasma with ion inertial scales between 1 and 10 cm and ion gyroscales between 0.1 and 1.0 cm. From these measurements, a magnetic Reynolds number can be computed. Since Taylor scale values are determined using multi-point correlations and a Richardson extrapolation technique, an estimate of the magnetic Reynolds number can be found without the added complication of specifying a model of microscopic diffusivity, a parameter often difficult to obtain experimentally.

Properties of Ion-Inertial Scale Plasmoids Observed by the Juno Spacecraft in the Jovian Magnetotail.

We expand on previous observations of magnetic reconnection in Jupiter's magnetosphere by constructing a survey of ion‐inertial scale plasmoids in the Jovian magnetotail. We developed an automated detection algorithm to identify reversals in the Bθ component and performed the minimum variance analysis for each identified plasmoid to characterize its helical structure. The magnetic field observations were complemented by data collected using the Juno Waves instrument, which is used to estimate the total electron density, and the JEDI energetic particle detectors. We identified 87 plasmoids with “peak‐to‐peak” durations between 10 and 300 s. Thirty‐one plasmoids possessed a core field and were classified as flux‐ropes. The other 56 plasmoids had minimum field strength at their centers and were termed O‐lines. Out of the 87 plasmoids, 58 had in situ signatures shorter than 60 s, despite the algorithm's upper limit being 300 s, suggesting that smaller plasmoids with shorter durations were more likely to be detected by Juno. We estimate the diameter of these plasmoids assuming a circular cross section and a travel speed equal to the Alfven speed in the surrounding lobes. Using the electron density inferred by Waves, we contend that these plasmoid diameters were within an order of the local ion‐inertial length. Our results demonstrate that magnetic reconnection in the Jovian magnetotail occurs at ion scales like in other space environments. We show that ion‐scale plasmoids would need to be released every 0.1 s or less to match the canonical 1 ton/s rate of plasma production due to Io.

Relating the Solar Wind Turbulence Spectral Break at the Dissipation Range with an Upstream Spectral Bump at Planetary Bow Shocks

Abstract At scales much larger than the ion inertial scale and the gyroradius of thermal protons, the magnetohydrodynamic (MHD) theory is well equipped to describe the nature of solar wind turbulence. The turbulent spectrum itself is defined by a power law manifesting the energy cascading process. A break in the turbulence spectrum develops near-ion scales, signaling the onset of energy dissipation. The exact mechanism for the spectral break is still a matter of debate. In this work, we use the 20 Hz Mercury Surface, Space Environment, Geochemistry, and Ranging (MESSENGER) magnetic field data during four planetary flybys at different heliocentric distances to examine the nature of the spectral break in the solar wind. We relate the spectral break frequencies of the solar wind MHD turbulence, found in the range of 0.3–0.7 Hz, with the well-known characteristic spectral bump at frequencies ∼1 Hz upstream of planetary bow shocks. Spectral breaks and spectral bumps during three planetary flybys are identified from the MESSENGER observations, with heliocentric distances in the range of 0.3–0.7 au. The MESSENGER observations are complemented by one Magnetospheric Multiscale observation made at 1 au. We find that the ratio of the spectral bump frequency to the spectral break frequency appears to be r- and B-independent. From this, we postulate that the wavenumber of the spectral break and the frequency of the spectral bump have the same dependence on the magnetic field strength ∣B∣. The implication of our work on the nature of the break scale is discussed.

MMS Observations of Double Mid-Latitude Reconnection Ion Beams in the Early Non-Linear Phase of the Kelvin-Helmholtz Instability

The MMS satellites encountered a Kelvin-Helmholtz instability (KHI) period in the early non-linear phase at the post-noon flank magnetopause on 8 Sep 2015. The adjacent magnetosheath was characterized by a pre-dominantly northward Bz > 0 magnetic field with weakly positive in-plane components in a GSM coordinate system. Ion velocity distribution functions indicate at least 17 KH vortex intervals with two typically D-shaped ion beam distributions, commonly associated with reconnection exhausts, that stream in both directions along a mostly northward magnetic field at 350–775 km/s with a median 525 km/s ion beam speed. The counter-streaming ion beams are superposed on a core population of slowly drifting magnetosheath ions with a field-aligned 50–200 km/s speed. Each interval lasted no more than 5.25 s with a median duration of 1.95 s corresponding to in-plane spatial scales 3 < ΔS < 22 di assuming a constant 1 di = 61 km ion inertial scale and a tailward VKH∼258 km/s KH vortex propagation speed along the MMS trajectory. The counter-streaming ions are predominantly observed in the warm KH vortex region between the cold magnetosheath proper and the hot isotropic ion temperature of a low-latitude boundary layer as the MMS constellation traverses a KH vortex. The in-plane spatial scales and the locations of the observed counter-streaming ion beams generally agree with the predictions of twice-reconnected magnetic fields at two mid-latitude reconnection (MLR) regions in a two-fluid three-dimensional numerical simulation previously reported for this KH event. MMS typically recorded a higher phase space density of the fast parallel ion beam that we associate with a tailward reconnection exhaust from the southern MLR (SMLR) and a lower phase space density of the fast anti-parallel ion beam that we associate with a tailward reconnection exhaust from the northern MLR (NMLR) of similar speed. This is either consistent with MMS being closer to the SMLR region than the NMLR region, or that the KHI conditions may have favored reconnection in the SMLR region for the observed in-plane magnetosheath magnetic field as predicted by a two-fluid three-dimensional numerical simulation.

Stability of superthermal strahl electrons in the solar wind

ABSTRACT We present a kinetic stability analysis of the solar wind electron distribution function consisting of the Maxwellian core and the magnetic-field aligned strahl, a superthermal electron beam propagating away from the sun. We use an electron strahl distribution function obtained as a solution of a weakly collisional drift-kinetic equation, representative of a strahl affected by Coulomb collisions but unadulterated by possible broadening from turbulence. This distribution function is essentially non-Maxwellian and varies with the heliospheric distance. The stability analysis is performed with the Vlasov–Maxwell linear solver leopard. We find that depending on the heliospheric distance, the core-strahl electron distribution becomes unstable with respect to sunward-propagating kinetic-Alfvén, magnetosonic, and whistler modes, in a broad range of propagation angles. The wavenumbers of the unstable modes are close to the ion inertial scales, and the radial distances at which the instabilities first appear are on the order of 1 au. However, we have not detected any instabilities driven by resonant wave interactions with the superthermal strahl electrons. Instead, the observed instabilities are triggered by a relative drift between the electron and ion cores necessary to maintain zero electric current in the solar wind frame (ion frame). Contrary to strahl distributions modelled by shifted Maxwellians, the electron strahl obtained as a solution of the kinetic equation is stable. Our results are consistent with the previous studies based on a more restricted solution for the electron strahl.

Identification and characterization of current sheets in collisionless plasmas turbulence

Current sheets forming in an ion-kinetic turbulent collisionless plasma are characterized by utilizing the results of two-dimensional hybrid-kinetic numerical simulations (ions treated as kinetic species using particle-in-cell method and electrons treated as mass-less fluid). For this sake, the algorithm, proposed by Zhdankin et al. [Astrophys. J. 771, 124 (2013)] for the identification and characterization of current sheets forming in magnetohydrodynamic plasma turbulence, is extended to current sheets forming in a much noisier kinetic plasma turbulence. Effects of the algorithm parameters on the identification of current sheets in plasma turbulence are analyzed. Based on this analysis, appropriate values of algorithm parameters for the current sheet identification are chosen. Current sheets are statistically characterized in terms of their main properties—the peak current density, the peak current carrier velocity (mainly electrons), the thickness, the length, and also the aspect ratio (length/thickness). Results of the characterization show that the current sheets have a tendency to thin down to below ion inertial length scale until numerical (grid-resolution based) dissipation stops any further thinning. Simulations that include the physics at the electron inertial scales, e.g., hybrid simulations with electron inertia, are required to fully address the current sheet thinning processes and determination of the final thicknesses of current sheets.

Identification and characterization of current sheets in collisionless plasma turbulence

Current sheets forming in an ion-kinetic turbulent collisionless plasma are characterized by utilizing the results of two-dimensional hybrid-kinetic numerical simulations (ions treated as kinetic species using particle-in-cell method and electrons treated as mass-less fluid). For this sake, the algorithm, proposed by Zhdankin et al. [Astrophys. J. 771, 124 (2013)] for the identification and characterization of current sheets forming in magnetohydrodynamic plasma turbulence, is extended to current sheets forming in a much noisier kinetic plasma turbulence. Effects of the algorithm parameters on the identification of current sheets in plasma turbulence are analyzed. Based on this analysis, appropriate values of algorithm parameters for the current sheet identification are chosen. Current sheets are statistically characterized in terms of their main properties—the peak current density, the peak current carrier velocity (mainly electrons), the thickness, the length, and also the aspect ratio (length/thickness). Results of the characterization show that the current sheets have a tendency to thin down to below ion inertial length scale until numerical (grid-resolution based) dissipation stops any further thinning. Simulations that include the physics at the electron inertial scales, e.g., hybrid simulations with electron inertia, are required to fully address the current sheet thinning processes and determination of the final thicknesses of current sheets.

Subproton-scale Intermittency in Near-Sun Solar Wind Turbulence Observed by the Parker Solar Probe

Abstract High time-resolution solar wind magnetic field data are employed to study statistics describing intermittency near the first perihelion (∼35.6 R ⊙) of the Parker Solar Probe mission. A merged data set employing two instruments on the FIELDS suite enables broadband estimation of higher-order moments of magnetic field increments, with five orders established with reliable accuracy. The duration, cadence, and low noise level of the data permit evaluation of scale dependence of the observed intermittency from the inertial range to deep subproton scales. The results support multifractal scaling in the inertial range, and monofractal but non-Gaussian scaling in the subproton range, thus clarifying suggestions based on data near Earth that had remained ambiguous due to possible interference of the terrestrial foreshock. The physics of the transition to monofractality remains unclear but we suggest that it is due to a scale-invariant population of current sheets between ion and electron inertial scales; the previous suggestion of incoherent kinetic-scale wave activity is disfavored as it presumably leads re-Gaussianization that is not observed.

MMS Observations of Reconnection Separatrix Region in the Magnetotail at Different Distances From the Active Neutral X‐Line

AbstractThe region surrounding the reconnection separatrix consists of many particle and wave transient features (electron, cold and hot ion beams, Hall E&B fields, kinetic Alfvén, LH, etc. waves) whose pattern and parameters may vary depending on the distance from active neutral line. We study nine quick MMS entries into the plasma sheet boundary layer (PSBL) from the tail lobe to address the meso‐scale pattern and other characteristics of phenomena for active separatrix crossings as deduced from particle observations. The outermost thin layer (a fraction of ion inertial scale, di) of low‐density plasma consists of accelerated electron beams and lobe cold ions and displays density depletions (EBL region). It is followed by hot proton beam (PBL region) in which the plasma density grows from lobe‐like towards plasma sheet‐like values; the beam energy‐dispersion is used to estimate the distance from the active neutral line. Thin (usually ≤ di) region containing intense Hall‐like Ez perturbations (HR) usually overlaps with EBL and PBL regions. It often includes correlated B perturbations suggesting the Alfvén wave‐related transport from the reconnection source; the estimated Alfvénic ratio δE/(VA δB) varied between 0.3 and 1.3 in studied examples. The HR is associated with profound plasma property changes, including the heating of cold ion beams in its innermost part, it hosts intense structured field‐aligned currents and intense E‐field fluctuations. Surprisingly, most of abovementioned findings are valid for crossings observed at large distances from the reconnection region (exceeding a few tens Re or >100 di) except for longer time‐scales and larger spatial scales of the pattern.

Evolution of Three-dimensional Coherent Structures in Hall Magnetohydrodynamics

Abstract This work extends the computational model EULAG-MHD to include Hall magnetohydrodynamics (HMHD)—important to explore physical systems undergoing fast magnetic reconnection at the order of the ion inertial length scale. Examples include solar transients along with reconnections in magnetosphere, magnetotail, and laboratory plasmas. The paper documents the results of two distinct sets of implicit large-eddy simulations in the presence and absence of the Hall forcing term, initiated with an unidirectional sinusoidal magnetic field. The HMHD simulation while benchmarking the code also emphasizes the complexity of three-dimensional (3D) evolution over its two-dimensional counterpart. The magnetic reconnections onset significantly earlier in HMHD. Importantly, the magnetic field generated by the Hall term breaks any inherent symmetry, ultimately making the evolution 3D. The resulting 3D reconnections develop magnetic flux ropes (MFRs) and magnetic flux tubes. Projected on the reconnection plane, the ropes and tubes appear as magnetic islands, which later break into secondary islands, and finally coalesce to generate an X-type neutral point. These findings are in agreement with the theory and contemporary simulations of HMHD, and thus verify our extension of the EULAG-MHD model. The second set explores the influence of the Hall forcing on generation and ascend of an MFR from sheared magnetic arcades—a novel scenario instructive in understanding the coronal transients. The rope evolves through intermediate complex structures, ultimately breaking locally because of reconnections. Interestingly, the breakage occurs earlier in the presence of the Hall term, signifying faster dynamics leading to magnetic topology favorable for reconnections.

MMS Observations of Intense Whistler Waves Within Earth's Supercritical Bow Shock: Source Mechanism and Impact on Shock Structure and Plasma Transport

AbstractThe properties of whistler waves near lower‐hybrid frequencies within Earth's quasi‐perpendicular bow shock are examined using data from the Magnetospheric Multiscale (MMS) mission. These waves appear as right‐hand polarized wave packets propagating upstream obliquely to the magnetic field and shock normal with phase speeds from a few hundred up to 1,600 km/s. The wavelengths are near the ion inertial length scale (λ∼ 0.3–1.3 λi). Detailed analysis finds characteristics consistent with the modified two‐stream instability mechanism driven by the reflected ion and electron drift. Correlations between wave and electron anisotropy variations reveal that the whistlers are affecting electron dynamics and thus their perpendicular and parallel temperatures. The electron signatures are explainable via the interaction of magnetized electrons in the whistler induced nonmonotonic magnetic fields. These waves have intense magnetic fields (δB/B∘∼ 0.1–1) and carry sizable currents that are a significant fraction of the thermal current (|J/Jvte|∼ 0.1–0.5). The whistler‐induced currents and the electron anisotropies are sufficiently large to respectively excite high‐frequency (HF) electrostatic (>100 Hz) and HF whistler waves (f∼ 0.1–0.5 fce). Energy dissipation J·E from whistlers at 30 Hz and below range from a few thousandths to few hundredths of μW/m3. Comparisons reveal that plasma energy is converted to wave energy in the foot, whereas wave energy gets dissipated into the plasma in the ramp, where irreversible heating occurs. These observed features are indicative of an intricate coupling between small‐scale interaction processes and larger‐scale structure transpiring within the layer. Such a characterization is only made possible now with the MMS high‐time‐resolution measurements.

Electron-only Reconnection in Kinetic-Alfvén Turbulence

Abstract We study numerically small-scale reconnection events in kinetic, low-frequency, quasi-2D turbulence (termed kinetic-Alfvén turbulence). Using 2D particle-in-cell simulations, we demonstrate that such turbulence generates reconnection structures where the electron dynamics do not couple to the ions, similarly to the electron-only reconnection events recently detected in the Earth’s magnetosheath by Phan et al. Electron-only reconnection is thus an inherent property of kinetic-Alfvén turbulence, where the electron current sheets have limited anisotropy and, as a result, their sizes are smaller than the ion inertial scale. The reconnection rate of such electron-only events is found to be close to 0.1.