Articles | Volume 36, issue 3
Ann. Geophys., 36, 695–703, 2018
Ann. Geophys., 36, 695–703, 2018

Regular paper 03 May 2018

Regular paper | 03 May 2018

Plasma flow patterns in and around magnetosheath jets

Plasma flow patterns in and around magnetosheath jets
Ferdinand Plaschke1,2 and Heli Hietala3 Ferdinand Plaschke and Heli Hietala
  • 1Institute of Physics, University of Graz, Graz, Austria
  • 2Space Research Institute, Austrian Academy of Sciences, Graz, Austria
  • 3Department of Earth, Planetary, and Space Sciences, University of California Los Angeles, Los Angeles, CA, USA

Correspondence: Ferdinand Plaschke (


The magnetosheath is commonly permeated by localized high-speed jets downstream of the quasi-parallel bow shock. These jets are much faster than the ambient magnetosheath plasma, thus raising the question of how that latter plasma reacts to incoming jets. We have performed a statistical analysis based on 662 cases of one THEMIS spacecraft observing a jet and another (second) THEMIS spacecraft providing context observations of nearby plasma to uncover the flow patterns in and around jets. The following results are found: along the jet's path, slower plasma is accelerated and pushed aside ahead of the fastest core jet plasma. Behind the jet core, plasma flows into the path to fill the wake. This evasive plasma motion affects the ambient magnetosheath, close to the jet's path. Diverging and converging plasma flows ahead and behind the jet are complemented by plasma flows opposite to the jet's propagation direction, in the vicinity of the jet. This vortical plasma motion results in a deceleration of ambient plasma when a jet passes nearby.

Keywords. Magnetospheric physics (magnetosheath; MHD waves and instabilities; solar wind–magnetosphere interactions)

1 Introduction

Plasma jets defined as transient increases in dynamic pressure are a very common phenomenon in the subsolar magnetosheath (Němeček et al.1998; Savin et al.2008; Karlsson et al.2012; Plaschke et al.2016). In particular, the magnetosheath region downstream of the quasi-parallel shock is permeated by magnetosheath jets. Correspondingly, subsolar jet occurrence is found to be controlled – almost exclusively – by the cone angle of the interplanetary magnetic field (IMF), i.e., the angle between the IMF and the Earth–Sun line, while other solar wind parameters or their variability only play a minor role (Archer et al.2012; Archer and Horbury2013; Plaschke et al.2013). A substantial fraction of jets is believed to originate from bow shock ripples or undulations (Hietala and Plaschke2013), which are common at the reforming, patchy, quasi-parallel bow shock (Schwartz and Burgess1991; Omidi et al.2005). At inclined shock surfaces, the solar wind plasma is compressed, but less decelerated and thermalized, leading to entities of dense and fast plasma inside the magnetosheath (Hietala et al.2009, 2012). Other generation mechanisms are related to IMF discontinuities, discontinuity-related hot flow anomalies (HFAs), and spontaneous HFAs originating from foreshock cavitons (Lin et al.1996a, b; Savin et al.2012; Zhang et al.2013; Kajdič et al.2013; Omidi et al.2013; Archer et al.2014; Chu et al.2017).

Due to their mechanisms of generation, jets are more prevalent closer to the bow shock (Plaschke et al.2013). Nevertheless, large-scale jets with a cross-sectional diameter of over 2 RE (Earth radii) have been found to impact the subsolar magnetopause at relatively high rates (with respect to other transients) of once every 21 min, on average, and once every 6 min when the IMF cone angle is low (under 30). These impact rates are even more remarkable when taking into account that jets of smaller scale occur more often and that typical jet scales are on the order of 1 RE (Karlsson et al.2012; Archer et al.2012; Hietala et al.2012; Plaschke et al.2013, 2016; Gunell et al.2014). When jets impact the magnetopause, the consequences can be substantial. Due to their excess dynamic pressure, they will indent the magnetopause locally, generating surface waves on the boundary and inner magnetospheric compressional waves (Glassmeier and Heppner1992; Plaschke et al.2009; Shue et al.2009; Amata et al.2011; Plaschke and Glassmeier2011; Archer et al.2013a, b). In addition, local reconnection may perhaps be triggered at the magnetopause (Hietala et al.2012). In the magnetosphere, the radiation belt electron population may be modified by magnetopause shadowing (Elkington et al.2003; Turner et al.2012). Consequences may even be seen on the ground, in the form of increased ionospheric convection, geomagnetic field variations, and possibly “throat” aurora observations (Hietala et al.2012; Dmitriev and Suvorova2012; Archer et al.2013b; Han et al.2017).

Consequences of jets are not only restricted to the magnetopause and the magnetosphere. Also within the magnetosheath, jets are expected to alter ambient plasma due to their excess velocity. Simulation results by Karimabadi et al. (2014) show that fast jets push slower ambient magnetosheath plasma out of their way. As a result, that ambient plasma performs an evasive motion around jets, such that it is further slowed down or even pushed in a sunward direction in the vicinity of jets. Thereby, jets may create anomalous flows, stir the plasma in the magnetosheath, and hence be a source of additional turbulence. Note, however, that these simulations were 2-D, which may affect the flow patterns.

A first observational indication for this plasma motion has recently been reported by Plaschke et al. (2017). Within an interval of repeated jet observations downstream of the quasi-parallel bow shock by the Magnetospheric Multiscale (MMS) spacecraft (Burch et al.2016), a high plasma velocity variability is seen, which includes sunward plasma flows in the subsolar magnetosheath. Unfortunately, due to the close MMS spacecraft configuration on the order of a few tens of kilometers (i.e., much smaller than the typical jet scale sizes), it could not be directly proven that the sunward flows were indeed caused by the nearby passage of high-speed jets. This is the primary aim of this paper, to ascertain whether or not jets cause evasive motion of the plasma in the magnetosheath.

2 Data sets

To achieve this goal, we need observations of jets by multiple spacecraft (at least spacecraft pairs) that feature separations on the order of typical jet cross-sectional scales, so that the jet and the plasma outside can be monitored simultaneously. The orbits of the five Time History of Events and Macroscale Interactions during Substorms (THEMIS) spacecraft turn out to be ideally suited (Angelopoulos2008). In particular, the inner three THEMIS spacecraft (THA, THD, and THE) regularly traverse the dayside magnetosheath at the required distances from one another when their orbit apogees are located in the subsolar local time sector.

A data set of 2859 subsolar magnetosheath high-speed jets already exists that we can use for this study; ion velocity, density and dynamic pressure measurements of one of them are shown as an example in Fig. 1. The data set was created for a statistical study by Plaschke et al. (2013), based on 4 years (2008 to 2011) of five THEMIS spacecraft measurements (THA to THE, each spacecraft treated separately as a single spacecraft). Here we briefly recall the steps that led to the compilation of the jet data set; the selection process is explained in more detail in section 2 of Plaschke et al. (2013).

First, THEMIS measurements at geocentric distances between 7 and 18 RE, taken inside a 30 wide cone with the tip at Earth and open towards the Sun, were preselected. From these measurements, magnetosheath intervals were selected: (i) where THEMIS ion density measurements exceeded twice the solar wind density (see panel b of Fig. 1) as given by the OMNI high-resolution solar wind data set (King and Papitashvili2005); (ii) where the differential ion energy flux of 1 keV ions exceeded that of 10 keV ions; and (iii) where THEMIS magnetometer (Auster et al.2008), electrostatic analyzer (McFadden et al.2008), and OMNI IMF, ion density, ion velocity, and plasma beta data were available. Note that the OMNI data are already propagated to the bow shock nose. Averages over 5 min of OMNI measurements preceding any time of interest were used to account for the additional propagation delays to the locations of the THEMIS spacecraft. Therewith, 6960 intervals with durations exceeding 2 min were obtained, comprising in total 2736.9 h of magnetosheath measurements.

Figure 1Figure showing a magnetosheath high-speed jet observed by THC. From top to bottom: (a) ESA ion velocity measurements in GSE, (b) ESA ion density measurements in black and (twice the) solar wind density measurements by OMNI in red (blue), (c) dynamic pressure pd,x measured by THC (black) as well as (1∕4 and 1∕2) the solar wind value in red (green and blue). On top, the pre-jet, jet, and post-jet intervals are marked; normalized times are given in red. Figure modified from Fig. 1 in Plaschke et al. (2013).


Figure 2Cases where second spacecraft observed jets. The top panel (a) shows v and the bottom panel (b) shows Δv values. Median and (upper and lower) quartile values are depicted with solid and dotted lines, respectively. Red lines correspond to reference spacecraft measurements; purple and blue lines correspond to second spacecraft measurements for d<0.4 RE (〈d〉=0.20 RE) and d>0.4 RE (〈d〉=0.51 RE), respectively. Initial median values pertaining to normalized times tn=0…0.4 are depicted by bold circles.


Second, within these subsolar magnetosheath data intervals, jets were identified as follows: (i) the dynamic pressure in the geocentric solar ecliptic (GSE) x-direction (pd,x=ρvx2) shall be larger than half the respective solar wind value pd,sw (see panel c of Fig. 1). The interval where pd,x>pd,sw/4 shall be denoted as a jet interval of length tjet. Start and end times of that interval shall be denoted with tjet,s and tjet,e, respectively, and the time of maximum pd,x/pd,sw shall be denoted with t0 (all these times are marked with vertical lines in Fig. 1). (ii) The observing spacecraft shall be located in the magnetosheath (as defined above) between tjet,s-1 min and tjet,e+1 min. Intervals between tjet,s-1 min and tjet,s and between tjet,e and tjet,e+1 min are called pre-jet and post-jet intervals, respectively (see the top of Fig. 1). (iii) The ion velocity vx shall be negative within the jet interval, and shall surpass vx(t0)∕2 at least once within both pre-jet and post-jet intervals (see panel a of Fig. 1). Therewith we ensure that jets propagate towards the magnetopause and are associated with significant, localized enhancements in velocity. By applying these criteria, Plaschke et al. (2013) identified 2859 jets in total in the preselected magnetosheath data.

Plaschke et al. (2016) go on to identify jet observations by one THEMIS spacecraft (denoted as the reference spacecraft at position rref) for which additional observations by another THEMIS spacecraft (second spacecraft at rsec) are available in a plane perpendicular to the jet propagation direction. That latter direction shall be given by the ion velocity v measured by the reference spacecraft at time t0. In detail, it is ascertained whether there is a second spacecraft such that the angle between

(1) d = r sec - r ref

and v is between 80 and 100 at times t0. The second spacecraft shall be located in the magnetosheath between t0±(tjet+1 min). These criteria are fulfilled for 561 jets. For 101 of those, there are even two “secondary” spacecraft measurements available (i.e., all three THEMIS spacecraft were located appropriately in the magnetosheath). Hence, we obtain 662 observations by pairs of jet observing reference spacecraft and context providing second spacecraft in the subsolar magnetosheath. Plaschke et al. (2016) call the data set comprising these 662 cases the two spacecraft (2SC) data set (see Sect. 3 of that paper for more details).

Figure 3Visualization of velocities depicted in Fig. 2; cases where second spacecraft observed jets. Left panel (a): arrows show median velocity components v (upwards) and Δv (to the right) corresponding to different normalized times tn. Red arrows are based on reference spacecraft measurements. Purple and blue arrows show median v and Δv for second spacecraft at two different distances 〈d〉=0.20 RE and 〈d〉=0.51 RE. Initial values pertaining to normalized times tn=0…0.4 are depicted by bold arrows/circles. Right panel (b): derived from (a) by subtracting those initial flow values, i.e., the respective background flows in the magnetosheath.


3 Analysis and results

A jet is observed by the reference spacecraft in each of the 662 cases; i.e., there are defined times tjet,s-1 min, tjet,s, t0, tjet,e, and tjet,e+1 min. We identify normalized times tn=0…4, such that tn=0 and tn=4 correspond to tjet,s-1 min and tjet,e+1 min, respectively. Times t0 are associated with tn=2. For the jet shown in Fig. 1, the normalized times are given on top of the figure in red.

In order to characterize the flow patterns in and around jets, we use the ion velocities v measured by the pairs of THEMIS spacecraft: (i) in the direction of v as seen by the reference spacecraft (vref), which at tn=2 corresponds to the jet direction as observed by the reference spacecraft, and (ii) in the direction d, which is approximately perpendicular to vref, at least at tn=2. We denote ion velocity measurements along vref with v, and along d as v. The difference of the latter measurements

(2) Δ v = v , sec - v , ref

indicates whether the flow along the inter-spacecraft line is divergent (Δv positive) or convergent (Δv negative).

3.1 Flow pattern inside jets

First, we look into the plasma flow within jets. Therefore, we select those cases of the 2SC data set where second spacecraft observed simultaneously a jet, i.e., where times t0 (tn=2) as seen by the reference spacecraft were within a jet interval observed by second spacecraft. This holds for 230 cases. We divide this group into cases where |d|=d<0.4RE and cases where d>0.4 RE. There are 84 and 146, respectively, and their median distances 〈d〉 are 〈d〉=0.20 RE and 〈d〉=0.51 RE. Furthermore, we divide the normalized time sequence (between tn=0 and 4) into 10 intervals of 0.4 time units each. For each group of cases and time interval, we compute median v and Δv values, as well as upper and lower quartiles thereof. These are displayed in Fig. 2.

In the top panel (a), v is shown as a function of normalized time tn. The median velocity v as seen by the reference spacecraft (red solid line) increases from just over 100 to about 160 km s−1. That velocity is computed from all 662 2SC cases, as it does not depend on second spacecraft observations; the same applies to the reference spacecraft observations shown in Figs. 34, and 5. In normalized time, the maximum in v is symmetric around tn=2. Second spacecraft at different distances d (but still within the jets that the reference spacecraft are in) see corresponding increases in v (purple and blue lines). However, the peak median v values are lower (less than 140 km s−1). This is expected, since we use the ion velocity vectors measured by the reference spacecraft as reference directions for v. Hence, reference/second spacecraft measurements of v are necessarily equal to/lower than the full modulus of the ion velocity at the respective spacecraft locations.

In the bottom panel (b) of Fig. 2, Δv values are shown. Almost all median values (solid lines) are positive. Hence, in general, we observe diverging flows, as expected in the subsolar magnetosheath due to the velocity deviation imposed on the solar wind plasma when passing the curved bow shock. In agreement with expectations, this effect is larger for larger spacecraft separations d, as evidenced by the initial (tn=0…0.4) median Δv values that are depicted by bold circles. With respect to the variations in Δv, a clear bipolar pattern is apparent. Divergence of plasma flows first increases ahead of tn=2. After tn=2, the divergence becomes notably smaller, before returning to pre-jet levels. At tn=2.2 even slightly negative median Δv values (converging flows) are observed, for 2SC cases with d<0.4 RE (〈d〉=0.20 RE, solid purple line).

Figure 4Cases where second spacecraft did not observe jets. The top panel (a) shows v and the bottom panel (b) shows Δv values. Median and (upper and lower) quartile values are depicted with solid and dotted lines, respectively. Red lines correspond to reference spacecraft measurements; purple, blue, and green lines correspond to second spacecraft measurements for 〈d〉=0.44 RE, 〈d〉=0.63 RE, and 〈d〉=1.10 RE, respectively. Initial median values pertaining to normalized times tn=0…0.4 are depicted by bold circles.


The data from Fig. 2 are visualized in Fig. 3a. The median velocities v and Δv yield the vertical and horizontal components of the arrows. Red arrows correspond to reference spacecraft median v observations, purple and blue arrows to the median velocity measurements by second spacecraft. Note that the normalized time increases from the top to the bottom of the figure. Hence, earlier observations corresponding to the front side of the jets are on top, while later observations (rear side) are shown at the bottom.

Figure 3b is directly derived from Fig. 3a by subtracting the background magnetosheath flows, i.e., the initial median values pertaining to normalized times tn=0…0.4. Hence, it depicts and emphasizes the jet-induced changes in plasma flow with respect to the magnetosheath background. Clearly, the change from flow divergence to convergence around tn=2 can be observed (arrows pointing right, then left). This change coincides with an overall flow velocity increase in the direction of vref (upward pointing arrows), not only at the reference spacecraft, but also at the jet-crossing second spacecraft locations.

Figure 5Visualization of velocities depicted in Fig. 4; cases where second spacecraft did not observe jets. Left panel (a): arrows show median velocity components v (upwards) and Δv (to the right) corresponding to different normalized times tn. Red arrows are based on reference spacecraft measurements. Purple, blue, and green arrows show median v and Δv for the second spacecraft at three different distances 〈d〉=0.44 RE, 〈d〉=0.63 RE, and 〈d〉=1.10 RE. Initial values pertaining to normalized times tn=0…0.4 are depicted by bold arrows/circles. Right panel (b): derived from (a) by subtracting those initial flow values, i.e., the respective background flows in the magnetosheath. Note that the arrow scales in panel (b) are different for reference and second spacecraft measurements.


3.2 Flow pattern outside jets

Similar to Fig. 2, Fig. 4 shows velocities v and Δv for median distances 〈d〉=0.44 RE (spacecraft distances d<0.55 RE, 202 cases), 〈d〉=0.63 RE (0.55RE<d<0.85RE, 122 cases), and 〈d〉=1.10 RE (0.85RE<d<1.7RE, 67 cases), derived from spacecraft pair observations where second spacecraft did not observe jets while the reference spacecraft did. Apparently, results for v are very different when comparing Figs. 2a and 4a. Second spacecraft outside jets do not observe an increase in v, but rather a decrease close to tn=2. This decrease is slight in the v medians, but quite pronounced in their upper quartiles.

Again, in general, median v are smaller and median Δv values are larger for larger distances 〈d〉, as evidenced for instance by the values pertaining to tn=0…0.4 that are depicted by bold circles in both panels of Fig. 4.

In the median Δv data pertaining to the smallest distances d (solid purple line in bottom panel b), the bipolar pattern of increasing/decreasing flow divergence over jet passages can be recognized well. This bipolar pattern is also clearly visible in the corresponding upper and lower quartile data, similar to what is seen inside the jet (see Fig. 2b). However, in comparison, that pattern seems to be slightly shifted towards larger normalized times: the minimum in median Δv is reached at tn=3, where it becomes slightly negative (converging flows). With larger distances d from the reference spacecraft, a clear bipolar pattern is not apparent any more. However, we still see an increase in Δv from tn=0 to about tn=1…1.4. That increase starts earlier for more distant spacecraft. Thereafter, Δv essentially decreases to lower levels with respect to tn=0 (bold circles). Similar to Fig. 3, the data from Fig. 4 are visualized in Fig. 5.

Figure 5b depicts how the magnetosheath flow is altered by the nearby passage of a jet. As in Fig. 3b, the divergence and later convergence of plasma flows can be clearly seen. However, in the vref direction, magnetosheath plasma outside of jets is decelerated rather than accelerated around tn=2, as evidenced by slightly downward pointing arrows.

4 Discussion

The results obtained in the previous section allow for the following interpretation. When a high-speed jet penetrates slower magnetosheath plasma, it acts as a plough. The plasma immediately ahead of the fastest jet (core) region (i) will be accelerated in jet propagation direction and thereby contribute to the jet, and (ii) will be pushed to the side (increased divergence of plasma) to make way for the faster jet plasma. This is consistent with increases in both v and Δv around tn=1 as shown in Figs. 2 and 3. The passage of the core jet plasma creates a wake region to which relatively slower plasma is taken in, leading to a minimum in Δv (less diverging or even converging plasma) and a decrease from the maximum in v between tn=2 and 3 (again, Figs. 2 and 3). The situation is illustrated in an exaggerated and simplified manner in Fig. 6: see the thin and thick arrows along the trajectory of the “second spacecraft inside jet” that resemble those in Figs. 3a and b, respectively. Note that Fig. 6, similar to Figs. 3b and 5b, does not show any non-jet-related, i.e., general, differences in v between reference and second spacecraft, as well as generally positive Δv levels that increase with d, which should be attributed to the common divergence of plasma flows in the subsolar magnetosheath.

Figure 6Sketch of plasma flow patterns inside (red) and outside (blue) a magnetosheath jet. The flows are depicted by thin solid arrows resembling those in Figs. 3a and 5a. Differences from the background magnetosheath flow are illustrated by thick arrows and dots in lighter colors (compare to Figs. 3b and 5b). The jet plasma itself is shown to be layered, with slower plasma (in yellow) outside a faster core region (in orange). The upper part corresponds to the magnetosheath ahead of the jet (tn=0), whereas the lower part of the sketch depicts the trailing magnetosheath region (tn=4).


Pushing plasma out of the region ahead of the core jet region needs to have repercussions also on the plasma that is not directly in the jet's way, but in the vicinity of the propagation path. It will also be pushed away from that path ahead of the jet, consistent with an increase in Δv before tn=2 in Fig. 4b. Likewise, behind the jet, plasma streaming into the wake region will also lead ambient plasma outside of the jet's way to follow suit, thereby decreasing Δv after tn=2, in agreement with observations. The plasma motion out/into the jet's way ahead/behind it needs to be closed by plasma motion from the region ahead to the region behind the jet. That plasma motion is opposite to the jet's direction of propagation and to the regular magnetosheath plasma motion, to which it is superposed. As a result, the ambient plasma motion in vref direction (i.e., v) should exhibit a local minimum on the jet passage at tn=2, as illustrated by the thin arrows in Fig. 6 along the trajectory “second spacecraft outside jet”. This can be seen in Fig. 4a. Note, however, that this minimum in median v is very shallow and far from exhibiting the sunward plasma motion (negative v) in the ambient magnetosheath, seen in simulations by Karimabadi et al. (2014) as a reaction to jet penetration. Not even the lower quartiles of second spacecraft v measurements become negative.

In the frame of reference of the background magnetosheath flow (Figs. 3b, 5b, and thick arrows in Fig. 6), the plasma motion in and around jets is clearly vortical. In that respect, jets may be comparable to bursty bulk flows, as those cause somewhat similar plasma flows while ploughing through slower ambient plasma in the tail of the magnetosphere (e.g., Kauristie et al.2000; Birn et al.2004; Keiling et al.2009; Panov et al.2010).

There are also similarities in plasma motion with the laminar flow of an incompressible medium around a sphere or cylinder, in the sub-(magneto)-sonic regime applicable to most jets (Plaschke et al.2013). In agreement with this interpretation, the change around tn=2 from more to less diverging (even converging) flows, between maximum and minimum in Δv, occurs on longer timescales when observing ambient plasma further away from the jet (see Fig. 4b); these changes are fastest inside the jet (see Fig. 2b). In fact, observations of the second spacecraft not traversing the jet and being distant from the reference spacecraft hardly exhibit the full bipolar Δv signature between tn=0 and tn=4 (blue and green lines/arrows in Figs. 4b and 5b). Increases in divergence seem to start earlier than at tn=0. And at tn=4, median Δv values are still not recovering from the decreasing trend. Hence, the plasma motion of jet evasion starts earlier and is concluded later for plasma elements that are further away from the jet's path.

Lastly, it should be mentioned that the evasive plasma motion and, in general, the flow patterns in and around jets may be dependent on particular properties of those jets. For instance, it seems reasonable to assume that the cross-sectional scale size of jets, perpendicular to their propagation direction, may affect the evasive plasma motion, as it correlates with the amount of plasma that needs to be displaced. Unfortunately, it is impossible from two-spacecraft measurements to determine that perpendicular scale size for individual jets. Hence, any possible effect on plasma flow patterns associated therewith cannot be addressed in this paper. These and other effects dependent on jet properties remain to be studied in the future.

5 Summary and conclusions

We have performed a statistical study of plasma flows in and around magnetosheath high-speed jets, based on 662 cases in which one THEMIS spacecraft observed a jet and another nearby THEMIS spacecraft provided context to that observation. We find direct evidence for the following behavior: magnetosheath high-speed jets accelerate slower plasma ahead of them in their propagation path and push it to the side. After passage, plasma that is slower than the fastest jet plasma fills the wake region of the jet along its path. This evasive motion of slower plasma to give way to the jet core plasma is found to occur similarly in the ambient magnetosheath, near but not on a jet's path. In the frame of reference of the background magnetosheath flow, the plasma clearly performs a vortical motion. Divergent and convergent flows out of and towards the jet's path are complemented by flows in the opposite direction to the jet's propagation direction. These flows are superposed to the usual background flow of ambient plasma in the magnetosheath (mostly in the direction of propagation), thereby slowing plasma down. That deceleration is, however, small on average, so that no sunward flows become apparent in the statistical (median) results, contrary to simulations by Karimabadi et al. (2014) and case study observations by Plaschke et al. (2017).

Data availability

Data from the THEMIS mission including level 2 FGM and ESA data are publicly available from the University of California Berkeley and can be obtained from (THEMIS, 2018). The solar wind data from NASA's OMNI high-resolution data set (1 min cadence) are also publicly available and can be obtained from (NASA, 2018).

Competing interests

The authors declare that they have no conflict of interest.


We acknowledge NASA contract NAS5-02099 and Vassilis Angelopoulos for use of data from the THEMIS Mission. Specifically, we acknowledge Charles W. Carlson and James P. McFadden for use of ESA data and Karl-Heinz Glassmeier, Hans-Ulrich Auster and Wolfgang Baumjohann for the use of FGM data provided under the lead of the Technical University of Braunschweig and with financial support through the German Ministry for Economy and Technology and the German Center for Aviation and Space (DLR) under contract 50 OC 0302. We acknowledge valuable discussions within the International Space Science Institute (ISSI) team 350 called “Jets downstream of collisionless shocks” led by the two authors of this paper (Ferdinand Plaschke and Heli Hietala). The work of Ferdinand Plaschke was supported by the Austrian Science Fund (FWF): P 28764-N27. The work of Heli Hietala was supported by NASA grant NNX17AI45G and contract NAS5-02099.
The topical editor, Christopher Owen, thanks two anonymous referees for help in evaluating this paper.


Amata, E., Savin, S. P., Ambrosino, D., Bogdanova, Y. V., Marcucci, M. F., Romanov, S., and Skalsky, A.: High kinetic energy density jets in the Earth's magnetosheath: A case study, Planet. Space Sci., 59, 482–494,, 2011. a

Angelopoulos, V.: The THEMIS Mission, Space Sci. Rev., 141, 5–34,, 2008. a

Archer, M. O. and Horbury, T. S.: Magnetosheath dynamic pressure enhancements: occurrence and typical properties, Ann. Geophys., 31, 319–331,, 2013. a

Archer, M. O., Horbury, T. S., and Eastwood, J. P.: Magnetosheath pressure pulses: Generation downstream of the bow shock from solar wind discontinuities, J. Geophys. Res., 117, A05228,, 2012. a, b

Archer, M. O., Hartinger, M. D., and Horbury, T. S.: Magnetospheric ”magic” frequencies as magnetopause surface eigenmodes, Geophys. Res. Lett., 40, 5003–5008,, 2013a. a

Archer, M. O., Horbury, T. S., Eastwood, J. P., Weygand, J. M., and Yeoman, T. K.: Magnetospheric response to magnetosheath pressure pulses: A low-pass filter effect, J. Geophys. Res., 118, 5454–5466,, 2013b. a, b

Archer, M. O., Turner, D. L., Eastwood, J. P., Horbury, T. S., and Schwartz, S. J.: The role of pressure gradients in driving sunward magnetosheath flows and magnetopause motion, J. Geophys. Res., 119, 8117–8125,, 2014. a

Auster, H. U., Glassmeier, K. H., Magnes, W., Aydogar, O., Baumjohann, W., Constantinescu, D., Fischer, D., Fornacon, K. H., Georgescu, E., Harvey, P., Hillenmaier, O., Kroth, R., Ludlam, M., Narita, Y., Nakamura, R., Okrafka, K., Plaschke, F., Richter, I., Schwarzl, H., Stoll, B., Valavanoglou, A., and Wiedemann, M.: The THEMIS Fluxgate Magnetometer, Space Sci. Rev., 141, 235–264,, 2008. a

Birn, J., Raeder, J., Wang, Y. L., Wolf, R. A., and Hesse, M.: On the propagation of bubbles in the geomagnetic tail, Ann. Geophys., 22, 1773–1786,, 2004. a

Burch, J. L., Moore, T. E., Torbert, R. B., and Giles, B. L.: Magnetospheric Multiscale Overview and Science Objectives, Space Sci. Rev., 199, 5–21,, 2016. a

Chu, C., Zhang, H., Sibeck, D., Otto, A., Zong, Q., Omidi, N., McFadden, J. P., Fruehauff, D., and Angelopoulos, V.: THEMIS satellite observations of hot flow anomalies at Earth's bow shock, Ann. Geophys., 35, 443–451,, 2017. a

Dmitriev, A. V. and Suvorova, A. V.: Traveling magnetopause distortion related to a large-scale magnetosheath plasma jet: THEMIS and ground-based observations, J. Geophys. Res., 117, A08217,, 2012. a

Elkington, S. R., Hudson, M. K., and Chan, A. A.: Resonant acceleration and diffusion of outer zone electrons in an asymmetric geomagnetic field, J. Geophys. Res., 108, 1116,, 2003. a

Glassmeier, K.-H. and Heppner, C.: Traveling magnetospheric convection twin vortices – Another case study, global characteristics, and a model, J. Geophys. Res., 97, 3977–3992,, 1992. a

Gunell, H., Stenberg Wieser, G., Mella, M., Maggiolo, R., Nilsson, H., Darrouzet, F., Hamrin, M., Karlsson, T., Brenning, N., De Keyser, J., André, M., and Dandouras, I.: Waves in high-speed plasmoids in the magnetosheath and at the magnetopause, Ann. Geophys., 32, 991–1009,, 2014. a

Han, D.-S., Hietala, H., Chen, X.-C., Nishimura, Y., Lyons, L. R., Liu, J.-J., Hu, H.-Q., and Yang, H.-G.: Observational properties of dayside throat aurora and implications on the possible generation mechanisms, J. Geophys. Res.-Space, 122, 1853–1870,, 2017. a

Hietala, H. and Plaschke, F.: On the generation of magnetosheath high-speed jets by bow shock ripples, J. Geophys. Res., 118, 7237–7245,, 2013. a

Hietala, H., Laitinen, T. V., Andréeová, K., Vainio, R., Vaivads, A., Palmroth, M., Pulkkinen, T. I., Koskinen, H. E. J., Lucek, E. A., and Rème, H.: Supermagnetosonic Jets behind a Collisionless Quasiparallel Shock, Phys. Rev. Lett., 103, 245001,, 2009. a

Hietala, H., Partamies, N., Laitinen, T. V., Clausen, L. B. N., Facskó, G., Vaivads, A., Koskinen, H. E. J., Dandouras, I., Rème, H., and Lucek, E. A.: Supermagnetosonic subsolar magnetosheath jets and their effects: from the solar wind to the ionospheric convection, Ann. Geophys., 30, 33–48,, 2012. a, b, c, d

Kajdič, P., Blanco-Cano, X., Omidi, N., Meziane, K., Russell, C. T., Sauvaud, J.-A., Dandouras, I., and Lavraud, B.: Statistical study of foreshock cavitons, Ann. Geophys., 31, 2163–2178,, 2013. a

Karimabadi, H., Roytershteyn, V., Vu, H., Omelchenko, Y., Scudder, J., Daughton, W., Dimmock, A., Nykyri, K., Wan, M., Sibeck, D., Tatineni, M., Majumdar, A., Loring, B., and Geveci, B.: The link between shocks, turbulence, and magnetic reconnection in collisionless plasmas, Phys. Plasmas (1994–present), 21, 062308,, 2014. a, b, c

Karlsson, T., Brenning, N., Nilsson, H., Trotignon, J.-G., Vallières, X., and Facsko, G.: Localized density enhancements in the magnetosheath: Three-dimensional morphology and possible importance for impulsive penetration, J. Geophys. Res., 117, A03227,, 2012. a, b

Kauristie, K., Sergeev, V. A., Kubyshkina, M., Pulkkinen, T. I., Angelopoulos, V., Phan, T., Lin, R. P., and Slavin, J. A.: Ionospheric current signatures of transient plasma sheet flows, J. Geophys. Res., 105, 10677–10690,, 2000. a

Keiling, A., Angelopoulos, V., Runov, A., Weygand, J., Apatenkov, S. V., Mende, S., McFadden, J., Larson, D., Amm, O., Glassmeier, K.-H., and Auster, H. U.: Substorm current wedge driven by plasma flow vortices: THEMIS observations, J. Geophys. Res., 114, A00C22,, 2009. a

King, J. H. and Papitashvili, N. E.: Solar wind spatial scales in and comparisons of hourly Wind and ACE plasma and magnetic field data, J. Geophys. Res., 110, A02104,, 2005. a

Lin, Y., Lee, L. C., and Yan, M.: Generation of dynamic pressure pulses downstream of the bow shock by variations in the interplanetary magnetic field orientation, J. Geophys. Res., 101, 479–493,, 1996a. a

Lin, Y., Swift, D. W., and Lee, L. C.: Simulation of pressure pulses in the bow shock and magnetosheath driven by variations in interplanetary magnetic field direction, J. Geophys. Res., 101, 27251–27269,, 1996b. a

McFadden, J. P., Carlson, C. W., Larson, D., Ludlam, M., Abiad, R., Elliott, B., Turin, P., Marckwordt, M., and Angelopoulos, V.: The THEMIS ESA Plasma Instrument and In-flight Calibration, Space Sci. Rev., 141, 277–302,, 2008. a

NASA: Solar wind data from NASA's OMNI high resolution data set, available at:, last access: 23 April 2018. 

Němeček, Z., Šafránková, J., Přech, L., Sibeck, D. G., Kokubun, S., and Mukai, T.: Transient flux enhancements in the magnetosheath, Geophys. Res. Lett., 25, 1273–1276,, 1998. a

Omidi, N., Blanco-Cano, X., and Russell, C. T.: Macrostructure of collisionless bow shocks: 1. Scale lengths, J. Geophys. Res., 110, A12212,, 2005. a

Omidi, N., Sibeck, D., Blanco-Cano, X., Rojas-Castillo, D., Turner, D., Zhang, H., and Kajdič, P.: Dynamics of the foreshock compressional boundary and its connection to foreshock cavities, J. Geophys. Res., 118, 823–831,, 2013. a

Panov, E. V., Nakamura, R., Baumjohann, W., Angelopoulos, V., Petrukovich, A. A., Retinò, A., Volwerk, M., Takada, T., Glassmeier, K.-H., McFadden, J. P., and Larson, D.: Multiple overshoot and rebound of a bursty bulk flow, Geophys. Res. Lett., 37, L08103,, 2010. a

Plaschke, F. and Glassmeier, K.-H.: Properties of standing Kruskal-Schwarzschild-modes at the magnetopause, Ann. Geophys., 29, 1793–1807,, 2011. a

Plaschke, F., Glassmeier, K.-H., Sibeck, D. G., Auster, H. U., Constantinescu, O. D., Angelopoulos, V., and Magnes, W.: Magnetopause surface oscillation frequencies at different solar wind conditions, Ann. Geophys., 27, 4521–4532,, 2009. a

Plaschke, F., Hietala, H., and Angelopoulos, V.: Anti-sunward high-speed jets in the subsolar magnetosheath, Ann. Geophys., 31, 1877–1889,, 2013. a, b, c, d, e, f, g, h

Plaschke, F., Hietala, H., Angelopoulos, V., and Nakamura, R.: Geoeffective jets impacting the magnetopause are very common, J. Geophys. Res., 121, 3240–3253,, 2016. a, b, c, d

Plaschke, F., Karlsson, T., Hietala, H., Archer, M., Vörös, Z., Nakamura, R., Magnes, W., Baumjohann, W., Torbert, R. B., Russell, C. T., and Giles, B. L.: Magnetosheath high-speed jets: internal structure and interaction with ambient plasma, J. Geophys. Res., 122, 10157–10175,, 2017. a, b

Savin, S., Amata, E., Zelenyi, L., Budaev, V., Consolini, G., Treumann, R., Lucek, E., Safrankova, J., Nemecek, Z., Khotyaintsev, Y., Andre, M., Buechner, J., Alleyne, H., Song, P., Blecki, J., Rauch, J. L., Romanov, S., Klimov, S., and Skalsky, A.: High energy jets in the Earth's magnetosheath: Implications for plasma dynamics and anomalous transport, J. Exp. Theor. Phys. Lett., 87, 593–599,, 2008. a

Savin, S., Amata, E., Zelenyi, L., Lutsenko, V., Safrankova, J., Nemecek, Z., Borodkova, N., Buechner, J., Daly, P. W., Kronberg, E. A., Blecki, J., Budaev, V., Kozak, L., Skalsky, A., and Lezhen, L.: Super fast plasma streams as drivers of transient and anomalous magnetospheric dynamics, Ann. Geophys., 30, 1–7,, 2012. a

Schwartz, S. J. and Burgess, D.: Quasi-parallel shocks – A patchwork of three-dimensional structures, Geophys. Res. Lett., 18, 373–376,, 1991. a

Shue, J.-H., Chao, J.-K., Song, P., McFadden, J. P., Suvorova, A., Angelopoulos, V., Glassmeier, K. H., and Plaschke, F.: Anomalous magnetosheath flows and distorted subsolar magnetopause for radial interplanetary magnetic fields, Geophys. Res. Lett., 36, L18112,, 2009. a

THEMIS: THEMIS mission including level 2 FGM and ESA data, available at:, last access: 23 April 2018. 

Turner, D. L., Shprits, Y., Hartinger, M., and Angelopoulos, V.: Explaining sudden losses of outer radiation belt electrons during geomagnetic storms, Nat. Phys., 8, 208–212,, 2012. a

Zhang, H., Sibeck, D. G., Zong, Q.-G., Omidi, N., Turner, D., and Clausen, L. B. N.: Spontaneous hot flow anomalies at quasi-parallel shocks: 1. Observations, J. Geophys. Res., 118, 3357–3363,, 2013. a

Short summary
Fast jets of solar wind particles drive through a slower environment in the magnetosheath, located sunward of the region dominated by the Earth’s magnetic field. THEMIS multi-spacecraft observations show that jets push ambient particles out of their way. These particles flow around the faster jets into the jets’ wake. Thereby, jets stir the magnetosheath, changing the properties of this key region whose particles and magnetic fields can directly interact with the Earth’s magnetic field.