Issue |
A&A
Volume 635, March 2020
|
|
---|---|---|
Article Number | A100 | |
Number of page(s) | 11 | |
Section | The Sun and the Heliosphere | |
DOI | https://doi.org/10.1051/0004-6361/201936664 | |
Published online | 16 March 2020 |
Asymmetric expansion of coronal mass ejections in the low corona
1
Universidad Tecnológica Nacional, Facultad Regional Mendoza, Centro de Estudios para el Desarrollo Sustentable, Rodriguez 243, M5502AJE Mendoza, Argentina
e-mail: hebe.cremades@frm.utn.edu.ar
2
Consejo Nacional de Investigaciones Científicas y Técnicas, Godoy Cruz 2290, C1425FQB Ciudad Autónoma de Buenos Aires, Argentina
Received:
10
September
2019
Accepted:
5
January
2020
Aims. Understanding how magnetic fields are structured within coronal mass ejections (CMEs), and how they evolve from the low corona into the heliosphere, is a major challenge for space weather forecasting and for solar physics. The study of CME morphology is a particularly auspicious approach to this problem, given that it holds a close relationship with the CME magnetic field configuration. Although earlier studies have suggested an asymmetry in the width of CMEs in orthogonal directions, this has not been inspected using multi-viewpoint observations.
Methods. The improved spatial, temporal, and spectral resolution, added to the multiple vantage points offered by missions of the Heliophysics System Observatory, constitute a unique opportunity to gain insight into this regard. We inspect the early evolution (below ten solar radii) of the morphology of a dozen CMEs occurring under specific conditions of observing spacecraft location and CME trajectory, favorable to reduce uncertainties typically involved in the 3D reconstruction used here. These events are carefully reconstructed by means of a forward modeling tool using simultaneous observations of the Solar-Terrestrial Relations Observatory (STEREO) Extreme Ultraviolet Imager and the Solar Dynamics Observatory Atmospheric Imaging Assembly as input when originating low in the corona, and followed up in the outer fields of view of the STEREO and the Solar and Heliospheric Observatory coronagraphs. We then examine the height evolution of the morphological parameters arising from the reconstructions.
Results. The multi-viewpoint analysis of this set of CMEs revealed that their initial expansion – below three solar radii – is considerably asymmetric and non-self-similar. Both angular widths, namely along the main axes of CMEs (AWL) and in the orthogonal direction (AWD, representative of the flux rope diameter), exhibit much steeper change rates below this height, with the growth rate of AWL found to be larger than that of AWD, also below that height. Angular widths along the main axes of CMEs are on average ≈1.8 times larger than widths in the orthogonal direction AWD. The ratios of the two expansion speeds, namely in the directions of CMEs main axes and in their orthogonal, are nearly constant in time after ∼4 solar radii, with an average ratio ≈1.6. Heights at which the width change rate is defined to stabilize are greater for AWL than for AWD.
Key words: Sun: coronal mass ejections (CMEs) / Sun: corona / Sun: activity
© ESO 2020
1. Introduction
Coronal mass ejections (CMEs) are bright outward-moving structures ejected from the Sun in all directions, detected in the white-light corona by means of coronagraphs. It is well established that magnetic fields associated with interplanetary CMEs, particularly those involving a southward component, can interact with Earth’s magnetosphere under certain conditions. Furthermore, from in situ detections of interplanetary CMEs it is known that there are particular cases, dubbed magnetic clouds (MCs), in which magnetic fields are organized following a helical structure (e.g., Bothmer & Schwenn 1998). As a consequence, from the practical point of view of space weather, there have been a number of recent efforts devoted to understanding how magnetic fields within CMEs are organized. This necessity has inspired several pieces of work. For instance, Yurchyshyn et al. (2007) looked for a correspondence between the orientation angle of elongated CMEs and of the axis angles of MCs. Marubashi et al. (2015) investigated the connection between the orientation and handedness of interplanetary magnetic flux ropes detected in situ and those of their solar source regions along magnetic polarity inversion lines. Following a similar reasoning, Savani et al. (2015) proposed a means of predicting the magnetic field structure of interplanetary CMEs arriving at Earth on the basis of their solar sources and the early characteristics of CMEs.
The study of CME morphology is of particular interest because it holds a close relationship with the magnetic field configuration. Both morphology and its evolution from the earliest stages of CMEs are key factors in order to associate their structure with that of their source region fields. In turn, these two factors can give hints on the involved initiation processes, and can benefit space weather forecasts on the basis of pre-eruptive field topologies in the low corona or developing CMEs in the corona.
A fundamental challenge in understanding the morphology of CMEs resides in the limitations imposed by coronagraphic observations. The coronagraph, the prime instrument that enables us to visualize CMEs, is capable of detecting the Thomson scattering brightness produced by coronal electrons, which are shaped by the prevailing magnetic fields. Despite the dramatic improvement in sensitivity and spatial resolution in the past few decades, a fundamental restriction persists: coronagraphs record the integrated brightness along the line of sight, and thus provide a 2D view of a 3D entity, as are CMEs.
Opposite to the early general belief that CMEs were spherical bubbles holding rotational symmetry (e.g., Crifo et al. 1983) supported by observations of Earth-directed circular halo CMEs (e.g., Howard et al. 1982), a number of studies have indicated that many CMEs (if not all, see Vourlidas et al. 2013) are organized along a major axis. Moreover, they present a topology consistent with that of a twisted magnetic flux rope in agreement with in situ detections of MCs. Evidence of helical magnetic fields can be found in several CMEs (e.g., Dere et al. 1999; Wood et al. 1999; Krall et al. 2001), while a configuration of arcades with cylindrical symmetry has also been proposed as consistent with CME observations (Moran & Davila 2004). Cremades & Bothmer (2004) noted that the projected views of a set of highly structured CMEs depended on the location and inclination of the neutral lines of their associated source regions, also suggestive of a helical flux rope configuration consistent with a cylindrical symmetry. According to this picture, the differing appearances of CMEs may be partly the result of different orientations of their main axis of symmetry with respect to the observer’s line of sight (LOS). This configuration leads to two extreme projected views: axial, detected when the LOS is aligned with the main axis of the cylinder, and lateral, seen when the LOS is perpendicular to this axis. These views have been also referred to respectively as edge-on and face-on views in the literature (Thernisien et al. 2006, 2009, 2011), and also axial and broadside (e.g., Krall & St. Cyr 2006). The axial (edge-on) view enables detection of the archetypal three-part structure, as in the left panel of Fig. 1, with the dark void at times being outlined by circular threads. The lateral (face-on; broadside) view, however, exhibits a more elongated but ragged and diffuse structure (see right panel of Fig. 1). Such cases were investigated by Cremades & Bothmer (2005) from the Earth-perspective images provided by the Large Angle Spectroscopic Coronagraph (LASCO; Brueckner et al. 1995) on board Solar and Heliospheric Observatory (SOHO; Domingo et al. 1995). They found that CMEs exhibiting their axial view in coronagraph images were in general much narrower than those exhibiting their lateral perspective.
![]() |
Fig. 1. Two CMEs exhibiting extreme cases of projection: axial view, with the LOS aligned with the CME main axis (left); and lateral view, with the LOS nearly perpendicular to the main axis (right). |
Until the mid-2000s, the only feasible view of the solar corona was that provided by instruments on the Sun-Earth line, thus hindering visualization of the same CME from different perspectives. This changed with the Solar-Terrestrial Relations Observatory Mission (STEREO; Kaiser et al. 2008), whose twin spacecraft enabled simultaneous imaging of the solar corona from multiple viewpoints. In fact, a proper and uncommon combination of factors, namely spacecraft separation, CME trajectory, and main axis orientation, enabled the first report of the simultaneous observation of the axial and lateral views of the same CME (Cabello et al. 2016). As pointed out by Mierla et al. (2009), the three-dimensional (3D) reconstruction of CME morphology from available coronagraphic images is an intrinsically ill-posed problem, given that a proper tomographic reconstruction would require a large number of images of the same CME from numerous different viewpoints. Although this is clearly a limitation to the 3D reconstruction of CMEs, images of the solar corona provided by the twin STEREO spacecraft, in combination with the terrestrial view provided by SOHO, constitute a triplet that better estimates the 3D morphology of CMEs, enabling the study of its evolution from their initiation in the low corona.
Inspired by the results of Cremades & Bothmer (2004, 2005) and Cabello et al. (2016) concerning the asymmetry in the axial versus lateral perspectives of CMEs, we examined the early evolution of the morphology of a distinct set of CMEs. The particular conditions of spacecraft location and CME trajectory under which the analyzed events occur reduce uncertainties typically involved in the 3D reconstruction via forward modeling. In Sect. 2.1 we explain why and how these events were analyzed. Section 2.2 describes the method used to reconstruct the overall 3D morphology of the CMEs under study. Section 3 presents the main findings, while Sect. 4 puts them in the context of previous results and discusses their implications.
2. Methodology
2.1. Data and event selection
White-light coronal images captured from multiple vantage points constitute the main data set used in this analysis. The Earth view is provided by LASCO C2 on board the legendary SOHO Mission, while the other perspectives are given by the COR1 and COR2 coronagraphs of the STEREO/Sun-Earth Connection Coronal and Heliospheric Investigation instrument suite (SECCHI; Howard et al. 2008). The early evolution of CMEs in the low corona was investigated by means of images in the extreme-ultraviolet (EUV) provided by the Atmospheric Imaging Assembly (AIA; Lemen et al. 2012) on board the Solar Dynamics Observatory (SDO; Pesnell et al. 2012) and by the SECCHI Extreme-Ultraviolet Imagers (EUVI) on board STEREO.
To minimize uncertainties due to projection effects, and with the fundamental goal of characterizing the morphology of a number of CMEs, we selected events according to the criteria described below. Otherwise, three views of the same event may not have been sufficient to correctly discern the orientation of its main axis of symmetry or the overall configuration of its magnetic field. First of all, the events under study were identified during times of spacecraft quadrature: November 2010–July 2011 (STEREO–SOHO quadrature) and December 2012–June 2013 (quadrature between the two STEREO spacecraft).
The next and most crucial criterion refers to the direction of propagation of the CMEs, which ideally should be perpendicular to the plane containing the observing spacecraft (STEREO, SOHO, and SDO) and thus to the ecliptic. Given that most CMEs originate at the two activity belts, those propagating at high latitudes are relatively rare. The reason for this requirement regarding the propagation direction has been addressed in detail by Cabello et al. (2016). In brief, three distinct views of the same CME are not enough to constrain the morphological parameters if all vantage points and the propagation direction lie on the same plane. This also agrees with the findings of Wood et al. (2009), who noted that the best results from forward modeling (see Sect. 2.2) are obtained when the CME’s appearance is very distinct from the different viewpoints, and analogous to the remark from Thernisien et al. (2009), who noted that imaging instruments placed away from the ecliptic plane would also help resolve ambiguities in the determination of CME parameters.
As a first approach to identifying events meeting these criteria, the SOHO LASCO CME Catalog1 (Yashiro et al. 2004) was inspected for CMEs with central position angle (CPA) within ±25° with respect to a CPA = 0° (north pole) and a CPA = 180° (south pole) during the two mentioned quadrature periods. We found a total of 35 events that met the criteria. However, as one of the fundamental objectives of this work is to study the evolution of the 3D morphology of CMEs from their beginnings in the low corona, only those whose eruption could be visualized by the STEREO/SECCHI EUVI and SDO/AIA instruments were considered, in total 12. One of the selected events at an early stage of its evolution is shown in Fig. 2.
![]() |
Fig. 2. Incipient CME on 14 December 2010 at 15:20 as seen in the low corona from STEREO EUVI-B 195 Å (left), SDO AIA 193 Å (center), and STEREO EUVI-A 195 Å (right). A previous image has been subtracted in all cases to enhance visibility of the event. The bottom row adds the GCS model fit superimposed as a green mesh, as described in Sect. 2.2. |
2.2. Determination of 3D parameters
As mentioned above, the general configuration of CMEs has proven to be in agreement with that of a twisted magnetic flux rope having cylindrical symmetry. The Graduated Cylindrical Shell (GCS) model proposed by Thernisien et al. (2009) is a modeling tool capable of adjusting an ad hoc 3D hollow structure, representative of such a flux rope, to two or three nearly simultaneous views of a CME recorded from different vantage points. This forward model is an upgrade to the version developed by Thernisien et al. (2006), based on the findings of Cremades & Bothmer (2004). The ad hoc geometrical surface that is fit to the images looks like a hollow croissant, consisting of two conical legs and a front similar to a section of a torus whose cross section increases with height. Despite its widespread use, the thoughtful application of the model is not trivial. As mentioned, three views of the same CME may not be enough to unambiguously find the true parameters describing the geometrical figure that best fits the CME, since various solutions may seem appropriate, especially if its propagation direction lies close to the plane containing the observing spacecraft.
To perform this investigation, nearly simultaneous triplets of images at various points in time describing the evolution of the selected events were fit by means of the GCS model. This means that fits were performed for each point in time when the CME was observed, from its beginning in the low corona (see, e.g., Patsourakos et al. 2010) as in the example shown in the bottom panel of Fig. 2, and during its propagation through the field of view (FOV) of the coronagraphs STEREO COR1 and COR2, and SOHO/LASCO C2. Initially, proxy values of latitude (θ), longitude (ϕ), and tilt (γ) are set to those of the source region and/or associated Hα filament. They are then adjusted by visual inspection on a best educated guess basis together with the three remaining parameters of the model, namely height (H), half width (α; half of the angle between the axes of the CME conical legs), and aspect ratio (κ; given by a(r)/r with a(r) being the variable radius of each of the conical legs). For each event, the triplets of images corresponding to the first and last analyzed times, respectively in EUVI 195 Å–AIA 193 Å and in COR2–LASCO C2, are adjusted first. When performing fits for the times comprised in between, the value of each parameter is allowed to vary between those corresponding to the first and last frames. These parameters are iteratively fit until the best visual match of the CME structure in each image and a smooth time evolution of the fit parameters are found. In the first iterations of the fitting process, some events showed a slight jump in parameter values when going from the COR1 to the COR2 FOV. To diminish this effect a compromise fit was performed. Naturally, this may result in the masking of slight changes; however, we are interested in the overall temporal evolution of CME widths. Parameters at the last fitted times usually show little to no variation; however, none of them is forced to be constant at a certain height. Figure 3 shows GCS model fits for all events superimposed to CME snapshots as viewed from the three perspectives. Images of the white-light corona are running-difference to enhance contrast, where special care was taken to avoid misinterpretations between flux rope and shock.
![]() |
Fig. 3. GCS model fits (green meshes) superimposed on snapshots of the 12 analyzed high-latitude CMEs. Each row corresponds to the three views of a particular event, while the left, middle, and right columns are the views from coronagraphs on board STEREO-B, SOHO, and STEREO-A, respectively. |
![]() |
Fig. 3. continued. GCS model fits (green meshes) superimposed on snapshots of the 12 analyzed high-latitude CMEs. |
For each event, plots like those in Fig. 4 were produced to show the evolution of the GCS parameters. CME height H is displayed as a function of time, while all other parameters (α, ϕ, θ, κ, γ) are shown as a function of height. Given that there is no merit function that quantitatively estimates the goodness of the fit, and thus the GCS parameters uncertainties, and that the interpretation of the geometry may differ from event to event and from observer to observer, precise quantitative errors for the fitted GCS parameters cannot be determined (see, e.g., Mierla et al. 2009), only roughly estimated (see Sect. 3). Moreover, we would like to emphasize that we are more interested in the temporal evolution of CME parameters than in their absolute values.
![]() |
Fig. 4. Evolution of the GCS parameters deduced for the analyzed event on 14 December 2010. From left to right and top to bottom: CME height H, aspect ratio κ, half angle α, longitude ϕ, latitude θ, and tilt γ. The evolution is shown as a function of height in all panels except for the height, which is shown against time. |
The main 3D attributes found for the investigated CMEs are summarized in Table 1. The first column displays the date and time of the first set of reconstructed parameters for each event. The second and third columns refer to the final central latitude and longitude arising from the GCS reconstructions. The fourth and fifth columns show the radial propagation speed at 6 R⊙ and the mean acceleration within the analyzed time interval. The remaining columns will be addressed in the following section.
Analyzed events and main attributes.
3. Results
Since we would like to understand how CMEs expand as they propagate outward in the corona, first we need to examine their height-time behavior. Figure 5 displays the height-time plots of the 12 analyzed events, where the symbols represent the data points arising from the GCS fitting. The dashed lines are second-order functions that best fit the time evolution of each event. To normalize time and avoid negative velocities, each fit curve starts at a height of 1.25 R⊙, which is a typical height for quiescent coronal cavities before being triggered to eruption (Gibson et al. 2006). In this sample, all events exhibit either an accelerated or a linear profile. To follow a smooth general trend, the quantities analyzed below rely on the height values that arise from the second-order fits rather than on the measurements themselves.
To understand how CMEs expand in the direction of their main axis and orthogonal to it, the angular widths in the axial (AWD) and lateral (AWL) directions are of particular interest. They correspond to the views shown in Fig. 1, and are determined as AWD = 2δ and AWL = 2α + δ; with δ being the arcsin(κ) (Thernisien et al. 2011). Both AWs are plotted as a function of height in Fig. 6, along with fit power laws extending back to the earliest data point of each event. A rapid increase in the AW in both axial and lateral directions is evident in the first 3 R⊙, followed by a stabilization phase at larger heights. The final values of AWD and AWL are listed in Cols. 6 and 7 of Table 1, followed by a measure of the corresponding uncertainties arising from the half angle α and the aspect ratio κ. These values have been estimated, as was done by Patsourakos et al. (2010), by verifying that the fits become unacceptable when varying α and κ beyond 20% of the chosen values for the frame at the lowest height. There is no correlation between the final values of AWL and AWD: events with large final AWL do not necessarily imply a large final AWD. It is thus worth noting that AWL and AWD depend on δ, but only AWL also depends on α. Since α and δ are free parameters of the model, AWL and AWD can vary independently. The initial change rate (deg h−1) reaches a few hundred for some events, and drops off to a few tens of degrees per hour by 3 R⊙.
![]() |
Fig. 5. Height-time plots arising from the GCS fits to the 12 analyzed events. The various symbol-color combinations correspond to the different events, and the dashed lines are constrained second-order fits. |
![]() |
Fig. 6. Height evolution of the CME AW in the axial (AWD, top panel) and lateral (AWL, bottom panel) directions. The symbols represent the data points, while the dashed lines are fit power laws extending back to the first data point of each event (see legend in the bottom panel, which also applies to the top plot). |
Figure 7 shows the difference between the derivatives with respect to height of the fit AWL and AWD curves in Fig. 6. This difference is always positive, denoting that the AW grows faster in the direction of the main axis of CMEs (lateral direction), than perpendicular to it (axial direction). Moreover, the difference is much larger in the first ∼3 R⊙, surpassing the 10° per solar radius. After that height both expansion rates and their differences drop to become negligible.
Expansion speeds in km s−1 have also been computed along both directions; their values at 6 R⊙ are listed in the last two columns of Table 1. The top panel of Fig. 8 displays the ratio of the speed along the main axes of CMEs (lateral speed) to the perpendicular (axial speed) as a function of height. It should be noted that both speeds are directly proportional at nearly all heights, with the lateral speed being on average ≈1.6 higher than the axial value (more precisely 1.56 ± 0.22). The bottom panel shows the ratio of lateral to radial (propagation) speeds for all analyzed events. In the first few solar radii, some events show a much more height-dependent ratio than others; after ∼4 R⊙ the ratio becomes nearly constant, with the average among all events being close to one (0.96 ± 0.63). In general, it can be said that the lateral speeds tend to be higher than the axial and similar to the radial (propagation) speeds, although this is not the case for a few events.
![]() |
Fig. 7. Difference between the fit height change rate of the angular widths in the lateral direction and in the axial direction for each event (dashed lines, see the legend). The black solid line is the average curve, starting at the first height where all events could be measured. The symbols here only indicate the location over the curves where the event was measured. The curve corresponding to 20110317 exceeds the scale and has been intentionally cut off. |
![]() |
Fig. 8. Speed ratios as a function of height. Top panel: lateral to axial speed. Bottom panel: lateral to radial (propagation) speed. The dashed colored lines correspond to individual events (see legend in bottom panel), which also applies to the top plot. The black solid line is their average, starting at the first height where all events could be measured. |
In order to find a means of evaluating the height at which the angular expansion of CMEs stabilizes, we define the settling height as the height at which AWL and AWD change less than 5% per solar radius. The settling heights of AWL and AWD are plotted against each other in Fig. 9, where each symbol represents a particular event. In general, events with larger settling heights in one direction tend to have larger settling heights in the other. From the figure, a slight trend toward larger settling heights for AWL can also be inferred.
4. Discussion and conclusions
Here we analyzed twelve polar CME events that took place during epochs of quadrature between spacecraft, and whose eruption could be tracked from their inception in the EUV low corona. The evolution of their morphology was characterized by means of multi-viewpoint coronal observations and a forward-modeling tool that was applied to several time instants. Even though 3D modeling of the morphology of CMEs is a problem with multiple possible solutions, a proper combination of spacecraft location and CME propagation direction allows a better constrained fit. More specifically, conditions for the good assessment of morphological parameters of CMEs are given in this study by (i) the stereoscopic views provided by the STEREO spacecraft located at the ecliptic while away from the Sun-Earth line, plus the view from Earth’s perspective provided by SOHO; (ii) a propagation direction of the analyzed CMEs nearly perpendicular to the plane containing these spacecraft; and (iii) the consideration of several time instants for the analysis, including the initial moments in the low corona.
One of the main findings of this study refers to the amount of asymmetry in the extents and expansion of CMEs in the lateral and axial directions of their embedded flux ropes. In fact, AWL is greater than AWD at all heights and for all events. This is imposed by the model, which defines α and δ as always positive, on the basis of the widths of observed face-on and edge-on CMEs by Cremades & Bothmer (2004). If any of the analyzed CMEs had a wide flux rope diameter, but a narrow face-on width, i.e., AWD > AWL, the model would not have been able to successfully fit it.
This can be also visualized in Fig. 10, which shows the ratio AWL/AWD as a function of height. This ratio presents steeper variations for individual events in the first few solar radii. It then stabilizes at differing values ranging from about 1.3 to 2.7, depending on the event. The average AWL/AWD ratio over height for all events is 1.84 ± 0.37, and is represented in the figure by the black solid line. Different angular widths depending on the flux rope orientation with respect to the observer had been suggested by St. Cyr et al. (2004) and Cremades & Bothmer (2005), who found a ratio AWL/AWD ≈ 1.6. We would like to note that the latter study measured widths of specific features – circular void and prominence extent – rather than full AWs. Moreover, these ratios arose from the comparison of the average lateral (broadside, AWL) and axial (AWD) widths found for different events from the perspective of SOHO, the single available view at the time. The first case study for which it was possible to determine the AWL/AWD ratio for the same CME was presented by Cabello et al. (2016), who found a value ≈1.6 as well. The analysis of Krall & St. Cyr (2006) performed on a parameterized model CME yielded a ratio of 1.9 for the best considered fit.
![]() |
Fig. 9. Settling height of AWL vs. that of AWD for all analyzed events. The solid line indicates a slope equal to 1. |
![]() |
Fig. 10. Ratio AWL/AWD. The dashed colored lines correspond to individual events. The black solid line is their average, starting at the first height where all events could be measured. |
Although significant work has been devoted to understanding how flux-rope cross sections evolve, little has been reported on the behavior of the lateral versus axial extents. In the case study by Wood et al. (2010) the CME flux rope was interpreted to be very fat at early times, and to become narrower and flatter as it evolved outward from the Sun. Janvier et al. (2014) discuss the unfeasibility of comparing distributions of CME apparent AWs and of the radius of magnetic clouds, given that CME apparent widths involve projection and selection effects that hinder the proper measurement of the extension of a flux rope’s axis and radius. Balmaceda et al. (2018), in turn, classified CMEs according to their morphological appearance, where the L (loop) and F (flux rope) types are attributed to CMEs seen face-on and edge-on, respectively, and found an average AW of L-CMEs nearly two times larger than that of F-CMEs.
Unlike the lateral versus axial asymmetry, the aspect ratio of CMEs, defined as the quotient of the CME center height over the CME’s half width (center of the flux rope/diameter), has been addressed by a fair number of studies in the literature. For instance, Krall et al. (2001) noticed that when looplike features could be measured in EIT, the loop expanded in width faster than it moved outward, implying a decreasing aspect ratio at low heights. Patsourakos et al. (2010) analyzed the 3D expansion of a CME and found two distinct phases, with the first of them also exhibiting smaller aspect ratios and the following phase showing an almost constant aspect ratio. A constant aspect ratio is commonly regarded as a sign of self-similarity, meaning that CMEs maintain their shape as they expand (Wood et al. 2009). Veronig et al. (2018) also find a fast decrease in the CME aspect ratio in the first stages of the CME genesis. A similar concept of aspect ratio has been investigated much farther in the heliosphere, in relation with the cross section of magnetic clouds regarded as cylindrical objects. It has been proposed that these cross sections are greatly flattened (e.g., Riley & Crooker 2004; Nieves-Chinchilla et al. 2018). However the “pancaking” may cease at some point before the arrival of these interplanetary structures to Earth (Savani et al. 2011).
We also find that the early times of the analyzed eruptions are the most important in terms of morphology, as demonstrated by the steep slope of the aspect ratio and of the difference between the AW change rates in the first few solar radii. This result is consistent with previous findings, like those of Zhang et al. (2004) and Patsourakos et al. (2010). The former noticed a superexpansion in the lateral direction during the acceleration phase in the lower corona, in contrast to the self-similar expansion pattern of constant AW observed in the outer corona. Patsourakos et al. (2010) also proposed two phases of evolution, with the initial one exhibiting a faster expansion and rise below 2 R⊙, followed by a nearly constant aspect ratio implying a self-similar evolution. Other observational studies revealing that the first instants of CME evolution present rapid changes and are decisive in their morphological development, report the acceleration of most events peaking at heights below 0.5 R⊙ (Bein et al. 2011) or 2 R⊙ (Joshi & Srivastava 2011).
The main results arising from the multi-viewpoint analysis of these 12 high-latitude CMEs can be summarized as follows:
-
Angular widths along the main axes of CMEs (AWL, face-on view) are larger than angular widths in the orthogonal direction (AWD, edge-on view) by an average factor of ≈1.8 (Figs. 6 and 10).
-
Both angular widths change considerably with height below ∼3 R⊙ (Fig. 6).
-
The growth rate of AWL is higher than that of AWD below ∼3 R⊙ (Fig. 7).
-
The ratio of the two expansion speeds, namely in the lateral (AWL) and axial (AWD) directions, is nearly constant after ∼4 R⊙, with an average ≈1.6 (Fig. 8, top).
-
On average, the expansion speed in the lateral (AWL) direction is similar to the outward radial propagation speed (Fig. 8, bottom).
-
The settling heights of AWL tend to be larger than those of AWD, showing some degree of proportionality.
According to the obtained results, we conclude that the expansion of the studied CMEs below ∼3 R⊙ cannot be considered self-similar and is evidently asymmetric, with higher expansion speeds in the lateral direction. This result has important implications for understanding and modeling not only CME morphology, but also shock generation, EUV wave formation, and their properties, such as profile and shape, among others.
As a follow-up to this work, we plan to understand our findings by means of MHD-based models. Although most of them do not consider the asymmetric expansion of CMEs and of their embedded flux ropes, they do address the mechanisms governing the overall expansion. For instance, Byrne et al. (2010) attribute the expansion of a modeled CME from 30° to 60° over a height range of 2–46 R⊙ to the convection with the ambient solar wind in a diverging geometry, and to the pressure gradient between the flux rope and solar wind. Mishra & Wang (2018) treat CMEs as axisymmetric cylinders assuming self-similar expansion, with the net force direction in agreement with the expansion acceleration, suggesting that the thermal pressure force is the internal driver of CME expansion and the Lorentz force the restraining agent. Liu et al. (2008) also found that that thermal pressure accounts for most of the acceleration of a CME, with the magnetic pressure contributing only to the acceleration early in the evolution up to ∼3 R⊙. In turn, Patsourakos et al. (2010) argue that the initial overexpansion of an erupting flux rope is caused neither by the decreasing ambient pressure as the flux rope rises nor by photospheric motions. Instead, they suggest that it results from an expansion of the flux surfaces of the poloidal flux external to the flux rope itself (driven by flux conservation) and/or from the rapid addition of flux by reconnection in a current sheet under a growing pre-existing flux rope in growth or a new one in formation. Expansion in the lateral direction is not addressed in this interpretation.
Out-of-the-ecliptic missions like Solar Orbiter will enable similar analyses that minimize uncertainties due to projection effects, but on events propagating close to the ecliptic plane, particularly Earth-directed. All currently existing coronagraphs lie on this plane, thus hindering proper 3D reconstruction of potentially geoeffective CMEs.
Acknowledgments
HC is member of the Carrera del Investigador Científico (CONICET). FAI is a postdoctoral fellow of CONICET and LAM holds a scholarship from the National Inter-University Council (Argentina). The authors are grateful to an anonymous referee for useful comments and suggestions. Authors appreciate financial support from UTN grants UTI4915TC and ME5445 and acknowledge use of data from the STEREO (NASA), SDO (NASA), and SOHO (ESA/NASA) missions, produced by the SECCHI, AIA, and LASCO international consortia. This work uses data from the SOHO/LASCO CME Catalog generated and maintained at the CDAW Data Center by NASA and the CUA in cooperation with NRL.
References
- Balmaceda, L. A., Vourlidas, A., Stenborg, G., & Dal Lago, A. 2018, ApJ, 863, 57 [NASA ADS] [CrossRef] [Google Scholar]
- Bein, B. M., Berkebile-Stoiser, S., Veronig, A. M., et al. 2011, ApJ, 738, 191 [Google Scholar]
- Bothmer, V., & Schwenn, R. 1998, Ann. Geophys., 16, 1 [Google Scholar]
- Brueckner, G. E., Howard, R. A., Koomen, M. J., et al. 1995, Sol. Phys., 162, 357 [NASA ADS] [CrossRef] [Google Scholar]
- Byrne, J. P., Maloney, S. A., McAteer, R. T. J., Refojo, J. M., & Gallagher, P. T. 2010, Nat. Commun., 1, 74 [Google Scholar]
- Cabello, I., Cremades, H., Balmaceda, L., & Dohmen, I. 2016, Sol. Phys., 291, 1799 [NASA ADS] [CrossRef] [Google Scholar]
- Cremades, H., & Bothmer, V. 2004, A&A, 422, 307 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cremades, H., & Bothmer, V. 2005, in Coronal and Stellar Mass Ejections, eds. K. Dere, J. Wang, & Y. Yan, IAU Symp., 226, 48 [NASA ADS] [Google Scholar]
- Crifo, F., Picat, J. P., & Cailloux, M. 1983, Sol. Phys., 83, 143 [NASA ADS] [CrossRef] [Google Scholar]
- Dere, K. P., Brueckner, G. E., Howard, R. A., Michels, D. J., & Delaboudiniere, J. P. 1999, ApJ, 516, 465 [NASA ADS] [CrossRef] [Google Scholar]
- Domingo, V., Fleck, B., & Poland, A. I. 1995, Sol. Phys., 162, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Gibson, S. E., Foster, D., Burkepile, J., de Toma, G., & Stanger, A. 2006, ApJ, 641, 590 [NASA ADS] [CrossRef] [Google Scholar]
- Howard, R. A., Michels, D. J., Sheeley, N. R., & Koomen, M. J. 1982, ApJ, 263, L101 [NASA ADS] [CrossRef] [Google Scholar]
- Howard, R. A., Moses, J. D., Vourlidas, A., et al. 2008, Space Sci. Rev., 136, 67 [NASA ADS] [CrossRef] [Google Scholar]
- Janvier, M., Démoulin, P., & Dasso, S. 2014, Sol. Phys., 289, 2633 [NASA ADS] [CrossRef] [Google Scholar]
- Joshi, A. D., & Srivastava, N. 2011, ApJ, 739, 8 [NASA ADS] [CrossRef] [Google Scholar]
- Kaiser, M. L., Kucera, T. A., Davila, J. M., et al. 2008, Space Sci. Rev., 136, 5 [NASA ADS] [CrossRef] [Google Scholar]
- Krall, J., & St Cyr, O. C. 2006, ApJ, 652, 1740 [NASA ADS] [CrossRef] [Google Scholar]
- Krall, J., Chen, J., Duffin, R. T., Howard, R. A., & Thompson, B. J. 2001, ApJ, 562, 1045 [NASA ADS] [CrossRef] [Google Scholar]
- Lemen, J. R., Title, A. M., Akin, D. J., et al. 2012, Sol. Phys., 275, 17 [NASA ADS] [CrossRef] [Google Scholar]
- Liu, Y. C. M., Opher, M., Cohen, O., Liewer, P. C., & Gombosi, T. I. 2008, ApJ, 680, 757 [NASA ADS] [CrossRef] [Google Scholar]
- Marubashi, K., Akiyama, S., Yashiro, S., et al. 2015, Sol. Phys., 290, 1371 [Google Scholar]
- Mierla, M., Inhester, B., Marqué, C., et al. 2009, Sol. Phys., 259, 123 [NASA ADS] [CrossRef] [Google Scholar]
- Mishra, W., & Wang, Y. 2018, ApJ, 865, 50 [NASA ADS] [CrossRef] [Google Scholar]
- Moran, T. G., & Davila, J. M. 2004, Science, 305, 66 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Nieves-Chinchilla, T., Linton, M. G., Hidalgo, M. A., & Vourlidas, A. 2018, ApJ, 861, 139 [NASA ADS] [CrossRef] [Google Scholar]
- Patsourakos, S., Vourlidas, A., & Kliem, B. 2010, A&A, 522, A100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pesnell, W. D., Thompson, B. J., & Chamberlin, P. C. 2012, Sol. Phys., 275, 3 [NASA ADS] [CrossRef] [Google Scholar]
- Riley, P., & Crooker, N. U. 2004, ApJ, 600, 1035 [Google Scholar]
- Savani, N. P., Owens, M. J., Rouillard, A. P., et al. 2011, ApJ, 731, 109 [NASA ADS] [CrossRef] [Google Scholar]
- Savani, N. P., Vourlidas, A., Szabo, A., et al. 2015, Space Weather, 13, 374 [NASA ADS] [CrossRef] [Google Scholar]
- St Cyr, O. C., Cremades, H., Bothmer, V., Krall, J., & Burkepile, J. T. 2004, AGU Fall Meet. Abstr., 2004, SH22A-04 [Google Scholar]
- Thernisien, A. F. R., Howard, R. A., & Vourlidas, A. 2006, ApJ, 652, 763 [Google Scholar]
- Thernisien, A., Vourlidas, A., & Howard, R. A. 2009, Sol. Phys., 256, 111 [Google Scholar]
- Thernisien, A., Vourlidas, A., & Howard, R. A. 2011, J. Atm. Solar-Terr. Phys., 73, 1156 [Google Scholar]
- Veronig, A. M., Podladchikova, T., Dissauer, K., et al. 2018, ApJ, 868, 107 [Google Scholar]
- Vourlidas, A., Lynch, B. J., Howard, R. A., & Li, Y. 2013, Sol. Phys., 284, 179 [NASA ADS] [Google Scholar]
- Wood, B. E., Karovska, M., Chen, J., et al. 1999, ApJ, 512, 484 [Google Scholar]
- Wood, B. E., Howard, R. A., Thernisien, A., Plunkett, S. P., & Socker, D. G. 2009, Sol. Phys., 259, 163 [NASA ADS] [CrossRef] [Google Scholar]
- Wood, B. E., Howard, R. A., & Socker, D. G. 2010, ApJ, 715, 1524 [NASA ADS] [CrossRef] [Google Scholar]
- Yashiro, S., Gopalswamy, N., Michalek, G., et al. 2004, J. Geophys. Res. (Space Phys.), 109, 7105 [CrossRef] [Google Scholar]
- Yurchyshyn, V., Hu, Q., Lepping, R. P., Lynch, B. J., & Krall, J. 2007, Adv. Space Res., 40, 1821 [NASA ADS] [CrossRef] [Google Scholar]
- Zhang, J., Dere, K. P., Howard, R. A., & Vourlidas, A. 2004, ApJ, 604, 420 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
All Figures
![]() |
Fig. 1. Two CMEs exhibiting extreme cases of projection: axial view, with the LOS aligned with the CME main axis (left); and lateral view, with the LOS nearly perpendicular to the main axis (right). |
In the text |
![]() |
Fig. 2. Incipient CME on 14 December 2010 at 15:20 as seen in the low corona from STEREO EUVI-B 195 Å (left), SDO AIA 193 Å (center), and STEREO EUVI-A 195 Å (right). A previous image has been subtracted in all cases to enhance visibility of the event. The bottom row adds the GCS model fit superimposed as a green mesh, as described in Sect. 2.2. |
In the text |
![]() |
Fig. 3. GCS model fits (green meshes) superimposed on snapshots of the 12 analyzed high-latitude CMEs. Each row corresponds to the three views of a particular event, while the left, middle, and right columns are the views from coronagraphs on board STEREO-B, SOHO, and STEREO-A, respectively. |
In the text |
![]() |
Fig. 3. continued. GCS model fits (green meshes) superimposed on snapshots of the 12 analyzed high-latitude CMEs. |
In the text |
![]() |
Fig. 4. Evolution of the GCS parameters deduced for the analyzed event on 14 December 2010. From left to right and top to bottom: CME height H, aspect ratio κ, half angle α, longitude ϕ, latitude θ, and tilt γ. The evolution is shown as a function of height in all panels except for the height, which is shown against time. |
In the text |
![]() |
Fig. 5. Height-time plots arising from the GCS fits to the 12 analyzed events. The various symbol-color combinations correspond to the different events, and the dashed lines are constrained second-order fits. |
In the text |
![]() |
Fig. 6. Height evolution of the CME AW in the axial (AWD, top panel) and lateral (AWL, bottom panel) directions. The symbols represent the data points, while the dashed lines are fit power laws extending back to the first data point of each event (see legend in the bottom panel, which also applies to the top plot). |
In the text |
![]() |
Fig. 7. Difference between the fit height change rate of the angular widths in the lateral direction and in the axial direction for each event (dashed lines, see the legend). The black solid line is the average curve, starting at the first height where all events could be measured. The symbols here only indicate the location over the curves where the event was measured. The curve corresponding to 20110317 exceeds the scale and has been intentionally cut off. |
In the text |
![]() |
Fig. 8. Speed ratios as a function of height. Top panel: lateral to axial speed. Bottom panel: lateral to radial (propagation) speed. The dashed colored lines correspond to individual events (see legend in bottom panel), which also applies to the top plot. The black solid line is their average, starting at the first height where all events could be measured. |
In the text |
![]() |
Fig. 9. Settling height of AWL vs. that of AWD for all analyzed events. The solid line indicates a slope equal to 1. |
In the text |
![]() |
Fig. 10. Ratio AWL/AWD. The dashed colored lines correspond to individual events. The black solid line is their average, starting at the first height where all events could be measured. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.