Testing the asymmetric cone model for halo CMEs using STEREO/SECCHI coronagraphic observations

Space weather is significantly controlled by halo coronal mass ejections (HCMEs) originating close to the central meridian and directing toward the Earth. Unfortunately, coronagraphic observations (especially for HCMEs) are subject to a projection effect which makes it impossible to determine the true radial velocity and width of CMEs. However, these parameters can be estimated by correcting for the projection effect using the asymmetric cone model (Michalek, 2006). A set of 20 CMEs, observed as halo events in the LASCO field of view and simultaneously as limb events in the STEREO/SECCHI field of view, are used to check the accuracy of the asymmetric cone model. For this purpose, characteristics of the considered CMEs (angular widths and radial speeds) measured in STEREO/SECCHI images are compared with those obtained by the asymmetric cone model. We demonstrate that the widths and speeds determined by both methods are very similar. Correlation coefficients for speeds and angular widths are 0.99 and 0.96, respectively. We have also shown that the projection effect is unpredictable and could sometimes be very significant (up to 100% of the velocity measured in the LASCO field of view). On average, the SOHO/LASCO projected speeds for the HCMEs are 23% smaller than the radial velocities obtained from the STEREO/SECCHI images.     

Solar and interplanetary parameters of CMEs with and without type II radio bursts

We have analyzed 101 CMEs, and their associated ICMEs and interplanetary (IP) shocks observed during the period 1997–2005. The main aim of the present work is to study the interplanetary characteristics of metric and DH type II associated CMEs such as, shock strength, IP shock speed, ICME speed, stand off distance and transit time. Among these 101 CMEs, 38 events show both metric and DH type II bursts characteristics. There are no metric and DH type II association for 52 events. While DH type II alone is found in 7 cases, metric type II alone is found in 4 events. It is found that the mean speeds of CMEs increase progressively from CMEs without type II events to CMEs associated with metric and DH type IIs as suggested by Gopalswamy et al. (2005). In addition, we found that the speeds of ICMEs and IP shocks progressively increase in the following order: events without metric and DH type IIs, events with metric alone, events with DH alone and events with both metric and DH type IIs. Similarly the Mach number is found to increase in the same order. While there is not much change in the stand-off distance among these cases, it is minimum (∼18 R_⊙) for CMEs with speed greater than 2200 km/s. The above results confirm that more energetic CMEs can produce both metric and DH type IIs for which the interplanetary parameters such as mean values of ICME speed and IP shock speed and Mach number are found to be higher.     

Investigating the Coronal Mass Ejections associated with DH type-II radio bursts and solar flares during the ascending phase of the solar cycle 24

We studied a set of 74 CMEs, with shedding the light on the halo-CMEs (HCMEs), that are associated with decametric – hectometric (DH) type-II radio bursts (1–16?MHz) and solar flares during the period 2008–2014. The events were classified into 3 groups (disk, intermediate, and limb events) based on their longitudinal distribution.We found that the events are mostly distributed around 15.32° and 15.97° at the northern and southern solar hemispheres, respectively. We found that there is a clear dependence between the longitude and the CME’s width, speed, acceleration, mass, and kinetic energy. For the CMEs’ widths, most of the events were HCMEs (～62%), while the partial HCMEs comprised ～35% and the rest of events were CMEs with widths less than 120°. For the CMEs’ speeds, masses, and kinetic energies, the mean values showed a direct proportionality with the longitude, in which the limb events had the highest speeds, the largest masses, and the highest kinetic energies. The mean peak flux of the solar flares for different longitudes was comparable, but the disk flares were more energetic. The intermediate flares were considered as gradual flares since they tended to last longer, while the limb flares were considered as impulsive flares since they tended to last shorter.A weak correlation (R?=?0.32) between the kinetic energy of the CMEs and the duration of the associated flares has been noticed, while there was a good correlation (R?=?0.76) between the kinetic energy of the CMEs and the peak flux of the associated flares. We found a fair correlation (R?=?0.58) between the kinetic energy of the CMEs and the duration of the associated DH type-II radio bursts.     

Height of shock formation in the solar corona inferred from observations of type II radio bursts and coronal mass ejections

Employing coronagraphic and EUV observations close to the solar surface made by the Solar Terrestrial Relations Observatory (STEREO) mission, we determined the heliocentric distance of coronal mass ejections (CMEs) at the starting time of associated metric type II bursts. We used the wave diameter and leading edge methods and measured the CME heights for a set of 32 metric type II bursts from solar cycle 24. We minimized the projection effects by making the measurements from a view that is roughly orthogonal to the direction of the ejection. We also chose image frames close to the onset times of the type II bursts, so no extrapolation was necessary. We found that the CMEs were located in the heliocentric distance range from 1.20 to 1.93 solar radii (Rs), with mean and median values of 1.43 and 1.38 Rs, respectively. We conclusively find that the shock formation can occur at heights substantially below 1.5 Rs. In a few cases, the CME height at type II onset was close to 2 Rs. In these cases, the starting frequency of the type II bursts was very low, in the range 25–40 MHz, which confirms that the shock can also form at larger heights. The starting frequencies of metric type II bursts have a weak correlation with the measured CME/shock heights and are consistent with the rapid decline of density with height in the inner corona.     

Modeling of “gradual” solar energeticparticle events using a stochastic differential equation method

We have modeled “gradual” solar energetic particle events through numerical simulations using a StochasticDifferential Equation (SDE) method. We consider that energetic particle events are roughly divided into two groups: (1) where the shock was driven by coronal mass ejections (CMEs) associated with large solar flares, and (2) where they have no related solar events apart from the CMEs. (The detailed classification of energetic particle events was discussed in our previous paper.) What we call “gradual” solar energetic particle events belong to the former group. Particles with energies greater than 10 MeV are observed within several hours after the occurrence of flares and CMEs in many gradual events. By applying the SDE method coupled with particle splitting to diffusive acceleration, we found that an injection of high energy particles is necessary for early enhancement of such a high-energy proton flux and that it should not be presumed that the solar wind particles act as the seed population.     

Properties and relationship between solar eruptive flares and Coronal Mass Ejections during rising phase of Solar Cycles 23 and 24

Statistical relationship between major flares and the associated CMEs during rising phases of Solar Cycles 23 and 24 are studied. Totally more than 6000 and 10,000 CMEs were observed by SOHO/LASCO (Solar and Heliospheric Observatory/Large Angle Spectrometric Coronagraph) during 23rd [May 1996–June 2002] and 24th [December 2008–December 2014] solar cycles, respectively. In particular, we studied the relationship between properties of flares and CMEs using the limb events (longitude 70–85°) to avoid projection effects of CMEs and partial occultation of flares that occurred near 90°. After selecting a sample of limb flares, we used certain spatial and temporal constraints to find the flare-CME pairs. Using these constraints, we compiled 129 events in Solar Cycle 23 and 92 events in Solar Cycle 24. We compared the flare-CME relationship in the two solar cycles and no significant differences are found between the two cycles. We only found out that the CME mean width was slightly larger and the CME mean acceleration was slightly higher in cycle 24, and that there was somewhat a better relation between flare flux and CME deceleration in cycle 24 than in cycle 23.     

Remote sensing of the ring current using energetic neutral atoms

Energetic neutral atom (ENA) images of the storm-time ring current obtained from the ISEE-1 spacecraft provide information for a “zero-order” global model of the energetic ion distribution. With the assumption of isotropic pressure and magnetostatic, non-convective pressure balance, the global system of electrical currents driven by the ion pressure can be calculated using Euler potentials for the divergenceless current density. Radial pressure gradients drive azimuthal currents, and azimuthal pressure gradients drive radial currents. The radial currents cause current lines in the inner magnetosphere to close in the ionosphere, forming a “partial” ring current. The intensities and locations of these field-aligned currents driven into and out of the ionosphere resemble those of the observed Region 2 current system, but not all observed properties of the Region 2 system are reproduced by the “zero-order” model.     

CME dynamics using coronagraph and interplanetary ejecta data

In this work, we present a study of the coronal mass ejection (CME) dynamics using LASCO coronagraph observations combined with in-situ ACE plasma and magnetic field data, covering a continuous period of time from January 1997 to April 2001, complemented by few extreme events observed in 2001 and 2003. We show, for the first time, that the CME expansion speed correlates very well with the travel time to 1 AU of the interplanetary ejecta (or ICMEs) associated with the CMEs, as well as with their preceding shocks. The events analyzed in this work are a subset of the events studied in Schwenn et al. (2005), from which only the CMEs associated with interplanetary ejecta (ICMEs) were selected. Three models to predict CME travel time to Earth, two proposed by Gopalswamy et al. (2001) and one by Schwenn et al. (2005), were used to characterize the dynamical behavior of this set of events. Extreme events occurred in 2001 and 2003 were used to test the prediction capability of the models regarding CMEs with very high LASCO C3 speeds.     

Study of CME transit speeds for the event of 07-NOV-2004

Several methods for CME speed estimation are discussed. These include velocity derivation based on the frequency drifts observed in metric and decametric radio wave data using a range of coronal density models. Coronagraph height–time plots allow measurement of plane-of-sky and expansion speeds. These in turn can enable propagation speeds to be derived from a range of empirical relations. Simple geometric e.g., cone, models can provide propagation velocity estimates for suitable halo or partial halo events. Interplanetary scintillation observations allow speed estimates at large distances from the Sun detecting in particular the deceleration of the faster CMEs. Related interplanetary shocks and the arrival times and speeds of the associated magnetic clouds at Earth can also be considered. We discuss the application of some of these methods to the transit to Earth of a complex CME that originated earlier than 16:54 U.T. on 07-NOV-2004. The difficulties in making velocity estimates from radio observations, particularly under disturbed coronal conditions, are highlighted.     

A statistical study of post-flare-associated CME events

We present a statistical study of post-flare-associated CMEs (PFA-CMEs) during the period from 1996 to 2010. By investigating all CMEs and X-ray flares, respectively, in the LASCO and GOES archives, we found 15875 CMEs of which masses are well measured and 25112 X-ray flares of which positions are determined from their optical counterparts. Under certain temporal and spatial criteria of these CMEs and solar flare events, 291PFA-CMEs events have been selected. Linking the flare fluxes with CME speeds of these paired events, we found that there is a reasonable positive linear relation between the CME linear speed and associated flare flux. The results show also the CME width increases as the flux of its associated solar flare increases. Besides we found that there is a fine positive linear relation between the CME mass and its width. Matching the flare fluxes with CME masses of these paired events, we find the CME mass increases as the flux of its associated solar flare increases. Finally we find the PFA-CME events are in regular more decelerated than the other CMEs.     

Properties of narrow coronal massejections observed with LASCO

We report the statistical properties of narrow coronal mass ejections (CMEs, angular width < 20°) withparticular emphasis on comparison with normal CMEs. We investigated 806 narrow CMEs from an online LASCO/CME catalog and found that (1) the fraction of narrow CMEs increases from 12% to 22% towards solar maximum, (2) during the solar maximum, the narrow CMEs are generally faster than normal ones, (3) the maximum speed of narrow CMEs (1141 km s⁻¹) is much smaller than that of the normal CMEs (2604 km s⁻¹). These results imply that narrow CMEs do not form a subset of normal CMEs and have a different acceleration mechanism from normal CMEs.     

Local-area helioseismology as a diagnostic tool for solar variability

Dynamical and thermal variations of the internal structure of the Sun can affect the energy flow and result in variations in irradiance at the surface. Studying variations in the interior is crucial for understanding the mechanisms of the irradiance variations. “Global” helioseismology based on analysis of normal mode frequencies, has helped to reveal radial and latitudinal variations of the solar structure and dynamics associated with the solar cycle in the deep interior. A new technique, - “local-area” helioseismology or heliotomography, offers additional potentially important diagnostics by providing three-dimensional maps of the sound speed and flows in the upper convection zone. These diagnostics are based on inversion of travel times of acoustic waves which propagate between different points on the solar surface through the interior. The most significant variations in the thermodynamic structure found by this method are associated with sunspots and complexes of solar activity. The inversion results provide evidence for areas of higher sound speed beneath sunspot regions located at depths of 4–20 Mm, which may be due to accumulated heat or magnetic field concentrations. However, the physics of these structures is not yet understood. Heliotomography also provides information about large-scale stable longitudinal structures in the solar interior, which can be used in irradiance models. This new diagnostic tool for solar variability is currently under development. It will require both a substantial theoretical and modeling effort and high-resolution data to develop new capabilities for understanding mechanisms of solar variability.     

Evolution of Dst index,cosmic rays and global temperature during solar cycles 20–23

We have studied conditions in interplanetary space, which can have an influence on galactic cosmic ray (CR) and climate change. In this connection the solar wind and interplanetary magnetic field parameters and cosmic ray variations have been compared with geomagnetic activity represented by the equatorial Dst index from the beginning 1965 to the end of 2012. Dst index is commonly used as the solar wind–magnetosphere–ionosphere interaction characteristic. The important drivers in interplanetary medium which have effect on cosmic rays as CMEs (coronal mass ejections) and CIRs (corotating interaction regions) undergo very strong changes during their propagation to the Earth. Because of this CMEs, coronal holes and the solar spot numbers (SSN) do not adequately reflect peculiarities concerned with the solar wind arrival to 1 AU. Therefore, the geomagnetic indices have some inestimable advantage as continuous series other the irregular solar wind measurements. We have compared the yearly average variations of Dst index and the solar wind parameters with cosmic ray data from Moscow, Climax, and Haleakala neutron monitors during the solar cycles 20–23. The descending phases of these solar cycles (CSs) had the long-lasting solar wind high speed streams occurred frequently and were the primary contributors to the recurrent Dst variations. They also had effects on cosmic rays variations. We show that long-term Dst variations in these solar cycles were correlated with the cosmic ray count rate and can be used for study of CR variations. Global temperature variations in connection with evolution of Dst index and CR variations is discussed.     

CMEs and active regions on the sun

The relationship between active regions (ARs) and coronal mass ejections (CMEs) is studied. For this purpose a statistical analysis of 694 CMEs associated with ARs was carried out. We considered the relationship between properties of the CMEs and ARs characterized using the McIntosh classification. We demonstrated that CMEs are likely to be launched from ARs in the mature phase of their evolution when they have complex magnetic field. The fastest and halo CMEs can be ejected only from the most complex ARs (when an AR is a bipolar group of spots with large asymmetric penumbras around the main spot with many smaller spots in the group). We also showed that the wider events have a tendency to originate from uncomplicated magnetic structures. This tendency was used for estimation of the real angular widths of the halo CMEs. The probability of launching of fast CMEs increases together with increase of the complexity and size of ARs. The widest, but slow, CMEs originate from the simplest magnetic structure which are still able to produce CMEs. Our results could be useful for forecasting of space weather.     

日冕物质抛射的磁绳正演法及源区投影修正方法的讨论

利用磁绳正演法(forward modeling)研究了2007---2008年STEREO卫星观测到的21个结构清晰的日冕物质抛射(CME)事件, 得到这些CME事件在三维空间的运动速度; 与根据STEREO A星和B星观测速度进行源区投影修正得到的结果进行比较, 发现这种源区投影修正方法存在很大的局限性. 统计结果表明, 当CME源区在日面上的日心角距(CME源区和日心连线与观测点和日心连线的夹角)大于50o时, 修正速度与三维速度之间的差别不大; 当日心角距小于这个值时, 修正速度与三维速度差别明显, 尤其对于小日心角距, 相差甚大. 统计结果还表明, 当日心角距大于65o时, A星和B星得到的投影速度与三维速度相近, 投影速度可以作为CME的三维速度.      

Determining CME parameters by fitting heliospheric observations: Numerical investigation of the accuracy of the methods

Transients in the heliosphere, including coronal mass ejections (CMEs) and corotating interaction regions can be imaged to large heliocentric distances by heliospheric imagers (HIs), such as the HIs onboard STEREO and SMEI onboard Coriolis. These observations can be analyzed using different techniques to derive the CME speed and direction. In this paper, we use a three-dimensional (3-D) magneto-hydrodynamic (MHD) numerical simulation to investigate one of these methods, the fitting method of  and . Because we use a 3-D simulation, we can determine with great accuracy the CME initial speed, its speed at 1 AU and its average transit speed as well as its size and direction of propagation. We are able to compare the results of the fitting method with the values from the simulation for different viewing angles between the CME direction of propagation and the Sun-spacecraft line. We focus on one simulation of a wide (120–140°) CME, whose initial speed is about 800 km s⁻¹. For this case, we find that the best-fit speed is in good agreement with the speed of the CME at 1 AU, and this, independently of the viewing angle. The fitted direction of propagation is not in good agreement with the viewing angle in the simulation, although smaller viewing angles result in smaller fitted directions. This is due to the extremely wide nature of the ejection. A new fitting method, proposed to take into account the CME width, results in better agreement between fitted and actual directions for directions close to the Sun–Earth line. For other directions, it gives results comparable to the fitting method of Sheeley et al. (1999). The CME deceleration has only a small effect on the fitted direction, resulting in fitted values about 1–4° higher than the actual values.     

A study on radio-loud interacting/non-interacting CMEs-associated SEPs and solar flares

We have established a data set of 58 major hybrid SEP events associated with meter-to-decahectometer wavelength (m-to-DH) type II bursts, solar flares, and radio-load CMEs during the period of 1997–2014. The main focus of our study is to address the following two questions: Does the interaction of CMEs play a role in the enhancement of SEP intensity? Is there any difference in the seed population, and parent eruptions in the SEP events with and without CME interactions? Hence, the sample of 58 events is classified into two sets: (i) 35 non-interacting-CME-associated SEP events; (ii) 23 interacting-CME-associated SEP events. All the characteristics of SEPs, their associated CMEs/flares and the relationships between them are statistically analyzed and compared. Some of the basic attributes and relative elemental abundances (Fe/O ratios) of the both the sets are also compared. The results indicate that the seed particles in non-interacting-CME-associated SEP events are mostly from solar wind/coronal materials. But in the case of interacting-CME-associated SEP events, it may be associated with both flare material from preceding flares and coronal materials from solar wind/preceding CMEs. The correlation studies reveal that there are clear correlations between logarithmic peak intensity of SEP events and properties of CMEs (space speed: cc?=?0.56) and solar flares (peak intensity: cc?=?0.40; integrated flux: cc?=?0.52) for non-interacting-CME-associated SEP events. But these correlations are absent for the interacting-CME-associated events. In addition, the results suggest that interaction of primary CMEs with their preceding CMEs plays an important role in the enhancement of peak intensity of SEPs at least for a set of m-to-DH type II bursts associated SEP events.     

Coronal mass ejection speeds measured in the solar corona using LASCO C2 and C3 images

In this work we present height-time diagrams of 2 halo coronal mass ejections, observed on September 28th, 1997 and June 29th, 1999. The CMEs were observed by the Large Angle and Spectroscopic Coronagraph (LASCO), which observes the solar corona from 2 to 32 solar radii. To obtain these diagrams we divide the LASCO images of a given sequence in angular slices, transform them into rectangular slices (their width chosen proportional to the time distance to the next image) and place them side by side. Thus, the speed profile of any pattern moving in the particular latitudinal slice can be derived. With this method we were able to identify even minor speed changes in several angular positions for the chosen events. This technique is particularly appropriate to identify acceleration or deceleration of structures in halo CMEs.     

Acceleration and transport of energeticparticles at CME-driven shocks

Historically, solar energetic particle (SEP) events are classified in two classes as “impulsive” and “gradual”. Whether there is a clear distinction between the two classes is still a matter of debate, but it is now commonly accepted that in large “gradual” SEP events, Fermi acceleration, also known as diffusive shock acceleration, is the underlying acceleration mechanism. At shock waves driven by coronal mass ejections (CMEs), particles are accelerated diffusively at the shock and often reach > MeV energies (and perhaps up to GeV energies). As a CME-driven shock propagates, expands and weakens, the accelerated particles can escape ahead of the shock into the interplanetary medium. These escaping energized particles then propagate along the interplanetary magnetic field, experiencing only weak scattering from fluctuations in the interplanetary magnetic field (IMF). In this paper, we use a Monte-Carlo approach to study the transport of energetic particles escaping from a CME-driven shock. We present particle spectra observed at 1 AU. We also discuss the particle “crossing number” at 1AU and its implication to particle anisotropy. Based on previous models of particle acceleration at CME-driven shocks, our simulation allows us to investigate various characteristics of energetic particles arriving at various distances from the sun. This provides us an excellent basis for understanding the observations of high-energy particles made at 1 AU by ACE and WIND.