New Evidence that Magnetoconvection Drives Solar–Stellar Coronal Heating

, , , , and

Published 2017 July 5 © 2017. The American Astronomical Society. All rights reserved.
, , Citation Sanjiv K. Tiwari et al 2017 ApJL 843 L20 DOI 10.3847/2041-8213/aa794c

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

2041-8205/843/2/L20

Abstract

How magnetic energy is injected and released in the solar corona, keeping it heated to several million degrees, remains elusive. Coronal heating generally increases with increasing magnetic field strength. From a comparison of a nonlinear force-free model of the three-dimensional active region coronal field to observed extreme-ultraviolet loops, we find that (1) umbra-to-umbra coronal loops, despite being rooted in the strongest magnetic flux, are invisible, and (2) the brightest loops have one foot in an umbra or penumbra and the other foot in another sunspot's penumbra or in unipolar or mixed-polarity plage. The invisibility of umbra-to-umbra loops is new evidence that magnetoconvection drives solar-stellar coronal heating: evidently, the strong umbral field at both ends quenches the magnetoconvection and hence the heating. Broadly, our results indicate that depending on the field strength in both feet, the photospheric feet of a coronal loop on any convective star can either engender or quench coronal heating in the loop's body.

Export citation and abstract BibTeX RIS

1. Introduction

Active regions have the strongest magnetic field on the Sun and contain the brightest extreme-ultraviolet (EUV) and X-ray coronal loops (Withbroe & Noyes 1977; Vaiana & Rosner 1978; Golub et al. 1980; Reale 2014), heated to 2–6 MK (Zirker 1993; Reale 2014) by unknown magnetic processes (Fisher et al. 1998; Schrijver et al. 1998; Moore et al. 1999; Katsukawa & Tsuneta 2005; Klimchuk 2006; De Pontieu et al. 2007). As a rule, the stronger the photospheric magnetic flux, the brighter the corona—the EUV/X-ray corona in active regions is 10–100 times more luminous and 2–4 times hotter than that in quiet regions and coronal holes, which are heated to only about 1.5 MK (Withbroe & Noyes 1977; Vaiana & Rosner 1978; Zirker 1993; Reale 2014) and have fields that are 10–100 times weaker than active regions (Wiegelmann et al. 2014). The two most widely argued mechanisms for coronal heating are (i) magnetic-wave heating (e.g., De Pontieu et al. 2007; van Ballegooijen et al. 2011), which is possibly dominant for heating the coronal plasma in the quiet-Sun and coronal holes, and (ii) nanoflare heating (Parker 1983, 1988), which is possibly dominant in active region coronal loops (Cirtain et al. 2013). For (ii), the continuous shuffling of the field in the feet of coronal loops by the (sub-)photospheric convection supposedly braids the field in the coronal loops. This results in the corona being heated to EUV and X-ray temperatures by Ohmic dissipation of magnetic energy via reconnection in these braided loops, impulsively at many small-scale current sheets, at a rate of 107 erg cm−2 s−1 (Parker 1983, 1988). Recent observations and coronal magnetic field modeling have provided evidence for the presence of braided loop structures and nanoflaring in the active region corona (Cirtain et al. 2013; Winebarger et al. 2013; Brosius et al. 2014; Thalmann et al. 2014; Tiwari et al. 2014).

Active regions with stronger magnetic fields (often having sunspots) have brighter coronal emission than active regions with weaker magnetic fields (Withbroe & Noyes 1977; Vaiana & Rosner 1978; Golub et al. 1980; Reale 2014). Most coronal loop feet in sunspots are rooted in the umbra or in penumbral spines in the inner penumbra (these penumbral spines are intrusions of steeply inclined strong umbral field into the inner penumbra; Tiwari et al. 2015). Sunspots host one foot of some of an active region's brightest coronal loops observed in EUV and X-ray wavelengths, and some of these have that foot in umbra (e.g., Foukal 1975; Alissandrakis & Patsourakos 2013), consistent with sunspot umbrae having the strongest field in an active region. Even so, umbrae often have the darkest areas of an active region in coronal EUV and X-ray images (Pallavicini et al. 1979; Webb & Zirin 1981; Sams et al. 1992; Golub et al. 1994). Here, we resolve this paradox (of why some active region coronal loops rooted in an umbra are the brightest but others are the dimmest) by showing that the coronal EUV brightness of active region coronal loops depends on the convective freedom in both photospheric feet of the loops, and that considering only the loop ends rooted in an umbra and not knowing the rooting of the opposite ends, as was the case in previous studies (e.g., Foukal 1975; Pallavicini et al. 1979; Webb & Zirin 1981; Sams et al. 1992; Golub et al. 1994; Alissandrakis & Patsourakos 2013), resulted in this paradox.

2. Data and Modeling

2.1. Instrumentation and Data

We use Solar Dynamics Observatory (SDO) Atmospheric Imaging Assembly (Lemen et al. 2012; AIA) images, mainly from the 94 and 193 Å channels, for identifying and tracking active region EUV coronal loops. The 193 Å channel of AIA detects Fe xii emission at about 1.5 MK. The AIA 94 Å channel is mainly sensitive to hot emission centered on an Fe xviii line (6–8 MK), but it also detects emission from about 1 MK plasma (Martínez-Sykora et al. 2011; Testa & Reale 2012; Warren et al. 2012; Boerner et al. 2014). We also inspected our active regions in the images from all other AIA channels.

AIA UV 1600 Å channel images are used for identifying the locations of plage in coronal images. The AIA 1600 Å passband primarily passes lower-chromospheric continuum emission (characteristic temperature of 5000 K) but also transmits radiation from two CIV lines (at ∼1550 Å) formed at 0.1 MK in the transition region. Six AIA channels are used for differential-emission-measure (DEM) analysis. The pixel size of AIA images is 0farcs6. The time cadence for EUV AIA channels is 12 s and that for UV channels is 24 s. However, we have used three-minute-cadence AIA movies for the present study. We removed solar rotation by derotating all the SDO images to a particular time.

Line-of-sight magnetograms from the SDO/Helioseismic and Magnetic Imager (HMI; Schou et al. 2012) are used to map the magnetic field strength in active regions. HMI provides line-of-sight magnetograms every 45 s, with a pixel size of 0farcs5. The cadence of the line-of-sight magnetograms used in the movies of this paper is also three minutes. For nonlinear force-free (NLFF) field modeling, which is described next, SDO/HMI vector magnetograms at several selected times are used. HMI vector magnetograms of active regions are computed using the Very-Fast-Inversion of the Stokes-Vector algorithm (Borrero et al. 2011) and are available at a 12 m cadence.

2.2. Nonlinear Force-free Field Modeling

To date, there are no routine direct measurements of the coronal magnetic field. To compensate for this lack, NLFF magnetic field modeling based on photospheric vector magnetograms is presently the most accurate way to deduce the three-dimensional coronal magnetic field (De Rosa et al. 2009; Wiegelmann & Sakurai 2012). We first transform the HMI vector magnetograms to the solar disk center (Gary & Hagyard 1990), i.e., we deduce the vertical and horizontal vector components of the magnetic field from the image-plane measurements (i.e., from the line-of-sight and transverse field). Because the photospheric magnetic field of active regions is not completely force-free (Metcalf et al. 1995; Gary 2001; Moon et al. 2002; Tiwari 2012), a preprocessing technique is applied to the HMI vector magnetograms (Wiegelmann et al. 2006) to achieve suitable force-free boundary conditions (e.g., Low 1985) for the coronal NLFF field modeling. The optimization procedure is described in earlier publications (e.g., Wiegelmann et al. 2012).

Since the coronal plasma is frozen to the magnetic field, the bright coronal loops trace the coronal magnetic field. We find a great number of projected model magnetic field lines that trace the coronal loops observed in AIA 193 and 94 Å images of our active regions. This verifies that our NLFF model field sufficiently accurately approximates the true coronal magnetic field. For a coronal loop visible in the AIA images, the matching model magnetic field line gives the magnetic setting of the coronal loop's photospheric feet. We performed modeling for each active region every two hours during the observation period.

It is important to note that from the NLFF-modeled field we need to specify only the general location of the loop feet, i.e., whether a loop foot is in unipolar plage, in mixed-polarity plage, in penumbra, or in umbra. This method is sufficiently accurate for this purpose (De Rosa et al. 2009; Wiegelmann & Sakurai 2012; Wiegelmann et al. 2012).

3. Results

In Figure 1, we show an example active region (NOAA 12108) in UV, EUV, continuum-intensity images, and a line-of-sight magnetogram. This active region has a simple bipolar field configuration and contains sunspots of opposite magnetic polarity. The active region was flare-quiet—it produced no C-class or larger flares within 24 hr before and after our observations. We selected a non-flaring active region because we are interested in quasi-steady heating and not in the large bursts of heating by flares. During the 26 hr of our observations, the active region passed the central meridian 8° south of disk center, and thus was always close to disk center, thus avoiding projection errors. Coronal loops of different temperatures can be identified in the coronal EUV images of SDO/AIA. Noticeably, the sunspots of the active region appear dark in continuum-intensity (Figure 1(b)) and UV images (Figure 1(a)). The sunspots are surrounded by weaker-field plage regions, which are bright in the AIA 1600 Å image (Figure 1(a)). The loops are rooted in different photospheric magnetic features such as sunspots, unipolar plage, and mixed-polarity plage (compare Figures 1(c)–(e)). The temperature map (Figure 1(f)), derived from taking into account the emission across the temperature regime ∼0.6–6 MK (corresponding to AIA 171–94 Å channels) through DEM, following Aschwanden et al. (2013), shows the approximate temperatures of coronal loops. The temperature map shows that the brightest coronal loops in the 94 Å image are the hottest ones. We emphasize that, although the DEM analysis indicates that the brightest AIA 94 Å loops are the hottest, our main results of the paper do not depend on the DEM calculations.

Figure 1. NOAA active region 12108 as observed near central meridian by SDO on 2014 July 07 at 23:45UT. (a) Photospheric and transition region UV emission (AIA 1600 Å). (b) Photospheric visible-light emission (HMI continuum). (c) Coronal EUV emission (AIA 193 Å). (d) Photospheric line-of-sight magnetic field (HMI:6173 Å). (e) Coronal EUV emission from the hottest plasma (AIA 94 Å). The white arrow indicates the brightest coronal loop at this time (see also Movie 1, in which arrows indicate five typical examples of the brightest loops in the 94 Å movie, at 02:36:01, 06:33:01, 12:36:01, 19:00:01, and 23:42:01 UT, and five typical examples of less-bright loops in the 193 Å movie, at 00:09:06, 03:30:06, 13:30:06, 14:15:06, and 00:48:06 UT). All images in Movie 1 are derotated to 22:00:06 UT on 2014 April 01. (f) A temperature map of the active region obtained by DEM analysis. Blue/red (blue/white in f) contours, which outline sunspot umbrae, are contours of ±1000 G from the smoothed (by 10 pixels) line-of-sight magnetogram.

(An animation of this figure is available.)

Video Standard image High-resolution image

We use sequences of AIA 193 and 94 Å images, such as those shown in Figure 1, to study the evolution of bright active region coronal loops (Movie 1). In these sequences (i) bright EUV coronal loops are episodic, with lifetimes of 1–2 hr; (ii) many of the brightest EUV loops are apparently rooted within or at the edge of a sunspot umbra at one foot, with their other foot located in opposite-polarity plage or penumbra; (iii) some of the other brightest loops have one foot rooted in umbra or penumbra, and the other foot in mixed-polarity flux; (iv) the bright coronal loops connecting opposite unipolar plage regions are never as bright as those with one foot rooted in a sunspot or in mixed-polarity plage; (v) coronal loops with each foot in a sunspot umbra of opposite polarity ("umbra-to-umbra" loops) are never visible, neither in the wavelength bands presented here, nor in any of the other AIA passbands.

To examine whether umbra-to-umbra coronal magnetic loops are present, we use a NLFF field model of the active region coronal magnetic field, based on a SDO/HMI photospheric vector magnetogram. In Figure 2(a), selected model field lines are shown overlaid on the nearest-in-time AIA 94 Å image. Each of these model field lines traces an observed loop closely, validating the model field. Despite the absence of any discernible connection between the two sunspot umbrae in EUV coronal images, the three-dimensional model magnetic field shows that there are field lines connecting the opposite-polarity umbrae of the active region (Movie 2). In Figures 2(b)–(c), we plot two sets of the model field lines that link the two opposite-polarity umbrae. Since the length of a coronal loop is an important factor in determining the loop's coronal brightness—shorter loops require a smaller amount of heating than longer ones for the same brightness (Mandrini et al. 2000; Klimchuk 2006; Kano et al. 2014; Reale 2014)—we display umbra-to-umbra loops of different apex heights (and therefore of different lengths). The two sets of umbra-to-umbra loops displayed in Figures 2(b) and (c) have apex heights of ∼90–100 mega-meters (hereafter Mm), "low-apex" fields, and ∼150–175 Mm, "high-apex" fields, respectively. Neither low-apex (short) nor high-apex (long) loops that connect the sunspot umbrae are visible in the AIA 94 and 193 Å images. We also verified the invisibility of these loops in the movies from all other AIA channels. This indicates that the coronal loops connecting opposite-polarity umbrae have the lowest coronal-temperature plasma density (and thus the lowest coronal brightness), and hence are the least heated (Vaiana & Rosner 1978). It is worth mentioning here that a dim umbra-to-umbra loop seen in AIA 171 Å images of the active region (NOAA 12108) two days after our observations (Chitta et al. 2016), is reasonable because a light bridge, which is a signature of magnetoconvection in the umbra, was present then in the leading-polarity sunspot, and was at the foot of that loop.

Figure 2. Model coronal magnetic field of NOAA active region 12108 overplotted on an AIA 94 Å image. The NLFF field model is based on an HMI photospheric vector magnetogram recorded at 23:36:00 UT. The background AIA 94 Å image is the nearest in time, at 23:36:01 UT. The blue/red contours here and in Movie 1 were made the same way as those in Figure 1. The height at each point along each model field line is color-coded. (a) Sample model magnetic field lines matching the coronal loops in 94 Å and/or 193 Å images. (b) A set of lower-altitude model umbra-to-umbra field lines having apex heights of 90–100 Mm. (c) A set of higher-altitude model umbra-to-umbra field lines having heights of 150–175 Mm. The higher/longer umbral field lines are naturally rooted near the outer edge of each umbra, in or close to penumbral spines in the inner penumbra. None of the umbra-to-umbra model field lines match any discernible coronal loops. See Movie 2 for a three-dimensional view of these umbra-to-umbra model field lines.

(An animation of this figure is available.)

Video Standard image High-resolution image

In Figure 3, we show images and the temperature map of another flare-quiet active region (NOAA 12021), the model field lines of which are shown in Figure 4. We find the same pattern for the EUV brightness of coronal loops with respect to their magnetic rooting within the active region. Again, the modeled umbra-to-umbra field lines do not match any discernible EUV coronal loops within an extended 24 hr time span (Movie 3). Although here we show ∼24–26 hr of data for only two examples, we have used JHelioviewer Software to find  similar active regions characterized by a large sunspot for each polarity with no light bridges in the umbrae. These active regions show no exception to our results, e.g., we do not find any visible umbra-to-umbra coronal loops.

Figure 3. Similar to Figure 1 but for a second example, NOAA active region 12021, observed on 2014 April 01/02, when it was about as near disk center as the first active region was during its observations. The contours have the same levels and meaning as in Figure 1. An arrow in e again points to the brightest and hottest loop at this time (see also Movie 3, in which arrows indicate five typical examples of the brightest loops in the 94 Å movie, at 12:30:01, 13:18:01, 20:57:01, 21:21:01, and 04:39:01 UT, and five typical examples of less-bright loops in the 193 Å movie, at 10:00:06, 01:12:06, 02:30:06, 06:21:06, and 09:54:06 UT). All images in the Movie 3 are derotated to 11:30:40 UT on 2014 July 07.

(An animation of this figure is available.)

Video Standard image High-resolution image

Figure 4. Model active region coronal magnetic field of NOAA active region 12021. The NLFF model field is based on the HMI vector data observed on 2014 April 02 at 03:35:25UT, and the AIA 94 Å background image is the nearest in time at 03:36:01 UT. See Movie 4 for a three-dimensional view of these umbra-to-umbra model field lines.

(An animation of this figure is available.)

Video Standard image High-resolution image

4. Discussion and Conclusions

We interpret our observational and modeling results as follows. The convective freedom at both photospheric feet is the primary determinant of the brightness of active region coronal loops. Umbra-to-umbra loops are the dimmest (invisible) because the convection in each foot of these loops is so strongly suppressed by the strong umbral fields that the magnetic energy input to those loops is quenched. That is, field braiding for powering heating by magnetic reconnection is quenched in these coronal loops. On the one hand, the EUV invisibility of umbra-to-umbra loops and the observed magnetic rooting of the bright EUV loops are consistent with Parker's idea (Parker 1983) that convection shuffles field lines and the subsequent small-scale reconnection (nanoflares; Parker 1988) in the resulting braided coronal loops (Cirtain et al. 2013; Winebarger et al. 2013; Brosius et al. 2014; Tiwari et al. 2014; Thalmann et al. 2014; Pontin et al. 2017) releases magnetic energy to heat the corona, which implies that the strong suppression of convection by the strong field in umbrae should quench coronal heating in umbra-to-umbra loops. On the other hand, any appreciable heating from magnetic waves generated by convection-churning of the loop feet is evidently also suppressed in umbra-to-umbra coronal loops.

Because the brightest coronal loops are the ones linking an umbra or inner penumbra to unipolar plage, mixed-polarity plage, or penumbra, we infer that they are the brightest because they have one foot rooted in umbra or inner penumbral spines and therefore have stronger magnetic field in the loop body than do plage-to-plage loops of the same length (which never become as bright as the brightest loops). Plausibly, because of that and because the other foot has vigorous convection, these loops with only one end in umbra have a higher rate of magnetic energy injection (higher Poynting flux), and a higher dissipation rate via reconnection, and thus brighter coronal emission than the plage-to-plage loops. All previous studies of the brightness of coronal loops stemming from an umbra considered only the loop feet in that umbra and did not know the magnetic settings of the other ends of the loops. Thus, they did not discern why some active region coronal loops rooted in umbra are the brightest (Foukal 1975; Alissandrakis & Patsourakos 2013) and others are the dimmest (Pallavicini et al. 1979; Webb & Zirin 1981; Sams et al. 1992; Golub et al. 1994). The sketch in Figure 5 depicts what we infer to be the dependence of the brightness of coronal loops on their magnetic rooting in an active region. Our observations imply that as long as the field can be braided by convection in a loop foot, the stronger the field in the loop, the stronger the coronal heating. Our results qualitatively support the models of Hurlburt et al. (2002) and Chen et al. (2014).

Figure 5.

Figure 5. A schematic drawing depicting the dependence of the coronal EUV brightness of active region coronal loops on the photospheric magnetic rooting of the loops. Yellow/blue  represent positive/negative polarity of surface magnetic flux. "U" and "P" stand for sunspot umbra and penumbra, respectively. The different loop colors indicate different brightnesses, with red/dark-gray being the brightest/dimmest, and bright orange being dimmer than red and brighter than the pale orange. Each of the two foreground positive and negative magnetic areas is a unipolar plage region. A mixed-polarity plage region is present on the left in the background, where one of the brightest loops has one foot, the other foot being in a sunspot penumbra of positive polarity. The taller of the two penumbra-to-penumbra coronal loops is the dimmer, presumably because it is longer.

Standard image High-resolution image

Thus, magnetoconvective driving at both feet of a loop is the key. Whether the subsequent transport and dissipation mechanism is braiding/nanoflare or wave heating is beyond this investigation. Irrespective of whether the active-region corona is heated more by dissipation of magnetic waves or more by nanoflares, our results indicate that, depending on the field strength in both feet, the photospheric feet of any coronal loop on the Sun—or on any star with a convective envelope—can either engender or quench coronal heating in the loop's body. Our results are new evidence for the widely held view that the quasi-steady coronal heating in active regions is driven by magnetoconvection in the feet of the loops (Parker 1983, 1988). The alternative possibility, that the heating mostly comes from dissipation of free magnetic energy already in the loop field when it emerged from inside the Sun (Fisher et al. 1998), is unlikely in view of the EUV invisibility of umbra-to-umbra loops.

We observed, as have a few other researchers (Falconer et al. 1997; Chitta et al. 2017), that the active region coronal loops that are most persistently the brightest are rooted in regions of mixed-polarity magnetic flux at one end or at both ends. This indicates that the evolutionary interaction of the opposite-polarity flux, which is driven by magnetoconvection in and below the photosphere, somehow enhances the rate of injection of free magnetic energy into the accessed coronal loops, in excess of what the rate would be were the magnetic flux unipolar at both ends.

We thank Allen Gary, Alphonse Sterling, David McKenzie, Davina Innes, and P. Venkatakrishnan for useful comments on the manuscript. S.K.T., R.L.M., and A.R.W. were supported by funding from the LWS-TRT program and the HGI Program of the Heliophysics Division of NASA's SMD. N.K.P.'s research was supported by an appointment to the NASA Postdoctoral Program at the NASA/MSFC, administered by USRA under a contract with NASA. J.K.T. acknowledges support from the Austrian Science Fund (FWF): P25383-N27. We acknowledge using NASA/SDO-AIA-HMI data and the software JHelioviewer.

Please wait… references are loading.
10.3847/2041-8213/aa794c