{'fraction_non_alphanumeric': 0.02719553886678249, 'fraction_numerical': 0.015454239757888247, 'mean_word_length': 4.6663226421086055, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The formation and disintegration of magnetic bright points observed by Sunrise/IMaX

The evolution of the physical parameters of magnetic bright points (MBPs) located in the quiet Sun (mainly in the interwork) during their lifetime is studied. First we concentrate on the detailed description of the magnetic field evolution of three MBPs. This reveals that individual features follow different, generally complex, and rather dynamic scenarios of evolution. Next we apply statistical methods on roughly 200 observed MBP evolutionary tracks. MBPs are found to be formed by the strengthening of an equipartition field patch, which initially exhibits a moderate downflow. During the evolution, strong downdrafts with an average velocity of 2.4 km/s set in. These flows, taken together with the concurrent strengthening of the field, suggest that we are witnessing the occurrence of convective collapses in these features, although only 30% of them reach kG field strengths. This fraction might turn out to be larger when the new 4 m class solar telescopes are operational as observations of MBPs with current state of the art instrumentation could still be suffering from resolution limitations. Finally, when the bright point disappears (although the magnetic field often continues to exist) the magnetic field strength has dropped to the equipartition level and is generally somewhat weaker than at the beginning of the MBP's evolution. Noteworthy is that in about 10% of the cases we observe in the vicinity of the downflows small-scale strong (exceeding 2 km/s) intergranular upflows related spatially and temporally to these downflows.


INTRODUCTION
Magnetic bright points (MBPs) are manifestations of small-scale solar magnetic fields (for a recent review see, e.g., de Wijn et al. 2009) and are counted among the most interesting features on the Sun. Only isolated, small-scale magnetic field concentrations of kG strengths appear bright in the solar photosphere, as was shown by Criscuoli & Uitenbroek (2014) and Riethmüller et al. (2014). But the influence of the telescope PSF leads to a smoothing of the "true" magnetic field strength to lower values, so that recent observations display a distribution of field strengths from small fractions of kG up to 1.5 kG (see, e.g., Beck et al. 2007;Utz et al. 2013b, and references therein). Such concentrations of magnetic field appear as bright point-like features in filtergrams and are termed MBPs (e.g., Wiehr et al. 2004). The increase in brightness can be explained due to a partial evacuation of the magnetic field structure leading to a decreased density within the structure. This decreased density turns into an effective lowering of the optical unity layer showing then deeper and hotter plasma. Besides, it allows more radiation and hence heat to penetrate into the feature from the surrounding hot walls (see e.g. Keller et al. 2004;Vögler et al. 2004).
MBPs have been studied for several decades with the first observations dating back to the 1970s (see, e.g., Dunn & Zirker 1973;Muller 1983). Since then their dynamical (e.g., Berger & Title 1996;de Wijn et al. 2005;Chitta et al. 2012;Jafarzadeh et al. 2013) as well as their statistical and magnetic properties (e.g., Beck et al. 2007;Bovelet & Wiehr 2008;Utz et al. 2013a;Giannattasio et al. 2013) have been characterised in increasing detail. In recent years such filtergram studies were done mostly in the Fraunhofer G-band (Berger et al. 1998(Berger et al. , 2004Bovelet & Wiehr 2008) but also in other molecular bands such as the CN band-head (Zakharov et al. 2005) or the TiO (Abramenko et al. 2010). Besides the analysis of pure observational data also the comparison between characteristics in observations and simulations attract more and more interest (e.g., Keys et al. 2011). Moreover, spectropolarimetric studies have contributed to our knowledge of MBPs (e.g., Bellot Rubio et al. 2001;Nagata et al. 2008;Viticchié et al. 2010). Lately, combinations of filtergram (e.g., G-band and white light) data and spectropolarimetric observations have been performed (see, e.g., Riethmüller et al. 2010;Romano et al. 2012).
The interest in MBPs arises because they are easy to observe proxies of kG magnetic flux concentrations. These in turn are of great interest for several reasons, such as the not yet completely understood chromospheric and coronal heating problem (e.g., Walsh & Ireland 2003;Klimchuk 2006). By connecting the solar interior to the outer solar atmosphere magnetic flux concentrations play an important role in most proposed heating mechanisms, irrespectively of whether the actual heating takes place via magnetic reconnection, ohmic dissipation, or wave heating processes. Another important aspect is the contribution of MBPs to variations in total solar irradiance, which may influence the Earth's climate (see Solanki et al. 2013, for a review). Furthermore, they influence the granulation itself and the energy transport. A number of authors including Title et al. (1989), Muller et al. (1989), Narayan & Scharmer (2010), and Andić et al. (2011) have reported on smaller granular cells around ensembles of magnetic flux concentrations and concluded that the presence of strong magnetic fields suppresses to a certain extent normal convective flows. Furthermore there is recently a debate going on about a possible subpopulation of granules -so-called mini granules (e.g., Abramenko et al. 2012;Lemmerer et al. 2014). Hence the question arises if and how this kind of granules is related to small-scale magnetic fields? Studies of the evolution of the number of network magnetic elements and MBPs in the quiet Sun over the solar cycle have given results that depend on the magnetic flux per feature and on the length of time considered (see Hagenaar et al. 2003). Another recent attempt (for the last solar minimum) of tracking the number and variation of the number of individual MBPs was, unfortunately, not conclusive (see Muller et al. 2011).
Since kG fields are needed to produce MBPs, their appearance, disappearance, and evolution teach us how the vigorous photospheric convection interacts with magnetic fields, strengthens them and keeps them in equilibrium until they dissolve or disintegrate. The most feasible theory for the creation of small-scale kG magnetic features is that of convective collapse proposed by Parker (1978) or similar by Spruit (1979). The convective collapse model is based on the idea that a concentrated magnetic field suppresses efficiently the transport of energy by convection. Therefore the plasma locally cools down and a downdraft sets in. This downdraft evacuates the flux tube, so that the excess pressure of the surrounding gas compresses it (or it "collapses"). After the collapse, a new stable state can emerge (depending on the initial conditions) with magnetic field strengths above 1 kG. The convective instability acting in small-scale magnetic features and their subsequent concentration have been verified observationally (e.g., Bellot Rubio et al. 2001;Nagata et al. 2008;Narayan 2011;Requerey et al. 2014) as well as by numerical experiments (see, e.g, Grossmann-Doerth et al. 1998;Danilovic et al. 2010;Hewitt et al. 2014).
While the convective collapse model serves us well in understanding the processes leading to kG flux tubes, much less is known about the processes leading to the dissolution and disintegration of the strong kG tubes. The aforementioned work of Bellot Rubio et al. (2001) states that disintegration might be caused by upward moving shocks in the flux tubes produced by fast downflowing material hitting the base of the flux tube and bouncing back. Such a mechanism was observed in simulations by Steiner et al. (1998). This scenario has been studied by computer simulations such as the one of Takeuchi (1999). Another possibility for the disintegration of thin magnetic flux tubes is the interchange instability, which can be overcome by whirl flows around the flux tube (see Parker 1975;Schüssler 1984). Alternatively, the magnetic field may leak out of the features by diffusion, or may be torn out of them by turbulent convection, or may reconnect with turbulent fields in the immediate surroundings of the magnetic features. Hence investigating the evolution of the plasma parameters of magnetic flux concentrations during the very end of their lifetime is of great interest. How is the flux concentration maintained and stabilized during their life, and why and how does it dissolve?
To shed more light on these questions it is necessary to track individual magnetic flux concentrations during their lifetime and obtain as detailed and highly resolved information about the interesting quantities as possible. This is the main aim of the present study, where we employ MBPs as proxies of kG magnetic concentrations. These are used to distinguish such concentrations from other, weaker field, features. In the first part a few specific evolutionary tracks will be discussed in detail while later on a statistical approach is applied to all of the roughly 200 tracked MBPs. This work complements the recent work of Requerey et al. (2014) in which the evolution of a single magnetic patch was studied in full detail and in addition the study of Hewitt et al. (2014), who investigated in detail the formation, evolution, and dissolution of MBP features in MHD simulations.
2. DATA The data were acquired with the IMaX instrument (see Martínez Pillet et al. 2011) on board the Sunrise balloon-borne mission Barthol et al. 2011;Berkefeld et al. 2011) which was flown from 2009 June 8 to 13 along a route from ESRANGE/Kiruna (balloon launching facility in northern Sweden) to Somerset island (northern Canada). During this period it ascended to the stratosphere to a height of about 37 km above the ground where it took practically seeing-free images of the Sun. The IMaX instrument is an imaging magnetograph making use of the Fe I 5250.2Å line, which is the most magnetically sensitive line in the visible with an effective Landé factor of three. According to its instrumental paper (Martínez Pillet et al. 2011): "Gauss equivalent sensitivities are 4 G for longitudinal fields and 80 G for transverse fields per wavelength sample. The lineof-sight (LOS) velocities are estimated with statistical errors of the order of 5-40 ms −1 ." This achieved performance of the instrument and used spectral line makes it the best bet to detect weak fields in quiet Sun when one does not have access to the near infrared part of the spectrum. In operation it takes at each of the N λ wavelengths, N p differently polarised (by using liquid crystal variable retarders as modulators) single exposures that are repeated N A times in order to increase the signalto-noise ratio (S/N). Therefore, a single line scan data set consists of N λ × N p images that are taken in a time period equal to N λ × N p × N A × t exp , where t exp is the exposure time of a single shot plus overhead for readout, etc. (about 250 ms).
The analysed data consist of 100 such single data sets obtained by the IMaX instrument in the so-called V5-6 mode. In this mode, V means that N p = 4 (vector mode), and 5 and 6 stand for N λ and N A , respectively. The Fe I spectral line at 5250.2Å was sampled at −80 mÅ, −40 mÅ, 40 mÅ, 80 mÅ and 227 mÅ off the line centre. The last sample is located in the nearby continuum. With this sampling scheme, the V5-6 mode reaches a typical S/N of 1000 and consists of all four Stokes parameters (I, Q, U , V ).
The obtained sequence is broken in two parts: the first, lasting from 0:35:49 UT until 0:58:32 UT and consisting of 42 sets of images taken about every 33 s, and the second, with the same temporal resolution taken from 1:30:41 UT until 2:02:16 UT comprising 58 sets of Stokes profile maps. Both were acquired on 2009 June 9 close to the solar disc centre under exceptional quiet Sun conditions. The spatial sampling was 0.055 arcsec per pixel. Considering the detector area of 936 × 936 pixels 2 , the field of view (FOV) is 51.48 × 51.48 arcsec 2 . Due to the applied image reconstruction algorithms (resulting in apodised images; see Martínez Pillet et al. 2011 for details) the useful size of the FOV was reduced to about 780 × 780 pixels 2 . The application of the phase diversity reconstruction algorithm yielded images with a typical signal to noise ratio of about 350 in the Stokes I continuum and a spatial resolution of about 0.15 -0.18 arcsec.
A snapshot of a part of the field of view and some of the data products can be seen in Fig. 1. From left to right and from top to bottom, we display images of the continuum intensity, of the blue-wing intensity (−40 mÅ), of the total absolute circular polarization, of the magnetic field strength and inclination, of the LOS velocity, and of the temperature at log τ = 0, −1, and −2. The six latter images come from inversions with the SIR code (Stokes Inversion based on Response functions; Ruiz Cobo & del Toro Iniesta 1992).

ANALYSIS
We used an already reduced data set publicly available on the mission Web page. The full reduction comprises flat fielding, dark current reduction, polarimetric demodulation (necessary to retrieve the Stokes parameters from the measured linear combinations), cross-talk removal, and image reconstruction by using the point-spread function of the instrument as obtained from specifically devoted phase-diversity observations. All the detailed information about data reduction and reconstruction can be found in the instrument paper of Martínez Pillet et al. (2011) and references therein. To stay as close as possible to the original data we have abstained from p-mode filtering. Hence, the p-modes may have an influence on the measured LOS velocities.
The inversions were carried out using two nodes for the temperature (i.e., modifying the initial guess atmosphere linearly) and one node for magnetic field strength, inclination, azimuth, LOS velocity, and microturbulence (assumed to be height independent). Due to the high spatial resolution, a zero macroturbulence velocity and a magnetic filling factor of unity were assumed. This means that for unresolved structures, the parameters returned by the inversion might represent averages over several atmospheric components. In total, seven free parameters were determined from the 5 × 4 = 20 observables. The average LOS velocity in the complete FOV was set to zero. Therefore LOS velocity values in this study have to be regarded as difference values to the average (quiet Sun) LOS velocity.
Several different algorithms for small-scale feature identification were developed in the past and presented in literature among with results on MBPs. To give at least two examples of automated and more sophisticated algorithms: a "blob finding" algorithm was successfully applied by Berger et al. (1995)   algorithm based on the idea of identifying MBPs within intergranular lanes. In a second step they use the fact that MBPs should be surrounded by pixels belonging to intergranular lanes and hence a "compass" search (investigation of the surrounding of a pixel) should yield only dark intergranular pixels. In addition other known properties of MBPs like their brightness and brightness gradient can be used.
The MBP identification in this study was done by using the automated procedure developed by Utz et al. (2009Utz et al. ( , 2010 on the blue-wing wavelength sample maps (−40 mÅ from the line core). The high contrast of these images and the relatively high formation region of the line at that wavelength makes it ideal for MBP detection among the available IMaX data. The interesting features are situated in the intergranular lanes and are commonly associated with downflows, these downflows might further facilitate the detection by enhancing the contrast to the surroundings via line shifts (shifting the Fe line to longer wavelengths and thus giving larger measured intensities in the line wing). In contrast, the direct usage of the continuum intensity maps leads to identification difficulties, with a significant fraction of wrongly identified non-magnetic features (Utz et al. 2013c). We have refrained from using the high-contrast images provided by the Sunrise Filter Imager (SuFI Gandorfer et al. 2011), which are not affected by downflows and display MBPs clearly, because of its limited FOV, which would have restricted the statistical sample at our disposal too strongly.
After the identification of the MBPs the tracking can be done (for an overview and in depth discussion of several available algorithms; see DeForest et al. 2007). The identification algorithm leaves us with a set of x and y positions of the intensity barycentre of the identified features (plus parameters such as brightness or size) for every image. The tracking procedure has now the purpose of connecting these single positional measurements of individual MBPs detected in each image to generate tracks or time-series of their evolutions. Such an algorithm must consider three constraints: a) the distance of connected features from one image to the next has to be smaller than a maximum distance given as v max × ∆t where v max is the maximum horizontal velocity of the MBP (about 4 km/s; see, e.g., Nisenson et al. 2003, Jafarzadeh et al. 2013 who found MBPs to occasionally reach velocities of up to 16 km/s) and ∆t is the temporal resolution of the data set (in our case 33 s); b) the set of connections between the feature positions should have a minimum total length, where the total length of a possible set of connections is defined as the sum of the distances between connected features; and c) no two paths of MBPs should cross each other from one time step to the next one as this would mean that one would have "tunneled" through the other one.
The approach to the tracking problem adopted in this paper does not describe the transition from one image to the next (executing a direct feature track assignment for every MBP), but instead, adds up additional pointer information. These pointers drive the programme to a follow-up and predecessor realisation. To obtain the pointers, we compare the current feature positions with the positions in the previous and subsequent image. The spatially closest following and preceding structure, regarding also the maximum distance rule, gets assigned to the MBP in question via two pointer variables which are added to the data structure. This structure contains then all the information of the single measurements, such as a unique identifier number (for the single measurement), brightness, size and so on, plus the pointer information to the follow up measurement and preceding one. After this tracking, the creation of time series of MBPs is a simple task. All features which have no pointer assigned to a previous realisation are obviously the starting points of MBP sequences. These structures can be taken and the evolution can be followed by jumping via the pointers from one measurement to the next.
The retrieval of discussed physical quantities is now easily done by taking the plasma parameters at the barycentre position of the MBP. For a better comparison with the surrounding area and since the barycentre position did not always agree with the brightest or strongest pixel, we also computed the mean, maximum and minimum values in the vicinity of the feature's barycentre positions (keeping in mind that those evaluated positions must not coincide with the same pixel in all instants of time). The vicinity is defined as a square of 5 × 5 pixels 2 centered around the pixel of interest.

RESULTS
Prior to a statistical study of some 200 evolutionary tracks, we start this section with a detailed discussion of some specific, even though not classifying, cases (i.e., that each MBP track is too distinct from each other to be assigned to a specific class or type of MBP evolution). The first shown example is fairly representative of a canonical convective collapse. This theoretical concept is believed to represent an important process leading to the magnetic field amplification for small-scale magnetic structures, but it may not be the only possible amplification process. The second case illustrates that instead of the presence of a downflow per se, the interaction of the strong downflows with nearby small-scale upflows may be responsible for the magnetic field amplification, in agreement with recent findings by Requerey et al. (2014). In the third example we show an evolution in which plasma flows play only a minor role, if at all, for the evolution and creation of a strong magnetic field concentration.

Case I: double convective collapse event
In Fig. 2 we show the tracked feature in the plasma parameter maps. The first row gives the blue line-wing sample. The brightening is clearly seen in the maps and can still be visually identified in the first two images after the automated tracking has lost the MBP. This is probably due to a better human flexibility and ability in identifying structures and features compared to the employed automated algorithm. The next row illustrates the continuum sample. The MBP is formed at the merging site of several granules. The same feature can be identified also in the map of temperature at log τ = −2 shown in the third row. The next row gives the magnetic field strength map. Here, the MBP is formed in a region of weak magnetic field that first strengthens then weakens again, and persists after the identification ceases. In the LOS velocity maps (fifth row) we see a strong downflow appearing at 66 s and again at 132 s. To probe the details of the MBP's evolution, we display in Fig. 3 plots of the physical parameters within a square of 5 × 5 pixels centred on the barycentre of the structure in the blue wing maps. The evolution of the barycentre is plotted in black. Orange lines correspond to these pixels presenting maximum values of the corresponding quantity in the selected area. Blue lines correspond to the pixel having the minimum value. Note that lines with the same colour may correspond to different pixels in each panel. Finally, dashed, green lines give the evolution of the average in the area. Error bars indicate the rms dispersion of each quantity within the selected box. The dashed vertical lines mark the beginning and the end of the automated MBP tracking.
Noteworthy, the brightness barycentre does not coincide spatially with the strongest field pixel, neither with the fastest downflow one. Except for a couple of instants where the latter is almost located at the brightness barycentre it seems clear that the maximum intensity and the largest velocities are located in the surroundings of the strongest field (tentatively, the magnetic structure's core). This is in agreement with the results found by Requerey et al. (2014), Hewitt et al. (2014), and the ones outlined in the thesis of Bühler (2013) using a new inversion technique developed by van Noort (2012). The authors of the latter analysed a plage region and came to a similar conclusion, namely that the downflows avoid the points of the strongest magnetic field.
Before the automated tracking starts it is clear that a magnetic field with a strength of a small fraction of a kG already exists. About 1 minute later, a first strong downflow of 3 kms −1 sets in that is accompanied by a clear strengthening of the magnetic field up to some 1.4 kG. Both, the downflow and the magnetic field strength, significantly decrease in just 30 s to 1 kms −1 and 0.5 kG just to rapidly increase again in the next half-minute. This time the field strength reaches 1.3 kG (0.9 kG for the brightness barycentre) and the downflow accelerates to almost 5 kms −1 . Such rapid variations and fluctuations of plasma parameters were already reported previously by Jurčák et al. (2013). While we see in these observations the magnetic field and LOS velocity maxima to occur co-temporal it was reported by Hewitt et al. (2014) that at the point of maximum field strength the LOS velocity cedes. The difference in observations and simulations might be explained by the lower temporal resolution in observations (about 30 s) compared to simulations (about 2 s).
The size varies in a very nicely anti-correlated fashion with the field strength, as expected from the conservation of magnetic flux. At τ = 1, the temperature does not show clear correlations relative to the magnetic and dynamic evolution of the MBP (i.e., that the maximum field strength does not vary in phase with the temperature maximum), which is in agreement with the difficulties in MBP selection in the continuum maps. However, in higher layers (mostly at log τ = −2), clear signatures of temperature enhancements can be seen at the same time as the convective collapse episodes happen. These observational findings agree mostly with recent simulations (see Fig. 1, fourth panel, of Hewitt et al. 2014) where an increase in temperature in all three discussed atmospheric layers have been seen.
Still within the tracking period and only 33 s after the second collapse, both, the magnetic field strength and the plasma flow, relax to the values they had prior to the first convective collapse. Therefore, the kG state does not seem to be a stable configuration of this feature that has nevertheless experienced a double convective collapse event.
4.2. Case II: illustration of the importance of shear flows/upflows The maps of the plasma parameters and the evolution of a second MBP are depicted in Fig. 4. While the track of the feature in the blue line-wing, the continuum, the temperature, and also in the magnetic field strength map does not show any special behaviour with respect to the other studied cases, one can see an interesting detail in the velocity maps. Starting in the second image after the automated tracking started, and vanishing again in the fourth image, there is a very small inter-granular upflow in the middle of the downflowing region, very close to the detected barycentre of the MBP. The connection of the upflow with the neighbouring granule and its coincidence with strong downflows, have recently been investigated by Requerey et al. (2014). In addition such small-scale upflows can be seen in simulations like the one of Hewitt et al. (2014, e.g., Fig. 4, middle panel and lower panel) where one can see that the barycentre of the found MBP is situated actually just in between a strong downflow and an upflow.
Interestingly, within the same evolutionary track, we find a second such occurrence in the last two tracked images. Again, a small upflow evolves very close to the studied downflow and disappears just after one more image.
As for the former case, Fig. 5 illustrates the evolution of the brightness barycentre, and of the maximum and minimum quantities as well as the average over 5 by 5 pixels 2 centred on the barycentre. Again, this MBP shows two conspicuous magnetic strengthenings. The strongest pixel reaches 1.2 kG and later 1 kG, while the barycentre itself has some 0.8 kG on both occasions. Peaks in the downflow of 2.5 and 2 kms −1 , respectively, appear simultaneously with the magnetic field strength maxima. Noteworthy is the fact that at the time of the strong downflows, two upflows are detected within the defined vicinity box around the MBP barycentre. These upflows cover only a rather small region of approximately 10.000 km 2 (estimated by visual inspection) and are weaker than the downflows, but strong enough to be reliable with values between −0.5 to −1 kms −1 . These striking coincidences of up-and downflows, while the magnetic field strengthens, point to the idea that both processes might be coupled somehow and that the canonical convective collapse onset could be accompanied by some mechanisms (e.g., Steiner et al. 1998;Steiner 1999;Kato et al. 2011;Requerey et al. 2014). The other plasma parameters show no behaviour worth mentioning except for the evolution of the temperature in the higher photospheric layers. Interestingly, this temperature displays a double peak, as seen for the magnetic field strength and velocity evolution, but only the second peak is co-temporal to the velocity and magnetic field strength peaks.
A manual inspection of all the tracked MBPs show that about 10% of all the tracked features show stronger spatially located upflows close to the downflows related to the magnetic feature. Therefore the discussed case might still be quite a special one but definitely not the only one in which a spatially correlated upflow plays an important role in the evolution of the feature. 4.3. Case III: magnetic field evolution without strong LOS velocities The third presented case complements the picture and shows the complexity and variety of possible evolutionary paths. This case differs from the others by the large extent of the magnetic field patch at whose border the MBP feature was detected, as can be seen from Fig. 6. This patch also harbours a number of further bright points and filigree and is probably a network magnetic field concentration which is a major difference to the other, more or less well isolated, small and compact, MBPs. Figure 7 illustrates the evolution of the different The FOV displayed is about 1.6 arcsec × 1.6 arcsec. The contours in the velocity maps correspond to the intersection between upflow and downflow region while the contour in the magnetic field strength maps outlines significant magnetic fields with strengths above 0.45 kG. The images have been clipped in such a way as to enhance the contrast with the clipping range given by the colourbars.   plasma parameters. The magnetic field strength shows a well defined maximum, which exceeds the kG threshold, 132 s after the start of tracking. Also, the size follows the commonly accepted picture of a fast shrinkage with a co-temporal increase of magnetic field strength (see e.g. Spruit 1979;Fischer et al. 2009). But the most interesting curve for this case is the evolution of the LOS velocity component. There is a strong downflow of about 2.6 kms −1 at the beginning of the tracking witnessed by the orange line (maximum in the FOV) that gets weaker (dropping to below 0.3 kms −1 for the barycentre -black line -and average values -green line) when the magnetic field strength maximum is reached. It increases again with weaker magnetic fields when the tracking is, in fact, stopped. This is in contrast to the other studied cases, where the LOS velocity component behaved more or less co-temporal with the changing magnetic field strength. On the other hand this is in perfect agreement with the simulations done by Hewitt et al. (2014) who stated that the LOS velocity cedes or at least strongly decreases at the moment of intensity and magnetic field maximum. In their simulations the LOS velocity peaked always before the intensity and magnetic field strength reached its maximum. An increase in the temperature of the higher layers can also be seen, although weaker than in the two former cases.
The major difference to the two former cases is probably that this feature was not created out of a weak internetwork magnetic field patch but is on the contrary related to an extended network magnetic field region and second that the formed bright point was never really isolated. Even though it has been partly dislocated and split up from the network patch. Furthermore the final field strength at the point of dissociation is much higher with roughly 0.6 kG. This makes it more compatible to the two cases discussed in the Appendix of the paper of Hewitt et al. (2014), where the MBPs were formed in a magnetic field patch already clearly exceeding at the point of formation the 1 kG threshold. Besides, the tracking in such "network" case is more elaborate as nearby features constantly interact with the feature under investigation. We will go into more details on this in the discussion section.
From the three presented examples, it is clear that while one can learn a lot from single (partly special) cases, it is quite hard to understand the evolution of MBPs on average or in a general way. Different cases seem to be just too distinct from each other to generalise based on such singular observations. To achieve the goal of a more general characterisation we applied statistical methods on all of the tracked and investigated evolutions at hand.

MBP evolutionary statistics
In this subsection we outline the major statistical results derived from the roughly 200 single MBP tracks, incorporating a few thousand single MBP features identified in the course of this study. For this purpose, we have sorted the tracks with increasing initial magnetic field strength, whereby the field strength prior to the Beginning of tracking t=0 s Posterior to tracking t=198 s  first tracked position is used. In the top panel of Fig. 8 we plot this field strength (black asterisks) together with the maximum field strength reached during the evolution of the tracked MBPs (top of the vertical bars) and their dissolution magnetic field strength (bottom of the vertical bars; that is extracted at the last known position in the first image after the tracking).
MBPs are most likely distinctive features that should have a different magnetic field strength distributionand, hence, also onset and dissolution strengths-than the pure background magnetic field in the FOV due to noise in the polarisation signals. To estimate such a background noise field strength we have averaged that of pixels having Stokes Q, U , and V signals below 3 σ of the continuum intensity and found 0.14±0.08 kG. This noise level is displayed in the top panel of Fig. 8 with a blue, dash-dotted line. Since 97 % of the tracks show at least at one instance of time field strengths larger than 0.38 kG, that is, the average background noise field strength plus three times the rms, we make sure that the overwhelming majority of the automatic MBP detections correspond to distinctive magnetic features. The previously mentioned threshold strength of 0.38 kG is almost coincident with the upper limit (0.4 kG) for the equipartition field strength range (see, e.g., Wiegelmann et al. 2010, and references therein) which is marked with a green, dashed line in the upper panel of the figure. A glance at the figure reveals that the onset field strengths of all the MBPs are distributed around that equipartition value. Indeed, the median of these onsets is 0.42 kG. The vertical bars show us that most of the MBP field strengths stay below a deduced field strength of one kG with an amplification of two to three times their onset value. Most of them dissolve with field strengths below the equipartition value.
The LOS velocity component (second panel of Fig. 8) displays in general downflows, for both, the onset and the dissolution, i.e., there are only a few cases of detected flow reversals at the very end (24 upflows out of more than 200 features resulting in about 10% of measured flow reversals). Important to note is that the statistics is done for the barycentre position and not for the maxima or minima in the vicinity of the MBP (those quantities have been shown previously in the evolutionary track Figs. 3, 5, and 7 in addition to the barycentre), i.e., it is still possible that there are even more upflows in the vicinity. Also, because MBPs are located almost exclusively in downflowing intergranular lanes, it may well be that upflows are hidden by the effects of scattered light. Illustrations on how scattered light can hide small-scale flows of opposite sign to the dominant ones are given by (Joshi et al. 2011;Scharmer et al. 2011) for sunspot penumbrae. The initial phase is characterised by weak downflows on the order of 1 km/s (0.8 ± 0.6 km/s), which is also the case for the dissolution stage when plasma flows are nearly back to the initial ones (only slightly changed to an even smaller average value of 0.6 ± 0.6 km/s). In a few cases it can be observed that the LOS velocity component changes sign, meaning that there is an upflow at the very end of the evolution. In about half of the cases, the downflows only reach a maximum value of about 2 km/s and only a few cases (about 5 %) show very strong downflows exceeding 4 km/s, with the strongest reaching 6-7 km/s.

DISCUSSION
In the following we discuss the various aspects of this work, evaluate it critically and compare the methods and results with the literature. We structure the discussion by breaking it into subsections that consider the various steps of the analysis and the results.

Identification of MBPs in the blue line-wing:
Narrowband filtergrams taken in the G band, which is a CH absorption band centred around 430.5 nm, have become a standard observable in recent MBP studies (see, e.g., Berger et al. 2004;Sánchez Almeida et al. 2010). Since the formation height of the G band is at mid-upper photospheric layers, the closest available proxy within the data set at hand is one of the wavelength samples closer to the line core which is therefore more sensitive to the temperature in such layers. The very high temperature dependence of the Fe I 525.0 nm line strength makes this a sensitive parameter. Hence, we have chosen the blue line-wing sample at −40 mÅ off the line centre for the identification of MBPs. While this facilitates the detection of MBPs as they show most of the time downflows which contributes to an even higher intensity into this line sample (due to line shifts), it bears the risk of masking upflows as those velocities might decrease the brightness in the blue line-wing sample.

Feature tracking:
Two contrasting tracking approaches are illustrated in Fig. 9 on the basis of a simple example covering three time steps, with a maximum of two features per time step. On the left-hand side we have a feature assigning tracking approach. Imagine that we have two features (A and B) at a certain instant. If at the next one we only find one feature (C), which is closer to the old position of A than to the old position of B, then we will replace the MBP identifier C by A (always under the prerequisite that the features are close enough, determined by a maximum distance variable). The same happens in the following time step with feature E, which is most likely not a new structure but the continuation of A. On the right-hand side of Fig. 9 we have the pointer solution to the tracking problem. Here we are not replacing the feature identifier with a general track identifier and therefore assigning the measured MBP to a certain track. Rather, we add two pointers to every feature pointing to its predecessor and its follower. In this case one measurement can be used by several tracks (several other MBPs point to it) but overall the algorithm remains fairly simple by only allowing a given MBP to point to another one, namely, the spatially closest. For a more general discussion and the implications raised by other tracking approaches we refer to the publication of DeForest et al. (2007).
The feature assigning tracking has the advantage of creating a time reversible solution but the assignment of the feature track can be very complicated if there are several features close by in one or more instances giving rise to plentiful possible combinations and hence potential feature tracks. This is particularly true if the only criterion the decision is based on are the locations of the barycentres as in this work. More sophisticated tracking algorithms (but slower in runtime) would also incorporate the shape of the tracked features and look for example on the overlap between two realised instances of the feature (e.g., Hagenaar et al. 1999). On the other hand the pointer tracking approach is fast and easy to realize but not reversible as the example, depicted in Fig. 9, shows. The thin lines on the right-hand side subfigure show the realised MBP tracks. Clearly one of the solid line tracks (B-C-E) is different to the dashed line track (D-C-A). This means, if the time would be reversed, different features would get connected to each other. At the end, especially for our case, it is not such a serious problem as we are i) in quiet Sun and hence can expect that there are mostly isolated MBPs and ii) we are not going to estimate parameters which could be crucially influenced by this behaviour. Such parameters would be number densities, fractional occupation rates of the surface (here double counting of features would lead directly to errors), flux emergence or decay rates and of course lifetimes.

Inversions:
The used data are among the best currently available spectropolarimetric data sets with regard to their temporal and spatial resolution and in particular the consistently high data quality over the whole time series. Furthermore, the spectral resolution is sufficient to provide enough information for standard inversions such as Milne-Eddington (M-E) inversions. In our case the inversions were carried out with the SIR code (see Ruiz Cobo & del Toro Iniesta 1992) under the assumption of an easy atmospheric model stratification, namely assuming two nodes for the temperature (that is, linear perturbations), and one node (constant perturbations) for the three components of the magnetic field, the LOS velocity, and the microturbulence. Therefore, only seven free parameters are used, even less than with M-E inversions, but gaining additional information on the temperature stratification. Due to weak Stokes Q and U signals and the general tendency of inversion codes to overestimate magnetic field inclination in the case of such weak Stokes signals (see, e.g., Khomenko et al. 2003;Borrero & Kobel 2011, 2012 we refrained in this work from analysing and discussing the evolution of magnetic field inclinations. A closer look on the results of the inversion shows that the formal error of the inversions for the inclinations is in most cases about 25 degrees but in some very poorly defined cases (practically there is only noise in the Q and U component) this error can be well above 40 degrees. Some pixels on the other hand, with higher signal strengths, show formal inclination errors of less than 8 degrees.
The reliability of the SIR inversion results on these data for other plasma parameters was already assessed and reported by Requerey et al. (2014) who estimated typical uncertainties (rms) of 0.1 kG and 100 ms −1 from several different realizations of added noise to the data at a level of 10 −3 I c . In detail they repeated the inversions with 100 different realisations of noise and were thus able to assess the impact of noise on the quality of the inversions with the previously stated outcome. Hence we take these estimates also as valid for our analysis as we work with the same data set as well as with the same inversion code set up in the same fashion. Nevertheless, one could legitimately wonder whether there are limits above which the inferences are wrong, based on the typ-ical behavior expected for Stokes V in the two extreme cases of weak and strong fields. These doubts could be especially addressed to the strongest fields in our observations that might be too strong. The errors might come, e.g., from the quite large instrumental broadening on a fairly reduced sampling. We are, however, confident on the Requerey et al.'s estimates.
To digress a little on this, let us point out that, first, SIR already takes the broadening into account and, second, that arguments based just on Stokes V may fail because the inversion code makes full use of the information contained in all four Stokes parameters. Information from Stokes I is also important and even the fact that Stokes Q and U are typically small can be relevant to the final fit. Hence, our only intuition that is built on Stokes V may not be of good help. Last but not least, the most reassuring confirmation is already present in our statistical results. Note that field strengths in Fig.  8 barely exceed 1.5 kG. Among the approximately 2000 analyzed sets of profiles, only two of them display clearly erroneous values above 2 kG. We have neither eliminated nor analyzed in detail these two singular cases. On the one hand, showing them on the plot stresses the good quality of the results. On the other hand, a manual repair after a specific analysis is indeed irrelevant for our statistical purposes.

Comparison with previous studies:
Bellot Rubio et al. (2001) were among the first ones to report on the observational detection of a convective collapse (earlier observations were done already by, e.g., Sigwarth et al. 1999;Lin & Rimmele 1999). They studied exhaustively the observed Stokes profiles and performed also an inversion of the spectropolarimetric data. The occurrence of a convective collapse with magnetic field strengths increasing from 0.4 kG to 0.6 kG and downflow velocities of up to 6 km/s were confirmed by them. The overall process lasted for about 4 minutes before a seemingly rebounding plasma created a strong flow and shock which finally seemed to disintegrate the structure. One might ask at this point if it really takes a full-fledged instability to produce this moderate increase in magnetic field strength or if other mechanisms could also be responsible for moderate increases in the magnetic field? Shimizu et al. (2008) and coworkers were among the first to use the recent Hinode spacecraft and the onboard spectropolarimeter to investigate high-speed local mass downflows. They prepared a detailed study of the Stokes profiles taken in several interesting regions from the sunspot moat to the quiet Sun. In the case of the quiet Sun they were able to identify co-temporal brightenings in the G-band and in the Ca II-H filtergram occurring together with strong transient downflows. Nagata et al. (2008) was focusing then on such downflows in the quiet Sun and reported a detailed case study of a convective collapse and the formation of an MBP from the Hinode spectropolarimeter data. The inversion of the data yielded the magnetic field strength and the LOS velocity component. Furthermore, the brightness evolution was followed and, from the brightness difference, an estimate of the temperature difference was obtained. As in the study of Bellot Rubio et al. (2001) they found an increasing plasma downflow which went along with in-creasing magnetic field strength and brightness. The maximum field strength was reached about 3 minutes after the onset of the downflow with a value of 2 kG representing an amplification factor of four. The maximum downflow velocity was reached at the same time, with a value of 6 km/s. Afterwards, the velocity decreased while the brightness still increased for about 100 more seconds, with the magnetic field strength remaining nearly constant. The final point of the observed evolution was reached when the velocity reversed and an upflow of about 2 km/s set in, causing the magnetic field to reduce to about 1.5 kG along with a reduction of intensity. Unfortunately, the authors did not follow the evolution any further. We note that in this exceptionally well-studied evolution of an MBP the maximum field strength reached is near the upper limit of observed field strengths for MBPs. It did not take the community long to reproduce the results of Nagata et al. (2008) on a broader statistical base which was done by Fischer et al. (2009), who studied a total of 49 such convective collapse events with spectro-polarimetric data and filtergrams of Mg I-b2 as well as Ca II-H (higher photosphere and lower chromosphere). They were able to confirm the major implications of the convective collapse scenario; A shrinkage in size with observed downflows and the forming of bright points in the photosphere as well as in the chromosphere for most of the observed features.
The previous studies of Bellot Rubio et al. (2001) and Nagata et al. (2008) are congruent in reporting strong downflows (of about 6 km/s), followed later on by an upflow, which disintegrated in the first case the magnetic field, while it weakened the field in the second. Indeed we were able to identify such strong downflows in several cases, e.g. in 14 tracks the downflow exceeds 4 km/s (about 7% of the total amount of studied evolutions) and in five cases (about 3% of the total number of tracks) the value of 6 km/s is exceeded (or at least closely reached). Nevertheless this did not lead to an observable velocity reversal with upflows and a weakening of the magnetic field strength of MBPs. All the same, some upflows have been observed very close to our MBPs at given moments of their evolution. Our identification algorithm-with different criteria to that of Nagata et al. (2008)-or different amounts of magnetic flux contained in the structures could possibly explain the discrepancy. The identification in the blue line-wing could be ruled out as an explanation for the discrepancy as upflows would increase the probability of ending the tracking of a feature (causing lower brightnesses in the blue line-wing) and hence be in favour for detecting upflows at the end of the tracking period. Another possibility would be that the better spatial resolution reached in our study reveals spatially separated up and downflows where previous studies only found a mixture of both, resulting in an effective up or downflow depending on what quantity had a stronger influence on the average.
In a recent study, using the Imaging Spectro Polarimeter (CRISP; see Scharmer et al. 2008) at the Swedish Solar Telescope, Narayan (2011) presented results about strong transient downflows. These downflows were obviously associated with the creation of strong magnetic fields (flux tubes) and hence interpreted as manifestations of convective collapse. The author reports eight cases with downflows peaking at a maximum between 3.0 km/s and 5.2 km/s and showing peak magnetic field strengths from 0.48 kG to 1.32 kG. As in our study, the author does not find a reversal of the downflow into an upflow to cause the dispersion and weakening/disintegration of the flux concentration as reported by others. Furthermore the magnetic concentrations found by Narayan are of a transient nature with short lifetimes and do not lead to a permanent increase in magnetic field strength, in contrast to the MBP followed by Nagata et al. (2008). Overall our findings in the current study agree very well with the ones presented in the paper of Narayan. Furthermore, we are able to identify small concentrated strong upflows close to downflows in Narayan's study (as reported here for our second case) for his case g (Fig. A.6, first row fifth and sixth image within the marking circle), even though it looks in the particular case as if this would be the starting point of a newly born granular cell, and for his case h (Fig.  A.7, first row, second to fourth image left of the marked region).
In another recent study, Requerey et al. (2014) investigated in detail the evolution of a single magnetic structure from its formation due to flux concentration by granular advection, the further amplification by a process interpreted as convective collapse, till an oscillatory-like behavior at the final state of a collapsed flux tube is reached. During the convective collapse process they were able to observe the formation of two MBPs which were located close to the edge of the original magnetic structure which undergoes the collapse. Furthermore they showed that the plasma parameters of these features are correlated in a complex and not straightforward manner to each other (undergoing rapid fluctuations without a clear correlation to the brightness of the MBP) while the core of the magnetic structure (which undergoes the collapse) shows a nice correlation between magnetic field strength and LOS velocity. These findings may explain the complexity seen in the evolutionary tracks of the plasma parameters (Figs. 3, 5, and 7).
Finally, our study agrees with recent computer simulations such as the one of Danilovic et al. (2010) where the authors compare the magnetic field amplification process as seen in MHD simulations with those observed by Hinode/SOT data. They found downflow velocities between 5 km/s and 10 km/s in the simulations. For the strongest of the studied downflows they observed a velocity reversal after 200 s. However, the upflow did not disintegrate the magnetic field concentration, in contrast to the findings of Bellot Rubio et al. (2001), but led to a considerable weakening of the field as was observed by Nagata et al. (2008). An interesting recent study shedding light on the evolution of MBPs is the one of Hewitt et al. (2014), which shows nicely the creation of strong kG magnetic field brightenings in numerical simulations. The major difference to most observational papers is the non-cotemporality between the magnetic field strength, brightness, and LOS velocity maximum. Only our third observational case agrees with their general finding of preceding velocity maxima before field strength and intensity maxima are reached. Besides, they see moderate downflows with 3 to 4 km/s leading to really strong magnetic fields exceeding even the 2 kG threshold. In one of their three cases they are able to identify a very small flow reversal but generally speaking the evolution of the MBPs is dominated by downflows. Interestingly the brightness maxima of their first case lies between a downflow and a weaker upflow like in the second presented case of the study at hand.
If we compare the properties of our features, which are mainly located in the internetwork quiet Sun, with the results of high-resolution studies of magnetic elements in active region plage, we see clear differences. For example, our magnetic features display often sizable downflows over most of their generally short lifetime (at least as defined as their visibility as bright points), see Figs. 8,11,and 12. The magnetic features in plage regions are surrounded by strong downflows, but in their interiors they tend to be relatively free of downflows, as found by Rouppe van der Voort et al. (2005), Langangen et al. (2007) and Bühler (2013). Obviously, these plage magnetic elements have a rather different internal structure than the highly dynamic internetwork structures that we are studying.

Evolutionary MBP track statistics:
An important part of the current study considered the statistics of evolutionary MBP tracks and the results gained by such statistics.
Let us now discuss in more detail the statistics of the magnetic field strength evolution. As we found in Sect. 4.4, the distribution of initial magnetic field strengths is centred around the equipartition field strength (actually the mean value of all the initial fields is 0.46 ± 0.17 kG). MBPs are clearly distinct from the background magnetic field. This can be seen in Fig. 10 too. Here we plot the distribution of the initial magnetic field strengths (yellow triangles), the maximum magnetic field strength reached during the evolution (orange crosses), the dissolution field strength (black asterisks; average value of 0.41 ± 0.18 kG), and the background magnetic field strength of the complete FOV (blue diamonds). The complete magnetic field strength distribution, considering all sin-  gle measurements (i.e., every barycentre magnetic field strength of the identified MBPs in all frames), is represented by green circles. The correspondingly coloured lines are the results of best fits of log-normal distributions to the data. Such log-normal distributions were already found for other characteristic MBP parameters in the past (e.g., for the feature size; Abramenko et al. 2010;Crockett et al. 2010) and can be related to the fragmentation processes MBPs are undergoing (splitting and merging; see, e.g., Bogdan et al. 1988;Limpert et al. 2008). The distributions turned out to be fair represen-tations of the distributions except for the maximum field strengths that tend to show a bimodal distribution with values either below or above 1 kG. Such bi-modality was already reported in Utz et al. (2013b), who speculated that the two maxima might be related to internetwork fields and network fields.
This coincides with results of Beck et al. (2007) who investigated the magnetic field strength distribution (and other parameters) of MBPs around a sunspot moat and found a continuous distribution of field strengths starting with values around 0.5 kG up to 1.5 kG. On the other hand Lagg et al. (2010) showed (using the same data sets as in the current study) that at least some of the magnetic features in the internetwork have kG strength. Also, bear in mind the work of Riethmüller et al. (2014) who show that only magnetic features with kG fields are associated with a clear MBP, but that many of these appear as weaker fields in the observations, due to the insufficient spatial resolution. Furthermore straylight might affect the Stokes measurements and lead to a reduced detected magnetic field strength. A similar conclusion was also drawn by Criscuoli & Uitenbroek (2014), who used simulated data and investigated the behaviour of the field strength distribution with changing spatial resolutions and instrument effects.
The dashed lines (FOV distribution and distribution of all MBP measurements) were normalised to the same total frequency (number of measurements) as for the other distributions, i.e., to the number of analysed MBP tracks. This was done to allow a better comparison between the distributions. Else the FOV distribution and the distribution of all MBP measurements would be on a different absolute scale (far away from the others). The vertical dash-dotted lines in Fig. 10 indicate the average background noise magnetic field strength over the FOV (0.14 kG; see Sect. 4.4) and an upper boundary of the equipartition field (0.4 kG; see Solanki et al. 1996, who gives a range for the equipartition field from 0.2 to 0.4 kG). Note that the mean value indeed coincides with the expectation value of the (log-normal) distribution of magnetic field strengths over the full FOV. Such an expectation value is given for a log-normal distribution in general by where µ and σ are the mean and standard deviation of the natural logarithm of X. We see that the initial magnetic field strength distribution is separated by two times the rms value from the distribution of the complete FOV. We stated before that the expectation value (generally also known as mean value) for a log-normal distribution is situated at the right slope of the distribution due to its skewness. Very interestingly, the expectation value of the initial magnetic field strength distribution is 0.424 kG, which lies within the range of reasonable values for the equipartition field strength. This discussion supports the previously stated result that MBPs, which resemble the strongest part of magnetic structures (the core of magnetic structures) seem to be formed out of larger patches of equipartition magnetic fields and possess a distinct magnetic field strength distribution that is distinct from the background magnetic fields. This result is in agreement with the suggestion made by Venkatakrishnan (1986) that for the convective collapse to be successful the magnetic features must have a minimum size, i.e., a minimum amount of magnetic flux (as was observationally confirmed by Solanki et al. 1996). Indeed, 30 % (63 out 207) tracked MBPs reach kG field strengths at a given moment of their lifetime (although this fraction may be influenced by the limited spatial resolution).
As we have seen in Sect. 4.4 and in Fig. 8, initial downflows are on the order of 1 to 2 km/s, which increase twoto threefold at some point during their evolution, reaching values from 4 to 6 km/s. We show in Fig. 11 the LOS velocity distributions of the initial measurement (yellow triangles), the distribution of maximum LOS velocities (orange crosses) reached during the evolution and the dissolution distribution (black asterisks). Furthermore, as before, we depict the distribution of all LOS velocities in the FOV for the sake of comparison (blue diamonds) and of all single MBP measurements (green circles). The latter two distributions were rescaled to the number of tracks (corresponding to the measurements entering the previous three distributions) to enable an easier comparison. The lines are best-fit Gaussians. Although the different fits may not be equally good for all the distributions, their Gaussian mean parameters can easily be used to see that MBPs statistically start to be detected when downflows have already developed (0.8 km/s), that maximum velocities are clearly higher (2.4 km/s), and that dynamic effects fade away at the same time that the MBPs disappear (0.6 km/s). Interestingly, the green curve shows a maximum very close to that of the orange curve, i.e., below 1 km/s. The main difference is that the green curve is somewhat broader, with larger, but also a few smaller velocities. This suggests that for most of their time, the MBPs have velocities similar to that at the beginning.
To study in more detail the pre-and post-collapse phase we investigate the distribution of LOS velocities prior and posterior to the moment that the maximum field strength is reached, let us call this time t max . Figure 12 shows the histogram of the LOS velocities prior to t max as black solid line (downward-pointing solid triangles) and posterior to t max as yellow line (solid upward-pointing triangles). Moreover the green dashed curve gives the ascending phase of the MBP evolution and the descending phase of the evolution after t max is shown in blue colour. The latter four curves are slightly shifted relative to each other suggesting that the velocities prior to the maximum field strength are slightly higher compared to the posterior phase. Additionally the LOS distributions just prior and posterior to t max seem to be more skewed and hence also deviate more from Gaussian fits (illustrated by the corresponding lines) than the distributions for the full ascending/descending phase. The distribution at t max is plotted in violet colour and for comparison we display again the distribution of maximum LOS velocities during the features' lifetime. Interestingly the distribution at t max is clearly shifted to lower values compared to the maximum LOS velocity distribution (red colour) suggesting that the highest velocities during the MBP evolution are not reached at exactly t max . Nevertheless the velocities during the maximum are on average higher than at the other phases of the evolution (green/blue curve). The green and the blue curve have been again normalised to the number of tracks, as in Fig. 11. Special interest is often drawn to the question of flow reversals or upflows during the evolution of the MBP. An investigation of the found distribution yields values of 13 % of upflows prior to t max and 14% just posterior to this moment, which is therefore not significantly different. Interestingly this picture changes if we have a look at the distribution right at t max . Here we can only find in 10 % of the cases an upflow. Additionally the evolution seems to be not symmetric, which we can see if we consider all the measurements during the two phases of the MBP evolution (the phase prior to t max and posterior). Here we find only 11 % of the LOS measurements giving upflows in the ascending phase of the magnetic field evolution but 17 % of the LOS measurements give upflows in the descending phase of the magnetic field. Therefore it is probably safe to say that even if we cannot find direct flow reversals and strong upflows after the creation of the MBPs, the phase after t max has generally speaking weaker downflows and exhibits also more often some upflows.
The last detail we focus on is the relative timing of the maximum seen in different plasma quantities of the MBP. The lifetime of the MBPs has been normalised for that purpose. Negative values mean in that respect that the maximum occurred already before the tracking and a value greater of one that the maximum happened in the post-tracking phase. The resulting histograms of the maximum field strength (marked by red crosses) and the maximum LOS velocity (green asterisks are plotted in Fig. 13). Interestingly both quantities show their maximum just in the beginning of the evolution. The coincidence in the timing of the maxima of both quantities can be considered in agreement with former studies having less spatial and temporal resolution. It is not at odds with the previous finding of asynchrony between field strength and LOS velocity maxima in some specific cases since the histogram only gives statistical information. It does not say that for a particular feature the maximum velocity must occur at the same time as the maximum field strength.
In this study we tracked the evolution of small-scale magnetic fields that manifest themselves as MBPs in spectropolarimetric data observed with Sunrise/IMaX. The used data set was taken during a period of exceptionally quiet Sun close to the solar disc centre in 2009 June. The identification and tracking of the features were done automatically in the blue Fe line wing (−40 mÅ off the line centre) which appears by visual inspection to have similar properties to the G-band for the identification of MBPs although this proxy has to be used carefully as the intensity might be influenced by vertical motions along the LOS. Finally, plasma parameters of interest were extracted from fully processed maps obtained by the inversion code SIR.
The three case studies of MBP evolution presented in detail illustrate the variety of evolutionary histories and indicate that a variety of processes may be involved in the creation of MBPs as well as for stabilising them during their lifetime. In fact the observed plasma parameters are in a complicated way related to each other and drawing conclusions from one single observed track would not be generally applicable and may be misleading. One conclusion which can be drawn, however, is that MBPs, at least the ones in the quiet Sun, are rapidly evolving structures with constant, significant fluctuations, i.e., most of them are transient features and the longer lived ones show several brightness and field strength peaks. A similar conclusion was drawn by Jurčák et al. (2013), Narayan (2011), and by Fischer et al. (2009), who did a statistical analysis of convective collapse events in a Hinode data set.
We have in addition obtained more generally applicable results by considering all the observed tracks: • The magnetic field strength starts weak with values near the upper end of estimates of the equipartition field strength of about 0.4 kG, reaches maximum values two to three times the initial value, which is in most cases still below 1 kG probably also caused by resolution effects, and declines to values consistent with the equipartition field strength before the feature dissolves.
• The LOS velocity shows a downflow practically all the time. It starts with a weak value of about 0.8 km/s, reaches a clear maximum (in a few cases of up to 6 km/s but on average 2.4 km/s) before declining again to about 0.6 km/s on average, i.e., a value smaller than the starting value.
These observational facts are partly congruent with the expectations one would have from the convective collapse scenario and confirm on a statistical basis some of the expected implications of the theory such as an increase of the magnetic field strength coupled with strong, transient downflows. The presence of significant downflows co-located with the bright point throughout its lifetime of most of our features is in complete contrast with the absence of such downflows in bright points in active region plage. Moreover, we have outlined in the second studied case the importance of near small-scale upflows. The importance, evolution, and driving force of such spatially co-located up and downflows (shear flows) have to be studied in more detail in the future. Another point to be investigated in the future is the decay of the features. Are MBPs really disintegrated by a downflow reversal leading to an upwards drifting shock front? If so, then why are upflows recorded so rarely in our statistical sample? Only 33 out of a total of 207 investigated MBPs displayed upflows just before they disappeared as bright points. An even smaller fraction of MBPs display upflows just before or at the time their B values start to decrease. These values suggest that another mechanism is responsible for the weakening of the field than strong upflows. Naturally, we cannot exclude the possibility that upflows (maybe caused by upwards propagating shocks as discussed in Bellot Rubio et al. 2001) are too quickly evolving so that they cannot be detected within the temporal resolution of about 33 s at hand. Some of the previous questions have been partially addressed by Requerey et al. (2014), but it would be worthwhile to search for longer and more stable MBP evolutions (such as the one reported by Nagata et al. 2008) and to investigate what physical processes lead to the stabilisation. Due to their lateral expansion with height, the surface of flux tubes is in principle susceptible to the hydromagnetic exchange instability (e.g., Parker 1975). Therefore several possibilities were proposed in literature how to overcome the instabilities and stabilise flux tubes. This could be achieved by suitable environmental conditions such as the right temperature fields (Bünte et al. 1993b,a) or swirl flows (Schüssler 1984) around the magnetic flux tube. The latter is definitely worthwhile to be studied in more detail as recently the interest in swirls and vortices (e.g. Bonet et al. 2008;Shelyag et al. 2011) in the solar atmosphere has gained a lot of interest again due to new observations showing clearly such "tornado-like" structures (e.g. Wedemeyer-Böhm et al. 2012;Wedemeyer et al. 2013).
Another rather interesting follow-up study would be to check if the downflows (and hence the MBPs) studied here were preceded by magnetic flux emergence (e.g. as seen in linear polarization), i.e., could these downflows really be the signature of convective collapse, or rather are they simply the draining of rising loops that emerged recently? Possibly there is a component of both in the observed downflows associated with the enhancement of the magnetic field strength?
Deutsches Zentrum für Luft-und Raumfahrt e.V. (DLR), Grant No. 50OU 0401, and by the Innovationsfond of the President of the Max Planck Society (MPG). The High Altitude Observatory (HAO) contribution was partly funded through NASA grant NNX08AH38G. This work was partly supported by the BK21 plus program through the National Research Foundation (NRF) funded by the Ministry of Education of Korea. Finally the authors wish to express their sincere gratitude to the anonymous referee whose comments and remarks helped significantly to improve the manuscript. Moreover, the whole department is grateful for the English proof reader for creating a heated debate about the possibility for introducing new words into the English language. To put it into the words of the proof reader: "The words "till" and "still" have been changed to "suntil" here and elsewhere in the text. Please confirm if it is correct."