Open Access
Issue
A&A
Volume 686, June 2024
Article Number A112
Number of page(s) 9
Section The Sun and the Heliosphere
DOI https://doi.org/10.1051/0004-6361/202348764
Published online 03 June 2024

© The Authors 2024

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

As discussed by Bello González et al. (2019), the basic necessary conditions of a minimum required flux and the presence of strong horizontal fields for the formation of penumbrae are known in general. However, there are a number of aspects that are still not known and/or understood. Observational analyses of pores and sunspots imply that a minimum magnetic flux on the order of 1020 Mx is necessary for the penumbra to form (Zwaan 1987; Leka & Skumanich 1998). This number is consistent with the theoretical model by Rucklidge et al. (1995), who investigated the stability of magnetic flux tubes and their properties based on their total flux. These models and observations describe pores and sunspots created by rather symmetrical magnetic flux tubes. As the magnetic flux within the flux tube increases, the magnetic field inclination on its edges increases, and the horizontal component of the magnetic field is strong enough to shape the rising convective cells to form the elongated penumbral filaments. The magnetohydrodynamic (MHD) simulations of sunspots by Rempel et al. (2009) as well as analysis of the observations by Jurčák et al. (2014a) suggest that the critical magnetic field inclination necessary for penumbra formation is around 45°. This value is also in agreement with the simpler models of Rucklidge et al. (1995).

Another property of penumbrae in sunspots that might play a role in sustaining horizontal fields resulting in penumbra formation is the presence of an overlaying magnetic field that traps the emerging horizontal field in the photosphere in the form of penumbral filaments. Solanki & Montavon (1993) proposed the uncombed structure of the magnetic field in sunspot penumbra where horizontal flux tubes or tops of elongated convective cells (as suggested by MHD simulations) are interleaved by a stronger and more vertical magnetic field known as the background magnetic field or intraspines. The actual configuration of the magnetic field is very similar to the proposed model, as confirmed by Tiwari et al. (2013), and it was shown by Borrero et al. (2007, 2008) that the background magnetic field closes above the penumbral filaments.

Another form of penumbra is the orphan penumbra. These are penumbrae not coupled to an umbra. Jurčák et al. (2014a) reported no background magnetic field present in these structures, while Zuccarello et al. (2014) found indications of it. Orphan penumbrae are in most aspects (e.g. Evershed flows and penumbral grains) similar to sunspot penumbrae (Jurčák et al. 2014a). Lim et al. (2013) and Zuccarello et al. (2014) found orphan penumbrae originating from an emerging magnetic field trapped at the photospheric level by the pre-existing overlying fields in areas of complex topology in active regions, typically near the polarity inversion lines.

Details of penumbra formation were studied by Schlichenmaier et al. (2010a,b, 2011, 2012), who used unique observations of sunspot formation over a period of 4.5 h. The formation of penumbra was in this case caused by the accumulation of magnetic flux on one side of the protospot where a stable segment of a penumbra formed on the opposite side. Prior to the formation of a stable penumbra, Schlichenmaier et al. (2011) observed temporal filaments with what they introduced as counter-Evershed flows. Using the same data, Rezaei et al. (2012) found signatures of a low lying (observed in the Fe I 617.3 nm line) magnetic canopy in the regions where penumbra later formed. This scenario was subsequently confirmed by the numerical simulations of MacTaggart et al. (2016), who studied the flux emergence of a twisted magnetic flux tube and the magnetic canopy formed preferentially away from the emergence site. Indication of a chromospheric magnetic canopy was found also by Shimizu et al. (2012) and Romano et al. (2013), who observed a chromospheric halo in regions where penumbra later formed. The presence of a chromospheric halo preceding the formation of the penumbra was interpreted by Romano et al. (2013, 2014), and Murabito et al. (2016) as a necessary configuration of the magnetic field that allows for downward bending of the chromospheric field lines and the creation of a shallow loop that is also associated with temporal counter-Evershed flows caused by siphon flows. The recent work by Lindner et al. (2023) favours a scenario in which emerging flux blocked by an overlying (chromospheric) canopy provides the strong horizontal fields necessary for penumbra formation.

Li et al. (2018) investigated the formation of a penumbra sector and found that in a first stage, the area and magnetic flux of the forming penumbral sector and the umbra both increased and that in a subsequent second stage, the area and magnetic flux of the penumbra increased at the expense of the umbra area and the flux. They also found a persistent blueshift in the penumbra formation area that they interpreted as a signature of magnetic flux emergence. However, they did not find any chromospheric counterpart coupled to the penumbra formation unlike the annular region reported by Shimizu et al. (2012). Yet, they found a correlation between the penumbral sector formation and sunspot rotation. Li et al. (2019) also found indications of penumbral areas formed by emerging flux trapped in the photosphere by the overlying pre-existing field in a complex active region.

Jurčák (2011) discovered the magnetic property defining the umbra-penumbra boundary of stable sunspots, a critical value of the vertical magnetic field (1.86 kG, Jurčák et al. 2018). Jurčák et al. (2015) found this canonical value to be a natural inner boundary for the expansion of penumbral magneto-convection in umbral areas. Their result is in agreement with the findings of Li et al. (2018), as umbral areas are originally occupied by developing penumbra. Therefore, penumbra can only set in and operate in regions with a vertical field below the critical value. The observation of such a vertical magnetic field boundary against vigorous modes of magneto-convection (e.g. the penumbral mode) is supported by the analysis of the Gough and Tayler criterion (Gough & Tayler 1966) of stability by vertical fields against overturning convection by Schmassmann et al. (2021).

From a modelling point of view, there are a number of realistic simulations of sunspots with penumbra (see, e.g., Rempel et al. 2009; Rempel 2011, 2012; Chen et al. 2017; Panja et al. 2021). Yet only a few works address the onset and formation of penumbra. In accordance with observations, Hurlburt & Rucklidge (2000), Simon & Weiss (1970), and Jahn & Schmidt (1994) found that an increasing inclination (with respect to the normal to the solar surface) of the magnetic field and an increasing flux in sunspot models are key factors for the origin of penumbra. Wentzel (1992) proposed the “fallen magnetic flux tubes” model for penumbra formation. There, the upwelling of a mass flow in the inner footpoint of the field lines within the umbra provokes the flux tubes to fall onto the photosphere and become submerged in the surroundings of the sunspot (i.e. the penumbra forms from top to bottom). On local scales, Rempel (2011) presented a more elaborated model where the flux tubes bend due to the mass loads originated from pressure-driven upflows (see Fig. 19 of the cited work). On the other hand, Chen et al. (2017) simulated the formation of an asymmetric sunspot formed by the coalescence of magnetic features with the same polarity from small-scale dipoles embedded in granular cells (i.e. the penumbra forms from bottom to top). Regarding penumbral radial flows (e.g. Evershed flow or counter-Evershed flow), Chen et al. (2017) found that penumbral filaments were dominated by inflows in the region where the flux emergence was more active (between the two main spots). More interestingly, some penumbral filaments developed intermittent outflows even when they were originally dominated and surrounded by inflows. Transient inflows embedded in outflows were obtained by Siu-Tapia et al. (2018) from simulations based on the MURaM radiative MHD code.

In this work, we describe the observational properties of magneto-convection at the outskirts of a protospot before the onset of a penumbra and during the first stages of the penumbral formation process. The observations and data analysis are described in Sect. 2. Three distinct penumbral formation processes are respectively presented and discussed in Sects. 3 and 4. This work closes with conclusions obtained from the results.

2. Observations and data analysis

The data were recorded on 9 July 2009 with the “Göttingen” Fabry-Pérot Interferometer (GFPI, Bendlin et al. 1992; Puschmann et al. 2006) in the full-Stokes polarimetric mode (Bello González et al. 2008) at the VTT Telescope (Tenerife). The observations consisted of a time series of 4:40 h following the evolution of the leading sunspot in NOAA 11024 located at a heliocentric angle θ ∼ 28°. Here, we make use of the imaging data in the G-band and full-Stokes polarimetry in the Fe I 617.3 nm photospheric line. We note that these data were previously used in other investigations (Schlichenmaier et al. 2010a,b; Rezaei et al. 2012; Bello González et al. 2012; Jurčák et al. 2014b), and therefore we defer to the referenced works for further details about the observations.

We restricted our analysis to a 2 h time span (08:32 UT–10:32 UT) of uninterrupted observation of penumbra formation consisting of 109 GFPI scans. We focused on selected regions of the protospot periphery in which the formation of penumbral segments can be followed over time. We ran Milne-Eddington inversions of the Fe I 617.3 nm spectropolarimetric observations with the VFISV code (Borrero et al. 2011), which provides the components of the magnetic field vector and line-of-sight (LOS) velocity maps. We used the code AMBIG (Leka et al. 2009) to solve the 180° ambiguity in the component of the magnetic field perpendicular to the observer’s LOS. We also made use of routines of the AZAM code (Lites et al. 1995) to transform the magnetic field vector from the LOS frame to the local reference frame (LRF). During the time series, the field of view varies. This resulted in a limited effective field of view.

3. Results

Figure 1 shows a snapshot of the field of view under analysis at the beginning of the observations. An animation of the development and evolution of the protospot can be found online. At this early stage, the protospot already shows segments of stable penumbra around (22″, 20″) and (16″, 11″). The penumbra around (11″, 19″) is highly distorted by the continuous emergence of the new flux in the developing active region (Rezaei et al. 2012). We note that the protospot is the leading polarity of NOAA 11024, and the following polarity is directed towards the east (not visible within the field of view).

thumbnail Fig. 1.

Snapshot of NOOA 11024 sunspot at the early stages of penumbra development (08:40:52 UT) as seen in G-band intensity (upper left), LOS velocity (upper right), total magnetic field strength (lower left), and magnetic field inclination in the LRF (lower right). The pink iso-contour outlines the umbral boundary using a threshold of 50% of the quiet-Sun intensity. The three pink boxes in the G-band intensity maps mark the regions studied in detail. The maps have been rotated according to the solar north-south direction. A movie showing the temporal evolution is available online.

We studied the magnetic field and velocity properties in regions where penumbrae later form. We note that the Milne-Eddington inversion scheme does not allow for disentangling of the different physical properties of the spot background component from the penumbral filaments. Instead, the Milne-Eddington inversion provides mean values of LOS velocity, magnetic field strength, and inclination (Borrero et al. 2014) from regions where the 617.3 nm line forms and is most sensitive to changes of the physical properties. This implies that the magnetic canopy affects the plasma properties inferred at the external boundary of the umbral core and beyond. In this case, the influence is asymmetric due to the heliocentric angle and the resulting projection effects. As seen in the bottom maps of Fig. 1 (B and γLRF), the magnetic boundary oriented towards the disc centre is sharper than the magnetic boundary oriented towards the solar limb (i.e. the magnetic radial gradient at the umbral boundary is larger at the centre side than at the limb side). Figure 2 illustrates how the projection of the magnetopause on the quiet-Sun regions induces a larger magnetic gradient at the centre side of the umbral boundary, while it induces a smaller magnetic gradient at the limb side, which translates into a sharper magnetic boundary at the centre side of the spot. We note that this is a simplified sketch where the complex nature of the magnetopause is not shown and the line formation region, spot size, and Wilson depression are not to scale.

thumbnail Fig. 2.

Sketch of the projection effects on the derived magnetic properties around a spot at θ = 28°. The influence of the magnetopause (i.e. the most external red lines) on the derived magnetic properties reaches areas farther away from the umbral boundary at the limb side compared to the centre side.

In the following three subsections, we identify three distinct types of magneto-convection: a bi-directional flow in elongated granules located at the centre side of the spot, a transient filament with counter-Evershed flow, and a granular pattern located at the limb side. All of them are followed by the formation of the penumbra and the onset of the Evershed flow.

3.1. Bi-directional flow in elongated granules

On the centre-side region of the protospot, the penumbra starts to form at approximately 9:20 UT (i.e. some 50 min after the start of our observations). In Fig. 3, we show a 2 h time slice of a 5 pixel-wide selected region from 08:32 UT until 10:30 UT for intensity, inclination, field strength, and LOS velocity.

thumbnail Fig. 3.

Bi-directional flow in elongated granules. Top panels: time slices of selected penumbra-formation region on the sunspot centre side. From top to bottom: LOS velocities, magnetic field strength, inclination (LRF), and G-band intensity. Bottom panels: context maps containing the areas (marked by purple boxes) selected for the time slices as seen in the G-band intensity and LOS velocity at 8:35 UT (two left panels) and at 10:24 UT (two right panels). The violet arrow points the direction along the slice length. The purple and red iso-contours highlight velocities of ±130 m s−1, respectively.

The time slice shows that from the start of the observations and for about 50 min, the velocity in the selected region at the very outskirts of the umbra is characterised by a persistent redshift-blueshift pattern (i.e. the flow is directed towards the umbra in its close vicinity and outwards further out). This bi-directional velocity pattern occurs in regions with weak magnetic fields of less than 500 G. Only the innermost end of the redshifted region is rooted in regions with 1000 G field strength inside the umbra, as seen in the continuum intensity slices. The time slice shows that before the penumbra sets in, the (radial) red-blue pattern is associated with inclinations that go from positive polarity (blue-green) to horizontal (yellow) to opposite polarity (orange-red). The G-band intensity images show that the bi-directional flow pattern is associated with one intensity structure that resembles an elongated granule. The mean intensity of this granule is comparable to the quiet-Sun intensity.

We note that in the region under study, the LOS inclination of the magnetic field is around 30° and towards the centre of the solar disc. The strong magnetic field observed close to the umbral boundary is part of the magnetic body of the spot. Further away from the umbral boundary, the LOS does not cross the magnetic canopy of the protospot, as illustrated in Fig. 2 (centre side), and thus we observe the actual magnetic field configuration of the bi-directional flow pattern.

From about 9:10 UT onward, the red-blue patch is displaced away from the umbra, and the displacement increases in time. The reason for this displacement is the formation of a new convective cell featuring an outwards-oriented plasma motion (blueshifted region) that is initially most visible in the G-band continuum images as faint, bright filaments that do not have a strong enough LOS flow to be encircled by the red contour (i.e. −130 m s−1). This outwards flow can be considered the first appearance of the Evershed flow. As time progresses, the length of the region occupied by the incipient Evershed flow expands. In the continuum intensity slices, one can see multiple bright structures whose inner footpoint, which is also the brightest, moves inwards into the umbra in a manner similar to penumbral grains in stable penumbrae (e.g. Sobotka et al. 1999). From the intensity time slices it becomes obvious that the forming penumbra with Evershed flow does not heat the solar surface as effectively as the convection pattern characterised by the bi-directional flow. The mean G-band intensity in the forming penumbra decreases to 70% of the quiet-Sun intensity (0.7 IQS). This value is comparable to the average penumbral intensity measured in fully developed sunspot penumbrae.

Along with the expansion of the region occupied by the Evershed flow, we observed an expansion of the region with a strong magnetic field. The initially sharp gradient of the magnetic field strength at the protospot boundary became more diffuse. At the end of the analysed period, the magnetic field strength increased to approximately 600 G at the outer edge of the slice. A similar trend was also observed for the magnetic field inclination and can be explained by the expansion of the magnetic field that is associated with the penumbra formation. In the beginning, we had a sharp boundary between the nearly vertical magnetic field at the protospot boundary and the horizontal field in the weakly magnetised bi-directional flow patch. Later on, we found LRF inclinations around 80° at the end of the slice, and such values are comparable to those found at the outer penumbral boundaries of stable sunspots.

3.2. Transient filament with counter-Evershed flow

On the limb-side boundary of the protospot, we observed one clear example of the formation of a transient filament carrying a counter-Evershed flow. In Fig. 4, we show the time slice evolution of this region. We note that the LOS inclination in this region is around 70° and points towards the limb. Therefore, the magnetic boundary of the protospot is not as sharp as on the centre side, and our results are influenced by looking through the magnetic canopy of the protospot. Also, the viewing angle is such that upflows located at the penumbral bright grains are seen as blueshifts, and a regular Evershed flow is therefore not represented only by redshifts.

thumbnail Fig. 4.

Analogous to Fig. 3 but showing the limb-side region where a transient filament with a counter-Evershed flow formed. The bottom panels correspond to observations at 9:23 UT (two left panels) and at 10:01 UT (two right panels).

Prior to the formation of the transient filament, we observed within the slice an edge of the region with a formed penumbra. The innermost part of the penumbra had intensities around 0.7 IQS, but the second row of penumbral bright grains (located around 2″ of the slice length) reached intensities around IQS. As mentioned before, at this location, the bright grains are associated with blueshifts; these blue patches are visible until 8:52 UT. The magnetic boundary of the protospot was diffuse due to the presence of the penumbra, while the magnetic field strength and inclination steadily decreased and increased along the slice, respectively.

Around 9:18 UT, a filament with counter-Evershed flow developed. Within 10 min, it reached its full length and maximum LOS velocity and then decayed for approximately 15 min. It reached its maximum continuum intensity in later stages, around 9:34 UT. Its formation can also be identified in the time-slice plots of the magnetic field strength and inclination. The magnetic field is weaker and more horizontal than that of regular penumbra observed before and afterward. A similar behaviour has been reported for some of the cases in a statistical study of counter-Evershed flows by Castellanos Durán et al. (2021). After the disappearance of the transient filament, we observed a similar magnetic field configuration as before its formation. There are no strong LOS velocities observed, as the LOS inclination is close to 90° within the slice. We note that the Evershed flow observed before 9:12 UT is associated with LRF inclinations around 60°, and such inclinations are not found within the slice between minutes 9:42 UT and 10:02 UT. Starting at 9:52 UT, we observed weak blueshifts associated with bright penumbral grains of very short penumbral filaments that have LOS inclinations close to 90°. After 10:02 UT, the filaments became longer and more inclined, so the redshifted regions associated with Evershed flow appeared. As in the case of Sect. 3.1, we observed steady lengthening of the filaments carrying the Evershed flow in time.

3.3. Granular pattern

On a different segment of the limb-side boundary of the protospot, we analysed the granular pattern that precedes the formation of the Evershed flow. The evolution of this region is shown in Fig. 5. As in the case of Sect. 3.2, we were looking through the magnetic canopy of the protospot, and the physical parameters resulting from the inversion are influenced by this.

thumbnail Fig. 5.

Analogous to Fig. 3 but showing the limb-side region where the granular pattern was observed before the penumbra formation. The bottom panels correspond to observations at 9:00 UT (two left panels) and at 10:24 UT (two right panels).

At the displayed segment of the protospot boundary, in the continuum intensity maps we observed for the first 80 min granular patterns with granules smaller than those observed in the quiet-Sun regions. Despite being smaller, the mean intensity is comparable to the IQS. During this time, we observed blueshifts. The structure of the LOS motions does not correspond to the granular pattern seen in the continuum intensity. In the dopplergram, the motions have a filamentary structure, and their lifetimes are longer than those of granular upflows. Such flows would correspond to the counter-Evershed flow described in Sect. 3.2.

During the two studied hours, we did not detect any significant nor sudden changes in the configuration of the magnetic field along the slice in Fig. 5. The gradients of the magnetic field strength and inclination along the slice length steadily decreased in time. During the first 80 min until 9:52 UT, we did not detect any filaments with Evershed flows represented by redshifts. However, during the whole two hours, we observed short filaments on the protospot boundary that moved inwards and contained weak blueshifts. As explained in Sect. 3.2, these can be associated with upflows located in bright grains of such short filaments, where the outflow motion characterising the Evershed flow cannot be detected due to the viewing angle.

The filaments carrying the Evershed flow started to be apparent in the intensity slices around 09:44 UT as bright structures. We could not detect any LOS velocity due to the viewing angle. As the filaments lengthened in time, we started to detect redshifted LOS velocities around 09:52 UT in parts of the filaments where the flow is horizontal and oriented towards the limb side. It is noteworthy that the Evershed flow sets in rapidly, within five to ten minutes. Until the end of the analysed period, the filaments lengthened and the Evershed flow became stronger.

4. Discussion and conclusions

We have identified three distinct types of magneto-convection-related motions prior to penumbra formation. These have been reported previously in other solar contexts unrelated to penumbra formation.

The bi-directional flows observed on the centre side and contained in one elongated granule, as seen in the intensity maps, resemble in many aspects the flux emergence event described in detail by Centeno et al. (2017). In the case of the flux emergence described by Centeno et al. (2017), the transversal magnetic field strength reached a maximum of 800 G. In our case, the horizontal field component is, at maximum, around 200 G. Higher field strength values can be found at the footpoints of the flow pattern, but the inner footpoint is influenced by the magnetic canopy of the protospot (see Fig. 2), and the outer footpoint is situated in the vicinity of another patch of positive polarity. In our case, the bright elongated granule is not as long and thin as in the case described by Centeno et al. (2017). Based on this comparison and on the observations of orphan penumbra, where the horizontal field component is even higher (above 1000 G; Jurčák et al. 2014a) and where intensity structures are visually comparable to penumbral filaments, we surmise that the horizontal component of the magnetic field determines the width and length of the convective cell.

Counter-Evershed flows have already been reported on the boundaries of the studied protospot by Schlichenmaier et al. (2011). There are other studies that have described observations of counter-Evershed flows in penumbra since then (Louis et al. 2014; Siu-Tapia et al. 2017). The recent statistical analysis of Castellanos Durán et al. (2021) shows that counter-Evershed flows are common phenomena in the majority of active regions, although they are more frequently observed when the filaments carrying the counter-Evershed flow are related to sunspot light bridges (see also Kleint & Sainz Dalda 2013; Louis et al. 2020). However, the typical lifetimes found by Castellanos Durán et al. (2021) are significantly longer (most commonly around 7 h) than the case studied in Sect. 3.2. Similar lifetimes in counter-Evershed flows have been found in sunspot MHD simulations (i.e. of a few hours) by Siu-Tapia et al. (2018). The big disparity between the counter-Evershed flow found in Sect. 3.2 that lasts only approximately 30 min and the observed and simulated inflows that last a few hours suggests that the penumbral formation process may benefit from the existence of short-lived counter-Evershed flows. Based on a more detailed analysis of the counter-Evershed flows reported by Schlichenmaier et al. (2011) at the centre side of the spot, we found that they actually correspond to the inner end of a bi-directional flow pattern.

It is well known that the magnetic boundaries of pores exceed those defined by intensity (Keppens & Martinez Pillet 1996) and that the size of the granules is influenced by the magnetic field strength (Hirzberger et al. 2002; Lagg et al. 2014; Jurčák et al. 2017; Vigeesh et al. 2017). Therefore, it is common to observe small granules near the boundaries of pores and protospots. Moreover, Falco et al. (2017) showed that the effects of strong magnetic fields on inhibiting the convection are visible in the granules observed in light bridges and in umbral dots, where granular sizes are smaller than those reported in quiet-Sun regions and their mean continuum intensity is darker than in quiet-Sun granules. Yet granules spanning from small to quiet-Sun-like sizes are observed in plage regions with mean magnetic fields weaker than 700 G and quiet-Sun-like mean continuum intensities. As shown in Jurčák et al. (2017), the intensity of granules is not influenced significantly by the magnetic field strength. We observed small but bright granules in regions with magnetic field strengths from 200 G to 700 G, depending on the distance from the protospot boundary, and with LRF inclinations of around 70° (the retrieved values are influenced by the LOS crossing the magnetic field canopy of the protospot; see Fig. 2). Thus, the horizontal component is not strong enough (yet) to shape the granules into more elongated convective cells. However, LOS velocity maps of the same regions show a filamentary blueshifted structure that does not match a typical granular motion. The disagreement between the granular structure and the filamentary flow is possibly caused by the height difference between the formation regions of the continuum intensity and the Fe I line. The lack of significant or sudden changes in the configuration of the magnetic field during the convective transformation from small granules to penumbral filament might be a combination of the effects of the LOS crossing the magnetic field canopy and the viewing angle.

Irrespective of the location where we observed the initial phases of the penumbra formation, the process always follows a similar pattern. We note that due to the viewing angle and LOS effects, the situation is best seen in Fig. 3, where the LOS is well aligned with the magnetic field on the umbral boundary, and we thus detected very weak flows along magnetic field lines. On the other hand, on the limb side, the LOS is close to being perpendicular to the magnetic field orientation, and the LOS velocities are thus hard to detect.

The process of penumbra formation starts within umbral areas with vertical fields below the critical threshold (Jurčák et al. 2015) and the appearance of short, bright filaments. Schmassmann et al. (2021) has described the role of vertical fields in the inhibition of vigorous modes of magneto-convection in umbrae. These filaments carry weak flows oriented outwards, and the bright penumbral grains (heads) of these filaments move further into the umbra. As time progresses, the newly appearing filaments become longer and extend radially from the original location of the umbra–quiet-Sun boundary. Thus the sunspot expands to become larger in size. To some extent, the umbra-penumbra boundary also moves inwards to the umbra, as explained in Jurčák et al. (2015). As the penumbra forms, the radial gradient of the magnetic field strength and inclination decreases, eventually reaching values typically observed in stable sunspots (i.e. field strength and inclination around 2500 G and 30°, respectively, on the inner penumbral boundary and around 600 G and 70° on the outer penumbral boundary).

The initial appearance of the penumbral filaments can be explained by the MHD simulation of Rempel (2011), which describes how the mass loading along inclined field lines bends the magnetic field and creates filaments. The length of the filaments is given by the magnetic field inclination, so the more inclined the field initially is, the longer are the filaments that can be produced. However, it has not yet been observed if the actual process of filament formation is identical to this scenario. It is also possible that the short filaments are created by small-scale flux emergence in which the horizontal field is carried up to the photosphere by convective motions, as obtained with a numerical model of a forming sunspot by Chen et al. (2017). Yet, these authors found that most of the penumbral filaments form in the flux emergence region, while we find that the penumbral filaments mostly form in the other side of the spot. The main difference between these scenarios is that in the simulations, the penumbral filaments form from the pre-existing field, while in the second scenario, the penumbral filaments form due to additional flux emergence. Only in the case of the first scenario would we observe the re-emergence of the bent field lines (sea-serpent configuration). Given the complex configuration of the sunspot magnetic field, we cannot be sure if such a sea-serpent structure is indeed related to the re-emergence of the field lines that were bent by the mass loading or if it is just an emerging flux loop and part of the so-called background magnetic field that is observed around filaments everywhere in sunspot penumbra. We note that we are comparing the observations of forming filaments in a developing penumbra to the Rempel (2011) MHD simulation of filaments in a stable penumbra.

According to Rempel (2011), Evershed flows develop from the deflection of the (convective) upflow into a radial outflow by the action of Lorentz forces in a thin layer right below the τ = 1 surface. We surmise that in the observations presented here, the bi-directional flows result from the elongation of (magneto-) convective cells by horizontal fields building up at the outskirts of the protospot. At this stage, the Lorentz force is not strong enough to drive an Evershed flow. During the following flux accumulation process (Rezaei et al. 2012), the horizontal fields develop the necessary gradient (with height) to develop a Lorentz force that eventually deflects the emerging upflows into an Evershed flow.

We find that the Evershed flow onsets in the umbra, and as it develops, it expands outwards into actual penumbral filaments. Physical conditions allowing for bi-directional flows and their evolution into Evershed flows will be investigated in future sunspot simulations.

Movie

Movie 1 associated with Fig. 1 (global_evolution_review1) Access here

Acknowledgments

This work was supported by the Czech-German common grant, funded by the Czech Science Foundation under the project 23-07633K and by the Deutsche Forschungsgemeinschaft under the project BE 5771/3-1 (eBer-23-13412), and the institutional support ASU:67985815 of the Czech Academy of Sciences. The Vacuum Tower Telescope is operated by the Leibniz-Institut für Sonnenphysik (KIS) in Freiburg, at the Spanish Observatorio del Teide of the Instituto de Astrofísica de Canarias.

References

  1. Bello González, N., Okunev, O., & Kneer, F. 2008, A&A, 490, L23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Bello González, N., Kneer, F., & Schlichenmaier, R. 2012, A&A, 538, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Bello González, N., Jurčák, J., Schlichenmaier, R., & Rezaei, R. 2019, in Solar Polariation Workshop 8, eds. L. Belluzzi, R. Casini, M. Romoli, & J. Trujillo Bueno, ASP Conf. Ser., 526, 261 [Google Scholar]
  4. Bendlin, C., Volkmer, R., & Kneer, F. 1992, A&A, 257, 817 [NASA ADS] [Google Scholar]
  5. Borrero, J. M., Bellot Rubio, L. R., & Müller, D. A. N. 2007, ApJ, 666, L133 [NASA ADS] [CrossRef] [Google Scholar]
  6. Borrero, J. M., Lites, B. W., & Solanki, S. K. 2008, A&A, 481, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Borrero, J. M., Tomczyk, S., Kubo, M., et al. 2011, Sol. Phys., 273, 267 [Google Scholar]
  8. Borrero, J. M., Lites, B. W., Lagg, A., Rezaei, R., & Rempel, M. 2014, A&A, 572, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Castellanos Durán, J. S., Lagg, A., & Solanki, S. K. 2021, A&A, 651, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Centeno, R., Blanco Rodríguez, J., Del Toro Iniesta, J. C., et al. 2017, ApJS, 229, 3 [NASA ADS] [CrossRef] [Google Scholar]
  11. Chen, F., Rempel, M., & Fan, Y. 2017, ApJ, 846, 149 [Google Scholar]
  12. Falco, M., Puglisi, G., Guglielmino, S. L., et al. 2017, A&A, 605, A87 [EDP Sciences] [Google Scholar]
  13. Gough, D. O., & Tayler, R. J. 1966, MNRAS, 133, 85 [Google Scholar]
  14. Hirzberger, J., Bonet, J. A., Sobotka, M., Vázquez, M., & Hanslmeier, A. 2002, A&A, 383, 275 [EDP Sciences] [Google Scholar]
  15. Hurlburt, N. E., & Rucklidge, A. M. 2000, MNRAS, 314, 793 [NASA ADS] [CrossRef] [Google Scholar]
  16. Jahn, K., & Schmidt, H. U. 1994, A&A, 290, 295 [NASA ADS] [Google Scholar]
  17. Jurčák, J. 2011, A&A, 531, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Jurčák, J., Bellot Rubio, L. R., & Sobotka, M. 2014a, A&A, 564, A91 [Google Scholar]
  19. Jurčák, J., Bello González, N., Schlichenmaier, R., & Rezaei, R. 2014b, PASJ, 66, S3 [CrossRef] [Google Scholar]
  20. Jurčák, J., Bello González, N., Schlichenmaier, R., & Rezaei, R. 2015, A&A, 580, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Jurčák, J., Lemmerer, B., & van Noort, M. 2017, in Fine Structure and Dynamics of the Solar Atmosphere, eds. S. Vargas Domínguez, A. G. Kosovichev, P. Antolin, & L. Harra, 327, 34 [Google Scholar]
  22. Jurčák, J., Rezaei, R., González, N. B., Schlichenmaier, R., & Vomlel, J. 2018, A&A, 611, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Keppens, R., & Martinez Pillet, V. 1996, A&A, 316, 229 [NASA ADS] [Google Scholar]
  24. Kleint, L., & Sainz Dalda, A. 2013, ApJ, 770, 74 [Google Scholar]
  25. Lagg, A., Solanki, S. K., van Noort, M., & Danilovic, S. 2014, A&A, 568, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Leka, K. D., & Skumanich, A. 1998, ApJ, 507, 454 [NASA ADS] [CrossRef] [Google Scholar]
  27. Leka, K. D., Barnes, G., Crouch, A. D., et al. 2009, Sol. Phys., 260, 83 [Google Scholar]
  28. Li, Q., Yan, X., Wang, J., et al. 2018, ApJ, 857, 21 [NASA ADS] [CrossRef] [Google Scholar]
  29. Li, Q., Yan, X., Wang, J., et al. 2019, ApJ, 886, 149 [NASA ADS] [CrossRef] [Google Scholar]
  30. Lim, E.-K., Yurchyshyn, V., Goode, P., & Cho, K.-S. 2013, ApJ, 769, L18 [NASA ADS] [CrossRef] [Google Scholar]
  31. Lindner, P., Kuckein, C., González Manrique, S. J., et al. 2023, A&A, 673, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Lites, B. W., Low, B. C., Martinez Pillet, V., et al. 1995, ApJ, 446, 877 [Google Scholar]
  33. Louis, R. E., Beck, C., Mathew, S. K., & Venkatakrishnan, P. 2014, A&A, 570, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Louis, R. E., Beck, C., & Choudhary, D. P. 2020, ApJ, 905, 153 [CrossRef] [Google Scholar]
  35. MacTaggart, D., Guglielmino, S. L., & Zuccarello, F. 2016, ApJ, 831, L4 [NASA ADS] [CrossRef] [Google Scholar]
  36. Murabito, M., Romano, P., Guglielmino, S. L., Zuccarello, F., & Solanki, S. K. 2016, ApJ, 825, 75 [NASA ADS] [CrossRef] [Google Scholar]
  37. Panja, M., Cameron, R. H., & Solanki, S. K. 2021, ApJ, 907, 102 [Google Scholar]
  38. Puschmann, K. G., Kneer, F., Seelemann, T., & Wittmann, A. D. 2006, A&A, 451, 1151 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Rempel, M. 2011, ApJ, 729, 5 [Google Scholar]
  40. Rempel, M. 2012, ApJ, 750, 62 [Google Scholar]
  41. Rempel, M., Schüssler, M., Cameron, R. H., & Knölker, M. 2009, Science, 325, 171 [CrossRef] [Google Scholar]
  42. Rezaei, R., Bello González, N., & Schlichenmaier, R. 2012, A&A, 537, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Romano, P., Frasca, D., Guglielmino, S. L., et al. 2013, ApJ, 771, L3 [NASA ADS] [CrossRef] [Google Scholar]
  44. Romano, P., Guglielmino, S. L., Cristaldi, A., et al. 2014, ApJ, 784, 10 [NASA ADS] [CrossRef] [Google Scholar]
  45. Rucklidge, A. M., Schmidt, H. U., & Weiss, N. O. 1995, MNRAS, 273, 491 [NASA ADS] [CrossRef] [Google Scholar]
  46. Schlichenmaier, R., Bello González, N., Rezaei, R., & Waldmann, T. A. 2010a, Astron. Nachr., 331, 563 [NASA ADS] [CrossRef] [Google Scholar]
  47. Schlichenmaier, R., Rezaei, R., Bello González, N., & Waldmann, T. A. 2010b, A&A, 512, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Schlichenmaier, R., González, N. B., & Rezaei, R. 2011, in Physics of Sun and Star Spots, eds. D. Prasad Choudhary, & K. G. Strassmeier, IAU Symp., 273, 134 [NASA ADS] [Google Scholar]
  49. Schlichenmaier, R., Rezaei, R., & González, N. B. 2012, in 4th Hinode Science Meeting: Unsolved Problems and Recent Insights, eds. L. Bellot Rubio, F. Reale, & M. Carlsson, ASP Conf. Ser., 455, 61 [Google Scholar]
  50. Schmassmann, M., Rempel, M., Bello González, N., Schlichenmaier, R., & Jurčák, J. 2021, A&A, 656, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Shimizu, T., Ichimoto, K., & Suematsu, Y. 2012, ApJ, 747, L18 [NASA ADS] [CrossRef] [Google Scholar]
  52. Simon, G. W., & Weiss, N. O. 1970, Sol. Phys., 13, 85 [Google Scholar]
  53. Siu-Tapia, A., Lagg, A., Solanki, S. K., van Noort, M., & Jurčák, J. 2017, A&A, 607, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Siu-Tapia, A. L., Rempel, M., Lagg, A., & Solanki, S. K. 2018, ApJ, 852, 66 [Google Scholar]
  55. Sobotka, M., Brandt, P. N., & Simon, G. W. 1999, A&A, 348, 621 [NASA ADS] [Google Scholar]
  56. Solanki, S. K., & Montavon, C. A. P. 1993, A&A, 275, 283 [NASA ADS] [Google Scholar]
  57. Tiwari, S. K., van Noort, M., Lagg, A., & Solanki, S. K. 2013, A&A, 557, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Vigeesh, G., Jackiewicz, J., & Steiner, O. 2017, ApJ, 835, 148 [NASA ADS] [CrossRef] [Google Scholar]
  59. Wentzel, D. G. 1992, ApJ, 388, 211 [NASA ADS] [CrossRef] [Google Scholar]
  60. Zuccarello, F., Guglielmino, S. L., & Romano, P. 2014, ApJ, 787, 57 [NASA ADS] [CrossRef] [Google Scholar]
  61. Zwaan, C. 1987, ARA&A, 25, 83 [NASA ADS] [CrossRef] [Google Scholar]

All Figures

thumbnail Fig. 1.

Snapshot of NOOA 11024 sunspot at the early stages of penumbra development (08:40:52 UT) as seen in G-band intensity (upper left), LOS velocity (upper right), total magnetic field strength (lower left), and magnetic field inclination in the LRF (lower right). The pink iso-contour outlines the umbral boundary using a threshold of 50% of the quiet-Sun intensity. The three pink boxes in the G-band intensity maps mark the regions studied in detail. The maps have been rotated according to the solar north-south direction. A movie showing the temporal evolution is available online.

In the text
thumbnail Fig. 2.

Sketch of the projection effects on the derived magnetic properties around a spot at θ = 28°. The influence of the magnetopause (i.e. the most external red lines) on the derived magnetic properties reaches areas farther away from the umbral boundary at the limb side compared to the centre side.

In the text
thumbnail Fig. 3.

Bi-directional flow in elongated granules. Top panels: time slices of selected penumbra-formation region on the sunspot centre side. From top to bottom: LOS velocities, magnetic field strength, inclination (LRF), and G-band intensity. Bottom panels: context maps containing the areas (marked by purple boxes) selected for the time slices as seen in the G-band intensity and LOS velocity at 8:35 UT (two left panels) and at 10:24 UT (two right panels). The violet arrow points the direction along the slice length. The purple and red iso-contours highlight velocities of ±130 m s−1, respectively.

In the text
thumbnail Fig. 4.

Analogous to Fig. 3 but showing the limb-side region where a transient filament with a counter-Evershed flow formed. The bottom panels correspond to observations at 9:23 UT (two left panels) and at 10:01 UT (two right panels).

In the text
thumbnail Fig. 5.

Analogous to Fig. 3 but showing the limb-side region where the granular pattern was observed before the penumbra formation. The bottom panels correspond to observations at 9:00 UT (two left panels) and at 10:24 UT (two right panels).

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.