Open Access
Issue
A&A
Volume 695, March 2025
Article Number A45
Number of page(s) 16
Section The Sun and the Heliosphere
DOI https://doi.org/10.1051/0004-6361/202453322
Published online 03 March 2025

© The Authors 2025

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

Estimates of the emergence rate of the global magnetic flux show a staggering amount of magnetic flux, in some cases more than 1025 Mx, emerging through the solar surface each day (Thornton & Parnell 2011). One important contribution to the global magnetic flux budget is the small-scale magnetic flux in the quiet Sun, on spatial scales ranging from 5″ down to the current observational limit. It supplies between 1023 to 1024 Mx of the total magnetic flux at anytime (Simon et al. 2001; Hagenaar et al. 2003; Zhou et al. 2013; Gošić 2015), which is similar to the amount present in active regions during a solar maximum (Jin et al. 2011). Much of this magnetic flux is provided by network regions (for a review, see Bellot Rubio & Orozco Suárez 2019, and references therein), defined as regions of relatively strong and vertical magnetic fields occupying the outer boundaries of supergranular cells. Gošić et al. (2014) reported that 40% of the magnetic flux that emerged inside of the supergranular cells, in the internetwork regions (see Bellot Rubio & Orozco Suárez 2019, and references therein), ended up in the network, transferring flux at a rate of 1.5 × 1024 Mx per day. This highlights the importance of the extremely dynamic internetwork regions in supplying magnetic flux to the network.

In internetwork regions, the emergence of the magnetic flux mainly occurs in the form of small-scale magnetic loops within granules, which are advected towards the surface by the upward motion of the plasma. Spectro-polarimetric measurements of Hinode allowed researchers to trace small-scale magnetic loops from their initial appearance (Martínez González et al. 2007; Centeno et al. 2007; Gömöry et al. 2010; Ishikawa et al. 2010). A linear polarisation patch appears when the loop emerges. Simultaneously or shortly after, the footpoints are visible as longitudinal magnetic field patches of opposite polarity connected to the linear polarisation patch. Finally, the footpoints move apart and are dragged to the intergranular lanes, while the apex disappears. Small-scale loops can rise and eventually reach the chromosphere transferring energy into upper layers of the atmosphere (Martínez González & Bellot Rubio 2009; Ishikawa et al. 2010; Smitha et al. 2017).

Another mode of flux emergence within granules, namely as magnetic flux sheets, was discovered in the framework of 3D magnetoconvection simulations by Moreno-Insertis et al. (2018) in addition to the emergence of small-scale loops. Those authors found that these two emergence modes differ, first of all, in the fraction of granular surface they cover along their evolution, and in the granular evolutionary stage in which they appear. On the one hand, small-scale magnetic loops are concentrated, tube-like magnetic structures that emerge in well-formed granules. On the other hand, flux sheet emergence results from the stretching of a magnetic flux structure near the photosphere in expanding phases of the host granule, leading to an organised horizontal magnetic ‘mantle’ as was obtained in the simulation of Moreno-Insertis et al. (2018) and also in the observations of Fischer et al. (2019). In the latter paper, as the sheet expands towards the granule’s edges, it eventually fragments, leaving only the footpoints, which can be observed as opposite-polarity patches in circular polarisation maps. In fact, the magnetic flux sheet that emerged exhibited larger magnetic flux than the standard loop emergence and showed small clusters of mixed-polarity elements surrounding it, similar to those later analysed by Gošić et al. (2022). At an even stronger emergent magnetic flux, the magnetic structure might interact with the convective plasma during emergence, shaping the flux appearance. In special cases such as active regions, this interaction leads to elongated granulation, indicating the arrival of strong transverse magnetic flux (Centeno 2012; Centeno et al. 2017; Wang et al. 2020).

High-resolution observations of the flux emergence process at granular scales can help to validate the mechanisms of small-scale magnetic field generation identified in numerical simulations (Rempel et al. 2023). Statistical characterisation of magnetic flux sheets is a way to shed light on how they emerge at the surface and interact with the photospheric flows, how large their contribution is to the energy budget of the solar atmosphere, and on the conditions under which they can reach chromospheric heights. This work aims to provide a statistical analysis of photospheric magnetic flux sheets and quantify the amount of magnetic flux they contribute at the photospheric level. We study the plasma dynamics underlying flux sheet emergence in the host granule and the relation with different stages and structural features of the granules such as exploding granules (Rast 1995; Palacios et al. 2012; Fischer et al. 2017), granules in their initial stages of evolution or nascent granules (Bray et al. 1984; Stix 2002; Moreno-Insertis et al. 2018), and granular lanes (Steiner et al. 2010; Fischer et al. 2020).

2. Observations and methods

2.1. Calibrated database

To investigate the emergence of magnetic flux sheets associated with a growing granule seen in the continuum, we utilised data from two observatories: the Solar Optical Telescope (SOT, Tsuneta et al. 2008) of the Hinode satellite (Kosugi et al. 2007) and the ground-based Swedish 1-m Solar Telescope (SST, Scharmer et al. 2003). Spectropolarimetric observations were acquired by the Hinode Spectropolarimeter (SP, Lites et al. 2013) and the Imaging Spectropolarimeter (CRISP, Scharmer et al. 2008), respectively. All datasets were recorded during very low solar activity, near solar minima. Tables A.1 and A.2 list the specific datasets analysed in this study. Figure 1 presents two representative snapshots, highlighting two distinct cases of emerging magnetic flux sheets.

thumbnail Fig. 1.

Representative snapshots of two distinct cases of emerging magnetic flux sheets. Left: Hinode SP full-sized and reconstructed image from slit data in the continuum of the Fe I 630.2 nm line pair with a FoV of 4 . 5 × 41 $ 4{{\overset{\prime\prime}{.}}}5 \times 41{{\prime\prime}} $ recorded on October 22 2007 at 18:34 UT. Right: Image of the CRISP instrument in the continuum of the Fe I 630.2 nm line recorded on August 06 2011. The FoV is about 41″ × 41″ cut from the original size of 57″ × 57″ to fit with the Hinode SP image. The scale is presented in megametres [Mm]. Yellow boxes mark the location of a detected event in each particular dataset: event ID 523 for the Hinode SP image and event ID 53961 for the CRISP image, respectively. Green, red, and blue contours indicate the apparent transverse magnetic flux, the negative longitudinal magnetic flux, and the positive longitudinal magnetic flux, respectively.

2.1.1. Hinode data

To study the magnetic flux-sheet emergence within a growing granule seen in the continuum, we favoured high-spatial-resolution (no binning) datasets of SP, which have a pixel size of 0 . 16 $ 0{{\overset{\prime\prime}{.}}}16 $. These datasets have a spectral coverage of 2.39 Å centred in the Fe I 630.15 and 630.25 nm lines (Landé-factor of geff = 1.66 and geff = 2.5, respectively) and a spectral sampling of 21.55 mÅ in its normal operation mode. We selected the datasets with a cadence of less than or equal to a minute, to catch the events which sometimes last less than 5 minutes and to identify the different stages of the event. A disadvantage of the SP datasets is that one obtains a relatively small field of view (FoV). In addition, one needs to exclude events which only take place partially within the FoV, and therefore the effective FoV is further reduced. We try to alleviate the situation by collecting several datasets taken at different times and days, thereby, expanding the viewable area by combining multiple datasets.

We obtained suitable data from one of the online repositories1 choosing datasets close to the disc centre. In total, we have obtained over 15 hours of a FoV of approximately 2 . 7 × 41 $ 2{{\overset{\prime\prime}{.}}}7 \times 41{{\prime\prime}} $ and about 7 hours of a FoV of approximately 4 . 5 × 41 $ 4{{\overset{\prime\prime}{.}}}5 \times 41{{\prime\prime}} $. The details of the chosen datasets can be found in Table A.1.

All datasets were calibrated using IDL routines provided by the instrument scientists and publicly available through the SolarSoft distribution. The SP calibration provides us with the apparent transverse and longitudinal magnetic flux density calculated by the sp_prep.pro routine (see Lites et al. 2008), which we mostly use in this work. This magnetic flux corresponds to the magnetic field obtained for the area of the resolution element with filling factor one as if the magnetic field covered the entire resolution element.

2.1.2. SST data

We used science-ready CRISP data from August 2011. The dataset consists of spectropolarimetric data taken in the Fe I 630.15 and 630.25 lines sampled with 44 mÅ steps at 30 wavelength points. The dataset additionally includes spectropolarimetric data at the Fe I 557.6 nm line, which is not used in this work. CRISP recorded spectral images for 47 minutes with a cadence of 28 seconds. The spatial sampling was at 0 . 059 × 0 . 059 $ 0{{\overset{\prime\prime}{.}}}059 \times 0{{\overset{\prime\prime}{.}}}059 $ per pixel with a FoV of about 53″ × 53″. The dataset was calibrated (dark-, flat-, wavelength-, and polarisation calibration) and reconstructed with the multi-object multi-frame blind deconvolution (MOMFBD, Löfdahl 2002; van Noort et al. 2005) algorithm. The output includes Stokes-I, Q, U, and V images at the 30 wavelength points (for more information about the dataset see also Stangalini et al. 2015, 2017; Ledvina et al. 2022). We produced total linear polarisation and circular polarisation maps taking into account both lines of the Fe I 630 nm line pair. We used the resulting time series of polarisation and continuum maps for further analysis.

In addition, we used a dataset collected during a SOLARNET campaign executed in April 2019, which includes spectropolarimetric observations in the Fe I 617.3 nm line (Landé-factor of geff = 2.5) sampled with 15 spectral points symmetrically distributed around the line core (λ0) in 35 mÅ steps. This dataset also includes spectropolarimetric information from Ca II 854.2 nm and spectroscopic information from Hα line, which are not used in this work. The FoV covered by the CRISP images was approximately 50″ × 50″ sampled at a pixel size of 0 . 059 × 0 . 059 $ 0{{\overset{\prime\prime}{.}}}059\times 0{{\overset{\prime\prime}{.}}}059 $. The duration of the time sequence was ∼2 hours with a cadence 28 seconds. The data from CRISP was calibrated (dark-, flat-, wavelength-, and polarisation calibration) using the SSTRED package (Löfdahl et al. 2018) as well as reconstructed with the MOMFBD algorithm. The science-ready data was then made available by the Institute for Solar Physics of Stockholm University (for more information about the dataset see also Fischer et al. 2020; Díaz-Castillo et al. 2024). We produced total linear and circular polarisation maps and aligned the datasets with each other taking the Fe I 617.3 nm data as a reference. In addition to the polarisation maps, we calculated the apparent longitudinal magnetic flux density from the Stokes V profiles using the weak field approximation with the same method applied in Martínez González & Bellot Rubio (2009).

2.2. Magnetic flux sheet detection and characterisation methods

The main criteria used for identifying events of magnetic flux sheet emergence is the evidence of the simultaneous evolution of a transverse magnetic flux patch and a host granule. Therefore, the first task is to find patches of transverse magnetic flux. For the detection and tracking of the transverse magnetic flux, we used the Yet Another Feature Tracking Algorithm (YAFTA) code (Welsch & Longcope 2003), which was applied to the apparent transverse magnetic flux densities calculated using both lines of the Hinode SP data and the total linear polarisation maps of the Fe I 630.25 and 617.3 nm of the SST CRISP datasets of 2011 and 2019, respectively. The minimum size of around 0.1 Mm2 for the transverse magnetic flux patches and threshold parameters were set separately for the Hinode SP or SST CRISP datasets. This step can then be performed fully automatically without visual guidance and allows us to obtain a group of potential events. For the Hinode SP, we set a threshold of 130 Mx cm−2 for the apparent transverse magnetic flux density corresponding to three times the noise level (∼40 Mx cm−2). For the SST CRISP dataset, we set a threshold of 5 × 10−3Ic of the total linear polarisation corresponding to three times the noise level (∼1.7 × 10−3Ic).

The second task is to link the evolution of the transverse magnetic flux to the evolution of a host granule co-located with the transverse magnetic flux patch and evolving simultaneously with it. We explored the granular pattern associated with each event by comparing the temporal evolution of the patches of transverse magnetic flux with spatially aligned continuum images. During our analysis, we found that the automated assignment of a single granule based solely on the variation in the intensity was in some cases rather difficult due to the complexity and fine structure of the convection pattern. Hence, to effectively associate the evolution of a unique granule in a semi-automatic procedure, we implemented a method that combines the horizontal velocity field and its divergence, and the line-of-sight (LoS) or vertical velocity field.

To estimate the horizontal velocity distribution of the plasma flow, we used the DeepVel neural network model (Asensio Ramos et al. 2017) applied to continuum images. DeepVel is a deep learning model based on a fully convolutional neural network architecture that predicts instantaneous horizontal velocities at different heights in the lower solar atmosphere. It directly estimates these velocities pixel-by-pixel from pairs of consecutive continuum images, without time or space averaging. This distinguishes DeepVel from other methods that require such averaging for accurate results such as the standard local correlation tracking methods (November & Simon 1988; Fisher & Welsch 2008). Synthetic data generated by suitable magnetohydrodynamics (MHD) models of the solar atmosphere are used to train DeepVel. These models provide images of the solar surface and the corresponding true velocities at three different heights. DeepVel learns to associate these images with the underlying velocities through a deep learning process. The accuracy of DeepVel’s predictions is closely tied to the specific characteristics of the training data, thus to estimate velocities at different spatial resolutions, observing cadences, or locations on the Sun, the DeepVel model should be retrained using synthetic data that reflects these specific conditions (Tremblay et al. 2018).

For the Hinode SP dataset, we performed a training and testing procedure using two different CO5BOLD simulations described in Vigeesh et al. (2017), one hydrodynamic (non-magnetic) and the other magnetohydrodynamic simulation (initial uniform vertical magnetic field of 50 G and advanced over a magnetic field redistribution timescale). The computational domain has a horizontal size of around 38 × 38 Mm with a cell size of 80 km and discretised on 480 × 480 horizontal grid cells. We used snapshots of intensity at the different cadences of the Hinode SP datasets, which had been degraded with the corresponding PSF profile. We used the bolometric intensity at τ = 1 of the simulations, which reproduces the behaviour of the observed continuum, and matched the histogram distribution of the simulated intensity with the observed continuum intensity. For training and validation of the model, we randomly extracted 64000 pairs of consecutive patches of 40 × 40 pixels from the PSF-degraded continuum snapshots together with the corresponding horizontal velocity field at three different τ values: τ = 1, 0.1, 0.01 from the simulation. As in Asensio Ramos et al. (2017), all inputs were normalised to the median intensity of the quiet Sun and velocities to the interval [0, 1] using training set bounds. The model architecture was implemented using PyTorch modules and the training was carried out by optimising the squared difference between the output of the network and the velocities in the training set with the Adam stochastic first-order gradient-based optimisation algorithm (Kingma & Ba 2014) using a learning rate of 0.001 and batches of 64 samples during 50 epochs.

For the SST dataset, we used the default and pre-trained model from Asensio Ramos et al. (2017). This model was trained on synthetic magneto-convection simulations and tested on SUNRISE I/IMAX (Barthol et al. 2011; Martínez Pillet et al. 2011); it is applicable to CRISP data successfully as there is only a slight difference between the SUNRISE I/IMAX data pixel size and the CRISP pixel size. We noted, however, that the model was trained considering snapshots with a cadence of 30 seconds, which may lead to an underestimation of velocities compared to the 28-second cadence of SST CRISP data used in this work.

While the horizontal velocities inferred by the DeepVel neural network are model-dependent, it accurately reconstructs the flow orientation and amplitude at subgranular and granular scales when using a training and validation input generated from different numerical simulations (Tremblay et al. 2018).

To estimate vertical velocity fields, we calculated the bisector velocity as a proxy for LoS velocities at specific line depths thus sampling different heights in the photosphere. Due to the different spectral resolutions of the Hinode SP and CRISP and the different formation heights of the lines, we selected lines and defined three line depths that roughly characterise the low, middle and high photosphere. For the Hinode SP spectra, we selected depth positions at 70%, 35% and core (0%) of the Fe I 630.15 nm line for the low, middle and high photosphere, respectively. For SST CRISP, we selected depth positions at 60%, 30% and core (0%) of the Fe I 630.15 nm line (2011 dataset) and of the Fe I 617.3 nm line (2019 dataset) for the low, middle and high photosphere, respectively. We also performed p-mode correction to the full sequence of the datasets.

Based on the description of Moreno-Insertis et al. (2018), the emergence of magnetic flux sheets can take place in the early stages of the host granule, which are characterised by rapid and expanding plasma flows. We therefore used the positive transverse divergence of the flow to determine the area of the host granule. In addition to the horizontal plasma flow, the vertical plasma flow is expected to be upward-directed inside expanding granules and gradually weaken at higher altitudes.

Considering the above, we defined two initial criteria which characterise a granule hosting a magnetic flux sheet: (1) closed regions of positive divergence and (2) closed regions of upward velocities, both co-located with the patch of transverse magnetic flux. To apply such criteria, we applied the DeepVel models to the time series of continuum intensity images of all datasets, to retrieve the apparent horizontal velocity field, vh, and then calculate its divergence, ∇h ⋅ vh(x, y). We selected the plasma flows at the lowest height in the photosphere to identify the host granule; that is, the apparent horizontal velocity field at τ = 1 provided by DeepVel and the vertical velocity field at the 70% and 60% of the line depth provided by the bisector analysis for Hinode SP and SST CRISP, respectively.

Our method applies the first criteria automatically as follows: (1) locate the superposition of regions of positive divergence of the apparent horizontal velocity field with the patch of transverse magnetic flux and (2) define the host granule boundaries using a single closed contour of zero-divergence. We noted that some events show a complex granulation (CG) pattern in their surroundings and that spatial resolution makes the zero-divergence criteria too sensitive to contrast changes resulting in merged granules. Thus, we manually inspected the distribution of the vertical velocity field inside the mask to correct merging granules. In some cases, we applied both criteria strictly even though this may result in a smaller granule. In addition, many of the events were associated with exploding granules and granular lanes (see Sect. 3 for details on these phenomena and findings). This leads to cases in which a part of the granule (a newly split or sectioned area) starts to expand again and develops its own convective pattern. In some cases, the expanding section is viewed as the host granule when is responsible for surfacing the associated magnetic flux.

Finally, the third task was to determine if the newly arriving transverse magnetic flux patch is flanked by appearing and diverging opposite polarity patches of longitudinal magnetic flux density, making this a true emerging event. We defined the longitudinal magnetic flux patches associated with the footpoints by proximity to the sheet. We selected the longitudinal magnetic flux patches within a radial distance corresponding to a circle centred on the sheet, with a maximum area of two times the maximum area of the sheet, and we defined the boundaries of the longitudinal magnetic flux patches with a threshold. To establish the threshold values, we performed a similar calculation of the correspondent magnetic flux density retrieved by a particular photon noise as presented by Martínez González & Bellot Rubio (2009). We set 10 Mx cm−2 and 15 Mx cm−2 for Hinode SP and SST CRISP datasets, respectively, corresponding to three times the photon noise in each case.

In both the Hinode SP and SST CRISP datasets, we visually inspected the results and, when necessary, manually re-adjusted boundaries. This visual inspection was crucial for confirming the effective presence of magnetic flux sheet cases. Feature tracking methods such as YAFTA identified and tracked all transversal magnetic flux patches in the considered datasets, encompassing all potential cases of loop-like magnetic flux emergence. Hundreds of such patches were detected. However, to accurately identify sheet emergence, in addition to the established criteria, we visually verified all cases to minimise the occurrence of false detections.

Figure 2 shows an example of a magnetic flux sheet emergence event in the Hinode SP dataset characterised by our method. All panels correspond to a single time step where the magnetic flux sheet has maximum extent over its host granule.

thumbnail Fig. 2.

Magnetic flux sheet event with ID 352 recorded by Hinode SP on September 26 2007 at 08:15 UT. From left to right, the first panel shows the continuum intensity in the background with contours. Pine green contours mark a threshold of 130 Mx cm−2 for the apparent transverse magnetic flux density patch ϕTapp. The blue (red) contours mark the positive +ϕLapp (negative −ϕLapp) polarity of the longitudinal magnetic flux density with thresholds of 20 Mx cm−2, 40 Mx cm−2 and 60 Mx cm−2. The black contours mark a single closed contour of zero-divergence defined as the boundary of the host granule. The second panel shows the divergence of the apparent horizontal velocity field marking the corresponding location of the sheet represented by the transverse magnetic flux density patch ϕTapp (pine green contour) and the boundary of the host granule (black contour). Finally, the third, fourth, and fifth panels show the bisector velocity field for the three different depth positions: 70%, 35%, and core (0%) of the Fe I 630.15 nm line corresponding roughly to the low, middle, and high photosphere, respectively.

3. Results

3.1. Occurrence rate, description, and properties of events

After detecting events as was described in Sect. 2.2 and discarding questionable candidates (mostly those located at the edges of the FoV and those with incomplete evolution information), we confirm 42 events as positive identifications (see full databases in Tables B.1 and B.2). This number likely represents a lower limit of the actual occurrence rate of magnetic flux sheets. The numerical simulations by Moreno-Insertis et al. (2018) estimate 0.3 to 1 events per day and Mm2. However, as the authors also note, this estimate is based on a small sample found in the simulations and might not be representative. Assuming the lower estimate of 0.3 events per day and Mm2, we would expect to detect approximately 45 events in our database, which is a number comparable to the 42 events we identified.

As an initial classification that helps to characterise the magnetic sheet events, we define four types of granular phenomena linked to the host granule, where the emergence of the magnetic flux sheet occurs: exploding granules (EG), nascent granules (NG), granules with granular lanes (GL), and CG. We found that the first three types of granular phenomena associated with magnetic flux sheet appearance are roughly equally represented in our sample.

The most frequent are magnetic flux sheet events associated with exploding granules, with 15 cases identified. Exploding granules appear in continuum images as rapidly expanding, often large granules that develop a dark core or dot (Hirzberger et al. 2001; Roudier et al. 2001; Ellwarth et al. 2021). This dark core indicates a developing downflow within the granule, which eventually evolves into a new intergranular lane. In addition to their morphological features, exploding granules have been associated with nearby magnetic flux dynamics (De Pontieu 2002; Zhang et al. 2009). In particular, Palacios et al. (2012) and Guglielmino et al. (2020) studied the emergence of extended magnetic structures in exploding granules observed by SUNRISE I/IMAX, carrying total magnetic fluxes of the order of 1018 Mx. Using numerical simulations, Rempel (2018) reported on the process of amplification of magnetic flux by turbulent shear in the newly forming downflow within the exploding granule. This process results in the formation of extended magnetic structures nearby, with lengths comparable to the granule scale.

The next phenomena in the ranking are nascent granule events, with 12 cases, and granules with granular lanes, with 11 cases. In the nascent granule events, the magnetic flux sheet appears nearly simultaneously with the host granule itself and evolves with it during its fast expansion phase. These are the cases mentioned in the numerical paper of Moreno-Insertis et al. (2018) as a probable cause of the formation of magnetic flux sheets. The average size of nascent granules in this sample is approximately 0.8 Mm2, which is smaller than the average size of host granules classified as exploding granules and host granules with granular lanes (1.7 Mm2). However, the average maximum area reached by nascent granules in the sample is around 1 Mm2, a characteristic size of standard granules. This suggests that these sheets tend to disappear (or become undetectable) in the advanced phases of the granular evolution.

Granular lanes, on the other hand, appear in continuum images as dark lanes moving within a granule, preceded by a bright rim. They have an arch-like shape, with their endpoints still connected to the intergranular lane as they move apart. Numerical simulations by Steiner et al. (2010) suggest that these features are the visible signatures of horizontal vortex tubes beneath the solar surface, which can be associated with magnetic fields in the photosphere, as confirmed observationally by Fischer et al. (2020).

We identified four cases with CG patterns where the host granule could not be clearly determined. As a result, we excluded these cases from the statistics on the host granule coverage and the plasma dynamics of host granules.

Figure 3 summarises the temporal and spatial properties of the events concerning the associated granular type. Columns 4−5 of Table B.1 and Col. 4 of Table B.2 show the absolute value for each quantity analysed. We found no clear relationship between the duration of magnetic flux sheet events and the associated granular type, as is shown in Fig. 3a. The duration of the events ranges from 3 to 10 minutes, with an average of 6 minutes. Similarly, the maximum area reached by the transversal magnetic flux density patch, from now on termed magnetic flux sheet area, does not show a clear correlation with the granular type, as is shown in Fig. 3b, with most magnetic flux sheets exhibiting areas between 0.25 and 0.75 Mm2.

thumbnail Fig. 3.

Duration of each event (left panel), the apparent maximum area reached by the magnetic flux sheet (central panel) and the proportion of a host granule’s area covered by its associated magnetic flux sheet (right panel). The granulation type is coded by different colours, as follows: exploding granules (EG): yellow bars and lines; granular lanes (GL): blue bars and lines; nascent granules (NG): green bars and lines; and complex granulation (CG): red bars and lines. The CG is only included in panels a and b. All histograms include the kernel density estimate (KDE) providing complementary information about the possible shape of the distribution. The actual values are included in Cols. 4−5 of Table B.1 and Col. 4 of Table B.2. (a) Event duration. (b) Magnetic flux sheet maximum area. (c) Host granule coverage.

However, the proportion of a host granule’s area covered by its associated magnetic flux sheet (host granule coverage) appears to be linked to the granular type (see Fig. 3c). Nascent and smaller granules (< 1 Mm2) are generally covered to a larger extent by the magnetic flux sheet compared to larger granules (> 1 Mm2) (green histogram in Fig. 3c). Magnetic flux sheets associated with granular lanes typically cover a smaller area in the host granule (areas between 1 Mm2 to 2 Mm2), as they tend to be elongated and located near the edges of granules, which prevents them from covering it completely (blue histogram in Fig. 3c). The distribution of the host granule coverage for events associated with exploding granules seems broader (yellow histogram in Fig. 3c). This could be due to two factors: (1) a wider range of the sheet sizes (blue histogram in Fig. 3b) or a wide range of host granule sizes (areas from 1 Mm2 to 5 Mm2), which will depend on whether the host granule is the entire exploding granule or a fragment of it when the sheet reaches its maximum. We found that the maximum coverage measured in our event list reached only 70%, which could be related to several factors: the criteria used to define the boundaries of the host granule or the intrinsic properties of the granule or sheet.

In the following subsections, we examine the evolution of selected cases for each granulation phenomenon. The simplest, standard ‘textbook-like’ cases of flux sheet emergence in our sample are found for the category of nascent granules (Sect. 3.1.2 and, in particular, Fig. 7). Yet, we start describing individual flux sheet emergence events with examples of exploding granules since they are the most frequent category in the sample. We also illustrate the diverse characteristics of emerging magnetic flux sheets and highlight, from the beginning, various difficulties in their identification.

3.1.1. Exploding granule events

Figure 4 shows the emergence of a magnetic flux sheet associated with an exploding granule. This event, labelled ID 801, was detected by Hinode SP during the September 25 2007 at 12:59 UT run. The magnetic flux sheet (pine green contour) becomes visible at the boundary of the exploding granule (black contour) at t = 121.2 min. The characteristic dark dot is observed around (x, y) = [1.0, 0.9] at t = 122.4 min and t = 123.6 min while the magnetic flux sheet expands. Along with the sheet, both footpoints of longitudinal magnetic flux are visible as small patches flanking the sheet (red contours for the negative polarity and blue contours for the positive polarity). The main fragment of the exploding granule is identified as the host granule (black contour), characterised by upflows of −2 km s−1 (blue pixels on the colour map within the host granule). Over the five-minute evolution, the magnetic flux sheet quadruples in size, the footpoints expand and intensify, and the host granule appears to grow in area. In this event, the positive-polarity footpoint remains anchored to the intergranular region from the start of the emergence, while the negative-polarity footpoint appears within the host granule. After t = 124.8 min, the negative-polarity footpoint fragments into two patches and a second patch of positive polarity appears surrounding the left flank of the magnetic flux sheet patch. This observation is similar to the organisation of the longitudinal magnetic flux described by Moreno-Insertis et al. (2018) during the later stages of magnetic flux sheet emergence. In their study, the two circular polarisation patches, representing the footpoints, gradually move apart and stretch over time, eventually reaching the intergranular lanes.

Figure 5 depicts the complete evolution of a magnetic flux sheet event within an exploding granule observed by SST CRISP. This event, identified as ID 91772, is detected in the 2011 dataset and is associated with a ‘double-event’ apparently related to a large-scale magnetic structure accumulating beneath the photosphere and emerging in sections. A small transverse magnetic flux patch appears below the typical central dark dot in the exploding granule at t = 42.9 min (pine green contour). Around the main transverse magnetic flux patch, non-negligible signals of linear polarisation located within neighbouring granules and associated with the distribution of the longitudinal magnetic flux are detected. At t = 43.9 min, the granule begins to disintegrate, forming intergranular lanes associated with downflows of 1 km s−1 (red pixels on the bisector velocity map within the host granule). At the same time, the transverse magnetic flux patch remains in a region dominated by upflows (blue pixels on the bisector velocity map within the host granule). One intergranular lane appears on the granule’s left side far from the central dot, while the other emerges on the right connected to the central dot. By t = 44.8 min, the transverse magnetic flux patch reaches its maximum extent and a fragment of the original granule becomes the host granule (black contour). The longitudinal magnetic flux patches linked to the footpoints of the sheet are located within the intergranular lane on the host granule surroundings, where they undergo splitting and merging with nearby longitudinal magnetic flux concentrations over the next three minutes. From t = 46.2 min onwards, a secondary event, identified as ID 94010, appears within a neighbouring granule around (x, y) = [3.1, 2.1] (orange contour in Figure 5). The surrounding longitudinal magnetic flux in the FoV shows a complex distribution associated with this event, evidencing a potential fragmentation of a large-scale magnetic structure while emerging from below the solar surface.

thumbnail Fig. 4.

Magnetic flux sheet emergence event alongside the expansion of its host granule, which fragments as an exploding granule. The image series shows the continuum maps reconstructed from Hinode SP slit data. The times are indicated above the images. Contours mark the apparent transverse (pine green) and longitudinal (red contours for the negative polarity and blue contours for the positive polarity) magnetic flux density associated with the event. The magnetic flux density thresholds are defined as in Fig. 2. Black contours mark the host granule border. Colour maps within the host granule show the distribution of the bisector velocity at 70% of the line depth.

thumbnail Fig. 5.

Second case of magnetic flux sheet emergence in an exploding granule. The image series shows the continuum close to the Fe I 617.3 nm line recorded by CRISP. The plot set-up is defined as in Fig. 4. The blue (red) contours mark the positive (negative) polarity of the longitudinal magnetic flux density with thresholds of 40 Mx cm−2, 60 Mx cm−2 and 80 Mx cm−2. Black contours mark the host granule border. A secondary magnetic flux sheet, depicted by orange contours, was detected near and concurrently with the event within the exploding granule. Colour maps within the host granule show the distribution of the bisector velocity at 60% of the line depth.

3.1.2. Nascent granule events

The event depicted in Fig. 6 showcases a magnetic flux sheet that appears over a nascent, isolated and relatively small granule. This event, labelled ID 458, was detected by Hinode SP during the October 22 2007 at 18:34 UT run. The maximum area of the transverse magnetic flux patch of 0.4 Mm2 is reached at t = 89.0 min (pine green contour), just a minute after the initial detection and persisting for the following three minutes. The area of the host granule (black contour) and the area of the magnetic flux sheet (pine green contour) are similar and do not change after t = 89.0 min. The vertical velocity within the host granule at t = 89.0 min and t = 90.0 min is characterised by upflows between −0.5 km s−1 to −1.25 km s−1 (blue pixels on the colour map within the host granule). These upflows transitioned to descending flows during the final two minutes, ultimately leading to the disappearance of both the sheet and the host granule, giving rise to an intergranular region. This event represents one of the smallest cases detected, exhibiting a relatively weak magnetic flux contribution. The spatial resolution of the Hinode SP dataset, coupled with the faint appearance of the event, significantly complicates the identification of its footpoints. As a result, we only consider a single footpoint with negative polarity, detected just above the detection threshold at the onset of the event.

thumbnail Fig. 6.

Case of a magnetic flux sheet emergence in a nascent and relatively small granule. The image series shows the continuum maps constructed from Hinode SP slit data. The plot set-up and information are as shown in Fig. 4.

Figure 7 depicts the evolution of a magnetic flux sheet event within a uniform granule in its fast-expanding phase observed by SST/CRISP. This event, identified as ID 58352 and detected in the 2019 dataset, is a clear example of a magnetic sheet covering a considerable portion of its host granule. Although the host granule appears to be well established in the initial continuum images, the vertical velocity field exhibits characteristics reminiscent of the early expansion phase of a nascent granule. The transverse magnetic flux patch appears within the host granule with an area below 0.1 Mm2 at t = 66.7 min, later on growing up to 0.5 Mm2 at t = 67.7 min. The longitudinal magnetic flux patches corresponding to the footpoints appear in the flanks of the sheet, and the host granule develops an intergranular lane at the bottom of the sheet co-located with the positive-polarity footpoint (blue contours at t = 67.7 min and t = 68.6 min). On the opposite side, the negative-polarity footpoint appears within the granule in a region of apparent downflow. At t = 69.5 min, the magnetic flux sheet reaches its maximum extent of 0.68 Mm2 and covers around 60% of the host granule. One can observe that during that time, the host granule has a localised region of upflow fully covered by the sheet (blue pixels on the colour map in the bottom part of the host granule) and regions of downflow mostly localised at the edges on the upper part of the granule. In particular, the region with larger downflow velocities is co-located with the negative-polarity footpoint that moves towards the upper zone of the host granule until it reaches the intergranular region at t = 71.4 min.

thumbnail Fig. 7.

Second case of magnetic flux sheet emergence in a nascent granule. The image series shows the continuum close to the Fe I 617.3 nm line recorded by CRISP. The plot set-up and information are as shown in Fig. 5.

3.1.3. Granular lane events

We observed that events associated with granular lanes often coincide with granule fragmentation or in the formation of a new convective pattern. Figures 8 and 9 exemplify these cases, the first one identified as ID 440 and detected by the Hinode SP during the October 06 2007 08:01 UT run and the second one identified as ID 39266 detected in the SST CRISP dataset of 2019, respectively.

thumbnail Fig. 8.

Case of a magnetic flux sheet emergence associated with a granular lane. The image series shows the continuum maps constructed from Hinode SP slit data. The plot set-up and information are as is shown in Fig. 4. The dashed yellow curves in t = 69.6 min and t = 70.8 min highlight the granular lane signature.

thumbnail Fig. 9.

Case of a magnetic flux sheet emergence associated with a granule with granular lanes. The image series shows the continuum close to the Fe I 617.3 nm line recorded by CRISP. The plot set-up and information are as is shown in Fig. 5. The dashed yellow curve in t = 49.0 min highlights a granular lane signature. The yellow arrow marks the direction of displacement of the magnetic flux sheet.

Figure 8 depicts the evolution of a magnetic flux sheet associated with a granular lane signature. The patch of transverse magnetic flux appears at the edge of an emerging granule at t = 66.6 min. It expands over the next two minutes, reaching a size of 0.58 Mm2 but only covering 20% of the total area of the host granule. Signatures of the granular lane are visible at t = 69.6 min and t = 70.8 min (dashed yellow curves). The sheet is co-located with a dark area of relatively small upflow velocities (see the dashed yellow curve at t = 69.6 min), preceded by a bright rim around (x, y) = [1.7, 2.0], that eventually lead to the formation of a new convective pattern after t = 71.4 min. This event occurs near a strong network element, which does not seem to be associated with the sheet, therefore, we take as the effective footpoint the longitudinal magnetic flux patch of positive polarity closest to the sheet, which is only observed in the first two minutes of evolution in (x, y) = [2.0, 1.5]. In this case, the footpoint of negative polarity is not observed.

Figure 9 shows an example of a magnetic flux sheet with granular lane signatures evolving within a complex-shaped granule that eventually fragments. The patch of transverse magnetic flux appears at t = 49.0 min in an arch-like border of the host granule in its upper part similar to an extended granular lane, around (x, y) = [3.0, 3.2] (dashed yellow curve). At t = 50.4 min, the host granule starts to fragment while the patch of transverse magnetic flux expands. The patch of transverse magnetic flux reaches its maximum extent of 0.2 Mm2 at t = 51.8 min when the host granule fragments in three smaller granules. At that moment, the patch of transverse magnetic flux is localised over a fragment characterised by upflow velocities of −1.2 km s−1, a flow likely contributing to the sheet’s rise. Over the next four minutes, the area of the patch of transverse magnetic flux decreases and moves to the right edge of the new host granule (towards the direction marked by the yellow arrow). The footpoints associated with this sheet are observed at t = 51.8 min, the negative-polarity one localised on an intergranular region, around (x, y) = [3.3, 3.4], and the positive-polarity one within the host granule. Starting from t = 53.2 min, the footpoints extend around the new host granule, until it either fragments (negative-polarity footpoint) or continues to expand (positive-polarity footpoints). The magnetic flux sheet exhibits a displacement towards the direction indicated by the yellow arrow. This displacement, likely driven by plasma flows, encompasses both the sheet itself and its footpoints, which extend along the direction of displacement. Notably, the pair of footpoints appears to be oriented perpendicular to the sheet’s displacement direction (towards the yellow arrow marks).

3.2. Emerged flux

To quantify the amount of magnetic flux carried by magnetic flux sheets, we measured the total magnetic flux emerging from each studied event and compared it with the values reported in the literature. We used the tracked longitudinal magnetic flux patches associated with the events. The patches were defined considering a threshold of three times the photon noise calculated for each dataset, as was described in Sect. 2.2. In addition, we determined by proximity which of the detected longitudinal magnetic flux patches were involved in the particular emergence event to select the corresponding footpoints, which would often be several of those. The total magnetic flux, ΦLapp, was calculated by averaging the magnetic flux density of the pixels within the patch and then multiplying by the area covered by the patch limited by the defined threshold. For the calculation, we used the apparent magnetic flux density output of the sp_prep.pro routine for the Hinode SP dataset and the magnetic flux densities calculated by the weak field approximation for the SST CRISP dataset. We analysed the total magnetic flux during each event’s lifetime, as well as for both associated footpoints. For events for which only one footpoint is detected, we used the maximum value reached in the observed footpoint during the event’s lifetime (see examples in Sects. 3.1.2 and 3.1.3). In the case of events for which both footpoints are observed, we noticed that most of them share similar values and evolution patterns, thus we selected the maximum value reached in time and footpoint. For cases for which the values and evolution patterns are not similar, we visually determined the footpoint that best fits the evolution of the magnetic flux sheet.

The histogram in Fig. 10 exhibits the distribution of total longitudinal magnetic flux ΦLapp calculated for each magnetic sheet event. The absolute values are listed in Col. 7 of Table B.1. The average value obtained for the 42 studied events is 1.6 × 1017 Mx with a slight tendency towards higher values. This average is 1.8 times larger than the mean value reported by Martínez González & Bellot Rubio (2009) for Ω-loop type emergence and is comparable to the values of high magnetic flux events in their sample.

thumbnail Fig. 10.

Total longitudinal magnetic flux measured in the footpoints of the magnetic flux sheet events. The histogram includes the kernel density estimate providing complementary information about the shape of the distribution.

To understand the relationship between the emergence of magnetic flux sheets and the magnetic flux carried by their footpoints, we compared the total magnetic flux for each event with three properties of the associated magnetic flux sheet: (i) its maximum area Amax (left panel in Fig. 11), (ii) the rate of change of the magnetic flux sheet area ΔA (central panel in Fig. 11) and (iii) the divergence of the horizontal flow field within the sheet measure at the surface ∇h ⋅ vh(x, y)τ = 1 (right panel in Fig. 11). We find that the total magnetic flux is weakly correlated with these properties, as is indicated by Pearson correlation coefficients that are below 0.5 in all cases. In particular, the correlations between total magnetic flux with the change rate of the magnetic flux sheet area and the flow divergence are very low (correlation coefficient of 0.14), suggesting a negligible impact of total magnetic flux on these horizontal expansion metrics of the sheet. The correlation between total flux and maximum area is slightly stronger (coefficient of 0.41), but remains modest.

thumbnail Fig. 11.

Scatter plots comparing the total longitudinal magnetic flux of the sheet, ΦLapp, to the apparent maximum area of the sheet Amax (left panel), the rate of change of the magnetic flux sheet area ΔA (central panel), and the divergence of the horizontal flow field within the sheet ∇h ⋅ vh(x, y)τ = 1 (right panel). The Pearson correlation coefficient is shown in each panel. See Table B.1 (Cols. 5, 6, 9) for numerical values.

3.3. Plasma dynamics

Magnetic flux sheets are coupled to the evolution of their host granule, thus the dynamics of the plasma within these granules may play a fundamental role in the occurrence and emergence process of the magnetic structure. For instance, studying the sub-photospheric origin of the magnetic flux emergence within granules, Moreno-Insertis et al. (2018) found that a representative amount of the non-weakly magnetised plasma is located in upflow zones within the convection cell. In the present paper, we have analysed the vertical velocity of the granules that host flux sheets to check whether they are emerging during the lifetime of the sheet. Figure 12 shows a box plot of the distribution of the average of the bisector velocities at three positions of the spectral line that sample roughly three different heights in the photosphere: low (closest point to the continuum intensity), middle (roughly half maximum of the depth of the line) and high photosphere (core of the line). The box plots display the full distribution, the percentiles and the median value for each height. First, we notice that most of the distribution has negative velocities, that is, upflows. Going from the upper to the lower level, the median value becomes increasingly negative and the dispersion grows. The actual values go from −0.18 km s−1 in the high photosphere to −0.45 km s−1 in the low photosphere, confirming that most of the host granules are emerging. We identify two outlier events. The first is characterised by positive velocities in the photosphere, indicative of a host granule with a weak upflow. Inspecting the event in detail, its associated sheet emerges near an intergranular region characterised by downward flows. The second outlier involves a host granule with unusually high upward velocities (approaching −1 km s−1) in the middle and high photosphere, significantly exceeding the typical velocities in the sample.

thumbnail Fig. 12.

Average vertical velocity at different heights of the photosphere of the host granules during the emergence of their associated magnetic flux sheet. The box plots display the full distribution (scatter points), maximum and minimum values (whiskers lines), the interquartile range (box), upper and lower quartiles (box boundaries), and the median value (line within the box) for each height: low photosphere (blue box), middle photosphere (white box), and high photosphere (red box).

Measuring the expansion rate of the magnetised elements that rise towards the photosphere, Moreno-Insertis et al. (2018) found that their expansion rate at depths such as 0.5 or 1 Mm was low enough (or negative) that they could maintain their field strength with no large weakening for several minutes. For instance, the expansion rate of the majority of elements with B > 300 G at depths of 0.5 and 1 Mm below the photosphere was below 3 × 10−3 s−1. For upflowing plasma, at the photosphere we expect the expansion rate to be substantially larger than at depths of half or one Mm; yet, it may be of interest to compare the values obtained here with those quoted by those authors.

We analysed two potential metrics for the horizontal expansion of the magnetic flux sheet, namely the rate of change of the magnetic flux sheet area and the divergence of the horizontal flow field within the sheet measured. For the first metric, the average rate of change of the transverse magnetic flux patch area, calculated across the sample, is 1.8 × 10−3 Mm2 s−1 (see central panel in Fig. 11). Assuming a linear, uniform, and circular expansion, the corresponding horizontal velocity is approximately 1.3 km s−1. This value is comparable to the velocities of the convective plasma in granules (∼1 km s−1, Nordlund et al. 2009) and aligns with the joint evolution of the magnetic flux sheet and its host granule. For the second metric, we analysed the distribution of the divergence of the plasma at the location of the magnetic flux sheet while it emerges at the surface; that is, at τ = 1 (see right panel in Fig. 11). Figure 13 shows the distribution of the mean values of the divergence of the horizontal flow field at τ = 1. We find that most of the values in the distribution are positive, with a distinct peak at around 8 × 10−3 s−1, which provides further evidence for the emergent nature of the event. The mean value of the distribution is 7.5 × 10−3 s−1 with considerable contribution of lower values. These numbers are in line with those obtained in the numerical simulation for deeper levels: the present values are larger, as was expected, but still within the same order of magnitude.

thumbnail Fig. 13.

Divergence of the horizontal plasma flow co-located with the magnetic flux sheet area. The horizontal velocity flow is estimated at τ = 1. The histogram includes the kernel density estimate providing complementary information about the shape of the distribution.

4. Summary and discussion

We have presented a statistical analysis of 42 magnetic flux sheet emergence events observed in the solar photosphere. We have characterised the magnetic properties of these events and the associated plasma dynamics of their host granules. To identify these events, we developed a systematic two-step method. First, we detected patches of transverse and longitudinal magnetic flux density based on polarisation signals of photospheric absorption lines. Second, we associated a single convective cell as a host granule for each event by analysing the estimated horizontal and vertical velocity fields at different photospheric heights. With our method, we confirm the distinct properties of magnetic flux sheets in the photosphere compared to subgranular-scale Ω-loop type emergence (Martínez González & Bellot Rubio 2009). We observe that magnetic flux sheets exhibit expansive, long-lasting transverse magnetic flux patches, unlike the concentrated, nearly constant-size, and short-duration structures of Ω-loop events. Our calculations of the expansion of the magnetic flux sheet, as well as the values of the divergence of the flow within the sheet and the lifetime of the transversal magnetic flux patch itself support this distinction. Additionally, we find that flux sheet footpoints often exhibit structures that expand, stretch, and fragment over time, eventually forming clusters of mixed-polarity elements. In contrast, Ω-loop footpoints are typically compact, opposite-polarity patches that separate during emergence.

Our observed occurrence rates, namely 0.27 events per day and Mm2, align well with those predicted by the numerical simulations of Moreno-Insertis et al. (2018), namely between 0.3 and 1 events per day and Mm2. However, given the small sample size in both the observational and numerical studies, it is difficult to draw definitive conclusions about their representativeness. Based on our study, we report the amount of magnetic flux that magnetic flux sheets contribute to the solar photosphere at an average value of 1.6 × 1017 Mx per sheet (see Fig. 10), around 1.8 times larger than the mean value reported by Martínez González & Bellot Rubio (2009) for Ω-loop type emergence (0.91 × 1017 Mx). Considering the estimated average magnetic flux density per sheet and appearance rate, magnetic flux sheet events could potentially contribute an average of 2.6 × 1023 Mx per day to the small-scale magnetic flux that feeds the quiet Sun atmosphere. This corresponds to 24% of the magnetic flux reported by Martínez González & Bellot Rubio (2009) for Ω-loop type emergence and 2% of the amount reported by Jin et al. (2009) for HIFs. We suggest that magnetic flux sheets of the type analysed in this paper could contribute a non-negligible amount of magnetic flux to the total quiet-Sun magnetic flux budget.

To delve into the process of emergence of the magnetic flux sheets, we have explored the relationship between the total longitudinal magnetic flux at their footpoints with two horizontal expansion metrics of the magnetic flux sheet, namely the rate of change of the transverse magnetic flux patch area (magnetic flux sheet area) during the emergence and the divergence of the horizontal velocity field within the sheet (see Sects. 3.2 and 3.3). Our findings suggest that none of these two expansion metrics are influenced by the total longitudinal flux calculated at the footpoints of the magnetic flux sheet within our sample. However, we cannot exclude the possibility that the amount of magnetic flux influences the dynamics of the three-dimensional magnetic flux sheet. Our analysis is limited to the formation heights of lines that roughly estimate different zones in the photosphere and the estimation of horizontal flows not completely in agreement with those zones, leading to two-dimensional views of the magnetic flux sheet, obscuring its full three-dimensional nature and the potential impact of the amount of magnetic flux on its dynamics. Future studies could probe the influence of magnetic flux on plasma dynamics by analysing these magnetic flux volumes at different atmospheric layers using various spectral lines.

Moreover, we have investigated the relationship between the emergence of magnetic flux sheets and granular phenomena. We identified cases where the magnetic flux emerges directly coupled to nascent and rising granules, consistent with the numerical simulations of Moreno-Insertis et al. (2018). The distribution of the vertical velocity of the host granules at various photospheric heights (see Fig. 12) confirms the emergent nature of these granules, especially in the lower photosphere. This distribution exhibits an asymmetric interquartile range, with a larger proportion of smaller velocities, potentially indicating host granules with weaker upflow velocities. While the influence of magnetic flux on convective flows is well-established in cases of strong magnetic flux (Centeno 2012; Wang et al. 2020), the influence on the flows of small, subgranular features like tiny magnetic flux tubes or flux sheets is probably weak, since their plasma beta is likely to be well above unity. However, it may be of interest to investigate this point in future studies.

Additionally, our study revealed that many magnetic flux sheets emerge in association not only with nascent granules but also with exploding granules and granules with granular lanes, exhibiting the diverse characteristics of emerging magnetic flux sheets. In the case of exploding granules, the rapid expansion and fragmentation of these structures might facilitate the emergence of magnetic flux accumulated beneath the photosphere, as is shown by Guglielmino et al. (2020). Numerical simulations of exploding granules demonstrate that plasma flow in the fragmentation zone can amplify the local magnetic field a few hundred kilometres below the surface, which is later transported to the surface by shallow recirculation (Rempel 2018). This could potentially contribute to magnetic flux near the fragmentation zones, in addition to the underlying emerging flux. In the case of granules with granular lanes, numerical simulations suggest that granular lanes can transport recycled magnetic flux from nearby intergranular lanes, regardless of the underlying magnetic flux (Fischer et al. 2020). While the significance of this recycled flux in magnetic flux emergence remains uncertain – as it could potentially occur in any granule – our findings and those of Moreno-Insertis et al. (2018) offer insights. Firstly, sheets associated with granular lanes exhibit significantly lower mean horizontal flow divergence (left tail of the distribution in Fig. 13 and values in Col. 9 of Table B.1) compared to other granular phenomena, indicating their proximity to intergranular regions. Secondly, Moreno-Insertis et al. (2018) observed sub-photospheric, strongly magnetised volumes in upflows, often located near downflow regions. These observations collectively suggest that recycled magnetic flux may accompany the emergence of magnetic flux sheets.

By studying this observational sample of emerging magnetic flux sheets in the photosphere, we gain insights into the subsurface structure of emerging magnetic fields. Similar to Martínez González & Bellot Rubio (2009) for Ω-loop type emergence, we identify cases of ‘double events’, two events occurring in nearby granules and almost simultaneously; for instance, the event shown in Figure 5 and the events ID 58352 and ID 58366 (see Tables B.1 and B.2). These simultaneous events may indicate a large-scale magnetic structure accumulating beneath the photosphere and emerging in sections. However, further research is needed to investigate the convective plasma conditions that facilitate the emergence of small-scale magnetic flux from the subsurface to the photosphere, particularly through numerical simulations.

In our future work, we shall analyse numerical simulations and compare their results to the observational analysis presented here. Additionally, we shall focus on specific events that emerge into the chromosphere, analysing their observational characteristics and comparing them to numerical simulations.


Acknowledgments

Hinode is a Japanese mission developed and launched by ISAS/JAXA, collaborating with NAOJ as a domestic partner, NASA and STFC (UK) as international partners. Scientific operation of the Hinode mission is conducted by the Hinode science team organised at ISAS/JAXA. This team mainly consists of scientists from institutes in the partner countries. Support for the post-launch operation is provided by JAXA and NAOJ (Japan), STFC (U.K.), NASA, ESA, and NSC (Norway). 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 2021-00169). This project has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No. 824135, the Trans-National Access Programme of SOLARNET. This research has been supported as well by the European Research Council (ERC) through the Synergy grant No. 810218 (“The Whole Sun”, ERC-2018-SyG) and by the Spanish Ministry of Science, Innovation and Universities through project PGC2018-095832-B-I00. We thank Philip Lindner and Anjali J. Kaithakkal for recording and providing us with this dataset from their observational campaign in April 2019. We thank Oleksii Andriienko for performing the data reconstruction with SSTRED. We thank Serena Criscuoli for providing the SST dataset recorded in August 2011 and Shahin Jafarzadeh and Luc Rouppe van der Voort for reducing this dataset. We thank Vigeesh Gangadharan for providing us with a time series of snapshots of his CO5BOLD simulation and for his fruitful discussion and clarification regarding granular lanes. We thank Nazareth Bello-González for her comments and suggestions for improving our manuscript. S.M.D.C. and C.E.F. were funded by the Leibniz Association grant for the SAW-2018-KIS-2-QUEST project. S.L.G. acknowledge support from INAF. The National Solar Observatory is operated by the Association of Universities for Research in Astronomy, Inc. (AURA) under a cooperative agreement with the National Science Foundation. This research has made use of NASA’s Astrophysics Data System Bibliographic Services. 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), PyTorch (pytorch.org) and TensorFlow (tensorflow.org).

References

  1. Asensio Ramos, A., Requerey, I. S., & Vitas, N. 2017, A&A, 604, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Barthol, P., Gandorfer, A., Solanki, S. K., et al. 2011, Sol. Phys., 268, 1 [Google Scholar]
  3. Bellot Rubio, L., & Orozco Suárez, D. 2019, Liv. Rev. Sol. Phys., 16, 1 [Google Scholar]
  4. Bray, R., Loughhead, R., & Durrant, C. 1984, The Solar Granulation (Cambridge: Cambridge University Press) [Google Scholar]
  5. Centeno, R. 2012, ApJ, 759, 72 [NASA ADS] [CrossRef] [Google Scholar]
  6. Centeno, R., Socas-Navarro, H., Lites, B., et al. 2007, ApJ, 666, L137 [CrossRef] [Google Scholar]
  7. Centeno, R., Blanco Rodríguez, J., Del Toro Iniesta, J. C., et al. 2017, ApJS, 229, 3 [NASA ADS] [CrossRef] [Google Scholar]
  8. De Pontieu, B. 2002, ApJ, 569, 474 [CrossRef] [Google Scholar]
  9. Díaz-Castillo, S. M., Fischer, C. E., Rezaei, R., Steiner, O., & Berdyugina, S. 2024, A&A, 691, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Ellwarth, M., Fischer, C. E., Vitas, N., Schmiz, S., & Schmidt, W. 2021, A&A, 653, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Fisher, G. H., & Welsch, B. T. 2008, ASP Conf. Ser., 383, 373 [NASA ADS] [Google Scholar]
  12. Fischer, C. E., Bello González, N., & Rezaei, R. 2017, A&A, 602, L12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Fischer, C. E., Borrero, J. M., Bello González, N., & Kaithakkal, A. J. 2019, A&A, 622, L12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Fischer, C. E., Vigeesh, G., Lindner, P., et al. 2020, ApJ, 903, L10 [NASA ADS] [CrossRef] [Google Scholar]
  15. Gošić, M. 2015, Ph.D. Thesis, University of Granada, Spain [Google Scholar]
  16. Gošić, M., Bellot Rubio, L. R., Orozco Suárez, D., Katsukawa, Y., & del Toro Iniesta, J. C. 2014, ApJ, 797, 49 [Google Scholar]
  17. Gošić, M., Rubio, L. R. B., Cheung, M. C. M., et al. 2022, ApJ, 925, 188 [CrossRef] [Google Scholar]
  18. Gömöry, P., Beck, C., Balthasar, H., et al. 2010, A&A, 511, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Guglielmino, S. L., Martínez Pillet, V., Ruiz Cobo, B., et al. 2020, ApJ, 896, 62 [NASA ADS] [CrossRef] [Google Scholar]
  20. Hagenaar, H. J., Schrijver, C. J., & Title, A. M. 2003, ApJ, 584, 1107 [NASA ADS] [CrossRef] [Google Scholar]
  21. Hirzberger, J., Koschinsky, M., Kneer, F., & Ritter, C. 2001, A&A, 367, 1011 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Ishikawa, R., Tsuneta, S., & Jurčák, J. 2010, ApJ, 713, 1310 [CrossRef] [Google Scholar]
  23. Jin, C., Wang, J., & Zhou, G. 2009, ApJ, 697, 693 [CrossRef] [Google Scholar]
  24. Jin, C. L., Wang, J. X., Song, Q., & Zhao, H. 2011, ApJ, 731, 37 [Google Scholar]
  25. Kingma, D. P., & Ba, J. 2014, ArXiv e-prints [arXiv:1412.6980] [Google Scholar]
  26. Kosugi, T., Matsuzaki, K., Sakao, T., et al. 2007, Sol. Phys., 243, 3 [Google Scholar]
  27. Ledvina, V. E., Kazachenko, M. D., Criscuoli, S., et al. 2022, ApJ, 934, 38 [CrossRef] [Google Scholar]
  28. Lites, B. W., Kubo, M., Socas-Navarro, H., et al. 2008, ApJ, 672, 1237 [Google Scholar]
  29. Lites, B., Akin, D., Card, G., et al. 2013, Sol. Phys., 283, 579 [NASA ADS] [CrossRef] [Google Scholar]
  30. Löfdahl, M. G. 2002, SPIE Conf. Ser., 4792, 146 [Google Scholar]
  31. Löfdahl, M. G., Hillberg, T., de la Cruz Rodriguez, J., et al. 2018, ArXiv e-prints [arXiv:1804.03030] [Google Scholar]
  32. Martínez González, M. J., & Bellot Rubio, L. R. 2009, ApJ, 700, 1391 [CrossRef] [Google Scholar]
  33. Martínez González, M. J., Collados, M., Ruiz Cobo, B., & Solanki, S. K. 2007, A&A, 469, L39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Martínez Pillet, V., del Toro Iniesta, J. C., Álvarez-Herrero, A., et al. 2011, Sol. Phys., 268, 57 [CrossRef] [Google Scholar]
  35. Moreno-Insertis, F., Martinez-Sykora, J., Hansteen, V. H., & Muñoz, D. 2018, ApJ, 859, L26 [Google Scholar]
  36. Nordlund, Å., Stein, R. F., & Asplund, M. 2009, Liv. Rev. Sol. Phys., 6, 2 [Google Scholar]
  37. November, L. J., & Simon, G. W. 1988, ApJ, 333, 427 [Google Scholar]
  38. Palacios, J., Blanco Rodríguez, J., Vargas Domínguez, S., et al. 2012, A&A, 537, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Rast, M. P. 1995, ApJ, 443, 863 [NASA ADS] [CrossRef] [Google Scholar]
  40. Rempel, M. 2018, ApJ, 859, 161 [Google Scholar]
  41. Rempel, M., Bhatia, T., Bellot Rubio, L., & Korpi-Lagg, M. J. 2023, Space Sci. Rev., 219, 36 [NASA ADS] [CrossRef] [Google Scholar]
  42. Roudier, T., Eibe, M. T., Malherbe, J. M., et al. 2001, A&A, 368, 652 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Scharmer, G. B., Bjelksjo, K., Korhonen, T. K., Lindberg, B., & Petterson, B. 2003, SPIE Conf. Ser., 4853, 341 [NASA ADS] [Google Scholar]
  44. Scharmer, G. B., Narayan, G., Hillberg, T., et al. 2008, ApJ, 689, L69 [Google Scholar]
  45. Simon, G. W., Title, A. M., & Weiss, N. O. 2001, ApJ, 561, 427 [CrossRef] [Google Scholar]
  46. Smitha, H. N., Anusha, L. S., Solanki, S. K., & Riethmüller, T. L. 2017, ApJS, 229, 17 [Google Scholar]
  47. Stangalini, M., Giannattasio, F., & Jafarzadeh, S. 2015, A&A, 577, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Stangalini, M., Giannattasio, F., Erdélyi, R., et al. 2017, ApJ, 840, 19 [NASA ADS] [CrossRef] [Google Scholar]
  49. Steiner, O., Franz, M., Bello González, N., et al. 2010, ApJ, 723, L180 [Google Scholar]
  50. Stix, M. 2002, The Sun: An Introduction (Berlin: Springer-Verlag) [Google Scholar]
  51. Thornton, L. M., & Parnell, C. E. 2011, Sol. Phys., 269, 13 [Google Scholar]
  52. Tremblay, B., Roudier, T., Rieutord, M., & Vincent, A. 2018, Sol. Phys., 293, 57 [NASA ADS] [CrossRef] [Google Scholar]
  53. Tsuneta, S., Ichimoto, K., Katsukawa, Y., et al. 2008, Sol. Phys., 249, 167 [Google Scholar]
  54. van Noort, M., Rouppe van der Voort, L., & Löfdahl, M. G. 2005, Sol. Phys., 228, 191 [Google Scholar]
  55. Vigeesh, G., Jackiewicz, J., & Steiner, O. 2017, ApJ, 835, 148 [NASA ADS] [CrossRef] [Google Scholar]
  56. Wang, J., Liu, C., Cao, W., & Wang, H. 2020, ApJ, 900, 84 [NASA ADS] [CrossRef] [Google Scholar]
  57. Welsch, B. T., & Longcope, D. W. 2003, ApJ, 588, 620 [NASA ADS] [CrossRef] [Google Scholar]
  58. Zhang, J., Yang, S.-H., & Jin, C.-L. 2009, RAA, 9, 921 [NASA ADS] [Google Scholar]
  59. Zhou, G., Wang, J., & Jin, C. 2013, Sol. Phys., 283, 273 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Dataset details

This section includes the characteristics of the datasets used for the present analysis. Context data information is included for completeness but is not used in this work.

Table A.1.

Database Hinode

Table A.2.

Database SST

Appendix B: Magnetic flux sheet emergence events

Table B.1.

Magnetic flux sheet properties

Table B.2.

Host granule properties

All Tables

Table A.1.

Database Hinode

Table A.2.

Database SST

Table B.1.

Magnetic flux sheet properties

Table B.2.

Host granule properties

All Figures

thumbnail Fig. 1.

Representative snapshots of two distinct cases of emerging magnetic flux sheets. Left: Hinode SP full-sized and reconstructed image from slit data in the continuum of the Fe I 630.2 nm line pair with a FoV of 4 . 5 × 41 $ 4{{\overset{\prime\prime}{.}}}5 \times 41{{\prime\prime}} $ recorded on October 22 2007 at 18:34 UT. Right: Image of the CRISP instrument in the continuum of the Fe I 630.2 nm line recorded on August 06 2011. The FoV is about 41″ × 41″ cut from the original size of 57″ × 57″ to fit with the Hinode SP image. The scale is presented in megametres [Mm]. Yellow boxes mark the location of a detected event in each particular dataset: event ID 523 for the Hinode SP image and event ID 53961 for the CRISP image, respectively. Green, red, and blue contours indicate the apparent transverse magnetic flux, the negative longitudinal magnetic flux, and the positive longitudinal magnetic flux, respectively.

In the text
thumbnail Fig. 2.

Magnetic flux sheet event with ID 352 recorded by Hinode SP on September 26 2007 at 08:15 UT. From left to right, the first panel shows the continuum intensity in the background with contours. Pine green contours mark a threshold of 130 Mx cm−2 for the apparent transverse magnetic flux density patch ϕTapp. The blue (red) contours mark the positive +ϕLapp (negative −ϕLapp) polarity of the longitudinal magnetic flux density with thresholds of 20 Mx cm−2, 40 Mx cm−2 and 60 Mx cm−2. The black contours mark a single closed contour of zero-divergence defined as the boundary of the host granule. The second panel shows the divergence of the apparent horizontal velocity field marking the corresponding location of the sheet represented by the transverse magnetic flux density patch ϕTapp (pine green contour) and the boundary of the host granule (black contour). Finally, the third, fourth, and fifth panels show the bisector velocity field for the three different depth positions: 70%, 35%, and core (0%) of the Fe I 630.15 nm line corresponding roughly to the low, middle, and high photosphere, respectively.

In the text
thumbnail Fig. 3.

Duration of each event (left panel), the apparent maximum area reached by the magnetic flux sheet (central panel) and the proportion of a host granule’s area covered by its associated magnetic flux sheet (right panel). The granulation type is coded by different colours, as follows: exploding granules (EG): yellow bars and lines; granular lanes (GL): blue bars and lines; nascent granules (NG): green bars and lines; and complex granulation (CG): red bars and lines. The CG is only included in panels a and b. All histograms include the kernel density estimate (KDE) providing complementary information about the possible shape of the distribution. The actual values are included in Cols. 4−5 of Table B.1 and Col. 4 of Table B.2. (a) Event duration. (b) Magnetic flux sheet maximum area. (c) Host granule coverage.

In the text
thumbnail Fig. 4.

Magnetic flux sheet emergence event alongside the expansion of its host granule, which fragments as an exploding granule. The image series shows the continuum maps reconstructed from Hinode SP slit data. The times are indicated above the images. Contours mark the apparent transverse (pine green) and longitudinal (red contours for the negative polarity and blue contours for the positive polarity) magnetic flux density associated with the event. The magnetic flux density thresholds are defined as in Fig. 2. Black contours mark the host granule border. Colour maps within the host granule show the distribution of the bisector velocity at 70% of the line depth.

In the text
thumbnail Fig. 5.

Second case of magnetic flux sheet emergence in an exploding granule. The image series shows the continuum close to the Fe I 617.3 nm line recorded by CRISP. The plot set-up is defined as in Fig. 4. The blue (red) contours mark the positive (negative) polarity of the longitudinal magnetic flux density with thresholds of 40 Mx cm−2, 60 Mx cm−2 and 80 Mx cm−2. Black contours mark the host granule border. A secondary magnetic flux sheet, depicted by orange contours, was detected near and concurrently with the event within the exploding granule. Colour maps within the host granule show the distribution of the bisector velocity at 60% of the line depth.

In the text
thumbnail Fig. 6.

Case of a magnetic flux sheet emergence in a nascent and relatively small granule. The image series shows the continuum maps constructed from Hinode SP slit data. The plot set-up and information are as shown in Fig. 4.

In the text
thumbnail Fig. 7.

Second case of magnetic flux sheet emergence in a nascent granule. The image series shows the continuum close to the Fe I 617.3 nm line recorded by CRISP. The plot set-up and information are as shown in Fig. 5.

In the text
thumbnail Fig. 8.

Case of a magnetic flux sheet emergence associated with a granular lane. The image series shows the continuum maps constructed from Hinode SP slit data. The plot set-up and information are as is shown in Fig. 4. The dashed yellow curves in t = 69.6 min and t = 70.8 min highlight the granular lane signature.

In the text
thumbnail Fig. 9.

Case of a magnetic flux sheet emergence associated with a granule with granular lanes. The image series shows the continuum close to the Fe I 617.3 nm line recorded by CRISP. The plot set-up and information are as is shown in Fig. 5. The dashed yellow curve in t = 49.0 min highlights a granular lane signature. The yellow arrow marks the direction of displacement of the magnetic flux sheet.

In the text
thumbnail Fig. 10.

Total longitudinal magnetic flux measured in the footpoints of the magnetic flux sheet events. The histogram includes the kernel density estimate providing complementary information about the shape of the distribution.

In the text
thumbnail Fig. 11.

Scatter plots comparing the total longitudinal magnetic flux of the sheet, ΦLapp, to the apparent maximum area of the sheet Amax (left panel), the rate of change of the magnetic flux sheet area ΔA (central panel), and the divergence of the horizontal flow field within the sheet ∇h ⋅ vh(x, y)τ = 1 (right panel). The Pearson correlation coefficient is shown in each panel. See Table B.1 (Cols. 5, 6, 9) for numerical values.

In the text
thumbnail Fig. 12.

Average vertical velocity at different heights of the photosphere of the host granules during the emergence of their associated magnetic flux sheet. The box plots display the full distribution (scatter points), maximum and minimum values (whiskers lines), the interquartile range (box), upper and lower quartiles (box boundaries), and the median value (line within the box) for each height: low photosphere (blue box), middle photosphere (white box), and high photosphere (red box).

In the text
thumbnail Fig. 13.

Divergence of the horizontal plasma flow co-located with the magnetic flux sheet area. The horizontal velocity flow is estimated at τ = 1. The histogram includes the kernel density estimate providing complementary information about the shape of the distribution.

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.