Issue 
J. Space Weather Space Clim.
Volume 6, 2016
Statistical Challenges in Solar Information Processing



Article Number  A27  
Number of page(s)  8  
DOI  https://doi.org/10.1051/swsc/2016022  
Published online  27 June 2016 
Research Article
Tracking of magnetic flux concentrations over a fiveday observation, and an insight into surface magnetic flux transport
Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency, 311 Yoshinodai, Chuo, Sagamihara, Kanagawa
2525210, Japan
^{*} Corresponding author: iida@solar.isas.jaxa.jp
Received:
31
January
2015
Accepted:
3
June
2016
The solar dynamo problem is the question of how the cyclic variation in the solar magnetic field is maintained. One of the important processes is the transport of magnetic flux by surface convection. To reveal this process, the dependence of the squared displacement of magnetic flux concentrations on the elapsed time is investigated in this paper via a featurerecognition technique and a continual fiveday magnetogram. This represents the longest time scale over which a satellite observation has ever been performed for this problem. The dependence is found to follow a power law and differ significantly from that of diffusion transport. Furthermore, there is a change in the behavior at a spatial scale of 10^{3.8} km. A superdiffusion behavior with an index of 1.4 is found at smaller scales, while changing to a subdiffusion behavior with an index of 0.6 on larger ones. We interpret this difference in the transport regime as coming from the networkflow pattern.
Key words: Sun / Magnetic field / Photosphere / Surface / Solar cycle
© Y. Iida, Published by EDP Sciences 2016
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
Magnetic fields on the solar surface are a fundamental driver of solar activity, with its effects on the SunEarth system. The solar dynamo problem, whose main manifestation is the 11year cycle in solar activity, and the 22year cycle in the solar magnetic field, is of great concern to solar physics and geoscience.
The fluxtransport model is one of the plausible scenarios to account for this cyclic variation. The basic theory behind the model was proposed by Babcock and Leighton (Babcock 1961; Leighton 1964). One important process in this model is poleward transport of magnetic flux at the solar surface. Sunspots appear near 30° latitude with a certain lean angle from the longitudinal direction and are dispersed within one or two months; in lineofsight observations of the photospheric magnetic field, the Sun reveals a patchy structure at small scales (≤10^{3} km), called magnetic flux concentration. The magnetic field can strengthen until it becomes an equipartition field, at which stage the magnetic energy density is equal to that of the kinetic energy; it may strengthen even further by the process of convective collapse. However, most of the magnetic field remains weaker than the equipartition field, and the motion of the concentration is thought to be dominated by that of the plasma. Hence plasma flows play an important role in magnetic field transport. These flows can be divided into two categories, globalscale and smallscale flows.
There are two kinds of globalscale flows at the solar surface. One is the differential rotation, which is the variation in the rotational velocity as a function of latitude and expanding the surface magnetic region along the longitudinal direction. The other one is the meridional flow, which is a flow pattern from the equator to the pole in each hemisphere. It transports the magnetic concentration to the pole and is important for the length of the solar cycle. The amplitude of the meridional flow is of the order of 10 m s^{−1} (Wang et al. 1989; Beck et al. 2002); it takes ∼1 year for the magnetic flux concentrations to travel the distance of the solar radius (∼7.0 × 10^{5} km). However, if only these globalscale flows were considered, then all the flux of positive and negative polarities in the sunspots would reach the pole and there would be no polarity reversal, e.g., no cyclic variation in the magnetic field.
The surface convective flow plays an important role here. Convective flow patterns at the solar surface are multiscale (Nordlund et al. 2009). There are three kinds of convections on the solar surface; the smallest but strongest convective flow pattern is called granulation. This can be observed directly in the intensity images of photospheric continuum. The spatial scale of these granules is ∼10^{3} km, and their amplitude is ∼1 km s^{−1}), which is larger than that of the meridional flow by two orders of magnitude. Mesogranulation occurs on a larger scale than granulation November et al. (1981) and has a spatial scale of ∼7 × 10^{3} km and a horizontal velocity of 500 m s^{−1}. However, it is difficult to see the significant signature corresponding to mesogranulation. Yelles Chaouche et al. (2011) investigated relationships between horizontal flow pattern from intensity images and magnetic field by the use of high spatial and temporal datasets obtained by the Imaging Magnetograph eXperiment (IMaX). Although there was a high correlation between the horizontal flow field and motion of magnetic concentrations, they found no significant scales within the range of <10 Mm in probability density functions of footpoint separation distance. Berrilli et al. (2013) investigated length scale of void structure of lineofsight magnetic field but they found no preferred scales of organization in the range of 2–10 Mm and suggested that the multiscale nature of flows on the solar surface masks a mesogranular scale. The last and largest one type of solar convective flow is supergranulation. It has a spatial size of ~1.6 × 10^{4} km, and a typical horizontal speed of 300–500 m s^{−1} (see Rieutord & Rincon 2010 and references therein).
Because these surface convections are much faster than the meridional flow, the magnetic flux concentration rapidly prevails during transport to the pole. Considering the inclination of sunspots and antipolarity between the hemispheres, this prevails across the hemisphere, resulting in imbalances of the cancellation occurrence between the polarities in each hemisphere. Hence, the flux to the pole is nonzero and changes in the magnetic polarity take place.
Wang et al. (1989) investigated this model by using numerical simulation. They calculated the time evolution of the surface magnetic field by changing the meridional flow speed and the diffusion coefficient for smallerscale convections. The simulation with 10–20 m s^{−1} as meridional flow and 600 km^{2} s^{−1} as convective diffusion was found to produce the results that are the most consistent with the observations.
However, recent observations reveal that the magnetic field transport has a character that is different from the normal diffusion regime. The relationship between the elapsed time and the squared displacement of the magnetic flux concentration is therefore investigated there because their proportionality is a necessary condition for the diffusion transport. The critical difference in the magnetic field observation between Leighton’s era and the present is that the tracking of each magnetic flux concentration is possible now, allowing us to directly investigate the magnetic flux transport. The tracking of the patches requires a stable spatial and temporal resolution better than 10^{3} km and 1 min, respectively, which could not be accomplished in Leighton’s era. Abramenko et al. (2011) summarize the recent situation and the superdiffusion scaling with an index varying from 1.27 to 1.67 is reported in recent work based on ground observation and satellite observations (Hagenaar et al. 1999; Lawrence et al. 2001; Giannattasio et al. 2014b; Caroli et al. 2015).
However, there is still a large gap between the spatial scale of tracking and the global scale. This is because of limitations in the observation period. Previous papers have mainly been devoted to time scales below one day. We need longer observations to investigate transport over larger scales.
Here, we investigate the relationship between the elapsed time and square of the displacement of the magnetic concentration in the longer time scale. Recent developments in technology have facilitated this kind of investigation in two ways; first, featurerecognition and tracking techniques of the magnetic flux concentration have advanced over the last few decades (Hagenaar et al. 1999; DeForest et al. 2007; Parnell et al. 2009; Thornton & Parnell 2011). The analysis of huge amounts of events is needed for investigating solar surface transport. This is because the magnetic elements are much smaller than the global scale, and hence, the statistical character is crucial; featurerecognition techniques are plausible solutions to this difficulty. Secondly, continual and uniform magnetogram data are now being provided by the satellites. The uniformity of the dataset is important for the featurerecognition method. Moreover, while the maximum groundobservational period is limited by the Earth’s rotation, satellite observation is free from this limitation and is thus preferable for this kind of study.
In Section 2, we briefly summarize the basic concept connecting the global transport regime and the motion of each element in 1D randomwalk modeling. Descriptions of the observational data and tracking algorithm are shown in Sections 3 and 4, respectively. The result is shown in Section 5. Section 6 is devoted to a comparison between the present data and previous observations and theoretical requirements.
2. Global transport regime and motion of each element
In this work we find that the dataset of magnetic field concentration exhibits properties of nonFickian diffusion. We will focus on analyzing the relation between mean squared displacements and time and will not explore concept related to the fractal dimensionality of the space. Therefore to introduce the main concepts we present here the case of the 1D randomwalk model leading to normal diffusion.
A 1D random walk is defined as a motion that has constant travel distance during one walk and constant waiting time between walks. Figure 1 shows a schematic example of a 1D random walk. The expected position of the particle is clearly 〈x(t)〉 = 0 at any time under random motion. However, the expected dispersion of the particle position, e.g., the squared displacement from t = 0, increases with time and is calculated as$$\langle {x}^{2}\left(t\right)\rangle =\left(\frac{{\delta}^{2}}{\tau}\right)t,$$(1)where τ and δ are the temporal and spatial scales for one jump. This result shows that the squared displacement is proportional to the elapsed time in the random walk of a particle. The analytical solution of the diffusion transport has the same character. The diffusion transport of the particle density in 1D is written as$$\frac{\mathrm{\partial}n\left(x,t\right)}{\mathrm{\partial}t}=D{\nabla}^{2}n\left(x,t\right),$$(2)where n(x,t) is the particle density at position of x and time of t and D is a diffusion coefficient. The expectation of 〈x^{2}(t)〉 can be calculated as$$\langle {x}^{2}\left(t\right)\rangle ={\int}_{\mathrm{\infty}}^{\mathrm{\infty}}\mathrm{}{x}^{2}n(x,t)\mathrm{d}x=2\mathrm{Dt}.$$(3)
Fig. 1. Schematic picture of a 1D randomwalk motion. τ and δ are, respectively, temporal and spatial scales for one walk. 
Again, we see that the expectation value of the squared displacement is proportional to the elapsed time. Diffusion transport gives the same expectation value of the dispersion as the randomwalk particle with a condition D = δ^{2}/2τ. Note that this proportionality is a necessary condition for the description by diffusion transport on a global scale.
Another simple case is ballistic transport. In this case, the particle motion is written as dx(t)/dt = v_{0} and hence . Here, the squared displacement is proportional to the square of the elapsed time.
This concept, the relationship between the elapsed time and the squared displacement, can be expanded to the more general transport regime, which is called nonFickian diffusion (Balescu 1988), e.g. anomalous diffusion. Figure 2 shows a schematic picture of this concept. The subdiffusion regime is where the powerlaw index is less than 1, meaning that the spreading of the particles slows down as compared to what would happen in a diffusive regime. A typical mechanism of this regime is trapping with random walks. When there are trapping points, some of the particles are captured and the average spreading velocity becomes slower than in the diffusion regime. In contrast, there is a superdiffusive regime between diffusion and ballistic motion. The powerlaw index then has a value between 1 and 2, indicating faster spread than under ballistic motion. A typical case of superdiffusion is the Lévy flight, wherein the velocity of the particle has a powerlaw distribution rather than a Gaussian one. Above ballistic motion, we find the hyperdiffusion regime. It has a powerlaw index larger than 2 and the pervading speed increases as it spreads. A selfavoiding random walk in which the paths of the elements never interact with each other is a typical case of the hyperdiffusion. However, this is not the case for solar surface transport. There is another idea to explain the nonFickian diffusion mechanism. Schrijver & Martin (1990) investigated motions of magnetic concentrations in the core and surroundings of active regions, and found different amplitudes for the diffusion coefficient. They discussed that the scatter in the amplitudes may be caused by differences in step length because there were no significant differences in velocity distributions. If the diffusion coefficient varies with space, then the diffusive transport may be globally anomalous, while being diffusive within each region.
Fig. 2. The transport regime on a global scale, based on the relationship between the elapsed time and the squared displacement of each element. Note that the axes are logarithmically scaled. 
Unlike the ballistic and diffusion cases, a fractal differential equation is necessary to include these nonFickian processes. Please see Bakunin (2008) for a more systematic and mathematical treatment.
Thus, we can see that the relationship between the elapsed time and squared displacement contains important information about the form of the global transport, although it is only a necessary condition. We investigate this relationship in the transport of magnetic flux concentration on the solar surface.
3. Instruments and data
For globalscale transport, the most important point is the duration of the observation. Because the observational duration of a ground observatory is limited by the Earth’s rotation, it is not suitable for this kind of analysis. We use satellite observations in this study. The second important point is spatial resolution. It is still difficult to detect small magnetic flux concentrations from recent satellite data although their spatial resolution has been improved. Higher spatial resolution makes feature recognition much easier. For these reasons, we select the magnetograms obtained by the Solar Optical Telescope (SOT) onboard the Hinode satellite (Kosugi et al. 2007).
The Hinode satellite was launched in September 2006. SOT is one of the three telescopes onboard the Hinode. It observes the Sun in the visible spectrum and provides spectropolarimetric data, from which we can calculate the magnetic information on the solar surface. The SOT has a high spatial resolution of 0.2–0.3″, although it cannot cover the entire Sun. There are two instruments on SOT: a filtergram (FG) and a spectropolarimeter (SP). We use the magnetogram obtained by the FG in this study. The FG consists of two filtergram imagers: the broadband filtergram imager (BFI) and the narrowband filtergram imager (NFI). NFI obtains the Stokespolarization signal in several lines, and the magnetic field information on the solar surface is calculated from the polarization signal based on the Zeeman effect. We use the magnetograms of the Na I 589.6 nm absorption line obtained by NFI.
The longest observation for which magnetogram data is available, taken between September 2006 and June 2014, is used in this study. SOT observed the quiet region from 10:24 UT on December 30th, 2008 to January 5th, 2009. The total duration of the observation was 115 h 13 min, which is ∼5 days. The temporal cadence and field of view were limited to ∼5 min and ∼110″ × 110″ = 6 × 10^{9} km^{2}, respectively. 1,642 magnetograms were obtained in total. The pixel scale was 0.16″ ∼ 120 km.
Data corrections are needed before the analysis. We performed the same correction as Iida et al. (2012). Dark and flatfield corrections of the CCD camera are done with fg_prep.pro, contained in the SolarSoftWare (SSW) package. The columnwise medianoffset of the CCD camera is known and the median value of each column is subtracted from all pixels in the same column (Lamb et al. 2010). Due to solar rotation, the position of the observing center moves from solar coordinates of (−520.9″, −9.3″) to (695.9″,−6.9″). All images are derotated to 8:01 UT January 2nd, 2009, which is the middle of the observing period, by the procedure drot_map.pro in the SSW package. Although the pointing stability of Hinode is excellent, thanks to a correlation tracker, a small discrepancy still exists. We investigate this discrepancy by determining a correlation between the consecutive images. In this long data set, the discrepancy reaches a value of 15″ ∼ 10^{4} km. Only the region observed over the whole period is analyzed, which has an angular size of 94.9″ × 91.4″. We need the conversion factor between the polarization signal and the actual magnetic field, although a preliminary conversion is done onboard. This conversion factor is derived from the linear fitting between the circular polarization (CP) of SOT/NFI and the lineofsight magnetic field derived from the Michelson Doppler Imager (MDI) onboard the Solar and Heliospheric Observatory (SoHO). The factor is derived as 1.83 G DN^{−1}. Figure 3 shows an example of the magnetogram after all the above corrections.
Fig. 3. An example of a magnetogram after the preprocessing described in the text. 
4. Featuretracking algorithm
We use the same algorithm for the recognition and tracking of the magnetic concentration as Iida et al. (2012). We briefly summarize their method here.
In the recognition process, we use the clumping method by Parnell et al. (2009), where we employ a signal threshold and the pixels above it are recognized as valid pixels. The threshold is determined in each magnetogram by fitting the histogram of the magnetic field strength with a Gaussian function. For our purposes, we employ two σ as the threshold. The σ varies around 5 G in this dataset, which is a typical value for quiet Sun regions, and hence we set ∼10 G as a magnetic field strength threshold in this study. The size threshold, which is used as the smallest size of the recognized concentration, is set to be 81 pixels, or ∼5 × 10^{5} km^{2}, corresponding to the granular size.
After the recognition, the concentrations between the consecutive magnetograms are tied by comparing the spatial overlaps, which are checked with an extra 5 pixel margin. This is necessary because of the 5min interval between the magnetograms. In some cases, more than one concentration overlaps with that in the previous magnetogram. The concentration that has the most similar flux content is related to the previous concentration in these cases. Note also that tracking is done from the concentration with a larger flux content. This is because the smaller concentrations may have a tendency to disappear or stay smaller than concentrations with a larger flux content. With the tracking method explained above, the concentration with the largest flux content survives after coalescence and splitting.
Figure 4 shows the schematic picture of the method. In this case, a concentration of negative polarity starts from the bottom left of the figure. It moves and splits into two concentrations at a certain time. During this splitting, the concentration with larger flux content in the latter magnetogram is treated as the same as the originally tracked concentration. The concentration with smaller flux content is treated as one that has newly emerged. After the splitting, the originally tracked concentration collides with one having the opposite polarity. In this case, three concentrations exist in the tracking method, which are shown by the orange arrows in Figure 4. The lifetime and displacement from birth to death for each concentration are recorded.
Fig. 4. Schematic picture of tracking and squared displacement of the magnetic flux concentration. Three concentrations are tracked. The orange arrows show the displacements for each concentration. Three sets of the elapsed time and the squared displacement are obtained in this example. 
After tracking all concentrations, the elapsed time and squared displacement are obtained for each. Figure 5 summarizes all possible cases of the birth and death of the magnetic flux concentrations in the analysis. The concentrations surrounded by the green square are the newly emerged or vanished ones in the left and right columns, respectively.
Fig. 5. Schematic pictures of births and deaths of magnetic concentrations in our analysis. The left and right columns show the birth and death events, respectively. The concentration surrounded by the green square is that which is produced or dies in each case. 
5. Results
In total, 21,823 positivepolarity and 19,544 negativepolarity concentrations are tracked. First, we investigated the typical lifetime of the magnetic concentrations. The average lifetimes are estimated to be 22.5 and 23.9 min for the positive and negativepolarity concentrations, respectively. Further, the lifetime distribution is investigated. Figure 6 shows the number histograms of the lifetimes, using logarithmic axes. The thin solid and dashed histograms, respectively, correspond to positive and negative concentrations. The thick solid histogram shows the total of both polarities.
Fig. 6. Number histograms of the lifetime. The thin solid and dashed histograms show those of positive and negative polarities, respectively. The thick solid one shows the total of both polarities. 
Next, we investigate the relationship between the elapsed time and the squared displacement of the magnetic concentrations, as displayed in Figure 7. The horizontal axis is the lifetime from birth to death of the concentrations and the vertical axis is the squared displacement. The bins are set [10^{x}, 10^{x+0.1}]. Diamonds indicate the average squared displacement in each bin; bars indicate the standard deviations in each bin as an error of this analysis. We see the change in the dependence around scales of 2 × 10^{4} s (corresponding to a squared displacement of 10^{7.5} km^{2}, or a length scale of 10^{3.8} km, and hence we fit the results into two domains. The dashed and dotted lines show the results of fitting in the range below and above 2 × 10^{4} s, respectively. The dashed one has a powerlaw index of 1.4 ± 0.1 (superdiffusion scaling) and the dotted one has a powerlaw index of 0.6 ± 0.2 (subdiffusion scaling).
Fig. 7. Relationship between the elapsed time and the squared displacement of the patches. Diamonds show the average of the squared displacement and bars show the onesigma error in each bin. The dashed and dotted lines show the fitting results below and above the time scale of 2 × 10^{4} s, respectively. 
Recent papers have investigated the separation distance of the paired magnetic flux concentrations (Lepreti et al. 2012; Giannattasio et al. 2014a). Here, we perform the same analysis on our dataset. Figure 8 shows the dependence of the squared separation distance on the elapsed time, again with logarithmic axes. Diamonds indicate the average squared separation in each bin. We see a linear relationship in this plot and hence fit the result with a line with the range below 10^{5} s. The result of the fitting is shown by the dashed line. We obtain the powerlaw index of 1.46 with the fitting. This value is consistent with recent papers. Lepreti et al. (2012) reported an index of 1.5 over spatial scales from 10^{1.2} km to 10^{2.5} km and temporal scales from 10 s to 400 s. Giannattasio et al. (2014a) reported 1.55 for a quiet region over spatial scales from 10^{2} km to 10^{3.6} km and temporal scales from 10 s to 10^{4} s. Our result extends their result to larger temporal and spatial scales.
Fig. 8. Relationship between the elapsed time and the squared separation distance of the paired patches. Diamonds show the average of the squared separation. The dashed line corresponds to the powerlaw fit. 
6. Discussion
We have investigated the relationship between the elapsed time and the squared displacement of magnetic flux concentrations over the longest observation by the Hinode satellite and found a different behavior above and below the scale of 10^{3.8} km. On short temporal and spatial scales, the superdiffusion regime has an index of 1.4. This scaling is consistent with those of previous studies (Lawrence et al. 2001; Abramenko et al. 2011; Giannattasio et al. 2014b), where the powerlaw index varies from 1.27 to 1.67.
On the other hand, we have newly found in the quiet Sun that the scaling becomes subdiffusive above the scale of 10^{3.8} km. Lawrence & Schrijver (1993) reported subdiffusion scaling in the active region for temporal scales longer than half a day, and spatial scales larger than 6 × 10^{3} km. These scales are similar to the ones found here.
What causes this change of the scaling near 10^{3.8} km? We suggest the supergranulation as a plausible candidate; although the strongest convective flow pattern is a granular flow as mentioned, the magnetic field also forms magnetic network due to being transported by supergranular flow, whose mechanism remains unknown. The typical scale of the network field is ∼10^{4} km, which is slightly larger than that of the change found in this paper. Our interpretation is as follows: inside the network boundary, the magnetic concentration is affected not only by the granular flow patterns but also by the networkflow pattern. The granular flow pattern is expected to cause a diffusivelike transport of the magnetic field because its characteristic size is smaller than that of the network. On the other hand, the networkflow pattern is expected to transport magnetic field ballistically below the scale of 10^{4} km because its cell size is larger; hence, it is natural to assume ballistic transport by the networkflow pattern. With diffusive transport by the granulation and ballistic transport by the networkflow pattern, we expect an intermediate scaling between them, namely superdiffusive transport. In a larger scale than the network field, it is expected that the magnetic concentration will reach the conjunction point of the network. At this point, the magnetic field will be trapped by the network flow. This trapping is expected to result in the subdiffusion scaling above 10^{4} km. To support the hypothesis that the magnetic field is trapped by network downward flows, dopplergrams can be used to identify the cospatiality of the stagnation points and intense downflows in the future study. Figure 9 shows a schematic view of this speculation. Numerical simulation is also a powerful tool to resolve this matter.
Fig. 9. Schematic of the interpretation of the magnetic concentration transport on the network boundary. In and on the network boundary, the motion of the magnetic concentration behaves as a superdiffusion by the network flow and the granular diffusive flow. When it reaches at the conjunction point of the network, it is trapped by the network field and the motion becomes subdiffusive. 
A recent paper by Giannattasio et al. (2014b) found a change in the powerlaw index from 1.34 below the scale of granules (1.5 × 10^{3} km) to 1.25 above it. The spatial scale is smaller by one order of magnitude than the findings in this paper and it remains in the superdiffusion region. The reason we cannot see this may be that the temporal resolution is poorer in the dataset of this study than in theirs. The time interval of our dataset is 5 min, which is insufficient for investigating granular scales of ∼10 min. Thus, it is difficult to see the change in the dependence around the granular scale. Why does the superdiffusion dependence continue at temporal scales longer than granulation? We believe the answer lies in the different balance between the magnetic and kinetic energies of convection. The conjunction point of the network is known as a place where strong magnetic fields exist; it is reasonable to expect such magnetic fields to prevent convective flow from prevailing. On the other hand, the magnetic field at the edge of a granule is weaker and may not be able to stop the convective flow. We believe that this difference may cause the difference in their behaviors.
Next, we will discuss the inconsistency with the theoretical modeling. The requirement for diffusion coefficient from kinetic dynamo simulation is 600 km^{2} s^{−1} (Wang et al. 1991); in previous papers, this value has been reached approximately at the network scale (∼10^{4} km). If the dependence of squared travel distance on time has subdiffusion scaling on longer scales, the effective diffusion coefficient should decrease at such large scales. Considering the 1.5 or 2 orders of spatial difference between the network and the solar global scale, the effective diffusion coefficient becomes less than the requirement by one order of the magnitude from the present result.
Although the subdiffusion scaling obtained in this study covers only one order of magnitude, the investigation of longer time scales is difficult because Hinode has a geocentric orbit and the observational duration is limited by the Sun’s rotation. We need a heliocentric satellite for longer time scales. However, the accuracy of the result can be improved with the accumulation of the dataset. The dependences among quantities such as latitude, longitude, and magnetic flux amount are an interesting question for future research.
At the last part of the Discussion, we shall suggest plausible processes to fill the gap between the measured value of the diffusion coefficient and that required by numerical simulation. We assumed rigid particles without interactions and evaluated the diffusion coefficient of the horizontal motion on the plane, e.g., the photosphere. Some processes are not considered in this modeling. One is the merging and splitting of the magnetic flux concentrations. When magnetic concentrations merge together and then split, the magnetic flux is transported. Hence, the merging and consequent splitting increase the transport coefficient. Because recent studies reveal that merging and splitting of magnetic flux concentrations frequently occur on the actual solar surface (Lamb et al. 2008, 2010; Iida et al. 2012; Gošić et al. 2014), we expect that this process is dominant in global transport. Iida et al. (2015) investigated the same dataset with this study and found the frequent merging and splitting of magnetic concentrations with a time scale of ∼30 min. Other processes that are not considered in this study are cancellation and emergence. The magnetic field is transported in the vertical direction by these processes, making it possible for the magnetic field to go through a conjunction point, increasing the effective transport coefficient. Again Iida et al. (2015) investigated cancellation events in this dataset but found them less frequent than merging and splitting. However, the dependence of cancellation occurrence on magnetic flux has a powerlaw index of −2.48 ± 0.26. The powerlaw index less than −2 implies that smaller cancellation is important in terms of total amount of magnetic flux transport. Hence there is a possibility that cancellation events smaller than those investigated in Iida et al. (2015) transport significant magnetic flux. We expect either or both of these processes to fill the gap between the theoretical requirement and the result in this study. However, there have been no theoretical models taking these processes into account so far. The matter of including them in the global transport equation is left for future work.
Acknowledgments
First of all, we would also like to express our thanks to the Hinode team for providing the continuous observational data. Hinode is a Japanese mission developed and launched by ISAS/JAXA, collaborating with NAOJ as a domestic partner, NASA and STFC (UK) as international partners. Scientific operation of the Hinode mission is conducted by the Hinode science team 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 (U.K.), NASA, ESA, and NSC (Norway). The author would also like to thank Enago (www.enago.jp) for the English language review. The editor thanks two anonymous referees for their assistance in evaluating this paper.
References
 Abramenko, V.I., V. Carbone, V. Yurchyshyn, P.R. Goode, R.F. Stein, F. Lepreti, V. Capparelli, and A. Vecchio. Turbulent diffusion in the photosphere as derived from photospheric bright point motion. Astrophys. J., 743, 133, 2011. [Google Scholar]
 Babcock, H.W. The topology of the Sun’s magnetic field and the 22year cycle. Astrophys. J., 133, 572, 1961. [Google Scholar]
 Bakunin, O.G. Diffusion Turbulence: Scaling Versus Equations, Springer Verlag, Berlin, 2008. [Google Scholar]
 Balescu, R. Transport Processes in Plasmas: Classical Transport Theory, Elsevier Science Ltd., Amsterdam, 1988. [Google Scholar]
 Beck, J.G., L. Gizon, and T.L. Duvall Jr. A new component of solar dynamics: northsouth diverging flows migrating toward the equator with an 11 year period. Astrophys. J., 575, L47, 2002. [CrossRef] [Google Scholar]
 Berrilli, F., S. Scardigli, and S. Giordano. Multiscale magnetic underdense regions on the solar surface: Granular and mesogranular scales. Sol. Phys., 282, 379, 2013. [NASA ADS] [CrossRef] [Google Scholar]
 Caroli, A., F. Giannattasio, M. Fanfoni, D. Del Moro, G. Consolini, and F. Berrilli. Turbulent convective flows in the solar photospheric plasma. J. Plasma Phys., 81, 5, 2015. [Google Scholar]
 DeForest, C.E., H.J. Hagenaar, D.A. Lamb, C.E. Parnell, and B.T. Welsch. Solar magnetic tracking. I. Software comparison and recommended practices. Astrophys. J., 666, 576, 2007. [Google Scholar]
 Giannattasio, F., F. Berrilli, L. Biferale, D. Del Moro, M. Sbragaglia, L. Bellot Rubio, M. Gošić, and D. Orozco Su′arez. Pair separation of magnetic elements in the quiet Sun. A&A, 569, A121, 2014a. [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Giannattasio, F., M. Stangalini, F. Berrilli, D. Del Moro, and L. Bellot Rubio. Diffusion of magnetic elements in a supergranular cell. Astrophys. J., 788, 137, 2014b. [CrossRef] [Google Scholar]
 Gošić, M., L.R. Bellot Rubio, D. Orozco Suárez, Y. Katsukawa, and J.C. del Toro Iniesta. The solar internetwork. I. Contribution to the network magnetic flux. Astrophys. J., 797, 49, 2014. [CrossRef] [Google Scholar]
 Hagenaar, H.J., C.J. Schrijver, A.M. Title, and R.A. Shine. Dispersal of magnetic flux in the quiet solar photosphere. Astrophys. J., 511, 932, 1999. [Google Scholar]
 Iida, Y., H.J. Hagenaar, and T. Yokoyama. Detection of flux emergence, splitting, merging, and cancellation of network field. I. Splitting and merging. Astrophys. J., 752, 149, 2012. [Google Scholar]
 Iida, Y., H.J. Hagenaar, and T. Yokoyama. Detection of flux emergence, splitting, merging, and cancellation of network fields. II. Apparent unipolar flux change and cancellation. Astrophys. J., 814, 134, 2015. [Google Scholar]
 Kosugi, T., K. Matsuzaki, T. Sakao, T. Shimizu, Y. Sone, et al. The Hinode (SolarB) mission: an overview. Sol. Phys., 243, 3, 2007. [NASA ADS] [CrossRef] [Google Scholar]
 Lamb, D.A., C.E. DeForest, H.J. Hagenaar, C.E. Parnell, and B.T. Welsch. Solar magnetic tracking. II. The apparent unipolar origin of quietsun flux. Astrophys. J., 674, 520, 2008. [Google Scholar]
 Lamb, D.A., C.E. DeForest, H.J. Hagenaar, C.E. Parnell, and B.T. Welsch. Solar magnetic tracking.III. Apparent unipolar flux emergence in highresolution observations. Astrophys. J., 720, 1405, 2010. [Google Scholar]
 Lawrence, J.K., and C.J. Schrijver. Anomalous diffusion of magnetic elements across the solar surface. Astrophys. J., 411, 402, 1993. [Google Scholar]
 Lawrence, J.K., A.C. Cadavid, A. Ruzmaikin, and T.E. Berger. Spatiotemporal scaling of solar surface flows. Phys. Rev. Lett., 86, 5894, 2001. [NASA ADS] [CrossRef] [Google Scholar]
 Leighton, R.B. Transport of magnetic fields on the Sun. Astrophys. J., 140, 1547, 1964. [Google Scholar]
 Lepreti, F., V. Carbone, V.I. Abramenko, V. Yurchyshyn, P.R. Goode, V. Capparelli, and A. Vecchio. Turbulent pair dispersion of photospheric bright points. Astrophys. J., 759, L17, 2012. [Google Scholar]
 Nordlund, Å., R.F. Stein, and M. Asplund. Solar surface convection. Living Rev. Sol. Phys., 6, 2, 2009. [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 November, L.J., J. Toomre, K.B. Gebbie, and G.W. Simon. The detection of mesogranulation on the Sun. Astrophys. J., 245, L123, 1981. [Google Scholar]
 Parnell, C.E., C.E. DeForest, H.J. Hagenaar, B.A. Johnston, D.A. Lamb, and B.T. Welsch. A powerlaw distribution of solar magnetic fields over more than five decades in flux. Astrophys. J., 698, 75, 2009. [Google Scholar]
 Rieutord, M., and F. Rincon. The Sun’s supergranulation. Living Rev. Sol. Phys., 7, 2, 2010. [CrossRef] [Google Scholar]
 Schrijver, C.J., and S.F. Martin. Properties of the large and smallscale flow patterns in and around AR 19824. Sol. Phys., 129, 95, 1990. [NASA ADS] [CrossRef] [Google Scholar]
 Thornton, L.M., and C.E. Parnell. Smallscale flux emergence observed using Hinode/SOT. Sol. Phys., 269, 13, 2011. [NASA ADS] [CrossRef] [Google Scholar]
 Wang, Y.M., A.G. Nash, and N.R. Sheeley Jr. Magnetic flux transport on the Sun. Science, 245, 712, 1989. [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Wang, Y.M., N.R. Sheeley Jr., and A.G. Nash. A new solar cycle model including meridional circulation. Astrophys. J., 383, 431, 1991. [Google Scholar]
 Yelles Chaouche, L., F. MorenoInsertis, V. Martínez Pillet, T. Wiegelmann, J.A. Bonet, et al. Mesogranulation and the solar surface magnetic field distribution. Astrophys. J., 727, L30, 2011. [Google Scholar]
Cite this article as: Iida Y. Tracking of magnetic flux concentrations over a fiveday observation, and an insight into surface magnetic flux transport. J. Space Weather Space Clim., 6, A27, 2016, DOI: 10.1051/swsc/2016022.
All Figures
Fig. 1. Schematic picture of a 1D randomwalk motion. τ and δ are, respectively, temporal and spatial scales for one walk. 

In the text 
Fig. 2. The transport regime on a global scale, based on the relationship between the elapsed time and the squared displacement of each element. Note that the axes are logarithmically scaled. 

In the text 
Fig. 3. An example of a magnetogram after the preprocessing described in the text. 

In the text 
Fig. 4. Schematic picture of tracking and squared displacement of the magnetic flux concentration. Three concentrations are tracked. The orange arrows show the displacements for each concentration. Three sets of the elapsed time and the squared displacement are obtained in this example. 

In the text 
Fig. 5. Schematic pictures of births and deaths of magnetic concentrations in our analysis. The left and right columns show the birth and death events, respectively. The concentration surrounded by the green square is that which is produced or dies in each case. 

In the text 
Fig. 6. Number histograms of the lifetime. The thin solid and dashed histograms show those of positive and negative polarities, respectively. The thick solid one shows the total of both polarities. 

In the text 
Fig. 7. Relationship between the elapsed time and the squared displacement of the patches. Diamonds show the average of the squared displacement and bars show the onesigma error in each bin. The dashed and dotted lines show the fitting results below and above the time scale of 2 × 10^{4} s, respectively. 

In the text 
Fig. 8. Relationship between the elapsed time and the squared separation distance of the paired patches. Diamonds show the average of the squared separation. The dashed line corresponds to the powerlaw fit. 

In the text 
Fig. 9. Schematic of the interpretation of the magnetic concentration transport on the network boundary. In and on the network boundary, the motion of the magnetic concentration behaves as a superdiffusion by the network flow and the granular diffusive flow. When it reaches at the conjunction point of the network, it is trapped by the network field and the motion becomes subdiffusive. 

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.