Characteristics of plasma boundaries with large density gradients and their effects on Kelvin–Helmholtz instability

Abstract

Boundaries between space plasmas occur in numerous contexts and scales, from astrophysical jets to planetary magnetospheres. Mass and momentum transport across boundaries poses a fundamental problem in magnetospheric physics. Kelvin–Helmholtz instability (KHI) is a promising mechanism to facilitate transport. Although previous studies have suggested KHI occurrence in various space plasmas, theory predicts that compressibility prevents KHI excitation at boundaries with large density gradients because of previously considered boundary structures where density varies with velocity. Based on the observations of a large density gradient boundary by MAVEN at Mars, where we can observe an extreme case, in this study, we show that it is the entropy, instead of the previously considered density, that varies with the velocity in the real velocity-sheared boundary. The entropy-based boundary structure places the velocity shear in a lower-density region than the traditional density-based structure and weakens the compressibility effect. This new boundary structure thus enables KHI excitation even at large density gradient boundaries, such as at the ionopause of unmagnetized planets and the plasmapause of magnetized planets. The result suggests the ubiquitous occurrence of KHI in the plasma universe and emphasizes its important role in planetary cold plasma escape from unmagnetized planets.

1 Introduction

Momentum and mass transport mechanisms at a velocity-sheared boundary in the space plasma are important processes that affect the supply of magnetospheric plasma at various planets. In various space plasmas, collisions are often negligible (Matthaeus et al., 2003), and thus, momentum and mass transport primarily occurs through couplings between constituent charged particles and the electric/magnetic fields (Seki et al., 2003). When two different plasma regimes come into contact, instabilities facilitate mass and momentum transport across the ambient magnetic field. Such boundaries occur in numerous contexts and scales from astrophysical plasmas such as at neutron star mergers (Price and Rosswog, 2006), star-forming molecular clouds (Berne et al., 2010), nova explosions (Casanova et al., 2011), and solar corona (Hillier and Arregui, 2019), in addition to various planetary magnetospheres (Hasegawa et al., 2004; Slavin et al., 2008; Delamere et al., 2013). The Kelvin–Helmholtz instability (KHI) (Fujimoto and Terasawa, 1994) is considered a promising candidate to facilitate momentum and mass transfers across boundaries.

One of the important parameters for assessment of the occurrence of KHI is compressibility. The stability analysis of compressional plasma predicts that KHI excitation is suppressed when the density ratio across the boundary () becomes larger, due to the stabilizing effects of the compressibility with a high Mach number (, where represents the magnetosonic Mach number) (Miura and Pritchett, 1982). An experimental study indicates a reduction in the growth rate of KHI and supports the theoretical prediction (Wan et al., 2015). In theoretical considerations, both the velocity shear and density gradient are often expressed with the hyperbolic tangent (Miura and Pritchett, 1982; Amerstorfer et al., 2010; Matsumoto and Seki, 2010). This traditional density-based model of the velocity-sheared boundary has a high region in the middle of the boundary, and the KHI suppression by compressibility can become significant in large density gradient cases.

In Earth, which has a strong intrinsic magnetic field, KHI has been identified by in situ satellite observations at the boundary between the shocked solar wind and terrestrial magnetosphere (Hasegawa et al., 2004), where the density gradient is usually small (with the density ratio between the two sides of the boundary being ) and compressibility is not very important. The in-plane magnetic field (magnetic field parallel to velocity) also stabilizes KHI in the low regime, and the magnetic field geometry effects are pronounced enough to excite KHI in the Earth’s magnetosphere (Kavosi et al., 2023). Thus, the KHI excitation at the magnetopause depends on the dipole tilt angle and interplanetary magnetic field (IMF) (Hwang et al., 2023). KHI can also be unstable in the inner magnetosphere, associated with localized flow channels (Chaston, 2022).

In the case of unmagnetized planets such as Mars and Venus, that do not have a global intrinsic magnetic field, the velocity-sheared boundary appears between the shocked solar wind and the ionized upper atmosphere (ionosphere), i.e., around the ionopause. The magnetic field geometry effects are also important at Venus and Mars, and KHI is mainly excited in the perpendicular plane to the IMF (Dang et al., 2022). Due to the large density in the ionosphere, the compressibility effects are expected to be important and can suppress the occurrence of KHI. Despite this expectation from theory, KHI vortex-like signatures (boundary waves) have been observed on Mars (Ruhunusiri et al., 2016; Poh et al., 2021; Wang et al., 2022) and Venus (Chong et al., 2018), where is larger than in the terrestrial case.

Particularly on Mars, tends to be large due to the low solar wind density, compared to Venus, at the Martian orbit (), and KHI is expected to be stabilized by compressibility (Miura and Pritchett, 1982). Thus, the observation of KHI-like signatures around the Martian ionopause poses a fundamental question about the ubiquity of KHI in space plasmas, where a large density gradient across velocity-sheared boundaries is sometimes expected: What enables the KHI excitation even at a plasma boundary with a large density gradient? Understanding of the characteristics of KHI around the ionopause is also important for comprehending the ionospheric cold ion outflows from the Martian atmosphere (Inui et al., 2019)

Using a comprehensive in situ plasma observation at Mars by MAVEN and MHD simulations of KHI, in this study, we investigate the characteristics of plasma boundaries with large density gradients and their effects on Kelvin–Helmholtz instability. Following the description of the instrumentation and model in Section 2, the results of the comparison between the observations and simulations are shown in Section 3. Based on observations, we propose a new realistic boundary structure for the velocity-sheared layer with a large density gradient based on entropy and conduct MHD simulations of KHI. In Section 4 and Section 5, discussions and conclusion are provided, respectively.

2 Instrumentation and model description

A comprehensive in situ observation of plasma at Mars by MAVEN (Andersson et al., 2015; Connerney et al., 2015; Halekas et al., 2015; McFadden et al., 2015) provides an ideal opportunity to investigate the fundamental question described in Section 1. We used the data obtained in the MAVEN orbit #118 on 20 October 2014 from the MAG instrument (Connerney et al., 2015) for the vector magnetic fields, STATIC (McFadden et al., 2015)/SWIA (Halekas et al., 2015) for ion moments with/without mass discrimination, and LPW (Andersson et al., 2015) for electron density in the ionosphere. We used SWIA density and velocity data only for the time period when the SupraThermal and Thermal Ion Composition (STATIC) energy coverage was limited due to its observation mode, as described in Section 3.1.

A high-resolution, numerically robust magnetohydrodynamic (MHD) simulation code CANS+ (Matsumoto et al., 2019) is used to simulate KHI in velocity-sheared boundary layers with large density gradients. The nonlinear evolution of the physical quantities—plasma number density, pressure, velocity, and magnetic field B—is followed on the basis of MHD equations. We numerically solve the equations by using the fifth-order MP5 scheme (Suresh and Huynh, 1997) and the approximate Rieman solver of the HLLD scheme (Miyoshi and Kusano, 2005). In order to compare traditional density-based and new entropy-based initial structures of the velocity-sheared boundary, we conducted two sets of two-dimensional MHD simulations. Two models of the boundary structure will be described in detail, together with observations, in Section 3.1. Based on the observations, the number density, velocity, temperature, and magnetic field strength in the shocked solar wind (ionosphere) side are set to 2 (10⁴) cm⁻³, 177 (0) km/s, 275 (0.055) eV, and 8 (8) nT, respectively. It corresponds to . The total number of grid points used is 504 600. The spatial resolution is set as . Thus, the velocity shear layer is resolved with about 48 grids, and the resolution is good enough to reproduce the linear growth of KHI (Matsumoto et al., 2019). The settings for the two simulation runs other than the initial conditions are identical.

3 Results

3.1 Boundary observations by MAVEN and the entropy-based boundary structure model

Since Mars does not possess a global magnetic field, the solar wind directly interacts with the ionosphere. As it approaches Mars, the supersonic solar wind forms a bow shock in front of the planet (Figure 1), and the decelerated, shocked solar wind interacts with the dense ionosphere. At the boundary between the shocked solar wind and the ionosphere, the density gradient is large, with a typical density ratio of . A combination of mechanical and electrostatic attenuators of the STATIC instrument onboard MAVEN (McFadden et al., 2015) enables continuous ion density measurements over four orders of magnitude, which is essential in order to reveal the detailed structure inside the boundary.

FIGURE 1

Trotignon et al., 2006) locations, respectively. The white line shows the MAVEN orbital projection onto the - coordinates. For the period from 09:45 to 11:00 UT, which corresponds to the time period shown in Figure 2, the total ion density is color-coded on the orbit. Horizontal gray thin lines show the center of the mass ion velocity component parallel to the perpendicular flow in the shocked solar wind.

MAVEN observed an example of a velocity-sheared boundary with a large density gradient of between the shocked solar wind and dayside Martian ionosphere, as shown by the color codes and gray bars in Figure 1. The magnetic field in the ionosphere (light blue bar at the bottom in Figure 2) and shocked solar wind (orange) has a similar value (−8 nT) with some enhancement inside the boundary (the top panel in Figure 2). The high magnetic field strength at low altitudes is due to the localized crustal magnetic fields and will not have major effects on the global boundary structure. Alternating high and low energy populations are seen inside the boundary (second and third panels from the top of Figure 2), and increases in heavy ions at masses of −16 (O⁺) and 32 (O₂⁺) coincide with the low energy populations (not shown). It is consistent with the density variations of each ion species shown in the fifth panel in Figure 2, i.e., O₂⁺ is dominant in the ionosphere, while proton contribution is important in the shocked solar wind side. As the spacecraft traverses from the ionosphere to the shocked solar wind, the velocity increases from −0 to −177 km/s (sixth panel in Figure 2), while the density decreases from −6570 to 1.3 cm⁻³ (5th panel), i.e., , by taking averages for each region.

FIGURE 2

Connerney et al., 2015), the ion energy–time, mass–time spectrograms, density of each species, combined electron/ion density, velocity, and region identification bar, respectively. The black line in the sixth panel of the ion density shows the total ion density obtained both from STATIC (McFadden et al., 2015) and SWIA (Halekas et al., 2015) in addition to electron density from LPW (Andersson et al., 2015) in the ionosphere (green) to confirm the accuracy of the density observation. The ion velocity panel displays the bulk velocity component in the direction of the shocked solar wind flow averaged over the orange bar time interval at the bottom. Note that we used SWIA density and velocity data for the time period shown in the gray bar at the bottom of the third panel due to the limited energy coverage of STATIC. The bottom color bar indicates time intervals used in the data analysis. Light blue (orange) corresponds to the interval used to calculate ionospheric (shocked solar wind) parameters, while pink indicates the boundary region in which the velocity shear and density gradient are observed.

In order to investigate the location of the velocity shear relative to the density gradient, the observed relation between the density and velocity is determined and plotted in Figure 3A. In theoretical considerations (Miura and Pritchett, 1982; Amerstorfer et al., 2010; Matsumoto and Seki, 2010), the velocity shear, taking the X-axis in the original shocked solar wind direction and Y across the shear, is often expressed with the hyperbolic tangent asin the rest frame of the ionosphere, where and indicate velocity in the shocked solar wind and the half-width of the shear layer, respectively. In the case of a density gradient, previous studies (Amerstorfer et al., 2010; Matsumoto and Seki, 2010) have assumed that the density gradient is expressed aswhere and show number densities in the ionospheric and shocked solar wind sides, respectively. This traditional density-based model of the velocity-sheared boundary is displayed in Figure 4A and also by the blue line in Figure 3A. The observation therefore reveals that the observed density–velocity relation (black dots in Figure 4A) at the boundary with a large density gradient differs largely from the density-based model (blue line).

FIGURE 3

FIGURE 4

Cabral and Leedom, 1993) in the simulation domain at 209, where a KHI vortex is developed in panel (D), but not in panel (C). The |V| profile with white noise was used as a texture for LIC visualization.

In order to describe the observed boundary structure, we here introduce the entropy per unit mass, (, where indicates the thermal pressure), instead of the density, to construct an entropy-based model:where and show the entropy in the ionospheric and shocked solar wind sides, respectively. Figure 3A shows that the density–velocity relation with this entropy-based model (red line) agrees much better with observations than the density-based model (blue). It should be noted that the entropy is calculated under the assumption that the pressure is constant across the boundary. We investigated other MAVEN orbits that had comprehensive plasma data, including low-energy range observations by STATIC around the periapsis, which were available between 1 December 2014 and 26 January 2015. The observations show that all of the investigated boundary structure observations agree better with the entropy-based model than the density-based model. Compared to other physical parameters, such as momentum, to construct a new model, the entropy-based model shows the best agreement with the observation. We think it is natural to consider the boundary structure in which the entropy varies with the velocity since the density perturbation is theoretically transported by the entropy mode wave in the MHD approximation. In the next section, we use the traditional density-based and new entropy-based models as the initial conditions of the KHI simulations to observe the effects of the boundary model on momentum and mass transfers across the boundary by KHI.

3.2 MHD simulation results

We here conducted two-dimensional MHD simulations of KHI with a transverse magnetic field for the two initial conditions, i.e., the density-based model (Figure 4A) and entropy-based model (Figure 4B). The density, velocity, temperature, and magnetic field strength in the shocked solar wind (ionosphere) side are set to 2 (10⁴) cm⁻³, 177 (0) km/s, 275 (0.055) eV, and 8 (8) nT based on observations shown in Section 3.1, respectively. These parameters correspond to . The settings for the two simulation runs other than the initial conditions are identical. Figure 3B shows the time evolution with the entropy-based initial condition (red line) at the early (light blue dots) and fully developed (black dots) stages of the KHI time evolution. The observed density–velocity relation (black dots in Figure 3A) resembles that of the early-stage KHI (light blue in Figure 3B), suggesting that KHI had already started to grow in the real Martian boundary.

The center upper part of Figure 4D shows a well-developed vortex with its center at approximately . The velocity of the vortex along the boundary was approximately 62 km/s. In contrast to the KHI vortex development with the entropy-based model in Figure 4D, we did not see any clear growth of KHI in the case of the density-based model (Figure 4C). In the density-based boundary structure, the magnetosonic Mach number (defined by the speed of sound and Alfvén velocity as ), in most parts of the velocity shear, becomes large enough () to suppress the KHI growth due to the compressibility (Miura and Pritchett, 1982). A remarkable difference between the density-based and entropy-based models is the relative location of the velocity shear to the density gradient. As shown in Figure 4, the entropy-based model places the velocity shear in the lower-density region, i.e., the speed of sound at the center of the velocity shear can be larger, i.e., is smaller than in the density-based model. It thus weakens the compressibility suppression effect and enables KHI excitation.

4 Discussion

The real boundary structure with velocity shear and density gradient in space plasma reported here requires alteration of the initial boundary structure in theoretical studies of KHI to estimate mass and momentum transport across the boundary. In collisionless space plasmas, plasma mixing across the magnetic field lines occurs at small-scale structures, caused by secondary instabilities excited in the rolled-up vortex of KHI, i.e., at a well-developed KHI. To understand the momentum and mass transport mechanisms across the velocity-sheared boundary, high-resolution kinetic simulations of KHI have been conducted (Matsumoto and Seki, 2010; Dargent et al., 2019). Some of the previous studies have identified the importance of including electron inertial effects when assessing mass transport mechanisms across the boundaries. The observed entropy-based boundary structure presented here has a fundamental difference from the traditional initial conditions of the various simulation studies of KHI (density-based model) (Suresh and Huynh, 1997; Casanova et al., 2011; Andersson et al., 2015; Matsumoto et al., 2019; Hwang et al., 2023). The change in the initial condition largely alters the time evolution of KHI, even in the small and moderate density gradient cases. For example, the usage of the proposed entropy-based model based on the observations will be important for understanding the interplay between KHI and the lower-hybrid drift instability in plasma mixing (Dargent et al., 2019).

Although kinetic particle-in-cell (PIC) simulations are important to investigate plasma mixing processes, MHD (Michael et al., 2021; Dang et al., 2022) or hybrid (Wang et al., 2023) simulations have the advantage of conducting three-dimensional global simulations of the interaction between the solar wind and planetary magnetosphere/ionosphere. The comparison between MHD and PIC simulations indicates that the MHD approximation can capture important features of the plasma mixing layer in certain conditions when the grid resolution and numerical dissipation are small enough (Matsumoto and Seki, 2010). Thus, this study encourages investigations of dayside velocity-sheared boundary structures in the global simulation, which determine the initial condition for KHI growth, to compare them with the entropy-based boundary structure obtained from observations in this study.

The specific entropy used to formulate the realistic boundary structure proposed in this study has been utilized to explore transport of plasma from the solar wind to various regions in the Earth’s magnetosphere (Borovsky and Cayton, 2011). is a conserved quantity in the adiabatic transport of isotropic plasma, and it is often used to track plasma transport, i.e., there is often a gap of across the boundaries. Non-conservation of is expected when nonadiabatic processes such as plasma mixing and wave–particle interaction occur. In the MHD approximation, the density perturbation is transported by the entropy mode. We thus think that it is a reasonable choice for the velocity-sheared boundary model to assume that varies with velocity in the boundary.

The entropy-based model obtained in this study has no artificial peak in the magnetosonic Mach number and momentum inside the boundary, in contrast to the density-based model, as shown by the red lines in Figures 4A,B. It prevents the overestimation of mass and momentum transport. The observation also indicates that the realistic boundary structure in space plasma can facilitate KHI over a wide range of density gradients by evading compressibility suppression, and thus enabling the ubiquitous occurrence of Kelvin–Helmholtz instability in the plasma universe. It might explain the KHI occurrence reported in various space and astrophysical plasmas (Hasegawa et al., 2004; Price and Rosswog, 2006; Slavin et al., 2008; Berne et al., 2010; Casanova et al., 2011; Delamere et al., 2013; Hillier and Arregui, 2019). Specific to Mars, KHI excitation can cause cold ions to escape from the Martian ionosphere (Inui et al., 2019), and it might provide a new escape channel to answer the mystery of the fate of its ancient carbon dioxide (Barabash et al., 2007) and understand the ancient Martian climate change (Jakosky et al., 2018).

5 Conclusion

These results suggest that the suppression of KHI can be weakened in the realistic velocity-sheared boundary structure even when the density gradient across the boundary is large, and it can facilitate the ubiquitous occurrence of Kelvin–Helmholtz instability in various planetary magnetospheres and astrophysical plasmas. It will be useful to use the entropy-based boundary model for future simulation studies of velocity-sheared boundaries with a density gradient.

References

• 1

  AmerstorferU. V.ErkaevN. V.TaubenschussU.BiernatH. K. (2010). Influence of a density increase on the evolution of the Kelvin–Helmholtz instability and vortices. Phys. Plasmas17, 072901. 10.1063/1.3453705

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  AnderssonL.ErgunR. E.DeloryG. T.ErikssonA.WestfallJ.ReedH.et al (2015). The Langmuir probe and waves (LPW) instrument for MAVEN. Space Sci. Rev.195, 173–198. 10.1007/s11214-015-0194-3

  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  BarabashS.FedorovA.LundinR.SauvaudJ. A. (2007). Martian atmospheric erosion rates. Science315, 501–503. 10.1126/science.1134358

  • CrossRef
  • Google Scholar

• 4

  BerneO.MarcelinoN.CernicharoJ. (2010). Waves on the surface of the Orion molecular cloud. Nature466, 947–949. 10.1038/nature09289

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  BorovskyJ. E.CaytonT. E. (2011). Entropy mapping of the outer electron radiation belt between the magnetotail and geosynchronous orbit. J. Geophys. Res.116, A06216. 10.1029/2011JA016470

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  CabralB.LeedomL. C. (1993). “Imaging vector fields using line integral convolution,” in Proceedings of the 20th annual conference on Computer graphics and interactive techniques (SIGGRAPH 93), 263–270.

  • Google Scholar
  • View reference in article

• 7

  CasanovaJ.JoséJ.García-BerroE.ShoreS. N.CalderA. C. (2011). Kelvin–Helmholtz instabilities as the source of inhomogeneous mixing in nova explosions. Nature478, 490–492. 10.1038/nature10520

  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  ChastonC. C. (2022). Flow channels and the generation of Alfvenic turbulence along storm-time inner magnetospheric field-lines. Geophys. Res. Lett.49, e2022GL101321. 10.1029/2022GL101321

  • CrossRef
  • Google Scholar

• 9

  ChongG. S.PopeS.WalkerS. N.FrahmR. A.ZhangT. L.FutaanaY. (2018). A statistical study of ionospheric boundary wave formation at Venus. J. Geophys. Res.123, 7668–7685. 10.1029/2018JA025644

  • CrossRef
  • Google Scholar

• 10

  ConnerneyJ. E. P.EspleyJ.LawtonP.MurphyS.OdomJ.OliversenR.et al (2015). The MAVEN magnetic field investigation. Space Sci. Rev.195, 257–291. 10.1007/s11214-015-0169-4

  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  DangT.LeiJ.ZhangB.ZhangT.YaoZ.LyonJ.et al (2022). Oxygen ion escape at Venus associated with three-dimensional Kelvin-Helmholtz instability, Geophys. Res. Lett.49, e2021GL096961. 10.1029/2021GL096961

  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  DargentJ.LavorentiF.CalifanoF.HenriP.PucciF.CerriS. S. (2019). Interplay between Kelvin Helmholtz and lower-hybrid drift instabilities. J. Plasma Phys.85, 805850601. 10.1017/S0022377819000758

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  DelamereP. A.WilsonR. J.ErikssonS.BagenalF. (2013). Magnetic signatures of Kelvin-Helmholtz vortices on Saturn's magnetopause: global survey. J. Geophys. Res.118, 393–404. 10.1029/2012JA018197

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  FujimotoM.TerasawaT. (1994). Anomalous ion mixing within an MHD scale Kelvin-Helmholtz vortex. J. Geophys. Res.99, 8601–8613. 10.1029/93ja02722

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  HalekasJ. S.TaylorE. R.DaltonG.JohnsonG.CurtisD. W.McFaddenJ. P.et al (2015). The solar wind ion analyzer for MAVEN. Space Sci. Rev.195, 125–151. 10.1007/s11214-013-0029-z

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  HasegawaH.FujimotoM.PhanT. D.RèmeH.BaloghA.DunlopM. W.et al (2004). Transport of solar wind into Earth’s magnetosphere through rolled-up Kelvin–Helmholtz vortices. Nature430, 755–758. 10.1038/nature02799

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  HillierA.ArreguiI. (2019). Coronal cooling as a result of mixing by the nonlinear kelvin–helmholtz instability. ApJ885, 101. 10.3847/1538-4357/ab4795

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  HwangK.-J.WangC.-P.NykyriK.HasegawaH.TapleyM. B.BurchJ. L.et al (2023). Kelvin-Helmholtz instability-driven magnetopause dynamics as turbulent pathway for the solar wind-magnetosphere coupling and the flank-central plasma sheet communication. Front. Astron. Space Sci.10, 1151869. 10.3389/fspas.2023.1151869

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  InuiS.SekiK.SakaiS.BrainD. A.HaraT.McFaddenJ. P.et al (2019). Statistical study of heavy ion outflows from Mars observed in the Martian‐induced magnetotail by MAVEN. J. Geophys. Res.124, 5482–5497. 10.1029/2018JA026452

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  JakoskyB. M.BrainD.ChaffinM.CurryS.DeighanJ.GrebowskyJ.et al (2018). Loss of the Martian atmosphere to space: present-day loss rates determined from MAVEN observations and integrated loss through time. Icarus315, 146–157. 10.1016/j.icarus.2018.05.030

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  KavosiS.RaederJ.JohnsonJ. R.NykyriK.FarrugiaC. J. (2023). Seasonal and diurnal variations of Kelvin-Helmholtz Instability at terrestrial magnetopause. Nat. Comm.14, 2513. 10.1038/s41467-023-37485-x

  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  MatsumotoY.AsahinaY.KudohY.KawashimaT.MatsumotoJ.TakahashiH. R.et al (2019). Magnetohydrodynamic simulation code CANS+: assessments and applications. Publ. Astron. Soc. Jpn.71. 10.1093/pasj/psz064

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  MatsumotoY.SekiK. (2010). Formation of a broad plasma turbulent layer by forward and inverse energy cascades of the Kelvin-Helmholtz instability. J. Geophys. Res.115, A10231. 10.1029/2009JA014637

  • CrossRef
  • Google Scholar

• 24

  MatthaeusW. H.QinG.BieberJ. W.ZankG. P. (2003). Nonlinear collisionless perpendicular diffusion of charged particles. Astrophys. J.590, L53–L56. 10.1086/376613

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  McFaddenJ. P.KortmannO.CurtisD.DaltonG.JohnsonG.AbiadR.et al (2015). MAVEN SupraThermal and thermal ion compostion (STATIC) instrument. Space Sci. Rev.195, 199–256. 10.1007/s11214-015-0175-6

  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  MichaelA. T.SorathiaK. A.MerkinV. G.NykyriK.BurkholderB.MaX.et al (2021). Modeling Kelvin-Helmholtz instability at the high latitude boundary layer in a global magnetosphere simulation. Geophys. Res. Lett.48, e2021GL094002. 10.1029/2021GL094002

  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  MiuraA.PritchettP. L. (1982). Nonlocal stability analysis of the MHD Kelvin-Helmholtz instability in a compressible plasma. J. Geophys. Res.87, 7431–7444. 10.1029/ja087ia09p07431

  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  MiyoshiT.KusanoK. (2005). A multi-state HLL approximate Riemann solver for ideal magnetohydrodynamics. J. Comput. Phys.208, 315–344. 10.1016/j.jcp.2005.02.017

  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  PohG.EspleyJ. R.NykyriK.FowlerC. M.MaX.XuS.et al (2021). On the growth and development of nonlinear Kelvin–Helmholtz instability at Mars: MAVEN observations. J. Geophys. Res.126, e2021JA029224. 10.1029/2021JA029224

  • CrossRef
  • Google Scholar

• 30

  PriceD. J.RosswogS. (2006). Producing ultrastrong magnetic fields in neutron star mergers. Science312, 719–722. 10.1126/science.1125201

  • CrossRef
  • Google Scholar

• 31

  RuhunusiriS.HalekasJ. S.McFaddenJ. P.ConnerneyJ. E. P.EspleyJ. R.HaradaY.et al (2016). MAVEN observations of partially developed Kelvin-Helmholtz vortices at Mars. Geophys. Res. Lett.43, 4763–4773. 10.1002/2016gl068926

  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  SekiK.HiraharaM.HoshinoM.TerasawaT.ElphicR. C.SaitoY.et al (2003). Cold ions in the hot plasma sheet of Earth’s magnetotail. Nature422, 589–592. 10.1038/nature01502

  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  SlavinJ. A.AcunaM. H.AndersonB. J.BakerD. N.BennaM.GloecklerG.et al (2008). Mercury's magnetosphere after MESSENGER's first flyby. Science321, 85–89. 10.1126/science.1159040

  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  SureshA.HuynhH. T. (1997). Accurate monotonicity-preserving schemes with Runge–Kutta time stepping. J. Comput. Phys.136, 83–99. 10.1006/jcph.1997.5745

  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  TrotignonJ. G.MazelleC.BertucciC.AcuñaM. (2006). Martian shock and magnetic pile-up boundary positions and shapes determined from the Phobos 2 and Mars Global Surveyor data sets. Planet. spa. Sci.52, 357–369. 10.1016/j.pss.2006.01.003

  • CrossRef
  • Google Scholar

• 36

  WanW. C.MalamudG.ShimonyA.Di StefanoC.TranthamM.KleinS.et al (2015). Observation of single-mode, kelvin-helmholtz instability in a supersonic flow. Phys. Rev. Lett.115, 145001. 10.1103/PhysRevLett.115.145001

  • CrossRef
  • Google Scholar

• 37

  WangL.HuangC.DuA.GeY.ChenG.YangZ.et al (2023). Kelvin–helmholtz instability at Mars: in situ observations and kinetic simulations. Astrophys. J.947, 51. 10.3847/1538-4357/acc655

  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  WangX.XuX.YeY.WangJ.WangM.ZhouZ.et al (2022). MAVEN observations of the Kelvin-Helmholtz instability developing at the ionopause of Mars. Geophys. Res. Lett.49, e2022GL098673. 10.1029/2022GL098673

  • CrossRef
  • Google Scholar
  • View reference in article