Abstract
Magnetosonic (MS) waves, i.e., ion Bernstein mode waves, are one of the common plasma waves in the Earth’s magnetosphere, which are important for regulating charged particle dynamics. How MS waves propagate in the magnetosphere is critical to understanding the global distribution of the waves, but it remains unclear. Although previous studies present that MS waves can be reflected by fine-scale density structures, the dissipation of waves by background plasma has long been neglected. In this study, we perform one-dimensional (1-D) particle-in-cell (PIC) simulations to study the propagation of MS waves through density structures, where both absorption and reflection have been included. We find that absorption is as important as reflection when considering the propagation of MS waves through density structures, and both of them are strongly dependent on the shape of density structures. Specifically, the reflectivity of MS waves is positively and negatively correlated with the height and width of density structures, respectively, while the absorptivity of MS waves has a positive correlation with both the height and width of density structures. Our study demonstrates the significance of absorption during the propagation of MS waves, which may help better understand the distribution of MS waves in the Earth’s magnetosphere.
1 Introduction
Magnetosonic (MS) waves, also known as ion Bernstein mode waves, are one of the intense electromagnetic emissions observed in the Earth’s inner magnetosphere []. These waves were first detected by the OGO 3 satellite and named “equatorial noise” due to their occurrence within about of the magnetic equator [–]. Recent observations made by Cluster and THEMIS satellites [,] have shown that MS waves can occur both inside and outside the plasmasphere near the magnetic equator. The waves are excited at harmonics of the proton gyrofrequency [] and at large (∼) wave normal angles [,] by ring velocity distributions of ring current protons [,]. MS waves play a significant role in regulating the dynamics of charged particles in the Earth’s magnetosphere [–]. They have been proposed as a candidate for accelerating ∼100 keV electrons up to relativistic energies in the outer radiation belt []. Additionally, scattering by MS waves may explain the formation of butterfly distributions of radiation belt electrons [,,–]. Furthermore, MS waves can effectively energize the background cold protons and electrons [,].
Understanding the propagation of MS waves in the Earth’s magnetosphere is crucial in comprehending the global distributions of these waves and their impact on energy transfer among different particle populations. Satellite observations indicate that the occurrence rate of MS waves strongly depends on the magnetic local time (MLT) outside the plasmapause, but remains nearly uniform inside the plasmapause []. This coincides with the scenario that MS waves are initially generated outside the plasmasphere in the noon and dusk sectors and then propagate both outward and inward, crossing the plasmapause and migrating globally over MLT [,]. Moreover, the occurrence rate and intensity of MS waves outside the plasmapause are higher than inside it [] and the majority of MS waves inside the plasmapause have lower frequencies than the local proton cyclotron frequency [,], making radial propagation the most plausible explanation [,,].
The propagation of MS waves is strongly influenced by the inhomogeneous background plasma density. By performing one-dimensional (1-D) full wave simulations with the finite difference time domain (FDTD) method, Liu et al. [] have found that MS waves can propagate deep into the plasmasphere with only a small fraction of the MS wave power being reflected by the plasmapause. Instead, the fine-scale density structures near the outer edge of the plasmapause can effectively reflect MS waves. Such fine-scale density structures have been widely observed in the Earth’s magnetosphere [,–]. However, previous simulations have also revealed that MS waves can be significantly damped by the background cold plasma [,,], which was neglected in the study of Liu et al. [] due to the limitation of their model. Therefore, we utilize the 1-D PIC model to simulate the propagation of MS waves through the fine-scale density structures, where both absorption and reflection have been considered. We have also quantified the reflectivity and absorptivity of MS waves passing through the density structure and investigated their dependences on the shape of the density structure.
2 Simulation model and method
In this study, we employ a 1-D PIC simulation model to investigate the effects of density structures on the propagation of MS waves. The background magnetic field is directed along the z-axis, and the wave vector of MS waves is lying in the - plane. Here, the simulation box is along the wave vector (or the propagating direction), which is defined as the r direction. This model includes full three-dimensional electromagnetic fields and velocities but only allows spatial variations in the r direction. The periodic boundary conditions are adopted. The units of time and space are the inverse of the proton gyrofrequency , and the proton inertial length , respectively. The plasma system only consists of background protons and electrons which are denoted by subscripts “p” and “e” hereafter. Both protons and electrons satisfy a Maxwellian velocity distribution and have the same temperature . To reduce computational costs, the mass ratio of proton to electron is set to 1600, and the ratio of light speed to the Alfven speed is set to 20. The simulation domain with a length of is divided equally into 30000 grids. The average number of superparticles in each grid is approximately 100 for each species, and the time step is set to .
The angle between the wave vector and the background magnetic field is defined as , i.e., the wave normal angle. Here we will consider two categories of MS waves: perpendicular () and quasi-perpendicular waves (). For each run, we initially pump the monochromatic MS wave from the left boundary to ( is the wavelength) by assigning fluctuating wave fields on each grid and fluctuating bulk velocity to each particle in the form of ( is the related parameter and is the wave number, respectively) along the r direction. Based on the dispersion relation of MS waves in a cold plasma, we can obtain the wave fields by the following relations:where P, D, R, L, and S are the Stix parameters [,] and is the refraction index in Eqs 1–6. Besides, the corresponding bulk velocities of protons and electrons are given by Eqs 7–9:where qj, mj, and Ωcj denote the charge, mass, and cyclotron frequency of the j-component of plasma (j indicates p or e), respectively. In each run, the is set to , and other parameters can be calculated according to the above relations.
For convenience, the density structure is assumed as the sinusoidal variation of density, so the plasma density as a function of r is given by Eq. 10:where n0 is the ambient number density (outside the density structure), and H is the height of the density structure which is the ratio of the peak density to the ambient density . and denote the location and width of the density structure, respectively. The location of the density structure is fixed at in each run.
In the following section, we will present the simulation results of three runs in detail: Run 1 with , , and , Run 2 with and no structure, and Run 3 with , , and . To show how we initialize the simulation model, we present the spatial profiles of (a) plasma density, (b) wave fields, bulk velocities of (c) electrons, and (d) protons at for Run 1 in Figure 1. There is a density structure located at with the width and the height (Figure 1A). The MS waves are launched within the region of (i.e., ). For the perpendicular MS wave, there are only one component of fluctuating magnetic fields (δBz) and two components of fluctuating electric fields (δEx and δEy) (Figure 1B). The bulk velocities of protons and electrons are shown in Figures 1C, D, respectively. Although the corresponding fluctuating density is not initialized, the density fluctuation will be self-consistently coupled to the MS wave very quickly in the PIC model
FIGURE 1
3 Simulation results
3.1 Perpendicular waves:
Figure 2 displays the propagation of the MS waves in Run 1 with a density structure and in Run 2 without it. Without the density structure, the MS waves propagate along the r direction (x-axis) with a constant speed of , which is consistent with the linear dispersion relation (Figures 2C, D). As expected, we can find that the wave amplitude is gradually decreasing during the propagation due to the dissipation caused by background plasma. In Run 1, the MS wave encounters the density structure at about , and then there appears the weak backward propagating MS wave emitting from the left boundary of the density structure. Such a wave is just reflected by the density structure. Besides, the phase velocity inside the density structure becomes smaller than that outside the structure due to the enhanced plasma density (Figure 2A). Meanwhile, the wave fields and become larger and smaller inside the density structure (Figures 2A, B), respectively, consistent with the results of Liu et al. [].
FIGURE 2
Figure 3A shows the time evolution of the energy of MS wave fields (orange line), the kinetic energy of charged particles (green line), and their sum (blue line) for Run 1. First of all, the total energy in our model is well conserved within a 0.06% margin of error. Secondly, it is clearly shown that the wave energy decreases as the particle kinetic energy increases with time, suggesting the wave dissipation takes effect during the propagation. We further display the time evolution of kinetic energies of protons and electrons inside the density structure in Figures 3B, C, respectively, where the vertical dashed lines mark the tie when the wave enters and leaves this structure. For both protons and electrons, their parallel kinetic energy remains nearly unchanged, since the perpendicular MS wave does not have the parallel electric field and cannot energize particles in the parallel direction. The perpendicular kinetic energy of protons rapidly increases as the wave enters the density structure mainly due to the pick-up process. After the wave leaves the density structure, the rapidly drops, but there is still the obvious net increase of kinetic energy (Figure 3B), which is caused by the perpendicular heating by the MS wave (also see Supplementary Material). However, the net increase of kinetic energy for electrons is not significant (Figure 3C).
FIGURE 3
To quantify the effects of the density structure on the propagation of MS waves, we need to calculate the reflectivity , absorptivity , and transmissivity of MS waves. Here, the , , are estimated based on the following formulas:where , , and are the Poynting flux of the incident, reflected, and transmitted waves, respectively, and the corresponding amplitudes are , , and . The waveform of transmitted waves is shown in Figure 4C, which is directly obtained by recording the time series of at the right boundary of the density structure. The eight clearest wavelengths are selected by marking the peaks and troughs with red dots, and then the amplitude of transmitted wave is given by their average. However, the reflected and incident MS waves cannot be directly obtained, since the incident and reflected waves are mixed at the left boundary in Run 1. Instead, we use the recorded waveform at the left boundary of the density structure in Run 2 to represent the incident wave (Figure 4A), and the amplitude is the average of those marked points. Then, the waveform of reflected waves is obtained by subtracting the reflected waveform shown in Figure 4A from the mixed waveform recorded at the left boundary in Run 1, and the amplitude is calculated as above. For Run 1, the , , and are estimated as 0.01783, 0.00392, and 0.01683, so the , , and are 4.83%, 6.07%, and 89.1% according to Eqs 11–13. It is clearly found that the absorption of MS wave by the density structure is even more significant than the reflection in this case.
FIGURE 4
Furthermore, we have also studied the effects of the shape of density structure on the reflection and absorption of MS wave by considering the various widths and heights . Figure 5 summarizes the simulation results (32 additional runs). First of all, there is a clear trend that both the and are positively correlated with the height of density structure (Figures 5A, B). Just as expected, the corresponding transmissivity must be decreasing with the height (Figure 5C). However, with the increasing , the decreases while the increases (Figures 5D, E). Their opposite correlation with the width of density structure results in the particular dependence of on the . With the increasing , the first increases and then decreases, leading to a maximum on a certain value of (Figure 5F). Moreover, comparing the with , we can find that the absorptivity is comparable to or even larger than the reflectivity in most cases, suggesting the absorption of waves should not be neglected when modeling the propagation of MS wave through the density structure.
FIGURE 5
3.2 Quasi-perpendicular waves:
Besides the strictly perpendicular waves, we have also studied the propagation of MS waves with through the density structure due to the different interactions with charged particles []. Figure 6 shows the spatial-temporal evolution of wave fields for Run 3. Similar to Run1, the reflected wave immediately shows up when the MS wave encounters the left boundary of the density structure at , and the transmitted waves become weaker after the MS wave leaves the density structure.
FIGURE 6
Figure 7A displays the time history of the wave energy and particle kinetic energy for Run 3. The total energy (blue line) is well conserved in this simulation model, and the energy exchange between MS wave and plasma can be easily found. Figures 7B, C present the time evolution of the kinetic energies of protons and electrons inside the density structure, respectively. The perpendicular kinetic energy of protons and electrons inside the structure first increases as the wave enters the density structure and then decreases after the wave leaves, with a weak net increase. However, the parallel kinetic energy of electrons gradually increases when the wave passes through the structure and then remains nearly constant after the wave leaves. This is due to the efficient parallel heating of electrons by quasi-parallel MS waves (Supplementary Figure S2).
FIGURE 7
We further investigate the effects of the shape of density structure on the reflection and absorption of the MS wave with . The calculation method of , , and is the same as above. The simulation results marked by red dots and solid lines are summarized in Figure 8. Similar to the perpendicular MS wave, the and increase with the increasing , resulting in the decreasing . Then, the and the have a positive correlation and an anti-correlation with the , respectively, leading to a maximum of at a certain . Although the trends of , , and with the and are very similar to the perpendicular cases, the reflectivity and the absorptivity of the quasi-perpendicular MS waves are smaller than those of the MS waves with , which may imply that the quasi-perpendicular MS waves can propagate further than those perpendicular waves in the magnetosphere.
FIGURE 8
4 Conclusion and discussion
The effects of density structures on MS waves are important to understand the distribution and propagation of MS waves in the Earth’s magnetosphere, which are attracting more and more attention. However, previous simulations and theoretical models [,] only include the reflection of MS waves caused by the density structure. To include both the absorption and reflection of waves, we utilize a self-consistent model, i.e., PIC model, to study the propagation of MS waves across density structures. We find that both perpendicular and quasi-perpendicular propagating MS waves can be effectively reflected and absorbed by the fine-scale density structure. Generally, the absorption of MS waves is as important as the reflection when MS waves propagate through the density structure, and they are strongly dependent on the shape of the density structure. The reflection of MS waves is positively correlated with the height but is inversely related to the width of a density structure. While the absorption of MS waves is positively correlated with the height and width of a density structure. Our simulation results reveal that the absorption also plays an important role in the propagation of MS waves in the Earth’s magnetosphere, which can help better understand the properties and distribution of MS waves.
To obtain the reliable reflectivity, absorptivity, and transmissivity of the MS waves, we must ensure that the total energy of this system is conserved. As shown in Figures 3A, 7A, it is clear shown that the total energy is well conserved within a margin of error below 0.1%, which is much lower than the energy change (>5%) of charged particles or wave fields. This is true for all simulation runs in this study. Thus, the dependences of reflectivity and absorptivity of MS waves on the shape of density structure as shown in Figures 5, 8 are quite reliable. Since the reflectivity of MS wave is strongly dependent on the density gradient, so the will increase with the increase of the height or the decrease of width, i.e., steep density structure. While, the absorptivity should be positively correlated with the number of particles inside the density structure, so the increases with the increase of the height or width, i.e., large density structure. However, because the corresponding transmissivity relies on the sum of and , the dependence of on the shape of density structure is somehow unpredictable.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary materials, further inquiries can be directed to the corresponding author/s. Figure data are available at https://cstr.cn/14804.41.sciencedb.space.00813.00D00D7D.
Author contributions
TS: Writing–original draft. XG: Writing–review and editing. YK: Writing–review and editing. QL: Writing–original draft, Writing–review and editing, Supervision. XW: Writing–original draft, Writing–review and editing, Validation.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article. This work was supported by the NSFC grant 41774151, 42104155, B-type Strategic Priority Program of the Chinese Academy of Sciences, Grant No. XDB41000000.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2024.1254024/full#supplementary-material
References
1.
SantolíkONěmecFGereováKMacúšováEde ConchyYCornilleau-WehrlinN. Systematic analysis of equatorial noise below the lower hybrid frequency. Ann Geophys (2004) 22:2587–95. 10.5194/angeo-22-2587-2004
2.
RussellCHolzerRSmithE. OGO 3 observations of ELF noise in the magnetosphere: 1. Spatial extent and frequency of occurrence. J Geophys Res (1969) 74:755–77. 10.1029/JA074i003p00755
3.
RussellCTHolzerRESmithEJ. OGO 3 observations of ELF noise in the magnetosphere 2. The nature of the equatorial noise. J Geophys Res (1970) 75:755–68. 10.1029/JA075i004p00755
4.
GurnettDA. Plasma wave interactions with energetic ions near the magnetic equator. J Geophys Res (1976) 81:2765–70. 10.1029/JA081i016p02765
5.
XiaoFZhouQHeYYangCLiuSBakerDNet alPenetration of magnetosonic waves into the plasmasphere observed by the Van Allen Probes. Geophys Res Lett (2015) 42:7287–94. 10.1002/2015GL065745
6.
MaQLiWThorneRMAngelopoulosV. Global distribution of equatorial magnetosonic waves observed by THEMIS. Geophys Res Lett (2013) 40:1895–901. 10.1002/grl.50434
7.
BoardsenSAHospodarskyGBMinKAverkampTFBoundsSRKletzingCAet alDetermining the wave vector direction of equatorial fast magnetosonic waves. Geophys Res Lett (2018) 45(16):7951–9. 10.1029/2018GL078695
8.
KasaharaYKenmochiHKimuraI. Propagation characteristics of the ELF emissions observed by the satellite Akebono in the magnetic equatorial region. Radio Sci (1994) 29(4):751–67. 10.1029/94RS00445
9.
KeYChenLGaoXLuQWangXChenRet alWhistler-mode waves trapped by density irregularities in the earth's magnetosphere. Geophys Res Lett (2021) 48:e2020GL092305. 10.1029/2020GL092305
10.
MeredithNHorneRAndersonR. Survey of magnetosonic waves and proton ring distributions in the Earth’s inner magnetosphere. J Geophys Res (2008) 113:A06213. 10.1029/2007JA012975
11.
MinKNěmecFLiuKDentonREBoardsenSA. Equatorial propagation of the magnetosonic mode across the plasmapause: 2-D PIC simulations. J Geophys Res Space Phys (2019) 124:4424–44. 10.1029/2019JA026567
12.
BortnikJThorneRM. Transit time scattering of energetic electrons due to equatorially confined magnetosonic waves. J Geophys Res (2010) 115:A07213. 10.1029/2010JA015283
13.
BortnikJThorneRMNiBLiJ. Analytical approximation of transit time scattering due to magnetosonic waves. Geophys Res Lett (2015) 42:1318–25. 10.1002/2014GL062710
14.
HorneRBThorneRMGlauertSAMeredithNPPokhotelovDSantolıkO. Electron acceleration in the Van Allen radiation belts by fast magnetosonic waves. Geophys Res Lett (2007) 34(17):L17107. 10.1029/2007GL030267
15.
HorneRBWheelerGVAlleyneHSCK. Proton and electron heating by radially propagating fast magnetosonic waves. J Geophys Res (2000) 105(27):27597–610. 10.1029/2000JA000018
16.
LiXTaoX. Validation and analysis of bounce resonance diffusion coefficients. J Geophys Res Space Phys (2018) 123:104–13. 10.1002/2017JA024506
17.
ChenLSunJLuQWangXGaoXWangDet alTwo-dimensional particle-in-cell simulation of magnetosonic wave excitation in a dipole magnetic field. Geophys Res Lett (2018) 45:8712–20. 10.1029/2018GL079067
18.
LiJBortnikJThorneRMLiWMaQBakerDNet alUltrarelativistic electron butterfly distributions created by parallel acceleration due to magnetosonic waves. J Geophys Res Space Phys (2016) 121:3212–22. 10.1002/2016JA022370
19.
XiaoFLYangCSuZPZhouQHHeZGHeYHet alWave-driven butterfly distribution of Van Allen belt relativistic electrons. Nat Commun (2015) 6(1):8590. 10.1038/ncomms9590
20.
NiBBZouZYFuSCaoXGuXDXiangZ. Resonant scattering of radiation belt electrons by off-equatorial magnetosonic waves. Geophys Res Lett (2018) 45(3):1228–36. 10.1002/2017GL075788
21.
YangCSuZXiaoFZhengHWangYWangSet alA positive correlation between energetic electron butterfly distributions and magnetosonic waves in the radiation belt slot region. Geophys Res Lett (2017) 44:3980–90. 10.1002/2017GL073116
22.
MaQLLiWThorneRMBortnikJKletzingCAKurthWSet alElectron scattering by magnetosonic waves in the inner magnetosphere. J Geophys Res Space Phys (2016) 121(1):274–85. 10.1002/2015JA021992
23.
SunJGaoXLuQChenLLiuXWangXet alSpectral properties and associated plasma energization by magnetosonic waves in the Earth’s magnetosphere: particle-in-cell simulations. J Geophys Res Space Phys (2017) 122:5377–90. 10.1002/2017JA024027
24.
YuanZYuXHuangSQiaoZYaoFFunstenHO. Cold ion heating by magnetosonic waves in a density cavity of the plasmasphere. J Geophys Res Space Phys (2018) 123:1242–50. 10.1002/2017JA024919
25.
ZhouQJiangZYangCHeYLiuSXiaoF. Correlated observation on global distributions of magnetosonic waves and proton rings in the radiation belts. J Geophys Res Space Phys (2021) 126:e2020JA028354. 10.1029/2020JA028354
26.
ChenLThorneRM. Perpendicular propagation of magnetosonic waves. Geophys Res Lett (2012) 39:L14102. 10.1029/2012GL052485
27.
XiaoFZhouQHeZTangL. Three-dimensional ray tracing of fast magnetosonic waves. J Geophys Res (2012) 117:A06208. 10.1029/2012JA017589
28.
SantolíkOParrotMNěmecF. Propagation of equatorial noise to low altitudes: decoupling from the magnetosonic mode. Geophys Res Lett (2016) 43:6694–704. 10.1002/2016GL069582
29.
LiuXChenLYangLXiaZMalaspinaDM. One-dimensional full wave simulation of equatorial magnetosonic wave propagation in an inhomogeneous magnetosphere. J Geophys Res Space Phys (2018) 123:587–99. 10.1002/2017JA024336
30.
MaukBHFoxNJKanekalSGKesselRLSibeckDGUkhorskiyA. Science objectives and rationale for the radiation belt storm probes mission. Space Sci Rev (2013) 179(1):3–27. 10.1007/s11214-012-9908-y
31.
MaQLiWChenLThorneRMAngelopoulosV. Magnetosonic wave excitation by ion ring distributions in the Earth’s inner magnetosphere. J Geophys Res Space Phys (2014) 119:844–52. 10.1002/2013JA019591
32.
BoardsenSAHospodarskyGBKletzingCAPfaffffRFKurthWSWygantJRet alVan Allen Probe observations of periodic rising frequencies of the fast magnetosonic mode. Geophys Res Lett (2014) 41:8161–8. 10.1002/2014GL062020
33.
SunJLuQWangXLiuXGaoXYangH. Modulation of magnetosonic waves by background plasma density in a dipole magnetic field: 2-D PIC simulation. J Geophys Res Space Phys (2021) 126:e2021JA029729. 10.1029/2021JA029729
34.
Stix. Waves in plasmas. College Park, MD: American Institute of Physics (1992).
35.
StreltsovAVGoyalR. Whistlers in micro ducts. J Geophys Research:SpacePhysics (2021) 126:e2021JA029868. 10.1029/2021JA029868
36.
YuXYuanZOuyangZYaoF. Effects of the plasmapause on the radial propagation of fast magnetosonic waves: an analytical approach. J Geophys Res Space Phys (2021) 126:e2020JA028330. 10.1029/2020JA028330
Summary
Keywords
magnetosonic wave, wave propagation, density structure, magnetosphere, paticlein-cell simulation, wave-particle interaction
Citation
Shao T, Gao X, Ke Y, Lu Q and Wang X (2024) The effects of plasma density structure on the propagation of magnetosonic waves: 1-D particle-in-cell simulations. Front. Phys. 12:1254024. doi: 10.3389/fphy.2024.1254024
Received
06 July 2023
Accepted
15 April 2024
Published
28 May 2024
Volume
12 - 2024
Edited by
Chuanfei Dong, Boston University, United States
Reviewed by
Mehran Shahmansouri, Arak University, Iran
Xu Liu, The University of Texas at Dallas, United States
Hongyang Zhou, University of Helsinki, Finland
Updates
Copyright
© 2024 Shao, Gao, Ke, Lu and Wang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xinliang Gao, gaoxl@ustc.edu.cn; Yangguang Ke, keyg@ustc.edu.cn
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.