Accuracy and Limitations of Fitting and Stereoscopic Methods to Determine the Direction of Coronal Mass Ejections from Heliospheric Imagers Observations

Using data from the Heliospheric Imagers (HIs) onboard STEREO, it is possible to derive the direction of propagation of coronal mass ejections (CMEs) in addition to their speed with a variety of methods. For CMEs observed by both STEREO spacecraft, it is possible to derive their direction using simultaneous observations from the twin spacecraft and also, using observations from only one spacecraft with fitting methods. This makes it possible to test and compare different analyses techniques. In this article, we propose a new fitting method based on observations from one spacecraft, which we compare to the commonly used fitting method of Sheeley et al. (1999). We also compare the results from these two fitting methods with those from two stereoscopic methods, focusing on 12 CMEs observed simultaneously by the two STEREO spacecraft in 2008 and 2009. We find evidence that the fitting method of Sheeley et al. (1999) can result in significant errors in the determination of the CME direction when the CME propagates outside of 60deg \pm 20 deg from the Sun-spacecraft line. We expect our new fitting method to be better adapted to the analysis of halo or limb CMEs with respect to the observing spacecraft. We also find some evidence that direct triangulation in the HI fields-of-view should only be applied to CMEs propagating approximatively towards Earth (\pm 20deg from the Sun-Earth line). Last, we address one of the possible sources of errors of fitting methods: the assumption of radial propagation. Using stereoscopic methods, we find that at least seven of the 12 studied CMEs had an heliospheric deflection of less than 20deg as they propagated in the HI fields-of-view, which, we believe, validates this approximation.Comment: 17 pages, 6 figures, 2 tables, accepted to Solar Physic

Using an Ellipsoid Model to Track and Predict the Evolution and Propagation of Coronal Mass Ejections

We present a method for tracking and predicting the propagation and evolution of coronal mass ejections (CMEs) using the imagers on the STEREO and SOHO satellites. By empirically modeling the material between the inner core and leading edge of a CME as an expanding, outward propagating ellipsoid, we track its evolution in three-dimensional space. Though more complex empirical CME models have been developed, we examine the accuracy of this relatively simple geometric model, which incorporates relatively few physical assumptions, including i) a constant propagation angle and ii) an azimuthally symmetric structure. Testing our ellipsoid model developed herein on three separate CMEs, we find that it is an effective tool for predicting the arrival of density enhancements and the duration of each event near 1 AU. For each CME studied, the trends in the trajectory, as well as the radial and transverse expansion are studied from 0 to ~.3 AU to create predictions at 1 AU with an average accuracy of 2.9 hours.Comment: 18 pages, 11 figure

Reconstructing the 3-D Trajectories of CMEs in the Inner Heliosphere

A method for the full three-dimensional (3-D) reconstruction of the trajectories of coronal mass ejections (CMEs) using Solar TErrestrial RElations Observatory (STEREO) data is presented. Four CMEs that were simultaneously observed by the inner and outer coronagraphs (COR1 and 2) of the Ahead and Behind STEREO satellites were analysed. These observations were used to derive CME trajectories in 3-D out to ~15Rsun. The reconstructions using COR1/2 data support a radial propagation model. Assuming pseudo-radial propagation at large distances from the Sun (15-240Rsun), the CME positions were extrapolated into the Heliospheric Imager (HI) field-of-view. We estimated the CME velocities in the different fields-of-view. It was found that CMEs slower than the solar wind were accelerated, while CMEs faster than the solar wind were decelerated, with both tending to the solar wind velocity.Comment: 17 pages, 10 figures, 1 appendi

Effect of Solar Wind Drag on the Determination of the Properties of Coronal Mass Ejections from Heliospheric Images

The Fixed-\Phi (F\Phi) and Harmonic Mean (HM) fitting methods are two methods to determine the average direction and velocity of coronal mass ejections (CMEs) from time-elongation tracks produced by Heliospheric Imagers (HIs), such as the HIs onboard the STEREO spacecraft. Both methods assume a constant velocity in their descriptions of the time-elongation profiles of CMEs, which are used to fit the observed time-elongation data. Here, we analyze the effect of aerodynamic drag on CMEs propagating through interplanetary space, and how this drag affects the result of the F\Phi and HM fitting methods. A simple drag model is used to analytically construct time-elongation profiles which are then fitted with the two methods. It is found that higher angles and velocities give rise to greater error in both methods, reaching errors in the direction of propagation of up to 15 deg and 30 deg for the F\Phi and HM fitting methods, respectively. This is due to the physical accelerations of the CMEs being interpreted as geometrical accelerations by the fitting methods. Because of the geometrical definition of the HM fitting method, it is affected by the acceleration more greatly than the F\Phi fitting method. Overall, we find that both techniques overestimate the initial (and final) velocity and direction for fast CMEs propagating beyond 90 deg from the Sun-spacecraft line, meaning that arrival times at 1 AU would be predicted early (by up to 12 hours). We also find that the direction and arrival time of a wide and decelerating CME can be better reproduced by the F\Phi due to the cancellation of two errors: neglecting the CME width and neglecting the CME deceleration. Overall, the inaccuracies of the two fitting methods are expected to play an important role in the prediction of CME hit and arrival times as we head towards solar maximum and the STEREO spacecraft further move behind the Sun.Comment: Solar Physics, Online First, 17 page

Interplanetary and Geomagnetic Consequences of Interacting CMEs of 13-14 June 2012

We report on the kinematics of two interacting CMEs observed on 13 and 14 June 2012. Both CMEs originated from the same active region NOAA 11504. After their launches which were separated by several hours, they were observed to interact at a distance of 100 Rs from the Sun. The interaction led to a moderate geomagnetic storm at the Earth with Dst index of approximately, -86 nT. The kinematics of the two CMEs is estimated using data from the Sun Earth Connection Coronal and Heliospheric Investigation (SECCHI) onboard the Solar Terrestrial Relations Observatory (STEREO). Assuming a head-on collision scenario, we find that the collision is inelastic in nature. Further, the signatures of their interaction are examined using the in situ observations obtained by Wind and the Advance Composition Explorer (ACE) spacecraft. It is also found that this interaction event led to the strongest sudden storm commencement (SSC) (approximately 150 nT) of the present Solar Cycle 24. The SSC was of long duration, approximately 20 hours. The role of interacting CMEs in enhancing the geoeffectiveness is examined.Comment: 17 pages, 5 figures, Accepted in Solar Physics Journa

3D Coronal Density Reconstruction and Retrieving the Magnetic Field Structure during Solar Minimum

Measurement of the coronal magnetic field is a crucial ingredient in understanding the nature of solar coronal phenomena at all scales. We employed STEREO/COR1 data obtained during a deep minimum of solar activity in February 2008 (Carrington rotation CR 2066) to retrieve and analyze the three-dimensional (3D) coronal electron density in the range of heights from 1.5 to 4 Rsun using a tomography method. With this, we qualitatively deduced structures of the coronal magnetic field. The 3D electron density analysis is complemented by the 3D STEREO/EUVI emissivity in the 195 A band obtained by tomography for the same CR. A global 3D MHD model of the solar corona was used to relate the reconstructed 3D density and emissivity to open/closed magnetic field structures. We show that the density maximum locations can serve as an indicator of current sheet position, while the locations of the density gradient maximum can be a reliable indicator of coronal hole boundaries. We find that the magnetic field configuration during CR 2066 has a tendency to become radially open at heliocentric distances greater than 2.5 Rsun. We also find that the potential field model with a fixed source surface (PFSS) is inconsistent with the boundaries between the regions with open and closed magnetic field structures. This indicates that the assumption of the potential nature of the coronal global magnetic field is not satisfied even during the deep solar minimum. Results of our 3D density reconstruction will help to constrain solar coronal field models and test the accuracy of the magnetic field approximations for coronal modeling.Comment: Published in "Solar Physics

The random phase property and the Lyapunov Spectrum for disordered multi-channel systems

A random phase property establishing in the weak coupling limit a link between quasi-one-dimensional random Schrödinger operators and full random matrix theory is advocated. Briefly summarized it states that the random transfer matrices placed into a normal system of coordinates act on the isotropic frames and lead to a Markov process with a unique invariant measure which is of geometric nature. On the elliptic part of the transfer matrices, this measure is invariant under the unitaries in the hermitian symplectic group of the universality class under study. While the random phase property can up to now only be proved in special models or in a restricted sense, we provide strong numerical evidence that it holds in the Anderson model of localization. A main outcome of the random phase property is a perturbative calculation of the Lyapunov exponents which shows that the Lyapunov spectrum is equidistant and that the localization lengths for large systems in the unitary, orthogonal and symplectic ensemble differ by a factor 2 each. In an Anderson-Ando model on a tubular geometry with magnetic field and spin-orbit coupling, the normal system of coordinates is calculated and this is used to derive explicit energy dependent formulas for the Lyapunov spectrum