# High-resolution observations of recurrent jets from an arch filament system

Reetika Joshi, Luc Rouppe van der Voort, Brigitte Schmieder, Fernando Moreno-Insertis, Avijeet Prasad, Guillaume Aulanier, Daniel Nóbrega-Siverio

► To cite this version:

Reetika Joshi, Luc Rouppe van der Voort, Brigitte Schmieder, Fernando Moreno-Insertis, Avijeet Prasad, et al.. High-resolution observations of recurrent jets from an arch filament system. Astronomy & Astrophysics - A&A, 2024, 691, pp.A198. 10.1051/0004-6361/202449715. hal-05207447

HAL Id: hal-05207447

https://hal.science/hal-05207447v1

Submitted on 13 Aug 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# High-resolution observations of recurrent jets from an arch filament system

Reetika Joshi<sup>1,2,*</sup>①, Luc Rouppe van der Voort<sup>1,2</sup>①, Brigitte Schmieder<sup>3,4,5</sup>①, Fernando Moreno-Insertis<sup>6,7</sup>, Avijeet Prasad<sup>1,2</sup>①, Guillaume Aulanier<sup>8,1</sup>①, and Daniel Nóbrega-Siverio<sup>6,7,1,2</sup>①

Received 23 February 2024 / Accepted 29 August 2024

## ABSTRACT

**Context.** Solar jets are collimated plasma ejections along magnetic field lines observed in hot (extreme-ultraviolet (EUV) jets) and cool (chromospheric surges) temperature diagnostics. Their trigger mechanisms and the relationship between hot and cool jets are still not completely understood.

**Aims.** We aim to investigate the generation of a sequence of active-region solar jets and their evolution from the photospheric to the coronal heights using multithermal observations from ground-based and space-borne instruments.

**Methods.** Using the synergy of high-spatial-resolution and high-temporal-resolution observations by the Swedish 1-m Solar Telescope (SST), along with the Solar Dynamics Observatory (SDO), we analyzed a sequence of solar jets originating in a mixed-polarity region between the leading and following sunspots of an active region. We investigated the kinematics of these jets using the spectra from the SST observations. We used a non-force-free field (NFFF) extrapolation technique to derive the magnetic field topology of the active region.

**Results.** A mixed-polarity region is formed over a long period (24 hours) with persistent magnetic flux emergence. This region has been observed as an arch filament system (AFS) in chromospheric SST observations. In this region, negative polarities surrounded by positive polarities create a fan surface with a null point at a height of 6 Mm detected in the NFFF extrapolation. SST observations in the Hβ spectral line reveal a large flux rope over the AFS moving from north to south, causing successive EUV and cool jets to move in the east-west direction and later towards the south along the long open loops.

**Conclusions.** The high-resolution SST observations (0.038 per pixel) resolve the dark area observed at the jet base and reveal the existence of an AFS with an extended cool jet, which may be the result of a peeling-like mechanism of the AFS. Based on the combined analysis of SST and AIA observations along with extrapolated magnetic topology, it is suggested that the magnetic reconnection site may move southward by approximately 20 Mm until it reaches a region where the open magnetic field lines are oriented north-south.

**Key words.** Sun: chromosphere – Sun: corona – Sun: flares – Sun: magnetic fields

## 1. Introduction

High-resolution observations from ground and space have revealed that the Sun is continuously generating ejections starting in the low atmosphere and reaching the corona. The generic solar jets are observed as multithermal collimated plasma flows along the magnetic field lines (see e.g., Shibata et al. 1992; Canfield et al. 1996; Nisticò et al. 2009; Chandra et al. 2015; Joshi et al. 2017). Solar jets are observed in a wide range of wavelengths in the solar spectrum and have a great variety of sizes, from the very small scale, similar to spicules (Tian et al. 2014), to large-scale eruptions leading to coronal mass ejections (Wang et al. 2002; Shen et al. 2012; Liu et al. 2015; Joshi et al. 2020c). These jets are observed in different solar domains,

such as coronal holes (Panesar et al. 2018), active regions, penumbra (Liu et al. 2022), light bridges (Asai et al. 2001; Robustini et al. 2016), and even in the quiet sun or at the pole (Nisticò et al. 2009; Raouafi & Stenborg 2014). Up- and down-flows are commonly observed in solar jets (Schmieder et al. 1997; Cheung et al. 2015; Joshi et al. 2024). Several reviews summarize the observations and simulations, both past and present (Raouafi et al. 2016; Schmieder 2022).

Solar coronal jets have been classified as standard and blowout jets (Moore et al. 2010). The triggering mechanism behind both of these classes of jets is magnetic reconnection, which allows the stored magnetic energy to be released and converted into kinetic and thermal energy flux in the reconnection outflows. Magnetic reconnection in tiny regions showing small-scale jetting activity (called jetlets) is also proposed to be a source of heating for the solar corona and to act as a driver of

<footnote n="1">¹ Rosseland Centre for Solar Physics, University of Oslo, PO Box 1029, Blindern N-0315, Oslo, Norway</footnote>

<footnote n="2">² Institute of Theoretical Astrophysics, University of Oslo, PO Box 1029, Blindern N-0315, Oslo, Norway</footnote>

<footnote n="3">³ LIRA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Université de Paris, 5 place Jules Janssen, F-92195 Meudon, France</footnote>

<footnote n="4">⁴ Centre for mathematical Plasma Astrophysics, Dept. of Mathematics, KU Leuven, 3001 Leuven, Belgium</footnote>

<footnote n="5">⁵ School of Physics and Astronomy, University of Glasgow, Glasgow G12 8QQ, UK</footnote>

<footnote n="6">⁶ Instituto de Astrofísica de Canarias, E-38200 La Laguna, Tenerife, Spain</footnote>

<footnote n="7">⁷ Departamento de Astrofísica, Universidad de La Laguna, E-38206 La Laguna, Tenerife, Spain</footnote>

<footnote n="8">⁸ Sorbonne Université, École Polytechnique, Institut Polytechnique de Paris, Observatoire de Paris – PSL, CNRS, Laboratoire de physique des plasmas (LPP), 4 place Jussieu, F-75005 Paris, France</footnote>

<footnote>* Corresponding author; reetika.joshi@astro.uio.no</footnote>
the solar wind (Raouafi & Stenborg 2014; Panesar et al. 2018). There are different theoretical explanations with regard to the drivers of the jets, which may differ depending on the presence of emerging magnetic field, magnetic flux cancelation, or erup- tion of a mini-filament.

The theoretical basis of the standard jet eruption model was first proposed by Heyvaerts et al. (1977), further devel- oped by Shibata et al. (1992), and later numerically simulated by several authors in 2D (Yokoyama & Shibata 1995, 1996) and 3D (Moreno-Insertis et al. 2008; Moreno-Insertis & Galsgaard 2013). In these models, closed magnetic arch-like bipolar flux emerges into the ambient unipolar field and then reconnection occurs between the unipolar field and the branch of the arch with the polarity opposite to the unipolar field. This reconnec- tion outflow leads to a fast-mode shock with the plasma flow upward along the open magnetic field lines and creates a jet. Observations of anemone jets made with Hinode (Tsuneta et al. 2008; Suematsu et al. 2008) belong to this class of standard jets, as explained by (Nishizuka et al. 2008, 2011). In chromo- spheric Hα observations, these regions of flux emerging can be observed as arch filament systems (AFSs; Bruzek 1967). An AFS is defined as a system of very low and dark fila- ments (Mandrini et al. 2002; Ma et al. 2015). These systems have a lifetime of several days; however, individual arches show a lifetime of a few tens of minutes. It has been observed that the hot flare loops overlie these AFSs at coronal heights (Alissandrakis et al. 1990; Malherbe et al. 1998).

X-ray jets are frequently accompanied by cool plasma jets called solar surges; those emerge as roughly straight threads of chromospheric dark material (Canfield et al. 1996; Schmieder et al. 1996; Wang et al. 1998; Chae et al. 1999; Mandrini et al. 2002; Uddin et al. 2012). Surges can be detected in hot coronal channels, such as 171 Å, 193 Å, and 211 Å, due to the fact that they block the extreme-ultraviolet (EUV) emis- sion coming from layers below them (Anzer & Heinzel 2005; Heinzel et al. 2008). Therefore, surges appear as dark areas alongside the hot jets in these AIA EUV channels (Nelson et al. 2019; Joshi et al. 2020a).

The other class of jets, the so-called blowout jets, are identi- fied by their width, which is as large in the ejection part as at the base (Moore et al. 2010, 2013). Blowout jets have been observed with two phases; the first phase is similar to the standard model, but is followed by a second phase characterized by the pres- ence of an elongated darkening in the EUV images, interpreted as a mini-filament (Sterling et al. 2015, 2022; Panesar et al. 2016, 2022; Kumar et al. 2019). These rising filaments are assumed to bring free energy due to shear or twist and recon- nect with the ambient magnetic field, thereby initiating blowout jets (Wyper et al. 2019; Sterling et al. 2022; Kumar et al. 2023). Many jets have helical motions and their twist can be trans- ported from a nearby flux rope to the jet (Joshi et al. 2020b, 2021). The role of twist in unfurling the blowout jet mate- rial has also been revealed by several numerical experiments (Pariat et al. 2009, 2010, 2015; Moreno-Insertis & Galsgaard 2013; Archontis & Hood 2013; Fang et al. 2014; Wyper et al. 2017).

Using multi-instrument observations, several studies on jets have been carried out to explain the twist in mini-filaments at the bases of blowout jets (Cheung et al. 2015; Tiwari et al. 2018), bilateral flows at the reconnection points (Ruan et al. 2019), reconnection in null points or in bald patch regions (Joshi et al. 2020b), and the formation of plasmoids (blob-like plasma ejections or bright kernels) (Singh et al. 2012; Zhang & Ji 2014; Rouppe van der Voort et al. 2017; Bogdanova et al. 2018;

Chen et al. 2022). These bright kernels trace the bilateral flows forming at the junction of the arch and the spire of the jet and commonly have a higher density and temperature than the jet. Numerical experiments of the jet process propose the tearing mode instability to form the plasma islands associated with the bright kernels (Guo et al. 2013; Ni et al. 2016; Li et al. 2023).

In the present study, we use observations of high spatial and temporal resolution made with the Swedish 1-m Solar Telescope (SST; Scharmer et al. 2003), together with observations from the SDO satellite, to analyze a sequence of jets originating from an active region in which emerging magnetic flux (EMF) is continu- ally appearing. The moving magnetic patches associated with the EMF episodes form a favorable site for jet launching in the close vicinity of the leading sunspot. Small B-class flares, with one of them reaching the B8.8-class as recorded by the GOES satellite, were related to recurrent ejections of plasma blobs and jets from these emerging regions. In Sect. 2, we present an overview of the jet observations. In Sect. 3, we present the recurrent jet envi- ronment and the jet dynamics. In Sect. 4, we present the NFFF extrapolation results. In Sect. 5, we discuss the trigger and occur- rence of the jets according to the magnetic configuration and the drifting of the null-point location.

## 2. Observations

### 2.1. Ground-based observations

Observations of the AR NOAA 13038 (N16 W61) were taken during the coordinated observational campaign involving the SST and SDO in June 2022. The observations were recorded using two instruments at the SST: the CRisp Imaging Spec- troPolarimeter (CRISP; Scharmer et al. 2008) and the CHRO- Mospheric Imaging Spectrometer (CHROMIS; Scharmer 2017). CRISP and CHROMIS are Fabry-Pérot tunable filter instru- ments equipped with dual etalon systems (Scharmer 2006). The CRISP instrument has the capability to measure the linear and circular polarization in photospheric and chromospheric spectral lines. With CRISP, the observations were taken in Hα (31 line positions from −1500 to +1500 mÅ), Ca II 8542 Å (20 line posi- tions from −1680 to +2380 mÅ), and Fe I 6173 Å (14 line posi- tions from −320 to +680 mÅ) at a cadence of 40 s. The Ca II and Fe I observations were recorded in spectropolarimetric mode. The CHROMIS data were taken in Hβ 4861 Å (27 line positions between ±2100 mÅ) at a cadence of 7 s. The spatial sampling for CRISP and CHROMIS is 0.′058 per pixel and 0.′038 per pixel, respectively. This SST dataset was reduced using the SSTRED pipeline (Löfdahl et al. 2021; de la Cruz Rodríguez et al. 2015). High image quality was achieved with the aid of the SST adap- tive optics system (Scharmer et al. 2024) and image restora- tion with the Multi-Object Multi-Frame Blind Deconvolution (MOMFBD; Van Noort et al. 2005) method. The spectropolar- metric Fe I 6173 Å observations were inverted using a Milne- Eddington inversion code (de la Cruz Rodríguez 2019) to pro- duce maps of the magnetic field strength along the line of sight ($B_{LOS}$).

### 2.2. Space-borne observations

For the EUV and magnetic field observations, we used Atmo- sphere Imaging Assembly (AIA; Lemen et al. 2012) and Helio- seismic Magnetic Imager (HMI; Scherrer et al. 2012) respec- tively, on board the Solar Dynamics Observatory (SDO; Pesnell et al. 2012). AIA observes the full solar disk, which is
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** General overview of AR 13038, where the location of the AR is marked with a white rectangular box in the full disk AIA 171 Å image (panel a) and the focused FOV is shown in panels (b) and (c). The white box in panel (b-c) is the FOV for the middle and last row, showing the maximum FOV, which was covered during the coordinated observations by SST. The straight enhancement of brightening in 171 Å (yellow contours of AIA 171 in panel c) is along with the polarity inversion line between the main positive polarity and the emerged negative polarity (N1, N2, N3 in Fig. 2). The white and black contours label the magnetic field strength of ±500 G. The middle row presents the AIA observations aligned to SST observations at a specific time and the AIA spatial resolution is blown up to the SST spatial resolution. An animation is attached online with this figure to show the evolution of the jet region from 09:47 to 11:30 UT.</figcaption></figure>

ideal for studying the large-scale jet evolution in coronal diagnostics. The AR was observed with all EUV (94, 131, 171, 193, 211, 304, and 335 Å) and UV (1600 and 1700 Å) channels of AIA with cadences of 12 sec and 24 s, respectively. The pixel size of AIA is 0.76. HMI observes the photospheric magnetic field with a cadence of 45 s and a pixel size of 0.5. We used these magnetograms for the AR evolution, starting one day prior to the beginning of the jet activity. In addition, we also used the Spaceweather HMI Active Region Patch (SHARP;

Bobra et al. 2014) data to make the magnetic field extrapolations. This SHARP dataset is available for each AR in the solar disk with a cadence of 12 min. To compare ground-based observations with SDO observations, we aligned the AIA and HMI datasets with the SST observations by matching the AIA and HMI spatial sampling to that of the SST. This alignment procedure was carried out using SSWIDL<fnref n="1" />.

<footnote n="1">¹ https://robrutten.nl/Recipes_IDL.html</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Evolution of the emerging flux on the east side of the leading positive sunspot over a duration of 12 h. This sunspot acts as a source of flux emergence, and emerged polarity patches move to form the regions marked by the two ovals in panel (c). The negative polarity on the west side of the sunspot is an artifact due to the projection effect. The threshold of the magnetic field is ±500 G.</figcaption></figure>

**3. Recurrent jet environment**

On 24 June 2022, AR 13038 was the most highly jet-producing region in the Sun with several GOES B and C class solar flares. We focus on the jets and their environment occurring during the SST observations from 09:47 UT to 10:55 UT.

**3.1. Active region evolution and magnetic flux emergence**

Figure 1 presents a global picture of AR 13038 located near the west limb (N 15, W60). A part of this AR was observed with SST, which is shown in the last row of Fig. 1, with a rotated white square. The basic magnetic configuration of this region consists of a leading sunspot of positive polarity and an extended region of negative polarity trailing on the east side. The negative polarity patch on the upper west side (or limb side) of the leading positive polarity sunspot is due to the projection effect from its close proximity to the limb.

Long loops are observed in coronal temperatures (AIA 171 Å), joining the two main polarities in the east-west direction in the southern end of the sunspot. Towards the southern direction of the AR, nearly open field lines (long open loops) are observed. On the east side of the sunspot, an extensive region with EMF is found (Fig. 2). The large positive polarity is surrounded by small negative polarities (see the cyan oval), which were formed due to flux emergence over time. Looking at the evolution of magnetic field over the duration of 12 h, from June 23 22:00 UT to June 24 11:00 UT, the small negative polarities are moving to the northeast, indicating a strong shear between the negative and positive magnetic environment. The motion of the small negative polarities stops when reaching the former positive magnetic network located between the two main polarities. The zoomed-in image of this area with SST/CRISP (Fig. 3 panels c and f) shows more detail of the small negative-polarity patches.

Two systems can be identified with labels N1, N2, and N3 along with P1, P2, and P3 encircled in the cyan oval, and N4 and N5 along with P4 and P5 in the yellow oval. The latter small-polarity patches are not detectable in the HMI magnetograms. As a consequence of the continuous magnetic flux emergence for an extended time, the motion of small new magnetic polarities has been observed in the longitudinal magnetic

field maps, between the negative polarities (N1, N2, and N3) and the positive polarities (P1, P2, and P3). Due to the location of the AR towards the limb, the longitudinal magnetic field measurements have a strong horizontal component. Therefore, horizontal fibril-like structures between small positive and negative polarities are visible (an example is shown with a white arrow in panel c of Fig. 3). Their direction is inclined by an angle of about ~120° with respect to the PIL, indicating a strong shear.

**3.2. Overview of the jets**

We observed the AR for a duration of ~4 hours (from 08:00 to 12:00 UT) with SST. The observations by CRISP in $H\alpha$ focus on the EMF region in two successive time sequences. A GOES C1.5 class flare started at 08:49 UT, peaked at 08:58 UT, and decayed at 09:13 UT, and was followed by a jet towards the south at 09:17 UT and a jet towards the east at 09:32 UT. Soon after, at 09:59 UT, a GOES B8.8 class flare started to occur, peaked at 10:20 UT, and decayed at 10:30 UT in the GOES record. In the present paper, we focus on the later time series observations within 09:47–10:55 UT. The continuous ejection of jets at 09:53 UT, 10:12 UT, and 10:57 UT was observed to be associated with a small amount of flaring activity. Emergence of magnetic flux was a continuous process for many hours before the activity. New flux emergence creates favorable conditions for surges, jets, and flares to occur during this long process and we observed jets, surges, and bright kernels issued from a flare ribbon of ellipsoidal shape (encircled in Fig. 4 panels e-f) around the flux-emergence region at multiple times (see the attached movie associated with Fig. 1).

In Fig. 4 we present the jets observed in AIA 171 Å, 304 Å, and $H\alpha$ (Jet 1, Jet2, Jet3). Jet 1 and Jet 2 are going to the east, while Jet 3 is going to the south. After Jet 1 at 09:47 UT, which is apparent in the EUV filters as well as in the blue wing of $H\alpha$, we observed the formation of an AFS in the $H\alpha$ line core (see movie attached with Fig. 1). We believe that, during the flux-emergence process, dense material is lifted up, forming an AFS (as explained in Schmieder & Pariat 2010; Ma et al. 2015). Comparing the images in 171 Å, 304 Å, and $H\alpha$ (Fig. 4), we observe that the dark area surrounded by brightenings corre-
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Magnetic field configuration resulting from the emerged magnetic flux. The left column is observed with CRISP and right column is observed with HMI. Two systems of negative and positive polarities can be identified: in the north, N1, N2, and N3 with P1, P2, and P3, and in the south N4 and N5 with P4 and P5. We can see the weak negative polarities N4 and N5 in the CRISP observation (panel a). These are moving continuously and the HMI spatial resolution is not sufficient to resolve them. The right column zooms onto the blue rectangular area shown in the middle row. The white arrow in panel (c) shows the magnetic polarities moving between the negative (N3) and positive (P3-P4) boundaries observed in the longitudinal magnetic field maps. An animation is attached online with this figure to show the magnetic field evolution from 09:47 to 11:30 UT. The threshold of the magnetic field is ±500 G.</figcaption></figure>

sponds to the AFS region between the two flare ribbons in $H\alpha$. The AFS is well visible in the $H\alpha$ images (see panel f) as dark strands at the top of the ellipsoid covering the flare region. In $H\alpha$, we observe what seems to be a progressively peeling-like mechanism of the AFS followed by the ejection of elongated surges (see Fig. 5). In SST observations, we detected numerous stripes of dark material peeling away from the AFS just before the jet eruption. As the jet begins, fewer stripes are visible. From the animation associated with Fig. 5, we can follow the peeled material, which is moving towards the south (see the blue arrow in Fig. 5 last column) and going towards the east in the jet. This resembles the peeling mechanism of the external layers of the emerged domes, which is frequently reported in flux-emergence experiments, and precedes the surge ejection (see, e.g., the animation of Figure 10 by Moreno-Insertis & Galsgaard 2013). In the AIA filters, the dense and cool material absorbs a part of the EUV emission and appears as an unresolved dark structure (see panel (a) of Fig. 5), which could have been interpreted as a mini filament eruption, in the absence of high-resolution SST observations.

There are many blobs being ejected alongside the Jet2 in AIA 171 and 304 Å (shown with the black arrows in panel d of Fig. 4). Another arcade occurs for a second time after Jet 2

between 10:22 UT and 11:00 UT covering the former AFS with a different direction (NW to SE). These arcades could correspond to post flare loops after reconnection and cooling hotter loops visible in AIA 171 Å. Flare-ribbon brightenings are apparent in Fig. 4 in the negative and positive magnetic field. Around 10:18 UT (Jet 2), a set of east-west bright loops are observed in AIA 171 Å (panel d). The west footpoints of the bright loops are identified as a north-south (NS) elongated feature. This NS feature is co-spatial with the positive polarities P3, P4, and P5 included in the yellow oval with N4 and N5 in panel a of Fig. 2. Jet 3 appears as a fine shape in 171 Å and 304 Å (Fig. 4, panels g and h) and corresponds to a broad dark jet in $H\alpha$ (panel i).

High-temperature brightenings are detected in the hot AIA channels (the best is seen in the AIA 94 Å channel in Fig. 6). The light curves obtained in two selected boxes (marked in Fig. 6 first row) indicate the normalised intensity at different times (Fig. 6 bottom row). We have normalized the total intensity with the maximum value of the intensity of the time series. These two boxes are of ~10" in size. The intensity variation at the different box locations in six EUV channels of AIA are plotted in different colors. Box1 in the north of the EMF reaches maximum intensity between 10:20 and 10:25 UT, while box2 in the south reaches its
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** EUV (171 and 304 Å) and Hα observation of jets at their maximum height. Magnetic field contours of strength ±500 G from SST/CRISP Fe I 6173 Å are over-plotted in the second and third columns for all three jets. Jet 1, Jet 2, and Jet 3 are shown with white arrows and the ellipsoidal shaped flare ribbons around the AFS are shown with black dashed contour in panels (e–f) and the ejected blobs of Jet 2 are indicated by black arrows in panel (d). The bright loops (west to east) heated by the flare are shown in panel (d). The white rectangular box in panel (b) shows the FOV for Fig. 5.</figcaption></figure>

maximum between 10:30 and 10:40 UT. These brightenings are associated with the jets towards the east at 09:53 and 10:18 UT and towards the south at 10:58 UT; they only show a clear intensity peak in the hot AIA 94 Å channel. We believe that the plasma in the brightening reaches such a high temperature that it is more clearly visible in the hot AIA 94 Å wavelength. The brightening bounded by box1 in the north is observed to move towards the south (see box2) in ~18 minutes. This is a clear indication that the reconnection site is moving from north to south and the overlying magnetic field lines have different directions to the east or south, which explains the change of direction of the

jets. The distance between the two locations of the brightening in Fig. 6c and e is estimated to be 20 Mm.

### 3.3. Doppler shifts in Hα

On the basis of the spectroscopic observations from CRISP, one can calculate Doppler shifts using the Hα lines. We constructed Dopplergrams by subtracting two different Hα line positions at ±45 km s⁻¹ offset from line center. Figure 7 shows the resulting color maps for the vertical velocity component, with the downflows drawn in red and the upflows in blue, at different times
<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Zoom onto the AFS region, corresponding to the field of view displayed in Fig. 4 (panel b). The left column displays observations in the AIA 304 Å wavelength, the middle column shows SST observations in the Hα blue wing, and the right column is the same as the middle but in pale rendering. The right column was used to manually trace the peeled stripes (marked in red) from the AFS. The blue arrow indicates the direction of the peeling material flowing southward and contributing to the jet (highlighted in magenta color). The feature observed in AIA 304 Å appears as an unresolved dark structure, which is resolved in SST observations as an AFS. This figure is accompanied by an online animation.</figcaption></figure>

from 9:49:48 UT (top left panel) to 11:07:03 UT (bottom right panel). Negative and positive magnetic field contours (±500 G) are plotted in the second column (for all three jets) in black and green, respectively, using the Fe I 6173 Å magnetograms from SST. Here, upflowing jet material is visible between the negative polarities and the positive polarities $x = 760''-770''$, $y = 230''-240''$ (see panels b,e, and h). The maps clearly show the jet flows as elongated blue features (and blobs in the case of Jet 2); interestingly, one can also identify a perturbation that can be interpreted as a moving flux rope evolving from north to south. The initial location of this perturbation is shown with a black arrow in panel c; the motion can be discerned in the associated movie

as blueshifted plasma starting in between the N1, N2, and N3 and P1, P2, and P3 region and moving toward the south, with subsequent ejection of plasma material as Jet 2 at ~10:05 UT. We show that this motion can be interpreted as a flux rope (with a right-handed twist) in the following section (Sect. 3.4) with the help of Hβ observations.

The motion of the flux rope also causes collimated downflows oriented towards the south (panel e, f). For the southern region from P4 and P5, from 09:49 UT until 10:12 UT, there are mainly downflows (redshifted material) around $x = 770''-780''$ and $y = 220''-240''$. Then from 10:57 UT (Jet 3) there are upflows (blueshifted material) and Jet3 continues till 11:07 UT
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Drifting of the brightening from north (box1) to south (box2). Top row: AIA 94 Å middle row: AIA 171 Å bottom row: Temporal evolution of the normalized intensity in six AIA channels in box1 and box2. Three jets are marked with white arrows in panels (g), (i), and (l), respectively. Drifting of the hot material is quantified in the hot AIA 94 Å channel with the normalized intensity curves (blue curve in the bottom row).</figcaption></figure>

with up- and downflows. The base of Jet3 is approximately close to the negative-polarity N3 plotted in black contours. The downflows of Jet 3 are observed in many parallel structures covering a large area. The footpoints of these structures extend ~15" towards the south in the region of P5 (green contours in the last row).

### 3.4. Swirling flux rope in Hβ observations

Inspecting the high cadence (7 s) CHROMIS Hβ 4861 Å observations in different spectral line positions, similar to Hα Dopplermaps, we find evidence of a cool plasma cloud moving from the north of the EMF region (above the negative polarity N1) and ejecting the plasma material as Jet 2 towards the east from the location of negative polarities N2 and N3 (see Fig. 8 and the associated animation); we believe that this twisting material corresponds to the flux rope swirling over the AFS region and flowing away towards the east, as discussed in Sect. 3.3. To test the twisting motion of the flux rope, we performed spectral analysis of the Hβ profiles in the AFS region presented in Fig. 8. We chose two time steps at 10:00 UT and 10:14 UT, shown in the upper (panels a–c) and lower half (panels h–j) of Fig. 8, respectively; we selected two different points in each (shown with plus and asterisk symbols in red in the respective Hβ −0.4 Å panels) to calculate the spectral evolution by means of λ − t diagrams. The spectral profiles at the flux-rope location are plotted for various time steps. These specific time steps are indicated by blue, red, and orange dashes on the right axis of the spectra (see panels d, f, k, and m). The spectral profiles for different

locations within the flux rope, marked by red asterisks and plus symbols, are shown in panels e, g, l, and n. Exploring the spectra and profiles at different times, it is interesting to follow the dynamics of the plasma at the flux rope location as we see clear asymmetries and even separate components in both wings. For example, the red and orange profiles in panel e have a separated component at a Doppler shift of about −50 and −60 km s⁻¹. The blue profile in this panel e has a redshifted component at about +40 km s⁻¹. Panel g also shows examples from this region with similar Doppler-shifted components. The associated animation shows the temporal evolution of the AFS region in the Hβ line core (left panel), blue wing (−0.4 Å, central panel), and Doppler maps (±0.4 Å, right panel) in the interval 09:47 UT to 10:33 UT. From 09:55 UT, we see highly dynamical small dark threads that appear to move southward and some of them in a swirling pattern above the AFS region. These dynamical threads or stripes are shown in Fig. 5 in red color, explained as they are peeled out from the AFS. The swirling pattern of a large dark plasma cloud is particularly clear in the blue wing and Doppler movie and we interpret this as a swirling flux rope (similarly seen in Hα Doppler map observations). This swirling flux rope travels southward over this region and evolves into the ejection of Jet 2 towards the east at 10:06 UT.

We conjecture that these rapidly changing blueshifts and redshifts in the spectra correspond to the twisted flux rope passing from north over the AFS and contributing to the Jet 2 material. At 10:00 UT, the spectral profiles at different times at two locations (see the asterisk and plus symbol in panel b) show a blueshift and a redshift (panels e and g), where the spectral profiles have the
<figure><img src="image_8.png" /><figcaption>**Fig. 7.** SST/CRISP Hα Doppler maps showing the evolution of the three jets. The location of the dynamic swirling flux rope is shown in panel (c) (see the attached animation). These Doppler maps are constructed using two Hα line positions offset by 1 Å from the line center (±45 km s⁻¹). Magnetic field contours of ±500 G (black and green) from the SST Fe I 6173 Å magnetograms are overlaid in the first column. This figure is associated with an [online animation](#).</figcaption></figure>

red and blue asymmetry in comparison to the average profile at a quiet region (gray color). At the time of Jet 2 at 10:14 UT, the profiles show a clear blueshift of more than −60 km s⁻¹ (see the blue curve in panel l corresponding to *t* = 27 min in panel k). Similar blueshifts are also evident in other profiles from the jet material (location of the asterisk in panel i and spectra in panel m and n). In this way, these very highly resolved observations with a higher cadence (as the Hα observations has a cadence of 40 s in comparison to Hβ with 7 s) are showing the highly dynamical swirling motion of the flux rope over the AFS with alternate blueshifts and redshifts, and only blueshifts in the jet.

# 4. Non-force-free magnetic field extrapolation

## 4.1. Method

We apply the NFFF extrapolation technique, as developed by Hu & Dasgupta (2008) and Hu et al. (2010), to extrapolate the coronal magnetic field of AR 13038. The time-frames under investigation are from 10:00 to 11:12 UT on June 24, 2022, which include the time of the vector magnetograms portrayed in Figs. 9 and 10.

The NFFF extrapolation technique is built on the minimum energy dissipation rate principle (MDR), which assumes that the plasma system tends towards a state with the lowest possible
<figure><img src="image_9.png" /><figcaption>**Fig. 8.** CHROMIS Hβ observations of the swirling flux rope (panels a–g) and Jet 2 (panels h–n). For two time steps (10:00 UT and 10:14 UT), the row of images (a–c, and h–j) shows the Hβ line core, the blue wing, and a Dopplergram, respectively. In the Hβ wing images, a red plus symbol and an asterisk mark the locations for which the spectral evolution in λt diagrams is shown in panels (d), (f), (k), and (m). The spectral profiles for these locations for the two chosen time steps are shown in red in the profile panels (e), (g), (l), and (n). For each λt diagram, two more spectral profiles from different times are shown in orange and blue. The times for these profiles are indicated with short colored dashes in the λt diagrams. The thin gray profile is a reference Hβ profile, which is the average over a quiet region. The time along the t axis in the λt diagrams is in minutes after the start of the observation at 09:47:08 UT. An animation of the three image panels is provided in the online material.</figcaption></figure>

dissipation rate on a timescale shorter than the dynamical one. In that case, the magnetic field, denoted as *B*, fulfils the double-curl Beltrami equation (Bhattacharyya et al. 2007). The field *B* is consequently described as the sum of two linear-force-free fields and a potential field (Hu et al. 2010). This results in a magnetic field configuration that can have nonzero Lorentz force. An iterative method is then applied to the NFFF extrapolation that minimizes the average discrepancy between the observed and calculated horizontal field at the photospheric boundary. The

NFFF technique has recently been effectively used to interpret various transient phenomena in active regions, such as flares, coronal jets, and coronal dimmings (Prasad et al. 2018; Nayak et al. 2019; Prasad et al. 2020, 2023). We recognize the unreliability of vector magnetograms when the region of interest is too far from the disk center, as explained in Bobra et al. (2014) and Falconer et al. (2016). In our case, the region is at 58°, and so to miminize the problem, we employed SHARP data remapped from CCD coordinates to a heliographic cylindrical
<figure><img src="image_10.png" /></figure>

**Fig. 9.** Non-force free field extrapolation of AR 13038 (top view) at 10:12 UT. This rendering shows selected magnetic field lines with magenta field lines connecting to the region containing the N1, N2 and N3 with P1, P2, and P3 polarities. The yellow extrapolated field lines show the field-line geometry in the southern part of the AR and are co-spatial with the ejection of Jet 3.

equal-area (CEA) projection centered on the selected patch for the NFFF extrapolation. We also visually compared the SHARP data presented in Fig. 9 with the SST magnetograms and did not find any projection artifacts, such as the fictitious negative polarity on the limb side in the $B_{LOS}$ maps in Fig. 2. The magnetogram shown in panel a of Fig. 10 spans 640 × 320 pixels in the x and y directions of a Cartesian coordinate system, while the vertical extent of the domain is set to be 320 pixels. In physical terms, the dimensions of the computational domain are approximately 466 Mm × 233 Mm × 233 Mm in the x, y, and z directions, respectively.

## 4.2. Magnetic topology

The NFFF extrapolation using HMI vector magnetic field maps shows the existence of a null point over the region of negative polarities N1, N2, and N3 connecting with the positive polarities on the east side (Fig. 10). The second system of negative polarities (N4 and N5) is not very prominent in the HMI resolution, and so it is not possible to look for a second null point in that location using extrapolations from the HMI data. For extra guidance, we also calculated maps of the squashing factor $Q$ (Priest & Démoulin 1995; Demoulin et al. 1996; Titov & Hornig 2002; Liu et al. 2016) on the photospheric plane. Narrow regions with very large $Q$ values (plotted with colored contours in panels b and c of Fig. 10) indicate a strong gradient in the remote connectivity of the field lines issuing from them; that is, they mark clear candidates for the photospheric base of a separatrix or quasi-separatrix layer (QSL).

The result of the extrapolation is particularly interesting for the field lines starting in the region around the N1, N2, and N3 negative polarities (Fig. 10). Those negative magnetic concentrations act as parasitic polarities surrounded by positive patches (see, e.g., Fig. 10, panels b and c or Fig. 9); consequently, the presence of a null point above them is a natural feature of the magnetic configuration. In the figure, the isocontours for large values of the squashing factor $Q$, drawn on top of the photospheric magnetogram, show a simple closed-circuit shape; the field lines that leave the photosphere at the site of the widest

<figure><img src="image_11.png" /></figure>

**Fig. 10.** Different views of the NFFF extrapolation. Panel (a): side view. Panel (b): top view. Panel (c): zoomed-in FOV of the top view. The three negative polarities N1, N2, and N3 are connected to the positive polarities P1, P2, and P3 nearby with small inner spines, and globally to the plage region (large negative polarity region to the east). The squashing factor $Q$ is plotted at the photospheric level in panels (b) and (c) (colored contours) to show the separatrix at the periphery of the fan surface around N1, N2, and N3.

high-$Q$ contour are seen to delineate a fan surface converging toward a location at around $z = 6$ Mm, where the null point is expected; moreover, the field intensity at that point is seen to be a minimum compared to the surroundings. This is a standard fan-spine null configuration; the upper spine and neighboring field lines arch over the active region and link to the large negative polarity to the east of the image. These field lines provide a path for Jet 1 and Jet 2 towards the east.

Figure 11 shows a zoomed-in image of the fan-spine structure. The region harboring the null point is marked with a green arrow and a star; the figure also indicates the approximate location of the inner spine (linking the null point to the parasitic negative polarities below it) and the outer spine that joins the null point to the large negative polarity to the east of the active region. The SST observations show that flux emergence is continually taking place in the region underneath the fan surface; moving magnetic patches are detected in the space between the opposite polarities N1-N2-N3 and P1-P2-P3 (see Section 3.1 and the accompanying animation to Fig. 3). All of this is indicative of
<figure><img src="image_12.png" /></figure>

**Fig. 11.** Side view of the fan-spine topology with a clear view of the location of the null point. The null-point location at the top of the fan surface (the site of the widest high-Q contour indicated with an elongated red surface) at a height of 6 Mm above the photosphere is marked with a green arrow and star. The approximate location of the fan surface, the inner spine (linking the null point to the parasitic negative polarities), and the outer spine (moving away from the null point) are shown with black arrows.

reorganization of the magnetic structure in the low atmosphere underneath the fan surface, which is likely to lead to the production of jets 1 and 2.

Understanding the origin of jet 3 is more complicated. The weak visibility in the HMI vector magnetic maps of the magnetic field polarities to the south of the region analyzed in the previous paragraph prevents the identification of a second null point (or region with large connectivity gradient) in that area. On the other hand, the SST FOV contains that region and can resolve the negative polarities N4 and N5, but does not cover the full active region (see Fig. 2); and so the SST data do not provide a flux-balanced magnetic map on which an extrapolation could be based. Nevertheless, it is highly plausible that, after the launching of jets 1 and 2, there is a drift of the reconnection site toward the south, which is also supported by the moving brightening observed in the hot AIA 94 Å channel. Further, going back to Fig. 9, the long yellow extrapolated field lines were traced from a region to the south of (and near to) the fan-spine configuration discussed above; we note that they first have a southward-pointing orientation that roughly follows the open-loop path in AIA 171 Å, and are co-spatial with the path for the Jet 3. This could explain the abrupt change of orientation of Jet 3 compared with that of the two previous jets.

**5. Discussion and conclusions**

We present three recurrent solar jets observed in AR NOAA 13038 on June 24, 2022, using SDO and SST observations. The novelty of this work lies in the coordination of SDO observations with SST. It is exceptional that we were able to resolve the fine-scale chromospheric structures in the jet environment using observations of high spatial and temporal resolution. The main observational findings of the paper are as follows:

1. Due to the continuous emergence of small parasitic magnetic fluxes at the edge of the main positive sunspot in the active region over the course of a day, the negative polar-

ties are surrounded by positive polarities, as observed with CRISP/SST. The presence of these parasitic polarities leads to a magnetic configuration featuring a dome, a fan, and a null point. Due to the continuous emergence of new flux at the edge of the main positive sunspot in the active region over the course of a day, a collection of negative polarities appear next to the sunspot surrounded by positive flux patches, as observed with CRISP/SST.

2. The EUV images from AIA (171 Å, 304 Å) show a dark area surrounded by an oval of brightenings at this location. This bright oval corresponds to the base of the dome. The dark area might be interpreted as a mini filament; however, SST images in $H\alpha$ and $H\beta$ reveal an AFS instead of a filament. The filaments in this system are fine fibrils oriented from the positive polarities (P1, P2, P3) to the negative polarities (N1, N2, N3), similar to the positive threads seen in CRISP magnetic maps. They form loops or arcades with downflows at both ends, as illustrated in a sketch of emerging flux by Schmieder et al. (2014). These fibrils are not perpendicular to the PIL but exhibit a shear, which introduces free energy.

3. The Doppler shifts obtained using $H\alpha$ (CRISP) and $H\beta$ (CHROMIS) spectra clearly show that a large flux rope, with width comparable to the width of the dark area, is swirling over the AFS from north to south. This large twisted flux rope forms because of the shear during the flux-emergence phase.

4. The oval-shaped brightening observed in EUV lines corresponds to the flare ribbon seen in $H\alpha$ observations. Along this ribbon, several bright points in EUV lines are observed, which could be due to canceling flux. However, these points do not correspond to the bases of the jets.

5. In $H\alpha$, following the swirling flux rope, another system of AFS developed from north to south, but with a different orientation. These could either be cooling post-flare loops or another AFS. However, if reconnection were to occur at the null point, we would be likely to observe post-flare loops.

Magnetic flux emergence is crucial for injecting magnetic energy into a system and for stressing the pre-existing magnetic field through shearing. Magnetic reconnection can occur continuously and at different heights within the solar atmosphere. As a result of reconnection, flux cancelation can occur between the emerging parasitic polarities and the pre-existing region of opposite polarity. In our current observational case, we conjecture that the jets occurred due to long-term continuous magnetic flux emergence, which induced shear between the pre-existing magnetic field and the newly emerging magnetic patches. In such scenarios, magnetic flux emergence might inject Poynting flux into the system and drive magnetic flux cancelation and the ejection of jets. Initially, the continuous motion of negative and positive polarities created a sheared region, leading to the formation of a twisted flux rope. Two jets were ejected to the east, and later another jet (Jet3) was observed being ejected towards the south. Similarly, we also observed the drifting of a brightening at the jet base towards the south in the hot EUV channel, which could correspond to a moving reconnection site.

For this small and complex region, the NFFF extrapolation confirms the existence of a null point, corresponding to the jet flow towards the east. From the extrapolation, the height of the null point is around $z = 6$ Mm, which is consistent with the ranges given in the literature (Galsgaard et al. 2017; Joshi et al. 2020a). In this present observation, it is not possible to accurately establish the existence of a second null point in the south,
as the HMI magnetogram only shows weak polarities there, and
the area is outside the SST FOV.

The change of direction of the jets from east to south
could be a result of the displacement of the reconnection site
to the southern region where the initial sections of the mag-
netic loops are oriented north-south, as seen in projection in
the observation. From these EUV observations, we estimated
a possible drift of the reconnection site (brightening in the
hot EUV channels) of 20 Mm in 18 min, implying an aver-
age velocity of 18.5 km s⁻¹. This is larger than the velocities
inferred for null-point displacements related to coronal bright
points from observations (Galsgaard et al. 2017) and simulations
(Nóbrega-Siverio & Moreno-Insertis 2022). However, it is note-
worthy that the magnetic field topology analyzed in the present
paper is more complex. The direction of the jets would depend
on the site of the magnetic reconnection. The earlier reconne-
tion events (Jet1 and Jet2) potentially altered the position of the
reconnection site and the orientation of the magnetic field lines
passing through the reconnection region. Such a reconnection-
driven shifting can easily happen if the spine (in the case of a null
point) or the high-Q flux tube (in the case of a QSL) are located
within a region where all field lines strongly diverge from the
reconnection region. With such geometry, it is quite simple to
shift the spine or the high-Q tube to nearly any position within
the “fanning-out” diverging flux system.

Using the SST observations, we identified the presence of an AFS at the base of the jets, typically representing the chromospheric manifestations of emerging flux regions. This AFS appears as dark structures in EUV lines and can be identified as mini-filaments in the absence of high-resolution observations. In this regard, we look forward to high-resolution observations from both existing and upcoming ground-based instruments, like DKIST (Rimmele et al. 2020) and EST (Quintero Noda et al. 2022). Given the results of the extrapolation, and the changing direction of the jets, we conclude that, during the evolution of the jets, the reconnection site could have moved from one location to another. In general, during jet events or in confined flare events, the null point (reconnection site) may actually be able to move and eventually involve completely different magnetic loops in a new magnetic reconnection event (Masson et al. 2009) or sometimes may show a persistent nature (Cheng et al. 2023). This can also be interpreted as “interchange reconnection”, which gives rise to two reconnection products and sets a new footpoint connection for the coronal magnetic field lines (Sterling et al. 2015). Future high-resolution ground-based observations of the photosphere coordinated with space instruments for coronal observations (e.g., MUSE, De Pontieu et al. 2022; Cheung et al. 2022) will provide a clearer picture of the moving reconnection site scenario for solar jets.

## Data availability

Movies associated to Figs 1, 3, 5, 7, 8 are available at
https://www.aanda.org

Acknowledgements. We thank to the referee for the comments which help to improve the manuscript. This research is supported by the Research Council of Norway, project number 325491, through its Centres of Excellence scheme, project number 262622, and the European Research Council through the Synergy Grant number 810218 (“The Whole Sun”, ERC-2018-SyG). This work benefited from discussions at the Instituto de Astrofísica de Canarias under the “IAC early career visitor program” and at the International Space Science Institute (ISSI) in Bern, through ISSI International Team project #535 “Unraveling surges: a joint perspective from numerical models, observations, and machine learning”. 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). SDO observations are courtesy of NASA/SDO and the AIA, EVE, and HMI science teams. We made use of NASA’s Astrophysics Data System Bibliographic Service.

## References

Alissandrakis, C. E., Tsiropoula, G., & Mein, P. 1990, A&A, 230, 200

Anzer, U., & Heinzel, P. 2005, ApJ, 622, 714

Archontis, V., & Hood, A. W. 2013, ApJ, 769, L21

Asai, A., Ishii, T. T., & Kurokawa, H. 2001, ApJ, 555, L65

Bhattacharyya, R., Janaki, M. S., Dasgupta, B., & Zank, G. P. 2007, Sol. Phys., 240, 63

Bobra, M. G., Sun, X., Hoeksema, J. T., et al. 2014, Sol. Phys., 289, 3549

Bogdanova, M., Zhelyazkov, I., Joshi, R., & Chandra, R. 2018, New Astron., 63, 75

Bruzek, A. 1967, Sol. Phys., 2, 451

Canfield, R. C., Reardon, K. P., Leka, K. D., et al. 1996, ApJ, 464, 1016

Chae, J., Qiu, J., Wang, H., & Goode, P. R. 1999, ApJ, 513, L75

Chandra, R., Gupta, G. R., Mulay, S., & Tripathi, D. 2015, MNRAS, 446, 3741

Chen, J., Erdélyi, R., Liu, J., et al. 2022, Front. Astron. Space Sci., 8, 238

Cheng, X., Priest, E. R., Li, H. T., et al. 2023, Nat. Commun., 14, 2107

Cheung, M. C. M., De Pontieu, B., Tarbell, T. D., et al. 2015, ApJ, 801, 83

Cheung, M. C. M., Martínez-Sykora, J., Testa, P., et al. 2022, ApJ, 926, 53

de la Cruz Rodríguez, J. 2019, A&A, 631, A153

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

Demoulin, P., Henoux, J. C., Priest, E. R., & Mandrini, C. H. 1996, A&A, 308, 643

De Pontieu, B., Testa, P., Martínez-Sykora, J., et al. 2022, ApJ, 926, 52

Falconer, D. A., Tiwari, S. K., Moore, R. L., & Khazanov, I. 2016, ApJ, 833, L31

Fang, F., Fan, Y., & McIntosh, S. W. 2014, ApJ, 789, L19

Galsgaard, K., Madjarska, M. S., Moreno-Insertis, F., Huang, Z., & Wiegelmann, T. 2017, A&A, 606, A46

Guo, Y., Démoulin, P., Schmieder, B., et al. 2013, A&A, 555, A19

Heinzel, P., Schmieder, B., Fárník, F., et al. 2008, ApJ, 686, 1383

Heyvaerts, J., Priest, E. R., & Rust, D. M. 1977, ApJ, 216, 123

Hu, Q., & Dasgupta, B. 2008, Sol. Phys., 247, 87

Hu, Q., Dasgupta, B., Derosa, M. L., Büchner, J., & Gary, G. A. 2010, JASTP, 72, 219

Joshi, R., Schmieder, B., Chandra, R., et al. 2017, Sol. Phys., 292, 152

Joshi, R., Chandra, R., Schmieder, B., et al. 2020a, A&A, 639, A22

Joshi, R., Schmieder, B., Aulanier, G., Bommier, V., & Chandra, R. 2020b, A&A, 642, A169

Joshi, R., Wang, Y., Chandra, R., et al. 2020c, ApJ, 901, 94

Joshi, R., Aulanier, G., Radcliffe, A., et al. 2024, A&A, 687, A172

Joshi, R., Schmieder, B., Tei, A., et al. 2021, A&A, 645, A80

Kumar, P., Karpen, J. T., Antiochos, S. K., Wyper, P. F., & DeVore, C. R. 2019, ApJ, 885, L15

Kumar, P., Karpen, J. T., Antiochos, S. K., et al. 2023, ApJ, 943, 156

Lemen, J. R., Title, A. M., Akin, D. J., et al. 2012, Sol. Phys., 275, 17

Li, X., Keppens, R., & Zhou, Y. 2023, ApJ, 947, L17

Liu, J., Wang, Y., Shen, C., et al. 2015, ApJ, 813, 115

Liu, J., Wang, Y., Erdélyi, R., et al. 2016, ApJ, 833, 150

Liu, Y., Ruan, G. P., Schmieder, B., et al. 2022, A&A, 667, A24

Löfdahl, M. G., Hillberg, T., de la Cruz Rodríguez, J., et al. 2021, A&A, 653, A68

Ma, L., Zhou, W., Zhou, G., & Zhang, J. 2015, A&A, 583, A110

Malherbe, J. M., Schmieder, B., Mein, P., et al. 1998, Sol. Phys., 180, 265

Mandrini, C. H., Démoulin, P., Schmieder, B., Deng, Y. Y., & Rudawy, P. 2002, A&A, 391, 317

Masson, S., Pariat, E., Aulanier, G., & Schrijver, C. J. 2009, ApJ, 700, 559

Moore, R. L., Cirtain, J. W., Sterling, A. C., & Falconer, D. A. 2010, ApJ, 720, 757

Moore, R. L., Sterling, A. C., Falconer, D. A., & Robe, D. 2013, ApJ, 769, 134

Moreno-Insertis, F., & Galsgaard, K. 2013, ApJ, 771, 20

Moreno-Insertis, F., Galsgaard, K., & Ugarte-Urra, I. 2008, ApJ, 673, L211

Nayak, S. S., Bhattacharyya, R., Prasad, A., et al. 2019, ApJ, 875, 10

Nelson, C. J., Freij, N., Bennett, S., Erdélyi, R., & Mathioudakis, M. 2019, ApJ, 883, 115

Ni, L., Lin, J., Roussev, I. I., & Schmieder, B. 2016, ApJ, 832, 195

Nishizuka, N., Shimizu, M., Nakamura, T., et al. 2008, ApJ, 683, L83

Nishizuka, N., Nakamura, T., Kawate, T., Singh, K. A. P., & Shibata, K. 2011, ApJ, 731, 43

Nisticò, G., Bothmer, V., Patsourakos, S., & Zimbardo, G. 2009, Sol. Phys., 259, 87
Nóbrega-Siverio, D., & Moreno-Insertis, F. 2022, ApJ, 935, L21

Panesar, N. K., Sterling, A. C., Moore, R. L., & Chakrapani, P. 2016, ApJ, 832, L7

Panesar, N. K., Sterling, A. C., Moore, R. L., et al. 2018, ApJ, 868, L27

Panesar, N. K., Tiwari, S. K., Moore, R. L., Sterling, A. C., & De Pontieu, B. 2022, ApJ, 939, 25

Pariat, E., Antiochos, S. K., & DeVore, C. R. 2009, ApJ, 691, 61

Pariat, E., Antiochos, S. K., & DeVore, C. R. 2010, ApJ, 714, 1762

Pariat, E., Dalmasse, K., DeVore, C. R., Antiochos, S. K., & Karpen, J. T. 2015, A&A, 573, A130

Pesnell, W. D., Thompson, B. J., & Chamberlin, P. C. 2012, Sol. Phys., 275, 3

Prasad, A., Bhattacharyya, R., Hu, Q., Kumar, S., & Nayak, S. S. 2018, ApJ, 860, 96

Prasad, A., Dissauer, K., Hu, Q., et al. 2020, ApJ, 903, 129

Prasad, A., Kumar, S., Sterling, A. C., et al. 2023, A&A, 677, A43

Priest, E. R., & Démoulin, P. 1995, J. Geophys. Res., 100, 23443

Quintero Noda, C., Schlichenmaier, R., Bellot Rubio, L. R., et al. 2022, A&A, 666, A21

Raouafi, N. E., & Stenborg, G. 2014, ApJ, 787, 118

Raouafi, N. E., Patsourakos, S., Pariat, E., et al. 2016, Space Sci. Rev., 201, 1

Rimmele, T. R., Warner, M., Keil, S. L., et al. 2020, Sol. Phys., 295, 172

Robustini, C., Leenaarts, J., de la Cruz Rodriguez, J., & Rouppe van der Voort, L. 2016, A&A, 590, A57

Rouppe van der Voort, L., De Pontieu, B., Scharmer, G. B., et al. 2017, ApJ, 851, L6

Ruan, G., Schmieder, B., Masson, S., et al. 2019, ApJ, 883, 52

Scharmer, G. B. 2006, A&A, 447, 1111

Scharmer, G. 2017, SOLARNET IV: The Physics of the Sun from the Interior to the Outer Atmosphere, 85

Scharmer, G. B., Bjelksjo, K., Korhonen, T. K., Lindberg, B., & Petterson, B. 2003, SPIE Conf. Ser., 4853, 341

Scharmer, G. B., Narayan, G., Hillberg, T., et al. 2008, ApJ, 689, L69

Scharmer, G. B., Sliepen, G., Sinquin, J. C., et al. 2024, A&A, 685, A32

Scherrer, P. H., Schou, J., Bush, R. I., et al. 2012, Sol. Phys., 275, 207

Schmieder, B. 2022, Front. Astron. Space Sci., 9, 820183

Schmieder, B., & Pariat, E. 2010, Astrophys. Space Sci. Proc., 19, 505

Schmieder, B., Malherbe, J. M., Mein, P., et al. 1996, ASP Conf. Ser., 111, 43

Schmieder, B., Aulanier, G., Demoulin, P., et al. 1997, A&A, 325, 1213

Schmieder, B., Archontis, V., & Pariat, E. 2014, Space Sci. Rev., 186, 227

Shen, Y., Liu, Y., Su, J., & Deng, Y. 2012, ApJ, 745, 164

Shibata, K., Nozawa, S., & Matsumoto, R. 1992, PASJ, 44, 265

Singh, K. A. P., Isobe, H., Nishizuka, N., Nishida, K., & Shibata, K. 2012, ApJ, 759, 33

Sterling, A. C., Moore, R. L., Falconer, D. A., & Adams, M. 2015, Nature, 523, 437

Sterling, A. C., Moore, R. L., & Panesar, N. K. 2022, ApJ, 927, 127

Suematsu, Y., Tsuneta, S., Ichimoto, K., et al. 2008, Sol. Phys., 249, 197

Tian, H., DeLuca, E. E., Cranmer, S. R., et al. 2014, Science, 346, 1255711

Titov, V. S., & Hornig, G. 2002, Adv. Space Res., 29, 1087

Tiwari, S. K., Moore, R. L., De Pontieu, B., et al. 2018, ApJ, 869, 147

Tsuneta, S., Ichimoto, K., Katsukawa, Y., et al. 2008, Sol. Phys., 249, 167

Uddin, W., Schmieder, B., Chandra, R., et al. 2012, ApJ, 752, 70

Van Noort, M., Rouppe Van Der Voort, L., & Löfdahl, M. G. 2005, Sol. Phys., 228, 191

Wang, H., Denker, C., Spirock, T., et al. 1998, Sol. Phys., 183, 1

Wang, H., Gallagher, P., Yurchyshyn, V., Yang, G., & Goode, P. R. 2002, ApJ, 569, 1026

Wyper, P. F., Antiochos, S. K., & DeVore, C. R. 2017, Nature, 544, 452

Wyper, P. F., DeVore, C. R., & Antiochos, S. K. 2019, MNRAS, 490, 3679

Yokoyama, T., & Shibata, K. 1995, Nature, 375, 42

Yokoyama, T., & Shibata, K. 1996, Astrophys. Lett. Commun., 34, 133

Zhang, Q. M., & Ji, H. S. 2014, A&A, 567, A11