Issue 
A&A
Volume 569, September 2014



Article Number  A121  
Number of page(s)  6  
Section  The Sun  
DOI  https://doi.org/10.1051/00046361/201424380  
Published online  02 October 2014 
Pair separation of magnetic elements in the quiet Sun
^{1} Dipartimento di FisicaUniversità di Roma Tor Vergata and INFN, via della Ricerca Scientifica, 1, 00133 Rome, Italy
email: Fabio.Giannattasio@roma2.infn.it
^{2} Instituto de Astrofísica de Andalucía (CSIC), Apdo. de Correos 3004, 18080 Granada, Spain
^{3} Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife, Spain
Received: 11 June 2014
Accepted: 14 August 2014
The dynamic properties of the quiet Sun photosphere can be investigated by analyzing the pair dispersion of smallscale magnetic fields (i.e., magnetic elements). By using 25 hlong Hinode magnetograms at high spatial resolution (), we tracked 68 490 magnetic element pairs within a supergranular cell near the disk center. The computed pair separation spectrum, calculated on the whole set of particle pairs independently of their initial separation, points out what is known as a superdiffusive regime with spectral index γ = 1.55 ± 0.05, in agreement with the most recent literature, but extended to unprecedented spatial and temporal scales (from granular to supergranular). Furthermore, for the first time, we investigated here the spectrum of the mean square displacement of pairs of magnetic elements, depending on their initial separation r_{0}. We found that there is a typical initial distance above (below) which the pair separation is faster (slower) than the average. A possible physical interpretation of such a typical spatial scale is also provided.
Key words: Sun: photosphere / Sun: magnetic fields / Sun: granulation
© ESO, 2014
1. Introduction
In the outermost 30% of the solar radius, the transfer of energy towards the surface occurs via turbulent convection. To date, a comprehensive theory of solar turbulent convection from small up to global scales has not yet been formulated. In the last decades, magnetohydrodynamic (MHD) simulations (see, e.g., Nordlund & Stein 1997; Stein & Nordlund 1998, 2001; Beeck et al. 2012) have been extensively used to mimic the uppermost convection zone, as they match very well the observations of the solar photosphere. However, only tiny regions of the Sun can be realistically simulated, because of the wide range of temporal and spatial convective scales and current computer power.
A complementary approach to investigate the properties of convection on the solar quiet photosphere consists in the study of the interaction between convective flows and the smallscale magnetic fields (hereafter magnetic elements) in the interior of supergranular cells. These internetwork magnetic elements can be reasonably regarded as passive objects advected by the underlying flow, as the drag force due to plasma kinetic energy is greater than the magnetic force they exert on the surroundings. Under this assumption (discussed in Sect. 3), the dynamics of magnetic elements describes that of the plasma (see, e.g., Yelles Chaouche et al. 2011; Berrilli et al. 2013). Tracking magnetic elements also allows us to study the onset and amplification of magnetic fields in the quiet Sun, the scales on which they organize, and the rate of interaction between fields (see, e.g., Chitta et al. 2012; van Ballegooijen et al. 2014). This information is important in order to get insights, for example, on the mechanisms that contribute to heating the solar corona, such as magnetic reconnections (see, e.g., Parker 1983; Viticchié et al. 2006) and buffeting induced MHD waves (see, e.g., Stangalini et al. 2013a,b).
Previous studies have tracked Gband magnetic bright points and magnetic elements from magnetograms (we refer to both of them as magnetic features), regarding them as Lagrangian probes. Under this assumption, the mean square displacement of such single magnetic features, namely ⟨Δl^{2}⟩, has been measured and shown to follow a power law ⟨Δl^{2}(τ)⟩ ∝ τ^{γ}, where time τ is defined as starting from the first detection of the magnetic feature. In particular, a spectral index γ = 1 is associated with a normal diffusion (also known as random walk) with constant diffusivity K ∝ ⟨Δl^{2}⟩ /τ. In this case, ⟨Δl^{2}⟩ corresponds to the standard deviation of a Gauss function describing the distribution of displacements.
It is well known that the presence of the combined effects of the velocity field and a superposed diffusion can lead to a very large diffusive coefficient, the socalled eddydiffusivity, which is the only relevant parameter needed to predict the longtime, longspace diffusion scale properties in many applied cases (Moffatt 1983; Bouchard & Georges 1990; Crisanti et al. 1991; Biferale et al. 1995). On the other hand, when there is anomalous diffusion (i.e., γ ≠ 1) the diffusivity depends on both spatial and temporal scales, and superdiffusive (γ> 1) or subdiffusive (γ< 1) regimes can arise. The most recent works in the literature agreed that there was a superdiffusive regime in the quiet Sun (see, e.g., Abramenko et al. 2011; Giannattasio et al. 2013, 2014; Jafarzadeh et al. 2014). This implies that the effective diffusivity decreases when the temporal (spatial) scale is reduced (Abramenko et al. 2011), thus allowing magnetic fields to be enhanced on the very short (small) scales.
Diffusive (normal or anomalous) regimes are typically well defined only in the asymptotic limit of large time and large distance, something that is very difficult to achieve in our case. In many geophysical and astrophysical situations only the transient behavior is observable and/or relevant, hence the interest in discussing in a more quantitative way the separation of pairs of magnetic elements in our observational setup. While the diffusion of single Lagrangian probes is dominated by large scale motions, pair separation should be more universal, being affected only by the relative velocity fields on scales on the order of the distance between the magnetic elements. If we identify with the distance between two magnetic elements in our ensemble at a given time τ, we will be interested in quantifying the probability distribution function (PDF), p(r,τ  r_{0},τ_{0}), of observing a given separation r at time τ starting from an initial distance r_{0} at time τ_{0}. In general, we obtain (1)where δ_{r}v is the velocity difference among the magnetic elements, the average is meant over all the considered pairs in our ensemble (see Sect. 2) and we have assumed that the initial position is uncorrelated from the underlying velocity field (Sokolov 1999). In many cases, the integral in Eq. (1) is well behaved and in the limit of large time converges to ∝⟨v^{2}⟩ τ, i.e., we have asymptotically a normal diffusion process with diffusivity given by the onepoint velocity fluctuations along the trajectories of the magnetic elements. Nevertheless, there are many cases where this asymptotic regime is never reached, or it lies just outside the spatial and temporal limits of observation. For example, it is well known that in the presence of multiscale non trivial statistical properties for the underlying velocity field, anomalous diffusion might develop. This is for instance the case of the celebrated Richardson diffusion (Richardson 1926) for homogeneous and isotropic turbulent flows in the range of spatial scales where the velocity obeys a Kolmogorov 1941 turbulent cascade (Frisch 1995). In this case it is predicted and observed the presence of a superdiffusion with γ = 3 and selfsimilar Richardsonlike PDF (Jullien et al. 1999; Falkovich 2001; Yeung & Borgas 2004; Monin & Yaglom 2007; Salazar & Collins 2009; Bec et al. 2010). Similarly, it is known that in the presence of a spatially smooth velocity field with very long temporal correlation, there might be an anomalous (super or sub) diffusion (see, e.g., the pioneering works of Geisel et al. 1985; Avellaneda & Majda 1992; Zaslavsky et al. 1993). Finally, there might also be the possibility of observing strong anomalous diffusion, i.e., a PDF for pair separation that is not selfsimilar (Castiglione et al. 1999), something that has been recently detected in turbulent flows using a very highstatistical dataset (Scatamacchia et al. 2012). In all cases, we are in the presence of a sort of strong or weak failure of the central limit theorem, i.e., the righthand side of Eq. (1) cannot be trated as the sum of many uncorrelated variables, either because we are exploring spatial and temporal scales too small compared with the characteristic variations of the underlying velocity, or because the velocity field itself possesses nontrivial spatial and temporal asymptotic multiscale properties.
A local (in temporal and spatial scale) effective diffusion coefficient can be defined as (2)which gives us the typical separation speed of two magnetic elements found at separation r after a time τ and with initial separation r_{0}. When the central limit theorem holds, we have K ~ const., independently of the original separation r_{0}. Otherwise, different anomalous regimes arise. In particular, here we are interested in studying the effects of r_{0} on pair separation for the case of magnetic elements in the quiet Sun, where the observational and intrinsic physical limitations do not allow us to extend the observing regimes to couples with separations much larger than r_{0}, and therefore where only the preasymptotic regime is observable. For instance, for fully developed turbulent flows, it is known that the pair separation is only ballistic, namely ⟨Δr^{2}⟩ ~ t^{2}, up to a maximum time, the socalled Batchelor time (i.e., t_{B}, Bourgoin et al. 2006), wich depends on the initial separation and follows the law .
To our knowledge, only Berger et al. (1998) and Lepreti et al. (2012) have applied pair dispersion analysis to magnetic features in the photosphere. Berger et al. (1998) tracked 622Gband bright point pairs acquired at the SVST by observing a region 29′′ × 29′′ wide for 70 min, with a cadence of 25 s and a spatial resolution of ~. They measured a spectral index consistent with γ ≃ 1.3. Lepreti et al. (2012) used NST images (Goode et al. 2010) to track a maximum of 7912 magnetic bright points in a quiet Sun region, a coronal hole, and a plage region. They measured a spectral index γ ≃ 1.5 everywhere in the temporal range 10 ≲ τ ≲ 400 s. They interpreted this difference from known scaling as due to the imperfectly passive nature of bright points.
The limitation in the statistics prevented the aforementioned authors from carrying out pair dispersion analysis taking into account the initial separation r_{0}. In this work, for the first time the pair separation approach is used to investigate the dynamic behavior of magnetic elements in the quiet Sun for different initial pair separations. We believe this is an important point, since it is difficult to imagine that pair separation is not affected by the initial condition for those temporal and spatial scales accessible on the Sun. We present the results obtained with 68 490 tracked magnetic pairs.
2. Observations and data analysis
The data set used in this work was described in Gošić (2012) and analyzed by Giannattasio et al. (2013, 2014) to study the diffusion of single magnetic elements up to supergranular scales. It consists of 995 HinodeNFI magnetograms (Kosugi et al. 2007; Tsuneta et al. 2008) with a spatial resolution of and a noise of σ_{B} = 6 G for single frames. The magnetograms were coaligned, trimmed to the same field of view (FoV, which is ~50 Mm sized), and filtered out for oscillations at 3.3 mHz (Gošić 2012). As a consequence, the present analysis is free from effects of acoustic oscillations and atmospheric seeing, and aimed to magnetic elements and not magnetic proxies like Gband bright points. The large FoV, which encloses an entire supergranule, and the high spatial resolution enabled us to investigate a wide range of spatial scales and observe a large number of magnetic elements. The series, acquired on November 2, 2010, covers 25 h without interruption, with a cadence of 90 s. This allowed us to study the dynamics of magnetic elements on a wide range of temporal scales (from a minute to a day).
In Fig. 1, we show the 25 h timeaveraged magnetogram of the FoV saturated at 25 G. We focus on the region of interest (ROI) inside the green circle, which is centered in the center of the supergranule and has a radius of ~10 Mm, such that the ROI is completely enclosed within the supergranule itself, i.e., in the internetwork region (Giannattasio et al. 2014).
We applied the tracking algorithm described in Del Moro (2004). The algorithm uses a variable threshold in order to overcome the loss of weak fields and resolve the clustered peaks of the largest magnetic features (e.g., Berrilli et al. 2005). We discarded all the magnetic elements with speed v> 7 km s^{1}, which is roughly the speed of sound in the photosphere. We also discarded the magnetic elements passing close to the boundary of the ROI (at a distance ≲1.86 Mm). A total of 68 490 tracked magnetic pairs originated in the ROI has been detected. In Fig. 1 we also show, for the sake of visualization, the evolution of a few magnetic elements forming a subset of pairs.
3. Results and discussion
The detected pairs have been used to investigate the nature of turbulent convection under the hypothesis of magnetic elements passively transported by the supergranular flows in the ROI. By analyzing the same data set, Giannattasio et al. (2013) found that the equipartition magnetic flux density is B_{e} ≃ 255 G. Only less than 4% of the total magnetic elements have an average flux greater than that value, and are all located in the network. Therefore, the condition of passive magnetic elements is reasonably fulfilled in the internetwork regions, as demonstrated by the spectropolarimetric studies performed by Orozco Suárez et al. (2007); Orozco Suárez & Bellot Rubio (2012); and Bellot Rubio & Orozco Suárez (2012).
Fig. 1 Timeaveraged magnetogram saturated at 25 G. Only the magnetic elements inside the ROI (limited by the green circle) are considered for the analysis. A few trajectories of magnetic elements are shown in the ROI. The asterisks mark the first detection positions; the plus signs mark the last detection positions. 
Fig. 2 Mean square pair separation as a function of time since the first appearance. Only the data points up to ~5500 s have been used to make the fit. The error on γ was computed as the standard deviation of the values obtained after a random subsampling of magnetic pairs. The results described in Berger et al. (1998, dotted line) and Lepreti et al. (2012, dashed line), are superposed for comparison. 
As done in the previous works in the literature (see, e.g., Lepreti et al. 2012) the separation spectrum ⟨Δr^{2}⟩ was first computed for all the pairs of magnetic elements, regardless of their initial separation. We obtained the results shown in Fig. 2 (black diamonds). We found that the separation spectrum is best fitted by a power law with spectral index γ = 1.55 ± 0.05 (the red line in the same figure), in agreement with the results of Lepreti et al. (2012), but here extended to supergranular scales. The error on γ was computed as the standard deviation of the values obtained after a random subsampling of the magnetic pairs.
As mentioned in Sect. 1, the high number of magnetic elements tracked allowed us to perform for the first time the pair separation analysis for different values of the initial pair separation. By looking in greater detail at the distribution of initial separations, one easily recognizes that it is in general very broad. In Fig. 3 we show, for each initial separation r_{0}, the number of magnetic pairs found in the range r_{0} − 0.174 < Δr<r_{0} + 0.174 Mm, which is 3 HinodeNFI pixels wide. Vertical bars in the graph represent the errors, which are mainly due to the Poissonian contribution. The peak of the curve lies between 8 Mm and 9 Mm (r_{0,peak}). By comparing Figs. 2 and 3 one can see that the global mean displacement with respect to the initial separation, averaged over all pairs, is of the same order of magnitude as the spread in the initial distribution of r_{0} in our sample, i.e., we have not reached any asymptotic longtime regime. Hence, it is natural to ask the question: How robust is the observed superdiffusive behavior as a function of r_{0}?
As magnetic elements are not pointlike, but have diameters up to d_{min} ≃ 1.86 Mm, it is not possible to chose arbitrarily small mutual distances. Therefore, the minimum pair separation set is d_{min}. Moreover, the maximum achievable separation is given by the diameter of the ROI, which is d_{ROI} ~ 20 Mm. Thus, r_{0} must satisfy d_{min} ≤ r_{0} ≤ d_{ROI}. We choose bins of r_{0} large enough to collect a sufficiently high number of magnetic pairs within, and small enough to be able to study the variation of ⟨Δr^{2}⟩ with r_{0}. For this purpose, we set the bin size at 348 km (i.e., 3 HinodeNFI pixels), and computed ⟨Δr^{2}(τ,r_{0})⟩ for each bin. In Fig. 4 we plot seven of all the computed ⟨Δr^{2}(τ,r_{0})⟩. For comparison, we also overplot the power law behavior (corresponding to γ = 1.55). In order to emphasize the deviations from such a law, in the inset of the same figure we plot a compensated pair separation spectrum . The errors on data points were computed as the standard deviation of the values obtained after a random subsampling of the magnetic pairs. From these two plots we can deduce that 1) the smaller the r_{0}, the smaller the effective eddy diffusivity ⟨Δr^{2}⟩ /τ for any τ; and 2) the smaller the r_{0}, the smaller the value of the effective slope γ; in addition, there is a clear change in the trend for initial separation crossing the value r_{0} ~ 10 Mm, which roughly corresponds to the radius of the ROI. From an observational point of view, 1) and 2) could be interpreted by taking into account the recent results in Orozco Suárez et al. (2012) and Giannattasio et al. (2014). In those works, the authors showed that the horizontal velocity field within a supergranule is mostly radial and directed from the center to the boundaries. Following this sketch, we expect that magnetic elements starting close to each other will, on average, separate more slowly than magnetic elements starting farther away from one another. In fact, magnetic elements with a larger initial separation are most likely to be dragged on along very different directions, thus separating faster. This effect naturally introduces a dependence of ⟨Δr^{2}⟩ on r_{0}. In particular, the systematic increase of the effective γ from the granular to the supergranular scale suggests that the preasymptotic diffusion is a function of the probed spatial scale.
By comparing Figs. 3 and 4 we note that for initial separations smaller than r_{0,peak} the slopes in the pair separation spectrum decrease at longer times; while for larger initial separations the slopes increase. This trend is significant, as can be seen from the errors on the data points shown in the inset in Fig. 4. The pairs of magnetic elements with initial separation around r_{0,peak} (from ~6.5 to ~11.5 Mm), which are in number about half of the entire population of pairs, are characterized by a separation spectrum with γ around ≃1.55. This explains why in the separation spectrum in Fig. 2, which was retrieved by considering all the pairs of magnetic elements in the ROI (with any initial separation), there is an effective trend consistent with γ ≃ 1.55 even at longer temporal scales.
Fig. 3 Number of magnetic pairs as a function of the initial separation. Vertical bars represent errors (see the text). 
Fig. 4 Mean square separation ⟨Δr^{2}(τ,r_{0})⟩ for seven different and equally spaced values of r_{0}. The black solid line corresponds to the fitting curve of Fig. 2. In the inset the compensated mean square separation is shown. The errors (vertical bars) are shown only for a few data points. 
Fig. 5 Timedependent PDF for pair separation of magnetic elements (a), b)), and the same rescaled to unity rms c), d) for initial separations in the range 7.42 <r_{0}< 7.77 Mm a), c) and 10.21 <r_{0}< 10.56 Mm b), d). The shaded areas in panels a) and b) cover the initial separation bins. All PDFs are normalized to unit area. 
To further investigate the effects of r_{0} on the pair separation, we computed the timedependent PDF (normalized at unit area) of observing a given separation starting from the initial values of 7.42 <r_{0}< 7.77 Mm and 10.21 <r_{0}< 10.56 Mm, at which values the change in the trend shown in Fig. 4 is observed. The results are shown in panels (a) and (b) of Fig. 5. In that figure, the initial separation range is depicted as a shaded area. As we can see, the PDF broadens with time, its rms being σ(τ), and the peak moves to gradually increasing separations. At τ = 900 s the tails begin to be important, and affect substantially the pair separation spectrum, indicating that the largest separations begin to become dominant. We rescaled the timedependent PDF so that it is centered at zero, and its rms σ(τ) is unity (see, e.g., Jullien et al. 1999). To this end, we introduced the rescaled separation q = (Δr − ⟨Δr⟩ ) /σ, being ⟨Δr⟩ the mean separation value, and computed PDF(q) at each time. The correct renormalization required to consider σPDF(q) instead of PDF(q). The results are shown in panels (c) and (d) of Fig. 5. As we can see, the curves seem to collapse on each other, especially in proximity of ⟨Δr⟩ (which corresponds to q = 0). This means that the small deviations from the mean values seem to be selfsimilar. When it holds, selfsimilarity indicates that there is a single underlying distribution governing the process at any time (Jullien et al. 1999). However, in our case the statistics is still too low to come to a conclusion in this sense. More data points are nedeed to better sample the tails of the timedependent PDF, where any possible breaking of selfsimilarity can be detected.
4. Conclusions
Dynamic processes in the solar photosphere can be studied at spatial and temporal scales from granular to supergranular by measuring the pair separation rate of quiet Sun magnetic elements. By taking advantage of uninterrupted 25 h magnetograms acquired by Hinode at high resolution and imaging a whole supergranule, we computed the separation spectrum of 68 490 pairs of magnetic elements. When considering all the pairs within the supergranule, we found a spectral index γ = 1.55 ± 0.05, in agreement with the most recent literature, but extended at unprecedented spatial and temporal scales. Such a superdiffusive regime can be interpreted as being due to an underlying velocity field with either characteristic spatial (temporal) scales larger (longer) than the scales of observation; or nontrivial asymptotic multiscale properties. For the first time we investigated the separation spectrum for different values of the initial pair separation of magnetic elements, r_{0}. The main conclusion is that the rate of pair separation depends on the spatial scale under consideration.
The possibility that the preasymptotic diffusive behavior detected here possesses nontrivial multiscaling properties remains to be investigated; in other words, whether higher order moments do not scale proportionally to the second order moment, ⟨(Δr)^{2p}⟩ ≠ ⟨(Δr)^{2}⟩ ^{p}. This would indicate the presence of strong anomalous diffusion (Castiglione et al. 1999), possibly connected to the presence of intermittent properties of the advecting velocity field. This ambitious goal surely represents a great challenge for future research since it can only be achieved by extending by at least one order of magnitude the statistical ensemble and the temporal window of the observation.
Acknowledgments
This work was supported by a Ph.D. grant at the University of Rome, Tor Vergata. Part of this work was done while F.G. was a Visiting Scientist at Instituto de Astrofísica de Andalucía (CSIC). Financial support by the Spanish MEC through project AYA201239636C0605 (including European FEDER funds) is gratefully acknowledged. This work has also benefited from discussions in the Flux Emergence meetings held at ISSI, Bern in December 2011 and June 2012. The data used here were acquired in the framework of Hinode Operation Plan 151, entitled Flux replacement in the solar network and internetwork. L.B. acknowledges partial funding from the European Research Council under the European Community’s Seventh Framework Programme, ERC Grant Agreement No 339032. Hinode is a Japanese mission developed and launched by ISAS/JAXA, collaborating with NAOJ as a domestic partner, N.A.S.A. and STFC (UK) as international partners. Scientific operation of the Hinode mission is conducted by the Hinode science team organized at ISAS/JAXA. This team mainly consists of scientists from institutes in the partner countries. Support for the postlaunch operation is provided by JAXA and NAOJ (Japan), STFC (UK), NASA, ESA, and NSC (Norway).
References
 Abramenko, V. I., Carbone, V., Yurchyshyn, V., et al. 2011, ApJ, 743, 133 [NASA ADS] [CrossRef] [Google Scholar]
 Avellaneda, M., & Majda, A. 1992, J. Stat. Phys., 69, 385 [NASA ADS] [CrossRef] [Google Scholar]
 Bec, J., Biferale, L., Lanotte, A., Scagliarini, A., & Toschi, F. 2010, J. Fluid Mech., 645, 497 [NASA ADS] [CrossRef] [Google Scholar]
 Beeck, B., Collet, R., Steffen, M., et al. 2012, A&A, 539, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bellot Rubio, L. R., & Orozco Suárez, D. 2012, ApJ, 757, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Berger, T. E., Löfdahl, M. G., Shine, R. A., & Title, A. M. 1998, ApJ, 506, 439 [NASA ADS] [CrossRef] [Google Scholar]
 Berrilli, F., del Moro, D., Florio, A., & Santillo, L. 2005, Sol. Phys., 228, 81 [NASA ADS] [CrossRef] [Google Scholar]
 Berrilli, F., Scardigli, S., & Giordano, S. 2013, Sol. Phys., 282, 379 [NASA ADS] [CrossRef] [Google Scholar]
 Biferale, L., Crisanti, A., Vergassola, M., & Vulpiani, A. 1995, Phys. Fluids 7, 2725 [Google Scholar]
 Bouchaud, J. P., & Georges, A. 1990, Phys. Rep., 195, 127 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Bourgoin, M., Ouelette, N. T., Xu, H., Berg, J., & Bodenschatz, E. 2006, Science, 311 , 835 [Google Scholar]
 Castiglione, P., Mazzino, A., MurtoreGinanneschi, P., & Vulpiani, A. 1999, Physica D, 134, 75 [NASA ADS] [CrossRef] [Google Scholar]
 Chitta, L. P., van Ballegooijen, A. A, Rouppe van der Voort, L., DeLuca, E. E., & Kariyappa, R. 2012, ApJ, 752, 48 [NASA ADS] [CrossRef] [Google Scholar]
 Crisanti, A., Falcioni, M., Paladin, G., & Vulpiani, A. 1991, Nuovo Cimento 14, 1 [Google Scholar]
 Del Moro, D. 2004, A&A, 428, 1007 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Falkovich, G., Gawedzki, K., & Vergassola, M. 2001, Rev. Mod. Phys., 73, 913 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Frisch, U. 1995, Turbulence. The legacy of A. N. Kolmogorov (Cambridge University Press) [Google Scholar]
 Geisel, T., Nierwetberg, J., & Zachert, A. 1985, Phys. Rev. Lett., 54, 616 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Giannattasio, F., Del Moro, D., Berrilli, F., et al. 2013, ApJ, 770, L36 [NASA ADS] [CrossRef] [Google Scholar]
 Giannattasio, F., Stangalini, M., Berrilli, F., Del Moro, D., & Bellot Rubio, L. 2014, ApJ, 788, 137 [NASA ADS] [CrossRef] [Google Scholar]
 Goode, P. R., Coulter, R., Gorceix, N., Yurchyshyn, V., & Cao, W. 2010, Astron. Nachr., 331, 620 [NASA ADS] [CrossRef] [Google Scholar]
 Gošić, M. 2012, Master’s thesis, Univ. Granada, Spain [Google Scholar]
 Jafarzadeh, S., Cameron, R. H., Solanki, S. K., et al. 2014, A&A, 563, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Jullien, M.C., Paret, J., & Tabeling, P. 1999, Phys. Rev. Lett., 82 , 2872 [NASA ADS] [CrossRef] [Google Scholar]
 Kosugi, T., Matsuzaki, K., Sakao, T., et al. 2007, Sol. Phys., 243, 3 [NASA ADS] [CrossRef] [Google Scholar]
 Lepreti, F., Carbone, V., Abramenko, V. I., et al. 2012, ApJ, 759, L17 [NASA ADS] [CrossRef] [Google Scholar]
 Moffatt, H. K. 1983 Geophys. Res., Rep. Prog. Phys., 46, 621 [Google Scholar]
 Monin, A. S., & Yaglom, A. M. 2007, Statistical Fluid Mechanics, Vol. 2: Mechanics of Turbulence, MIT [Google Scholar]
 Nordlund, A., & Stein, R. 1997, in Proc. SCORe’96: Solar Convection and Oscillations and their Relationship, 225, 79 [Google Scholar]
 Orozco Suárez, D., & Bellot Rubio, L. R. 2012, ApJ, 751, 2 [NASA ADS] [CrossRef] [Google Scholar]
 Orozco Suárez, D., Bellot Rubio, L. R., del Toro Iniesta, J. C., et al. 2007, ApJ, 670, L61 [NASA ADS] [CrossRef] [Google Scholar]
 Orozco Suárez, D., Katsukawa, Y., & Bellot Rubio, L. R. 2012, ApJ, 758, L38 [NASA ADS] [CrossRef] [Google Scholar]
 Parker, E. N. 1983, ApJ, 264, 642 [NASA ADS] [CrossRef] [Google Scholar]
 Richardson, L. F. 1926, Proc. Roy. Soc. London, Ser., A, 110, 709 [Google Scholar]
 Salazar, J. P. L. C., & Collins, L. R. 2009, Ann. Rev. Fluid Mech, 41, 405 [NASA ADS] [CrossRef] [Google Scholar]
 Scatamacchia, R., Biferale, L., & Toschi, F. 2012, Phys. Rev. Lett. 109, 144501 [NASA ADS] [CrossRef] [Google Scholar]
 Sokolov, I. 1999, Phys. Rev. E, 60, 5528 [NASA ADS] [CrossRef] [Google Scholar]
 Stangalini, M., Solanki, S. K., Cameron, R., & Martínez Pillet, V. 2013a, A&A, 554, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Stangalini, M., Berrilli, F., & Consolini, G. 2013b, A&A, 559, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Stein, R. F., & Nordlund, A. 1998, ApJ, 499, 914 [Google Scholar]
 Stein, R. F., & Nordlund, Å. 2001, ApJ, 546, 585 [NASA ADS] [CrossRef] [Google Scholar]
 Tsuneta, S., Ichimoto, K., Katsukawa, Y., et al. 2008, Sol. Phys., 249, 167 [NASA ADS] [CrossRef] [Google Scholar]
 van Ballegooijen, A. A., AsgariTarghi, M., & Berger, M. A. 2014, ApJ, 787, 87 [NASA ADS] [CrossRef] [Google Scholar]
 Viticchié, B., Del Moro, D., & Berrilli, F. 2006, ApJ, 652, 1734 [NASA ADS] [CrossRef] [Google Scholar]
 Yelles Chaouche, L., MorenoInsertis, F., Martínez Pillet, V., et al. 2011, ApJ, 727, L30 [NASA ADS] [CrossRef] [Google Scholar]
 Yeung, P. K., & Borgas, M. S. 2004, J. Fluid Mech., 503, 125 [NASA ADS] [CrossRef] [Google Scholar]
 Zaslavsky, G. M., Stevens, D., Weitzener, H. 1993, Phys. Rev. E, 48, 1683 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
All Figures
Fig. 1 Timeaveraged magnetogram saturated at 25 G. Only the magnetic elements inside the ROI (limited by the green circle) are considered for the analysis. A few trajectories of magnetic elements are shown in the ROI. The asterisks mark the first detection positions; the plus signs mark the last detection positions. 

In the text 
Fig. 2 Mean square pair separation as a function of time since the first appearance. Only the data points up to ~5500 s have been used to make the fit. The error on γ was computed as the standard deviation of the values obtained after a random subsampling of magnetic pairs. The results described in Berger et al. (1998, dotted line) and Lepreti et al. (2012, dashed line), are superposed for comparison. 

In the text 
Fig. 3 Number of magnetic pairs as a function of the initial separation. Vertical bars represent errors (see the text). 

In the text 
Fig. 4 Mean square separation ⟨Δr^{2}(τ,r_{0})⟩ for seven different and equally spaced values of r_{0}. The black solid line corresponds to the fitting curve of Fig. 2. In the inset the compensated mean square separation is shown. The errors (vertical bars) are shown only for a few data points. 

In the text 
Fig. 5 Timedependent PDF for pair separation of magnetic elements (a), b)), and the same rescaled to unity rms c), d) for initial separations in the range 7.42 <r_{0}< 7.77 Mm a), c) and 10.21 <r_{0}< 10.56 Mm b), d). The shaded areas in panels a) and b) cover the initial separation bins. All PDFs are normalized to unit area. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext 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 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.