Free Access
Issue
A&A
Volume 644, December 2020
Article Number A137
Number of page(s) 14
Section The Sun and the Heliosphere
DOI https://doi.org/10.1051/0004-6361/202038781
Published online 11 December 2020

© ESO 2020

1. Introduction

Coronal mass ejections (CMEs) are expulsions of magnetised plasma from the solar atmosphere and represent the most energetic events in the Solar System. These eruptions are ultimately the result of a magnetic energy storage and release process that is the consequence of the emergence of magnetic flux from the solar interior combined with photospheric motions that further evolve the coronal field. The free magnetic energy stored in the corona is in the form of field-aligned electric currents, and it is this energy that powers CMEs.

Separate processes are required to describe the slow build-up of electric currents in the corona before a CME and the rapid onset and acceleration of the eruption. Following the convention of Aulanier et al. (2010), in this work, CME initiation mechanisms are classified as either triggers or drivers. Triggers create a stable structure before an eruption and the driving process is responsible for the fast expansion phase of the eruption. A wide range of mechanisms have been identified as CME triggers, including sunspot rotation (e.g. Yan et al. 2012; Török et al. 2013; James et al. 2017), magnetic reconnection (e.g. flux cancellation; van Ballegooijen & Martens 1989; Yardley et al. 2018, and tether-cutting; Moore et al. 2001), and the helical kink instability of a flux rope (Hood & Priest 1979; Török & Kliem 2005). There are two main groups of theories pertaining to how the rapid expansion of a CME is driven. One group assumes that CMEs are driven by reconnection, which may or may not involve a flare, for example, in the “breakout” scenario (Antiochos et al. 1999; Temmer et al. 2010; Karpen et al. 2012). The other group assumes that CMEs are driven by an ideal magnetohydrodynamic (MHD) instability or loss of equilibrium involving a flux rope, such as the torus instability (van Tend & Kuperus 1978; Kliem & Török 2006; Démoulin & Aulanier 2010). For a summary of triggers and drivers, see Table 1 of Green et al. (2018).

A substantial number of observational studies have been carried out to test the theories that involve flux ropes present at the time of eruption. Early observational studies tested the flux cancellation mechanism of van Ballegooijen & Martens (1989) whereby a low altitude flux rope forms via magnetic reconnection in the photosphere or chromosphere. This low altitude flux rope may have a so-called bald patch separatrix (BPS) configuration in which the underside of the rope (that runs along the photospheric polarity inversion line; PIL) is rooted in the dense, high plasma β environment of the lower atmosphere, as shown by Green & Kliem (2009). Such ropes can be formed in both emerging active regions in which opposite polarity fragments collide, or in decaying active regions where the separation, natural fragmentation, and dispersal of emerging bipoles can bring polarities into contact with neighbouring polarities of opposite sign in sufficiently complex active regions.

More recently, extreme-ultraviolet (EUV), X-ray, and radio observations have shown evidence of pre-eruption flux ropes containing 10 MK plasma (known as hot flux ropes) and that have their underside in the low plasma β environment of the corona. These flux ropes are seen best when they are close to the solar limb (e.g. Reeves & Golub 2011; Patsourakos et al. 2013; Nindos et al. 2015). Hot flux ropes also appear to be formed as a consequence of magnetic reconnection, but this time in the corona, where the reconnection not only leads to the formation of the flux rope, but also a confined flare (Patsourakos et al. 2013; James et al. 2017).

The formation of a flux rope in the corona necessitates quasi-separatrix layers (QSLs) separating the rope from the external field it is embedded within. James et al. (2018) showed that a hot flux rope had a specific magnetic configuration known as a hyperbolic flux tube (HFT; Titov et al. 2002). In an HFT, two intersecting QSLs are present and an X-line runs along the underside of the rope. In this configuration, perturbations to the coronal field caused by, for example, photospheric flows may preferentially lead to reconnection in the current sheet beneath the flux rope (Demoulin et al. 1996). This reconnection could build additional magnetic flux in to the flux rope over time, and perhaps even play a role in the eruption of the rope.

In order to fully understand the physics of CMEs and be able to predict them, relevant triggers must be identified so that pre-eruptive structures can be recognised, and knowledge of whether an eruption of these structures can be successfully driven must also be obtained. However, most observational studies that identified higher altitude coronal flux ropes at the limb, which are likely to be of the HFT configuration, were unable to observe the evolution of the photosphere beneath the forming flux ropes, hindering the identification of the trigger mechanisms involved in these CMEs.

This work follows on from two studies that investigated an HFT rope that formed and erupted near disc centre allowing both the photospheric and coronal evolution to be tracked (James et al. 2017; hereafter Paper I, and James et al. 2018; hereafter Paper II). The combination of observations and modelling in Papers I and II showed the presence of an HFT flux rope with its axis ∼0.2 R above the photosphere. In this event, the flux rope formed as a consequence of magnetic reconnection in the corona, driven as photospheric orbiting motions of flux fragments about each other led to a collision of flux systems. While magnetic reconnection in the photosphere associated with the formation of BPS ropes injects photospheric plasma into the structure (Baker et al. 2013), spectroscopic measurements in Paper I showed that the studied flux rope had a coronal plasma composition, supporting the conclusion that the orbiting motions had indeed driven magnetic reconnection in the corona transforming a sheared arcade into an HFT rope.

As well as likely playing a key role in flux rope formation, orbiting motions of magnetic flux systems about each other can inject magnetic energy and helicity in to the arcade field that overlies the rope causing the field to inflate and weakening its confining effect on the underlying flux rope (Török et al. 2013). This scenario is consistent with the findings of Paper II that suggested the flux rope eruption was driven by the torus instability due to the magnetic field gradient above the rope reaching a critical value. The motion of magnetic flux fragments around each other may also be framed in the broader observational context of sunspot rotation (spin). The causal connection between rotating sunspots and the onset of solar activity has been well-studied and linked to the formation of sigmoids, flaring, and CMEs (e.g. Gerrard et al. 2003; Brown et al. 2003; Yan et al. 2012; Török et al. 2013; Vemareddy et al. 2016). The sunspots in each of these studies feature multiple umbrae within one rotating sunspot, and may therefore be similar to the “orbiting” scenario of eruption triggering described above.

Papers I and II represent a case-study of a single event. A natural extension of these previous works is to investigate whether this process is observed in other eruptive active regions and could therefore be viewed as a CME trigger mechanism. In this study, we investigate the formation of five HFT flux ropes in four active regions, including the active region studied in Paper I and Paper II. Photospheric orbiting motions are quantified in each active region to identify whether the same processes that formed the flux rope in Papers I and II are systematic in HFT flux rope formation. The criteria used to select events are given in Sect. 2. Sections 3.13.3 describe the methods used to quantify photospheric orbiting motions, estimate the height of coronal structures, and determine the chiralities of each active region, respectively. The observations of each eruption and the measured motions are detailed in Sect. 4, the results are interpreted and discussed in Sect. 5, and conclusions are presented in Sect. 6.

2. Data and event selection

Events with HFT flux ropes present in the corona before erupting as CMEs were selected from the era of the Solar Dynamics Observatory (SDO; Pesnell et al. 2012). Observational signatures of HFT flux ropes include, for example, sigmoids and plasmoids with underlying arcade fields (Reeves & Golub 2011). Arcades can be identified in many of the EUV passbands of the SDO’s Atmospheric Imaging Assembly (AIA; Lemen et al. 2012), but sigmoids and plasmoids appear only in those bands with the hottest temperature responses ∼10 MK. The 131 Å channel of AIA was used to identify sigmoids and plasmoids in the corona because it has a peak in temperature response at ∼107.0 K, however, it also has another peak at the lower temperature of 105.6 K (see Fig. 11 of Boerner et al. 2012). This means AIA 131 Å images contain contributions from both relatively hot and cold plasma. To ensure that only hot plasma signatures were identified, images from the relatively cool 171 Å channel (temperature response peak at 105.8 K) were used for comparison. Any feature that appeared in the 131 Å images but did not appear in the 171 Å channel was assumed to correspond to hotter, ∼107.0 K plasma, and therefore may be a flux rope signature.

Additionally, observations from the 193 Å, and 211 Å channels of AIA (that image coronal plasma at ∼106.2 K and ∼106.3 K respectively) were used to identify EUV dimmings, and the 1600 Å channel (that images plasma in the lower solar atmosphere at ∼105.0 K) was used to locate flare ribbons. EUV dimmings and flare ribbons manifest at the footpoints of erupting flux ropes (Janvier et al. 2014), and are therefore critical for determining the locations where the photospheric motions associated with flux rope formation within the orbiting scenario should be studied. AIA data were processed to level 1.5 using the aia_prep routine available in SolarSoft.

The CMEs, and the solar hemisphere from which they originated, were identified using images from the Large Angle and Spectrometric Coronagraph (LASCO; Brueckner et al. 1995) onboard the Solar and Heliospheric Observatory (SOHO; Domingo et al. 1995) and the Solar Terrestrial Relations Observatory (STEREO; Kaiser et al. 2008) coronagraphs. The STEREO spacecraft were ∼ ± 120° from the Sun-Earth line at the times of the selected events (at first in 2012, and then having approximately swapped positions with each other by 2017, albeit with only STEREO A still functioning), giving complementary perspectives of the CMEs.

The evolution of the photospheric flux and horizontal motions in each CME source region were studied using data from the Helioseismic and Magnetic Imager (HMI; Scherrer et al. 2012) onboard SDO. White-light continuum images and radial magnetic field observations were obtained from the SHARP data series (Spaceweather HMI Active Region Patch; Bobra et al. 2014). In events where the SHARP magnetograms contained more than one region of strong magnetic flux, such as multiple active regions, the data were cropped to contain only the desired source active region. Pixels with magnetic flux densities lower than 30 G were excluded from the flux calculations to reduce the effect of noise. In order to study the evolution of the photosphere before the eruptions, observations of the photospheric magnetic field over a number of days were analysed. The radial magnetic field measurement is increasingly affected by noise outside ±60° from disc-centre (Liu et al. 2012), so only CME source regions that evolved within this longitude range were used. The heliographic latitudes and longitudes of the studied active regions were referenced from the Debrecen Photoheliographic Data sunspot catalogue1. Furthermore, because flux cancellation is not expected during the formation of HFT flux ropes, the radial field magnetograms were used to confirm that the active regions selected in this study did not exhibit significant flux cancellation along the CME-source polarity inversion lines.

Full-Sun integrated X-ray light curves were obtained from the Geostationary Operational Environmental Satellite X-ray Sensor (GOES XRS) system and used in combination with EUV images from SDO/AIA to identify flares and brightenings in the CME source active regions.

In summary, the selection criteria for events with hot flux ropes were: a white light CME was clearly identifiable in coronagraph data; the CME source region was on the Earth-facing side of the Sun within ±60° of disc-centre; a sigmoid or a plasmoid was observed before eruption in an imaging waveband that senses plasma at a temperature of ∼107 K; an arcade of loops was present under the identified flux rope before eruption; and no significant magnetic flux cancellation was observed along the polarity inversion line associated with the CME.

Using these criteria led to the selection of five CMEs from four active regions. The dates of the five CMEs are given in Table 1, along with information about the source active regions. In each of these four active regions, magnetic flux fragments were observed to orbit each other during the days before the CMEs.

Table 1.

CME source regions.

3. Methods

3.1. Quantifying the orbital motion of emerging magnetic flux fragments

This study builds on that presented in Paper I, in which the orbital motion of a magnetic flux fragment around a sunspot was estimated by-eye using two white-light continuum images that were taken 24 hours apart. However, in this work, a more systematic method is used to track the motion of individual magnetic flux fragments, as explained in this section.

HMI continuum images from the SHARP data series were used to study the motion of magnetic flux fragments in the photosphere. The full passage of each active region across the solar disc was followed at an image cadence of 6 hours. Contours were set on each HMI SHARP continuum image to encircle individual umbral fragments, and the flux-weighted centre of each fragment was found using the radial magnetic field strength inside each contour (see Fig. 1). Across different active regions and at different times during the evolution of a given active region, different contour values were chosen between 10 000–25 000 DN/s in the continuum images to best isolate specific fragments.

thumbnail Fig. 1.

Identifying magnetic flux fragments. Left: contours are drawn at fixed intensity levels on HMI continuum images to enclose individual umbral fragments. Right: the continuum contours are overlaid on the radial component of magnetic field. Positive polarity field is shown in white, and negative in black. The flux-weighted centroid of each fragment is shown in red and is assigned a unique identifying number.

In the first image of a sequence (at time t1), two fragments were selected: one as the “central” fragment and one as the “orbiting” fragment, and a vector (v1) was drawn between their flux-weighted centroids. In the next image in the sequence (at time t2), the orbiting fragment had moved relative to the central fragment, and a new vector (v2) was drawn between the newly-positioned flux-weighted centroids. By comparing the vectors in successive images (illustrated in Fig. 2), the angle the orbiting fragment had moved around the central fragment by was calculated using the relation

(1)

thumbnail Fig. 2.

Illustration of the method used to track the relative orbital motion of magnetic flux fragments.

This gives 0 ≤ θ21 ≤ 180, and therefore does not distinguish between rotation in the clockwise and anticlockwise directions. In order to obtain a signed orbit angle, −180 ≤ θ21 ≤ 180, the relations

(2)

were used, where vn is the vector normal to v1 and v2, oriented positively out of the image plane towards the observer. Here, the convention is chosen that clockwise (anticlockwise) motions correspond to positive (negative) orbit angles. This process was repeated for every image in the time sequence to track the orbital motion over time.

In most events, separate pairs of fragments were tracked at different times to best quantify the observed motions in each evolving active region. The vectors between different pairs of tracked fragments are colour-coded in, for example, Fig. 4, and the corresponding orbit rates of different fragments are colour-coded in, for example, Fig. 5.

At some time steps, it was not possible for the tracking routine to distinguish two fragments, for example after they had merged together. However, the seemingly-merged sunspots still often seemed to be comprised of distinct fragments, as evidenced by observed light bridges or the later re-separation in to similar fragments. In these cases, the solid-body rotation of the merged sunspot was quantified as a proxy for how the previously-individual fragments continued to move around each other. This was done by linearly fitting a vector to the major axis of the elliptical sunspot umbra to best describe the orientation of the sunspot and following the rotation of this vector through time in the same way as described by Eqs. (1) and (2). The effectiveness of this method is limited in the case of merged sunspots that are approximately circular.

It is difficult to quantify error in the orbit rates. The primary source of uncertainty arises from the computed flux-weighted centroids that can move over time as multiple fragments coalesce. This can lead to instances where the orbit rate appears to stop or even reverse direction briefly. It is hard to quantify this effect, but we attempt to account for it by showing three-point moving averages of the measured orbit rates, for example, in Fig. 5, in order to better show the qualitative trend of the orbiting with less of an effect from frame-to-frame error.

3.2. Estimating flux rope altitude

Determining the heights of the identified flux ropes above the photosphere is important for supporting the interpretation of an HFT configuration and for examining the overall stability of a flux rope. In this work, flux rope heights are estimated from observations using geometrical considerations and assumed symmetry.

The assumption is made that the mid-point of either an observed sigmoid (which is interpreted as representing plasma emission from the underside of a flux rope) or the bottom of a plasmoid lies radially above the centre of a photospheric PIL. That is, that the flux rope is not inclined relative to the perpendicular direction. In reality, this will not necessarily be true, as the symmetry of the system is affected by many factors (e.g. the spatial flux distribution at the photosphere, the distribution of currents in the active region, the dynamical evolution of the sigmoid, etc.). The asymmetries and projection effects make it difficult to implement a method to automatically determine the mid-points of the sigmoids, so we chose this point by-eye based on the geometry of the sigmoid. Specifically, we found the middle of the relatively straight central section between the curved ends of the fully-illuminated S-shaped loops. Due to the nature of the assumptions being made, any resulting heights should be taken as order-of-magnitude estimations only.

From plane-of-sky solar observations, two coordinates were taken: the coordinate of the centre of the coronal sigmoid or plasmoid-base (xcor, ycor) and the point at the middle of the photospheric PIL between sunspots (xpil, ypil). Then, the height of the coronal point above the photospheric PIL was estimated using either the difference in x coordinates or the difference in y coordinates as

(3)

If the assumption that the coronal point lies radially above the photospheric point is correct, then hx = hy, but otherwise these two equations will give different height estimates. Whenever this method was applied in this work, both hx and hy were computed and the average was taken to give a qualitative estimate of height. Errors in hx and hy were calculated from the uncertainty in selecting the coordinates of the PIL and sigmoid centres.

3.3. Chirality determination

As part of investigating the formation of the hot flux ropes, the chirality of each active region was determined in two ways. Firstly, the chirality of magnetic flux emerging in to the photosphere was inferred from observations of magnetic tongues (López Fuentes et al. 2000; Luoni et al. 2011) in radial field HMI observations. If the leading polarity in a bipolar emerging active region extends to the south (north) of the trailing polarity, the emerging magnetic flux has right-handed (left-handed) twist.

Secondly, the chirality of the coronal magnetic field was indicated by observations of forward-S or reverse-S sigmoids in EUV observations from the 131 Å channel of AIA. Sigmoids with forward-S (reverse-S) shapes are a signature of magnetic field with right-handed (left-handed) twist in the corona (Longcope et al. 1998).

4. Observations and results

4.1. 13 March 2012

4.1.1. Coronal evolution

On 13 March 2012, a CME and accompanying M7.9 flare occurred in NOAA AR 11429 at ∼17:10 UT. The CME was first observed in white-light by LASCO C2 at 17:36 UT (Fig. 3a) when the active region was 59° west of central meridian. An inverse-S (left-handed) sigmoid and underlying arcade brightened in the active region at 16:30 UT, shortly before the onset of the eruption (see Fig. 3b). The ends of the sigmoid were rooted in the east and west parts of the active region, and twin EUV dimmings and hooked flare ribbons were observed in these locations during the eruption itself (Figs. 3c and 3d respectively). The centre of the pre-eruptive sigmoid (at 16:47 UT) is estimated to be at an altitude of 56 ± 5 Mm above the photosphere, however, this is found by averaging two rather different estimates, with hx = 91 ± 8 Mm and hy = 22 ± 3 Mm. The relatively large difference between hx and hy in this estimate may be due to the strong projection effects arising from the active region’s position nearly 60° west of disc-centre. Due to the lack of overlap of the estimated heights, the height estimate here should be treated with caution. This is the only event in this study where there is no overlap between the hx and hy estimates. Together, the observations detailed above suggest the presence of an HFT flux rope with footpoints in the east and west sides of the active region that erupted during the CME.

thumbnail Fig. 3.

Observations of the 13 March 2012 eruption. (a) A white-light CME (indicated by the white arrow) seen by LASCO C2. (b) AIA 131 Å image of a sigmoid (highlighted by a dashed red line) brightening above a flaring arcade (indicated by a red arrow). (c) AIA 211 Å base difference image showing twin EUV dimmings (marked by two yellow arrows). (d) AIA 1600 Å image showing hooked flare ribbons (indicated by two black arrows).

The previous CME from NOAA AR 11429 occurred three days prior, on 10 March at 17:44 UT, in association with an M8.4 flare. Similarly to the CME on 13 March, the flare arcade produced during this event on 10 March spanned the full width of the active region, indicating that both eruptions originated from the active region’s internal polarity inversion line. Furthermore, a white-light CME with a bright front and cavity was seen at 18:18 UT in LASCO C2 images, suggesting that the CME had a flux rope structure (Vourlidas et al. 2013 and the references within). Assuming that the 10 March CME fully ejected any flux rope that was present in the active region, the flux rope that erupted on 13 March must have formed at some time during the 3 days between eruptions. The specific time(s) of flux rope formation can be inferred by identifying solar flares (and therefore magnetic reconnection episodes) during this time.

GOES detected four flares in NOAA AR 11429 between the CMEs on 10 and 13 March 2012. Of these four, two flares illuminated arcades that spanned across the whole active region in the same way as the flux rope eruptions and exhibited flare ribbons in the 1600 Å channel of AIA that were in the same locations as were seen during the 13 March CME. Furthermore, hot plasma structures in the AIA 131 Å channel appeared above the flare arcades with shapes similar to – and in the same location as – the sigmoid associated with the 13 March CME. These two flares were a C4.1 flare on 12 March at 22:20 UT and a C3.1 flare on 13 March at 06:55 UT. Observations from the SOHO/LASCO and STEREO coronagraphs confirm these were confined events, supporting the scenario that these flares were associated with the construction, but not eruption, of a flux rope. These observations suggest the flux rope had begun to form at least 19 hours before the onset of the CME on 13 March 2012.

4.1.2. Photospheric evolution

NOAA AR 11429 rotated on to the solar disc on 3 March 2012 and contained two sunspots at that time. The leading sunspot had positive polarity and the trailing spot was negative, meaning the active region did not follow Hale’s law for a northern-hemisphere region in solar cycle 24. A new flux emergence episode began on 4 March 2012 at two sites between the pre-existing sunspots. The orientation of the magnetic tongues of this newly emerging flux indicated that it had left-handed twist (negative chirality).

The positive and negative fragments of the emerging flux moved towards the pre-existing sunspots of the same polarities. As the fragments approached the sunspots, they orbited around them in an anticlockwise sense. The strongest orbiting was observed around the leading sunspot, and the EUV dimmings (shown in Fig. 3c) suggest that the flux rope had a footpoint in this leading spot. Therefore, the leading sunspot was chosen as the location to quantify orbiting motion.

The motion of positive magnetic flux around the leading sunspot was tracked by following distinct fragments using the method described in Sect. 3.1. Here, we refer to certain fragments of magnetic flux with letters A–C, as in Fig. 4. At first, emerging fragment B moved towards and anticlockwise around pre-exisiting fragment A, and then fragment C emerged and orbited anticlockwise around fragment B. The orbiting motion is quantified in Table 2 and also represented in Figs. 4 and 5. In addition, Fig. 5 also contains the GOES X-ray lightcurve over the duration of the orbiting and the evolution of the positive and negative magnetic flux of the active region for comparison.

thumbnail Fig. 4.

Anticlockwise motion of newly-emerged flux around the pre-existing positive (leading) sunspot in NOAA AR 11429. Tracked fragments of flux are labelled A–C. In each image, vectors are drawn either to connect the flux-weighted centroids of two orbiting fragments (panels a–d) or to best-fit the major axis of merged fragments (panels e and f). In the bottom two panels, red circles highlight the position of a fragment (C) whose ≈180° orbit from the south of the sunspot to the north between panels was not successfully tracked using the vector method. This figure is available as an online movie.

thumbnail Fig. 5.

Top: measured orbital motion of chosen magnetic flux fragments in 6-hour intervals. The colours correspond to different choices of fragments (see Fig. 4), except for the purple block which corresponds to the average orbiting estimated by-eye that was not detected by the automated method. The black curve represents the orbit angles smoothed using a three-point moving average. Middle: full-disc integrated GOES soft X-ray lightcurve. Bottom: the evolution of magnetic flux in NOAA AR 11429, made using the radial magnetic field component, Br, of the HMI SHARP data series and smoothed with a 24-hour moving average. Vertical dashed lines indicate the times of the CME onsets on 10 and 13 March 2012, and the thin, light-blue vertical bars show the timings of confined flares associated with flux rope formation described in Sect. 4.1.1.

Table 2.

Orbital motions of umbral flux fragments relative to previously-emerged umbrae and the observed active region chiralities.

Between 5 March and 7 March, the tracked fragments orbited anticlockwise by 74° (an average of 37° per day), and from 7 March until midday on 8 March, another pair orbited anticlockwise by 36° (an average of ≈29° per day). During these times, three X-flares occurred that were associated with halo CMEs: one on 5 March and two in quick succession on 7 March. The flux ropes that erupted on 7 March were not rooted in the leading, coalescing sunspot and formed low-down in the atmosphere as a result of shearing motions elsewhere in the active region (Chintzoglou et al. 2015). For more on the 7 March eruptions, see Wang et al. (2014), Syntelis et al. (2016), and Baker et al. (2019). The total magnetic flux of the active region decayed from 6 March until 10 March, and by 9 March, all of the emerging fragments had coalesced together forming one large sunspot umbra.

On 10 March, the magnetic flux content of the active region began to increase once again, albeit this time without any new orbiting fragments. Little orbital motion of the previously-merged fragments was measured from 9–15 March (an average of 3° per day in the anticlockwise direction). This measured motion appears to be noisy, consisting of small values ≲1° per hour that vary between clockwise and anticlockwise directions (Fig. 5). However, when examining the merged sunspot by eye, there appear to still be two distinct fragments moving closely around each other that are not resolved by the fragment-tracking method (fragments B and C in panels e and f of Fig. 4). The magnetic flux fragment that arrived at the south of the leading sunspot at 12:00 UT on 8 March (fragment C) orbited almost 180° anticlockwise up to the north of the sunspot by 12:00 UT on 11 March – an average rotation of 60° per day. The CME studied in this section that occurred on 13 March 2012 followed this second phase of flux emergence.

4.2. 13 & 14 June 2012

4.2.1. Coronal evolution

On 13 June 2012 at 13:00 UT a CME erupted from NOAA AR 11504 in association with an M1.2 GOES class flare. The CME was seen in white-light by LASCO C2 at 14:36 UT (Fig. 6a). A hot EUV plasma emission feature (plasmoid) started to develop around 11:30 UT and it then rose above the flare arcade (that developed into the arcade of the M1.2 flare) during the 2 hours leading up to its eruption as the CME (Fig. 6b). At 12:14 UT, the underside of the slowly-rising plasmoid is estimated to be at a height of 173 ± 26 Mm above the centre of the active region’s photospheric PIL (hx = 193 ± 32 Mm, hy = 153 ± 21 Mm). An EUV dimming and hooked flare ribbons that formed in the east and west parts of the active region during the eruption infer the footpoint locations of the flux rope (Figs. 6c and 6d).

thumbnail Fig. 6.

Observations of the 13 June 2012 eruption. (a) White-light CME (indicated by the white arrow) observed by LASCO C2. (b) A plasmoid (highlighted by a red dashed line) brightens above a flaring arcade (indicated by a red arrow) in the 131 Å channel of AIA. (c) An EUV dimming (marked by a yellow arrow) seen in base difference 211 Å images. (d) Hooked flare ribbons (indicated by two black arrows) seen in the 1600 Å AIA channel.

The previous CME from NOAA AR 11504 occurred on 10 June, seen from 07:30 UT in LASCO C2 images as a faint circular blob. An associated M1.3 flare began in the active region at ≈06:39 UT on 10 June. Between the CMEs on 10 June and 13 June, there were 11 flares or brightenings in NOAA AR 11504. All were fairly weak (≤C2.7 GOES class) and they illuminated loops that spanned the full width of active region as did the 13 June CME and flare. During seven of these flares, the 1600 Å channel of AIA showed homologous flare ribbons in similar locations to those seen at the onset of the CME on 13 June. These flares occurred on 11 June at 17:20 UT (no GOES class assigned), on 12 June: a C2.0 flare at 06:38 UT, a C1.6 at 09:00 UT, a C.1.1 at 12:15 UT, and on 13 June: a C2.2 flare at 03:15 UT, a C1.3 at 05:40 UT, and a C2.7 at 09:13 UT. We suggest these seven homologous confined flares were related to reconnection that built the pre-eruptive flux rope, with the first occurring at ≈17:20 UT on 11 June (≈42 hours before the CME on 13 June).

NOAA AR 11504 also produced a CME on 14 June in association with an M1.9 GOES class flare. As described in Paper I, this CME was seen in LASCO C2 on 14 June 2012 at 14:12 UT (Fig. 7a). A sigmoid brightened in the active region above the arcade (that became the loops of the M1.9 class flare) at least 2 hours before erupting (Fig. 7b). The middle of the sigmoid is estimated to be at an altitude of 89 ± 30 Mm above the photosphere at 12:24 UT (hx = 85 ± 37 Mm, hy = 94 ± 22 Mm). For comparison, the underside of the extrapolated flux rope in Paper II was 45 Mm above the photosphere at the same time, and its top was at 150 Mm. Whilst the estimated sigmoid height is greater than that of the modelled flux rope’s underside (where the sigmoid is expected to be), it is lower than the highest point, suggesting that the method used in this study can provide estimates of flux rope height that are at least of comparable orders of magnitude. EUV dimmings and hooked flare ribbons show the footpoint locations of the flux rope during the eruption (Figs. 7c, and 7d). The observations suggest that an HFT flux rope formed in the active region, which is further supported by spectroscopic measurements of coronal plasma and an extrapolated magnetic field model (Papers I and II).

thumbnail Fig. 7.

Observations of the 14 June 2012 eruption. (a) White-light CME (indicated by the white arrow) observed by LASCO C2. (b) A sigmoid (highlighted by the dashed red line) brightens above a flaring arcade (indicated by the red arrow) in the 131 Å channel of AIA. (c) Twin EUV dimmings (marked by two yellow arrows) seen in base difference 131 Å AIA images. (d) Hooked flare ribbons (indicated by two black arrows) seen in the 1600 Å channel of AIA.

There were six confined brightening events (flares and EUV flux rope signatures) that spanned the active region PIL between the CMEs on 13 and 14 June, with the flares ranging from GOES class C1.7–C5.0 (on 13 June at 19:17 UT and on 14 June at 03:10 UT, 04:44 UT, 07:31 UT, 10:45 UT, and 11:05 UT). In each case, flare ribbons were observed in the AIA 1600 Å channel that resemble those seen during the 14 June CME. Therefore, we infer that these flares evidence the formation of the flux rope that erupted on 14 June throughout an 18-hour period.

4.2.2. Photospheric evolution

NOAA AR 11504 was close to the centre of the solar disc when it erupted on 13 and 14 June 2012. This gives an excellent viewpoint from SDO for examining the corona and photosphere with minimal projection effects. As detailed in Paper I, the active region contained two existing sunspots as it rotated on to the solar disc on 8 June, and flux emergence was clearly taking place by 9 June. The emerging flux showed right-handed magnetic tongues, indicating that the emerging flux tube had right-handed twist (positive chirality). The magnetic flux emergence occurred in “episodes”, with a major episode that began on 11 June and continued until 15 June in which distinct fragments of magnetic flux emerged one after the other and moved towards and clockwise around the pre-existing sunspots of the same polarities (described in detail in Sect. 3.2 of Paper I).

The strongest orbital motion was observed around the leading sunspot (see Fig. 8), and in this section we refer to certain fragments of magnetic flux around this sunspot with letters A–D, as in Fig. 8. Little-to-no significant orbital motion was measured on 9 and 10 June (2° in 24 hours between the pre-existing fragment A and the recently-emerged B). Between 11 June 06:00 UT and 12 June 12:00 UT, a fragment of emerged positive magnetic flux (fragment C) moved clockwise around the pre-existing positive sunspot (the now-merged fragments A and B) by ≈60° (an average orbital motion of 48° per day). This emerging fragment (C) grew to become the dominant umbra in the leading sunspot, and between 12 June 18:00 UT and 15 June 18:00 UT, another fragment (D) travelled ≈186° around fragment C (an average orbital motion of ≈62° per day). These quantified orbits are given in Table 2 and also represented in the top panel of Fig. 9.

thumbnail Fig. 8.

Clockwise motion of newly-emerged flux around the pre-existing positive (leading) sunspot in NOAA AR 11504. Tracked fragments of flux are labelled A–D. In each panel, vectors are drawn to connect the flux-weighted centroids of two orbiting fragments. This figure is available as an online movie.

thumbnail Fig. 9.

Top: measured orbital motion of chosen magnetic flux fragments in 6-hour intervals. The colours correspond to different choices of fragments (see Fig. 8). The black curve represents the orbit angles smoothed with a three-point moving average. Middle: full-disc integrated GOES soft X-ray lightcurve. Bottom: the evolution of magnetic flux in NOAA AR 11504, made using the radial magnetic field component, Br, of the HMI SHARP data series and smoothed with a 24-hour moving average. Vertical dashed lines indicate the times of the CME onsets on 10, 13, and 14 June 2012, and the thin, light-blue vertical bars show the timings of confined flares associated with flux rope formation described in Sect. 4.2.1.

Figure 9 shows that the two eruptive M-class flares from the active region occurred during the second period of strongest orbiting. Between the CMEs on 10 and 13 June, a net clockwise orbiting of 74° was observed, and between the CMEs on 13 and 14 June, 145° of clockwise orbiting occurred. This correlation between stronger orbiting and shorter formation time between eruptions supports the hypothesis that the orbital motions bring coronal loops together and facilitate the magnetic reconnection that builds the flux ropes.

4.3. 8 October 2012

4.3.1. Coronal evolution

A CME erupted from NOAA AR 11585 at ≈20:30 UT on 8 October 2012. The eruption was linked to a slow white-light CME seen in STEREO-A coronagraphs: COR1 from 21:15 UT, and COR2 from 23:09 UT (Fig. 10a). A faint arcade brightened in the active region as the eruption proceeded, but it was not designated a GOES flare class. Loops that extended from the east and west of the active region appeared strongly curved in the hours before the eruption, and from 19:30 UT, both of the curved ends began to expand, suggesting an expanding sigmoid in the active region (Fig. 10b). The interpretation of a continuous sigmoid in the active region is supported by X-ray observations2 taken on 7 October 2012 by the Hinode X-Ray Telescope (Golub et al. 2007). Furthermore, hot, rising plasma emission was seen near the centre of the active region from 20:19 UT, suggesting that a flux rope had formed in the active region. The hot feature was at an estimated altitude of 98 ± 25 Mm above the centre of the photospheric PIL at 20:19 UT (hx = 112 ± 39 Mm, hy = 85 ± 12 Mm). During the eruption, twin EUV dimmings and a hooked flare ribbon developed in the active region, indicating that the erupting flux rope had its footpoints in the east and west sides of the region (Figs. 10c and 10d).

thumbnail Fig. 10.

Observations of the 8 October 2012 eruption. (a) White-light CME (indicated by the white arrow) observed by STEREO COR2-A. (b) An inferred sigmoid (indicated by the dashed red line) in the 131 Å channel of AIA. (c) Twin EUV dimmings (marked by the two yellow arrows) seen in base difference 193 Å images. (d) One hooked flare ribbon (indicated by the black arrow) seen in the 1600 Å AIA channel.

No previous CMEs were observed to have originated from NOAA AR 11585 during its passage across the solar disc. Therefore, the entire disc-passage of the region from its appearance on 1 October until the CME on 8 October was analysed to search for flux rope formation signatures. Five confined flares occurred in the centre of NOAA AR 11585 in the same place as the subsequent eruption arcade, ranging from GOES class B3.2–B5.3 (three on 6 October at 05:23 UT, 07:06 UT, and 11:19 UT, and two on 7 October at 03:16 UT and 12:22 UT), all with similar flare ribbons that resemble the ones seen during the 8 October CME. This suggests the confined flares were involved in the formation of the flux rope that erupted on 8 October, and therefore they indicate the times at which the flux rope formed. In this way, formation started with a confined flare at 05:23 UT on 6 October, ≈63 hours before the 8 October CME.

4.3.2. Photospheric evolution

The leading sunspot of NOAA AR 11585 was forming as it rotated on to the solar disc and a small trailing (negative) spot briefly coalesced on 3 October before dispersing. The region exhibited left-handed magnetic tongues, implying the emerging flux had left-handed twist (negative chirality). Flux emergence had ceased by 4 October, and the region began to decay (see Fig. 11).

thumbnail Fig. 11.

Top: measured orbital rotation of a chosen magnetic flux fragment around another in 6-hour intervals (see Fig. 12). The black curve represents the orbit angles smoothed with a three-point moving average. Middle: full-disc integrated GOES soft X-ray lightcurve. Bottom: the evolution of magnetic flux in NOAA AR 11585, made using the radial magnetic field component, Br, of the HMI SHARP data series and smoothed with a 24-hour moving average. The vertical dashed line indicates the time of the CME onset on 8 October 2012, and the thin, light-blue vertical bars show the timings of confined flares associated with flux rope formation described in Sect. 4.3.1.

The positive sunspot was comprised of three distinctly-emerged umbral fragments within a single penumbra that orbited around each other during the week before the eruption (labelled A–C in Fig. 12). The two fragments that showed the strongest orbiting relative to each other were chosen (A and B). It is worth emphasising that the majority of the orbiting motions observed in this event occurred when the active region was decaying, in other words, after flux emergence had finished (or at least after the horizontal part of the emerging flux tube had crossed the photosphere).

thumbnail Fig. 12.

Anticlockwise motion of newly-emerged flux around the pre-existing positive (leading) sunspot in NOAA AR 11585. Tracked fragments of flux are labelled A–C. In each image, vectors are drawn either to connect the flux-weighted centroids of two orbiting fragments (panels a-c) or to best-fit the major axis of merged fragments (panels d-f). This figure is available as an online movie.

Between 3 October 06:00 UT and 6 October 00:00 UT (2.75 days), the chosen fragments orbited by 111° (an average of ≈40° per day). Between 6 October 12:00 UT and 8 October 00:00 UT (1.5 days), the chosen fragment orbited around the other by 41° (an average of ≈27° per day) and merged into one umbra (labelled ABC in Fig. 12). Between 8 October 12:00 UT and 9 October 00:00 UT (0.5 days), the merged umbra rotated by 7° (an average of 14° per day). These motions are quantified in Table 2 and visualised in Fig. 11 along with the GOES X-ray activity of the Sun during these times.

From 3–7 October, the active region was decaying, but it entered a new phase of emergence from 7 October until past the time of eruption. 159° of anticlockwise orbiting was measured from 3–9 October – an average of ≈26.5° per day. The orbiting of fragments in this event was weaker than in some others, but occurred continuously over a long time. This could explain the lack of flux rope ejections from this active region, with weaker motions taking a longer time to build a flux rope.

4.4. 14 July 2017

4.4.1. Coronal evolution

On 14 July 2017, a CME erupted from NOAA AR 12665 at ∼01:00 UT, accompanied by an M2.4 flare. The corresponding white-light CME was first seen by LASCO C2 at 01:36 UT (Fig. 13a). The active region appeared faintly sigmoidal several hours before the eruption in the 131 Å channel of AIA (Fig. 13b). A relatively faint flare arcade began to appear beneath the sigmoid from ∼00:00 UT, suggesting an HFT flux rope was present in the active region before the CME. The sigmoid began to expand and erupt at 00:30 UT, and the flare arcade brightened and grew in to a cusp shape. At 00:30 UT, the peak of the flare arcade cusp was at an estimated 71 ± 20 Mm above the photosphere (hx = 59 ± 8 Mm, hy = 82 ± 32 Mm). During the eruption, observed twin EUV dimmings and flare ribbons imply that the erupting flux rope had footpoints in the east and west parts of the active region (Figs. 13c and 13d).

thumbnail Fig. 13.

Observations of the 14 July 2017 eruption. (a) White-light CME (indicated by the white arrow) observed by LASCO C2. (b) Curved loops (indicated by the dashed red lines) that may belong to a sigmoid above the observed flare arcade (indicated by the red arrow) in the 131 Å channel of AIA. (c) Twin EUV dimmings (marked by the two yellow arrows) seen in base difference 211 Å images. (d) Flare ribbons (indicated by two black arrows) seen in the 1600 Å channel of AIA.

There were no CMEs observed from NOAA AR 12665 before the one on 14 July, so we look back over the full disc-passage of the region for signs of flux rope formation. Shortly after a confined M1.3 flare on 9 July, there were two confined C-class flares: a C1.6 flare from 06:15–06:40 UT and a C1.2 flare from 07:28–07:49 UT. Around the times of these C-flares (between 05:35–06:31 UT and 07:51–08:47 UT), faint emission signatures were seen towards the eastern edge of the active region in the 131 Å and 94 Å AIA channels only (suggesting plasma at T ∼ 10 MK). We interpret these signatures as part of a hot flux rope, brightening in association with the confined flares. Since there were no successful eruptions from 9 July until the CME on 14 July, the hot flux rope evidenced on 9 July may be the same one that erupted on 14 July. This would suggest flux rope formation had begun 115 hours before the CME, although it is difficult to support this by comparing flare ribbons between 9 and 14 July due to the dramatic photospheric evolution that occurred between those times (the strong orbiting described in Sect. 4.4.2). After the two C-flares described above, eleven additional confined B- or C-class flares occurred on 9 July, followed by a confined C1.5 flare on 10 July. No hot flux rope signatures were observed in association with these flares, so we do not link them to flux rope formation. The earliest flare ribbons that appeared homologous to the eruptive flare were seen during a confined B4.9 flare at 02:00 UT on 11 July, with one bright globular ribbon visible near the leading sunspot (and flux rope footpoint) and a second ribbon near the trailing polarity. Two more confined B-flares (at 04:40 UT on 11 July and 05:15 UT on 13 July) also showed similar homologous flare ribbons, suggesting they, too, were involved in forming the flux rope. In summary, there is evidence to suggest that two C-flares on 9 July, two B-flares on 11 July, and a B-flare on 13 July were involved in the formation of the flux rope that erupted on 14 July 2017.

4.4.2. Photospheric evolution

As the region rotated on to the disc on 5 July, it featured two small pre-existing sunspots; a leading positive spot and a trailing negative. Observed magnetic tongues demonstrate that the magnetic flux that emerged in to the active region was left-handed (negative chirality). Some very weak anticlockwise orbiting can be observed around the trailing sunspot, but as with the other events, here we focus on the leading sunspot.

A fragment of magnetic flux emerged on 6 July (fragment B in Fig. 14) and moved towards the leading sunspot before rotating anticlockwise around the pre-existing sunspot umbra (fragment A). Between 8 July 00:00 UT and 10 July 00:00 UT (2 days), the emerging fragment orbited around the other by 92° (an average of 46° per day). The strongest orbiting during this period was seen on the evening of 9 July, where the fragments sheared past each other from 18:00–00:00 UT for a rotation of 54°. Between 10 July 06:00 UT and 14 July 00:00 UT (3.75 days), fragments A and B merged together, and an orbit of 47° is inferred from the continued rotation of the sunspot (an average of 12.5° per day). Between 14 July 00:00 UT and 15 July 00:00 UT (1 day), a fragment broke away from the merged umbra (labelled fragment C) and continued to orbit the sunspot by 43°. These values are collected in Table 2, and presented in Fig. 15 with the corresponding X-ray activity measured by GOES.

thumbnail Fig. 14.

Anticlockwise motion of newly-emerged flux around the pre-existing positive (leading) sunspot in NOAA AR 12665. Tracked fragments of flux are labelled A–C. In each image, vectors are drawn either to connect the flux-weighted centroids of two orbiting fragments (panels a, b, and f) or to best-fit the major axis of merged fragments (panels c, d, and e). This figure is available as an online movie.

thumbnail Fig. 15.

Top: measured orbital motion of chosen magnetic flux fragments in 6-hour intervals. The colours correspond to different choices of fragments (see Fig. 14). The black curve represents the orbit angles smoothed with a three-point moving average. Middle: full-disc integrated GOES soft X-ray lightcurve. Bottom: the evolution of magnetic flux in NOAA AR 12665, made using the radial magnetic field component, Br, of the HMI SHARP data series and smoothed with a 24-hour moving average. The vertical dashed line indicates the time of CME onset on 14 July 2017, and the thin, light-blue vertical bars show the timings of confined flares associated with flux rope formation described in Sect. 4.4.1.

From 8–14 July, 138° of anticlockwise orbiting was measured – an average of ≈20° per day. The average daily orbiting in this active region is weaker than in some events, but since there were no observed prior CMEs from the region, and the flux rope may have been forming since as early as 9 July, the flux rope may have formed gradually over several days. The main emerging sunspot fragment encountered the pre-existing sunspot on 8 July and showed strong shearing motion clockwise around it throughout 9 July (90° in 24 hours). Most of the strong flares from the active region occurred during this time, with an M2.4 flare, 4 C-flares, and 16 B-flares.

5. Discussion

Across Papers I and II, a methodology is detailed by which a combination of observations and modelling enables the identification of hot flux ropes when they are on the solar disc (as opposed to limb detection). This allowed for the investigation of the formation mechanism of hot flux ropes.

This work tests the hypothesis that the photospheric orbiting motions associated with emerging magnetic flux are the trigger for building HFT flux ropes via reconnection in the corona (see Fig. 16), as was suggested in Papers I and II. We have shown that this mechanism is responsible for building five hot flux ropes with the observational signatures of hyperbolic flux tube (HFT) configurations that then erupted as CMEs on 13 March 2012 (from NOAA active region 11429), 13 June 2012 (from NOAA active region 11504), 14 June 2012 (from NOAA active region 11504), 8 October 2012 (from NOAA active region 11585), and 14 July 2017 (from NOAA active region 12665). The radial magnetic flux of each active region (between ∼3 × 1021 Mx and 3 × 1022 Mx) places them in the large active region category according to Schrijver & Zwaan (2000, presented in Table 1 of van Driel-Gesztelyi & Green 2015).

thumbnail Fig. 16.

Orbiting motions of a magnetic flux fragment around a pre-existing sunspot. (a) A magnetic bipole (red) emerges beneath the pre-existing bipole field of sunspots (blue). (b) A fragment of the emerging flux moves towards and “orbits” around the pre-existing sunspot of the same polarity, wrapping field lines around each other. (c) Component magnetic reconnection occurs in the corona. (d) The products of this reconnection are sheared (twisted) field lines in the form of an overlying flux rope and an underlying arcade.

The classification as an HFT rope comes from the estimated height of the underside of the observed hot flux rope or sigmoid and the presence of lower-lying arcade field as determined from EUV images. As detailed in Table 2, the heights of the flux ropes that erupted as CME events on 13 March 2012, 13 June 2012, 14 June 2012, 8 October 2012 and 14 July 2017 were estimated to be 56 ± 5 Mm, 173 ± 26 Mm, 89 ± 30 Mm, 98 ± 25 Mm and 71 ± 20 Mm above the photosphere, respectively. However, the 13 March 2012 estimate may be strongly affected by projection effects, and the height determination in the 13 June 2012 case was made during the slow-rise phase of the structure, potentially leading to an overestimation of the height. Furthermore, the height of the underside of the flux rope that erupted on 14 July 2017 was determined from the highest point of the flare arcade associated with the event, which is expected to be just beneath the bottom of the flux rope, and the height measurement was also made during the very early stage of the eruption. The height estimates are of similar order to the height of the flux rope modelled in Paper II, which had its underside and top at 45 Mm and 150 Mm above the photosphere, respectively, and with the observational study of Patsourakos et al. (2013) who found a hot flux rope that formed at an altitude of 80 Mm during a confined flare, rose slowly, and had an altitude of 138 Mm at the start of the impulsive CME acceleration phase. The heights place the ropes in the β <  1 region of the solar atmosphere (see Fig. 3 of Gary 2001), implying that plasma does not contribute to the stability of these HFT flux ropes in the same way that it would if their undersides were line-tied to dense photospheric or chromospheric plasma in a BPS configuration or if they contained dense filament plasma (Jenkins et al. 2018, 2019). If perturbed, reconnection could occur at the HFT under these flux ropes, aiding the eruption of the structure.

The main finding in this work is that the combination of orbital motions of photospheric magnetic field and magnetic reconnection in the corona is able to act collectively as a CME trigger mechanism, forming pre-eruptive HFT ropes. In addition to the orbiting motions, all regions produced confined flares that could be termed homologous: the confined flares occurred along the same polarity inversion line, had broadly the same spatial extent as observed in the early stages of the eventual eruptive flares, and featured flare ribbons that resembled those seen later during the CMEs. This indicates that the formation process of coronal HFT flux ropes in this study is systematically related to magnetic reconnection driven by the photospheric orbiting motions.

The confined flare observations can therefore be used to probe the timescale over which each flux rope formed; each reconnection event transforming the configuration from sheared arcade to flux rope and feeding more axial flux into the structure. Homologous C-class flares before the 13 March 2012 CME (on 12 March at 22:20 UT and 13 March at 06:55 UT) suggest a flux rope was forming in the active region for ≈19 hours before it erupted. In the case of the 13 June 2012 CME, homologous, confined C-flares occurred on 11, 12, and 13 June, suggesting a flux rope formation time of up to 42 hours. Another CME occurred from the same active region on 14 June, with relevant C-class flares spread over an ≈18 hour period. Clear signatures of the EUV flux rope in the from of a sigmoid were observed 2 hours before the eruption (see Paper I) indicating its presence as a coherent structure by this time. No other CMEs from the same source region were observed before the one on 8 October 2012. B-class flares were seen in the active region from the morning of 6 October onwards with homologous flare ribbons to those seen during the CME, suggesting a flux rope formation timescale of 63 hours. However, all observational signatures in this active region were weak so there is uncertainty in this proposed timescale. Finally, faint flux rope signatures were observed in the hottest AIA channels (131 Å, 94 Å) during confined C-class flares on 9 July 2017, 115 hours before the CME on 14 July.

In summary, the timescale over which the HFT rope in each active region is likely being built through successive reconnection events varies significantly from 18 to 115 hours. The key point here is that several episodes of magnetic reconnection in the corona may be needed to build the HFT structure that becomes eruptive, in an analogous way to the ongoing episodes of magnetic reconnection in the photosphere or chromosphere that build BPS ropes at lower altitudes (van Ballegooijen & Martens 1989).

The observed photospheric motions were quantified using the method described in Sect 3.1. Across the studied intervals presented in Table 2, rates of averaged daily orbiting were found between 12–48° per day. The highest rate of orbiting measured between any two subsequent time steps was 54.1°, measured over a 6-hour interval beginning at 18:00 UT on 9 July 2017 in NOAA active region 12665 (see Fig. 15). We also note the strong rate of orbiting observed in NOAA active region 11504 between the CMEs on 13 and 14 June 2012. In this ≈24 hour period, 103° of orbiting was observed (measured from noon on 13 June until noon on 14 June). This high rate of orbiting may explain how a flux rope was able to form in NOAA AR 11504 so quickly after the previous CME the day before. However, when considering all of the events together, there are no apparent correlations between formation time and average orbit rate, total orbit angle, or active region flux. Future work could examine a larger sample of active regions to develop a clearer sense of possible relationships between such parameters.

The conjecture that the photospheric orbiting motions are ultimately responsible for the flux rope formation is further tested by comparing the direction of these motions with the handedness of the flux rope. The observations show that there is consistency between the handedness of the emerging magnetic flux, the direction of the orbiting motion, and the handedness of the flux ropes that form. Active regions with emerging right-handed (left-handed) twist exhibited dominantly clockwise (anticlockwise) orbiting motions in the photosphere, and right-handed (left-handed) flux ropes formed in the corona. As an aside, we note here that active region NOAA 11429 (March 2012) does not conform to the expected Hale orientation for the northern hemisphere in solar cycle 24, and there is still agreement between emerging twist, orbiting direction, and flux rope handedness in this region. Previous studies have found that sunspot rotation can play a role in building the twist of forming flux ropes (Vemareddy et al. 2016; Yan et al. 2018). Similarly, in this study, once orbiting fragments had coalesced, their continuing orbital motions at times resembled the scenario of a rotating sunspot, giving further support to the scenario that twisted magnetic flux ropes can be formed in these active regions. Furthermore, simulations show that rotating sunspots can cause the inflation of overlying magnetic fields in active regions, facilitating the onset of the ideal torus instability as a CME driver (Török et al. 2013).

The origin of the orbiting motions seen in the photosphere of the active regions under study here is an open and interesting question, which may be related to observations of bodily rotation of sunspots. Brown et al. (2003) found sunspot rotations of 40 ° −60° over 3–5 days, and posed two possible explanations for sunspot rotation: photospheric flows that lead to localised proper motions, or flux emergence. Extending these explanations to this study, the first scenario is that the observed orbiting in the photosphere may be caused by one magnetic flux tube being physically moved around another by sub-photospheric flows with sufficient energy. The second scenario is that the orbiting may be an apparent motion caused by the emergence of two intertwined sub-photospheric flux tubes. Min & Chae (2009) studied a case of sunspot rotation and they suggest that this was an apparent motion caused by the emergence of twisted flux tube, but were unable to rule out the effect of a torque force from the solar interior. Contrarily, simulations by Sturrock et al. (2015) support the torque scenario for sunspot rotation and rule out the possibility of the rotation being an apparent effect.

To investigate whether the orbiting in this study is a product of flux emergence, the evolution of magnetic flux in each active region is examined. The confined flares observed in the regions that produced the CMEs on 13 March 2012, 13 June 2012, and 14 June 2012 occurred when the active regions were in phases of flux emergence. In the case of the 8 October 2012 CME, the confined flares took place during a phase of decaying magnetic flux, and before the 14 July 2017 CME, confined flaring occurred during periods of both flux emergence and decay. In the cases of all five CMEs, active region magnetic flux was increasing for at least 24 hours before eruption.

The lack of an observed increase in active region flux content during some of the confined flaring (coronal flux rope formation) does not necessarily rule out the scenario that the observed orbiting motions are a manifestation of an emerging flux tube. If we consider a twisted Ω−shaped flux tube emerging through the photosphere, we would expect an increase in observed magnetic flux as the tube’s apex rises into the solar atmosphere. However, once the apex has finished crossing the photospheric boundary, the mostly-vertical legs of the flux tube may continue to emerge with no discernible change in magnetic flux. It is therefore difficult to comment on the origin of the observed orbiting motions in this study. Regardless, the main conclusion is that magnetic reconnection occurred in the corona in every event examined here.

6. Conclusions

Five hot flux ropes have been identified that formed via magnetic reconnection in the corona with estimated heights ranging from ∼56 − 176 Mm, where the plasma β <  1. The reconnection that built the flux ropes occurred in sporadic bursts, as evidenced by confined solar flares, and the timings of these flares suggests the flux ropes began to form from around 18–115 hours before they erupted. Each confined flare event transformed sheared arcade field in to flux rope field and fed more flux to the growing structure.

In searching for the photospheric process(es) that caused magnetic reconnection to occur in the corona, it was found that all of the active regions exhibited newly-emerged magnetic flux fragments that moved towards and then orbited around pre-existing sunspots at some point during the time period of study. The studied flux ropes each had one leg rooted in the sunspots where the strongest orbiting was observed, and there was consistency between the handedness of the emerging magnetic flux, the sense of orbiting, and the handedness of the flux ropes that formed. Three left-handed sigmoids formed that were rooted in left-handed active regions that showed anticlockwise orbiting and two right-handed sigmoids formed in a right-handed active region that showed clockwise orbiting. Furthermore, the event with the shortest amount of time between successive eruptions featured the strongest rate of orbiting, with ≈100° of motion measured between the CMEs on 13 and 14 June 2012 (although this is not observed to be a general correlation when considering the other regions studied).

It is inferred that the orbiting motions observed in the photosphere enable magnetic reconnection in the corona, forming hot magnetic flux ropes in HFT configurations above sheared, flaring arcades. We propose that this mechanism be considered as a new trigger for CMEs that originate from hot flux ropes.


Acknowledgments

A. W. J. is supported by a European Space Agency (ESA) Research Fellowship, and acknowledges the support of the Leverhulme Trust and the Royal Society. L. M. G. is grateful to the Royal Society for the award of a Royal Society Research Fellowship that has supported this work. L. v. D. G. is partially funded under STFC consolidated grant number ST/N000722/1. L. v. D. G. acknowledges the Hungarian National Research, Development and Innovation Office Grant OTKA K-131508. G. V. acknowledges the support of the Leverhulme Trust Research Project Grant 2014-051, the support from the European Union’s Horizon 2020 research and innovation programme under Grant agreement No. 824135, and of the STFC Grant number ST/T000317/. Data courtesy of NASA/SDO and the AIA and HMI science teams. This research has made use of SunPy v1.1.0, an open-source and free community-developed solar data analysis Python package (SunPy Project 2020).

References

  1. Antiochos, S. K., DeVore, C. R., & Klimchuk, J. A. 1999, ApJ, 510, 485 [Google Scholar]
  2. Aulanier, G., Török, T., Démoulin, P., & DeLuca, E. E. 2010, ApJ, 708, 314 [Google Scholar]
  3. Baker, D., Brooks, D. H., Démoulin, P., et al. 2013, ApJ, 778, 69 [NASA ADS] [CrossRef] [Google Scholar]
  4. Baker, D., van Driel-Gesztelyi, L., Brooks, D. H., et al. 2019, ApJ, 875, 35 [NASA ADS] [CrossRef] [Google Scholar]
  5. Bobra, M. G., Sun, X., Hoeksema, J. T., et al. 2014, Sol. Phys., 289, 3549 [NASA ADS] [CrossRef] [Google Scholar]
  6. Boerner, P., Edwards, C., Lemen, J., et al. 2012, Sol. Phys., 275, 41 [NASA ADS] [CrossRef] [Google Scholar]
  7. Brown, D. S., Nightingale, R. W., Alexander, D., et al. 2003, Sol. Phys., 216, 79 [NASA ADS] [CrossRef] [Google Scholar]
  8. Brueckner, G. E., Howard, R. A., Koomen, M. J., et al. 1995, Sol. Phys., 162, 357 [NASA ADS] [CrossRef] [Google Scholar]
  9. Chintzoglou, G., Patsourakos, S., & Vourlidas, A. 2015, ApJ, 809, 34 [NASA ADS] [CrossRef] [Google Scholar]
  10. Démoulin, P., & Aulanier, G. 2010, ApJ, 718, 1388 [NASA ADS] [CrossRef] [Google Scholar]
  11. Demoulin, P., Henoux, J. C., Priest, E. R., & Mandrini, C. H. 1996, A&A, 308, 643 [NASA ADS] [Google Scholar]
  12. Domingo, V., Fleck, B., & Poland, A. I. 1995, Sol. Phys., 162, 1 [NASA ADS] [CrossRef] [Google Scholar]
  13. Gary, G. A. 2001, Sol. Phys., 203, 71 [NASA ADS] [CrossRef] [Google Scholar]
  14. Gerrard, C. L., Brown, D. S., Mellor, C., Arber, T. D., & Hood, A. W. 2003, Sol. Phys., 213, 39 [NASA ADS] [CrossRef] [Google Scholar]
  15. Golub, L., Deluca, E., Austin, G., et al. 2007, Sol. Phys., 243, 63 [NASA ADS] [CrossRef] [Google Scholar]
  16. Green, L. M., & Kliem, B. 2009, ApJ, 700, L83 [NASA ADS] [CrossRef] [Google Scholar]
  17. Green, L. M., Török, T., Vršnak, B., Manchester, W., & Veronig, A. 2018, Space Sci. Rev., 214, 46 [Google Scholar]
  18. Hood, A. W., & Priest, E. R. 1979, Sol. Phys., 64, 303 [NASA ADS] [CrossRef] [Google Scholar]
  19. James, A. W., Green, L. M., Palmerio, E., et al. 2017, Sol. Phys., 292, 71 [NASA ADS] [CrossRef] [Google Scholar]
  20. James, A. W., Valori, G., Green, L. M., et al. 2018, ApJ, 855, L16 [NASA ADS] [CrossRef] [Google Scholar]
  21. Janvier, M., Aulanier, G., Bommier, V., et al. 2014, ApJ, 788, 60 [NASA ADS] [CrossRef] [Google Scholar]
  22. Jenkins, J. M., Long, D. M., van Driel-Gesztelyi, L., & Carlyle, J. 2018, Sol. Phys., 293, 7 [Google Scholar]
  23. Jenkins, J. M., Hopwood, M., Démoulin, P., et al. 2019, ApJ, 873, 49 [NASA ADS] [CrossRef] [Google Scholar]
  24. Kaiser, M. L., Kucera, T. A., Davila, J. M., et al. 2008, Space Sci. Rev., 136, 5 [NASA ADS] [CrossRef] [Google Scholar]
  25. Karpen, J. T., Antiochos, S. K., & DeVore, C. R. 2012, ApJ, 760, 81 [NASA ADS] [CrossRef] [Google Scholar]
  26. Kliem, B., & Török, T. 2006, Phys. Rev. Lett., 96, 255002 [Google Scholar]
  27. Lemen, J. R., Title, A. M., Akin, D. J., et al. 2012, Sol. Phys., 275, 17 [Google Scholar]
  28. Liu, Y., Hoeksema, J. T., Scherrer, P. H., et al. 2012, Sol. Phys., 279, 295 [Google Scholar]
  29. Longcope, D. W., Fisher, G. H., & Pevtsov, A. A. 1998, ApJ, 507, 417 [NASA ADS] [CrossRef] [Google Scholar]
  30. López Fuentes, M. C., Demoulin, P., Mandrini, C. H., & van Driel-Gesztelyi, L. 2000, ApJ, 544, 540 [Google Scholar]
  31. Luoni, M. L., Démoulin, P., Mandrini, C. H., & van Driel-Gesztelyi, L. 2011, Sol. Phys., 270, 45 [NASA ADS] [CrossRef] [Google Scholar]
  32. Min, S., & Chae, J. 2009, Sol. Phys., 258, 203 [NASA ADS] [CrossRef] [Google Scholar]
  33. Moore, R. L., Sterling, A. C., Hudson, H. S., & Lemen, J. R. 2001, ApJ, 552, 833 [Google Scholar]
  34. Nindos, A., Patsourakos, S., Vourlidas, A., & Tagikas, C. 2015, ApJ, 808, 117 [NASA ADS] [CrossRef] [Google Scholar]
  35. Patsourakos, S., Vourlidas, A., & Stenborg, G. 2013, ApJ, 764, 125 [NASA ADS] [CrossRef] [Google Scholar]
  36. Pesnell, W. D., Thompson, B. J., & Chamberlin, P. C. 2012, Sol. Phys., 275, 3 [Google Scholar]
  37. Reeves, K. K., & Golub, L. 2011, ApJ, 727, L52 [NASA ADS] [CrossRef] [Google Scholar]
  38. Scherrer, P. H., Schou, J., Bush, R. I., et al. 2012, Sol. Phys., 275, 207 [Google Scholar]
  39. Schrijver, C. J., & Zwaan, C. 2000, Solar and Stellar Magnetic Activity (Cambridge University Press) [Google Scholar]
  40. Sturrock, Z., Hood, A. W., Archontis, V., & McNeill, C. M. 2015, A&A, 582, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. SunPy Project (Barnes, W. T., et al.) 2020, ApJ, 890, 68 [CrossRef] [Google Scholar]
  42. Syntelis, P., Gontikakis, C., Patsourakos, S., & Tsinganos, K. 2016, A&A, 588, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Temmer, M., Veronig, A. M., Kontar, E. P., Krucker, S., & Vršnak, B. 2010, ApJ, 712, 1410 [NASA ADS] [CrossRef] [Google Scholar]
  44. Titov, V. S., Hornig, G., & Démoulin, P. 2002, J. Geophys. Res. (Space Phys.), 107, 1164 [NASA ADS] [CrossRef] [Google Scholar]
  45. Török, T., & Kliem, B. 2005, ApJ, 630, L97 [Google Scholar]
  46. Török, T., Temmer, M., Valori, G., et al. 2013, Sol. Phys., 286, 453 [NASA ADS] [CrossRef] [Google Scholar]
  47. van Ballegooijen, A. A., & Martens, P. C. H. 1989, ApJ, 343, 971 [Google Scholar]
  48. van Driel-Gesztelyi, L., & Green, L. M. 2015, Liv. Rev. Sol. Phys., 12, [Google Scholar]
  49. van Tend, W., & Kuperus, M. 1978, Sol. Phys., 59, 115 [NASA ADS] [CrossRef] [Google Scholar]
  50. Vemareddy, P., Cheng, X., & Ravindra, B. 2016, ApJ, 829, 24 [CrossRef] [Google Scholar]
  51. Vourlidas, A., Lynch, B. J., Howard, R. A., & Li, Y. 2013, Sol. Phys., 284, 179 [NASA ADS] [Google Scholar]
  52. Wang, R., Liu, Y. D., Yang, Z., & Hu, H. 2014, ApJ, 791, 84 [NASA ADS] [CrossRef] [Google Scholar]
  53. Yan, X. L., Qu, Z. Q., Kong, D. F., & Xu, C. L. 2012, ApJ, 754, 16 [NASA ADS] [CrossRef] [Google Scholar]
  54. Yan, X. L., Wang, J. C., Pan, G. M., et al. 2018, ApJ, 856, 79 [NASA ADS] [CrossRef] [Google Scholar]
  55. Yardley, S. L., Green, L. M., van Driel-Gesztelyi, L., Williams, D. R., & Mackay, D. H. 2018, ApJ, 866, 8 [Google Scholar]

Movies

Movie of Fig. 4 Access here

Movie of Fig. 8 Access here

Movie of Fig. 12 Access here

Movie of Fig. 14 Access here

All Tables

Table 1.

CME source regions.

Table 2.

Orbital motions of umbral flux fragments relative to previously-emerged umbrae and the observed active region chiralities.

All Figures

thumbnail Fig. 1.

Identifying magnetic flux fragments. Left: contours are drawn at fixed intensity levels on HMI continuum images to enclose individual umbral fragments. Right: the continuum contours are overlaid on the radial component of magnetic field. Positive polarity field is shown in white, and negative in black. The flux-weighted centroid of each fragment is shown in red and is assigned a unique identifying number.

In the text
thumbnail Fig. 2.

Illustration of the method used to track the relative orbital motion of magnetic flux fragments.

In the text
thumbnail Fig. 3.

Observations of the 13 March 2012 eruption. (a) A white-light CME (indicated by the white arrow) seen by LASCO C2. (b) AIA 131 Å image of a sigmoid (highlighted by a dashed red line) brightening above a flaring arcade (indicated by a red arrow). (c) AIA 211 Å base difference image showing twin EUV dimmings (marked by two yellow arrows). (d) AIA 1600 Å image showing hooked flare ribbons (indicated by two black arrows).

In the text
thumbnail Fig. 4.

Anticlockwise motion of newly-emerged flux around the pre-existing positive (leading) sunspot in NOAA AR 11429. Tracked fragments of flux are labelled A–C. In each image, vectors are drawn either to connect the flux-weighted centroids of two orbiting fragments (panels a–d) or to best-fit the major axis of merged fragments (panels e and f). In the bottom two panels, red circles highlight the position of a fragment (C) whose ≈180° orbit from the south of the sunspot to the north between panels was not successfully tracked using the vector method. This figure is available as an online movie.

In the text
thumbnail Fig. 5.

Top: measured orbital motion of chosen magnetic flux fragments in 6-hour intervals. The colours correspond to different choices of fragments (see Fig. 4), except for the purple block which corresponds to the average orbiting estimated by-eye that was not detected by the automated method. The black curve represents the orbit angles smoothed using a three-point moving average. Middle: full-disc integrated GOES soft X-ray lightcurve. Bottom: the evolution of magnetic flux in NOAA AR 11429, made using the radial magnetic field component, Br, of the HMI SHARP data series and smoothed with a 24-hour moving average. Vertical dashed lines indicate the times of the CME onsets on 10 and 13 March 2012, and the thin, light-blue vertical bars show the timings of confined flares associated with flux rope formation described in Sect. 4.1.1.

In the text
thumbnail Fig. 6.

Observations of the 13 June 2012 eruption. (a) White-light CME (indicated by the white arrow) observed by LASCO C2. (b) A plasmoid (highlighted by a red dashed line) brightens above a flaring arcade (indicated by a red arrow) in the 131 Å channel of AIA. (c) An EUV dimming (marked by a yellow arrow) seen in base difference 211 Å images. (d) Hooked flare ribbons (indicated by two black arrows) seen in the 1600 Å AIA channel.

In the text
thumbnail Fig. 7.

Observations of the 14 June 2012 eruption. (a) White-light CME (indicated by the white arrow) observed by LASCO C2. (b) A sigmoid (highlighted by the dashed red line) brightens above a flaring arcade (indicated by the red arrow) in the 131 Å channel of AIA. (c) Twin EUV dimmings (marked by two yellow arrows) seen in base difference 131 Å AIA images. (d) Hooked flare ribbons (indicated by two black arrows) seen in the 1600 Å channel of AIA.

In the text
thumbnail Fig. 8.

Clockwise motion of newly-emerged flux around the pre-existing positive (leading) sunspot in NOAA AR 11504. Tracked fragments of flux are labelled A–D. In each panel, vectors are drawn to connect the flux-weighted centroids of two orbiting fragments. This figure is available as an online movie.

In the text
thumbnail Fig. 9.

Top: measured orbital motion of chosen magnetic flux fragments in 6-hour intervals. The colours correspond to different choices of fragments (see Fig. 8). The black curve represents the orbit angles smoothed with a three-point moving average. Middle: full-disc integrated GOES soft X-ray lightcurve. Bottom: the evolution of magnetic flux in NOAA AR 11504, made using the radial magnetic field component, Br, of the HMI SHARP data series and smoothed with a 24-hour moving average. Vertical dashed lines indicate the times of the CME onsets on 10, 13, and 14 June 2012, and the thin, light-blue vertical bars show the timings of confined flares associated with flux rope formation described in Sect. 4.2.1.

In the text
thumbnail Fig. 10.

Observations of the 8 October 2012 eruption. (a) White-light CME (indicated by the white arrow) observed by STEREO COR2-A. (b) An inferred sigmoid (indicated by the dashed red line) in the 131 Å channel of AIA. (c) Twin EUV dimmings (marked by the two yellow arrows) seen in base difference 193 Å images. (d) One hooked flare ribbon (indicated by the black arrow) seen in the 1600 Å AIA channel.

In the text
thumbnail Fig. 11.

Top: measured orbital rotation of a chosen magnetic flux fragment around another in 6-hour intervals (see Fig. 12). The black curve represents the orbit angles smoothed with a three-point moving average. Middle: full-disc integrated GOES soft X-ray lightcurve. Bottom: the evolution of magnetic flux in NOAA AR 11585, made using the radial magnetic field component, Br, of the HMI SHARP data series and smoothed with a 24-hour moving average. The vertical dashed line indicates the time of the CME onset on 8 October 2012, and the thin, light-blue vertical bars show the timings of confined flares associated with flux rope formation described in Sect. 4.3.1.

In the text
thumbnail Fig. 12.

Anticlockwise motion of newly-emerged flux around the pre-existing positive (leading) sunspot in NOAA AR 11585. Tracked fragments of flux are labelled A–C. In each image, vectors are drawn either to connect the flux-weighted centroids of two orbiting fragments (panels a-c) or to best-fit the major axis of merged fragments (panels d-f). This figure is available as an online movie.

In the text
thumbnail Fig. 13.

Observations of the 14 July 2017 eruption. (a) White-light CME (indicated by the white arrow) observed by LASCO C2. (b) Curved loops (indicated by the dashed red lines) that may belong to a sigmoid above the observed flare arcade (indicated by the red arrow) in the 131 Å channel of AIA. (c) Twin EUV dimmings (marked by the two yellow arrows) seen in base difference 211 Å images. (d) Flare ribbons (indicated by two black arrows) seen in the 1600 Å channel of AIA.

In the text
thumbnail Fig. 14.

Anticlockwise motion of newly-emerged flux around the pre-existing positive (leading) sunspot in NOAA AR 12665. Tracked fragments of flux are labelled A–C. In each image, vectors are drawn either to connect the flux-weighted centroids of two orbiting fragments (panels a, b, and f) or to best-fit the major axis of merged fragments (panels c, d, and e). This figure is available as an online movie.

In the text
thumbnail Fig. 15.

Top: measured orbital motion of chosen magnetic flux fragments in 6-hour intervals. The colours correspond to different choices of fragments (see Fig. 14). The black curve represents the orbit angles smoothed with a three-point moving average. Middle: full-disc integrated GOES soft X-ray lightcurve. Bottom: the evolution of magnetic flux in NOAA AR 12665, made using the radial magnetic field component, Br, of the HMI SHARP data series and smoothed with a 24-hour moving average. The vertical dashed line indicates the time of CME onset on 14 July 2017, and the thin, light-blue vertical bars show the timings of confined flares associated with flux rope formation described in Sect. 4.4.1.

In the text
thumbnail Fig. 16.

Orbiting motions of a magnetic flux fragment around a pre-existing sunspot. (a) A magnetic bipole (red) emerges beneath the pre-existing bipole field of sunspots (blue). (b) A fragment of the emerging flux moves towards and “orbits” around the pre-existing sunspot of the same polarity, wrapping field lines around each other. (c) Component magnetic reconnection occurs in the corona. (d) The products of this reconnection are sheared (twisted) field lines in the form of an overlying flux rope and an underlying arcade.

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.