Free Access
Issue
A&A
Volume 647, March 2021
Article Number A188
Number of page(s) 13
Section The Sun and the Heliosphere
DOI https://doi.org/10.1051/0004-6361/202040111
Published online 02 April 2021

© ESO 2021

1. Introduction

The evolution of the photospheric magnetic field plays a key role in the energy transport into the chromosphere and the corona. In active regions, the emergence of new magnetic flux leads to interactions with the pre-existing magnetic field, releasing magnetic energy, heating the solar atmosphere, and accelerating the solar plasma. Signatures of such activity are commonly present in observations: from compact brightenings such as Ellerman bombs (EBs), ultraviolet bursts (UVBs), and microflares (MFs) to big eruptive flares or X-ray jets (Archontis 2012; Priest 2014; Young et al. 2018). These phenomena are thought to be manifestations of magnetic reconnection that occurred at different plasma conditions. Waves also have an important role in the transfer of energy as those jets develop into shock waves as they propagate into the corona (Yokoyama & Shibata 1995).

A particular set of these small-scale reconnection events are induced by the interaction of opposite-polarity magnetic elements, which can leave observational imprints in the photosphere and the chromosphere. Among these phenomena, EBs (Ellerman 1917) are observed as rapid flames visible in the wings of Hα or Ca II 8542 Å tracing magnetic reconnection in the lower photosphere (Vissers et al. 2015). Alternatively, UVBs (Peter et al. 2014) are compact transient brightenings observed in UV bands, with very complex Si IV profiles, and do not generally show a co-spatial brightening in the coronal channels of the Solar Dynamics Observatory (SDO) Atmospheric Imaging Assembly (AIA) instrument. Similar targets have been studied with recent simulations (Danilovic 2017; Hansteen et al. 2017, 2019; Libbrecht et al. 2020) and observations (Socas-Navarro et al. 2006; Toriumi et al. 2017; Libbrecht et al. 2017; Guglielmino et al. 2018; Smitha et al. 2018; Vissers et al. 2019; Yadav et al. 2019; Ortiz et al. 2020; da Silva Santos et al. 2020). Microflares or stronger events show a visible response in coronal channels and they would originate from or extend into higher layers. Bi-directional flows at these locations are commonly reported in both observations and simulations (Peter et al. 2014; Vissers et al. 2019; Hansteen et al. 2017). They usually have timescales of a few minutes (Vissers et al. 2019; Smitha et al. 2018), but can live up to several hours (Guglielmino et al. 2018; Chitta et al. 2017). Therefore, differences in their lifetime, size, location, and so on, are used to distinguish them. This classification is even more difficult when several of these phenomena can occur co-temporarily which is due to simultaneous reconnection at different heights (Ortiz et al. 2020; Hansteen et al. 2019). They are sometimes associated with a cool counterpart known as surges, which are observed as darkenings in images taken in the Hα blue and red wings with large velocities (Yokoyama & Shibata 1995; Nóbrega-Siverio et al. 2016).

Given their transient small-scale nature, high-spatiotemporal-resolution spectropolarimetric observations are also needed to discern the mechanisms involved in the reconnection. Observations and simulations of similar explosive events (e.g., Innes et al. 2015; Ni et al. 2015; Rouppe van der Voort et al. 2017; Nóbrega-Siverio et al. 2016; Peter et al. 2019) have revealed very small-scale plasmoids that are ejected at high speeds after the fragmentation of the current sheet (tearing instability) during magnetic reconnection that enables fast energy release. In contrast, other studies propose the Kelvin-Helmholtz instability to generate a turbulent velocity field which finally triggers fast reconnection (e.g., Lazarian & Vishniac 1999; Kowal et al. 2009; Jeffrey et al. 2018; Chitta & Lazarian 2020).

Chromospheric energy deposition can be estimated through the radiative losses, which account for the energy that must be replenished at any given time by heating mechanisms. Leenaarts et al. (2018) investigated the chromospheric heating in a strong flux-emerging region. They used the Ca II K brightness as a proxy of the radiative losses in the chromosphere, finding a good correlation with the total linear polarization, which originated from large-scale loop-like structures. However, they did not attempt to estimate radiative losses from their observations. To our knowledge, spatially resolved estimates of chromospheric radiative losses have not been carried out in emerging flux regions.

In this study, we present an analysis of chromospheric heating induced by magnetic flux cancelation rooted in the photosphere. Our target is larger in size than typical EBs or UVBs and more similar to that observed by Smitha et al. (2018) or Guglielmino et al. (2018). We reconstructed 3D models from non-local thermodynamic equilibrium (NLTE) inversions of high-resolution spectropolarimetric data acquired in the Ca II K, Ca II 8542 Å, and Fe I 6301 and 6302 Å lines. Our study focuses on the overall atmospheric stratification of temperature, line-of-sight (LOS) velocity, and magnetic field vector. In addition, we estimate for the first time the total radiative losses of such a strong magnetic reconnection event, allowing us to calculate the energy needed to sustain the high inferred temperatures and compare this with the magnetic energy available. We also present an analysis of the evolution of the polarization signals, as well as the detection of plasmoids during the time series.

2. Observations

2.1. Target and observational setup

Our target develops in the active region NOAA 12593, which started emerging on 2016-09-18 around 04:00 UT. We observed it on 2016-09-21 recording three time-series of around 11, 17, and 21 min at 11:00 h, 12:15 h, and 12:38 h UT, respectively. The observations were performed using the CRisp Imaging SpectroPolarimeter (CRISP; Scharmer et al. 2008) and the CHROMospheric Imaging Spectrometer (CHROMIS; Scharmer 2017) instruments at the Swedish 1-m Solar Telescope (SST; Scharmer et al. 2003). The active region was located close to solar disk center, with solar coordinates at the beginning of the observations (X, Y) = (400″, 25″) and an observing angle μ = 0.91.

The Fe I 6301.5 Å and the 6032.5 Å data were recorded with the CRISP instrument at 16 wavelength positions. The Ca II 8542 Å data were also recorded with CRISP at 21 line positions: two wavelength positions at ±1.7 Å from the line core, and 19 positions spaced evenly between ±0.765 Å. Both datasets were recorded in full-Stokes mode. The Hα line was sampled in 16 wavelength positions nonequidistantly spaced between −2.2 Å and +1.3 Å from the line core. Datasets in the He I 5875 Å (D3) line were also acquired, but they are not used in this study. A complete scan of the Fe I, Ca II, and H I lines required 37 s. In the first and third scan, He I D3 was acquired instead of the Fe I lines, with a cadence of 43 s. The Ca II K dataset was recorded using the CHROMIS instrument. The latter was sampled using two wavelength positions at ±1.33 Å from the line core, 19 positions spaced evenly between ±0.7 Å around the line center (i.e., 3933.682 Å), and a single wavelength point in the pseudo-continuum at 4001 Å. The cadence of the CHROMIS data is approximately 9 s. Unfortunately, this event has not been observed by RHESSI or IRIS.

2.2. Data reduction

Both data sets were reduced using the CRISPRED and CHROMISRED pipelines as described by de la Cruz Rodríguez et al. (2015) and Löfdahl et al. (2018) using the Multi-Object-Multi-Frame-Blind-Deconvolution image restoration technique (MOMFBD; van Noort et al. 2005) to remove the atmospheric distortion. Residual rubber-sheet seeing distortions were compensated as proposed by Henriques (2012).

The CRISP and CHROMIS datasets were accurately aligned by scaling up the CRISP pixel scale to that of CHROMIS (from 0.059″ to 0.038″). Given that the temporal cadence of the CHROMIS data is higher, CRISP images were interpolated to the CHROMIS cadence using nearest-neighbor interpolation in the temporal dimension. As a final step, we performed an intensity calibration of the data by scaling the intensity in a quiet-sun area in the observations to the values of a reference atlas taking into account the observed heliocentric angle following Neckel & Labs (1984) and Cox (2000)1.

After the standard data reduction, a final inspection showed some residual artifacts that were removed in the following way: polarized interference fringes of short period in the monochromatic images were isolated in the Fourier space and filtered out. Long-period fringes identified at nearly continuum wavelengths in the polarization signals were subtracted to the rest of the wavelength points by removing their contribution (scalar projection). For a more detailed explanation of those procedures, see Díaz Baso et al. (2019b). To reduce the noise in the polarization signals, we used the neural network presented in Díaz Baso et al. (2019a) on Stokes Q, U, and V. We estimate that the final data has a polarimetric accuracy of 6 × 10−4Ic.

2.3. Context data of Solar Dynamics Observatory

We used data provided by the Atmospheric Imaging Assembly (AIA; Lemen et al. 2012) and the Helioseismic and Magnetic Imager (HMI; Scherrer et al. 2012) to put the observations into context. HMI and AIA are both onboard the NASA Solar Dynamics Observatory (SDO; Pesnell et al. 2012). The active region at the time of our observations (around 12:15 UT) is displayed in Fig. 1. The field-of-view (FOV) of the SST (drawn as a black square) was surrounded by two small sunspots of different polarities (see the magnetogram at the bottom of Fig. 1). Inside the SST FOV, an emerging flux region was observed at that time near the location of some pre-existent pores, indicated with a green square; this is our region of interest (ROI).

thumbnail Fig. 1.

SDO/HMI continuum image and SDO/HMI LOS magnetogram of the active region NOAA 12593 at 12:15 h UT. The black square contour outlines the FOV covered by the SST and the smaller green square indicates the ROI.

We show the history of our ROI as seen in SDO observations in Fig. 2. The simultaneous continuum intensity and magnetograms show the evolution of the corresponding phenomena in the photosphere. Before 06:00 UT there was only a pore of negative polarity at the center of the ROI. A new magnetic concentration of positive polarity of around 5″ emerges at the surface around 06:00 UT (visible in the magnetogram of the same figure). After that, they converged towards each other and became smaller until the new magnetic concentration almost disappears after 13:00 UT (rightmost column). The total magnetic flux of both features decreases over time by a roughly similar amount. We estimated the total canceled flux to be around 4 × 1019 Mx, and the cancelation rate to be around 106 Mx s−1. The AIA panels (lower rows of Fig. 2) show the evolution of this event in the chromosphere and corona. After the magnetic concentrations are in contact, some brightenings are visible at the location of the flux cancelation (around 08:00 UT) that became brighter by one to two orders of magnitude at the end of the process around 11:00–13:00 UT (coinciding in time with our observations highlighted with a purple dashed square).

thumbnail Fig. 2.

Time evolution of the ROI (green region shown in Fig. 1) as seen from SDO-HMI and SDO-AIA channels. The purple dashed-line square indicates the temporal window where we have SST observations.

Some form of magnetic reconnection is expected in this situation of a new flux-emerging region inside a pre-existing active region. The energy released during this interaction may be responsible for the observed chromospheric and coronal brightening.

2.4. Swedish 1-m Solar Telescope observations

Although the brightenings are visible in the HMI and AIA channels, due to their moderate resolution (∼1.2″ in HMI and ∼1.6″ in AIA) it is not possible to study the internal structure of the region. In the much higher spatial resolution SST observations (close to 0.1″ at 4000 Å), the region shows a very complex morphology. The SST FOV is indicated in Fig. 3, where we highlight in green the same ROI defined in previous figures. The photospheric images at SST resolution contain elongated granules generated during the flux emergence as shown in the upper-left panel of Fig. 3. We also identify the complex shape of the different magnetic polarities during the interaction, as shown in the Stokes V at the wing of the Fe I line in the upper-right panel of the same figure.

thumbnail Fig. 3.

Overview images which overlap between the two instruments at four different wavelengths taken at 12:15 UT. The red arrows point to solar north and west, while the white arrows indicate the direction to the closest limb. The green square indicates our ROI which we analyzed in detail with the inversion code. The black dashed square is a close-up analyzed in Sect. 4.3.

In the chromosphere, the brightenings start to be visible in the wings and core of Ca II 8542 Å, Ca II K, and Hα as bright fibrils directed to the east, apparently emerging from the interface between both polarities with widths of 80–200 km (bottom panels of Fig. 3). These values are estimated using the Ca II K images by measuring the full-width at half-maximum (FWHM) of the intensity profile across the fibrils. Above those brightenings, we observe dark jets that are apparently escaping from the tail of the brightenings but they bend backward, moving towards the opposite direction (west). They are better visible in the blue wing of Hα with a projected length of approximately 10 Mm. Their detection in the blue wing and not in the red wing implies strong upflows of at least 50 km s−1 (points further in the wings also shows similar features). These dark jets could then be identified as surges, coexisting with the hot fibrils (Nóbrega-Siverio et al. 2017).

The intense chromospheric emission is found during the entire time-series of almost 2 hours. Although the evolution in the photosphere is very slow and the presence of chromospheric emission is relatively stable in time, the Ca II images also show a dynamic response with rapidly moving brightenings that evolve on smaller timescales than the CHROMIS cadence of 9 s.

3. Data analysis

We estimated the model atmosphere for the ROI through NLTE inversions. We performed an inversion of the Stokes profiles in the Fe I 6301.5 and 6302.5 Å, Ca II 8542 Å, and Ca II K lines simultaneously, using the parallel NLTE STockholm Inversion Code2 (STiC; de la Cruz Rodríguez et al. 2016, 2019). This code is built on top of an optimized version of the RH code (Uitenbroek 2001). It makes use of a fast angular approximation to include the effect of partial redistribution (PRD; Leenaarts et al. 2012) which is necessary for a realistic modeling of the Ca II K line formation. Electron densities were derived consistent with NLTE hydrogen ionization, by iteratively solving the statistical equilibrium equations whilst enforcing charge conservation. The polarized radiative transfer equation was integrated using an efficient cubic DELO-Bezier formal solver (de la Cruz Rodríguez & Piskunov 2013). A 3D model atmosphere is obtained assuming 1D plane-parallel atmospheres along each LOS. This approximation is usually referred to as 1.5D modeling.

STiC iteratively adjusts the physical parameters of a model atmosphere, such as the temperature, LOS velocity, magnetic field vector, and microturbulence to find a synthetic spectrum that reproduces the observed profile. The stratifications of density and gas pressure are computed by assuming hydrostatic equilibrium (HE). The physical parameters of the model atmosphere are given as functions of the optical depth scale at 5000 Å, hereafter log(τ). The physical parameters are modified at specific node locations, followed by an interpolation to all other depth points. The merit function that accounts for the likelihood between the synthetic and observed spectra includes an additional regularization term that ensures a smooth behavior in the atmosphere even when a physical parameter has many degrees of freedom (i.e., nodes in the stratification).

We treated the Ca II atom in NLTE with the Ca II K line in PRD, the Ca II 8542 Å line in complete frequency redistribution (CRD), and the Fe I lines also in NLTE in case the complex atmosphere affects their formation (Smitha et al. 2020). We initialize the model atmosphere from the FAL-C model (Fontenla et al. 1993) by interpolation to 54 depth points from log(τ) = − 7.0 to log(τ) = + 1.0. We used an enhanced gas pressure at the upper boundary that accounts for the typical values of active regions for all pixels (Ptop= 1.0 dyn cm−2; de la Cruz Rodríguez et al. 2019). The inversions were run in multiple cycles with spatial smoothing between cycles, with the general approach being to use fewer nodes in the firsts cycle to obtain the general trend and more nodes in the subsequent cycles to get a more detailed atmospheric structure. We initialized the magnetic field by using a Milne-Eddington inversion3 of the Fe I lines. The inversion strategy is summarized in Table 1. To obtain an estimate of the uncertainty of our model after the inversion, we used Eq. (42) of del Toro Iniesta & Ruiz Cobo (2016) which takes into account both the discrepancy between the observed and the synthetic profiles and the effective sensitivity of the spectral lines to the underlying parameters, except that we only introduce perturbations at the node location. The uncertainty is calculated at the node location and then interpolated to the whole stratification for better visualization.

Table 1.

Number of nodes used in each inversion cycle.

4. Results

4.1. One snapshot inversion

The inversion of the whole FOV and/or the complete temporal series is very computationally expensive. Therefore, we only inverted our ROI around 12:15 UT. This snapshot includes the observation of the Fe I lines (not included in the other two time-series), exhibits good seeing in all wavelength positions for all spectral lines, and the time difference between the observation of the central wavelength of Ca II K and Ca II 8542 Å is less than 1 s. As the analysis is shown with maps at a given optical depth, regions at equal τ values can sample different geometrical heights.

Temperature stratification. Figure 4 shows horizontal slices of the inferred temperature at three different optical depths: at log(τ) ∼ 0, associated with the continuum formation layer in the photosphere, log(τ) = − 3.2 which shows the upper photosphere (around the solar temperature minimum in quiet-sun models), and at log(τ) = − 4, which provides a view of the chromosphere. At log(τ) ∼ 0, the temperature shows the pattern of photospheric granulation with temperatures around 6000 K. In this layer, there is no indication of chromospheric activity. At log(τ) ∼  − 3.2, we find a region of about 6″  × 6″ with a temperature of around 7000 K compared with the ∼5000 K of the surrounding. It starts at the edge of the pore where magnetic flux cancelation is expected to occur and extends towards the right side coinciding in location with the bright fibrils (see Fig. 3).

thumbnail Fig. 4.

Atmospheric structure in the ROI as inferred from the inversion. The first and second rows show the temperature and the LOS velocity at three layers, and the third and fourth rows show the longitudinal magnetic field, the transverse magnetic field, and the microturbulent velocity at two layers. In the LOS velocity maps, the blue color represents motions toward the observer (upflows) while the red color represents motions from the observer (downflows). Crosses in the upper-left panel show the location of the pixel examples of Fig. 5 and the solid line the cross-cut shown in Fig. 12. Thin dashed and dotted contours are contours of B||(log(τ) = 0) = ± 200 G.

The atmosphere above (at log(τ) = − 4) shows a very different temperature distribution. The region of the fibrils have now temperatures similar to what is found above the pore. In both places, some areas exceed 10 kK. Outside these areas, temperatures are much lower. While the higher temperatures above the fibrils can be understood as a transfer of energy from below, the increase in temperature above the pore has no counterpart below and could be generated at higher altitudes.

Velocity stratification. The LOS velocity maps are shown below the temperature panels in Fig. 4. The velocity map at log(τ) ∼ 0 displays the typical granulation pattern, with upflows in the granules and downflows in the intergranular lanes. At log(τ) = − 3.2, the low-chromospheric velocity map shows downflows and upflows with peaks of the order of ±10 km s−1. Downward motions are found above the positive polarity whereas upward motions are right next to the downward motions. As the plasma flows tend to align along magnetic field lines, the apparent upward and downward motions located next to each other suggest the presence of bidirectional flows along inclined magnetic field lines with respect to the LOS.

There are also some downward motions near (14″, 7″) close the other small negative polarity, which could be hot plasma that is not able to escape and is redirected downwards again along low-lying field lines. Finally, one of the surges along y = 3″ appears to show upward motions. Higher in the atmosphere at log(τ) = − 4.0 the general trend is similar; the upflows persist while the downflows are less intense.

Magnetic field stratification. The inferred components of the magnetic field relative to the LOS (B and |B|) are also presented in Fig. 4. The longitudinal magnetic field B in the photosphere at log(τ) = − 1.0 shows the pre-existing pore (at the left side of the ROI) with a longitudinal magnetic field larger than −1000 G. The emergence of the magnetic bubble brings to the surface a positive polarity patch that interacts with the pore, and a negative polarity concentration visible in the lower-right part of the region (14″, 4″). In the chromosphere at log(τ) = − 3.2, the magnetic field is more diffuse and weaker although maintaining the general topology of the photospheric field. While in the pore the longitudinal magnetic field decreases to 800 G in the chromosphere, the longitudinal magnetic field in the location of the bright fibrils is very weak.

The transverse magnetic field |B| in the photosphere at log(τ) = − 1.0 correlates with the areas of a strong longitudinal magnetic field (highlighted in this panel with contours at ±200 G). The transverse field reaches more than 1.5 kG in the right edge of the pore, and decreases rapidly outside the latter. The direction of the inferred azimuth (ϕB) has also been drawn on this panel. Although it has not been corrected from its 180-degree ambiguity, it is sufficient to illustrate that the general direction of the field is radial to the interior of the pore. In the chromosphere at log(τ) = − 3.2, the transverse magnetic field is in general much weaker, except in the region between the pore and the new polarity around (7″, 10″) where the field increases up to 700 G compared to the very weak photospheric field at this location. At the location of the fibrils, the transverse magnetic field is of the order of 200 G. Given that the longitudinal component is around 120 G, the inclination of the field lines would be around 30 degrees. However, given the low field strength and the noise present in the signals, we cannot ensure the reliability of this estimation.

Microturbulent velocity. In some cases, STiC is not capable of finding a temperature stratification that is capable of explaining the width of these profiles without including a microturbulent contribution to the Doppler width. In the photosphere, the microturbulence has values below 1 km s−1 on average, reaching ∼2 km s−1 in the intergranules and near the boundaries of the pore. This has been associated with turbulent motions in the intergranular lanes (Borrero & Bellot Rubio 2002). In the chromosphere (between log(τ) = − 3 and log(τ) = − 4) the microturbulent velocity increases on average up to 2–3 km s−1 in places far from the cancelation site. In locations with strong downflows, it reaches higher values (3–4 km s−1). This increase could be due to the turbulent motions as the plasma is being compressed towards the denser lower atmosphere. On the other hand, we find that the plasma associated with upward motions does not require such additional broadening (see Fig. 4). This might indicate how the plasma can escape easier to higher layers of the atmosphere where the surrounding density is lower.

Spectral signatures. To show the shape of the typical observed Stokes profiles and the quality of the fits, some profiles have been extracted from the location marked in the temperature map at log(τ) ∼ 0 in Fig. 4 and shown in Fig. 5. Overall, the fits to the Stokes profiles follow the observations closely. The first pixel 5a was extracted from the region of the high transverse chromospheric magnetic field. In this location, the wings of the Ca II lines are in emission while the core is in absorption. Although weak, the Ca II 8542 Å signals are above the noise and have been reproduced by the inversion code.

thumbnail Fig. 5.

Examples of spectra (observed in black and synthetic in red) with the corresponding inferred atmosphere. Red shaded bands indicate node uncertainties. The location of each pixel is quoted in each panel. The gray dashed curve represents the average quiet-Sun temperature stratification. (a) Pixel (6.6″, 9.0″) is shown as blue ‘+’ symbol in Fig. 4. (b) Pixel (9.7″, 10.8″) is shown as black ‘x’ symbol in Fig. 4.

One more profile was extracted from the bright fibrils and it is shown in Fig. 5b. The line-core intensities of Ca II 8542 Å and Ca II K are high, up to 13 and 3 times compared to their typical values at the quiet Sun, even higher than during some strong flares (Kuridze et al. 2018; Yadav et al. 2021). The retrieved atmosphere shows a strong temperature enhancement in the chromosphere compared to the profile in Fig. 5a. One consequence of that is the loss of sensitivity higher up as the material ionizes. This behavior can be detected in the uncertainties calculated for this pixel compared with the previous pixel. We note that our estimate of the uncertainties does not include degeneracies between different parameters and does not allow us to estimate errors included in the assumptions of the models such as the 1.5D plane-parallel approximation, and the assumption of hydrostatic equilibrium, and so on. Our estimate for one parameter is purely based on the goodness of the fit and the sensitivity of the different spectral lines in all Stokes parameters to a few nodes in the atmosphere.

4.2. Evolution of the polarization signals in time

Although we do not discuss in detail the overall temporal evolution of our target, we briefly discuss the evolution of the polarization signals. Figure 6 shows the evolution of the linear polarization map of the Ca II 8542 Å line, and therefore the evolution of the chromospheric magnetic field indirectly. Some contours show the temporal evolution of the circular polarization signals in the photosphere (given by the Fe I 6301 Å line). The first panel resembles the transverse magnetic field map shown in Fig. 4, indicating that the inversion code was able to extract this chromospheric information.

thumbnail Fig. 6.

Linear polarization maps at the central wavelength of the Ca II 8542 Å line. The polarity is highlighted with contours at −3 × 10−3Ic and 5 × 10−4Ic of Stokes V in the red wing of the Fe I 6301 Å line. The time of the observations is shown in the lower-left corner of each panel.

Two main conclusions can be drawn from Fig. 6: the circular polarization is decreasing as the magnetic flux is canceled, while the linear polarization increases in amplitude and area. This enhanced linear polarization is consistent with a ∩-configuration as the field lines are being bent or reconfigured during the magnetic reconnection. These polarization signals indicate a transverse magnetic field of 1 kG in a large part of the region. These unusually strong transverse fields (more than 1 kG) have also been found in the chromosphere during magnetic reconnection events (Kuridze et al. 2018).

4.3. Plasmoid-mediated magnetic reconnection

During the evolution of this event, we detected the presence of many small-scale blobs with similar shapes and widths, between 0.15 and 0.3″, moving away from the cancelation site, identified in the blue wing (λ0 − 1.3 Å) of the Ca II K line. These blobs were also seen in the solar atmosphere (Innes et al. 2015; Rouppe van der Voort et al. 2017), and were identified as plasmoids. They are formed when the current sheet exceeds a certain aspect ratio, becomes unstable, and fragments into smaller elements (plasmoids). If the latter are detected during at least two time steps we can track them and infer their apparent speed in the plane of the sky. Figure 7 shows three examples of these blobs at different time frames with the inferred velocities. We detect fast transverse velocities between 70 km s−1 and 100 km s−1. Their high velocities allow us to distinguish them more easily from the surroundings at this wavelength where they are clearly visible over a more photospheric background. Given that both the LOS and plane-of-the-sky velocities are of the same order, the ejected plasma appears to be neither predominantly vertical nor horizontal, supporting the idea of a tilted current sheet assuming that the flows are parallel to that. The direction of the motions is pointing away from the main polarity (pore) that the small opposite polarity is moving into. This is consistent with other observations and simulations (Chitta et al. 2017; Peter et al. 2019) as the stronger vertical magnetic field of the dominant polarity might be more difficult to bend. We identified more than 25 similar blobs in the three time-series (detected in two consecutive frames). We also find hints of other blobs, but they are more difficult to identify as they are seen for less than 9 s (one frame). This means that on average we have more than 1 blob every 30 s. This value is compatible with the timescale of plasmoid formation found in simulations (Nóbrega-Siverio et al. 2016). The presence of plasmoids further supports the magnetic reconnection scenario.

thumbnail Fig. 7.

Monochromatic images in the blue wing of Ca II K (λ0 − 1.3 Å) at three temporal intervals. Each row shows a different blob. Red circles and arrows indicate the location and direction of the motion. FWHM inferred from the blobs are quoted in the right corner of the left panels. Transverse velocities are shown in the right panels.

4.4. Energy balance

We can estimate the radiative cooling rate per unit volume from our inferred models as the frequency-integrated radiative flux divergence: ∫∇Fνdν. We approximate this by summing the net cooling rates in bound–bound transitions for the spectral lines in our model atoms for Ca II and Mg II:

Q = h ν 0 ( n u R ul n l R lu ) , $$ \begin{aligned} Q = h \nu _0 (n_u R_{ul} - n_l R_{lu}), \end{aligned} $$(1)

where h is the Planck’s constant, ν0 the line center frequency, nu and nl the population of each level, and Rul the radiative rate coefficient from level u to level l. This is an approximation: in the photosphere we ignore contributions from many weak spectral lines and the H continuum, and in the upper chromosphere we ignore Lyα and the transitions of helium. In the low and mid-chromosphere, this approximation is expected to be rather accurate: the heating in the Balmer continuum and the cooling in the Hα line tend to cancel each other, and other transitions tend to contribute much less to the cooling than Ca II and Mg II (Vernazza et al. 1981).

The height distribution of the total radiative cooling rate per unit volume (Q) and per unit mass (Q/ρ) is shown in Fig. 8 for three different patches inside the ROI. We chose a quiet-sun region (A), a region close to the fibrils (B), and a region inside the brightest fibrils (C), all indicated in Fig. 9. As expected, higher temperature regions in the atmosphere correlate with enhanced radiative cooling. The radiative losses are largest at the bottom of the chromosphere because the density is also much larger deeper down.

thumbnail Fig. 8.

Height distribution of the total radiative cooling rate per unit volume Q (upper panel) and per unit mass Q/ρ (lower panel) in solid lines calculated as the average of the regions highlighted in Fig. 9 with their average temperature stratification (dashed lines).

thumbnail Fig. 9.

Total integrated radiative losses calculated from the semi-empirical model inferred in the inversion.

thumbnail Fig. 10.

Joint probability distributions of the wavelength-summed Ca II K intensity and the total integrated radiative losses computed from the data shown in Fig. 9. The color of each bin is proportional to the logarithm of the number of points, with darker color indicating more points.

The total radiative losses in each pixel have been integrated over geometrical height (which is derived assuming hydrostatic equilibrium in the inversion) between the height of the temperature minimum and the height where the temperature reaches 10 kK. Above that temperature, the Ca II lines lose sensitivity and the hydrogen Lyα line starts to contribute, and so our derived cooling rates are inaccurate there. The total integrated radiative losses for all the pixels of the FOV are shown in Fig. 9, and the contribution of each spectral line to these losses is given in Table 2.

Table 2.

Average contribution of each spectral line to the total integrated radiative losses for the regions highlighted in Fig. 9, in units of kW m−2.

In the quiet sun (region A) we obtain around 4.5 kW m−2, while in the intermediate area (region B) we obtain losses of 32 kW m−2. These values agree well with the previous estimates (Vernazza et al. 1981; Withbroe & Noyes 1977). We emphasize that our cooling rate estimates are more consistent with the real situation: they are computed for each spatially resolved atmosphere and then averaged, whereas previous studies first averaged observed line profiles, then inferred an atmosphere consistent with the average line profile, and finally computed the cooling rate from that atmosphere.

In region C, the average integrated radiative losses are 108 kW m−2 but Fig. 9 shows peak values up to 160 kW m−2. These values are roughly a factor five higher than the semi-empirical values for active regions, indicating that the averaging of the observations can lead to substantial underestimation of the local radiative cooling rate.

Finally, we show the correlation of the height-integrated radiative losses with the wavelength-integrated intensity in the Ca II K line. The wavelength integration was done between the K1 minima whose formation is close to the temperature minimum (Bjørgen et al. 2018). The correlation is strong and close to linear, with a Pearson coefficient of 0.97. This result supports the analysis and results of Leenaarts et al. (2018), who used wavelength-integrated Ca II K intensity as a proxy for chromospheric radiative losses.

As the energy lost during the whole observation needs to be continually replenished, the total integrated radiative losses in the ROI assuming a similar emission during the two hours is around ∼5 × 1027 erg. Motivated by the magnetic nature of this event, we propose that the energy required for this event can be extracted from magnetic reconnection in the chromosphere. We can use the magnetic field vector from the inversion to estimate the magnetic energy content of the atmosphere, by integrating the whole volume 1/8πB2dV including the stratification in height. In that case, the energy content of this volume is ∼5 × 1029 erg. We set the integration range from log(τ) = 0 to log(τ) = − 4 which corresponds approximately to a height range of Δz = 800 km where we have enough magnetic sensitivity.

4.5. Magnetic field topology

The analysis of the monochromatic images and the inversion results reveals several aspects that can be put together to propose a magnetic topology that is consistent with the following findings:

1. The emergence of a new small magnetic flux concentration causes the footpoints of opposite polarity to separate and the granulation is deformed accordingly. This indicates that the flux tube adopts the shape of an Ω loop.

2. The positive polarity of a newly emerged magnetic concentration interacts with a pre-existing pore of a fully developed active region.

3. In the location where these two polarities come into contact and cancelation of the magnetic flux occurs, strong emission is found in all chromospheric and coronal diagnostics.

4. The temperature inferred by the inversion code indicates that in the lower chromosphere and above, the temperature increases for more than 2000 K while no signs of activity are present in the photosphere.

5. Bidirectional jets (accelerated to ±10 km s−1) are found with downward motions located above the positive polarity and upward motions towards the negative polarities. This indicates that the plasma flows are very inclined with respect to the LOS. The apparent motion in the plane of the sky (detected in the monochromatic images) is also pointing away from the pore.

6. After reaching a certain height, given the topology of the field, the plasma above seems to move toward the left side of the FOV as indicated by the surges. According to Fig. 1, the magnetic field at higher layers might be connected to the sunspot of the left side of the active region.

All the above findings are suggestive of a model magnetic reconnection between an emerging magnetic field and the pre-existing field, as found in observations as well as simulations (Nóbrega-Siverio et al. 2016; Hansteen et al. 2019; Guglielmino et al. 2019). The topology of our event is expected to be very similar to the classical emerging flux model for solar reconnection (Heyvaerts et al. 1977; Yokoyama & Shibata 1995) but on a much smaller scale where the event occurs deeper in the lower chromosphere. In this scenario, a magnetic diffusion region would arise around the interface between the two magnetic polarities. Inside this diffusion region, free magnetic energy is released in several different ways: the electric field induced by the magnetic field creates a thin layer of an intense current sheet, which heats the local plasma through Joule heating. In addition, plasma is ejected away from the reconnection site by the magnetic tension force. Another difference with those models is that the magnetic field at higher layers in many of the simulated cases is nearly vertical, while here the numerous plasma flows are in the plane of the sky and elongated structures indicate a nearly horizontal canopy of the pre-existent magnetic field.

Figure 11 shows a simplified illustration of the proposed magnetic topology that is consistent with our results. This figure schematically shows: (a) the emerged magnetic loop and the motion of the footpoints (with orange arrows), (b) the pore with strong vertical fields (negative concentration on the left side), (c) the reconnection region (with a purple cross) with the corresponding heating (orange grid) and bidirectional flows (red and blue arrows), and (d) later plasma motions which escape from the region (gray arrow) or which fall again (right red arrow). To compare the similarities of our 2D sketch, we also show in Fig. 12 a cross-cut of the inversion results across the polarities. In these figures, some vertical lines highlight the different domains and their correspondence in both figures. In this figure, the extended heating occurring at log(τ) = − 3.5 is clearly visible, similarly to Leenaarts et al. (2018), which is located deeper in regions with downflows. The transition region is reconstructed at higher optical depth in the region showing temperature enhancements. This illustration also helps to understand how the reconnection site (dashed line) could take place to the west side of the cancelation site (dotted line) from a geometrical perspective. We also added to the bottom panel of Fig. 12 the location of maximum sensitivity of the linear polarization to the transverse magnetic field in the Ca II 8542 Å line. On average, it is placed around log(τ) = − 4 in quiet-sun areas and log(τ) = − 3 in the flaring region due to the temperature enhancement, as shown in Yadav et al. (2021).

thumbnail Fig. 11.

Schematic illustration of the proposed magnetic field topology. Suggested magnetic field lines are given in black. Upflows and downflows are colored in blue and red, respectively. Vertical dotted and dashed lines are drawn for comparison with Fig. 12.

thumbnail Fig. 12.

Vertical stratification across the solid line indicated in the upper-left panel of Fig. 4. From top to bottom: temperature, LOS velocity, and longitudinal magnetic field.

The dark jets observed in Hα at high velocity seem to match with the description of surges, that is, cool chromospheric jets that coexist with the hot jets (Takasao et al. 2013; Nóbrega-Siverio et al. 2016). Nevertheless, it has been shown that the simple reconnection as a driver seems insufficient to explain their higher altitude (barely detected in the Ca II lines), longitude (compared with the size of the hot fibrils), and velocities (compared with the reconnection outflow velocity). Numerical simulations suggest that shocks are the acceleration mechanism (Takasao et al. 2013; Shibata et al. 2007). The magnetic reconnection happening in the lower atmosphere can generate slow mode waves. These waves propagate into regions of lower density where eventually they form a shock as a result of the law of conservation of energy which accelerates the local plasma upwards (Takasao et al. 2013). These waves can also be generated when reconnection outflows hit the ambient magnetic field (right side of Fig. 11). Given the large timescale of our event, a recurrent generation of waves seems to match with our scenario, which was also suggested by Takasao et al. (2013). Observations at higher cadence are needed to confirm the wave and shocks scenario as the origin.

Given the dynamic nature and spectral signatures of this observation, an alternative explanation could be given by a flaring atmosphere scenario. In the standard flare model, the energy released in the corona can be transported and dissipated into the chromosphere and even the photosphere. This scenario seems unlikely because (a) the new small magnetic flux emergence should start interacting with the pre-existing canopy of the pore from the bottom of the atmosphere, (b) the shape of the light curves does not show a sudden peak and gradual decrease as in standard flares but rather a gradual evolution during the cancelation period, and (c) it would be a very small flare (as it is not detected in GOES) and small flares are unlikely to produce the deep-seated and persistent heating for almost two hours.

5. Summary and conclusions

Using multi-line NLTE inversions, we retrieved a 3D atmospheric model of a reconnection event caused by the emergence of new magnetic flux into a previously existing magnetic field configuration in the chromosphere. The addition of four lines with respect to single-line inversions and in particular the Ca II K, allows us to improve the reconstruction of the atmospheric parameters to greater heights, while better constraining the areas of the model where the sensitivities of several lines overlap. Based on the results of the inversion and the appearance of the observed target in Hα, we propose a likely magnetic field configuration that is compatible with all observational signatures and with the inversion results. In the region where the reconnection occurs, we infer a temperature increase in the low chromosphere that is about 2000 K higher than in the surrounding quiet Sun, whereas higher up in the chromosphere we find a 1000–3000 K increase in temperature in a more extended region surrounding the reconnection event. In both cases, we find a fine spatial structure down to subarcsecond scales. We also find bidirectional flows of the order of 10 km s−1 emanating from the region, as well as bright blobs of roughly 0.2 arcsec diameter in the Ca II K ejected from the same region at 100 km s−1.

Increased brightness in the 304 Å and 171 Å channels of SDO/AIA indicate heating associated with the reconnection event to transition region and coronal temperature. Our inversion is not sensitive to temperatures above approximately 10 000 K and so we cannot set constraints on the height where the heating observed in AIA is located. A coronal response is not normally seen in photospheric and low-chromospheric reconnection events such as Ellerman bombs, but heating to transition-region temperatures in reconnection events higher in the chromosphere can give a response in AIA if the column density of overlying cool material is low enough so that not all radiation is absorbed in the hydrogen and helium continua (Young et al. 2018; Guglielmino et al. 2019; Hansteen et al. 2019; da Silva Santos et al. 2020).

We estimated the chromospheric radiative losses from the inferred atmospheric model in the reconnection target. The values that we obtain for the quiet Sun surrounding the active region are of the order of ∼4.5 kW m−2, consistent with earlier estimates based on semi-empirical models of the quiet Sun (Vernazza et al. 1981; Withbroe & Noyes 1977). In the active region, but outside of the reconnection event, we find values of 20–50 kW m−2, which is comparable with earlier modeling. However, we also note substantial variation with location on the Sun at arcsecond scales. Inside the reconnection event, we find peak values of the radiative losses up to 160 kW m−2, with spatial variations on subarcsecond scales. These losses are a factor five or more than the typical value for the chromosphere in active regions.

The majority of the radiative losses occur at relatively low heights, around 300−500 km above the log(τ) = 0 surface, and are caused by a deep-seated chromospheric temperature rise. Whether these losses are counterbalanced by in-situ conversion from magnetic energy, or caused by relatively hot flows coming from larger heights is still an open question. The magnetic energy release estimated from the flux cancelation (∼5 × 1029 erg) is sufficient to sustain the chromospheric losses for more than two hours (∼5 × 1027 erg).

The spatial structure of, and the quantitative constraint on, the radiative losses will be of help to validate models that aim to reproduce reconnection in the lower solar atmosphere. Such characterization is important for understanding reconnection processes and the triggering role of turbulent versus plasmoid reconnection. The advent of next-generation solar telescopes such as the Daniel K. Inouye Solar Telescope (DKIST; Tritschler et al. 2015) and the European Solar Telescope (EST; Matthews et al. 2016) will be crucial in this regard.


1

These routines are available in the Python repository

https://github.com/ISP-SST/ISPy

Acknowledgments

We would like to thank the anonymous referee for their comments and suggestions. C.J.D.B. thanks Sanja Danilovic and Flavio Calvo for their comments. J.dl C.R. is supported by grants from the Swedish Research Council (2015-03994), the Swedish National Space Agency (128/15), and the Swedish Civil Contingencies Agency (MSB). This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (SUNMAG, grant agreement 759548). J.L. is supported through a project (CHROMATIC, 2016.0019) funded by the Knut and Alice Wallenberg foundation. The Swedish 1-m Solar Telescope is operated on the island of La Palma by the Institute for Solar Physics of Stockholm University in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias. The Institute for Solar Physics is supported by a grant for research infrastructures of national importance from the Swedish Research Council (registration number 2017-00625). We acknowledge the community effort devoted to the development of the following open-source packages that were used in this work: numpy (numpy.org), matplotlib (matplotlib.org), scipy (scipy.org), astropy (astropy.org) and sunpy (sunpy.org). Infrastructure for Computing (SNIC) at the PDC Centre for High Performance Computing (Beskow, PDC-HPC) at the Royal Institute of Technology in Stockholm as well as recourses at the National Supercomputer Centre (Tetralith, NSC) at Linköping University. This research has made use of NASA’s Astrophysics Data System Bibliographic Services.

References

  1. Archontis, V. 2012, Phil. Trans. Roy. Soc. London Ser. A, 370, 3088 [Google Scholar]
  2. Bjørgen, J. P., Sukhorukov, A. V., Leenaarts, J., et al. 2018, A&A, 611, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Borrero, J. M., & Bellot Rubio, L. R. 2002, A&A, 385, 1056 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Chitta, L. P., & Lazarian, A. 2020, ApJ, 890, L2 [CrossRef] [Google Scholar]
  5. Chitta, L. P., Peter, H., Young, P. R., & Huang, Y. M. 2017, A&A, 605, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Cox, A. N. 2000, Allen’s Astrophysical Quantities (New York: AIP Press), 700 [Google Scholar]
  7. da Silva Santos, J. M., de la Cruz Rodríguez, J., White, S. M., et al. 2020, A&A, 643, A41 [EDP Sciences] [Google Scholar]
  8. Danilovic, S. 2017, A&A, 601, A122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. de la Cruz Rodríguez, J., & Piskunov, N. 2013, ApJ, 764, 33 [NASA ADS] [CrossRef] [Google Scholar]
  10. de la Cruz Rodríguez, J., Löfdahl, M. G., Sütterlin, P., Hillberg, T., & Rouppe van der Voort, L. 2015, A&A, 573, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. de la Cruz Rodríguez, J., Leenaarts, J., & Asensio Ramos, A. 2016, ApJ, 830, L30 [NASA ADS] [CrossRef] [Google Scholar]
  12. de la Cruz Rodríguez, J., Leenaarts, J., Danilovic, S., & Uitenbroek, H. 2019, A&A, 623, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. del Toro Iniesta, J. C., & Ruiz Cobo, B. 2016, Liv. Rev. Sol. Phys., 13, 4 [Google Scholar]
  14. Díaz Baso, C. J., Martínez González, M. J., Asensio Ramos, A., & de la Cruz Rodríguez, J. 2019a, A&A, 623, A178 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Díaz Baso, C. J., de la Cruz Rodríguez, J., & Danilovic, S. 2019b, A&A, 629, A99 [CrossRef] [EDP Sciences] [Google Scholar]
  16. Ellerman, F. 1917, ApJ, 46, 298 [NASA ADS] [CrossRef] [Google Scholar]
  17. Fontenla, J. M., Avrett, E. H., & Loeser, R. 1993, ApJ, 406, 319 [Google Scholar]
  18. Guglielmino, S. L., Zuccarello, F., Young, P. R., Murabito, M., & Romano, P. 2018, ApJ, 856, 127 [NASA ADS] [CrossRef] [Google Scholar]
  19. Guglielmino, S. L., Young, P. R., & Zuccarello, F. 2019, ApJ, 871, 82 [NASA ADS] [CrossRef] [Google Scholar]
  20. Hansteen, V. H., Archontis, V., Pereira, T. M. D., et al. 2017, ApJ, 839, 22 [NASA ADS] [CrossRef] [Google Scholar]
  21. Hansteen, V., Ortiz, A., Archontis, V., et al. 2019, A&A, 626, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Henriques, V. M. J. 2012, A&A, 548, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Heyvaerts, J., Priest, E. R., & Rust, D. M. 1977, ApJ, 216, 123 [NASA ADS] [CrossRef] [Google Scholar]
  24. Innes, D. E., Guo, L. J., Huang, Y. M., & Bhattacharjee, A. 2015, ApJ, 813, 86 [NASA ADS] [CrossRef] [Google Scholar]
  25. Jeffrey, N. L. S., Fletcher, L., Labrosse, N., & Simões, P. J. A. 2018, Sci. Adv., 4, 2794 [CrossRef] [Google Scholar]
  26. Kowal, G., Lazarian, A., Vishniac, E. T., & Otmianowska-Mazur, K. 2009, ApJ, 700, 63 [NASA ADS] [CrossRef] [Google Scholar]
  27. Kuridze, D., Henriques, V. M. J., Mathioudakis, M., et al. 2018, ApJ, 860, 10 [NASA ADS] [CrossRef] [Google Scholar]
  28. Lazarian, A., & Vishniac, E. T. 1999, ApJ, 517, 700 [NASA ADS] [CrossRef] [Google Scholar]
  29. Leenaarts, J., Pereira, T., & Uitenbroek, H. 2012, A&A, 543, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Leenaarts, J., de la Cruz Rodríguez, J., Danilovic, S., Scharmer, G., & Carlsson, M. 2018, A&A, 612, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Lemen, J. R., Title, A. M., Akin, D. J., et al. 2012, Sol. Phys., 275, 17 [Google Scholar]
  32. Libbrecht, T., Joshi, J., Rodríguez, J. d. l. C., Leenaarts, J., & Ramos, A. A. 2017, A&A, 598, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Libbrecht, T., Bjørgen, J. P., Leenaarts, J., et al. 2020, A&A, submitted [arXiv:2010.15946] [Google Scholar]
  34. Löfdahl, M. G., Hillberg, T., de la Cruz Rodriguez, J., et al. 2018, ArXiv e-prints [arXiv:1804.03030] [Google Scholar]
  35. Matthews, S. A., Collados, M., Mathioudakis, M., & Erdelyi, R. 2016, in The European Solar Telescope (EST), SPIE Conf. Ser., 9908, 990809 [Google Scholar]
  36. Neckel, H., & Labs, D. 1984, Sol. Phys., 90, 205 [NASA ADS] [CrossRef] [Google Scholar]
  37. Ni, L., Kliem, B., Lin, J., & Wu, N. 2015, ApJ, 799, 79 [NASA ADS] [CrossRef] [Google Scholar]
  38. Nóbrega-Siverio, D., Moreno-Insertis, F., & Martínez-Sykora, J. 2016, ApJ, 822, 18 [NASA ADS] [CrossRef] [Google Scholar]
  39. Nóbrega-Siverio, D., Martínez-Sykora, J., Moreno-Insertis, F.& Rouppe van der Voort, L. 2017, ApJ, 850, 153 [NASA ADS] [CrossRef] [Google Scholar]
  40. Ortiz, A., Hansteen, V. H., Nóbrega-Siverio, D., & van der Voort, L. R. 2020, A&A, 633, A58 [CrossRef] [EDP Sciences] [Google Scholar]
  41. Pesnell, W. D., Thompson, B. J., & Chamberlin, P. C. 2012, Sol. Phys., 275, 3 [Google Scholar]
  42. Peter, H., Tian, H., Curdt, W., et al. 2014, Science, 346, 1255726 [Google Scholar]
  43. Peter, H., Huang, Y. M., Chitta, L. P., & Young, P. R. 2019, A&A, 628, A8 [CrossRef] [EDP Sciences] [Google Scholar]
  44. Priest, E. 2014, Magnetohydrodynamics of the Sun (Cambridge, UK: Cambridge University Press) [Google Scholar]
  45. Rouppe van der Voort, L., De Pontieu, B., Scharmer, G. B., et al. 2017, ApJ, 851, L6 [NASA ADS] [CrossRef] [Google Scholar]
  46. Scharmer, G. 2017, SOLARNET IV: The Physics of the Sun from the Interior to the Outer Atmosphere, 85 [Google Scholar]
  47. Scharmer, G. B., Bjelksjo, K., Korhonen, T. K., Lindberg, B., & Petterson, B. 2003, in Innovative Telescopes and Instrumentation for Solar Astrophysics, eds. S. L. Keil, & S. V. Avakyan, Proc. SPIE, 4853, 341 [Google Scholar]
  48. Scharmer, G. B., Narayan, G., Hillberg, T., et al. 2008, ApJ, 689, L69 [Google Scholar]
  49. Scherrer, P. H., Schou, J., Bush, R. I., et al. 2012, Sol. Phys., 275, 207 [Google Scholar]
  50. Shibata, K., Nakamura, T., Matsumoto, T., et al. 2007, Science, 318, 1591 [NASA ADS] [CrossRef] [Google Scholar]
  51. Smitha, H. N., Chitta, L. P., Wiegelmann, T., & Solanki, S. K. 2018, A&A, 617, A128 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Smitha, H. N., Holzreuter, R., van Noort, M., & Solanki, S. K. 2020, A&A, 633, A157 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Socas-Navarro, H., Martínez Pillet, V., Elmore, D., et al. 2006, Sol. Phys., 235, 75 [NASA ADS] [CrossRef] [Google Scholar]
  54. Takasao, S., Isobe, H., & Shibata, K. 2013, PASJ, 65, 62 [NASA ADS] [CrossRef] [Google Scholar]
  55. Toriumi, S., Katsukawa, Y., & Cheung, M. C. M. 2017, ApJ, 836, 63 [NASA ADS] [CrossRef] [Google Scholar]
  56. Tritschler, A., Rimmele, T. R., Berukoff, S., et al. 2015, in 18th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun, Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun, 18, 933 [NASA ADS] [Google Scholar]
  57. Uitenbroek, H. 2001, ApJ, 557, 389 [NASA ADS] [CrossRef] [Google Scholar]
  58. van Noort, M., Rouppe van der Voort, L., & Löfdahl, M. G. 2005, Sol. Phys., 228, 191 [Google Scholar]
  59. Vernazza, J. E., Avrett, E. H., & Loeser, R. 1981, ApJS, 45, 635 [Google Scholar]
  60. Vissers, G. J. M., Rouppe van der Voort, L. H. M., Rutten, R. J., Carlsson, M., & De Pontieu, B. 2015, ApJ, 812, 11 [NASA ADS] [CrossRef] [Google Scholar]
  61. Vissers, G. J. M., de la Cruz Rodríguez, J., Libbrecht, T., et al. 2019, A&A, 627, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Withbroe, G. L., & Noyes, R. W. 1977, ARA&A, 15, 363 [NASA ADS] [CrossRef] [Google Scholar]
  63. Yadav, R., de la Cruz Rodríguez, J., Díaz Baso, C. J., et al. 2019, A&A, 632, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Yadav, R., Baso, C. J. D., de la Cruz Rodríguez, J., Calvo, F., & Morosin, R. 2021, A&A, in press, https://doi.org/10.1051/0004-6361/202039857 [Google Scholar]
  65. Yokoyama, T., & Shibata, K. 1995, Nature, 375, 42 [NASA ADS] [CrossRef] [Google Scholar]
  66. Young, P. R., Tian, H., Peter, H., et al. 2018, Space Sci. Rev., 214, 120 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1.

Number of nodes used in each inversion cycle.

Table 2.

Average contribution of each spectral line to the total integrated radiative losses for the regions highlighted in Fig. 9, in units of kW m−2.

All Figures

thumbnail Fig. 1.

SDO/HMI continuum image and SDO/HMI LOS magnetogram of the active region NOAA 12593 at 12:15 h UT. The black square contour outlines the FOV covered by the SST and the smaller green square indicates the ROI.

In the text
thumbnail Fig. 2.

Time evolution of the ROI (green region shown in Fig. 1) as seen from SDO-HMI and SDO-AIA channels. The purple dashed-line square indicates the temporal window where we have SST observations.

In the text
thumbnail Fig. 3.

Overview images which overlap between the two instruments at four different wavelengths taken at 12:15 UT. The red arrows point to solar north and west, while the white arrows indicate the direction to the closest limb. The green square indicates our ROI which we analyzed in detail with the inversion code. The black dashed square is a close-up analyzed in Sect. 4.3.

In the text
thumbnail Fig. 4.

Atmospheric structure in the ROI as inferred from the inversion. The first and second rows show the temperature and the LOS velocity at three layers, and the third and fourth rows show the longitudinal magnetic field, the transverse magnetic field, and the microturbulent velocity at two layers. In the LOS velocity maps, the blue color represents motions toward the observer (upflows) while the red color represents motions from the observer (downflows). Crosses in the upper-left panel show the location of the pixel examples of Fig. 5 and the solid line the cross-cut shown in Fig. 12. Thin dashed and dotted contours are contours of B||(log(τ) = 0) = ± 200 G.

In the text
thumbnail Fig. 5.

Examples of spectra (observed in black and synthetic in red) with the corresponding inferred atmosphere. Red shaded bands indicate node uncertainties. The location of each pixel is quoted in each panel. The gray dashed curve represents the average quiet-Sun temperature stratification. (a) Pixel (6.6″, 9.0″) is shown as blue ‘+’ symbol in Fig. 4. (b) Pixel (9.7″, 10.8″) is shown as black ‘x’ symbol in Fig. 4.

In the text
thumbnail Fig. 6.

Linear polarization maps at the central wavelength of the Ca II 8542 Å line. The polarity is highlighted with contours at −3 × 10−3Ic and 5 × 10−4Ic of Stokes V in the red wing of the Fe I 6301 Å line. The time of the observations is shown in the lower-left corner of each panel.

In the text
thumbnail Fig. 7.

Monochromatic images in the blue wing of Ca II K (λ0 − 1.3 Å) at three temporal intervals. Each row shows a different blob. Red circles and arrows indicate the location and direction of the motion. FWHM inferred from the blobs are quoted in the right corner of the left panels. Transverse velocities are shown in the right panels.

In the text
thumbnail Fig. 8.

Height distribution of the total radiative cooling rate per unit volume Q (upper panel) and per unit mass Q/ρ (lower panel) in solid lines calculated as the average of the regions highlighted in Fig. 9 with their average temperature stratification (dashed lines).

In the text
thumbnail Fig. 9.

Total integrated radiative losses calculated from the semi-empirical model inferred in the inversion.

In the text
thumbnail Fig. 10.

Joint probability distributions of the wavelength-summed Ca II K intensity and the total integrated radiative losses computed from the data shown in Fig. 9. The color of each bin is proportional to the logarithm of the number of points, with darker color indicating more points.

In the text
thumbnail Fig. 11.

Schematic illustration of the proposed magnetic field topology. Suggested magnetic field lines are given in black. Upflows and downflows are colored in blue and red, respectively. Vertical dotted and dashed lines are drawn for comparison with Fig. 12.

In the text
thumbnail Fig. 12.

Vertical stratification across the solid line indicated in the upper-left panel of Fig. 4. From top to bottom: temperature, LOS velocity, and longitudinal magnetic field.

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.