Abstract
The Earth’s magnetosphere and its bow shock, which is formed by the interaction of the supersonic solar wind with the terrestrial magnetic field, constitute a rich natural laboratory enabling in situ investigations of universal plasma processes. Under suitable interplanetary magnetic field conditions, a foreshock with intense wave activity forms upstream of the bow shock. So-called 30 s waves, named after their typical period at Earth, are the dominant wave mode in the foreshock and play an important role in modulating the shape of the shock front and affect particle reflection at the shock. These waves are also observed inside the magnetosphere and down to the Earth’s surface, but how they are transmitted through the bow shock remains unknown. By combining state-of-the-art global numerical simulations and spacecraft observations, we demonstrate that the interaction of foreshock waves with the shock generates earthward-propagating, fast-mode waves, which reach the magnetosphere. These findings give crucial insight into the interaction of waves with collisionless shocks in general and their impact on the downstream medium.
Main
Collisionless super-critical shocks are highly efficient particle accelerators observed throughout the universe. When the shock geometry is quasi-parallel, that is, when the angle θBn between the upstream magnetic field and the shock normal is below 45°, shock-reflected particles can travel far upstream and excite instabilities, forming an extended foreshock or precursor hosting intense wave activity. These waves modulate particle reflection at the shock front1, can affect cosmic ray acceleration at astrophysical shocks2,3 and cause atmospheric particle escape at non-magnetized planets4. The interaction of these upstream waves with the shock, their influence on particle acceleration and their transmission into the downstream medium have received considerable attention (see, for example, refs. 1,2,5,6,7), but our understanding of these processes remains limited.
At Earth, ultra-low-frequency waves originating in the foreshock are considered the main source of magnetospheric Pc3 waves (22−100 mHz), suggesting wave transmission8,9,10,11. Compressional Pc3 waves are routinely observed in the dayside magnetosphere, where they couple with field-line resonances12, forming a remote diagnostic of magnetospheric density through magnetoseismology13. They also modulate energetic particle precipitation into the upper atmosphere14,15. Nevertheless, after decades of intensive research, it is still unclear how foreshock waves traverse through the bow shock and the downstream magnetosheath, populated by shocked solar wind plasma (Fig. 1a).
This paper addresses the important question of the interaction of low-frequency waves with a collisionless shock and presents the missing link in foreshock wave transmission. This discovery was sparked by large-scale numerical simulations, showing the global picture of wave transmission. This allows us to identify relevant observations from the Magnetospheric Multiscale (MMS) mission16, which reveal the presence of fast magnetosonic waves in the subsolar magnetosheath, with properties consistent with a foreshock source. Our findings provide compelling evidence of the process connecting foreshock and magnetospheric waves.
Foreshock ‘30 s waves’, named after their typical period at Earth, are fast magnetosonic waves generated by cyclotron-resonant instabilities driven by shock-reflected particles in the solar wind17,18,19,20. They play an important role in modulating the shape of the shock front5, affecting particle reflection at the shock1 and contributing to quasi-parallel shock reformation7. Two main observations support the connection between foreshock 30 s waves and Pc3 fluctuations: (1) in both regions, the waves have very similar frequencies showing the same dependency on the interplanetary magnetic field (IMF) strength8,9,21, and (2) magnetospheric Pc3 wave power recedes when the IMF cone angle, measured from the Sun–Earth line, increases, causing the foreshock to shift away from the subsolar point8,11,22.
Direct transmission of the waves through the bow shock and magnetosheath was initially envisioned22 and is still to date widely invoked23,24,25. However, the lack of observational evidence for fast-mode waves of foreshock origin in the magnetosheath, despite extensive surveys (see, for example, refs. 26,27), has cast doubt on this scenario. Furthermore, early numerical works suggested that foreshock waves mode-convert into Alfvén/ion cyclotron waves upon crossing the bow shock28,29 and thus could not transmit as fast-mode waves. More indirect pathways have been explored, for example, localized variations of magnetosheath dynamic pressure that would cause magnetopause motions30, or modulated precipitation into the ionosphere31, but no consensus has been reached, and this long-standing question remains unsolved.
Here, we explore this issue using state-of-the-art global ion kinetic simulations performed with Vlasiator, a hybrid-Vlasov model designed to simulate the near-Earth plasma environment (refs. 32,33 and Methods). Vlasiator has the unique capability of providing a global view of near-Earth space while including ion kinetic processes with correct scale separation, allowing for direct comparison with spacecraft observations. It has been extensively used to study foreshock processes, providing excellent agreement with observational works19,20,34,35,36,37,38, and we now apply it to the study of foreshock wave transmission. We analyse the same run as presented in ref. 38, with upstream conditions corresponding to those encountered at Earth on 20 July 2016 between 08:00 and 12:00 universal time. At that time, near-Earth space was engulfed in a large-scale solar wind structure, a magnetic cloud39, causing a long-lasting interval of steady solar wind parameters. Because of the magnetic cloud, the solar wind parameters deviate from their typical values (solar wind density n = 12 cm−3, velocity V = (−565, 0, 0) km s−1 and IMF vector B = (12.5, −6.5, 0) nT, with a 28° cone angle). The values of the magnetosonic and Alfvén Mach numbers (Mms = 5.5 and MA = 6.4, respectively) are, however, typical for Earth’s bow shock40. We select this run because the large amplitude of the foreshock waves38 facilitates their tracking across near-Earth space.
Large fluctuations of the magnetic field strength are observed throughout the simulation domain (Fig. 1a). Foreshock compressive fluctuations become stronger when approaching the shock, in agreement with spacecraft observations41. In the magnetosphere, the magnetic field variations decay when moving inwards, consistent with the attenuation of compressive waves propagating into the magnetosphere25,38.
Figure 1b shows the power spectral density (PSD) of the magnetic field strength fluctuations obtained from a wavelet transform of the time series extracted at the positions marked by the coloured circles in Fig. 1a. The PSD shows a clear peak at ~13 s in the foreshock (black curve) and in the magnetosheath (purple curve). The period of the peak agrees well with the expected wave period for these solar wind conditions (12 s according to the Takahashi et al.9 formula), below the typical 30 s because of the high IMF strength20. The wave power increases by a factor of ~5 from the foreshock to the magnetosheath, in agreement with ref. 42. A somewhat broader peak in the same period range, from 5 to 20 s, is found in the outer magnetosphere (green curve), where the wave power is further amplified by a factor of ~4.
Figure 2 shows the corresponding time series in the foreshock and magnetosheath. There is a large wave power near the predicted foreshock wave period (Fig. 2c,h), with some variability probably due to the large IMF strength causing a more complex wave field20. These magnetic field strength fluctuations are associated with density variations (Fig. 2a,f). Both are mostly positively correlated near the foreshock wave period (Fig. 2d,i), indicative of fast-mode fluctuations and ruling out mirror-mode waves. In the foreshock, the fluctuations near 12 s are only weakly compressional (Fig. 2e), consistent with 30 s wave properties43. In contrast, the wave power near 30−40 s is associated with more compressional waves and thus a different wave mode. In the magnetosheath, the fluctuations at the foreshock wave period are strongly compressional (Figs. 2j and 1c).
We apply multi-spacecraft timing analysis20,44 to the magnetic field measurements from a triplet of virtual spacecraft around x = 8RE (Earth radius = 6371 km) (Fig. 1a, purple dot) with a spacecraft separation of ~0.05RE, and find earthwards-oriented wavevectors lying at about 10° from the Sun–Earth line. The associated θkB angle, measured between the wavevector and the ambient magnetic field, is close to 90°. The wave velocity is within 20% of the local fast magnetosonic speed.
Figure 3a shows a time–position map of the magnetic field strength along the Sun–Earth line. The foreshock waves appear as alternating bands of purple and green in the right-hand part of the plot, moving earthwards as time progresses. Similar features are seen for the By and Bz components (Fig. 3b,c, blue and orange bands). Figure 3a,b shows waves in the magnetosheath with approximately the same slope as their upstream counterparts, extending to the magnetopause, located around x = 6RE. In contrast, the Bz fluctuations lose their coherency, suggesting a different wave polarization. In the magnetosheath, the waves are polarized along the magnetic field (By being the dominant component), different from their circular polarization in the foreshock (shown, for example, in refs. 17,19).
These waves propagate earthwards at the fast magnetosonic speed in the magnetosheath plasma rest frame, as indicated by their slope in the time–position map following the dashed pink lines, corresponding to the sum of the fast magnetosonic and plasma bulk speed. In addition to these fast-mode waves, Fig. 3a–c also shows structures advected by the magnetosheath flow (almost vertical lines, following the green lines), and disturbances propagating from the magnetopause to the bow shock at the fast magnetosonic speed (coloured stripes with a positive slope), which may be due to the reflection of incoming waves at the magnetopause.
As they approach the bow shock, foreshock waves participate in shock reformation, consistent with spacecraft observations7 (see also Extended Data Fig. 1). This is illustrated in Fig. 3a by the quasi-periodic motion of the density contour which serves as a proxy of the shock position45. The fast-mode structures identified by ref. 7 in the immediate downstream are the nearly vertical stripes within ~0.3RE from the shock in Fig. 3a–c (see Extended Data Fig. 1 for virtual spacecraft time series allowing direct comparison with MMS observations in ref. 7). Their propagation speed is close to the bulk flow speed, in good agreement with ref. 7. Although these structures are prominent in the immediate shock vicinity, they quickly dissipate, while those disturbances travelling earthwards at the fast magnetosonic speed become more visible further downstream.
The foreshock waves further modulate the magnetosonic Mach number upstream of the shock (Fig. 3d) and consequently the shock compression ratio. This results in total pressure variations in the downstream (Extended Data Fig. 2), which are probably the source of the fast-mode waves traversing the magnetosheath.
To further confirm the nature of these waves, Fig. 3e shows the wave power in (ω, kx) space obtained from the two-dimensional (2D) Fourier transform of Fig. 3a between x = 6.5RE and x = 9RE. We find a large wave power in fluctuations propagating earthwards (that is, with negative kx) with the fast magnetosonic speed. This lends further support to the foreshock waves traversing the magnetosheath as fast-mode waves. Again, we note that these waves are accompanied by structures travelling at the plasma bulk speed (along the dashed green line) and sunwards-propagating fluctuations at the fast magnetosonic speed (along the pink line at positive kx). The transverse wave power (Fig. 3f) suggests that there could be upstream-propagating Alfvén waves, confined to low frequencies (Extended discussion of Vlasiator 1D simulation section).
Similar dispersion plots are obtained when calculating the 2D Fourier transform along other cuts crossing the quasi-parallel magnetosheath further down on the flank (see the right-hand side of Extended Data Fig. 3 for an example), revealing that the wave transmission is not limited to the Sun–Earth line. However, the PSD peak associated with the magnetosheath fast-mode waves in virtual spacecraft data tends to disappear when moving away from the subsolar region, suggesting that the waves would not be identified in spacecraft measurements. This is likely due to other magnetosheath waves with comparable or higher power dominating the power spectrum and concealing the magnetosheath fast-mode waves.
Based on our simulation results, we analyse observations from the MMS mission16 in the subsolar region, downstream of the quasi-parallel shock, to test our numerical predictions. From 2015 to 2020, we identified intervals with compressional magnetic field fluctuations at periods consistent with those of foreshock waves. The determination of the wave properties (Methods) further required that high time resolution burst mode data were available, reducing our data set to seven intervals. For each interval, the period of the foreshock waves is obtained from a theoretical formula9. Three of the seven comprise direct foreshock wave observations shortly before or after the magnetosheath intervals, confirming the validity of the predicted wave period. The times of the events and the associated IMF conditions are listed in Extended Data Table 1. All intervals but one are associated with an IMF cone angle around 30−40° (which, in the subsolar region, approximates well the θBn angle at the shock upstream of the spacecraft) and are thus in a geometry similar to that in our numerical analysis.
Figure 4 shows MMS1 observations during one representative interval, when the spacecraft moved from the foreshock into the subsolar magnetosheath, downstream of the quasi-parallel shock, on the inbound leg of its orbit on 14 February 2020. In the foreshock, MMS1 observed weakly compressional waves near the expected foreshock wave frequency, with properties consistent with fast-mode 30 s waves (Fig. 4c–e). Shortly thereafter, waves at similar frequencies were encountered in the magnetosheath (see Extended Data Fig. 4 for the PSD), again associated with positively correlated density and magnetic field fluctuations and this time with a stronger compressive component (Fig. 4f–h), consistent with the properties of the fast-mode waves in our numerical analysis.
For this magnetosheath interval and for all other selected intervals (Extended Data Table 1), we calculate the wavevector and frequency in the plasma rest frame using a single-spacecraft method based on magnetic field and current density measurements (ref. 46 and Methods). We restrict our analysis to spacecraft frame frequencies below 0.1 Hz, that is, to the frequency range of the waves of interest. The results of this analysis for the seven magnetosheath intervals are displayed in Fig. 5. Figure 5a–d shows that the waves propagate at oblique to large angles with respect to the ambient magnetic field (θkB > 30°), and that there are both earthwards- (negative kx) and sunwards-oriented (positive kx) wavevectors in the plasma rest frame.
Figure 5e–h shows the experimental (ω, k) values in the plasma frame, separated into different ranges of θkB. We compare them with the linear solutions for the fast (red) and Alfvén (blue) modes obtained from linear Vlasov theory47. The percentages in Fig. 5e–h indicate which fraction of the data points are within the red and blue areas, which account for the different plasma conditions in the intervals. We find that a large fraction of the data points show good agreement with the fast-mode solution, while only fewer data points are closer to the Alfvén solution. This suggests that both modes co-exist in the magnetosheath, with the fast mode being predominant. Those data points that fall within the red/blue areas in Fig. 5e−h are marked with the same colours in Fig. 5a–d, suggesting that there are both earthwards- and sunwards-propagating fast-mode waves in the magnetosheath, as reported in our simulation (Fig. 3). Furthermore, the Alfvén solutions exhibit low frequency and, considering the error bars, may in fact be advected structures with no intrinsic frequency. Finally, the wave phase speed (Extended Data Fig. 5) is generally larger than the Alfvén velocity, bringing further support to the presence of fast-mode waves in the magnetosheath during these intervals.
In summary, MMS observations show the presence of earthwards-propagating fast-mode waves in the quasi-parallel subsolar magnetosheath, at frequencies matching those of foreshock waves, in agreement with our model predictions. Our numerical and observational results, therefore, provide strong evidence that foreshock waves traverse the magnetosheath as fast magnetosonic waves, as was first inferred to explain the occurrence of magnetospheric Pc3 waves (see, for example, ref. 22). Our findings regarding the wave frequency and their compressional nature in the magnetosheath are in excellent agreement with their entry into the magnetosphere as fast-mode Pc3 waves15,23,25.
However, despite these similarities with previous works, we also find that the wave propagation through the bow shock is more complex than the typically inferred direct transmission. The earthwards-oriented wavevector of the magnetosheath fast-mode waves is not consistent with that of directly transmitted foreshock waves, indicating that new waves are generated at the shock by a process modulated by the foreshock waves. These downstream waves transmit the information of the wave period through the magnetosheath, thus providing the missing link between foreshock and magnetosphere. We propose the following scenario for the downstream wave generation (see also Fig. 6): Foreshock waves modulate the magnetosonic Mach number upstream of the shock and consequently the shock compression ratio. An increased (decreased) compression ratio creates a zone of enhanced (reduced) pressure just downstream of the shock. It is then this pressure imbalance that generates fast-mode compressive/rarefaction waves travelling through the magnetosheath, in a process similar to that described by Wu et al.48 for discontinuities interacting with a shock in magnetohydrodynamic simulations, confirmed observationally49. The discontinuities in the Wu et al.48 study are comparable with the foreshock fast-mode waves in that they also cause a change in the upstream magnetosonic Mach number.
The hybrid simulations by Thomas et al.50 show that downstream-oriented fast-mode waves are also generated when pressure pulses, propagating upstream in the plasma rest frame as in the present work, hit the shock. These earthwards-propagating fast-mode waves are accompanied by mode-converted Alfvén waves in the downstream, indicating that both wave modes co-exist in their simulations as well50. As noted by the authors, the pressure pulses in their simulation correspond to another type of foreshock waves50. Our work expands upon their study in demonstrating that the interaction of 30 s waves with the shock produces similar downstream waves, and in unravelling the processes taking place at the shock. The generation of fast-mode waves in the downstream requires that foreshock waves are compressional just upstream of the bow shock, which is typically the case at Earth (see, for example, ref. 41). Although the scenario we propose is based on a fluid description, the downstream fast-mode waves were not predicted by previous theoretical works (see, for example, ref. 51), possibly because of their indirect generation, in the downstream of the shock, or because of the linear approximation used in these works, which implies that the waves are only weakly compressional.
To facilitate comparison with earlier works, we further performed one-dimensional (1D) shock simulations with different upstream conditions (Extended discussion of Vlasiator 1D simulation section and Extended Data Figs. 6 and 7). These simulations clearly show the co-existence of both the earthwards-propagating fast-mode disturbances reported here and the mode-converted Alfvén waves from earlier studies28,29,52. We also note that the fast-mode waves are not detectable in the magnetic field component transverse to both the upstream field and flow. This probably explains why the downstream fast-mode waves were not identified in previous works that focused solely on this component28,29,52.
Our numerical simulations show that these waves are more easily detected with the model’s global view (Fig. 3) than in time series mimicking actual spacecraft measurements. The downstream fast-mode waves appear as one of the dominant wave modes only in the subsolar magnetosheath, while they are masked by other modes further on the flanks (Extended Data Fig. 3, left). Large statistical surveys have been conducted using measurements away from the subsolar point, and often focused on the dominant wave mode27,53,54. This probably explains why these waves have remained elusive for so long.
The change of polarization of the fast-mode waves, from right-handed in the foreshock to linear in the magnetosheath, is due to the downstream waves being generated at the shock rather than being directly transmitted. Our global simulation demonstrates that the interaction of foreshock waves with the shock generates an array of waves in the downstream: reformation-related structures travelling with the flow (as in ref. 7), mode-converted Alfvén waves (as in earlier simulations28,29,52) and earthwards-propagating fast-mode waves (see Fig. 6 for a summary). Only the latter are responsible for the connection between the foreshock and the magnetosphere, and the generation of Pc3 magnetospheric waves.
Our results apply to super-critical collisionless shocks in general, showcasing the complexity of shock–upstream waves interactions. The consequences of our findings extend beyond near-Earth space physics, as collisionless shocks and foreshocks are ubiquitous in Solar System and astrophysical plasmas, in improving our understanding of shock processes that can affect particle acceleration2,3.
Methods
Vlasiator simulation
The numerical part of this work employs the hybrid-Vlasov model Vlasiator, targeted at global simulations of the interaction of the solar wind with the Earth’s magnetosphere32,33. In the hybrid-Vlasov formalism, ions are treated as velocity distribution functions, whose evolution is dictated by Vlasov’s equation, while electrons are considered as a massless charge-neutralizing fluid. Vlasov’s equation is coupled with Maxwell’s equations and Ohm’s law, with the Hall term. Vlasiator provides a self-consistent description of ion kinetic processes, such as the ion beam instabilities generating the foreshock waves of interest to the present study, in their global context.
The run presented here is performed in a 2D–3V space, that is, two dimensional in real space and three dimensional in velocity space. Vlasiator runs are computationally demanding, requiring millions of central processing unit hours and generating tens of terabytes of data even in a 2D setup. The simulation domain covers the equatorial plane of near-Earth space (x−y plane in the Geocentric Solar Ecliptic frame), extending from about −8RE to 76RE along x and from −60RE to 31RE along y. The simulation was run for 598 s. The solar wind is injected into the simulation from the +x boundary and can escape through the other edges of the simulation domain, at which Neumann conditions are applied. The circular inner boundary surrounding the Earth is located at 25,000 km from the Earth’s centre and is assumed to be a perfect conductor. The resolution of the simulation domain is 30 km s−1 in velocity space and 260 km in real space, the latter corresponding to about four times the ion inertial length in the solar wind. Previous works have shown and discussed in more detail that ion kinetic effects arise in Vlasiator even when not resolving the ion inertial length and lead to realistic foreshock dynamics19,55.
The solar wind is injected at the +x boundary as a Maxwellian distribution with density n = 12 cm−3, bulk velocity V = (−565, 0, 0) km s−1 and temperature T = 0.5 MK. The IMF vector B = (12.5, −6.5, 0) nT makes a 28° angle with the Sun–Earth line. The corresponding magnetosonic and Alfvén Mach number are Mms = 5.5 and MA = 6.4, respectively. All vector quantities are given in Geocentric Solar Ecliptic coordinates.
The Earth’s magnetic dipole is implemented with a realistic magnetic moment of 8.0 × 1022 A m2 and no tilt. The Earth’s dipolar field is therefore out of plane in this equatorial 2D run. Outside of the magnetosphere, in the magnetosheath and upstream of the shock, the magnetic field is dominated by the in-plane IMF. Out-of-plane components are self-consistently generated by the interaction of the IMF with the shock and the magnetosphere and by plasma instabilities.
Wavelet analysis
To determine the properties of the magnetic field and density fluctuations in both our numerical simulations and observations, we apply a Morlet wavelet transform on the time series56. Wavelet analysis allows the distribution of power in time and frequency space, revealing the temporal evolution of wave activity. The magnetic compressibility (Figs. 2 and 4) is calculated as the power of the magnetic fluctuations parallel to the mean magnetic field divided by the total power of the magnetic field fluctuations. The wavelet cross-correlation computes the common power between two time series, here the electron density and the magnetic field strength, in time–frequency space (see, for example, ref. 57).
Wave dispersion relation from MMS observations
To determine the wavevector from single-spacecraft observations, Bellan46 developed a method that uses the measured magnetic field and the plasma current derived from the density and the ion and electron velocity measurements. The wavevector, as a function of the wave frequency ω, is then given by
where J(ω) and B(ω) refer to the Fourier transforms of the current density and the magnetic field at a spacecraft frame frequency ω. This method assumes quasi-neutrality and that each spacecraft frame frequency maps to a single wavevector. Using the k-filtering method, which can resolve multiple wavevectors at a given spacecraft frame frequency, Gershman et al.58 demonstrated that this approach was justified and that both methods agreed well. The unique payload on the MMS spacecraft16 allows the simultaneous measurement of the magnetic field from the fluxgate magnetometers59 and the plasma current density J. The high time resolution capabilities from the Fast Plasma Investigation60 allow the plasma current density to be measured directly as
where ne is the electron density, q the elementary charge, and vi and ve the ion and electron bulk velocities, respectively. Here, we focus on intervals when the MMS spacecraft are operating in burst mode. The magnetic field, electron and ion distributions are sampled at 128, 33 and 6.7 Hz, respectively. Burst mode data are needed to have ion and electron measurements at a cadence comparable to that of the magnetic field measurements. Although the waves of interest have low frequencies, high-cadence measurements allow us to better reconstruct the wave dispersion relation, in providing a broad range of frequencies. Before using equation (1), all measured quantities are resampled onto the electron plasma time tags. We use Bellan’s method as the baseline sizes of MMS are too small when compared with the waves studied. Bellan’s method is applied to MMS1–3 as some of the heads of the electron spectrometer on MMS4 have failed. The wavevectors from the three MMS observatories are averaged, and we retain only those where the differences between the three individual wavevectors at a given Fourier mode are less than 35°. Using the obtained wavevectors k, the ion bulk velocity vi and the spacecraft frame frequency ω, the fluctuations are Doppler shifted to the plasma frame thus
Some Doppler shifts result in negative frequencies. This can be interpreted by considering the phase velocity vph = ωplak/k2 (see, for example, ref. 61). A negative frequency results in the direction of propagation of the wave reversing. To correct this, we reverse the sign of ωpla and k.
After Doppler shifting, we plot the relation of ωpla versus k and compare it with linear solutions of the Maxwell–Vlasov set of equations obtained from the New Hampshire Dispersion Relation Solver (NHDS)47. The NHDS solves the full hot-plasma dispersion relation for a plasma consisting of bi-Maxwellian background species for ions and electrons. For the NHDS solutions, we use the averaged observed plasma parameters as input. We calculate four dispersion relations for each branch in the given angle range, which correspond to the extremes of the ion and electron plasma β (that is, the dimensionless ratio of thermal to magnetic pressure) in our data set. The shaded areas in Fig. 5e–h are drawn between the highest and lowest values of these four solutions. To reduce the effects of magnetic nulls (which cause extremely large values of β), we calculate the median of β for each interval and use the limits βp = [5, 20], and βe = [1, 3], and isotropic ion and electron temperatures to calculate the dispersion relations.
For parallel propagation, the Alfvén branch describes the ion cyclotron wave and is heavily damped under these plasma conditions. The parallel fast-mode branch is not as heavily damped and can exist to larger wavenumbers. For quasi-perpendicular propagation, the Alfvén branch describes the kinetic Alfvén wave, which has a very low phase speed and can exist to larger wavenumbers. The fast branch transitions to the ion Bernstein wave at harmonics of the ion cyclotron frequency.
As mentioned above, the Bellan46 method assumes that there is only one k vector associated to each Fourier mode, which is probably not the case here (Fig. 3d). As a result, this method may alternately pick up wavevectors from the different co-existing modes at the same Fourier mode. The k-filtering method does not have this limitation, but cannot be used here because the spacecraft separation is too small compared with the wavelength, resulting in large errors in the determination of the wave properties62.
We compared the currents obtained from the plasma measurements with those from the curlometer method63 and found them to be in good agreement. We also calculated the wave properties using the Bellan method applied to the current estimates from the curlometer, which confirmed the presence of earthwards-propagating fast-mode waves in the magnetosheath. The wave phase speeds derived from the curlometer currents are shown in Extended Data Fig. 5.
Extended discussion of Vlasiator 1D simulation
In addition to the global simulation presented in the manuscript, we also carried out local 1D shock simulations to investigate the wave transmission in a set-up more similar to that used in previous numerical works28,29,52. The single spatial dimension is along the shock normal, while the velocity space remains three dimensional. The upstream conditions are set to n = 1 cm−3, V = (−750, 0, 0) km s−1 and T = 0.5 MK, with the fast solar wind flow ensuring a fast initialization. The IMF makes a 30° cone angle with the shock normal, similar to the subsolar shock in the global run, and the IMF strength is B = 5 nT in the run with Alfvén Mach number MA = 6.9 or B = 8.6 nT in the run with MA = 4. The shock is initialized at x = 0 in the de Hoffmann–Teller frame, and the downstream parameters are calculated using the Rankine–Hugoniot jump conditions, similarly to the shock simulations performed in ref. 55.
The top panels in Extended Data Figs. 6 and 7 show the density, magnetic field strength and magnetic field components in the simulation domain at t = 500 s from the beginning of these local runs. The shock position has moved towards negative x because the Rankine–Hugoniot jump conditions are based on magnetohydrodynamics and thus the shock is not in equilibrium when kinetic effects are at play. The upstream is filled with typical foreshock fast-mode waves that are carried towards the shock by the solar wind flow. The Bx component remains constant across the simulation domain to satisfy ∇⋅B = 0.
The dispersion plots in the bottom panels of Extended Data Figs. 6 and 7 are obtained by performing a 2D Fourier transform between x = −10RE and x = −5RE (black rectangles in panel b), and between t = 200 s and t = 500 s. As in the global simulation, a clear signal is observed in the magnetic field strength, consistent with downstream-oriented fast-mode waves (Extended Data Figs. 6d and 7d). These are accompanied by structures propagating with the bulk flow. Similar features are observed in the two local runs.
Regarding the transverse waves, appearing here on the By magnetic field component because the selected IMF is contained in the x−z plane, distinct signals for upstream- and downstream-oriented wavevectors are observed in the run with MA = 4. This is consistent with the findings by Krauss-Varban et al.28,29,52, who identified the upstream-propagating waves as foreshock waves having been mode-converted to Alfvén waves through the shock. At higher Mach number, however, the distinction between the two oppositely propagating modes tends to disappear. This is consistent with the work by Quest64, who showed that the downstream wave velocity of the Alfvén waves becomes much smaller than the Alfvén speed as the Mach number increases, resulting in the waves’ being nearly non-convective. We also note that the mode-converted Alfvén waves are restricted to lower frequencies than the other modes, probably due to wave damping29, which makes it more difficult to identify them.
This suggests that foreshock waves interacting with the bow shock give rise to both mode-converted Alfvén waves and downstream-oriented fast-mode waves. The latter were not reported in the studies by Krauss-Varban et al.28,29,52, possibly because they focused on the transverse magnetic field component, on which these waves are not detectable (Extended Data Figs. 6e and 7e), whereas they are clearly seen in the total magnetic field strength (Extended Data Figs. 6d and 7d).
On the other hand, although the local 1D run clearly demonstrates the presence of upstream-oriented Alfvén waves in the downstream, these waves are challenging to observe in the global simulation. We do, however, observe a signal consistent with those waves when performing a Fourier transform of the transverse, Bz, component (Fig. 3e, right), at low frequencies (\(\omega \le 0.05{{{\varOmega }}}_{{{{\rm{ci}}}}}^{-1}\)). This signal being much weaker in the global simulation than in the local runs could be due to the different properties of the downstream plasma: The subsolar magnetosheath flow is sub-Alfvénic in the global run, whereas the downstream remains super-Alfvénic in the local simulations. This could strongly affect the wave growth, as the sunwards-propagating Alfvén waves are effectively moving back towards the shock in the global simulation.
Data availability
The Vlasiator global run described here takes several terabytes of disk space and is kept in storage maintained within the CSC-IT Center for Science. It can be accessed through the following link: https://a3s.fi/swift/v1/AUTH_81f1cd490d494224880ea77e4f98490d/vlasiator-2d-afc. The data from the local 1D shock runs can be accessed through the following link: https://datacloud.helsinki.fi/index.php/s/NBFEj7TJf6oQjqN. Vlasiator uses a data structure developed in-house (https://github.com/fmihpc/vlsv/65), which can be read using the Analysator software https://github.com/fmihpc/analysator/66,67. Usage of Vlasiator data must comply with the data policy as described on the Vlasiator website (https://www.helsinki.fi/en/researchgroups/vlasiator/rules-of-the-road).
The MMS data are publicly available and can be found on https://lasp.colorado.edu/mms/sdc/public/. Source data are provided with this paper.
Code availability
Vlasiator http://www.helsinki.fi/en/researchgroups/vlasiator/68,69 is distributed under the GPL-2 open source license at http://github.com/fmihpc/vlasiator/70.
References
Wu, M. et al. The role of large amplitude upstream low-frequency waves in the generation of superthermal ions at a quasi-parallel collisionless shock: Cluster observations. Astrophys. J. 808, 2 (2015).
Haggerty, C. C. & Caprioli, D. Kinetic simulations of cosmic-ray-modified shocks. I. Hydrodynamics. Astrophys. J. 905, 1 (2020).
Caprioli, D., Haggerty, C. C. & Blasi, P. Kinetic simulations of cosmic-ray-modified shocks. II. Particle spectra. Astrophys. J. 905, 2 (2020).
Collinson, G. et al. Solar wind induced waves in the skies of Mars: ionospheric compression, energization, and escape resulting from the impact of ultralow frequency magnetosonic waves generated upstream of the Martian bow shock. J. Geophys. Res. Space Phys. 123, 7241–7256 (2018).
Burgess, D. Foreshock–shock interaction at collisionless quasi-parallel shocks. Adv. Space Res. 15, 159–169 (1995).
Vladimirov, A., Ellison, D. C. & Bykov, A. Nonlinear diffusive shock acceleration with magnetic field amplification. Astrophys. J. 652, 1246–1258 (2006).
Liu, T. Z., Hao, Y., Wilson, L. B., Turner, D. L. & Zhang, H. Magnetospheric Multiscale observations of Earth’s oblique bow shock reformation by foreshock ultralow-frequency waves. Geophys. Res. Lett. 48, e91184 (2021).
Troitskaya, V. A., Plyasova-Bakounina, T. A. & Gul’Elmi, A. V. The connection of Pc2–4 pulsations with the interplanetary magnetic field. Dokl. Akademiia Nauk SSSR 197, 1312–1314 (1971).
Takahashi, K., McPherron, R. L. & Terasawa, T. Dependence of the spectrum of Pc 3–4 pulsations on the interplanetary magnetic field. J. Geophys. Res. 89, 2770–2780 (1984).
Heilig, B., Lühr, H. & Rother, M. Comprehensive study of ULF upstream waves observed in the topside ionosphere by CHAMP and on the ground. Ann. Geophys. 25, 737–754 (2007).
Bier, E. A. et al. Investigating the IMF cone angle control of Pc3–4 pulsations observed on the ground. J. Geophys. Res. Space Phys. 119, 1797–1813 (2014).
Howard, T. A. & Menk, F. W. Ground observations of high-latitude Pc3–4 ULF waves. J. Geophys. Res. Space Phys. 110, A04205 (2005).
Del Corpo, A. et al. An empirical model for the dayside magnetospheric plasma mass density derived from EMMA magnetometer network observations. J. Geophys. Res. Space Phys. 125, e27381 (2020).
Rasinkangas, R. & Mursula, K. Modulation of magnetospheric EMIC waves by Pc 3 pulsations of upstream origin. Geophys. Res. Lett. 25, 869–872 (1998).
Motoba, T. et al. On the driver of daytime Pc3 auroral pulsations. Geophys. Res. Lett. 46, 553–561 (2019).
Burch, J. L., Moore, T. E., Torbert, R. B. & Giles, B. L. Magnetospheric Multiscale overview and science objectives. Space Sci. Rev. 199, 5–21 (2016).
Eastwood, J. P. et al. The foreshock. Space Sci. Rev. 118, 41–94 (2005).
Wilson, L. B. Low frequency waves at and upstream of collisionless shocks. Geophys. Monogr. Ser. 216, 269–291 (2016).
Turc, L. et al. Foreshock properties at typical and enhanced interplanetary magnetic field strengths: results from hybrid-Vlasov simulations. J. Geophys. Res. Space Phys. 123, 5476–5493 (2018).
Turc, L. et al. First observations of the disruption of the Earth’s foreshock wave field during magnetic clouds. Geophys. Res. Lett. 46, 12,644–12,653 (2019).
Le, G. & Russell, C. T. Solar wind control of upstream wave frequency. J. Geophys. Res. 101, 2571–2576 (1996).
Russell, C. T., Luhmann, J. G., Odera, T. J. & Stuart, W. F. The rate of occurrence of dayside Pc 3,4 pulsations: the L-value dependence of the IMF cone angle effect. Geophys. Res. Lett. 10, 663–666 (1983).
Clausen, L. B. N. et al. First simultaneous measurements of waves generated at the bow shock in the solar wind, the magnetosphere and on the ground. Ann. Geophys. 27, 357–371 (2009).
Francia, P., Regi, M., De Lauretis, M., Villante, U. & Pilipenko, V. A. A case study of upstream wave transmission to the ground at polar and low latitudes. J. Geophys. Res. Space Phys. 117, A01210 (2012).
Takahashi, K. et al. Propagation of ULF waves from the upstream region to the midnight sector of the inner magnetosphere. J. Geophys. Res. Space Phys. 121, 8428–8447 (2016).
Czaykowska, A., Bauer, T. M., Treumann, R. A. & Baumjohann, W. Magnetic field fluctuations across the Earth’s bow shock. Ann. Geophys. 19, 275–287 (2001).
Narita, Y. et al. Low-frequency wave characteristics in the upstream and downstream regime of the terrestrial bow shock. J. Geophys. Res. Space Phys. 111, A01203 (2006).
Krauss-Varban, D. Bow shock and magnetosheath simulations: wave transport and kinetic properties. Geophys. Monogr. Ser. 81, 121–134 (1994).
Krauss-Varban, D. Waves associated with quasi-parallel shocks: generation, mode conversion and implications. Adv. Space Res. 15, 271–284 (1995).
Engebretson, M. J. et al. A comparison of ULF fluctuations in the solar wind, magnetosheath, and dayside magnetosphere. 1. Magnetosheath morphology. J. Geophys. Res. 96, 3441–3454 (1991).
Engebretson, M. J. et al. The role of the ionosphere in coupling upstream ULF wave power into the dayside magnetosphere. J. Geophys. Res. 96, 1527–1542 (1991).
von Alfthan, S. et al. Vlasiator: first global hybrid-Vlasov simulations of Earth’s foreshock and magnetosheath. J. Atmos. Sol. Terrestrial Phys. 120, 24–35 (2014).
Palmroth, M. et al. Vlasov methods in space physics and astrophysics. Living Rev. Comput. Astrophys. 4, 1 (2018).
Palmroth, M. et al. ULF foreshock under radial IMF: THEMIS observations and global kinetic simulation Vlasiator results compared. J. Geophys. Res. Space Phys. 120, 8782–8798 (2015).
Kempf, Y. et al. Ion distributions in the Earth’s foreshock: hybrid-Vlasov simulation and THEMIS observations. J. Geophys. Res. Space Phys. 120, 3684–3701 (2015).
Pfau-Kempf, Y. et al. Evidence for transient, local ion foreshocks caused by dayside magnetopause reconnection. Ann. Geophys. 34, 943–959 (2016).
Blanco-Cano, X. et al. Cavitons and spontaneous hot flow anomalies in a hybrid-Vlasov global magnetospheric simulation. Ann. Geophys. 36, 1081–1097 (2018).
Takahashi, K. et al. Propagation of ultralow-frequency waves from the ion foreshock into the magnetosphere during the passage of a magnetic cloud. J. Geophys. Res. Space Phys. 126, e2020JA028474 (2021).
Kilpua, E. K. J., Balogh, A., von Steiger, R. & Liu, Y. D. Geoeffective properties of solar transients and stream interaction regions. Space Sci. Rev. 212, 1271–1314 (2017).
Winterhalter, D. & Kivelson, M. G. Observations of the Earth’s bow shock under high Mach number/high plasma beta solar wind conditions. Geophys. Res. Lett. 15, 1161–1164 (1988).
Le, G. & Russell, C. T. A study of ULF wave foreshock morphology - II: spatial variation of ULF waves. Planet. Space Sci. 40, 1215–1225 (1992).
Shevyrev, N. N., Zastenker, G. N., Eiges, P. E. & Richardson, J. D. Low frequency waves observed by Interball-1 in foreshock and magnetosheath. Adv. Space Res. 37, 1516–1521 (2006).
Eastwood, J. P., Balogh, A., Lucek, E. A., Mazelle, C. & Dandouras, I. Quasi-monochromatic ULF foreshock waves as observed by the four-spacecraft Cluster mission: 1. Statistical properties. J. Geophys. Res. Space Phys. 110, A11219 (2005).
Schwartz, S. J. Shock and discontinuity normals, Mach numbers, and related parameters. ISSI Sci. Rep. Ser. 1, 249–270 (1998).
Battarbee, M. et al. Non-locality of Earth’s quasi-parallel bow shock: injection of thermal protons in a hybrid-Vlasov simulation. Ann. Geophys. 38, 625–643 (2020).
Bellan, P. M. Revised single-spacecraft method for determining wave vector k and resolving space–time ambiguity. J. Geophys. Res. Space Phys. 121, 8589–8599 (2016).
Verscharen, D. & Chandran, B. D. G. NHDS: The New Hampshire Dispersion Relation Solver. Res. Notes AAS 2, 13 (2018).
Wu, B.-H., Mandt, M. E., Lee, L. C. & Chao, J. K. Magnetospheric response to solar wind dynamic pressure variations: interaction of interplanetary tangential discontinuities with the bow shock. J. Geophys. Res. 98, 21297–21312 (1993).
Maynard, N. C. et al. Interaction of the bow shock with a tangential discontinuity and solar wind density decrease: observations of predicted fast mode waves and magnetosheath merging. J. Geophys. Res. Space Phys. 112, A12219 (2007).
Thomas, V. A., Winske, D. & Thomsen, M. F. Simulation of upstream pressure pulse propagation through the bow shock. J. Geophys. Res. 100, 23481–23488 (1995).
McKenzie, J. F. & Westphal, K. O. Interaction of hydromagnetic waves with hydromagnetic shocks. Phys. Fluids 13, 630–640 (1970).
Krauss-Varban, D. & Omidi, N. Structure of medium mach number quasi-parallel shocks: upstream and downstream waves. J. Geophys. Res. 96, 17715–17731 (1991).
Matsuoka, A., Southwood, D. J., Kokubun, S. & Mukai, T. Propagation sense of low-frequency MHD waves in the magnetosheath observed by Geotail. J. Geophys. Res. 105, 18,361–18,376 (2000).
Schäfer, S. et al. Statistical phase propagation and dispersion analysis of low frequency waves in the magnetosheath. Ann. Geophys. 23, 3339–3349 (2005).
Pfau-Kempf, Y. et al. On the importance of spatial and velocity resolution in the hybrid-Vlasov modeling of collisionless shocks. Front. Phys. Plasma Phys. 6, 44 (2018).
Torrence, C. & Compo, G. P. A practical guide to wavelet analysis. Bull. Am. Meteorol. Soc. 79, 61–78 (1998).
Roberts, O. W. et al. Possible coexistence of kinetic Alfvén and ion Bernstein modes in sub-ion scale compressive turbulence in the solar wind. Phys. Rev. Res. 2, 043253 (2020).
Gershman, D. J. et al. Energy partitioning constraints at kinetic scales in low-β turbulence. Phys. Plasmas 25, 022303 (2018).
Russell, C. T. et al. The Magnetospheric Multiscale magnetometers. Space Sci. Rev. 199, 189–256 (2016).
Pollock, C. et al. Fast plasma investigation for Magnetospheric Multiscale. Space Sci. Rev. 199, 331–406 (2016).
Narita, Y. et al. Alfven waves in the foreshock propagating upstream in the plasma rest frame: statistics from Cluster observations. Ann. Geophys. 22, 2315–2323 (2004).
Sahraoui, F., Belmont, G., Goldstein, M. L. & Rezeau, L. Limitations of multispacecraft data techniques in measuring wave number spectra of space plasma turbulence. J. Geophys. Res. Space Phys. 115, A04206 (2010).
Robert, P., Dunlop, M. W., Roux, A. & Chanteur, G. Accuracy of current density determination. ISSI Sci. Rep. Ser. 1, 395–418 (1998).
Quest, K. B. Theory and simulation of collisionless parallel shocks. J. Geophys. Res. 93, 9649–9680 (1988).
Sandroos, A. Vlsv: file format and tools. GitHub https://github.com/fmihpc/vlsv/ (2022).
Hannuksela, O. & the Vlasiator team. Analysator: Python analysis toolkit. GitHub https://github.com/fmihpc/analysator/ (2022).
Battarbee, M. & the Vlasiator team. Analysator: Python analysis toolkit. Zenodo https://doi.org/10.5281/zenodo.4462515 (2021).
Palmroth, M. Vlasiator web site. Univ. of Helsinki https://www.helsinki.fi/en/researchgroups/vlasiator/ (2022).
Pfau-Kempf, Y. & the Vlasiator team. fmihpc/vlasiator: Vlasiator 5.1. Zenodo https://doi.org/10.5281/zenodo.4719554 (2021).
Palmroth, M. & the Vlasiator team. Vlasiator: hybrid-Vlasov simulation code. GitHub https://github.com/fmihpc/vlasiator/ (2022).
Acknowledgements
We acknowledge the International Space Science Institute for supporting the ISSI team 448 on ‘Global study of the transmission of foreshock ULF waves into the magnetosheath and the magnetosphere’ led by L. Turc and M. Palmroth, in which this study was initiated. We acknowledge the European Research Council for Starting grant 200141-QuESpace, with which Vlasiator was developed, and Consolidator grant 682068-PRESTISSIMO awarded to further develop Vlasiator and use it for scientific investigations. The CSC-IT Center for Science in Finland is acknowledged for the Grand Challenge award leading to the results shown here. We thank the Finnish Grid and Cloud Infrastructure (FGCI) and specifically the University of Helsinki computing services for supporting this project with computational and data storage resources. The work leading to these results has been carried out in the Finnish Centre of Excellence in Research of Sustainable Space (Academy of Finland grant numbers 336805 and 312390). We thank S. von Alfthan for his major contribution in the development of Vlasiator. We acknowledge financial support from the Academy of Finland (grant numbers 322544 to L.T., 328893 to L.T., 339756 to Y.P.-K., ICT_SUNVAC 335555 to M.B. and Profi4 318913 to U.G.), the University of Helsinki (three-year research grant 2020–2022 to L.T.), the STFC (Ernest Rutherford Fellowship ST/P003826/1 to D.V. and STFC Consolidated Grant ST/S000240/1 to D.V.), the Swedish Civil Contingencies Agency (grant 2016-2102 to A.P.D.), the Swedish National Space Agency (grant no. 2020-00111 to A.P.D.), the EU Horizon 2020 (project SHocks: structure, AcceleRation, dissiPation (SHARP) #101004131 to A.P.D.), the DGAPA/PAPIIT (grant IN105620 to P.K.), NASA (grant NNX17AD34G to K.T.), the NSF (grant 1744609 to K.T.) and the Grants-in-Aid for Scientific Research of the Japan Society for the Promotion of Science (grant numbers 16H06286 and 18KK0099 to N.T.).
Funding
Open Access funding provided by University of Helsinki including Helsinki University Central Hospital.
Author information
Authors and Affiliations
Contributions
L.T. initiated this study during the first meeting in Bern of the ISSI team 448, performed the numerical analysis and the identification of the MMS events, and wrote most of the first draft of the manuscript. O.W.R. carried out the analysis of MMS observations and the determination of the wave properties, and wrote parts of the manuscript. D.V. contributed to the comparison of MMS data with the theoretical wave dispersion solutions. A.P.D., P.K., M.P. and Y.P.-K. provided valuable inputs throughout the development of this study, both on the observational and the numerical analysis. A.J. and M.D. contributed to the analysis of MMS observations. E.K.J.K., J.S., K.T. and N.T. provided valuable input to the study during the ISSI team meeting. Y.P.-K., M.B. and U.G. ran the global Vlasiator simulation used in this study and provided guidance in analysing its outputs. Y.P.-K. ran the 1D Vlasiator simulations. All co-authors proofread and contributed to improving the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare that they have no competing interests.
Peer review
Peer review information
Nature Physics thanks Xuzhi Zhou, Takayuki Umeda and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data
Extended Data Fig. 1 Virtual spacecraft observations near the bow shock.
From top to bottom: (a) magnetic field strength (black) and ion density (red), (b) magnetic field components, (c) wavelet power spectrum of the magnetic field strength, and (d) cross-correlation of the magnetic field strength and density fluctuations. The pink dashed line in panel d shows the predicted foreshock wave period using the Takahashi et al.9 formula. The black contours in panels c and d delineate the 95% confidence interval of the wavelet power spectrum. The virtual spacecraft is positioned in (x = 9.6 RE; y = 0 RE) and is initially located in the foreshock, then crosses the shock shortly after t = 400 s and remains in the magnetosheath afterwards, due to the outward bow shock motion. Between t = 400 − 450 s, fast-mode oscillations at the foreshock wave period are observed just downstream of the shock, consistent with the observations reported by Liu et al.7.
Extended Data Fig. 2 Total pressure variations in the magnetosheath caused by the foreshock waves.
Time-position maps of the magnetosonic Mach number (a) and the total pressure (b) along the Sun-Earth line. The total pressure is calculated as the sum of the thermal pressure and the magnetic pressure. The white contour marks where the magnetosonic Mach number Mms = 1.
Extended Data Fig. 3 Magnetosheath wave activity away from the Sun-Earth line.
Left-hand side: time series of the magnetic field and plasma parameters extracted at a virtual spacecraft positioned in (x = 8 RE; y = -5 RE). The panels show, from top to bottom: (a) the magnetic field strength (in black) and the ion density (red), (b) the magnetic field components, (c) the wavelet power spectrum of the magnetic field strength, (d) the cross-correlation of the magnetic field strength and density fluctuations, and (e) the compressibility of the magnetic field fluctuations, defined as the wave power parallel to the mean magnetic field divided by the total wave power. The pink dashed line in panel d shows the predicted foreshock wave period using the Takahashi et al.9 formula. The black contours in panels c-e delineate the 95% confidence interval of the wavelet power spectrum and the hatched area marks the cone of influence. The time series used for panels c-e have been high-pass filtered to remove low frequency variations due to boundary motion (with a cutoff at 40 s), to better highlight the wave power in the relevant period range. Right-hand side: Dispersion plot obtained from the 2D Fourier transform of the magnetic field strength in the magnetosheath between x = 5 RE and x = 8 RE at y = - 5 RE using unfiltered data to which a Hann window has been applied along both dimensions. On the horizontal axis, the frequencies are normalised to the ion cyclotron frequency Ωci and the wave number, on the vertical axis, to the proton inertial length dp. The solid yellow lines show the Courant-Friedrichs-Lewy condition, which is the maximum speed at which information can travel in the simulation. The median bulk speed in the magnetosheath, at the locations used to calculate the 2D Fourier transform is indicated by the green dashed line. The blue dash-dotted lines and the pink dotted lines indicate sunward and earthward propagation at the median Alfvén speed and median fast-magnetosonic speed, respectively, in the plasma rest frame.
Extended Data Fig. 4 Power spectral density of the foreshock and magnetosheath fluctuations observed by MMS on 14 February 2020.
Power spectral density of the total magnetic field fluctuations calculated for the intervals marked by the vertical lines in Fig. 4 in the foreshock (black) and in the magnetosheath (red). The vertical dashed lines indicate the predicted foreshock wave frequency for each interval from the Takahashi et al.9 formula.
Extended Data Fig. 5 Wave phase speeds from MMS observations.
Total magnetic power as a function of the wave phase speeds vph = ω/k. The phase speeds correspond to the recovered points obtained from the Bellan method and are organised as a function of the orientation of the wavevector, with a negative speed indicating earthward propagation. The data points in black are the same as shown in Fig. 5, for which the current density was obtained from the particle measurements. Those in red are based on the current density calculated using the curlometer method. The results from both approaches are in excellent agreement. The vertical blue line denotes the Alfvén speed. The dashed lines denote the phase speeds for largest and smallest phase speeds \({v}_{{{{\rm{A}}}}}\cos ({\theta }_{{{{\rm{kB}}}}})\) and the shaded cyan area denotes the region where the speeds are consistent with Alfvén waves or advected structures.
Extended Data Fig. 6 Wave activity in a 1D local shock simulation with Alfvén Mach number MA = 6.9.
Proton density (a), magnetic field strength (b) and magnetic field components (c) along the x axis at t = 500 s from the beginning of the simulation. (d) and (e) Dispersion plots obtained from the 2D Fourier transform of the magnetic field strength (d) and By component (e) between x = -10 RE and x = -5 RE (marked by the black bar in panel (b)) and t = 250 - 500 s, using unfiltered data to which a Hann window has been applied along both dimensions. On the horizontal axis, the frequencies are normalised to the ion cyclotron frequency Ωci and the wave number, on the vertical axis, to the proton inertial length dp. The solid yellow lines show the Courant-Friedrichs-Lewy condition, which is the maximum speed at which information can travel in the simulation, the blue dash-dotted lines the median Alfvén speed, and the pink line the median fast-magnetosonic speed in the magnetosheath at the locations used to calculate the 2D Fourier transform. The median bulk speed is indicated by the green dashed line.
Extended Data Fig. 7 Wave activity in a 1D local shock simulation with Alfvén Mach number MA = 4.
Proton density (a), magnetic field strength (b) and magnetic field components (c) along the x axis at t = 500 s from the beginning of the simulation. (d) and (e) Dispersion plots obtained from the 2D Fourier transform of the magnetic field strength (d) and By component (e) between x = -10 RE and x = -5 RE (marked by the black bar in panel (b)) and t = 250 - 500 s, using unfiltered data to which a Hann window has been applied along both dimensions. On the horizontal axis, the frequencies are normalised to the ion cyclotron frequency Ωci and the wave number, on the vertical axis, to the proton inertial length dp. The solid yellow lines show the Courant-Friedrichs-Lewy condition, which is the maximum speed at which information can travel in the simulation, the blue dash-dotted lines the median Alfvén speed, and the pink line the median fast-magnetosonic speed in the magnetosheath at the locations used to calculate the 2D Fourier transform. The median bulk speed is indicated by the green dashed line.
Extended Data Table 1 List of all events used in the spacecraft data analysis, including important parameters for each interval.
Intervals during which the MMS satellites were located in the quasi-parallel subsolar magnetosheath. The upstream interplanetary magnetic field (IMF) vector and cone angle are obtained either from the OMNI data set or directly from measurements from the ACE or Wind spacecraft propagated to the bow shock, depending on the data availability for each event. The IMF cone angle provides a good estimate of the shock θBn angle upstream of MMS, because the spacecraft are located in the subsolar region. The magnetosheath wave period is obtained as the peak of the power spectral density of the magnetic field strength during the interval that is closest to the predicted foreshock wave period, given between parentheses. The last column provides the observed (predicted) foreshock wave period when MMS probed the foreshock shortly before or after the magnetosheath interval. Those events for which no data were available in the foreshock are marked with ‘-’. Note that the solar wind and IMF conditions somewhat differ between the magnetosheath and foreshock intervals, hence the slightly different wave periods, but which show good agreement with the predicted values.
Source data
Source Data Fig. 2
Time series of the virtual satellite measurements extracted at different positions along the Sun–Earth line, used to generate Fig. 2. Each line corresponds to a given time stamp (first column) and a given virtual satellite position (second to fourth columns). Plotted in Fig. 2 are the total magnetic field strength (ninth column, B.magnitude), the proton density (13th column, proton/rho) and the magnetic field components (5th to 8th columns, B.x, B.y, B.z). The other panels of Fig. 2 are obtained from a wavelet transform of the magnetic field strength (c and h) and the wavelet cross-correlation of the magnetic field strength and proton density (d and i) and the wavelet transform of the magnetic field components parallel and perpendicular to the mean magnetic field (e and j).
Source Data Extended Data Fig. 1
Time series of the virtual satellite measurements extracted at x = 9.6RE, y = 0RE, used to generate Extended Data Fig. 1. The time stamps of the virtual satellite measurements are given in the first column. Plotted in Extended Data Fig. 1 are the total magnetic field strength calculated as the norm of the three components (5th to 8th columns, B.x, B.y, B.z), the proton density (9th column, proton/rho) and the magnetic field components (5th to 8th columns, B.x, B.y, B.z). The other panels of Extended Data Fig. 1 are obtained from a wavelet transform of the magnetic field strength (c) and the wavelet cross-correlation of the magnetic field strength and proton density (d)
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Turc, L., Roberts, O.W., Verscharen, D. et al. Transmission of foreshock waves through Earth’s bow shock. Nat. Phys. (2022). https://doi.org/10.1038/s41567-022-01837-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41567-022-01837-z