P. Vemareddy

Heliophysicist
Indian Institute of Astrophysics

Research

Activity of the Sun is defined by eruptive events like Coronal mass ejections (CMEs) and flares occurring in the ARs. Solar flare is a sudden release of energy resulting in acceleration of particles and plasma flows and heat. Large flares can be accompanied by highly energetic particles that pose hazards ranging from excessive levels of harmful radiation to disruption and/or destruction of satellites in space. CMEs are energetic events in which the plasma that was initially entrapped in closed coronal magnetic field lines is ejected into interplanetary space. They propagate through the interplanetary medium with speeds up to a maximum of 3000 km/s. Sometimes CMEs precede flares and sometimes not, which are now referred to eruptive and confined flares, to bridge the gap between flares and CMEs.

It is believed that these events occur due to release of energy which is pre-stored in the magnetic field in stressed form by constantly straining plasma motions. When the AR magnetic system reaches a critical point of stored energy, it releases as an unavoidable eruption. Major fundamental questions in the field of MHD include how these events trigger and what mechanism drives them. On this line, my research interests on solar eruptions include flux rope formation and CME initiation, Sun-Earth connection of flux rope CMEs, quasi-static evolution of AR magnetic structure and topological study, data driven MHD simulations of AR magnetic structure and global Sun. Recent research projects are briefly described below.

Figure 1 Rising motion of flux rope as seen in AIA 304 (Top row), Hα (Middle row), and STEREO-A/EUVI 304Å (bottom row) observations, respectively. White thick curves represent the observed rising flux rope. In panel (d), the overlaid red (blue) contours represent positive (negative) polarities. Source: Vemareddy et al (2012)
During the eruptions, a commonly observed physical structure is magnetic flux rope (Burlaga et al 1988, Lepping et al 1990). It is a volumetric current channel with wrapped fieldlines around its central axis. Its preexistence, formation and role in eruptions are not well understood because it is believed that magnetic reconnection plays the main role and flux rope is treated as secondary.

Recent direct evidences of the flux rope come into existence as a conspicuous channel structure in the inner corona of the AR before and during a solar eruption (Zhang et al. 2012, Vemareddy et al 2012, and see Figure 1). This channel initially appears as a twisted and writhed sigmoidal structure in high temperature passbands. The channel evolves toward a semi-circular shape in the slow rise phase and then erupts upward rapidly in the impulsive acceleration phase, producing the front-cavity core components of the resulting CME (Zhang et al. 2012). The role of this hot channel in the eruption process is similar to that of flux rope in many numerical models (e.g., Aulanier et al. 2010; Fan & Gibson 2007) reproducing key features of eruptions. Filament channels are also regarded as flux ropes, which were recently recognized by Atmospheric Imaging Assembly onboard Solar Dynamics Observatory.
Figure 2 Variation of the total energy budget in the corona as a functional variation of helicity budget in AR 11560, AR 11726, and AR 11928. The dotted line having indicated slope and y-intercept in each panel refers to a straight line fit to the curve of variation of energy and helicity. The slope of this line indicates the constant torsion parameter α/8\π, and the y-intercept equals the potential energy (Ep).
Under the line-tied assumption, the field lines anchors from the photosphere. These field lines build stress by continuous slow plasma motions, which include vertical and horizontal motions. These motions are believed to be prime factors of complexity of magnetic field and its energy buildup.

Magnetic helicity is an important topological property of solar active regions (ARs) and is a measure of twist and writhe of the field lines (Berger & Field 1984; Finn & Antonsen 1985). To derive the time rate of helicity injection by these plasma motions, time sequence of photospheric magnetic field are being used. In a study of three emerging ARs (Vemareddy 2015), we found that the coronal helicity flux is dominantly coming from the shear term that is related to horizontal flux motions, whereas energy flux is dominantly contributed by the emergence term. Moreover, the ARs with same sign of Helicity-flux over time are found to launch CMEs but NOT in the one with helicity flux changing sign over time. And the non-linear curve of coronal energy vs helcity is approximately fitted with linear force-free model.
During the evolution of ARs, the sunspots exhibit proper as well as rotational motions, believed to be associated with the storage of energy by increasing the magnetic non-potentiality and eventual release in the form of observed eruptive events. We studied the rotating sunspots in AR 11158 and their probable role in the observed activity. Figure 2(a) shows a huge coronal mass ejection pointed by arrow on February 14, 2011. The rectangular region in top-left panel is shown in panel (top-right) in continuum intensity observed by Heliosiesmic Magnetic Imager with the labeled sunspots. We have noticed intrinsic rotation of sunspots SN1 and SP1 by scrutinizing animations of these images. The sunspot SN1 is remapped to polar coordinate system to have the information about these rotational parameters like angular distance (θ) and angular velocity (/dt) with respect to time. In bottom-right panel, we have plotted these parameters with respect to time, delineates that at maximum rotational speeds the CMEs occurred.
Figure 3 NOAA AR 11158 observed by (top-left) Atmospehric Imaging Assembly aboard Solar Dynamic Observatory in 304Å channel showing CME pointed along the arrow (top-right) in continuum intensity with labelled sunspots. SN1 and SP2 are observed to have intrinsic rotational motion, (bottom-left) movie showing rotation of sunspot SN1, (bottom-right) Plot of the feature displacement and angular velocity around the centroid with respect to time. Note that the observed CME occurred when SP1 has maximum rotational velocity marked by arrow.


We further estimated various physical parameters representing non-potential nature as the AR evolves. The temporal profiles of twist parameters, namely, average shear angle, aav, abest, derived from HMI vector magnetograms, and the rate of helicity injection, obtained from the horizontal flux motions of HMI line-of-sight magnetograms, correspond well with the rotational profile of the sunspot in the CME-prone region, giving predominant evidence of rotational motion causing magnetic non-potentiality. Moreover, the mean value of free energy from the virial theorem calculated at the photospheric level shows a clear step-down decrease at the onset time of the flares revealing unambiguous evidence of energy release intermittently that is stored by flux emergence and/or motions in pre-flare phases. Additionally, distribution of helicity injection is homogeneous in the CME-prone region while in the flare-prone region it is not and often changes sign. This study provides a clear picture that both proper and rotational motions of the observed fluxes played significant roles in enhancing the magnetic non-potentiality of the AR by injecting helicity, twisting the magnetic fields and thereby increasing the free energy, leading to favourable conditions for the observed transient activity. [Ref: Vemareddy et al, 2012, Astrophys. J, 761, 60]
Figure 4 Magnetic strcture around coronal null point in AR 11166. Fan separatrix surface intersects the photosphere along flare ribbons.
During the AR evolution, build-up of energy in the corona takes places through flux emergence and displacements for days (Schrijver 2009). Under these circumstances, the characteristic time of propagation of Alfven waves through the whole region is much less than the corresponding time of the global evolution of the AR, so that the dynamical effect is negligible. Moreover, the pressure and weight of plasma are also negligible in comparison with the magnetic pressure and tension, which therefore counterbalance each other to vanishing Lorentz force. Therefore, this slow change of AR evolution can be approximated by series of quasi-stationary, force-free equilibrium states in the low-plasma beta coronal environment.

Using HMI magnetic field measurements as bottom boundary conditions, we modelled the coronal magnetic field of AR 11166 around the time of X1.5 flare occurence. The magnetic structure reveals a null point topology with dome-like fan saperatrix surface and tail-like spine connecting remote ribbons (Vemareddy & Wiegelmann 2014). Fan-separatrix surface intersects the photosphere along the observed flare ribbons. It suggests that the reconnection triggers at the pre-existing null point, from where particles impinge along field lines to form chromospheric flare ribbons.
Electric currents in the astrophysical plasmas are generated by distortion of the magnetic field by external forces applied by field free plasma (Parker 1979). In the Sun, these currents enter the corona through embedded emerging magnetic fields from the convection zone to the corona through the photosphere. Despite their generation mechanisms, these currents are believed to play vital role in the magnetic solar eruptions. However, the Parker's theoretical prediction (Parker 1996) of neutralized net vertical current over a cross-section of flux tube at the photosphere is inconsistent to their role in eruptions. In many studied ARs (Ravindra et al 2012, Georgoulis et al 2012, Vemareddy et al 2015), however the observance of the non-neutralized net current is therefore needs proper explanation.
Figure 6 Evolution of |DC/RC| in north (Bz>0) and south (Bz<0) in AR 11166 and 12371. The ratio is also plotted for twist current ($I_ch$) in respective panels. The net magnetic flux is plotted in grey with y-axis scale on right of each panel. Horizontal dotted line marks current nuetralization level of unity.
Breakdown of net current neutralisation has implications to CME activity from the AR. The hoop-force of CME flux rope is proportional to square of the net current in the flux rope channel (Zakharov & Shafranov 1989), which is the algebraic sum of net direct current (DC) and return current (RC). The equilibrium and stability properties of the flux rope then requires the neutralisation of net current for which hoop force vanishes. Numerical simulations showed that a major breakdown occurs during emerging phase of AR (Torok et al 2014), during which opposite polarity regions are compact with a sheared PIL. Observational studies also found that the sheared PIL is the major contributor to DC and proposed the degree of current-neutralisation as the proxy for assessing the ability of ARs to produce CMEs (Liu et al 2017).

Using HMI photospheric vector magnetic field, we can only measure vertical component of electric current from Ampere's Law. In Figure 6, we plot the ratio of |DC/RC| (blue) in both polarities of flaring-only AR 11166 and CME-rich AR 12371. The evolution of this ratio is well above unity due to highly non-neutralized current in AR12371, whereas it is near unity indicating nearly neutralized currents in AR 11166. The red curve is the |DC/RC| for twist current, which follows above that of net current (blue). This is due to the fact first found in Vemareddy (2017) that the twist current is opposite in sign dominantly to that of shear current, which comes from gradient of field strength. From a study of 20 ARs of different activity, we conclude that the non-neutralised currents arise from sheared polarity inversion lines, whose segmented nature links to the flare/CME activity in the AR. (Ref: Vemareddy 2019, MNRAS)
Figure 5 Cylindrically symmetric constant linear force-free Lundquist model fitting (solid red) and non-linear force-free GH (dashed blue) model fitting are plotted against the in situ observations (black) in each panel. MC expansion is accounted in the models. Notice the better fitted rotating component with GH model.
Magnetic clouds (MCs) are large scale, organised magnetic structures in interplanetary space which are observed in situ as interplanetary coronal mass ejection (ICMEs).They are generally preceded by the occurrence of major coronal mass ejections (CMEs) at the Sun. It is now believed and shown from a variety of independent studies that MCs are magnetic flux ropes (MFRs) of locally straight cylindrical geometry (Burlaga 1988, Shodhan et al 2000, Liu et al 2008, Gopalaswamy et al 2013a). In this picture (see Figure 11), the MC is thought to be part of a large-scale bent flux rope extending from the Sun into interplanetary space with its feet possibly still connected to the Sun (Burlaga et al 1990, Bothmer and Schwenn 1998, Ferrugia et al 1993b, Schwenn 1998, Webb et al 2000). From the point of Sun-Earth connection, a major interesting, important question is how the solar source regions are connected to the in situ MCs, which should lead to important clues on how to forecast the internal magnetic field of CMEs around Earth and other planets.

A MC passage through the space-craft in in-situ is generally identified by a strong magnetic field strength, low plasma beta, and rotation of any magnetic field component (reversal of sign during passage). In the flux rope scenario, cylindrically symmetric force-free models fit the in-situ MC observations well (see Figure 12). Since the direction of axial component (southward/northward) in MC decides the occurence of geo-magnetic storms, it is crucial to determine the orientation of the flux rope and its axial field polarity.

In the above flux rope scenario, MC magnetic field observations are fitted (Figure 4) with Cylindrically symmetric constant linear force-free Lundquist model (solid red) and non-linear force-free Gold-Hoyle (dashed blue) models. A recent study (Vemareddy et al 2016) found that Gold-Hoyle models best fits the observations, in terms of chisquare, especially when MC expansion is accounded in the model.
Figure 6 Magnetic evolution in AR 12371 during its disk transit, (a) disk integrated soft X-ray GOES flux. Vertical lines mark the CME start times, (b) net magnetic flux, Φ, from north and south polarities, (c) coronal twist parameter av averaged over the AR, (d) helicity flux injection (in blue), and coronal helicity H normalized by Φ2 (in red), (e) energy flux injection dE/dt (in blue) and its coronal accumulation (in red). The observed decrease in net magnetic flux is due to flux cancellation at the polarity inversion line. Corresponding to the persistent strong shear motions, the H-flux is strong up to June 21 and accordingly |αav| shows an increasing trend.
Active regions (ARs) with sigmoidal structure, generally seen in soft X-rays and EUV, are prone to CME-eruptions and associated flares. Some ARs show rapid succession of CMEs and flares over a timescale of few hours. Because the timescale is too small compared to typical timescale for energy build-up in ARs, emergence of magnetic flux was suggested to be responsible for repeated CMEs of such short timescale. On the other hand, successive CMEs also occur from ARs in a timescale comparable to energy build-up by footpoint motions in the post emergence phase. The velocity field derived from tracked magnetograms indicates frequently persistent shear and converging motions of polarity regions about the polarity inversion line. These motions introduce huge amounts of magnetic energy and helicity in the AR magnetic system enabling flares and CMEs. While stored energy configuration is vital, the persistence of sigmoidal structure is the central feature for the successive CME occurrence from a source AR and had been the central subject for eruptive activity.

AR 12371 is in post-emerged phase producing successive fast CMEs in a span of its disk transit (18-25 June, 2015). EUV observations at the Sun register the initiation times of four eruptions as 15:05UT on 18 (CME1), 00:45UT on 21 (CME2), 16:15UT on 22 (CME3), and 07:30UT on 25 (CME4) of June 2015, respectively. Subsequently, the disk integrated GOES X-ray light curve (Figure 6a) delineates that the CMEs are associated with long duration M-class flares. The EUV double dimming, three-part CME structure, and fast propagation speeds in LASCO FOV (>1000kms-1) altogether characterize the CMEs as homologous events.

The magnetic evolution in this AR is summarised in Figure 6. The region of sigmoid is seen with persistent shearing and converging motions of opposite magnetic polarities over 7-days of evolution. At the same time, the net flux decreases over time (24%) with the build-up of magnetic twist (till 22nd June). The helicity flux, derived from footpoint motions, is negative and co-spatial with the coronal sigmoid. This implies qualitatively that the pumped helicity flux is being utilised in the repeated build-up of sheared core along the polarity inversion line. Besides injection helicity flux, the energy flux injection by flux motions is above 5x1027 ergs/s, which amounts to that required by a fast CME (~1033ergs) over about 50hours. Importantly, the normalised helicity flux is moderately high (0.15turns) compared to non-eruptive super AR 12192 (0.02 turns). A high value of magnetic flux normalized helicity flux typically suggests presence of a coherent twisted field.

Figure 7 First column: Pre-eruptive magnetic structure, for CME2, CME3, CME4 respectively, in AR 12371. Field lines are obtained from non-linear force-free field extrapolations. The background maps are the vertical field magnetic field component Bz (white: positive, black: negative). Major magnetic polarities are labelled with P*, N*. The white rectangle in 22T17:24UT panel marks the core region of the sigmoid shown in the right column. Second column: Magnetic field lines overplotted on AIA 94. Third column: Magnetic field lines in the AR core showing large scale sheared arcade and a weakly twisted core.
Force-free models of the AR are computed from the observed photospheric magnetic field. The computed magnetic structures reproduce the EUV sigmoid mostly with sheared arcades. Especially, the core-field is a combination of double inverse J-shaped and inverse S-shaped field lines with dips touching the photosphere. Such field-lines are formed by flux cancellation reconnection of opposite-J field-lines at bald-patch locations (where field-lines are tangent to the photosphere and curved upward (Figure 7). This study demonstrates the formation of a weakly twisted flux-rope from large scale sheared arcade field-lines. Further, our study infers a steep decrease of the background coronal field meeting the torus instability criteria at low height (~40Mm) in contrast to non-eruptive ARs. When the magnetic field becomes torus unstable, the weakly twisted flux rope is upward destabilised. This induces further reconnection below the flux rope and leads to a CME as the overlying magnetic field is weak. A new, but similar, stressed magnetic configuration is next reformed by photospheric motions leading to the repetition of the above scenario. This induces a series of homologous CMEs/flares. Finally, the modelled structure captured the major features sigmoid-to-arcade-to-sigmoid transformation, that is being cyclic under continuous photospheric flux motions. (also appeared as RHESSI sceince nugget)
Figure 8 Magnetic evolution in AR 12673: a) net flux in positive (north) and negative (south) magnetic polarities. Disk integrated GOES X-ray (1.0-8.0Å passband) flux is also shown with y-axis scale on the right. b) Evolution of net vertical current from north and south polarities. c) Neutralization of net current in individual polarities. Total length of all SPIL segments are also plotted with y-axis scale on right. d)time-rate of helicity flux normalized by averaged net flux of positive and negative polarities. Normalized accumulated helicity (H(t)/Φ2) is also plotted with y-scale on right. e) Energy flux injection and its accumulative quantity. Rapid flux emergence phase is marked with grey shade and the two major flaring phases are indicated with orange shade.
Magnetic helicity is believed to play a fundamental role in generating large-scale activity on the Sun. Magnetic helicity describes the magnetic field complexity, including twist, writhe, knot, and linkages of magnetic field. When coronal magnetic field is being pumped with helicity and energy, the magnetic complexity and non-potentiality increases. Using simultaneous magnetic field and coronal observations, several studies focused on the magnetic non-potential parameters to understand their roles in CMEs and flares. From those studies, it is now believed that the magnetic flux rope is built up by line-tied photospheric flux motions, such as magnetic flux emergence or horizontal flows. These processes inject magnetic helicity into the higher solar atmosphere, increasing the twist and kink of a flux rope (self-helicity) and the linkage between different flux ropes (mutual helicity). Magnetic helicity is conserved in an ideal MHD process and changes very slowly in a resistive process[1]. Thus, a flux rope with continuous injection of magnetic helicity inevitably erupts to remove the accumulated helicity, manifested as a CME.



Figure 9 Magnetic field extrapolation by force-free field approximation to the corona. a) Field lines showing the magnetic structure of AR 12673 on September 4, 2017 at 18:00UT. Background is radial magnetic field map. b) Magnetic structure on coronal 94Å image. c) Map of vertically integrated electric current. d) Coronal image of 304Å from AIA.
The above scenario of helicity role in the AR eruptivity is well portrayed in a recent study. Using the uninterrupted observations of space-borne Solar Dynamics Observatory, we studied the most violent AR 12673 that produced strongest flares in Solar Cycle 24. The AR emerged with two bipolar regions on September 2, 2017 in a pre-existing positive-polarity sunspot. The velocity field derived from tracked vector-magnetograms indicates persistent shear and converging motions of flux regions about the polarity inversion line (PIL). In Figure 8, the time evolution of the X-ray flux, net-electric current, helicity, energy flux injection are shown. A major helicity injection occurs during rapid flux emergence, consistent with the very fast flux emergence phase. While this helicity flux builds up the sigmoid by September 4, the helicity injection by the continued shear and converging motions in the later evolution contributes to the sigmoid sustenance and its core field twist, as a manifestation of the flux rope that erupts after exceeding critical value of twist. Moreover, the total length of the sheared PIL segments correlates with the non-neutralized currents and maintains a higher value in both polarity regions as a signature of eruptive capability of the AR, according to flux rope models. This result also implies that the non-neutralized currents arise in the vicinity of the sheared PIL. In Figure 9, an example of force-free extrapolation to the AR 12673 is displayed. The modeled magnetic structure well resembles the coronal observation of AIA 94Å, capturing major features like twisted core flux as flux rope, and hook-shaped sections connecting in the middle of the PIL. The emission in 304Å is also reproduced in vertical integration of volume electric current. The study of quasi-separatrix-layers reveals that the sheared arcade, enclosing the flux rope, is stressed to a critically stable state and its coronal height becomes doubled from September 4-6.

Figure 10 Comparison of normalized helicity flux injection (dH/dt/Φ2) in 4 flare/CME producing ARs. Starting times are 2017-08-31T19:12 UT, 2015-06-19T00:00 UT, 2014-10-21T00:00 UT, 2012-03-06T00:00 UT , for AR 12673, 12371, 12192, 11429 respectively. Normalized helicity flux injection in AR12673 is comparatively high with major injection phase co-temporal with rapid flux emergence.
In order to have more insight on the helicity flux input, we compare normalized helicity flux injection (dH/dt/Φ2) from different ARs (Figure 10). This value is a measure of non-potentiality per unit flux per unit time and indicates how fast the AR accumulates energy and helicity in the corona. As is clear from the plot, this parameter evolves at a higher rate by a factor of 3 in AR 12673 compared to other ARs. The AR 12192 is a flare-prolific region without CMEs in contrast to the rest of the ARs[4], and has small injection value per flux tube. Interestingly, this value in the AR 12673 is quite stronger by a factor of 2 than the strong CME-prolific ARs 12371 and AR11429. This suggests that the normalized helicity flux is a key parameter to distinguish eruptive and non-eruptive ARs and also predicts the severe space-weather events. (also appeared as HMI sceince nugget)