CME–HSS Interaction and Characteristics Tracked from Sun to Earth

In a thorough study, we investigate the origin of a remarkable plasma and magnetic field configuration observed in situ on June 22, 2011, near L1, which appears to be a magnetic ejecta (ME) and a shock signature engulfed by a solar wind high-speed stream (HSS). We identify the signatures as an Earth-directed coronal mass ejection (CME), associated with a C7.7 flare on June 21, 2011, and its interaction with a HSS, which emanates from a coronal hole (CH) close to the launch site of the CME. The results indicate that the major interaction between the CME and the HSS starts at a height of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$1.3~\mbox{R}_{\odot }$\end{document}1.3R⊙ up to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$3~\mbox{R}_{\odot }$\end{document}3R⊙. Over that distance range, the CME undergoes a strong north-eastward deflection of at least \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$30^{\circ }$\end{document}30∘ due to the open magnetic field configuration of the CH. We perform a comprehensive analysis for the CME–HSS event using multi-viewpoint data (from the Solar TErrestrial RElations Observatories, the Solar and Heliospheric Observatory and the Solar Dynamics Observatory), and combined modeling efforts (nonlinear force-free field modeling, Graduated Cylindrical Shell CME modeling, and the Forecasting a CME’s Altered Trajectory – ForeCAT model). We aim at better understanding its early evolution and interaction process as well as its interplanetary propagation and related in situ signatures, and finally the resulting impact on the Earth’s magnetosphere.


Introduction
Coronal mass ejections (CME) frequently pass over Earth at an average rate of 1 -2 events per month but with significant variations throughout the solar cycle (Richardson and Cane, 2010). They are the major cause of strong geomagnetic effects, especially during solar maximum (e.g., Farrugia, Burlaga, and Lepping, 1997;Liu et al., 2014). Besides CMEs, stream interaction regions (SIR) and, if persistent for several solar rotations so-called co-rotating interaction regions (CIR), structure interplanetary space. It is the interaction between solar wind high speed streams (HSS), emanating from coronal holes (CH), and the slow solar wind ahead, which forms compression regions, shocks and rarefaction regions, causing recurrent geomagnetic effects on Earth (see, e.g., Alves, Echer, and Gonzalez 2006, Verbanac et al. 2011, Vršnak et al. 2017, Yermolaev et al. 2018, Richardson 2018. While HSSs pose a continuous outflow, CMEs abruptly disrupt the rather steady solar wind structure, causing deviations (preconditioning) from the quiet solar wind conditions over the duration of several days (Temmer et al., 2017).
The interaction of CMEs with the solar wind, especially HSSs, may significantly change the CME properties en route through the heliosphere. The embedded flux rope may deform, kink or rotate (Manchester et al., 2004;Riley and Crooker, 2004;Wang et al., 2006;Yurchyshyn, 2008;Isavnin, Vourlidas, and Kilpua, 2013), erode due to reconnection (Dasso et al., 2006;Ruffenach et al., 2012Ruffenach et al., , 2015Lavraud et al., 2014;Wang et al., 2018), be deflected (Wang et al., 2004(Wang et al., , 2014a(Wang et al., , 2016Manchester et al., 2005;Evans, 2013, 2015;Zhuang et al., 2019) and may be related to increased turbulence in the sheath region (Lugaz et al., 2015;Kilpua, Koskinen, and Pulkkinen, 2017). Fast solar wind may also cause the CME to speed up, hence, shortening the propagation time between Sun and Earth. Nieves-Chinchilla et al. (2012) investigated in detail the rotation of a CME as it propagates through the heliosphere. It is understood that changes in CME properties differ strongly for processes taking place already low in the corona compared to those happening in interplanetary space (Wang et al., 2014b;Winslow et al., 2016). All these effects alter the initial CME properties observed close to the Sun making predictions of arrival time and geoeffectiveness a complex endeavor (Richardson, 2018). This also shows that the ambient solar wind plays an important role for CME propagation, and that it is necessary to study and understand the interaction processes already from its source on the Sun and related 1 AU signatures.
Previous studies of CME-HSS interaction events often focused on HSSs catching up with a CME deforming and compressing it (Winslow et al., 2016;He et al., 2018). In this study the CME appears to be propagating behind the SIR within the HSS, giving rise to a number of peculiarities in the observed in situ signatures. By studying in detail this CME-HSS interaction event we aim to unravel the complex physical processes related to a CME propagating in a HSS starting from the CME eruption site close to the HSS related CH. To do so, we investigate the Sun-Earth chain of a distinct and well-observed CME-HSS interaction event combining remote sensing observations, in situ measurements at 1 AU and modeling efforts. We further investigate the effects on the Earth's magnetosphere as a consequence of the CME-HSS interaction.

Motivation
Starting on June 22, 2011, we observed an intriguing configuration of plasma and magnetic field in situ signatures near L1. We find a SIR signature with a clear stream interface (SI) followed by a shock within the HSS on June 23, 2011 (Figure 1). After a short standoff distance, the magnetic ejecta (ME) can be identified propagating with the same speed as the ambient solar wind of the HSS. This gives rise to the two major questions which we address in this study: • How can these unique shock characteristics be explained? Usually, a shock within a high speed environment, like a HSS, would either quickly propagate through or would dissipate in the case it is not driven. The intuitive solution would be that the shock is driven by the ejecta. • Is the shock driven? Judging from the short standoff distance between the shock and the ejecta and that the shock signature is found in the middle of the HSS, one would suspect so. However, the speed of the ejecta is equal to the speed of the HSS, but the speed of the shock is higher, which indicates that it is not driven at that distance.
To answer these questions we track the CME from 1 AU back through the interplanetary space to its origin on the solar disk and examine what processes may produce such in situ signatures and what geomagnetic effects they cause. The paper is structured as follows: In Section 3 the in situ signatures are investigated. In Section 4 the CME launch and propagation is analyzed using remote sensing data and various modeling efforts. Section 5 presents the resulting geomagnetic effects. In Section 6 the results are discussed and summarized in Section 7.

In Situ Signatures
Using OMNI 1 data from the Advanced Composition Explorer (ACE; Stone et al. 1998) and the Global Geospace Science Wind satellite (Acuña et al., 1995), we investigate the properties of the solar wind and CME plasma at a distance of about 1 AU. Figure 1 shows a 24-hr period from 18 UT, June 22 to 18 UT, June 23, 2011. From top to bottom, we show the proton number density and the alpha-to-proton number density ratio, the bulk speed, the temperature and expected proton temperature for normal solar wind expansion (Lopez, 1987), the total perpendicular pressure (Russell, Shinde, and Jian, 2005) and dynamic pressure including the α particles, the components of the magnetic field in GSM coordinates, the total magnetic field, and the proton β and the Alfvèn Mach number. From the proton bulk speed one can identify a slow stream followed by a fast one. The stream interface occurs at around 20 UT, June 22 shown by the orange vertical guideline. The SI was identified based on the sharp drop of the density compression and the sharp rise in the temperature as shown in Figure 1. Additionally we observe a gradient in the flow, where the east-west flow changes sign (not shown). The following high speed solar wind stream is interrupted at 03:06 UT by simultaneous sharp rises in the density, temperature, speed and magnetic field, corresponding to a shock structure. Its strength is moderate, with a density ratio of ∼ 2 and in situ a field compression of ∼ 1.5. By using the coplanarity theorem (Colburn and Sonett, 1966;Abraham-Shrauner, 1972;Abraham-Shrauner and Yun, 1976), Figure 1 In situ measured Wind and ACE data from the OMNI database of a 24 h time interval during the HSS and CME arrival starting at June 22, 2011, 18 UT. From top to bottom: Proton density (black) with the α-particle ratio overlayed (blue); proton velocity; proton temperature (black) with the expected temperature in red; total perpendicular pressure (black) and the dynamic pressure (purple); magnetic field components (panel 5 -8); plasma-β (black) and the Alfvénic Mach number (red). The vertical guidelines represent the start times of the stream interface (orange), the shock (green) as well as the interval identified as the magnetic ejecta (blue lines) as derived from the data in this figure. which assumes that the magnetic field on both sides of the discontinuity and the shock normal all lie in the same plane, we calculate the shock normal. The shock normal speed we derive using the mass conservation equation across the discontinuity (e.g., see Paschmann and Daly, 1998, and the references therein). Using the magnetic coplanarity (Colburn and Sonett, 1966) we calculate the shock normal velocity to be ∼ 710 km s −1 . The angle between the upstream magnetic field and the shock normal (θ Bu,n ) is ∼ 30 • which makes it a quasi-parallel shock. When using the mixed mode coplanarity (Abraham-Shrauner, 1972) to calculate the shock normal, we derive v shock ∼ 678 km s −1 with θ Bu,n ∼ 34 • . Behind the shock the elevated high speed continues for ∼ 3.75 hrs. Then a structure with high magnetic field strength, low proton temperature, low beta and Alfvén Mach number, and increased alpha-to-proton number density ratio is present. Figure 2 presents the pitch angle distribution of suprathermal electrons in the ACE 272 eV channel showing isotropic flux during the shock arrival and the transition from the shock-sheath to the magnetic ejecta, which is expected. In this interval, we identify as the magnetic ejecta the period during which the electrons show bidirectionality, suggesting a closed magnetic structure (e.g., Montgomery et al. 1974, Gosling et al. 1987, Carcaboso Morales et al. 2018. This strongly supports the interpretation of a flux rope within the HSS. We note that the bulk speed before the shock and after the shock-sheath is the same, which suggests that the magnetic ejecta is currently not driving the shock. This can be interpreted as follows: the CME is embedded and dragged along with the co-rotating interaction region (CIR), which is an intriguing feature.

Remote Observations & Modeling
Using white-light as well as extreme ultra-violet (EUV) observations obtained by the Sun Earth Connection Coronal and Heliospheric Investigation (SECCHI; Howard et al. 2008) suite on-board the Solar TErrestrial RElations Observatories (STEREO-A/B; Kaiser et al. 2008) we are able to track the CME back to its origin on the solar disk. Figure 3 shows the propagation of the CME along the equatorial cut as seen by STEREO-A's Heliospheric Imager 1 and 2 (HI1, HI2) as well as its coronagraph (COR2) obtained in the form of a Figure 3 J-map of the ICME as seen in COR2, HI1 and HI2. The ICME kinematic is marked by the red dotted line and in situ arrival of the shock by the red dashed line. This feature is manually tracked back to the corresponding CME launch time at the Sun around 02 UT on June 21, 2011. J-map. We estimate the launch time at the Sun to be around June 21, 2011, 02 UT. This date coincides very well with an observed gradual C7.7 flare from the active region (AR) 11236 (location: N17/W12) starting at 01:22 UT on June 21, 2011. The flare can be associated with an Earth-directed halo CME whose signatures are subsequently observed by the Solar Dynamics Observatory (SDO; Pesnell, Thompson, and Chamberlin 2012), the Solar and Heliospheric Observatory (SOHO; Domingo, Fleck, and Poland 1995) and both STEREOs. A medium-size coronal hole (CH) is located in the south-west of the AR that we associate with the CIR observed in situ near 1 AU on June 22, 2011. Figure 4 shows the GOES soft X-ray flux around the time of the flare (top), which features a very gradual increase over nearly 2 hrs before the peak in the soft X-ray flux (1 -8 Å) of 7.7 × 10 −7 W/m 2 is reached at around 03:25 UT. The panels (from left to right) show the solar corona in the 211 Å EUV filter taken by the Atmospheric Imaging Assembly (AIA; Lemen et al. 2012) on-board SDO, before, during and after the flare. Before the flare starts, the AR and the CH to the south can be well observed. As the flare progresses, the post-eruptive loop system and coronal dimming regions evolve related to the magnetic field restructuring and density depletion due to the erupting flux rope structure (see e.g., Dissauer et al., 2019).

NLFFF Model
To qualitatively describe the configuration of the corona that leads to the non-radial evolution of the CME, we study the ambient magnetic field configuration of the Sun. We use synoptic vector magnetograms from the Heliospheric and Magnetic Imager (HMI, Schou et al. 2012, Couvidat et al. 2016) on-board SDO as input for a global nonlinear force-free field (NLFFF) model. The method was originally proposed in Cartesian geometry by Wheatland, Sturrock, and Roumeliotis (2000), but here we use the spherical optimization code developed by Wiegelmann (2007) and adjusted for the use of synoptic vector maps in Tadesse et al. (2014). As boundary condition we use a synoptic vector map for Carrington rotation 2111, which has been observed between June 05, 2011, and July 03, 2011.
The left panel of Figure 5 shows the resulting two-dimensional open field map at 1 R (f (θ, φ)) where yellow underlying contours corresponds to positive and blue underlying contours to negative polarity footpoints of open field lines. Overlayed is the NLFFF-model (Br) at a surface height of r = 2.5 R with red representing positive polarity and green representing negative polarity. We see that the northern hemisphere has a negative polarity and, vice versa, the southern hemisphere has a positive polarity. Near the flux center of the flare we find open field footpoints to the south which are well observed in EUV as the CH (right panel of Figure 5). The open field footpoints to the west, indicative of the extension of the CH, might be obscured as a consequence of nearby bright loops, and they are therefore not well visible as a dark region in EUV. The identified open field footpoints most likely form a magnetic potential barrier towards south and west, which the CME cannot easily cross due to the frozen-in condition in the coronal plasma. However, towards the north and east the CME can freely expand.

Coronal Dimmings and CME Launch
After establishing the conditions of the global solar magnetic configuration around the AR of interest, we analyze the evolution of coronal dimmings associated with the CME footpoints anchored at the surface. Coronal dimmings are regions of strongly reduced emission observed in EUV and SXRs (Hudson, Lemen, and Webb, 1996;Sterling and Hudson, 1997;Thompson et al., , 2000. They are a signature of the density depletion that is caused by the plasma expansion and evacuation during the early CME eruption. Bipolar coronal dimmings are generally interpreted to represent the footprints of CMEs in the low corona (e.g. Thompson et al. 2000, Dissauer et al. 2019) and can be used as proxies for studying the initial CME behavior.
In order to properly track the dimming regions, we use a thresholding technique applied on logarithmic base-ratio EUV images. This also allows us to identify so-called secondary dimmings, which are mapping the overlying magnetic field that is expanding and erupting. They are of special interest, since they mostly reflect the propagation direction of the CME in its early evolution phase. To capture the full extent of coronal dimmings over time, we derive cumulative dimming masks. These masks contain all dimming pixels that are detected below a certain threshold and over a given time range. In this way different parts of the dimming that may grow and recover, e.g. due to the associated flare, are also included in the detection. For further details on the method we refer to Dissauer et al. (2018Dissauer et al. ( , 2019. Figure 6 shows the evolution of the CME-associated dimming region, where each pixel is color-coded by the time of its first detection. Dark blue pixels are detected earlier than light blue pixels, respectively. The three contours represent the time evolution of the dimming region during its main impulsive phase, while the CME is propagating below 2 R , at 02:10 UT (red), 02:34 UT (green), and 02:58 UT (magenta), respectively. The dimming mainly grows/spreads towards the south-east and south until 03:00 UT (marked as black arrows) in the reference frame of the flux center of the flare (yellow asterisk). This clearly shows that the CME was launched from the southern part of the AR, evolving into the direction of the CH. The eastward and south-westward spread may be caused by the CME lateral expansion. Figure 7 gives the stereoscopic limb view of the CME using STEREO-A EUVI 195 Å running-difference images. We observe a clear southward propagation direction of the CME during its launch, which is in agreement with the early evolution of the associated coronal dimming. The red arrows indicate the observed propagation direction of the CME apex in the lower corona, which was estimated visually. It follows the trend shown in the latitudinal propagation profile of Figure 8. Following the CME evolution further on, one can clearly see that the CME gets deflected towards north over the distance range 1.3 -3 R , indicating an interaction with the open field from the CH.

Figure 6
Evolution of the coronal dimming caused by the CME. Each pixel is color-coded by the time of its first detection, where darker pixel represent an earlier detection time than lighter ones. The contours represent the size of the dimming region at three timesteps (as indicated by the colored lines in the colorbar).
The contours are at 02:10 UT, 02:34 UT and 02:58 UT, which represent the impulsive evolution phase of the CME (up to 1 R above the solar surface). The arrows indicate the major evolution direction of the dimming. The image is centered at the flux center of the flare (yellow asterisks).

ForeCAT and GCS
To further study the propagation of the CME we use two modeling approaches, the Forecasting a Coronal mass ejection's Altered Trajectory (ForeCAT; Evans 2013, 2015) and the Graduated Cylindrical Shell (GCS; Thernisien, Howard, and Vourlidas 2006, Thernisien, Vourlidas, and Howard 2009, Thernisien 2011 model. ForeCAT uses a Potential Field Source Surface model (PFSS; e.g. see Schatten, Ness, and Wilcox 1968) of the solar magnetic field to calculate the propagation, deflection and expansion of a CME on the basis of global magnetic pressure and tension. For further information on the model we refer to Kay, Opher, and Evans (2015).
Using the GCS model, we reconstruct the CME flux rope (FR) for 9 timesteps between 02:30 UT and 04:45 UT using COR1 and COR2 images from the SECCHI suite on both STEREO spacecraft in combination with LASCO C2 and C3 images (Brueckner et al., 1995). For the GCS reconstruction, which was fitted to best represent the white-light features of consecutive time steps, we derived large changes for latitude and longitude. This clearly shows the non-radial ejection and further deflection behavior of the CME due to its interaction with the open field of the HSS in the corona. Table 1 shows the fitted GCS parameters. Figure 8 shows results of ForeCAT ensemble modeling (Kay and Gopalswamy, 2018) in comparison to the GCS FR parameters. The two top panels show the latitude of the CME axis and the two bottom panels show the longitude in Stonyhurst coordinates. We run an ensemble of 100 ForeCAT simulations with slight variations in the initial CME position and orientation. The blue line shows the seed value of the ensemble -the initial value that determines the center of the ensemble parameter range and our best guess at the true initial  position and orientation. The dashed lines represent the median values, and the dark gray regions are one standard deviation about that. The red dots are the GCS values with error bars of ± 5 • for latitude and ± 10 • for longitude. The left panels show the results of the standard ForeCAT model, whereas in the right panels we included lateral overexpansion of the CME as well as artificial scaling of the magnetic field to simulate the compression of the field lines of the HSS. To simulate the compression we scaled the magnetic field strength in the direction of the CH by a scaling factor S = 1 + (R − R 0 )/0.15 R which was estimated empirically, where R 0 is the initial height of the CME nose. We find the model CME to propagate from the northern hemisphere southwards and from west to east. After the initial phase the model suggests a near constant propagation direction at a latitude of ∼ 7 • north and a longitude of ∼ 10 • east. Starting at a height of about 2.3 R GCS reconstructions are available. The GCS model parameters are very similar to the ForeCAT results. There is a slight eastward offset in the longitude but within the error bars. Only in the latitude we do derive a significant difference between the two results: the GCS reconstruction shows the values continuously rising (northward motion of the CME due to interaction with the HSS) which we cannot derive from the results of the standard ForeCAT model (Figure 8, left panel), which uses a static background magnetic field to model the deflection of a CME. The interaction of a CME and a HSS, however, involves dynamic processes in which, among others, the open magnetic field structure of the HSS will be compressed and the magnetic pressure increases to a point where the CME can no longer expand or propagate in that direction, which could lead to a change in the CME's trajectory. We see this northward motion is reproduced when we include the overexpansion and external compression in ForeCAT.

Flux Rope Evolution Using GCS Results
Using the last GCS reconstruction at 04:45 UT on June 21, and the near 1 AU measured in situ signatures of the magnetic ejecta (duration, average speed) we estimate the expansion factor of the FR in interplanetary space (IP). Assuming self-similar expansion at an arbitrary rate constrained by the initial (GCS) and final size (in situ). We neglect the effects of the CME-HSS interaction in IP space on the CME trajectory and shape due to low plasma densities and magnetic fields. Based on the GCS fitting results, we estimate that the initial FR radius at 13 R is 3.7 R . From the in situ measurement at 1 AU, we estimate the FR radius to be 31 R . The FR radius we estimate is the cross section along the Sun-Earth line, enabling us to investigate the geoeffective part of the magnetic structure. We note that this is the length along the observational path which means that it may not correctly represent the FR radius. The FR radius may be misrepresented in particular when the cross section is a skimming trajectory along the edge or a trajectory parallel to the FR axis. As we cannot clearly define the trajectory of the spacecraft through the magnetic structure (e.g. force-free flux tube fitting results are inconclusive due to the low magnetic field strength in the magnetic ejecta), we assume that the observational path is a reasonable estimate of the FR radius but we also have to consider that the uncertainties are large. Using a power-law equation for the increase of FR size with heliospheric distance given by Dumbović et al. (2018) adapted from a more general expression by Démoulin et al. (2008) and constrained by the derived GCS and in situ results we estimate the power-law index n A , i.e. the expansion factor to n A = 0.51 −0.13 +0.14 . We note that this resulting value lies on the lower end of the range obtained by statistical studies (e.g., see Gulisano et al., 2012; and indicates a rela-tive slow increase in FR size or that the in situ spacecraft passed through only a smaller part of the FR. Additionally, assuming a constant axial magnetic flux we can estimate the drop-rate of the central magnetic field strength, assuming that it follows a power-law behavior with the power-law index n B (see Dumbović et al., 2018). We expect n B to be 1.02 −0.26 +0.28 , which again is on the lower end of the range observed in statistical studies (Gulisano et al., 2012) and indicates a quite slow drop of the magnetic field strength within the FR. With a magnetic field strength of 10 nT at 1 AU we estimate the initial magnetic field strength at 13 R to be 0.002 G and the toroidal magnetic flux at 1 AU to be φ = 1.63 −0.63 +0.85 × 10 20 Mx. This value is at the lower end of the expected range of typical ICME fluxes of 10 20 to 10 22 Mx for C to X-class events (e.g., DeVore, 2000;Qiu et al., 2007;Temmer et al., 2017). With such small expansion, we would expect the B field at 1 AU to be high (since the CME did not expand much) but in fact, B 1AU 10 nT, which is low. So either the initial field was very low to begin with and the in situ spacecraft only encountered part of the FR or n A does not reflect n B .

CME Kinematics in the Sun-Earth Line
After analyzing the initial solar configuration, the launch and the start of the CME-HSS interaction using observations and models, as next step we investigate the propagation of the CME in the direction of Earth. We derived the height-time profile of the CME from multiple sources. We manually tracked the CME front in the equatorial plane using EUVI, COR1, and COR2 image data separately from STEREO-A and -B spacecraft. In addition, we obtained the height-time profile from the intersection of the front of the GCS reconstructed shell with the Sun-Earth line. Figure 9 shows the derived CME kinematics (top to bottom: height-time, velocity and acceleration profile). The dots are the measurements (height-time) and their direct numerical time derivatives (velocity, acceleration). To obtain robust estimates of the corresponding velocity and acceleration profiles, we first smooth the height-time curves and then derive the first and second time derivatives. The smoothing algorithm is based on the method presented in Podladchikova, Van der Linden, and Veronig (2017), extended toward non-equidistant data. The algorithm optimizes between two criteria in order to find a balance between data fidelity, i.e., the closeness of the approximating curve to the data, and smoothness of the approximating curve. From the acceleration profiles obtained in this way, we then interpolate to equidistant data points based on minimization of the second derivatives, and reconstruct the corresponding velocities and height profiles by integration (solid lines). The estimation errors of kinematic profiles are obtained by representing the reconstructed CME height, velocity and acceleration as an explicit function of original CME height-time data in the assumption that STEREO-A and STEREO-B height errors are 1.5% of height, and 3% of height for GCS data. In the height-time diagram we derive that the CME is propagating from east (closer to STEREO-A) towards west (getting closer to STEREO-B) because of the initial offset between the black and blue line which cross later. The GCS model results are in good agreement with the direct measurements. The fit to the GCS height-time profile was done assuming that the acceleration has already subsided.
The velocity and acceleration profiles show the impulsive phase around 02:20 UT followed by a slower acceleration, which subsides over time. Beyond a height of about 3 R the acceleration phase is finished and the CME propagates outwards at a speed of 800 -1000 km s −1 which then seems to adjust to the speed of the surrounding HSS. Near Earth, the in situ measured speed of the magnetic ejecta was ∼ 600 km s −1 which was roughly equal to the speed of the HSS (see Figure 1).

Effects on the Magnetosphere
Finally we analyze the geomagnetic response. In Figure 10, we show the sym-H storm index, the auroral electrojet AE and AL indices and the PCN-north index. The sym-H panel also shows the correction due to the magnetopause currents (blue trace). The sym-H index shows fairly quiet conditions. Interestingly, correcting for the contribution of the magnetopause currents (sym-H*) almost doubles the effects, leading to a moderate storm at the high speed stream before the shock. Basically this is the effect of the relatively high dynamic pressure (see Figure 1). Another moderate storm is caused by the CME sheath.
The AL and AE indices show signatures of a substorm at about 01 UT, June 23. At this time the ionospheric convection shows an increase (PCN index). This convection enhancement is likely being contributed from the nightside source of substorm activity (Sandholt, Andalsvik, and Farrugia, 2012). Some intermittent substorm activity is evident at the end of the interval. As Figure 1 shows, the dynamic pressure is higher than typical (i.e. 2 nPa). In addition, the B z component oscillates with large amplitudes (about 5 nT) in the sheath region of the CME. These two factors are bound to have an effect on the magnetopause (MP) shape when we consider the Shue et al. (1998) model which includes both the compression/rarefaction due to the dynamic pressure (P dyn ) and the erosion due to the negative B z component of the magnetic field. This result is confirmed in Figure 11, which shows from top to bottom the temporal profile of subsolar magnetopause position, the dawn-dusk terminator and their ratio, i.e. the aspect ratio of the magnetosphere. The subsolar distance (red) is closer to Earth than typical (i.e. 12 R E ). A large Earthward shift occurs at shock arrival. Thereafter, the magnetosphere generally expands slowly. The aspect ratio is typically 1.5. However, it changes significantly during the CME sheath passage visualizing the flaring of the magnetosphere. We conclude that the major effects were those produced by the sheath region.

Discussion
In this paper, we present a detailed case study of a CME ejected on June 21, 2011, using multi-viewpoint remote sensing and in situ observations supported by modeling to better understand and explain the intriguing in situ solar wind signatures observed as a consequence of the CME interacting with a nearby CH. We find that the local and global magnetic field configuration, especially large-scale open magnetic structures such as CHs, have a major effect on the early propagation direction of the CME. The low-lying, local magnetic field configuration of the AR and its vicinity first leads to a non-radial ejection southwards. The global magnetic field structure, which in the close proximity is dominated by the open field lines of the CH, causes the subsequent deflection of the CME north-eastwards.
From ForeCAT model results, we find that the open magnetic structures southwards and to the west of the AR seem to exert magnetic pressure on the CME, preventing it from propagating and expanding in these directions. The ForeCAT ensemble results show a larger spread in latitude than longitude, this is because the ensemble runs are based on how well the initial position within an active region can be identified. The spread in latitude and lon-gitude in the corona is simply a result of how sensitive the model is to the precise initial position, and in this case the magnetic forces consistently give a strong eastward deflection but the latitudinal motion is slightly more varied. The resulting propagation direction is also supported by the NLFFF model results. Cremades, Bothmer, and Tripathi (2006), Gopalswamy et al. (2009), Mäkelä et al. (2013 found that there is a correlation between a CME's direction and the properties of nearby CHs (distance, area and mean magnetic field strength). For the CME under study, we find that the direction in which the CME is ejected is mainly affected by the configuration of the local magnetic field at low heights and the position of the flux rope (as estimated from EUV observations) rather than the nearby CH. The initial propagation of the CME to the south can be well observed off-limb in STB EUV images. The on-disk observations from the evolution of the coronal dimming associated with CME show the same behavior and can therefore be used as a proxy to derive the CME propagation direction (see also Dissauer et al., 2019). A similar conclusion is drawn by Mandrini et al. (2007), who analyzed the coronal dimming of the X17 flare event on October 28, 2003.
Starting at a height of ∼ 1.3 R effects of the open field of the CH on the CME are clearly revealed. As the CME runs into the open field of the CH, the CME's further expansion and propagation in the direction of the open field is prohibited, resulting in a deflection towards north and east. This can be interpreted such that the southward motion is actually stopped and reversed by the compressed magnetic field of the HSS. While in the ForeCAT model this occurs as a net northward motion of the entire structure, in reality, where the CME is not forced to maintain a uniform shape, it could manifest as an asymmetric expansion in the northward direction. Separating these effects requires comparison of the external deflection forces and the internal forces that maintain a coherent CME structure, which is beyond the scope of this work. Wang et al. (2011) defined three groups of deflected CMEs: asymmetrical expansion, non-radial ejection, and deflected propagation. Our event can be classified as a combination of two different processes. Firstly, the CMEs initial direction deviates from the radial direction and can be classified as a non-radial eruption (second group). And secondly, due to the presence of an ambient magnetic structure the CME is deflected and as such belongs to the third group (deflected propagation).
The in situ measured data near and at Earth show interesting signatures of the interaction between a SIR and a CME that occurred in the interplanetary space between the Sun and Earth. On June 22, 2011, at around 20:00 UT a clear stream interface signature can be identified, followed on June 23, 2011, at 03:15 UT by the HSS that also reveals a clear shock signature. The shock coincides with typical CME signatures starting at 06:30 UT, hence, can be attributed to be CME related. This, however, raises the question how that shock signature was produced. As shown in Figure 1, the speed of the CME matches the bulk speed of the HSS before the shock arrival which suggests that the shock is not driven by the CME at this point. But the short duration of the shock-sheath region of about 3 hrs suggests otherwise. Due to the low compression of the magnetic field (∼ 1.5) and of the density (∼ 2) we can characterize the shock to have medium strength. We calculated the shock normal speed to be ∼ 700 km s −1 and found that the shock can be considered quasi-parallel.
Due to the inconclusive results in deriving the flux rope geometry using the Lundquist force-free flux rope model we are not able to reliably determine the part of the structure which the spacecraft intersected. However, we can estimate a trajectory from the modeled size and the propagation direction. As such, the most plausible interpretation is that the spacecraft only skims the outer edge of the CME and therefore the higher speed component, which drives the shock, is not observed in situ. This interpretation is consistent with the steadily decreasing magnetic strength profile (as opposed to peaking later). As shown before, due to the magnetic configuration and the interaction with the HSS, the CME does not follow a radial direction but deviates at least 30 • from it. This could explain the relative short standoff distance (e.g., shock-sheath duration) as the shock could still be driven at the apex but the measured bulk speed at the intersection trajectory is equal to the speed of the HSS. In other words, the shock at the spacecraft crossing is being driven but not by the part of the CME at this trajectory but rather by the plasma streaming from the apex.
Based on these results, we suspect that the CME is engulfed by the HSS during most of its propagation to 1 AU, which is consistent with the results from of the flux rope evolution. The flux tube expansion factor is with values of measured n A = 0.51 −0.13 +0.14 at the lower end of the range obtained by statistical studies (see DeVore, 2000;Gulisano et al., 2012;Qiu et al., 2007;Temmer et al., 2017). Possible interpretations are that the FRs expansion was hindered and/or that only part of the magnetic structure is encountered by the spacecraft. We suspect both to have happened in this case.
The effects of this interaction on the Earth's magnetosphere are weak-to-moderate but show some surprising features such as a strong compression of the magnetorsphere and significant flaring during the CME sheath passage. We find some normal storm as well as substorm activity caused by the CIR. Due to the unusual high dynamic pressure of the shock, for a weak CME, the subsolar position of the magnetopause is decreased up to 27% during the shock arrival, which causes significant magnetopause currents. This leads to the storm activity to be considered moderate.

Summary and Conclusions
We studied the origin of a peculiar in situ signature measured at 1 AU, caused by an Earthdirected ICME and its interaction with a nearby HSS with the aim to address two major questions: How such a unique in situ signature could be formed and whether the observed shock was still driven. Using a combination of observations and modeling efforts we were able to create a consistent interpretation of the event answering these questions: • The CME is launched non-radially in south-eastern direction due to the configuration of the local magnetic field. This is reflected in the evolution of the dimming regions, the CME loops low in the corona (as seen in EUV observations) and the ForeCAT results. • At a height after 1.3 R the CME runs into the HSS which slows down the southward propagation until it is reversed and the CME is propagating northward. This is indicated by both the GCS measurements and the ForeCAT model. However, to be consistent with the GCS results, the ForeCAT model input requires an artificially scaled background field and a wider CME (overexpansion of CME cross section). With this the ForeCAT model mimics the effects of compression of the magnetic field due to the CME, which compresses the CH's open field and is then deflected away from it. • Due to deflection and the HSS "wrapping" around the CME, we measure a relatively small FR size and possible slow expansion (owing to the estimated low-value n A = 0.51 −0.13 +0.14 ). • The short standoff distance as well as the shock signature may be interpreted such as that the CME, which is engulfed by the HSS, only skims the spacecraft so its high speed part might have been missed. This interpretation is also consistent with the shape of the magnetic field profile which is steadily decreasing and shows no clear peak as expected from a central flux rope crossing. Due to the magnetic configuration and the interaction with the HSS, the CME does not follow a radial direction but deviates at least 30 • from it, causing this geometric effect in the measurements.
Acknowledgements Open access funding provided by University of Graz. The SDO and STEREO image data and the Wind, ACE and STEREO in situ data is available by courtesy of NASA and the respective science teams. We also thank the OMNI website at Goddard Space Flight Center, USA for a significant portion of the data used. S.G.H., M.T., K.D. and A.M.V. acknowledge funding by the Austrian Space Applications Programme of the Austrian Research Promotion Agency FFG (ASAP-13 859729 SWAMI, ASAP-11 4900217 CORDIM and ASAP-14 865972 SSCME). S.J.H. acknowledges support from the JungforscherInnenfonds der Steiermärkischen Sparkassen. C.J.F. acknowledges support from NASA Wind NNX16AO04G and STEREO Quadrature grant. N.L. and C.J.F. acknowledge support from NSF grant AGS-1435785. F.C. acknowledges the support by FEDER/Ministerio de Ciencia, Innovación y Universidades -Agencia Estatal de Investigación/Proyecto (ESP2017-88436-R) and Predoctoral Research Grants 2016-MCIU/FSE (BES-2016-077267). T.W. acknowledges DLR-grant 50 OC 1701 and DFG-grant WI3211/4-1. We thank the anonymous referee for constructive comments, which helped to improve the manuscript.

Disclosure of Potential Conflicts of Interest
The authors declare that they have no conflicts of interest.
Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided 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.