A&A 419, 757-762 (2004)
DOI: 10.1051/0004-6361:20035739

Relation between families of granules, mesogranules and photospheric network

Th. Roudier - R. Muller

Laboratoire d'Astrophysique de Toulouse, Observatoire Midi-Pyrénées, 57 avenue d'Azereix 65008 Tarbes, France

Received 25 November 2003 / Accepted 13 February 2004

Abstract
The analysis of a solar granulation 3-h time sequence obtained at the Pic du Midi Observatory confirms the existence of Trees of Fragmenting Granules (TFG). The TFG lifetime histogram is fitted by a power law with an exponent equal to -1.68. The positive divergences (i.e. mesogranules) come from several families that are visible at different stages of the mesogranule evolution. We observe a good coincidence between the cork (passive scalars) distribution and the TFG boundaries. Thus, TFGs seem to play a role in the diffusion of the magnetic elements on the Sun surface.

Key words: Sun: photosphere - Sun: granulation

1 Introduction

The magnetic field in the quiet sun is composed in part of discrete flux tubes which form the well-known photospheric network. These flux tubes are buffeted and deformed continuously by granulation and diffuse towards the boundaries of the supergranule. The mechanism which drives the diffusion of the magnetic field towards the supergranule borders is not well understood. In the classical view there are three convective scales, granulation, mesogranulation and supergranulation, which are produced by the de-ionization of hydrogen or helium (Unsöld 1930; November et al. 1981; Simon & Weiss 1968). Simulations have shown that the combination of the motions generated by these three convective scales is a possible mechanism of diffusion of the magnetic elements towards the supergranule limits (Simon et al. 1991).

In recent work a different explanation has been presented on the origin of the meso and supergranulation. Rieutord et al. (2000) have proposed that, in all likelihood, mesogranulation is not a true scale of solar convection, but the combination of the effects of both highly energetic granules and data processing. According to these authors the supergranulation is the result of a nonlinear large-scale instability of the granular flow, triggered by exploding granules. Similarly, Rast (2003) claims that the scale of supergranulation has its origin in the interaction and merging of individual granular plumes. To understand the faster rotation of the supergranulation relative to the magnetic features, Gizon & Duvall (2003) suggest a wavelike nature of the supergranulation. It is therefore interesting to revisit the interaction of the various dynamical scales on the sun surface which can advect the flux tubes.

In this paper, we will limit our study to mesoscale evolution and interactions. Recently, Roudier et al. (2003) have shown that a significant fraction of the granules in the photosphere is organized in families of granules which are named "Trees of Fragmenting Granules'' (TFGs). A TFG consists of a family of repeatedly splitting granules, originating from a single granule at its beginning. The TFGs can live much longer (up to 8 h) than individual granules (10 min). When averaged in time, such long-lived TFGs produce coherent diverging flows which are identified as mesogranules. In the following, long-lived TFGs are defined as TFGs with lifetime greater than 1h30.

The aim of the present paper is to determine the relations between the divergence of the flow vector associated with mesogranulation, the TFGs and Network Bright Points (NBPs). For this purpose, we use the data that have been analyzed by Muller et al. (1992), where mesogranule identified as the positive divergences move towards supergranule boundaries.

2 Data

The granulation movie used for our investigation was obtained with the 50 cm refractor of the Pic du Midi Observatory, on September 20, 1988. It consists of 526 white light ( $5750~{\rm\AA}$) frames taken every 20 s. Each was selected in a burst of 50 frames. Most selected frames are of very good quality and near the diffraction limit (0 $.\!\!^{\prime\prime}$25). The field of view is $100\hbox{$^{\prime\prime}$ }$ $ \times$ $70\hbox{$^{\prime\prime}$ }$ and the pixel size is 0 $.\!\!^{\prime\prime}$14. The photographs were digitized with a $1024 \times 1024$ CCD camera at Lockheed Palo Alto Laboratory and then the frames were processed: aligned, destreched and finally filtered to remove 5-min p-mode oscillations. Our analysis is restricted to an area of $58\hbox{$^{\prime\prime}$ }$ $ \times$ $48\hbox{$^{\prime\prime}$ }$ common to all the frames, which includes a large part of a supergranule. The boundary of the supergranule is clearly outlined by the converging flows of granules and by the location of Network Bright Points (Muller et al. 1992).

3 TFGs detection and lifetime

The method of detecting families of splitting granules (TFGs) has been described in Roudier et al. (2003). To summarize, the intensity map are transformed into a binary maps using the local curvature of the intensity field (Strous 1994). When the second spatial derivatives have the same sign in the four directions, for example negative, the pixel belongs to a bright object (granule), in the opposite case it belongs to a dark intergranule. Once the segmentation is performed the objects are labelled. The evolution of the granules is tracked by identifying the temporal links in the 3D data box (2D in space and 1D in time). All the granules originating from a single progenitor, by splitting in sucessive generations, are labelled with the same label as the progenitor. The end of a family is determined when the last granules of the family dilute in most cases or merge in a granule with a larger area. This is the criterium to label in time the TFG in the case of merging.


  \begin{figure}
\par\includegraphics[width=7.7cm,clip]{0739fig1.ps}\end{figure} Figure 1: Lifetime histogram of individual granules and TFGs. It is well fitted by a power law with an exponent equal to -1.68.
Open with DEXTER

We constructed in this way 10 744 objects (the term object stands for a set of granules with the identical label) contained within the 3D (x,y,t) box. Among those, only 7471 had a lifetime greater than one frame. Of those 7471 objects, 512 had at least one splitting and formed the set of what we call "Trees of Fragmenting Granules" (TFGs). The remaining 6959 objects are called individual granules.

The lifetime histogram of individual granules and TFGs (Fig. 1) is well fitted by a power law $ N\propto t^\gamma$, with $\gamma =-1.68$ which is very close to the exponent found in Roudier et al. (2003) $\gamma =-1.72$. Such fit indicates that we cannot give a characteristic lifetime to the TFGs.

4 Relation between mesogranules and TFGs

It has been shown in Roudier et al. (2003) that horizontal flow fields of long-lived TFG computed by using only the granules belonging to that family, produce coherent diverging flows. These coherent diverging flows have been identified as mesogranules.

In this paper we extend the study of the relation between mesogranules and TFGs which was restricted up today to long-lived TFGs. The mesogranules detected as divergences by Muller et al. (1992) are used for that purpose. To detect the same features, we have computed here the divergence of the flow field with the same temporal and spatial windows, 17 min and 3 $.\!\!^{\prime\prime}$3 FWHM, respectively. They are shown in Fig. 2 at the same instant as in Fig. 2 (first top left) of Muller et al. (1992), and with the same labels. The windows parameters have been chosen to remove the flow pattern produced by individual granules. Thus, for each mesogranule identified in Fig. 2, we determined the granules and family numbers.


  \begin{figure}
\par\includegraphics[width=8.1cm,clip]{0739fig2.ps}\end{figure} Figure 2: Flow divergence (t=9 min) identical to Fig. 2 first image left column of Muller et al. (1992). Mesogranule numbers are indicated over the positive divergence selected by Muller et al. (1992).
Open with DEXTER

Figure 3 shows different families with the flow field at the same location. Some divergences visible in Fig. 2 are indicated by circles in Fig. 3. At these locations, we observe positive diverging flow associated with TFGs and confirm the result of Roudier et al. (2003).


  \begin{figure}
\par\includegraphics[width=14cm,clip]{0739fig3.eps}\end{figure} Figure 3: Families ( TFGs > 7200 s) together with the flow field at the same location. Circles indicate the location of some divergence visible in Fig. 2.
Open with DEXTER

In the following, we describe the most representative examples of the links between divergences and TFGs.

In the first divergence map of Fig. 4 (t=9 min; top left), the divergence labelled 5 corresponds to four families This is not surprising because the granule labeling of the first frame of our observation does not allow us to take care of the full granule history. So the detected divergence 5 at the beginning of the observation cannot rule out the link between the TFGs and the divergence. But during the time evolution of the divergence, new families are growing inside the mesogranule identified by the number 5. As shown in Fig. 4 at t=92 min, a new family is detected. This new family then pushes the old ones. This event repeats at t=171 min. Particular care has been taken to verify that the families are new ones. So Muller's mesogranule 5 which is identified as the same structure during the time sequence, corresponds in fact to the evolution of different families, unresolved by the spatial (3 $.\!\!^{\prime\prime}$3) and temporal (17 min) windows.


  \begin{figure}
\par\includegraphics[width=7.7cm,clip]{0739fi4a.eps}\hspace*{6mm}...
...m}\includegraphics[width=7.7cm,clip]{0739fi4e.eps}\hspace*{4.7cm}}
\end{figure} Figure 4: Time evolution of the divergence 5 of Muller et al. (1992), and the corresponding TFGs located in the area of this divergence. TFG numbers are indicated on each map.
Open with DEXTER

Table 1: Number of families (lifetime > 7200 s) found in each mesogranule during its life (i.e. Fig. 2).

Thus, the motion described in Muller et al. (1992) may appear as the combination of the evolution of different families and not as the motion of the same entity.

Another example is divergence number 15. It is quite stable during the time sequence; it is produced most of the time by only one family. An extension of Muller's mesogranule 15 is due to a close second family (t=117 min). This is also the case of the divergence number 2 which is composed of a family visible during all the time sequence and four adjacent families which appear at different times.

More generally, several long-lived families evolve in the divergences taken into consideration in Muller et al. (1992), which are among the largest and strongest detected in the field of view. Table 1 summarize the number of long-lived families for each divergence. In addition to the long-lived families, many shorter families appear and disappear, contributing to the divergence. A family of granules has its own identity. The spatial extension and the lifetime of the pattern depends of the characteritics (time averaging and spatial smoothing) of the windows to detect them. When 17 min and 3 $.\!\!^{\prime\prime}$3 windows are used, several families are associated with the largest divergences. This confirms previous results of Rieutord et al. (2000) who show that the increase of the size of windows tends to mix structures which are here the TFGs. The use of, divergence quantity is not suited to keep the temporal link of the positive divergence. On the contrary the temporal evolution of structure is followed by using the TFGs.

5 TFG proper motion


  \begin{figure}
\par\includegraphics[width=8.1cm,clip]{0739fi5a.eps}\hspace*{1cm}...
...c.eps}\hspace*{1cm}
\includegraphics[width=8.1cm,clip]{0739fi5d.eps}\end{figure} Figure 5: Time evolution of the TFGs overlapped with the displacement of their centre of gravity. Each TFG displacement represents the cumulative displacement between the beginning of the time sequence and the time of each plot. Encircled TFGs are those which change direction during their lifetime.
Open with DEXTER

To describe the evolution and the interactions of the long-lived TFGs (lifetime > 7200 s), we plotted the cumulative view of the TFGs at different times. The cumulative view has the advantage of showing directly the dominant structures, here TFGs, at each step of the time evolution in the field of view. The evolution is shown in Fig. 5 where the cumulative displacements of the centre of gravity of the families have been overplotted. Then depending of their locations, we observe a larger TFG displacement at mean distance between the centre and the limit of the supergranule outlined by the NBPs in Fig. 6. In the center of each panel of Fig. 5 we note that the motions of the TFGs are limited by their environment. Regarding the direction of their displacement, we observe that some of them change direction during their life (see encircled TFG in Fig. 5). However, displacements at the beginning of the sequence are more noisy due to the shorter time interval and they are probably influenced by older TFGs born before the beginning of the sequence. The movie of the long-lived TFGs shows that all families adjust their development relative to their neighbouring family, the most energetic ones dictating the evolution of others. This is similar to the description of the granule evolution given by Schriver et al. (1997).

At the end of our time sequence, we observe that long-lived families (lifetime > 7200 s, 2 h) cover pratically all the area and TFG displacements exhibit the large scale pattern which can be associated with the supergranular flows. The associated velocities of the centre of gravity of the TFGs lies between 0 and 0.8 km s-1with a peak around 0.3 km s-1. The divergence motions measured by Muller et al. (1992) are similar to those shown in Fig. 5. The similarity of their motions confirms the close relation between families and divergences described in the previous section and in Roudier et al. (2003). Then both long-lived TFGs and centres of divergence appear to be advected by supergranular flow.


  \begin{figure}
\par\includegraphics[width=8.1cm,clip]{0739fig6.ps}\end{figure} Figure 6: NBP location (cross) and the density of TFGs > 5400 s. The density of TFGs is proportional to the brightness in the image.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.1cm,clip]{0739fig7.eps}\end{figure} Figure 7: Cumulative view of the TFGs at the end of the time sequence and the final corks distribution after three hours of advection.
Open with DEXTER

Families of shorter lifetimes behave like the longer ones; they adjust their development relative to their environment and they move coherently with the long-lived family, toward the boundary of the supergranule.

6 Relation between Network Bright Points and TFGs

In Fig. 6, we plotted all the Network Bright Points (NBPs) identified in our 3-h time series and the sum of all the granules belonging to the families of granules living more than 5400 s. The NBPs outline roughly the boundary of a supergranule (Muller et al. 1992). In Fig. 6 the location of the NPBs coincides with the reduced density of long lived-families. Moreover, at the supergranule boundary are either short-lived or simply missing.


  \begin{figure}
\par\includegraphics[width=8.1cm,clip]{0739fig8.eps}\end{figure} Figure 8: Location of the NBPs (cross) and cork at the end of the time sequence.
Open with DEXTER

One way to know more about the role of TFGs in the diffusion of the magnetic elements is to compute the advection of floating corks by the granular flow. This technique was already used by Simon et al. (1991). It consists of following the trajectories of floating corks, initially uniformly distributed, and characterizing their spatial distribution after some time. Here we use the techniques described in Roudier et al. (1999). We determine the velocity field evolution of the Pic du Midi data set with a high spatial and temporal resolution (0 $.\!\!^{\prime\prime}$7 and 5 min). We then use this field to compute the cork trajectories during the three hours of recorded data. Figure 7 shows the cumulative view of the TFGs at the end of the time sequence and the final cork distribution after three hours of advection. As shown in this figure, corks are expelled from TFGs and concentrate in regions which coincide with limits of the TFGs. This result is quite remarkable because the computation of the velocity field which drives the corks motions and the TFG detection is independent. The coincidence of the corks distribution and the TFGs indicates that TFGs appear to play a role in the diffusion of the magnetic elements on the Sun surface. A previous result obtained with the same data, described in Rieutord et al. (2000), shows that corks concentrate in regions that coincide with those occupied by network bright points which delineated the supergranulation scale (Fig. 8). This suggests a key role of the TFGs in the building of the photospheric magnetic network in the quiet Sun.

7 Conclusion

From the analysis of the solar granulation 3-h time sequence obtained at the Pic du Midi Observatory, we confirm the existence of TFGs. We find a remarkable fitting of the TFG lifetime histogram of individual granules by a power law whose exponent is identical to the one found by Roudier et al. (2003). We confirm also the link between TFGs and positive divergences. The detailed analysis of this relation shows that divergences come from several families that are visible at different stages of the evolution. The spatial (3 $.\!\!^{\prime\prime}$3) and temporal (17 min) windows used in Muller et al. (1992) to identify divergences (i.e. mesogranules) were selected to remove the granulation scale. In this way they observed scales larger than granulation but with the tendency to mix structures, namely the TFGs.

The TFGs adjust their development relative to their family neighbours, the more energetic ones dictating the evolution of the others. The long-lived TFGs investigated here move like the divergences analysed by Muller et al. (1992): they move outward from the centre to the boundary of the supergranule, at supergranular speed, when they are located at mid-distances; they are more stationary when they evolve close to the centre or to the boundary. Thus both long-lived TFGs and divergence appear to be advected by supergranular flow.

The NBPs are located where the density of long-lived family of granules is in general low at the limit of the supergranule. The corks which are used as passive scalars are expelled from TFGs and concentrate in regions which coincide with the limits of the TFGs. This result is quite remarkable because the computation of the velocity field that drives the cork motion and the TFG detection are independent. This result must be compared with the persistent magnetic pattern found by Dominguez Cerdeña et al. (2003), which is visible at the spatial scale of the mesogranulation. The coincidence of the cork distribution and the TFG limits indicates that TFGs seem to play a role in the diffusion of the magnetic elements on the Sun surface.

While TFG's displacements exhibit large scale pattern that can be associated with the supergranule flows, some questions arise from our analysis, for instance:

- what is the degree of contribution of the TFGs to the diffusion of magnetic field?

- how do TFGs contribute to the building of the photospheric network?

The answer to such questions only will be possible by developing simulations in which the different sources, at different scales, will be easily identified.

Acknowledgements
We would like to thank very much Jose Antonio Bonet for very fruitful and helpful comments on the manuscript.We thank Michel Rieutord for useful advice and comments. This work was supported by the Centre National de la Recherche Scientifique (CNRS, UMR 5572) and the Programme National Soleil Terre (PNST). Special thanks are due to the Pic du Midi Observatory staff for their technical assistance.

References



Copyright ESO 2004