Terrestrial Magnetosphere Research Articles

Numerical simulations of temperature anisotropy instabilities stimulated by suprathermal protons

The new in situ measurements of the Solar Orbiter mission contribute to the knowledge of the suprathermal populations in the solar wind, especially of ions and protons whose characterization, although still in the early phase, seems to suggest a major involvement in the interaction with plasma wave fluctuations. Recent studies point to the stimulating effect of suprathermal populations on temperature anisotropy instabilities in the case of electrons already being demonstrated in theory and numerical simulations. Here, we investigate anisotropic protons, addressing the electromagnetic ion-cyclotron (EMIC) and the proton firehose (PFH) instabilities. Suprathermal populations enhance the high-energy tails of the Kappa velocity (or energy) distributions measured in situ, enabling characterization by contrasting to the quasi-thermal population in the low-energy (bi-)Maxwellian core. We use hybrid simulations to investigate the two instabilities (with ions or protons as particles and electrons as fluid) for various configurations relevant to the solar wind and terrestrial magnetosphere. The new simulation results confirm the linear theory and its predictions. In the presence of suprathermal protons, the wave fluctuations reach increased energy density levels for both instabilities and cause faster and/or deeper relaxation of temperature anisotropy. The magnitude of suprathermal effects also depends on each instability's specific (initial) parametric regimes. These results further strengthen the belief that wave-particle interactions govern space plasmas. These provide valuable clues for understanding their dynamics, particularly the involvement of suprathermal particles behind the quasi-stationary non-equilibrium states reported by in situ observations.

Field‐Aligned Current Structures During the Terrestrial Magnetosphere's Transformation Into Alfvén Wings and Recovery

AbstractOn 24 April 2023, a Coronal Mass Ejection event caused the solar wind to become sub‐Alfvénic, leading to the development of an Alfvén Wing configuration in the Earth's magnetosphere. Alfvén Wings have previously been observed as cavities of low flow around moons in Jupiter's and Saturn's magnetospheres, but the observing spacecraft did not have the ability to directly measure the Alfvén Wings' current structures. Through in situ measurements made by the Magnetospheric Multiscale spacecraft, the 24 April event provides us with the first direct measurements of current structures during an Alfvén Wing configuration. These structures are observed to be significantly more anti‐field‐aligned and electron‐driven than the typical diamagnetic magnetopause current, indicating the disruption caused to the magnetosphere current system by the Alfvén Wing formation. The magnetopause current is then observed to recover more of its typical, perpendicular structure during the magnetosphere's recovery from the Alfvén Wing formation.

The Relationship Between Large dB/dt and Field‐Aligned Currents During Five Geomagnetic Storms

AbstractDuring periods of increased geomagnetic activity, perturbations within the terrestrial magnetosphere are known to induce currents within conducting materials, at the surface of Earth through rapid changes in the local magnetic field over time (dB/dt). These currents are known as geomagnetically induced currents and have potentially detrimental effects on ground based infrastructure. In this study we undertake case studies of five geomagnetic storms, analyzing a total of 19 days of 1‐s SuperMAG data in order to better understand the magnetic local time (MLT) distribution, size, and occurrence of “spikes” in dB/dt, with 131,447 spikes in dB/dt exceeding 5 nT/s identified during these intervals. These spikes were concentrated in clusters over three MLT sectors: two previously identified pre‐midnight and dawn region hot‐spots, and a third, lower‐density population centered around 12 MLT (noon). The noon spike cluster was observed to be associated with pressure pulse impacts, however, due to incomplete magnetometer station coverage, this population is not observed for all investigated storms. The magnitude of spikes in dB/dt are determined to be greatest within these three “hot‐spot” locations. These spike occurrences were then compared with field‐aligned current (FAC) data, provided by the Active Magnetospheric Planetary Electrodynamic Response Experiment. Spikes are most likely to be co‐located with upward FACs (56%) rather than downward FACs (30%) or no FACs (14%).

Unveiling the 3D structure of magnetosheath jets

ABSTRACT Magnetosheath jets represent localized enhancements in dynamic pressure observed within the magnetosheath. These energetic entities, carrying excess energy and momentum, can impact the magnetopause and disrupt the magnetosphere. Therefore, they play a vital role in coupling the solar wind and terrestrial magnetosphere. However, our understanding of the morphology and formation of these complex, transient events remains incomplete over two decades after their initial observation. Previous studies have relied on oversimplified assumptions, considering jets as elongated cylinders with dimensions ranging from $0.1\, R_{\rm E}$ to $5\, R_{\rm E}$ (Earth radii). In this study, we present simulation results obtained from Amitis, a high-performance hybrid-kinetic plasma framework (particle ions and fluid electrons) running in parallel on graphics processing units (GPUs) for fast and more environmentally friendly computation compared to CPU-based models. Considering realistic scales, we present the first global, three-dimensional (3D in both configuration and velocity spaces) hybrid-kinetic simulation results of the interaction between solar wind plasma and the Earth. Our high-resolution kinetic simulations reveal the 3D structure of magnetosheath jets, showing that jets are far from being simple cylinders. Instead, they exhibit intricate and highly interconnected structures with dynamic 3D characteristics. As they move through the magnetosheath, they wrinkle, fold, merge, and split in complex ways before a subset reaches the magnetopause.

Field Line Curvature (FLC) Scattering in the Dayside Off‐Equatorial Minima Regions

AbstractMagnetic field line curvature (FLC) scattering is an effective mechanism for collisionless particle scattering. In the terrestrial magnetosphere, the FLC scattering plays an essential role in shaping the outer boundary of protons radiation belt, the rapid decay of ring current, and the formation of proton isotropic boundary (IB). However, previous studies have yet to adequately investigate the influence of FLC scattering on charged particles in the Earth's dayside magnetosphere, particularly in the off‐equatorial magnetic minima regions. This study employs T89 magnetic field model to investigate the impacts of FLC scattering on ring current protons in the dayside magnetosphere, with a specific focus on the off‐equatorial minimum regions. We analyze the spatial distributions of single and dual magnetic minima regions, adiabatic parameter, and pitch angle diffusion coefficients due to FLC scattering as functions of Kp. The results show that the effects of FLC scattering are significant not only on the dusk and dawn sides but also in the off‐equatorial minima regions on the noon. Additionally, we investigate the role of dipole tilt angle in the hemispheric asymmetry of FLC scattering effects. The dipole tilt angle controls the overall displacement of the dayside magnetosphere, resulting in different FLC scattering effects in the two hemispheres. Our study holds significance for understanding the FLC scattering effects in the off‐equatorial region of Earth's dayside magnetosphere and for constructing a more accurate dynamic model of particles.

Nonlinear electron-acoustic waves in non-Maxwellian plasma: application in terrestrial magnetosphere

Nonlinear electron-acoustic waves in non-Maxwellian plasma: application in terrestrial magnetosphere

High-energy particle observations from the Moon.

The Moon is a unique natural laboratory for the study of the deep space plasma and energetic particles environment. During more than 3/4 of its orbit around the Earth it is exposed to the solar wind. Being an unmagnetized body and lacking a substantial atmosphere, solar wind and solar energetic particles bombardthe Moon's surface, interacting with the lunar regolith and the tenuous lunar exosphere. Energetic particles arriving at the Moon's surface can be absorbed, or scattered, or can remove another particle from the lunar regolith by sputtering or desorption. A similar phenomenon occurs also with the galactic cosmic rays, which have fluxes and energy spectra representative of interplanetary space. During the remaining part of its orbit the Moon crosses the tail of the terrestrial magnetosphere. It then provides the opportunity to study in-situ the terrestrial magnetotail plasma environment as well as atmospheric escape from the Earth's ionosphere, in the form of heavy ions accelerated and streaming downtail. The lunar environment is thus a unique natural laboratory for analysing the interaction of the solar wind, the cosmic rays and the Earth's magnetosphere with the surface, the immediate subsurface, and the surface-bounded exosphere of an unmagnetized planetary body. This article is part of a discussion meeting issue 'Astronomy from the Moon: the next decades (part 2)'.

Signatures of Open Magnetic Flux in Jupiter's Dawnside Magnetotail

AbstractJupiter's magnetosphere exhibits notable distinctions from the terrestrial magnetosphere. The structure and dynamics of Jupiter's dawnside magnetosphere can be characterized as a competition between internally driven sunward flow and solar wind‐driven tailward flow. During the prime mission, Juno acquired extensive data from dawn to midnight, sampling the magnetodisc and higher latitude regions. Numerical moments from the Jovian Auroral Distributions Experiment (JADE‐I) plasma (ion) instrument revealed a mid‐latitude region of anticorotational (−vϕ) flow. While the magnitude of the flow is subject to uncertainty due to low count rates in these rarefied regions, we demonstrate in the raw JADE‐I data that the sign of vϕ is a robust measurement. Global Grid Agnostic Magnetohydrodyamics for Extended Research Applications simulations show a similar region of strongly reduced flow in proximity to open field lines. Additionally, we use Jupiter Energetic‐particle Detector Instrument integral moments to determine the Hen+/H+ ratio (where n refers to He+ or He++) and show that a transition to solar wind‐like composition occurs in the same region as the anticorotational flow. We conclude that the global simulations are consistent with the Juno data, where the simulations show a crescent of open magnetic flux that is bounded by the magnetodisc and a closed high‐latitude polar region (nominally the polar cap), which is never observed in the terrestrial magnetosphere. The distinct distribution of open flux in Jupiter's dawnside magnetosphere suggests the significance of planetary rotation and may represent a characteristic feature of rotating giant magnetospheres for future exploration.

Design and construction of the near-earth space plasma simulation system of the Space Plasma Environment Research Facility

Our earth is immersed in the near-earth space plasma environment, which plays a vital role in protecting our planet against the solar-wind impact and influencing space activities. It is significant to investigate the physical processes dominating the environment, for deepening our scientific understanding of it and improving the ability to forecast the space weather. As a crucial part of the National Major Scientific and Technological Infrastructure–Space Environment Simulation Research Infrastructure (SESRI) in Harbin, the Space Plasma Environment Research Facility (SPERF) builds a system to replicate the near-earth space plasma environment in the laboratory. The system aims to simulate the three-dimensional (3-D) structure and processes of the terrestrial magnetosphere for the first time in the world, providing a unique platform to reveal the physics of the 3-D asymmetric magnetic reconnection relevant to the earth's magnetopause, wave–particle interaction in the earth's radiation belt, particles’ dynamics during the geomagnetic storm, etc. The paper will present the engineering design and construction of the near-earth space plasma simulation system of the SPERF, with a focus on the critical technologies that have been resolved to achieve the scientific goals. Meanwhile, the possible physical issues that can be studied based on the apparatus are sketched briefly. The earth-based system is of great value in understanding the space plasma environment and supporting space exploration.

Betatron Acceleration of Suprathermal Electrons Upstream of the Martian Bow Shock

Betatron acceleration, a plasma process obtaining particle energy in the perpendicular direction but reserving energy in the field-aligned direction, is the consequence of magnetic strength enhancement when the first adiabatic invariant is conserved. Such process has been widely reported in the terrestrial magnetosphere but is barely reported in other planetary environments. Here, based on the in situ measurements from the Mars Atmosphere and Volatile Evolution mission, we report two events of betatron acceleration upstream of the Martian bow shock. In both events, betatron accelerations increase the fluxes of suprathermal electrons. The acceleration processes in these events are quantitatively reproduced with an analytical model. Gratifyingly, we find the acceleration factors derived from the analytical model are well consistent with the observations of magnetic strength enhancement. These results for the first time show that the betatron acceleration is an active upstream of the Martian bow shock and is very useful to help us understand the generation of energetic electrons in the Martian environment.

Three-dimensional magnetic reconnection in complex multiple X-point configurations in an ancient solar–lunar terrestrial system

Magnetic reconnection processes in three-dimensional (3D) complex field configurations have been investigated in different magneto-plasma systems in space, laboratory, and astrophysical systems. Two-dimensional (2D) features of magnetic reconnection have been well developed and applied successfully to systems with symmetrical property, such as toroidal fusion plasmas and laboratory experiments with an axial symmetry. But in asymmetric systems, the 3D features are inevitably different from those in the 2D case. Magnetic reconnection structures in multiple celestial body systems, particularly star–planet–Moon systems, bring fresh insights to the understanding of the 3D geometry of reconnection. Thus, we take magnetic reconnection in an ancient solar–lunar terrestrial magneto-plasma system as an example by using its crucial parameters approximately estimated already and also some specific applications in pathways for energy and matter transports among Earth, ancient Moon, and the interplanetary magnetic field (IMF). Then, magnetic reconnection of the ancient lunar–terrestrial magnetospheres with the IMF is investigated numerically in this work. In a 3D simulation for the Earth–Moon–IMF system, topological features of complex magnetic reconnection configurations and dynamical characteristics of magnetic reconnection processes are studied. It is found that a coupled lunar-terrestrial magnetosphere is formed, and under various IMF orientations, multiple X-points emerge at distinct locations, showing three typical magnetic reconnection structures in such a geometry, i.e., the X-line, the triple current sheets, and the A–B null pairs. The results can conduce to further understanding of reconnection physics in 3D for plasmas in complex magnetic configurations, and also a possible mechanism for energy and matters transport in evolutions of similar astrophysical systems.

Statistical Comparison of Various Dayside Magnetopause Reconnection X‐Line Prediction Models

AbstractReconnection at Earth's magnetopause drives magnetospheric convection and provides mass and energy input into the magnetosphere/ionosphere system thereby driving the coupling between solar wind and terrestrial magnetosphere. Despite its importance, the factors governing the location of dayside magnetopause reconnection are not well understood. Though a few models can predict X‐line locations reasonably well, the underlying physics is still unresolved. In this study we present results from a comparative analysis of 274 magnetic reconnection events as observed by the Magnetospheric Multiscale (MMS) mission to determine what quantities affect the accuracy of such models and are most strongly associated with the occurrence of dayside magnetopause reconnection. We also attempt to determine under what upstream solar wind conditions each global X‐line model becomes least reliable.

Electron Rolling‐Pin Distributions Preceding Dipolarization Fronts

AbstractElectron rolling‐pin distribution (RPDs), characterized by triple peaks at pitch angles 0°, 90°, and 180°, have recently been discovered in the terrestrial magnetosphere. Since the RPDs' formation is typically attributed to local betatron acceleration, RPDs have been believed to appear primarily inside strong magnetic regions, such as flux pileup regions (FPRs) behind dipolarization fronts (DFs). Different from such expectation, in this study we present unique observations of RPDs made by Cluster spacecraft in the terrestrial magnetotail, showing that RPDs are present inside weak magnetic field regions ahead of the DFs but absent in the strong magnetic field region behind them. The presence of RPD prior to the fronts may be atttributed to the combined effect of global betatron acceleration, global Fermi acceleration, and frontward transport driven by magnetic gradient drift. The atypical features of RPDs are important for fully understanding electron acceleration and transport in the magnetosphere.

Estimation of the Error in the Calculation of the Pressure‐Strain Term: Application in the Terrestrial Magnetosphere

AbstractCalculating the pressure‐strain terms has recently been performed to quantify energy conversion between the bulk flow energy and the internal energy of plasmas. It has been applied to numerical simulations and satellite data from the Magnetospheric MultiScale Mission. The method requires spatial gradients of the velocity and the use of the full pressure tensor. Here we present a derivation of the errors associated with calculating the pressure‐strain terms from multi‐spacecraft measurements and apply it to previously studied examples of magnetic reconnection at the magnetopause and the magnetotail. The errors are small in a dense magnetosheath event but much larger in the more tenuous magnetotail. This is likely due to larger counting statistics in the dense plasma at the magnetopause than in the magnetotail. The propagated errors analyzed in this work are important to understand uncertainties of energy conversion measurements in space plasmas and have applications to current and future multi‐spacecraft missions.

The effect of heavy ions on the dispersion properties of kinetic Alfvén waves in astrophysical plasmas

Context.Spacecraft measurements have shown Kinetic Alfvén Waves propagating in the terrestrial magnetosphere at lower wave-normal angles than predicted by linear Vlasov theory of electron-proton plasmas. To explain these observations, it has been suggested that the abundant heavy ion populations in this region may have strong, non-trivial effects that allow Alfvénic waves to acquire right-handed polarization at lower angles with respect to the background magnetic field, as in the case of typical electron-proton plasma.Aims.We study the dispersion properties of Alfvénic waves in plasmas with stationary phase-space distribution functions with different heavy ion populations. Our extensive numerical analysis has allowed us to quantify the role of the heavy ion components on the transition from the left-hand polarized electromagnetic ion-cyclotron (EMIC) mode to the right-hand polarized kinetic Alfvén wave (KAW) mode.Methods.We used linear Vlasov-Maxwell theory to obtain the dispersion relation for oblique electromagnetic waves. The dispersion relation of Alfvén waves was obtained numerically by considering four different oxygen ion concentrations ranging between 0.0 and 0.2 for all propagation angles, as a function of both the wavenumber and the plasma beta parameter.Results.The inclusion of the heavy O+ions is found to considerably reduce the transition angle from EMIC to KAW both as a function of the wave number and plasma beta. With increasing O+concentrations, waves become more damped in specific wavenumber regions. However, the inclusion of oxygen ions may allow weakly damped KAW to effectively propagate at smaller wave-normal angles than in the electron-proton case, as suggested by observations.

Parameter Study of Geomagnetic Storms and Associated Phenomena: CME Speed De-Projection vs. In Situ Data

In this study, we give correlations between the geomagnetic storm (GS) intensity and parameters of solar and interplanetary (IP) phenomena. We also perform 3D geometry reconstructions of geo-effective coronal mass ejections (CMEs) using the recently developed PyThea framework and compare on-sky and de-projected parameter values, focusing on the reliability of the de-projection capabilities. We utilize spheroid, ellipsoid and graduated cylindrical shell models. In addition, we collected a number of parameters of the GS-associated phenomena. A large variation in 3D de-projections is obtained for the CME speeds depending on the selected model for CME reconstruction and observer subjectivity. A combination of fast speed and frontal orientation of the magnetic structure upon its arrival at the terrestrial magnetosphere proves to be the best indicator for the GS strength. More reliable estimations of geometry and directivity, in addition to de-projected speeds, are important for GS forecasting in operational space weather schemes.

Multi‐Instrument Observations of the Effects of a Solar Wind Pressure Pulse on the High Latitude Ionosphere: A Detailed Case Study of a Geomagnetic Sudden Impulse

AbstractThe effects of a solar wind pressure pulse on the terrestrial magnetosphere have been observed in detail across multiple datasets. The communication of these effects into the magnetosphere is known as a positive geomagnetic sudden impulse (+SI), and are observed across latitudes and different phenomena to characterize the propagation of +SI effects through the magnetosphere. A superposition of Alfvén and compressional propagation modes are observed in magnetometer signatures, with the dominance of these signatures varying with latitude. For the first time, collocated lobe reconnection convection vortices and region 0 field aligned currents are observed preceding the +SI onset, and an enhancement of these signatures is observed as a result of +SI effects. Finally, cusp auroral emission is observed collocated with the convection and current signatures. For the first time, simultaneous observations across multiple phenomena are presented to confirm models of +SI propagation presented previously.

Ion Acceleration at the Quasi‐Parallel Shock: The Source Distributions of the Diffuse Ions

AbstractThe terrestrial bow shock is the boundary that slows and diverts the supermagnetosonic solar wind around the terrestrial magnetosphere by converting the kinetic energy of the solar wind into thermal and magnetic energy. Shock fronts are an important acceleration site for ions and electrons in collisionless plasmas, and are responsible for much of the particle acceleration in solar, planetary, and astrophysical regions. One of the fundamental outstanding questions of ion acceleration at shocks for which the upstream magnetic field is nearly aligned with the shock normal (i.e., quasi‐parallel shocks) is which portion of the incoming solar wind ion distribution ultimately becomes the seed population that is subsequently accelerated to high energies. This study discusses distribution functions of protons and alpha particles observed by the HPCA and FPI instruments onboard the MMS satellites during a crossing of the quasi‐parallel bow shock. The bow shock transition from the downstream region into the upstream solar wind shows the occasional presence of reflected ions and a population of 90° pitch angle ions in the shock ramp consistent with shock drift accelerated ions. Both populations contribute to the seed population of the shock accelerated ions known as the diffuse ion population.

Kinetic Scale Magnetic Reconnection with a Turbulent Forcing: Particle-in-cell Simulations

Turbulent magnetic reconnection has been observed by spacecraft to occur commonly in terrestrial magnetosphere and the solar wind, providing a new scenario of kinetic scale magnetic reconnection. Here by imposing a turbulent forcing on ions in particle-in-cell simulations, we simulate kinetic scale turbulent magnetic reconnection. We find formation of fluctuated electric and magnetic fields and filamentary currents in the diffusion region. Reconnection rate does not change much compared to that in laminar Hall reconnection. At the X-line, the electric and magnetic fields both exhibit a double power-law spectrum with a spectral break near local lower-hybrid frequency. The energy conversion rate is high in turbulent reconnection, leading to significant electron acceleration at the X-line. The accelerated electrons form a power-law spectrum in the high energy range, with a power-law index of about 3.7, much harder than one can obtain in laminar reconnection.

Response of Electric Field in Terrestrial Magnetosphere to Interplanetary Shock

Electric field impulses generated by interplanetary shocks can cause a series of dynamic processes in the Earth’s magnetosphere and were previously explained by either fast-mode wave propagation or flow related to compression of the magnetopause. Based on a Space Weather Modeling Framework simulation, we suggest a new scenario in which the evolution of the impulse is due to both the propagation of the fast-mode wave and the compression of the magnetopause, which can explain the simulation and observations in previous related studies. The onset of the electric field impulse is determined by the propagation of the fast-mode wave in the magnetosphere while the peak of the impulse is determined by the propagation of the compression of the magnetopause. The new understanding of the impulse is important for the generation of subsequent ultralow frequency waves through the coupling of the fast-mode to Alfvén waves and field line resonances and related radiation-belt electron acceleration.