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

Examining Periodic Solar Wind Density Structures Observed in the SECCHI Heliospheric Imagers

We present an analysis of small-scale, periodic, solar-wind density enhancements (length-scales as small as \approx 1000 Mm) observed in images from the Heliospheric Imager (HI) aboard STEREO A. We discuss their possible relationship to periodic fluctuations of the proton density that have been identified at 1 AU using in-situ plasma measurements. Specifically, Viall, Kepko, and Spence (2008) examined 11 years of in-situ solar-wind density measurements at 1 AU and demonstrated that not only turbulent structures, but also non-turbulent periodic density structures exist in the solar wind with scale sizes of hundreds to one thousand Mm. In a subsequent paper, Viall, Spence, and Kasper (2009) analyzed the {\alpha} to proton solar-wind abundance ratio measured during one such event of periodic density structures, demonstrating that the plasma behavior was highly suggestive that either temporally or spatially varying coronal source plasma created those density structures. Large periodic density structures observed at 1 AU, which were generated in the corona, can be observable in coronal and heliospheric white-light images if they possess sufficiently high density contrast. Indeed, we identify such periodic density structures as they enter the HI field of view and follow them as they advect with the solar wind through the images. The smaller periodic density structures that we identify in the images are comparable in size to the larger structures analyzed in-situ at 1 AU, yielding further evidence that periodic density enhancements are a consequence of coronal activity as the solar wind is formed.Comment: 15 pages, 12 figures. The final publication is available at http://www.springerlink.co

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

How Many CMEs Have Flux Ropes? Deciphering the Signatures of Shocks, Flux Ropes, and Prominences in Coronagraph Observations of CMEs

We intend to provide a comprehensive answer to the question on whether all Coronal Mass Ejections (CMEs) have flux rope structure. To achieve this, we present a synthesis of the LASCO CME observations over the last sixteen years, assisted by 3D MHD simulations of the breakout model, EUV and coronagraphic observations from STEREO and SDO, and statistics from a revised LASCO CME database. We argue that the bright loop often seen as the CME leading edge is the result of pileup at the boundary of the erupting flux rope irrespective of whether a cavity or, more generally, a 3-part CME can be identified. Based on our previous work on white light shock detection and supported by the MHD simulations, we identify a new type of morphology, the `two-front' morphology. It consists of a faint front followed by diffuse emission and the bright loop-like CME leading edge. We show that the faint front is caused by density compression at a wave (or possibly shock) front driven by the CME. We also present high-detailed multi-wavelength EUV observations that clarify the relative positioning of the prominence at the bottom of a coronal cavity with clear flux rope structure. Finally, we visually check the full LASCO CME database for flux rope structures. In the process, we classify the events into two clear flux rope classes (`3-part', `Loop'), jets and outflows (no clear structure). We find that at least 40% of the observed CMEs have clear flux rope structures. We propose a new definition for flux rope CMEs (FR-CMEs) as a coherent magnetic, twist-carrying coronal structure with angular width of at least 40 deg and able to reach beyond 10 Rsun which erupts on a time scale of a few minutes to several hours. We conclude that flux ropes are a common occurrence in CMEs and pose a challenge for future studies to identify CMEs that are clearly not FR-CMEs.Comment: 26 pages, 9 figs, to be published in Solar Physics Topical Issue "Flux Rope Structure of CMEs

The Impact of New EUV Diagnostics on CME-Related Kinematics

We present the application of novel diagnostics to the spectroscopic observation of a Coronal Mass Ejection (CME) on disk by the Extreme Ultraviolet Imaging Spectrometer (EIS) on the Hinode spacecraft. We apply a recently developed line profile asymmetry analysis to the spectroscopic observation of NOAA AR 10930 on 14-15 December 2006 to three raster observations before and during the eruption of a 1000km/s CME. We see the impact that the observer's line-of-sight and magnetic field geometry have on the diagnostics used. Further, and more importantly, we identify the on-disk signature of a high-speed outflow behind the CME in the dimming region arising as a result of the eruption. Supported by recent coronal observations of the STEREO spacecraft, we speculate about the momentum flux resulting from this outflow as a secondary momentum source to the CME. The results presented highlight the importance of spectroscopic measurements in relation to CME kinematics, and the need for full-disk synoptic spectroscopic observations of the coronal and chromospheric plasmas to capture the signature of such explosive energy release as a way of providing better constraints of CME propagation times to L1, or any other point of interest in the heliosphere.Comment: Accepted to appear in Solar Physics Topical Issue titled "Remote Sensing of the Inner Heliosphere". Manuscript has 14 pages, 5 color figures. Movies supporting the figures can be found in http://download.hao.ucar.edu/pub/mscott/papers/Weathe

Solar Stereoscopy with STEREO/EUVI A and B spacecraft from small (6 deg) to large (170 deg) spacecraft separation angles

We performed for the first time stereoscopic triangulation of coronal loops in active regions over the entire range of spacecraft separation angles ($\alpha_{sep}\approx 6^\circ, 43^\circ, 89^\circ, 127^\circ$, and $170^\circ$). The accuracy of stereoscopic correlation depends mostly on the viewing angle with respect to the solar surface for each spacecraft, which affects the stereoscopic correspondence identification of loops in image pairs. From a simple theoretical model we predict an optimum range of $\alpha_{sep} \approx 22^\circ-125^\circ$, which is also experimentally confirmed. The best accuracy is generally obtained when an active region passes the central meridian (viewed from Earth), which yields a symmetric view for both STEREO spacecraft and causes minimum horizontal foreshortening. For the extended angular range of $\alpha_{sep}\approx 6^\circ-127^{\circ}$ we find a mean 3D misalignment angle of $\mu_{PF} \approx 21^\circ-39^\circ$ of stereoscopically triangulated loops with magnetic potential field models, and $\mu_{FFF} \approx 15^\circ-21^\circ$ for a force-free field model, which is partly caused by stereoscopic uncertainties $\mu_{SE} \approx 9^\circ$. We predict optimum conditions for solar stereoscopy during the time intervals of 2012--2014, 2016--2017, and 2021--2023.Comment: Solar Physics, (in press), 22 pages, 9 figure

Automated Detection of Coronal Loops using a Wavelet Transform Modulus Maxima Method

We propose and test a wavelet transform modulus maxima method for the au- tomated detection and extraction of coronal loops in extreme ultraviolet images of the solar corona. This method decomposes an image into a number of size scales and tracks enhanced power along each ridge corresponding to a coronal loop at each scale. We compare the results across scales and suggest the optimum set of parameters to maximise completeness while minimising detection of noise. For a test coronal image, we compare the global statistics (e.g., number of loops at each length) to previous automated coronal-loop detection algorithms

Geomagnetic storm dependence on the solar flare class

Content. Solar flares are often used as precursors of geomagnetic storms. In particular, Howard and Tappin (2005) recently published in A&A a dependence between X-ray class of solar flares and Ap and Dst indexes of geomagnetic storms which contradicts to early published results. Aims. We compare published results on flare-storm dependences and discuss possible sources of the discrepancy. Methods. We analyze following sources of difference: (1) different intervals of observations, (2) different statistics and (3) different methods of event identification and comparison. Results. Our analysis shows that magnitude of geomagnetic storms is likely to be independent on X-ray class of solar flares.Comment: 3 pages, 1 tabl

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