Evolution and Consequences of Interacting CMEs of 2012 November 9-10 using STEREO/SECCHI and In Situ Observations

Wageesh Mishra 1 · Nandita Srivastava 1 · D. Chakrabarty 2

c © Springer

••••

Abstract Understanding of the kinematic evolution of Coronal Mass Ejections (CMEs) in the heliosphere is important to estimate their arrival time at the Earth. It is found that kinematics of CMEs can change when they interact or collide with each other as they propagate in the heliosphere. In this paper, we analyze the collision and post-interaction characteristics of two Earth-directed CMEs, launched successively on 2012 November 9 and 10, using white light imaging observations from STEREO/SECCHI and in situ observations taken from WIND spacecraft. We tracked two density enhanced features associated with leading and trailing edge of November 9 CME and one density enhanced feature associated with leading edge of November 10 CME by constructing J-maps. We found that the leading edge of November 10 CME interacted with the trailing edge of November 9 CME. We also estimated the kinematics of these features of the CMEs and found a significant change in their dynamics after interaction. In in situ observations, we identified distinct structures associated with interacted CMEs and also noticed their heating and compression as signatures of CMECME interaction. Our analysis shows an improvement in arrival time prediction of CMEs using their post-collision dynamics than using pre-collision dynamics. Estimating the true masses and speeds of these colliding CMEs, we investigated the nature of observed collision which is found to be close to perfectly inelastic. The investigation also places in perspective the geomagnetic consequences of the two CMEs and their interaction in terms of occurrence of geomagnetic storm and triggering of magnetospheric substorms.

Keywords: Coronal mass ejections, STEREO, Heliospheric Imagers

1 Udaipur Solar Observatory, Physical Research Laboratory,

Udaipur, India -313001

email:

wageesh@prl.res.in email:

nandita@prl.res.in

Space and Atmospheric Sciences Division, Physical

2

Research Laboratory, Ahmedabad-380009, India email:

dipu@prl.res.in

## 1. Introduction

Coronal Mass Ejections (CMEs) are huge magnetized plasma expulsions from the Sun into the heliosphere and are considered as primary drivers of geomagnetic storms and various other space weather effects (Gosling et al. , 1990; Echer et al. , 2008). It is therefore required to improve the physics and our understanding of evolution of CMEs in the solar wind, which is found to be dependent on heliospheric environment (Kilpua et al. , 2012; Vrˇnak s et al. , 2013). Various techniques based on stereoscopic observations which estimate the 3D kinematics of CMEs in the heliosphere have been developed (Mierla et al. , 2008; Thernisien, Vourlidas, and Howard, 2009; Thompson, 2009; Liu et al. , 2010a; Lugaz et al. , 2010; Davies et al. , 2013). Using these techniques, association of remote observations of CMEs with in situ observations has been carried out by Liu et al. (2010b); Temmer et al. (2011); Kilpua et al. (2012); M¨stl o et al. (2012); Mishra and Srivastava (2013); Mishra et al. (2014). Also, several models have been developed for predicting the arrival time of CMEs at 1 AU (Gopalswamy et al. , 2000; Vrˇ snak et al. , 2013). Sometimes, CMEs are observed to erupt in quick succession which under certain favorable initial conditions, can interact or merge with each other during their propagation in the heliosphere. In such a scenario, space weather prediction schemes may not be successful without taking into account the post-interaction CMEs characteristics. The compound streams (interaction of CME-CIR or CME-CME) were first inferred by Burlaga, Behannon, and Klein (1987) using observations from Helios and ISEEE 3 spacecraft. Their study suggested that a large geomagnetic storm may be associated with compound streams. The observational evidence for merging of CMEs (or cannibalism) was first given by Gopalswamy et al. (2001) based on the longwavelength radio and SOHO/LASCO observations. Thus, these interactions are of great importance from space weather point of view. Before the era of wide angle imaging far from the Sun, the understanding of involved physical mechanisms in CME-CME or CME-shock interaction was achieved mostly from magnetohydrodynamic (MHD) numerical simulations of the interaction of a shock wave with an magnetic cloud (MC) (Vandas et al. , 1997; Vandas and Odstrcil, 2004; Xiong et al. , 2006), the interaction of two ejecta (Gonzalez-Esparza, Santill´n, a and Ferrer, 2004; Lugaz, Manchester, and Gombosi, 2005; Wang et al. , 2005), and the interaction of two MCs (Xiong et al. , 2007; Xiong, Zheng, and Wang, 2009). Wang, Ye, and Wang (2003) have shown that a forward shock can cause an intense southward magnetic field of long duration in the preceding magnetic cloud. Such modifications in the preceding cloud are important for space weather prediction. Therefore, events involving CME-CME and CME-shock interactions are important candidates to investigate their kinematics and arrival time at near Earth.

After the launch of Solar TErrestrial RElation Observatory (STEREO) in 2006 (Kaiser et al. , 2008), it is possible to continuously image the CMEs and to determine the 3D locations of their features in the heliosphere and hence have direct evidence of CME interaction using images from Heliospheric Imager (HI) on Sun Earth Connection Coronal and Heliospheric Investigation (SECCHI) package (Howard et al. , 2008). However, immediately after the launch of

STEREO, during deep extended solar minimum, not many interacting CMEs were observed. As the current solar cycle progressed, interaction of CMEs now appears a fairly common phenomenon, in particular around solar maximum. The interacting CMEs of 2010 August 1 have been studied extensively by using primarily the STEREO white light imaging, near Earth in situ and, S-Waves radio observations (Harrison et al. , 2012; Liu et al. , 2012; M¨ ostl et al. , 2012; Temmer et al. , 2012; Mart´ ınez Oliveros et al. , 2012; Webb et al. , 2013). Also, Lugaz et al. (2012) have reported a clear deflection of 2010 May 23 - 24 CMEs after their interaction. As mentioned above, therefore, only a few studies on CME interaction have been reported. Several key questions regarding CME interaction need to be addressed, which are not well understood quantitatively,viz. (1) How the dynamics of CMEs change after interaction? What is the regime of interaction, i.e. elastic, inelastic or super-elastic? (Lugaz et al. , 2012; Shen et al. , 2012; Mishra and Srivastava, 2014) (2) What are the consequences of the interaction of CME-shock structure? How the overtaking shock changes the plasma and the magnetic field properties into preceding magnetic cloud? (Lugaz, Manchester, and Gombosi, 2005; Lugaz et al. , 2012; Liu et al. , 2012) (3) What are the favorable conditions for the CME cannibalism and the role of magnetic reconnection in it? (Gopalswamy et al. , 2001) (4) What is the possibility for production of reverse shock at CME-CME interaction site? (Lugaz, Manchester, and Gombosi, 2005) (5) Whether these interacted structures produce different geomagnetic consequences than individual CMEs, on their arrival to magnetosphere? (Farrugia et al. , 2006). (6) What are the favorable conditions for the CME's deflection and enhanced radio emission during CME-CME interaction? (Lugaz et al. , 2012; Mart´ ınez Oliveros et al. , 2012). In light of aforementioned questions and only a few studies reported, it is clear that the prediction of arrival time of interacting CMEs and association of remote observations of such CMEs with in situ is challenging.

In the present study, we focus on the identification, evolution and propagation of two CMEs launched on 2010 November 9 and 10, as they travel from the corona to the inner heliosphere. This study, in which the launch time of the selected CMEs are separated by about 14 hr, gives an opportunity to understand the CME-CME interaction unambiguously for simpler scenario, contrary to complex events around 2010 August 1, wherein four CMEs were launched in quick succession. We find that CMEs of November 9 and 10 are directed towards the Earth and their initial characteristics indicate a probability of their interaction. In Section 2.1, observations and estimated kinematics of these CMEs are presented to investigate if they interacted during their propagation in the heliosphere. In Section 2.2, nature of collision, momentum and energy transfer during collision phase of CMEs are described using their estimated true mass and kinematics. In Section 3.1, in situ observations of remotely tracked CME features are described. Using the derived kinematics combined with Drag Based Model (DBM) of Vrˇ snak et al. (2013), the arrival time and transit speed of tracked features of CMEs at L1 is predicted in Section 3.2. The geomagnetic response of interacted CMEs structures is described in Section 3.3. In Section 4, the results and discussions of the present study are emphasized and conclusions are given in Section 5.

## 2. Observations

For the present study of 2012 November 9-10, CME-CME interaction, we analyzed the white light images from SECCHI suite on-board NASA's twin STEREO (A &amp; B) mission. SECCHI suite consists of an Extreme Ultraviolet Imager (EUVI), two coronagraphs (COR1 and COR2) and, two Heliospheric Imagers (HI1 and HI2). These instruments can image the heliosphere from near the Sun to the Earth and beyond. The EUVI (Field-of-view (FOV): 1-1.7 R glyph[circledot] ), COR1 (FOV: 1.5-4 R glyph[circledot] ) and COR2 (FOV: 2.5-15 R glyph[circledot] ) FOV are centered on the Sun. The HI1 and HI2 have their optical axis aligned in the ecliptic plane and are off-centered from the Sun at solar elongation of 14 ◦ and 53 7 , respectively. . ◦ HI1 and HI2 have wide FOV of 20 ◦ and 70 , respectively. With this identical ◦ imaging package (SECCHI) on twin spacecraft, a CME can be tracked from 0.4 to 88 7 . ◦ elongations. At the time of observations around November 9, the STEREO A &amp; B were 127 ◦ westward and 123 ◦ eastward from the Sun-Earth line at a distance of 0.96 AU and 1.08 AU from the Sun, respectively. The in situ properties were investigated through data collected by instruments on-board Wind spacecraft (Lepping et al. , 1995; Ogilvie et al. , 1995).

## 2.1. Interaction of 2012 November 9-10 CMEs

In order to determine the true direction and speed of CMEs launched on November 9 and 10 in quick succession, 3D reconstruction technique is applied to CME features close to the Sun. We find that the CME of November 9 (hereinafter referred as, CME1) has a slow speed than fast speed November 10 CME (hereinafter referred as, CME2) and both move in a similar direction towards the Earth. These observations indicate the probability of their interaction in the heliosphere. Using SECCHI/HI observations, we attempt to investigate if these CMEs indeed interacted with each other in the heliosphere. We also attempt to identify various CME structures in the in situ observations taken at 1 AU.

## 2.1.1. 3D reconstruction in COR2 Field-of-view

Apartial halo CME1 (angular width of 276 ◦ ) was observed in SOHO/LASCO-C2 FOV at 15:12 UT on 2012 November 9 ( http://cdaw.gsfc.nasa.gov/CME list/ ). In the SECCHI/COR1-B and A, this CME was observed at 14:25 UT in the SW and SE quadrants, respectively. Another partial halo CME2 (angular width ∼ = 210 ) was observed in LASCO C2 FOV at 05:12 UT on November 10. This ◦ CME was detected by SECCHI/COR1-B and A in the SW and SE quadrants, respectively, at 05:05 UT on November 10, 2012.

Figure 1 shows true kinematics of CME1 and CME2 from 3D reconstruction of selected features along the leading edges of corresponding CMEs using tiepointing procedure (scc measure: Thompson 2009) on SECCHI/COR2 data. This tie-pointing procedure uses the concept of epipolar geometry (Inhester, 2006) and due to epipolar constraints, the identification of same feature in two images (from STEREO-A and B) is always possible along the same epipolar line in both images. Before implementing this technique, the processing of

SECCHI/COR2 images and creation of minimum intensity images and then its subtraction from the sequence of processed COR2 images was carried out following the approach described in Mierla et al. (2008); Srivastava et al. (2009). As the separation angle of the STEREO spacecraft is large, it is better to use two independent methods to confirm the results of 3D reconstruction. Therefore, we have also used the visual fitting of both CMEs in the SECCHI/COR2 FOV using the Graduated Cylindrical Shell (GCS) model (Thernisien, Vourlidas, and Howard, 2009; Thernisien, 2011) for the 3D reconstruction of these CMEs. This model represents the flux rope structure of CMEs with two shapes: the conical legs and the curved (tubular) fronts, also known as the 'hollow croissant'. We used contemporaneous image triplets of CME1 around 17:39 UT on November 9 from STEREO-A/COR2, STEREO-B/COR2 and SOHO/LASCO-C3. The best visual fit for CME1 is found in the direction of W02S14, with a half angle of 19 , tilt angle of 9 ◦ ◦ around the axis of symmetry of the model (i.e. rotated 9 ◦ anticlockwise out of the ecliptic plane), and an aspect ratio of 0.52. At this time, CME1 was at a distance of about 9.6 R glyph[circledot] from the Sun. Using the obtained fitted values of half angle and aspect ratio of CME1, its 3D face-on and edge-on angular width is estimated as 100 ◦ and 62 , respectively. We also carried out ◦ the visual fitting using GCS model for CME2 around 06:39 UT on November 10 by exploiting the concurrent image triplets of STEREO-A/COR2, STEREOB/COR2 and SOHO/LASCO-C3. The best fit for CME2 is obtained in the direction of W06S25, at a tilt angle of 9 ◦ around the axis of symmetry of the model, with a half angle of 12 ◦ , and an aspect ratio of 0.19. At this time CME2 was at a distance of 8.2 R glyph[circledot] from the Sun. Using the GCS model fitted values, the 3D face-on and edge-on angular width of CME2 is estimated as 46 ◦ and 22 , respectively. The GCS model fit for CME1 and CME2 is shown in the top ◦ and bottom panel of Figure 2, respectively. We find that the estimated latitude and longitude of both CMEs from GCS model and tie-pointing procedure are in agreement within few degrees. This shows that the results of both tie-pointing and GCS model are reliable and can be used for estimating kinematics of CMEs in the coronagraphic FOV.

From a study of 3D kinematics of the two CMEs using tie-pointing approach (scc measure: Thompson 2009), we found that CME1 is slow with a true speed of 620 km s -1 around 15 R glyph[circledot] while CME2 is faster with a speed of 910 km s -1 at a distance of approximately 15 R glyph[circledot] (Figure 1). Latitude and longitude of reconstructed features for both CMEs suggest that these were Earth-directed and could possibly interact in the interplanetary (IP) medium. If the speed of CMEs is assumed to be constant beyond the outer edge of COR2 FOV, then these CMEs should have collided at approximately 130 R glyph[circledot] on November 11 around 02:30 UT. However, previous studies have shown that the speed of a CMEcan significantly change after the coronagraphic FOV (Lindsay et al. , 1999; Gopalswamy et al. , 2000; Cargill, 2004; Manoharan, 2006; Vrˇ snak et al. , 2010). Therefore, further tracking of CME features in the heliosphere is required to determine the exact location and time of interaction of these CMEs.

## 2.1.2. Tracking of CMEs in HI Field-of-view

Evolution of CME1 and CME2 in the running difference images of COR2, HI1 and HI2 FOV is shown in Figure 3. To track and estimate the arrival times of the CMEs in the heliosphere using HI1 &amp; 2 images, we constructed J-maps based on the method developed by Sheeley et al. (1999); Davies et al. (2009) and derived the variation of elongation of selected features with time. The details about extraction of strip of constant position angle, selection of bin size and adopted procedure for construction of J-maps and derivation of elongation angle from it has been described in section 3.1.1 of Mishra and Srivastava (2013). The positively inclined bright features in the J-maps (Figure 4) correspond to enhanced density structure of the CMEs. We tracked the leading and trailing edge (marked in green and red) of bright features corresponding to slow CME1 and leading edge (blue) of bright feature corresponding to the fast CME2. Prior to interpreting the tracked features of CME1 in J-map as its leading and trailing edge, respectively, the derived elongations corresponding to all the three tracked features were overplotted on the base difference HI1 images. On careful inspection of the sequence of these images, we noticed that the two tracked features of CME1 correspond to first density enhancement at the front and second density enhancement at the rear edge of CME1, respectively. The base difference HI1-A images with overplotted contours of elongation corresponding to tracked features of CME1 and CME2, are shown in Figure 5. It is to be noted that the term 'trailing edge' of CME1 does not correspond to the rear most portion of CME1 but corresponds to the part behind the front and is not an image artifact in J-map constructed from running difference images.

Since, the tracked features of CME2 were not observed well in STEREO B ecliptic J-map, we could not implement the stereoscopic reconstruction technique to estimate the CME kinematics. Instead, we used Harmonic Mean (HM) approximation which requires single spacecraft observation to derive the elongation variations with time for the tracked features which could then be converted into heliocentric distances from the Sun (Lugaz, Vourlidas, and Roussev, 2009). In the HMapproximation, geometry of a CME is considered as a self-similar expanding sphere attached to the Sun and its centre moves along a fixed longitude. The choice of the HM method is based on our earlier study on comparison of various reconstruction techniques involving both stereoscopic and single spacecraft observations (Mishra et al. , 2014). It was found that amongst the single spacecraft techniques, the HM method is the most suitable for arrival time prediction. Implementing the 3D reconstruction technique in Section 2.1.1 (Figure 1), we have shown that the estimated longitude of CME1 and CME2 are approximately 10 ◦ and 2 ◦ east of the Sun-Earth line. We assume that beyond the COR2 FOV, CMEs continue to propagate along the same direction i.e. we ignore any heliospheric longitudinal deflection of CMEs. Therefore, the aforesaid estimated values of longitude were used in the HM approximation to convert the derived elongation from J-maps to radial distance from the Sun.

Figure 6 shows the distance and speed estimated for different tracked features of the two CMEs. The estimated speed is derived from the adjacent distances using a numerical differentiation with three point Lagrange interpolation and

therefore has systematic variations due to slight errors in the estimated distance. It is clear that the leading edge (LE) of CME1 has higher speed ( ≈ 500 km s -1 ) than its trailing edge (TE) speed ( ≈ 350 km s -1 ), averaged over few data points at the entrance of HI1 FOV. Also, the LE and TE of CME1 have lower speeds than LE of CME2. LE of CME2 shows a large radial speed of approximately 1100 km s -1 (ecliptic speed = 950 km s -1 ) in the COR2 FOV i.e.(2.5-15 R glyph[circledot] ) (Figure 1). Beyond 10 R glyph[circledot] distance, the LE of CME2 continuously decelerates for ≈ 10 hours up to 46 R glyph[circledot] where its speed reduces to 430 km s -1 . The extreme fast deceleration of LE of the CME2, starting from COR2 FOV, seems to be due to possible interaction with the preceding CME1. It is likely that the CME1 had large spatial scale due to which the trailing plasma and magnetic fields from CME1 created sufficiently dense ambient medium acting as a huge drag force for CME2 which resulted in its observed deceleration. The extreme fast deceleration of CME2 LE can also be due to closed magnetic structure of CME1 which may act like a magnetic obstacle for the CME2 (Temmer et al. , 2012).

From estimated kinematics of tracked features in heliosphere (top and bottom panel of Figure 6), it is clear that around November 10 at 11:30 UT speed of CME1 TE started to increase from 365 km s -1 with simultaneous decrease of CME2 LE speed from 625 km s -1 . Such an observation of acceleration of CME1 and deceleration of CME2 provides an evidence for the commencement of collision (Temmer et al. , 2012; Lugaz et al. , 2012; Shen et al. , 2012; Mariˇi´ c c et al. , 2014; Temmer et al. , 2014; Mishra and Srivastava, 2014). We define collision phase, as the interval during which the two CMEs come in close contact with each other and show opposite trend of acceleration relative to one another until they attain an approximately equal speed or their trend of acceleration is reversed. Carefully observing the variations of speed of tracked CME1 TE and CME2 LE, the start and end boundary of collision phase is drawn as vertical lines in the bottom panel of Figure 6. We noticed that at the end of the collision phase around November 10 at 17:15 UT, speed of CME1 TE is ≈ 450 km s -1 and speed of CME2 LE is 430 km s -1 . In the beginning of the marked collision phase, CME1 TE is at a distance of 37 R glyph[circledot] and CME2 is at 30 R glyph[circledot] . At the end of the marked collision phase, the CME1 TE and CME2 LE are at a distance of 50 R glyph[circledot] and 46 R glyph[circledot] , respectively. Here, we must highlight that throughout the collision phase, the estimated heliocentric distance of CME1 TE and CME2 LE is never found to be equal. This is because of the large scale structure of CMEs and we suggest that the tracked TE of CME1 and LE of CME2 using J-maps, are not strictly the rearmost trail of CME1 and outermost front of CME2, respectively. Therefore, we had defined 'collision' when tracked features of CMEs (i.e. time-elongation tracks) come in close contact (separation of ≈ 5 R ) with one another. Such small separation between observed collided features glyph[circledot] is possible if forward shock driven by CME2 interacts with the rear most portion of the CME1 as also suggested by Temmer et al. (2012); Mariˇi´ c c et al. (2014).

## 2.2. Momentum, Energy Exchange and Nature of Collision

To understand the momentum exchange during collision of CMEs, we require their true masses. Since, appearance of a CME is due to Thomson scattered

photospheric light from the electrons in the CME (Minnaert, 1930; Billings, 1966; Howard and Tappin, 2009), therefore, the recorded scattered intensity can be converted into the number of electrons and hence mass of a CME can be estimated, if composition of CME is known. But an observer from different vantage points receives different amount of scattered light by the electrons, therefore, the true location of electrons of CME, i.e. propagation direction of CME must be known to estimate the true mass of CME. For a long time, mass of a CME was calculated assuming the CME location in the observer's plane of sky (Munro et al. , 1979; Poland et al. , 1981; Vourlidas et al. , 2000). Although, the propagation direction of CME is calculated in our study using tie-pointing (Thompson, 2009) and Forward modeling (Thernisien, Vourlidas, and Howard, 2009) methods and described in Section 2.1.1, but to avoid any bias of the reconstruction methods, we use the method of Colaninno and Vourlidas (2009) which is based on Thomson scattering theory, to estimate the true propagation direction and true mass of both the CMEs. We used the simultaneous image pair of SECCHI/COR2 to estimate the true mass of CMEs. First, we estimated the mass of CME i.e. M A and M B , from STEREO-A and STEREO-B locations. Then we took the ratio of equations (7) and (8) of Colaninno and Vourlidas (2009) and obtained a new expression as M /M A B = B e ( θ A ) /B e ( θ A +∆), where θ A is the angle of propagation direction of CME measured from the plane-ofsky (POS) of STEREO-A , B ( e θ A ) is the brightness of a single electron at an angular distance of θ A from the POS and ∆ is the separation angle between both STEREO-A and STEREO-B from the Sun-Earth line. Using the measured M A and M B in aforementioned expression, the value of θ A can be estimated. However, we noticed that multiple values of θ A correspond to same ratio of M A and M B but the correct value of θ A among them can be selected by visual inspection of CME in COR2 images. We obtained the true mass of CME using the derived value of θ A in equation (4) of Colaninno and Vourlidas (2009). The estimated mass of CME1 at 18:39 UT on November 9 at ≈ 15 R glyph[circledot] , from STEREO-A and STEREO-B vantage point is 4.60 × 10 12 kg and 2.81 × 10 12 kg, respectively. For CME2, the estimated mass at 07:24 UT on November 10 at ≈ 15 R glyph[circledot] , from STEREO-A and STEREO-B locations is 2.25 × 10 12 kg and 1.31 × 10 12 kg, respectively. The estimated propagation direction of CME1 and CME2 is 19 ◦ and 21 ◦ West from the Sun-Earth line, respectively. We estimated the true mass of CME1 and CME2 as M 1 = 4.66 × 10 12 kg and M 2 = 2.27 × 10 12 kg, respectively.

We assume that after crossing the COR2 FOV and during collision of CME1 and CME2, their estimated true masses (M 1 and M ) remain constant as es2 timated earlier. From the Figure 6, the observed velocity of CME1 and CME2 before the collision is estimated as (u ,u ) = (365,625) km s 1 2 -1 and observed velocity of CME1 and CME2 after the collision is (v ,v ) = (450,430) km s 1 2 -1 . To understand the nature of collision of large scale magnetized plasmoids, we attempt to estimate the coefficient of restitution ( e ) of colliding CME1 and CME2. The coefficient of restitution measures the bounciness (efficiency to rebound) of a pair of objects in collision and is defined as ratio of their relative velocity of separation to relative velocity of approach. Hence, for e &lt; 1, e = 1, and e &gt; 1 the collision is termed as inelastic, elastic and super-elastic and consequently

the kinetic energy of the system after the collision is found to decrease, equal and increase than before the collision, respectively. We restrict ourselves not to estimate the value of e directly by using the pre and post collision velocity of CMEs. This is because the velocities estimated from reconstruction method have some errors and do not guarantee conservation of momentum, a necessary condition of collision, therefore, can lead to erroneous estimation of e value. In our approach, considering the condition of conservation of total momentum for the collision, the velocity of CME1 and CME2 after the collision can be estimated theoretically (v 1 th ,v 2 th ), provided that velocity of CME1 and CME2 before the collision and the coefficient of restitution ( e ) value are known.

Using the velocities (u 1 ,u 2 ) and masses (M ,M ), we iterate for a range of 1 2 e values and estimate the (v 1 th ,v 2 th ) each time. Using this approach, the most suitable value of e , using which the theoretically estimated (v 1 th ,v 2 th ) is found to be closest to the observed (v 1 ,v 2 ) i.e. for which the variance,

σ = √ ( v 1 th -v 1 ) 2 +( v 2 th -v 2 ) 2 , is minimum, can be obtained. The details of the above approach and the equations used, are explained in Mishra and Srivastava (2014). Using the aforementioned approach, the value of e = 0.1 is found with (v 1 th ,v 2 th ) = (458,432) km s -1 and σ = 9. With this theoretically estimated e value, total kinetic energy of CMEs is found to decrease by 6.7% of its value before the collision. Using the estimated velocities of the CMEs before and after the collision (Figure 6), the value for e is calculated equal to 0.08, which is approximately same as obtained from iterations described above. Therefore, our analysis suggests that the observed collision between the CMEs are close to perfectly inelastic in nature. The kinetic energy of CME1 and CME2 before the collision was 3.1 × 10 23 joules and 4.4 × 10 23 joules, respectively. After the collision, based on the observed speeds, we found that the kinetic energy of CME1 increased by 51% to its value before the collision while the kinetic energy of CME2 decreased by 54.5% to its value before the collision. We also noticed that after the collision, momentum of CME1 increased by 23% and momentum of CME2 decreased by 31% to their values before the collision. Such calculations support the claim that significant exchange of kinetic energy and momentum take place during CME-CME collision (Temmer et al. , 2012; Lugaz et al. , 2012; Mariˇ ci´ c et al. , 2014; Shen et al. , 2012).

Considering the uncertainties of ± 100 km s -1 in the speed, we repeated the above analysis. Taking the (v ,v ) = (550,530), we found 1 2 e = 0, i.e. perfectly inelastic. But in this case, (v 1 th ,v 2 th ) is estimated as (450,450) km s -1 which highlight that theoretically estimated velocity is not close to the observed velocity (v 1 ,v 2 ). If we consider the (v ,v ) = (350,330), then it implies that the 1 2 velocities of both CMEs after the collision are less than their values before the collision. Therefore, such velocities also cannot be considered for our case (as defined in Section 2.1.2). As, there are many sources of errors in the estimation of true mass of CMEs, therefore, the effect of uncertainty in mass must also be examined before confirming the nature of observed collision. In our case, the mass ratio (M /M ) is equal to 2.05, varying this ratio between 0.5 to 3.0 in 1 2 a step of 0.25, and repeating the iterative procedure for estimating the most suitable value of e corresponding to each mass ratio, we obtain e close to 0. Figure 7 shows the variation of e and σ with varying mass ratio. It is clear that

despite taking large uncertainties in the mass of CME1 and CME2, the nature of observed collision remains close to perfectly inelastic for the case of interacting CMEs of 2012 November 9-10.

## 3. In Situ observations, Arrival time and Geomagnetic Response of 2012 November 9-10 CMEs

## 3.1. In Situ Identification of Tracked CME Features

We analyzed the in situ data taken from WIND spacecraft located at L1 to identify the tracked density enhanced features of CMEs. Figure 8 shows magnetic field and plasma measurements during November 12, 12:00 UT to November 15, 12:00 UT. The arrival of a forward shock (labeled as S) marked by a sudden enhancement in speed, temperature and density is noticed at 22:20 UT on November 12. The region between the first and second vertical line represents the turbulent sheath region. Based on the CME identification criteria of Zurbuchen and Richardson (2006), the region bounded between second vertical line at 08:52 UT on November 13 and third vertical line at 02:25 UT on November 14, is identified as CME structure. In this region, we observed expansion of CME characterized by a monotonic decrease in proton speed and temperature. Based on the predicted arrival times which is derived using the estimated kinematics of the remotely observed tracked features of CME1 and CME2 as inputs in the DBM (explained in Section 3.2), this region is associated with Earth-directed CME1 launched on November 9. During the passage of CME1, the magnetic field is observed to be high ( ≈ 20 nT), plasma beta is less than unity ( β &lt; 1) with smooth rotation in magnetic field vector. Also, the latitude ( θ ) value of magnetic field vector decreased from 43 ◦ to -43 ◦ and its longitude ( φ ) decreased from 203 ◦ to 74 ◦ . Therefore, this region can be classified as a magnetic cloud (MC) and based on the observed arrival time it is associated with CME1. Due to the interaction of the CMEs, the region associated with CME1 is found to be at higher temperature than found generally for a normal isolated CME. After the third vertical line, magnetic field strength decreased reaching a minimum value of 6 nT around 4:00 UT on November 14. This interval of sudden drop in magnetic field is associated with a sudden rise in density, temperature and plasma β , suggestive of a possible magnetic hole (MH) (Tsurutani et al. , 2006) which is considered as a signature of magnetic reconnection (Burlaga and Lemaire, 1978). Another region of magnetic field depression from 08:05 UT to 10:15 UT on November 14, reaching a minimum value of 3 nT is also noticed. Corresponding to this minimum value of magnetic field, the plasma β and temperature is found to increase. The region during 03:45 UT - 08:05 UT on November 14, bounded between two distinct MH like structures, has enhanced magnetic field ( ≈ 15 nT), and plasma β less than unity. This region, between two MH, seems to be a magnetic field remnant of reconnecting CMEs structures. Based on the extremely elevated interval of plasma beta, temperature, observation of magnetic hole (MH) and probably sudden fast rotation in magnetic field vector during the region bounded between third and fourth vertical line, in Figure 8, this region

is identified as interaction region (labeled as IR) of CME1 TE with CME2 LE. Another structure is identified based on the elevated fluctuating magnetic field and temperature during 12:00 UT - 21:21 UT on November 14, bounded between fourth and fifth vertical line in Figure 8. During this interval, we noticed high magnetic field (9 nT) with no monotonic decrease in temperature and speed profile as well as plasma beta ( β ) is not less than unity. From these observations, i.e. lack of MC signatures and the short duration (9.5 hr) of the structure associated with CME2, we infer that the WIND spacecraft perhaps intersected the flank of CME2 (M¨ ostl et al. , 2010). This is also confirmed by the estimated latitude (-25 ◦ ) of CME2 using GCS model (described in Section 2.1.1). On examining the evolution of CME2 in HI1-A movie, we noticed that CME2 is directed towards the southern hemisphere which is a favorable condition for its flank encounter with in situ WIND spacecraft.

We observed that the magnitude of magnetic field in MC region is constant around ≈ 20 nT which may be due to the passage of shock from CME1. We also noticed that average temperature in first half of CME1 is high ( ≈ 10 5 K) than found in general ( ≈ 10 4 K). The high temperature of the CME1 may occur due to its collision with CME2, thereby resulting in its compression. Another possibility of high temperature of CME1 is due to the passage of forward shock driven by CME2 as reported in earlier studies by Lugaz, Manchester, and Gombosi (2005); Liu et al. (2012); Mariˇi´ c c et al. (2014); Mishra and Srivastava (2014). However, in the present case, the following CME2 is also observed with a high temperature ( ≈ 5 × 10 5 K), which has not been reported in earlier studies of interacting CMEs. We noticed significant high density at the front of CME1, possibly due to overall compression by sweeping of the plasma of CME1 at its leading edge by CME2 driven shock.

## 3.2. Arrival time of Tracked Features

We used the estimated speed at the last point of measurement (up to where CMEs could be tracked unambiguously in HI) and used it as input in the DBM developed by Vrˇ snak et al. (2013) to estimate the arrival time of tracked features at L1. Vrˇnak s et al. (2013) have shown that drag parameter lies in the range of (0.2 - 2.0) × 10 -7 km -1 . Using these values of the drag parameter in the DBM, the arrival time of CME at L1 can be predicted with a reasonably good accuracy (Vrˇ snak et al. , 2013; Mishra and Srivastava, 2013). The estimated speed, time and distance (v , t 0 0 and R glyph[circledot] ) of LE of CME1 (green track in J-map) are used as inputs in the DBM, corresponding to extreme range of drag parameter, to predict its arrival time and transit speed at L1. In situ observations show a peak in density, ≈ 0.5 hr after the shock arrival around 23:00 UT on November 12 with a transit speed of 375 km s -1 , which is expected to be the actual arrival of tracked LE feature corresponding to CME1. The predicted values of arrival time and transit speed of features and errors from the actual values are shown in Table 1.

For TE feature of CME1 and LE of CME2, we assume that they encounter the dense ambient solar wind medium created by the preceding CME1 LE. Therefore, their kinematics with the maximum value of the statistical range of the drag

parameter (2.0 × 10 -7 km -1 ), are used as inputs in the DBM. We consider that the TE of CME1 corresponds to the density enhancement at trailing front of CME1. At the rear edge of CME1, a density enhancement (12 particles cm -3 ) is observed around 23:30 UT in WIND observations on November 13, which is considered as the actual arrival of TE of CME1. The predicted arrival time and transit speed of this feature and errors therein, are given in Table 1. We further noticed that in situ data, an enhancement in density corresponding to LE of CME2 around 12:00 UT on November 14 (marked as arrival of CME2 with fourth vertical line in Figure 8) is observed which can be considered as the actual arrival time of LE of CME2. The predicted values of CME2 LE, using its kinematics with DBM, and errors therein are also listed in Table 1. The error for CME2 LE is large, but, several factors can lead to such large errors which are discussed in Section 4. It must be highlighted that if the 3D speed of CMEs estimated at final height in COR2 FOV is assumed as constant up to L1, then the predicted arrival time of CME1 and CME2 will be ≈ 10-16 hr and 44 hr earlier, respectively, than their predicted arrival times using post-collision speeds combined with DBM. This emphasizes the use of HI observations and post-collision speeds of CMEs as inputs in the DBM for improved arrival time prediction of interacting CMEs. Similarly, using HI observations, Colaninno, Vourlidas, and Wu (2013) have shown that linear fit of deprojected height-time above 50 R glyph[circledot] gives a half a day improvement over CME arrival time estimated using LASCO data.

## 3.3. Geomagnetic consequences of 2012 November 9-10 CMEs

As mentioned in the introduction section, very few studies have been dedicated to the study of the geomagnetic consequences of interacting CMEs in the past. The CMEs of 2012 November 9-10 resulted in a single strong geomagnetic storm with Dst index ≈ -108 nT at 8:00 UT on November 14, therefore it is important to investigate the impact of the interaction of the CME1 of November 9 with the CME2 of November 10 on the terrestrial magnetosphere-ionosphere system (MI) in details.

Figure 9(a-e) reveal the variations in the solar wind parameters in the Geocentric Solar Ecliptic (GSE) coordinate system during 2012 November 12-15. These parameters include solar wind proton density (in cc -1 ), velocity (in km s -1 , negative X-direction), ram pressure (in nPa), Z-component of the interplanetary magnetic field (IMF B z , in nT) and Y-component of the interplanetary electric field ( IEF y , in mV/m) respectively. These data with cadence of 1 minute are taken from the NASA GSFC CDAWeb ( www.cdaweb.gsfc.nasa. gov/istp public/ ). It is also important to note that the solar wind parameters presented in Figure 9(a-e) are corrected for propagation lag till the nose of the terrestrial bow shock. In order to compare the variations of these parameters with the magnetospheric and ionospheric parameters, additional time lags that account for the magnetosheath transit time and the Alfven transit time are calculated (Chakrabarty et al. , 2005). Therefore, the solar wind parameters presented in Figure 9(a-e) are corrected for the propagation lag, point by point, till ionosphere. Figure 9(f) represents the variation in the polar cap (PC) index.

The PC index is shown (Troshichev et al. , 2000) to capture the variations in the ionospheric electric field over polar region efficiently. Figure 9(g) shows the variation in the Sym-H (in nT) index which primarily represents the variation in the magnetospheric ring current (Iyemori and Rao, 1996). The Sym-H index is essentially the high temporal resolution (1 min) version of the Dst index. Figure 9(h) shows variations in westward (midnight sector) auroral electrojet (AL) current (in nT) which captures auroral substorm processes reasonably well. Further, in Figure 9, the arrival of the shock is denoted by S, different vertical lines mark the arrival of different features of CME1 and CME2 and are labeled in the same manner as in Figure 8.

Figure 9(a) reveals the arrivals of two distinctly enhanced density structures, first enhancement occurred during arrival of the shock-sheath region before CME1 LE (on November 12 20:00 UT to November 13 04:00 UT) and second enhancement occurred during arrival of the CME1 TE-IR region (November 13 12:00 UT to November 14 12:00 UT). During the arrival of shock-sheath region before CME1 LE, the peak density reached ≈ 68 cc -1 which is more than two times of the corresponding peak density ( ≈ 30 cc -1 ) observed during the arrival of CME1 TE-IR. Sharp enhancement in the solar wind velocity was also observed vis-a-vis the sharp density enhancement during arrival of shock-sheath region (Figure 9b) when the velocity reached from ≈ 300 km s -1 to ≈ 470 km s -1 . However, solar wind velocity did not change sharply and significantly during arrival of CME1 TE-IR. The changes in the density and the velocity resulted in changes in the solar wind ram pressure shown in Figure 9(c). The peak ram pressure during shock-sheath region before CME LE ( ≈ 18 nPa) was almost twice than that observed during CME1 TE-IR ( ≈ 9 nPa). Figure 9(d) reveals that IMF B z was predominantly southward in both density enhanced intervals although fast fluctuations were observed corresponding to the arrival of the shock-sheath region. The first peak in IMF B z is ≈ -8 nT before the arrival of the shock. Thereafter, IMF B z fluctuated sharply between ≈ -20 nT to ≈ +20 nT during arrival of shock-sheath region before arrival of CME1 LE. During CME1 TEIR, IMF B z reached a peak value of ≈ -18 nT. No significant change in IMF B z is observed during the passage of the CME2 on November 14-15, when its magnitude hovered around the zero line.

Figure 9(e) elicits the variation in IEF y during 2012 November 12-15. It is observed that the peak value of IEF y before the arrival of the shock was ≈ 2 mV/m. However, during the shock event, IEF y fluctuated between ±∼ 8 mV/m. During arrival of CME1 TE-IR, the peak value of IEF y reached ≈ 7.5 mV/m. In fact, similar to IMF B z , a sharp polarity change in IEF y was noticed at ∼ 09:40 UT on November 14. No significant change in IEF y is observed during the passage of the CME2 on November 14-15, when its magnitude hovered around the zero line. Figure 9(f) shows that the PC index increased during both density enhanced intervals and the peak value is ≈ 5 for both the intervals. Figure 9(g) shows variation in the Sym-H index which revealed the development of a geomagnetic storm during CME1 TE region. The Sym-H index reached a value of ≈ -115 nT during the passage of IR on November 14. The ring current activity was not significant during arrival of the shock-sheath region before CME1. Lastly,

Figure 9(h) elicits the variation in the AL index. It is seen that AL reached ≈ -600 nT during the arrival of shock-sheath before CME1 and ≈ -1400 nT during the arrival of CME1 TE-IR. Therefore, significant intensification of westward auroral electrojet occurred during the arrival of CME1 TE-IR.

The above observations (from Figure 9) reveal several interesting points. First, the magnitude of AL seems to remain unaffected by large amplitudes of fluctuations in IMF B z and IEF y during the shock-sheath region. Second and the most important point is that the duration of occurrence of the AL intensification during shock-sheath before CME1 LE and in CME1 TE-IR region were nearly identical with the duration of the southward IMF B z and positive IEF y phases. However, the AL amplitudes seem to be more closely correlated with the IEF y amplitudes during CME1 TE-IR region compared to shock-sheath region. Third, although the peak amplitudes of PC index were nearly same during shock-sheath region before CME1 LE and CME1 TE-IR region, the peak amplitudes of the AL index were significantly different during these two intervals. Fourth, the substorm activity seems to be over during the passage of the CME2, therefore, CME2 did not have any bearing on the triggering of substorms. Figure 9 demonstrates the direct role played by the positive IEF y (or southward IMF B z ) in the storm-time AL intensification particularly when the terrestrial magnetosphere encounters the CME1 TE and IR region. Further, the arrival of CME2 did not affect the terrestrial magnetosphere-ionosphere system.

As aforementioned, major geomagnetic response is noticed during the arrival of trailing edge of the preceding CME (CME1) and the interaction region (IR) of the two CMEs near the Earth. We also conclude that the following CME (CME2) failed to cause a significant geomagnetic activity possibly due to the fact that the in situ spacecraft encountered the flank of this CME as mentioned in Section 3.1. We understand that due to interaction and collision between the trailing edge of the preceding CME1 and the leading edge of the following CME2 as revealed in imaging observations (described in Section 2.1.2), the parameters responsible for geomagnetic activity were significantly intensified at CME1 rear edge and in interaction region found between CME1 and CME2. These results bring out the importance of CME-CME interaction in the formation of interaction region and its role in the significant development of geomagnetic disturbances.

## 4. Results and Discussion

In this paper, we have studied the interactions between two successively launched Earth directed CMEs from the Sun on 2012 November 9 and 10, respectively. The analysis involved remote sensing observations of CME propagation in inner heliosphere from SECCHI/COR2 and HI images and the associated in situ measurements from WIND spacecraft. The first CME (CME1) associated with a filament eruption is observed to be slower (600 km s -1 ) than the successive CME (CME2) launched on November 10 with a higher velocity (900 km s -1 ). Both the CMEs were launched in approximately the same direction towards the Earth, indicating a possibility of interaction because of their relative speeds and the same direction. Different features were tracked as brightness enhancements

in J-maps obtained from COR2 and HI images and subsequently HM technique was applied to estimate the kinematics of the tracked features of the CMEs. From the estimated kinematics, the site of collision of CMEs could be located at a distance of approximately 35 R glyph[circledot] at 12:00 UT on November 10, which is at least 85 R glyph[circledot] before and 15 hours earlier than as predicted by using the initial kinematics of CMEs in COR2 FOV. The estimated kinematics were then used as inputs in the DBM to determine the arrival times of the tracked features. After the interaction, CME2 transferred some of its velocity to CME1 and both were observed to propagate together, showing a distinct deceleration of CME2.

As described earlier, CME1 LE is found to be propagating with higher speed than CME1 TE. Therefore, CME1 LE (green track) did not collide with CME2 LE (blue track). With the start of collision phase, the LE of CME1 is also found to accelerate and after the collision phase, significant acceleration of CME1 LE is noticed. This may occur either because of sudden impact (push) from CME1 TE to CME1 LE during collision or due to the passage of a shock driven by the CME2 or a combined effect of both. Beyond the collision phase, up to an estimated distance of nearly 100 R glyph[circledot] , LE of CME2 was found moving behind the TE of CME1 and both these features propagated together decelerating slowly. For few hours after the observed collision phase, it is noticed that speed of the LE and TE of CME1 is slightly higher than LE of CME2, which can increase the observed separation of these structures. The J-maps show that both these features can be tracked for further elongations after interaction. However, due to the limitation of the HM method for higher elongations (Lugaz, Vourlidas, and Roussev, 2009), we restricted our measurements on the tracked features up to 100 R glyph[circledot] . Our analyses show an evidence that after the interaction, the two features did not merge. It appears that LE of CME2 interacted with TE of CME1 and continued to propagate with a reduced speed of 500 km s -1 . In the SECCHI/HI-B images these features cannot be well distinguished and therefore it is difficult to use any stereoscopic method to infer about a possible longitudinal deflection (based on stereoscopic methods) of these features after interaction.

During the collision between CME1 and CME2, we noticed a large deceleration of fast CME2 while relatively less acceleration of slow CME1, till both approached an equal speed. This is expected to occur if the mass of CME1 is larger than that of CME2, which is indeed the case. In Section 2.2 we have shown that the mass of CME1 is ≈ 2.0 times larger than the mass of CME2. This result is an important finding and is in agreement with the second scenario of interaction described in Lugaz, Vourlidas, and Roussev (2009). It also must be noted that we have estimated the mass of CMEs in COR2 FOV while interaction takes place in HI FOV. Further, we cannot ignore the possibility of increase in mass of CME due to mass accretion at its front via snowplough effect in the solar wind beyond COR FOV (DeForest, Howard, and McComas, 2013). We must also mention that in our calculation of momentum and energy transfer, although we use the total mass of CME1 but only a part of CME1 (i.e. TE of CME1) takes part in collision with CME2. Considering the uncertainties in the derived speeds and mass of the CMEs, our analysis reveals that the nature of collision remains close to perfectly inelastic. However, we also acknowledge that defining collision phase is crucial and seems to play an important role in this calculation. In this

context, it must be emphasized that marking of the start of collision is often difficult. This is because the following fast speed CME2 starts to decelerate, well before the actual merging or collision is revealed in imaging observations. This is also shown in earlier studies for other interacting CMEs (Temmer et al. , 2012; Shen et al. , 2012; Mariˇi´ c c et al. , 2014; Mishra and Srivastava, 2014). The interaction of CME2 with streamers or overlying coronal magnetic field lines (Temmer et al. , 2008) may be responsible for deceleration of following CME2 before its collision with CME1 takes place. However, the contribution of other factors, e.g., the solar wind drag acting on CME2, interaction of CME2 with trailing edge of preceding CME1, and magnetic interaction of CME2 with preceding CME1 cannot be completely ignored. Also, possibly a shock launched by the following CME2 can contribute in acceleration of the preceding CME1 before the actual merging of both CMEs. Therefore, different timing and large time-interval of acceleration of one CME and deceleration of other, prevent us to pinpoint the exact start and end of the collision phase.

We concluded that the collision of interacting CMEs of 2012 November 910 is close to perfectly inelastic in nature. Here, it is noted that we are dealing with large scale magnetized three-dimensional structures, however, in the present study we have reduced the calculations to one-dimensional collision which seems to be accurate as tracked features of both CME1 and CME2 are considered to move along approximately the same trajectory. This implies that there will be a small difference between velocity of tracked feature of one CME and its component along the direction of propagation (i.e. longitude) of tracked feature of other CME. We must admit that possibility of change in longitude of tracked features during their collision can partially contribute to some errors in the estimated speed. In light of earlier reported observations (Shen et al. , 2012; Temmer et al. , 2012; Mishra and Srivastava, 2014), it seems that collision of CMEs can occur in all i.e. elastic, inelastic or super-elastic regimes. An in-depth study is required further for understanding the nature of CME-CME collision which may depend on the characteristics of CMEs, locations and duration of collision phase in the heliosphere.

At the last point of measurement in HI FOV around 120 R glyph[circledot] , the estimated speed of TE of CME1 and LE of CME2 are approximately equal ( ≈ 470 km s -1 ). Therefore, both features are expected to arrive at L1 at approximately the same time. Comparing to in situ measured actual arrival time, the delayed arrival of LE of CME2 is possible due to higher drag force acting on it resulting in its deceleration. However, the propagation direction of CME2 in southward hemisphere, as noticed in HI1-A images and also estimated in COR2 FOV using 3D reconstruction, can also account for flank encounter of CME2 and thus its delayed sampling by WIND spacecraft at L1. This is consistent with our interpretation in Section 3.1. Since, in the HM method, the front of the CME is assumed spherical which may not be the real case for CME2, as it was moving with fast speed and has large probability to distort and flatten its front while interacting with structured solar wind (Odstrcil, Pizzo, and Arge, 2005). In this case, the estimated speed using HM method will lead to some error and can result in the delayed arrival of LE of CME2. Keeping all these issues in mind, we believe that it is quite probable that the tracked feature corresponding to LE of CME2

(blue track in J-map) is not sampled by in situ spacecraft. Therefore, erroneous predicted arrival time for LE of CME2 is related to incorrect identification of tracked CME2 LE in in situ data which is taken as a reference (actual) for remotely observed tracked CME2 LE.

The present study not only estimates the kinematics and arrival times based on HI observations and their application to the DBM model of CME propagation, it also attempts to associate tracked brightness enhancements as observed in remote sensing COR and HI (J-maps) with the in situ observations (M¨stl o et al. , 2010, 2011; Rollett et al. , 2012, 2013). The association of HI observations with in-situ measurements from WIND leads to many interesting results. In spite of two CMEs launched from the Sun in succession in the Earthward direction, we observe only one shock in in situ data which may suggest the merging of shocks driven by CME1 and CME2, if both CMEs would have driven shocks. However, such a claim cannot be made unless we are well familiar with the in situ signatures of merged shock and plasma structure following it. The sweeping of plasma to high density at the front of CME1 and its compressed heating is most likely due to passage of CME2 driven shock through the MC associated with CME1. Based on the predicted arrival time of tracked features, it seems that the CME2 driven shock and CME1-sheath region is tracked as CME1 LE in J-map. Therefore, we infer that CME1 LE is propagating probably into an unperturbed solar wind. Our study also provides a possibility of formation of interaction region (IR) at the junction of trailing edge of preceding CME and leading edge of following CME. We show that during collision of the CMEs, kinetic energy exchange up to 50% and momentum exchange between 23 to 30%. Our study also demonstrate that the arrival time prediction significantly improved using HIs on-board STEREO compared to COR2 observations, and also emphasizes the importance of understanding of post-collision kinematics in further improving the arrival time prediction for a reliable space weather prediction scheme.

Our study reveals clear signatures of interaction of these CMEs in remote and in situ observations and also helps in identification of separate structures corresponding to these CMEs. In spite of the interaction of the two CMEs in the interplanetary medium which generally results in complex structures as suggested by Burlaga, Behannon, and Klein (1987); Burlaga, Plunkett, and St. Cyr (2002), in our case, we could identify interacted CMEs as distinguished structures in WIND spacecraft data. Even after collision of these CMEs, they did not merge which may be possibly because of strong magnetic field and higher density of CME1 than CME2. This needs further confirmation, and therefore, it is worth to investigate what decides the formation of merged CME structure or complex ejecta during CME-CME interaction. Due to single point in situ observations of CME, we acknowledge the possibility of ambiguity in marking the boundaries of CMEs. In the present case, the boundaries for CME1 and IR are distinctly clear. Also, slight ambiguity in the boundary of CME2 (flank encounter) will not change our interpretation because the main geomagnetic response is caused due to enhanced negative Bz in trailing portion of CME1 and its extension in IR. Here, we also point out that temperature in IR is lesser than the temperature for CME2 region but temperatures in both regions are

elevated as compared to a normal non-interacting CME. The observations of unexpected larger temperature in CME2 region than IR region may be due to the possibility that sheath region of CME2 is intersected by in situ spacecraft, as we have interpreted.

The association of geomagnetic storms with isolated single CMEs has been carried out extensively for a long time (Tsurutani et al. , 1988; Gosling et al. , 1991; Gosling, 1993; Gonzalez et al. , 1989, 1994; Echer et al. , 2008; Richardson and Cane, 2011). However, only few studies have been dedicated to understand the role of interacting CMEs in the generation of geomagnetic storms (Burlaga, Behannon, and Klein, 1987; Farrugia et al. , 2006; Wang, Ye, and Wang, 2003; Xiong et al. , 2006). Our study is important as it focuses on the role of interacting CMEs in the generation of geomagnetic storms as well as substorms. In the context of substorms, our study highlights that persistence of IMF Bz in the southward direction is more important than the amplitude in driving the substorm activity as manifested by the AL intensification. Using WINDMI model, Mays et al. (2007) have shown that the IP shock and sheath features for CMEs contribute significantly to the development of storms and substorms. But, in our study, sharp and large southward excursions in the midst of fluctuating IMF Bz associated with shock (shock-sheath region before CME1) were found less effective in producing strong substorm activity. Therefore, further investigations are required regarding characteristics (geometry, intensity) of shock and preceding CME in the context of triggering of substorms, as has been shown in earlier studies (Jurac et al. , 2002; Wang et al. , 2003). Regarding geomagnetic storms, our study suggests that the trailing edge of preceding CME (CME1) and IR formed between two interacting CMEs are efficient candidates for intense geomagnetic storms.

Another interesting aspect regarding substorms noticed in our study is that the nearly equal amplitude responses of the PC index corresponding to the shocksheath region (IMF Bz sharply fluctuating between southward and northward directions) preceding the CME1 interval vis-`-vis the CME1 interval (IMF Bz a steadily turning southward). This is interesting as the responses of Sym-H and AL during these two intervals are quite different in terms of amplitudes of variations. Janzhura, Troshichev, and Stauning (2007) have inferred that the PC index for the sunlit polar cap (summer hemisphere, high polar ionospheric conductance) responds mainly to the geo-effective interplanetary electric field whereas the PC index for the dark winter cap (winter hemisphere, low polar ionospheric conductance) responds better to the particle precipitation in the auroral zone like the AE and AL indices. It is to be noted here that, in the present study, no distinction is made between the variations of PC index in the northern (PCN) and southern (PCS) hemispheres. On the other hand, it is known that AL index is constructed based on magnetometer observations around the northern auroral oval only. The discrepancy between the variations in the PC index and the AL index may, therefore, suggest towards the hemispherical asymmetry in the response of polar ionosphere corresponding to this event. This is not unexpected as asymmetric auroral intensities in the northern and southern hemispheres are reported in literature (Laundal and Østgaard, 2009). Further investigations are needed to understand this aspect.

## 5. Conclusions

The CMEs of November 9 and 10, provide us a rare opportunity to investigate the consequences of CME-CME interaction. A combination of heliospheric imaging and in situ observations are used for improving our understanding of CME kinematics, post-collision characteristics and nature of collision. Main results of our study can be concluded as follows:

- i) The analysis of propagation kinematics as obtained from J-maps provide evidence that the CME1 and CME2 collide at 35 R glyph[circledot] much earlier than by using the estimated kinematics in the COR2 FOV. This emphasizes the importance of heliospheric imaging particularly, for interacting CMEs and ascertaining their impact and arrival at the Earth.
- ii) Our analysis shows that post interaction and collision kinematics is required for a better prediction of CME arrival time at 1 AU. In the present case of November 9 and 10 CMEs, the speeds and momentum of CMEs changed from 23% to 30% compared to their values before the collision. Our results also highlight that estimated kinematics, in particular after collisions are important to combine with DBM for improving the estimation of arrival times of different features of CME which experience different drag forces during their propagation through the heliosphere. Based on this study, we conclude that CMEs cannot be treated as completely isolated magnetized plasma blobs, especially when they are launched in quick succession.
- iii) Using estimated mass and kinematics before and after the collision, we estimated the total kinetic energy of the system before and after the collision. We found that total kinetic energy of the system decreased by 6.7% to its value before the collision and the nature of collision is close to perfectly inelastic. The study also supports the idea that CME interaction or collision can lead to the heating and compression of both preceding and following CMEs.
- iv) It is clear that if the trailing edge of November 9 CME would not have been tracked using imaging observation, one could not have witnessed the collision. Further, trailing edge of this CME had strong negative component of magnetic field for long duration ( ≈ 13 hr) which may be due to its collision with the leading edge of November 10 CME. Therefore, our study reveals that tracking of different features of CME seems to be necessary for better understanding of CME-CME interaction.
- v) It is observed in the present case that the persistence of IMF B z in the southward direction is more important (rather than the amplitude) in driving the substorm activity as manifested in the AL intensification. Sharp and equally large southward excursions in the midst of fluctuating IMF B z associated with shock were found less effective in producing equally strong substorm activity. It is identified that the interaction region (IR), formed due to collision between CME1 TE and CME2 LE, has intensified plasma and magnetic field parameters which are responsible for major geomagnetic activity.

The present case study of interacting CMEs of 2012 November 9-10, highlight the importance of heliospheric imaging for estimating the kinematics of CME features before and after their collision and interactions. Our study support that

interacting CMEs can result in strong geomagnetic storms and also in substorms. Further statistical study of such interacting CMEs is required to understand their nature of collision and to investigate the characteristics of the CMEs (mass, strength and orientation of magnetic field, speed and direction of propagation, duration of collision phase) which are responsible for their interaction.

Acknowledgements We acknowledge the UK Solar System Data Centre for providing the processed Level-2 STEREO/HI data. The in situ measurements of solar wind data from ACE and WIND spacecraft were obtained from NASA CDAweb ( http://cdaweb.gsfc.nasa.gov/ ). We also acknowledge the use of DBM, developed by Bojan Vr˘nak and available at s http://oh.geof. unizg.hr/CADBM/cadbm.php , in this study. We also thank the reviewer for useful comments, which helped improve the manuscript. The work by N.S. partially contributes to the research for European Union Seventh Framework Programme (FP7/2007-2013) for the Coronal Mass Ejections and Solar Energetic Particles (COMESEP) project under Grant Agreement No. 263252.

Table 1. First column shows the tracked features of CMEs and second column shows their estimated kinematics (by HM technique) which is used as input to the DBM. The predicted arrival time and transit speed of the tracked features at L1, corresponding to the extreme range of the drag parameter used in the DBM, is shown in the third and fourth column. Errors in predicted arrival time and speed, based on comparison with in situ arrival time and speeds, are shown in column fifth and sixth. The errors in arrival time with negative (positive) sign indicate that predicted arrival is earlier (later) than the actual arrival time. The errors in transit speed with negative (positive) sign indicate that the predicted transit speed is lesser (greater) than the transit speed measured in situ.

| Tracked Features   | Kinematics as in- puts in DBM [t 0 , R 0 (R glyph[circledot] ), v 0 (km s - 1 )]   | Predicted arrival time (UT) using kinematics + DBM [ γ = 0.2 - 2.0 (10 - 7 km - 1 )]   | Predicted transit speed (km s - 1 ) at L1 [ γ =0.2 - 2.0 (10 - 7 km - 1 )]   | Error in predicted ar- rival time (hr) [ γ = 0.2 - 2.0 (10 - 7 km - 1 )]   | Error in predicted speed (km s - 1 ) [ γ = 0.2 - 2.0 (10 - 7 km - 1 )]   |
|--------------------|------------------------------------------------------------------------------------|----------------------------------------------------------------------------------------|------------------------------------------------------------------------------|----------------------------------------------------------------------------|--------------------------------------------------------------------------|
| CME1 LE            | Nov 11 13:42, 545, 137                                                             | Nov 12 18:10 - Nov 13 00:30                                                            | 490 to 380                                                                   | -5 to 1.5                                                                  | 115 to 5                                                                 |
| CME1 TE            | Nov 11 18:35, 120, 470                                                             | Nov 13 15:25                                                                           | 375                                                                          | -8                                                                         | -15                                                                      |
| CME2 LE            | Nov 12 00:30, 124, 455                                                             | Nov 13 19:40                                                                           | 375                                                                          | -16                                                                        | -35                                                                      |

Figure 1. From top to bottom panels, true height, velocity, acceleration, longitude and latitude of selected feature along the leading edge as derived from Tie-Pointing method have been plotted as a function of time for CME1 (left) and CME2 (right)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0352v1-figure-215b3a2a-3a0a-4b37-9037-999f0de4e04e.jpg)

Figure 2. The top and bottom rows show the GCS model wireframe overlaid on the CME1 and CME2 images, respectively. The triplet of concurrent images around 17:39 UT on November 9 and around 06:39 UT on November 10 corresponding to CME1 and CME2, respectively, are from STEREO/COR2-B (left), SOHO/LASCO-C3 (middle) and STEREO/COR2-A (right).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0352v1-figure-77bd80e2-a59c-4340-ac25-83e288b764e9.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0352v1-figure-e9968077-5a07-411e-8474-e899ba3c6fa6.jpg)

Earth PA:

2012-11-10T02.09.01.003 STEREO Hl1 Spacecraft Longitude (HEEQ): 127.190 Spacecraft Latitudee (HEEQ):

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0352v1-figure-d8de83ca-89e7-4306-886e-228a79efa732.jpg)

Earth PA:

91.260

2012-11-11722.09.21.003 STEREO A HI2 Longitude (HEEQ) 127,.290 STEREO\_A Latitude (HEEQ): 7.189

7.230

Earth PA: 91.280

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0352v1-figure-79a4ebda-39e4-4ad5-8cf8-95ec0dfbe754.jpg)

2012-11-10710.09.01.004 STEREO HI1 Spacecraft Longitude (HEEQ): 127.209 Spacecraft Latitudee (HEEQ):

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0352v1-figure-97d7a34a-30dd-44fe-a542-f6421dfd16a4.jpg)

Earth PA:

91.300

2012-11-12T04.09.21.004 STEREO HI2 Longitude (HEEQ) 127,31 STEREO HI2  Latitude (HEEQ):

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0352v1-figure-7d1fb200-aced-4b77-bdde-a993c9500bcb.jpg)

Earth PA:

91.500

Earth PA: 91.530

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0352v1-figure-0db74638-25ea-487c-ae35-bbf5dee355d7.jpg)

Figure 3. Evolution of CMEs in running difference images of COR2, HI1 and HI2 FOV is shown in upper, middle and lower panels respectively. Left and right panels show observations from STEREO/SECCHI Ahead spacecraft at two different times. Contours of elongation (green) and position angle (blue) is overlaid on the images. In each image, the postion angle is overlaid in interval of 10 ◦ and the horizontal red line is along the ecliptic at the position angle of Earth. In upper panel (both left and right) the vertical red line marks the 0 ◦ position angle. In each panels, the C1 and C2 corresponds to CME1 and CME2.

## STEREO HI and COR-2 Spacecralt A: Ecliptic

Figure 4. Time-elongation map (J-map) using the COR2 and HI observations of STEREO/SECCHI spacecraft during interval of 2012 November 9-14 is shown. The features corresponding to CME1 leading edge (LE), CME1 trailing edge (TE) and CME2 leading edge are (LE) tracked and over plotted on the J-map with green, red and blue, respectively. The red rectangle (rightmost) is an enlarged plot of the red rectangle (on the left) which clearly shows that the red and blue tracks meet in HI1 FOV.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0352v1-figure-588dbbda-de4c-4c1c-abc5-75d4d1034c94.jpg)

Figure 5. The contours of derived elongation of tracked CME1 LE (green), CME1 TE (red) and CME2 LE (blue) features from the J-map is overplotted on the base difference HI1-A images. In each image, the position angle (sky blue) is overlaid in interval of 10 ◦ and the horizontal red line is along the ecliptic at the position angle of Earth.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0352v1-figure-3844e1a7-ca69-4f04-bb4a-1f9a59e9956b.jpg)

Figure 6. In the top and bottom panels, distance and speed of the tracked features CME1 LE, CME1 TE and CME2 LE as marked in Figure 4 with green, red and blue, respectively, are shown. Speeds of these features are calculated from the differentiation of adjacent distance points using the three point Lagrange interpolation (shown with ∆ symbol in the bottom panel). Speed of these features are estimated by differentiating the first order polynomial fit to estimated distance points for an interval of approximately 5 hr (shown in solid line in the bottom panel). In the top and bottom panels, vertical lines show the error bars. We have assumed a fractional error of 5% in estimated distance which is used to determine the uncertainties in the speeds. From the left, the two vertical lines (black) mark the start and end of the collision phase.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0352v1-figure-91679bd3-8df2-48c2-b9bf-2feb4dd59f42.jpg)

Figure 7. The best suited coefficient of restitution ( e ) corresponding to different mass ratios of CME1 and CME2 are shown with + symbol and corresponding variance ( σ ) in velocity is shown with ∆ symbol.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0352v1-figure-19854495-5392-4b87-9151-77ee00b901e4.jpg)

Figure 8. From top to bottom, total magnetic field magnitude, z-component of magnetic field, proton density, proton temperature, proton speed, latitude and longitude of magnetic field vector, and plasma beta ( β ) is shown for the time interval of 12:00 UT on November 12 to 12:00 UT on November 15. From the left, first, second, third, fourth and fifth vertical (dashed) lines mark the arrival of shock, CME1 leading edge (LE), CME1 trailing edge (TE), CME2 LE and CME2 TE, respectively. Interaction region (IR) is shown during the interval of third and fourth vertical lines.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0352v1-figure-d4f2dcdd-0bf3-4a6a-b5f0-56d9f1d748f9.jpg)

Figure 9. Plot shows variations of geomagnetic field parameters during the interval 2012 November 12-15. Left panels show, top to bottom: Proton density (n/cc), Velocity (km/s), Ram Pressure (nPa) and B z component (nT). Right panels show, from top to bottom: interplanetary electric field's 'y' component IEF y (mV/m), polar cap (PC) index, Sym-H (nT) and AL index (nT). From the left, vertical lines and their labels are as defined in Figure 8.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0352v1-figure-f56b7fec-c10b-4cb8-97c4-9f7a431cbbb5.jpg)

## References

Billings, D.E.: 1966, A guide to the solar corona , San Diego, CA: Academic.

Burlaga, L.F., Lemaire, J.F.: 1978, Interplanetary magnetic holes - Theory. J. Geophys. Res. 83 , 5157 - 5160. doi: 10.1029/JA083iA11p05157 .

Burlaga, L.F., Behannon, K.W., Klein, L.W.: 1987, Compound streams, magnetic clouds, and major geomagnetic storms. J. Geophys. Res. 92 , 5725 - 5734. doi: 10.1029/JA092iA06p05725 .

Burlaga, L.F., Plunkett, S.P., St. Cyr, O.C.: 2002, Successive CMEs and complex ejecta. J. Geophys. Res. 107 , 1266. doi: 10.1029/2001JA000255 .

- Cargill, P.J.: 2004, On the Aerodynamic Drag Force Acting on Interplanetary Coronal Mass Ejections. Solar Phys. 221 , 135 - 149. doi: 10.1023/B:SOLA.0000033366.10725.a2 .
- Chakrabarty, D., Sekar, R., Narayanan, R., Devasia, C.V., Pathan, B.M.: 2005, Evidence for the interplanetary electric field effect on the OI 630.0 nm airglow over low latitude. Journal of Geophysical Research (Space Physics) 110 , 11301. doi: 10.1029/2005JA011221 .
- Colaninno, R.C., Vourlidas, A.: 2009, First Determination of the True Mass of Coronal Mass Ejections: A Novel Approach to Using the Two STEREO Viewpoints. Astrophys. J. 698 , 852 - 858. doi: 10.1088/0004-637X/698/1/852 .

Colaninno, R.C., Vourlidas, A., Wu, C.C.: 2013, Quantitative comparison of methods for predicting the arrival of coronal mass ejections at Earth based on multiview imaging. J. Geophys. Res. 118 , 6866 - 6879. doi: 10.1002/2013JA019205 .

Davies, J.A., Harrison, R.A., Rouillard, A.P., Sheeley, N.R., Perry, C.H., Bewsher, D., Davis, C.J., Eyles, C.J., Crothers, S.R., Brown, D.S.: 2009, A synoptic view of solar transient evolution in the inner heliosphere using the Heliospheric Imagers on STEREO. Geophys. Res. Lett. 36 , 2102. doi: 10.1029/2008GL036182 .

Davies, J.A., Perry, C.H., Trines, R.M.G.M., Harrison, R.A., Lugaz, N., M¨ ostl, C., Liu, Y.D., Steed, K.: 2013, Establishing a stereoscopic technique for determining the kinematic properties of solar wind transients based on a generalised self-similarly expanding circular geometry. Astrophys. J. 776 , 1.

DeForest, C.E., Howard, T.A., McComas, D.J.: 2013, Tracking Coronal Features from the Low Corona to Earth: A Quantitative Analysis of the 2008 December 12 Coronal Mass Ejection. Astrophys. J. 769 , 43. doi: 10.1088/0004-637X/769/1/43 .

Echer, E., Gonzalez, W.D., Tsurutani, B.T., Gonzalez, A.L.C.: 2008, Interplanetary conditions causing intense geomagnetic storms (Dst = -100 nT) during solar cycle 23 (1996-2006). J. Geophys. Res. 113 , 5221. doi: 10.1029/2007JA012744 .

Farrugia, C.J., Jordanova, V.K., Thomsen, M.F., Lu, G., Cowley, S.W.H., Ogilvie, K.W.: 2006, A two-ejecta event associated with a two-step geomagnetic storm. J. Geophys. Res. 111 , 11104. doi: 10.1029/2006JA011893 .

Gonzalez, W.D., Gonzalez, A.L.C., Tsurutani, B.T., Smith, E.J., Tang, F.: 1989, Solar windmagnetosphere coupling during intense magnetic storms (1978-1979). J. Geophys. Res. 94 , 8835 - 8851. doi: 10.1029/JA094iA07p08835 .

Gonzalez, W.D., Joselyn, J.A., Kamide, Y., Kroehl, H.W., Rostoker, G., Tsurutani, B.T., Vasyliunas, V.M.: 1994, What is a geomagnetic storm? J. Geophys. Res. 99 , 5771 - 5792. doi: 10.1029/93JA02867 .

Gonzalez-Esparza, A., Santill´n, a A., Ferrer, J.: 2004, A numerical study of the interaction between two ejecta in the interplanetary medium: one- and two-dimensional hydrodynamic simulations. Annales Geophysicae 22 , 3741 - 3749. doi: 10.5194/angeo-22-3741-2004 .

Gopalswamy, N., Lara, A., Lepping, R.P., Kaiser, M.L., Berdichevsky, D., St. Cyr, O.C.: 2000, Interplanetary acceleration of coronal mass ejections. Geophys. Res. Lett. 27 , 145 - 148. doi: 10.1029/1999GL003639 .

Gopalswamy, N., Yashiro, S., Kaiser, M.L., Howard, R.A., Bougeret, J.-L.: 2001, Radio Signatures of Coronal Mass Ejection Interaction: Coronal Mass Ejection Cannibalism? Astrophys. J. Lett. 548 , L91-L94. doi: 10.1086/318939 .

Gosling, J.T.: 1993, The solar flare myth. J. Geophys. Res. 98 , 18937 - 18950. doi: 10.1029/93JA01896 .

- Gosling, J.T., Bame, S.J., McComas, D.J., Phillips, J.L.: 1990, Coronal mass ejections and large geomagnetic storms. Geophys. Res. Lett. 17 , 901 - 904. doi: 10.1029/GL017i007p00901 .
- Gosling, J.T., McComas, D.J., Phillips, J.L., Bame, S.J.: 1991, Geomagnetic activity associated with earth passage of interplanetary shock disturbances and coronal mass ejections. J. Geophys. Res. 96 , 7831 - 7839. doi: 10.1029/91JA00316 .

Harrison, R.A., Davies, J.A., M¨stl, C., Liu, Y., Temmer, M., Bisi, M.M., Eastwood, J.P., de o Koning, C.A., Nitta, N., Rollett, T., Farrugia, C.J., Forsyth, R.J., Jackson, B.V., Jensen, E.A., Kilpua, E.K.J., Odstrcil, D., Webb, D.F.: 2012, An Analysis of the Origin and Propagation of the Multiple Coronal Mass Ejections of 2010 August 1. Astrophys. J. 750 , 45. doi: 10.1088/0004-637X/750/1/45 .

Howard, R.A., Moses, J.D., Vourlidas, A., Newmark, J.S., Socker, D.G., Plunkett, S.P., Korendyke, C.M., Cook, J.W., Hurley, A., Davila, J.M., Thompson, W.T., St Cyr, O.C., Mentzell, E., Mehalick, K., Lemen, J.R., Wuelser, J.P., Duncan, D.W., Tarbell, T.D., Wolfson, C.J., Moore, A., Harrison, R.A., Waltham, N.R., Lang, J., Davis, C.J., Eyles, C.J., Mapson-Menard, H., Simnett, G.M., Halain, J.P., Defise, J.M., Mazy, E., Rochus, P., Mercier, R., Ravet, M.F., Delmotte, F., Auchere, F., Delaboudiniere, J.P., Bothmer, V., Deutsch, W., Wang, D., Rich, N., Cooper, S., Stephens, V., Maahs, G., Baugh, R., McMullin, D., Carter, T.: 2008, Sun Earth Connection Coronal and Heliospheric Investigation (SECCHI). Space Sci. Rev. 136 , 67 - 115. doi: 10.1007/s11214-008-9341-4 .

Howard, T.A., Tappin, S.J.: 2009, Interplanetary Coronal Mass Ejections Observed in the Heliosphere: 1. Review of Theory. Space Sci. Rev. 147 , 31 - 54. doi: 10.1007/s11214-0099542-5 .

Inhester, B.: 2006, Stereoscopy basics for the STEREO mission. arXiv:astro-ph/0612649 .

Iyemori, T., Rao, D.R.K.: 1996, Decay of the Dst field of geomagnetic disturbance after substorm onset and its implication to storm-substorm relation. Annales Geophysicae 14 , 608 - 618. doi: 10.1007/s00585-996-0608-3 .

Janzhura, A., Troshichev, O., Stauning, P.: 2007, Unified PC indices: Relation to isolated magnetic substorms. Journal of Geophysical Research (Space Physics) 112 , 9207. doi: 10.1029/2006JA012132 .

Jurac, S., Kasper, J.C., Richardson, J.D., Lazarus, A.J.: 2002, Geomagnetic disturbances and their relationship to Interplanetary shock parameters. Geophys. Res. Lett. 29 , 1463. doi: 10.1029/2001GL014034 .

Kaiser, M.L., Kucera, T.A., Davila, J.M., St. Cyr, O.C., Guhathakurta, M., Christian, E.: 2008, The STEREO Mission: An Introduction. Space Sci. Rev. 136 , 5 - 16. doi: 10.1007/s11214007-9277-0 .

Kilpua, E.K.J., Mierla, M., Rodriguez, L., Zhukov, A.N., Srivastava, N., West, M.J.: 2012, Estimating Travel Times of Coronal Mass Ejections to 1 AU Using Multi-spacecraft Coronagraph Data. Solar Phys. 279 , 477 - 496. doi: 10.1007/s11207-012-0005-x .

Laundal, K.M., Østgaard, N.: 2009, Asymmetric auroral intensities in the Earth's Northern and Southern hemispheres. Nature 460 , 491 - 493. doi: 10.1038/nature08154 .

Lepping, R.P., Ac˜na, M.H., Burlaga, L.F., Farrell, W.M., Slavin, J.A., Schatten, K.H., Maru iani, F., Ness, N.F., Neubauer, F.M., Whang, Y.C., Byrnes, J.B., Kennon, R.S., Panetta, P.V., Scheifele, J., Worley, E.M.: 1995, The Wind Magnetic Field Investigation. Space Sci. Rev. 71 , 207 - 229. doi: 10.1007/BF00751330 .

Lindsay, G.M., Luhmann, J.G., Russell, C.T., Gosling, J.T.: 1999, Relationships between coronal mass ejection speeds from coronagraph images and interplanetary characteristics of associated interplanetary coronal mass ejections. J. Geophys. Res. 104 , 12515 - 12524. doi: 10.1029/1999JA900051 .

Liu, Y.D., Luhmann, J.G., M¨stl, C., Martinez-Oliveros, J.C., Bale, S.D., Lin, R.P., Harrison, o R.A., Temmer, M., Webb, D.F., Odstrcil, D.: 2012, Interactions between Coronal Mass Ejections Viewed in Coordinated Imaging and in situ Observations. Astrophys. J. Lett. 746 , L15. doi: 10.1088/2041-8205/746/2/L15 .

Liu, Y., Davies, J.A., Luhmann, J.G., Vourlidas, A., Bale, S.D., Lin, R.P.: 2010a, Geometric Triangulation of Imaging Observations to Track Coronal Mass Ejections Continuously Out to 1 AU. Astrophys. J. Lett. 710 , L82-L87. doi: 10.1088/2041-8205/710/1/L82 .

Liu, Y., Thernisien, A., Luhmann, J.G., Vourlidas, A., Davies, J.A., Lin, R.P., Bale, S.D.: 2010b, Reconstructing Coronal Mass Ejections with Coordinated Imaging and in Situ Observations: Global Structure, Kinematics, and Implications for Space Weather Forecasting. Astrophys. J. 722 , 1762 - 1777. doi: 10.1088/0004-637X/722/2/1762 .

Lugaz, N., Manchester, W.B. IV, Gombosi, T.I.: 2005, Numerical Simulation of the Interaction of Two Coronal Mass Ejections from Sun to Earth. Astrophys. J. 634 , 651 - 662. doi: 10.1086/491782 .

Lugaz, N., Vourlidas, A., Roussev, I.I.: 2009, Deriving the radial distances of wide coronal mass ejections from elongation measurements in the heliosphere - application to CME-CME interaction. Annales Geophysicae 27 , 3479 - 3488. doi: 10.5194/angeo-27-3479-2009 .

Lugaz, N., Hernandez-Charpak, J.N., Roussev, I.I., Davis, C.J., Vourlidas, A., Davies, J.A.: 2010, Determining the Azimuthal Properties of Coronal Mass Ejections from MultiSpacecraft Remote-Sensing Observations with STEREO SECCHI. Astrophys. J. 715 , 493 - 499. doi: 10.1088/0004-637X/715/1/493 .

Lugaz, N., Farrugia, C.J., Davies, J.A., M¨ ostl, C., Davis, C.J., Roussev, I.I., Temmer, M.: 2012, The Deflection of the Two Interacting Coronal Mass Ejections of 2010 May 23-24 as Revealed by Combined in Situ Measurements and Heliospheric Imaging. Astrophys. J. 759 , 68. doi: 10.1088/0004-637X/759/1/68 .

Manoharan, P.K.: 2006, Evolution of Coronal Mass Ejections in the Inner Heliosphere: A Study Using White-Light and Scintillation Images. Solar Phys. 235 , 345 - 368. doi: 10.1007/s11207006-0100-y .

Mariˇ ci´, c D., Vrˇ snak, B., Dumbovi´ c, M., Zic, T., Roˇa, D., Hrˇina, D., Luli´, S., Romˇtajn, ˇ s z c s I., Buˇ si´, c I., Salamon, K., Temmer, M., Rollett, T., Veronig, A., Bostanjyan, N., Chilingarian, A., Mailyan, B., Arakelyan, K., Hovhannisyan, A., Muji´, N.: 2014, Kinematics of c Interacting ICMEs and Related Forbush Decrease: Case Study. Solar Phys. 289 , 351 - 368. doi: 10.1007/s11207-013-0314-8 .

Mart´ ınez Oliveros, J.C., Raftery, C.L., Bain, H.M., Liu, Y., Krupar, V., Bale, S., Krucker, S.: 2012, The 2010 August 1 Type II Burst: A CME-CME Interaction and its Radio and White-light Manifestations. Astrophys. J. 748 , 66. doi: 10.1088/0004-637X/748/1/66 .

Mays, M.L., Horton, W., Kozyra, J., Zurbuchen, T.H., Huang, C., Spencer, E.: 2007, Effect of Interplanetary Shocks on the AL and Dst Indices. Geophys. Res. Lett. 34 , 11104. doi: 10.1029/2007GL029844 .

Mierla, M., Davila, J., Thompson, W., Inhester, B., Srivastava, N., Kramar, M., St. Cyr, O.C., Stenborg, G., Howard, R.A.: 2008, A Quick Method for Estimating the Propagation Direction of Coronal Mass Ejections Using STEREO-COR1 Images. Solar Phys. 252 , 385 396. doi: 10.1007/s11207-008-9267-8 .

Minnaert, M.: 1930, On the continuous spectrum of the corona and its polarisation. With 3 figures. (Received July 30, 1930). Z. Astrophys. 1 , 209.

Mishra, W., Srivastava, N.: 2013, Estimating the Arrival Time of Earth-directed Coronal Mass Ejections at in Situ Spacecraft Using COR and HI Observations from STEREO. Astrophys. J. 772 , 70. doi: 10.1088/0004-637X/772/1/70 .

Mishra, W., Srivastava, N.: 2014, Morphological and kinematic evolution of three interacting coronal mass ejections of 2011 february 13-15. Astrophys. J. , Under review.

Mishra, W., Srivastava, N., , Davies, J.A.: 2014, A comparison of reconstruction methods for the estimation of coronal mass ejections kinematics based on secchi/hi observations. Astrophys. J. , in press. doi: 10.1088/0004-637X/784/1/1 .

M¨ ostl, C., Temmer, M., Rollett, T., Farrugia, C.J., Liu, Y., Veronig, A.M., Leitner, M., Galvin, A.B., Biernat, H.K.: 2010, STEREO and Wind observations of a fast ICME flank triggering a prolonged geomagnetic storm on 5-7 April 2010. Geophys. Res. Lett. 37 , 24103. doi: 10.1029/2010GL045175 .

M¨ ostl, C., Rollett, T., Lugaz, N., Farrugia, C.J., Davies, J.A., Temmer, M., Veronig, A.M., Harrison, R.A., Crothers, S., Luhmann, J.G., Galvin, A.B., Zhang, T.L., Baumjohann, W., Biernat, H.K.: 2011, Arrival Time Calculation for Interplanetary Coronal Mass Ejections with Circular Fronts and Application to STEREO Observations of the 2009 February 13 Eruption. Astrophys. J. 741 , 34. doi: 10.1088/0004-637X/741/1/34 .

M¨ ostl, C., Farrugia, C.J., Kilpua, E.K.J., Jian, L.K., Liu, Y., Eastwood, J.P., Harrison, R.A., Webb, D.F., Temmer, M., Odstrcil, D., Davies, J.A., Rollett, T., Luhmann, J.G., Nitta, N., Mulligan, T., Jensen, E.A., Forsyth, R., Lavraud, B., de Koning, C.A., Veronig, A.M., Galvin, A.B., Zhang, T.L., Anderson, B.J.: 2012, Multi-point Shock and Flux Rope Analysis of Multiple Interplanetary Coronal Mass Ejections around 2010 August 1 in the Inner Heliosphere. Astrophys. J. 758 , 10. doi: 10.1088/0004-637X/758/1/10 .

Munro, R.H., Gosling, J.T., Hildner, E., MacQueen, R.M., Poland, A.I., Ross, C.L.: 1979, The association of coronal mass ejection transients with other forms of solar activity. Solar Phys. 61 , 201 - 215. doi: 10.1007/BF00155456 .

Odstrcil, D., Pizzo, V.J., Arge, C.N.: 2005, Propagation of the 12 May 1997 interplanetary coronal mass ejection in evolving solar wind structures. Journal of Geophysical Research (Space Physics) 110 , 2106. doi: 10.1029/2004JA010745 .

Ogilvie, K.W., Chornay, D.J., Fritzenreiter, R.J., Hunsaker, F., Keller, J., Lobell, J., Miller, G., Scudder, J.D., Sittler, E.C. Jr., Torbert, R.B., Bodet, D., Needell, G., Lazarus, A.J., Steinberg, J.T., Tappan, J.H., Mavretic, A., Gergin, E.: 1995, SWE, A Comprehensive Plasma Instrument for the Wind Spacecraft. Space Sci. Rev. 71 , 55 - 77. doi: 10.1007/BF00751326 .

Poland, A.I., Howard, R.A., Koomen, M.J., Michels, D.J., Sheeley, N.R. Jr.: 1981, Coronal transients near sunspot maximum. Solar Phys. 69 , 169 - 175. doi: 10.1007/BF00151264 .

Richardson, I.G., Cane, H.V.: 2011, Geoeffectiveness (Dst and Kp) of interplanetary coronal mass ejections during 1995-2009 and implications for storm forecasting. Space Weather 9 , 7005. doi: 10.1029/2011SW000670 .

Rollett, T., M¨stl, C., Temmer, M., Veronig, A.M., Farrugia, C.J., Biernat, H.K.: 2012, Cono straining the Kinematics of Coronal Mass Ejections in the Inner Heliosphere with In-Situ Signatures. Solar Phys. 276 , 293 - 314. doi: 10.1007/s11207-011-9897-0 .

Rollett, T., Temmer, M., M¨stl, C., Lugaz, N., Veronig, A.M., M¨stl, U.V.: 2013, Assessing o o the Constrained Harmonic Mean Method for Deriving the Kinematics of ICMEs with a Numerical Simulation. Solar Phys. 283 , 541 - 556. doi: 10.1007/s11207-013-0246-3 .

Sheeley, N.R., Walters, J.H., Wang, Y.-M., Howard, R.A.: 1999, Continuous tracking of coronal outflows: Two kinds of coronal mass ejections. J. Geophys. Res. 104 , 24739 - 24768. doi: 10.1029/1999JA900308 .

Shen, C., Wang, Y., Wang, S., Liu, Y., Liu, R., Vourlidas, A., Miao, B., Ye, P., Liu, J., Zhou, Z.: 2012, Super-elastic collision of large-scale magnetized plasmoids in the heliosphere. Nature 8 , 923 - 928. doi: 10.1038/nphys2440 .

Srivastava, N., Inhester, B., Mierla, M., Podlipnik, B.: 2009, 3D Reconstruction of the Leading Edge of the 20 May 2007 Partial Halo CME. Solar Phys. 259 , 213 - 225. doi: 10.1007/s11207009-9423-9 .

Temmer, M., Veronig, A.M., Vrˇ snak, B., Ryb´k, J., G¨m¨ry, P., Stoiser, S., Mariˇi´, D.: 2008, a o o c c Acceleration in Fast Halo CMEs and Synchronized Flare HXR Bursts. Astrophys. J. Lett. 673 , L95-L98. doi: 10.1086/527414 .

Temmer, M., Rollett, T., M¨stl, C., Veronig, A.M., Vrˇ snak, B., Odstrˇil, D.: 2011, Influence o c of the Ambient Solar Wind Flow on the Propagation Behavior of Interplanetary Coronal Mass Ejections. Astrophys. J. 743 , 101. doi: 10.1088/0004-637X/743/2/101 .

Temmer, M., Vrˇ snak, B., Rollett, T., Bein, B., de Koning, C.A., Liu, Y., Bosman, E., Davies, J.A., M¨stl, C., ˇ ic, T., Veronig, A.M., Bothmer, V., Harrison, R., Nitta, N., Bisi, M., Flor, o Z O., Eastwood, J., Odstrcil, D., Forsyth, R.: 2012, Characteristics of Kinematics of a Coronal Mass Ejection during the 2010 August 1 CME-CME Interaction Event. Astrophys. J. 749 , 57. doi: 10.1088/0004-637X/749/1/57 .

Temmer, M., Veronig, A.M., Peinhart, V., Vrˇ snak, B.: 2014, Asymmetry in the CME-CME Interaction Process for the Events from 2011 February 14-15. Astrophys. J. 785 , 85. doi: 10.1088/0004-637X/785/2/85 .

Thernisien, A.: 2011, Implementation of the Graduated Cylindrical Shell Model for the Threedimensional Reconstruction of Coronal Mass Ejections. ApJS 194 , 33. doi: 10.1088/00670049/194/2/33 .

Thernisien, A., Vourlidas, A., Howard, R.A.: 2009, Forward Modeling of Coronal Mass Ejections Using STEREO/SECCHI Data. Solar Phys. 256 , 111 - 130. doi: 10.1007/s11207009-9346-5 .

Thompson, W.T.: 2009, 3D triangulation of a Sun-grazing comet. Icarus 200 , 351 - 357. doi: 10.1016/j.icarus.2008.12.011 .

Troshichev, O.A., Lukianova, R.Y., Papitashvili, V.O., Rich, F.J., Rasmussen, O.: 2000, Polar cap index (PC) as a proxy for ionospheric electric field in the near-pole region. Geophys. Res. Lett. 27 , 3809 - 3812. doi: 10.1029/2000GL003756 .

Tsurutani, B.T., Smith, E.J., Gonzalez, W.D., Tang, F., Akasofu, S.I.: 1988, Origin of interplanetary southward magnetic fields responsible for major magnetic storms near solar maximum (1978-1979). J. Geophys. Res. 93 , 8519 - 8531. doi: 10.1029/JA093iA08p08519 .

Tsurutani, B.T., Gonzalez, W.D., Gonzalez, A.L.C., Guarnieri, F.L., Gopalswamy, N., Grande, M., Kamide, Y., Kasahara, Y., Lu, G., Mann, I., McPherron, R., Soraas, F., Vasyliunas, V.: 2006, Corotating solar wind streams and recurrent geomagnetic activity: A review. Journal of Geophysical Research (Space Physics) 111 , 7. doi: 10.1029/2005JA011273 .

Vandas, M., Odstrcil, D.: 2004, Acceleration of electrons by interacting CMEs. Astron. Astrophys. 415 , 755 - 761. doi: 10.1051/0004-6361:20031763 .

Vandas, M., Fischer, S., Dryer, M., Smith, Z., Detman, T., Geranios, A.: 1997, MHD simulation of an interaction of a shock wave with a magnetic cloud. J. Geophys. Res. 102 , 22295 22300. doi: 10.1029/97JA01675 .

Vourlidas, A., Subramanian, P., Dere, K.P., Howard, R.A.: 2000, Large-Angle Spectrometric Coronagraph Measurements of the Energetics of Coronal Mass Ejections. Astrophys. J. 534 , 456 - 467. doi: 10.1086/308747 .

Vrˇ snak, B., ˇ ic, Z T., Falkenberg, T.V., M¨stl, C., Vennerstrom, S., Vrbanec, D.: 2010, The o role of aerodynamic drag in propagation of interplanetary coronal mass ejections. Astron. Astrophys. 512 , A43. doi: 10.1051/0004-6361/200913482 .

Vrˇ snak, B., Zic, T., Vrbanec, D., Temmer, M., Rollett, T., M¨stl, C., Veronig, A., Calogovi´, ˇ o ˇ c J., Dumbovi´ c, M., Luli´, S., Moon, Y.-J., Shanmugaraju, A.: 2013, Propagation of Interc planetary Coronal Mass Ejections: The Drag-Based Model. Solar Phys. 285 , 295 - 315. doi: 10.1007/s11207-012-0035-4 .

Wang, Y.M., Ye, P.Z., Wang, S.: 2003, Multiple magnetic clouds: Several examples during March-April 2001. Journal of Geophysical Research (Space Physics) 108 , 1370. doi: 10.1029/2003JA009850 .

Wang, Y.M., Ye, P.Z., Wang, S., Xue, X.H.: 2003, An interplanetary cause of large geomagnetic storms: Fast forward shock overtaking preceding magnetic cloud. Geophys. Res. Lett. 30 , 1700. doi: 10.1029/2002GL016861 .

Wang, Y., Zheng, H., Wang, S., Ye, P.: 2005, MHD simulation of the formation and propagation of multiple magnetic clouds in the heliosphere. Astron. Astrophys. 434 , 309 - 316. doi: 10.1051/0004-6361:20041423 .

Webb, D.F., M¨ ostl, C., Jackson, B.V., Bisi, M.M., Howard, T.A., Mulligan, T., Jensen, E.A., Jian, L.K., Davies, J.A., de Koning, C.A., Liu, Y., Temmer, M., Clover, J.M., Farrugia, C.J., Harrison, R.A., Nitta, N., Odstrcil, D., Tappin, S.J., Yu, H.-S.: 2013, Heliospheric Imaging of 3D Density Structures During the Multiple Coronal Mass Ejections of Late July to Early August 2010. Solar Phys. 285 , 317 - 348. doi: 10.1007/s11207-013-0260-5 .

Xiong, M., Zheng, H., Wang, S.: 2009, Magnetohydrodynamic simulation of the interaction between two interplanetary magnetic clouds and its consequent geoeffectiveness: 2. Oblique collision. Journal of Geophysical Research (Space Physics) 114 , 11101. doi: 10.1029/2009JA014079 .

Xiong, M., Zheng, H., Wang, Y., Wang, S.: 2006, Magnetohydrodynamic simulation of the interaction between interplanetary strong shock and magnetic cloud and its consequent geoeffectiveness: 2. Oblique collision. Journal of Geophysical Research (Space Physics) 111 , 11102. doi: 10.1029/2006JA011901 .

Xiong, M., Zheng, H., Wu, S.T., Wang, Y., Wang, S.: 2007, Magnetohydrodynamic simulation of the interaction between two interplanetary magnetic clouds and its consequent geoeffectiveness. Journal of Geophysical Research (Space Physics) 112 , 11103. doi: 10.1029/2007JA012320 .

Zurbuchen, T.H., Richardson, I.G.: 2006, In-Situ Solar Wind and Magnetic Field Signatures of Interplanetary Coronal Mass Ejections. Space Sci. Rev. 123 , 31 - 43. doi: 10.1007/s11214006-9010-4 .
Wageesh Mishra Nandita Srivastava D. Chakrabarty
Evolution and Consequences of Interacting CMEs of 2012 November 9-10 using STEREO/SECCHI and In Situ Observations
Understanding of the kinematic evolution of Coronal Mass Ejections (CMEs) in
the heliosphere is important to estimate their arrival time at the Earth. It is
found that kinematics of CMEs can change when they interact or collide with
each other as they propagate in the heliosphere. In this paper, we analyze the
collision and post-interaction characteristics of two Earth-directed CMEs,
launched successively on 2012 November 9 and 10, using white light imaging
observations from STEREO/SECCHI and in situ observations taken from WIND
spacecraft. We tracked two density enhanced features associated with leading
and trailing edge of November 9 CME and one density enhanced feature associated
with leading edge of November 10 CME by constructing J-maps. We found that the
leading edge of November 10 CME interacted with the trailing edge of November 9
CME. We also estimated the kinematics of these features of the CMEs and found a
significant change in their dynamics after interaction. In in situ
observations, we identified distinct structures associated with interacted CMEs
and also noticed their heating and compression as signatures of CME-CME
interaction. Our analysis shows an improvement in arrival time prediction of
CMEs using their post-collision dynamics than using pre-collision dynamics.
Estimating the true masses and speeds of these colliding CMEs, we investigated
the nature of observed collision which is found to be close to perfectly
inelastic. The investigation also places in perspective the geomagnetic
consequences of the two CMEs and their interaction in terms of occurrence of
geomagnetic storm and triggering of magnetospheric substorms.