Issue 
A&A
Volume 569, September 2014



Article Number  A102  
Number of page(s)  6  
Section  The Sun  
DOI  https://doi.org/10.1051/00046361/201424221  
Published online  30 September 2014 
Observational evidence for buffetinginduced kink waves in solar magnetic elements
^{1} INAF – Osservatorio Astronomico di Roma, 00040 Monte Porzio Catone (RM), Italy
email: marco.stangalini@inaf.it
^{2} INAF – Istituto di Astrofisica e Planetologia Spaziali, 00133 Roma, Italy
^{3} Universita degli Studi di Roma Tor Vergata, 00133 Roma, Italy
^{4} Istituto Nazionale di Geofisica e Vulcanologia, 00143 Roma, Italy
Received: 16 May 2014
Accepted: 12 August 2014
The role of diffuse photospheric magnetic elements in the energy budget of the upper layers of the Sun’s atmosphere has been the recent subject of many studies. This was made possible by the availability of high temporal and spatial resolution observations of the solar photosphere, allowing large numbers of magnetic elements to be tracked to study their dynamics. In this work we exploit a long temporal series of seeingfree magnetograms of the solar photosphere to study the effect of the turbulent convection on the excitation of kink oscillations in magnetic elements. We make use of the empirical mode decomposition technique in order to study the transverse oscillations of several magnetic flux tubes. This technique permits analysis of nonstationary time series like those associated to the horizontal velocities of these flux tubes, which are continuously advected and dispersed by granular flows. Our primary findings reveal the excitation of low frequency modes of kink oscillations, which are subharmonics of a fundamental mode with a 7.6 ± 0.2 min periodicity. These results constitute observational proof of the excitation of kink waves by the buffeting of the convection cells in the solar photosphere, and they are discussed in light of their possible role in the energy budget of the upper Sun’s atmosphere.
Key words: Sun: photosphere / Sun: magnetic fields / Sun: oscillations
© ESO, 2014
1. Introduction
It is widely known that small magnetic elements in the solar atmosphere host a rich variety of magnetohydrodynamic (MHD) waves that may play an active role in the energy budget of the upper layers of the Sun’s atmosphere. This possibility appears even more attractive since high resolution observations have shown that small magnetic elements with diameters ~100−150 km cover a significant fraction of the solar photosphere (Lagg et al. 2010; Bonet et al. 2012), although it has been shown that magnetic flux emergence is not homogeneous (Stangalini 2014).
After being generated and amplified, possibly by the action of smallscale dynamo (SSD) acting in the solar convection zone (Lites 2011; Buehler et al. 2013), the smallscale magnetic fields are perturbed by the external photospheric plasma flows. Under this action the flux tubes are shaken and, at the same time, advected toward the boundaries of the supergranular cells (Abramenko et al. 2011; Orozco Suárez et al. 2012; Lepreti et al. 2012; Giannattasio et al. 2013, 2014) to form the socalled magnetic network (Simon & Weiss 1989). Orozco Suárez et al. (2012), employing a long time series of high resolution magnetograms (8 h), observed the small magnetic concentrations within a supergranule moving radially outward from the center toward the network, with a velocity aligned to the plasma flow. In this regard, Giannattasio et al. (2014) observed a change in the diffusion regime of small magnetic elements with a reduction of the diffusivity in the proximity of the magnetic network. This implies different dynamic regimes in supergranular cells.
At the same time many oscillations (kink, sausage, etc.) are expected in the flux tubes in response to the external plasma forcing (Nakariakov & Verwichte 2005). This scenario is supported by numerical simulations (Khomenko et al. 2008b; Fedun et al. 2011; Vigeesh et al. 2012; Nutto et al. 2012) and observations (Volkmer et al. 1995; Jess et al. 2009; Martínez González et al. 2011; Jess et al. 2012a,b; Jafarzadeh et al. 2013; Stangalini et al. 2013b), as well as by several theoretical studies (Roberts & Webb 1978; Spruit 1981; Edwin & Roberts 1983; Roberts 1983; Musielak et al. 1989; Steiner et al. 1998; Hasan et al. 2003; Musielak & Ulmschneider 2003a, to name a few).
Stangalini et al. (2013b) observed that kink waves coexist with longitudinal compressive waves in the same flux tube at the same time and are closely related, possibly through nonlinear interactions. Although solar convection is commonly believed to be the primary driver of waves in flux tubes, little is known about the excitation mechanisms of these waves. Hasan et al. (2003), Musielak & Ulmschneider (2003b,a), and Hasan & van Ballegooijen (2008) argue that the photospheric forcing of magnetic elements can be considered a viable mechanism through which energy is extracted from the photosphere and channelled to the upper layers of the atmosphere. It has been shown, in fact, that the rapid footpoint motion due to the turbulent granular buffeting can effectively excite kink waves that can propagate upward and couple with longitudinal waves (Kalkofen 1997; Hasan et al. 2003). By using high resolution observations and simulations of the solar photosphere, Keys et al. (2011) found that while the majority of smallscale magnetic elements have velocities between 0 and 1 km s^{1} and 6% of them have velocities in excess of 2 km s^{1}. The same authors argue that this significant fraction of magnetic elements may play an important role in the energy budget of the upper layers of the solar atmosphere, since velocities over 2 km s^{1} can drive kink waves effectively (Choudhuri et al. 1993). These velocities are consistent with the results found by many authors (see for example Möstl et al. 2006; Criscuoli et al. 2012; Chitta et al. 2012; Bodnárová et al. 2014; Stangalini et al. 2013b, to mention a few).
Jess et al. (2012b) detected slow upward propagating longitudinal waves in small magnetic elements, using high cadence broadband 2D data, visible as periodic intensity fluctuations in the range 110−600 s. This provided an observational proof of the propagation scenario.
Morton et al. (2013) have shown, by exploiting high spatial and temporal resolution observations, that torsional modes (Erdélyi & Fedun 2007) generated by vortices in the solar photosphere can effectively propagate to the chromosphere.
Very recently, Morton et al. (2014) have investigated the propagation of incompressive waves to the chromosphere, finding a good agreement between photospheric and chromospheric velocity power spectra at frequencies below 8 mHz.
Stangalini et al. (2013a) have observed the presence of high frequency peaks of up to 10−12 mHz in the power spectra associated to the horizontal velocity of individual magnetic elements. These high frequency oscillations are well above the cutoff frequency expected for kink waves, which is always smaller than the acoustic cutoff at ~5.3 mHz (Spruit 1981), therefore they are expected to propagate upwards. Although the physical properties of the magnetic elements analyzed were quite similar, Stangalini et al. (2013a) found that each individual magnetic element has its own signature in the power spectrum. This was interpreted by considering large scale physical parameters such as the connectivity of each magnetic element and the complexity of the magnetic carpet that make the power spectrum of the horizontal oscillations of a flux tube not only dictated by the local physical properties of the flux tube.
Stangalini et al. (2013a) estimated the response function of the magnetic elements subject to the forcing of the photospheric plasma flows by comparing the median power spectrum of kink oscillations in magnetic elements with the spectrum of the horizontal velocity associated to the granulation found by Matsumoto & Kitai (2010). This response revealed that magnetic elements are prone to oscillating at high frequency (ν> 5−6 mHz).
These results, as well as most of those present in the literature so far, were obtained by employing FFT^{1}based techniques. Since the process associated to the horizontal swaying of the flux tubes is far from being stationary, the use of FFTbased techniques may pose some limitation. Indeed, the modal basis used for the decomposition of the signals in FFT analysis is not well suited for analyzing nonstationary processes.
To overcome this shortcoming in this work we analyzed the horizontal fluctuations of several longlived photospheric flux tubes making use of empirical mode decomposition (EMD), a technique that is based on a modal basis that is defined on top of the signal itself (Huang et al. 1998). This allows us to highlight some peculiarities in flux tube oscillations as well as a method of providing direct observational proof of the nonlinear interaction between the external plasma and the magnetic elements, which excite MHD waves.
Fig. 1 Time averaged magnetogram obtained by averaging over the whole fourhours data set (saturated between −100 G and 100 G). The red dots indicate the mean position of the longest lived magnetic elements (≥4800 s) considered in the analysis. 

Open with DEXTER 
2. Data set and analysis
The data set used in this work consists of a sequence of high spatial resolution magnetograms acquired by SOT/NFI, the narrow band imager on board Hinode satellite (Tsuneta et al. 2008). The magnetograms were obtained in the Na I 589.6 nm spectral line from shutterless V and I Stokes filtergrams taken on 2008 August 18, close to disk center. The cadence of the data is 30 s. The diffraction limit in NFI filtergrams at 589.6 nm is 0.24 arcsec, while the pixel scale is set at 0.16 arcsec. This results in a slight spatial downsampling in the focal plane, providing a 2pixel Nyquist spatial sampling of 0.32 arcsec (or ~230 km on the solar photosphere).
The fieldofview (FoV) is approximately 80 × 80 arcsec, and the total duration of the time sequence is approximately 4 h. In Fig. 1 the average magnetogram is shown. The FoV encompasses a few supergranules whose borders are highlighted by the network patches visible in the figure.
Along with the standard SOT/FG calibration procedure (the IDL code prepfg available in the Hinode Solarsoft package), we also apply an additional calibration to limit the effects of jitter and tracking residuals. The visual inspection of the data in fact reveals a low frequency trend in the tracking and a few sudden shifts of the FoV which are not properly handled by the calibration code. Since we are mainly interested in the study of transverse oscillations of small magnetic elements, it is important to properly coalign the data set with subpixel accuracy. A registration procedure allowing subpixel accuracy must be applied. This procedure is based on FFT crosscorrelation and utilizes the whole FoV to estimate any misalignment between two images. We apply FFT registration iteratively until the mean residuals are minimized. This happens in three registration iterations.
To study the dynamics of small magnetic features we apply the YAFTA (Yet another Feature Tracking Algorithm) tracking algorithm (Welsch & Longcope 2003; DeForest et al. 2007). This algorithm identifies and tracks magnetic pixels belonging to the same local maximum. To ensure the reliability of the results, three constraints are applied. Each magnetic feature must lie above a threshold and must have an area slightly larger than the fullwidth at halfmaximum (FWHM) of the pointspread function (a 4 × 4 pixel box in our case) and a lifetime long enough to ensure a high spectral resolution. The threshold on the magnetic signal is chosen to be 2σ. Following Liu et al. (2012) we estimate the sigma of the magnetic signal by fitting a Gaussian to the lowfield pixels with absolute value of the magnetic flux density below 200 G. This resulted in a sigma of 11.8 G. Lastly, only the magnetic elements with a lifetime greater or equal to 4800 s are considered in the analysis (22 elements in total). This sets the frequency resolution at ~0.2 mHz. The threshold on the lifetime (4800 s) was chosen as a tradeoff between the number of elements and the frequency resolution. Thornton & Parnell (2011) and Stangalini et al. (2013b) have extensively tested this tracking code on high resolution data.
In Fig. 1 we show the average position of each magnetic element selected. From the position of each magnetic element given by the tracking algorithm the horizontal velocities are estimated. The position of the magnetic elements is estimated with subpixel accuracy as a fluxweighted average of the coordinates of the pixels constituting the magnetic element itself. The computation of the horizontal velocity using both the x and y components suffers from frequency doubling if the components themselves are oscillating around zero. This is a direct consequence of the definition of the total velocity where the squared value of the components is considered. In order to avoid this drawback in the spectral analysis we consider only the xcomponent of the horizontal velocity, hereafter v_{h} for simplicity.
2.1. Empirical mode decomposition
The EMD was originally introduced by Huang et al. (1998) as a precondition technique for the application of the Hilbert transform to a signal. The EMD moves from the assumption that most of the actual time series of natural systems are the result of the superposition of many different time scales, simultaneously active, and aims to decompose a signal into a set of intrinsic mode functions (IMFs), representing different oscillations at a local level. As opposed to the Fourier method, which decomposes a signal into a prefixed set of elementary functions (sine or cosine oscillations) on the basis of the energy content of each characteristic time scale and applies these to rigorously stationary time series, the characteristic time scales of the IMFs obtained by the EMD are based on the local distance between two successive extrema. In other words, an IMF represents a hidden oscillation mode, locally defined and thus not stationary. Indeed, an IMF is defined as a local monocomponent mode, which satisfies the following conditions: i) the number of extrema and zerocrossing should be equal or differ at the most by one; ii) at any time the mean value of the upper and lower envelopes, defined from the local maxima and minima respectively, is zero. It is important to point out that EMD technique decomposes a signal without any a priori assumption, providing a decomposition which is based on the data itself, i.e., a posteriori decomposition.
The EMD algorithm was employed to extract the IMFs of a time series x(t) and consists of a sequence of steps that can be summarized as follows:

1.
extract the local maxima and minima of the time series x(t);

2.
construct the upper (e_{u}(t)) and lower (e_{l}(t)) envelopes from local maxima and minima;

3.
calculate the mean m(t) = (e_{u}(t) + e_{l}(t))/2;

4.
define the new time series h_{i}(t) = x(t) − m(t);

5.
check if h_{i}(t) satisfies the condition required to be an IMF, if not repeat the steps 1−5 on the time series h_{i}(t);

6.
if h_{i}(t) satisfies the condition to be an IMF, define a new signal x^{i}(t) = x(t) − h_{i}(t) and start the procedure again to iteratively define the other IMFs;

7.
the iteration is stopped when the x^{k}(t) has less than two extrema. In this case x^{k}(t) is identified as the residue r(t).
The steps from (1) to (5) constitute the socalled sifting process, which is stopped on the application of one of the stoppage criteria, proposed in the wide range of literature about this technique (Huang et al. 1998, 2003, 2008; Flandrin et al. 2004). The application of a stoppage criterion is required to validate that the IMFs have a physical meaning. Here, we apply the stoppage criterion used by Rilling et al. (2003), fixing a maximal number of iterations of 300 steps and checking that the IMFs satisfy condition i). As a result of the EMD technique the original signal x(t) is decomposed into a set of IMFs and a residue, so that we can write: (1)where { h_{i}(t) } is the set of IMFs and r(t) is the residue. It is important to underline that the EMD method, being based on the extraction of the energy contained in the intrinsic time scales of the signal, is used to detect the rapid frequency changes that are typical features of nonstationary signals. Furthermore, as underlined in Dätig & Schlurmann (2004), compared to the usual Fourier transform, since the IMFs are not given in a closed analytical form, the mathematical orthogonality of the IMFs is only locally satisfied, instead of the usual theoretical sense (⟨ h_{i}(t)  h_{j}(t) ⟩ = δ_{ij}).
The EMD technique has already been used in many contexts of solar physics, heliospheric and magnetospheric physics such as, but not limited to, the study of coronal oscillations (Terradas et al. 2004), neutrino flux modulation (Vecchio et al. 2010, 2012), cosmic rays modulation (Laurenza et al. 2014), and magnetospheric studies (De Michelis et al. 2012).
3. Results
As described in the previous section, the EMD analysis provides a decomposition of the signals associated with the horizontal velocity of the selected magnetic elements in different intrinsic modes which are built on top of the time series itself. In our case, each v_{h} time series is decomposed into six IMFs (see for example Fig. 2). In the upper panel of Fig. 2 it is possible to see that the horizontal velocity of the flux tube shows strong fluctuations up to ± 4 km s^{1}. In addition, these strong fluctuations appear to be concentrated in short pulses in the time series, and the time series itself appears to be nonstationary. For each IMF we estimated the main frequency through an analysis of the autocorrelation function (ACF, see Fig. 3). This resulted in 132 frequency values in total (six values for each of the 22 magnetic elements). The ACF of each mode is fitted using an exponential function of the form (2)where A is a constant term, x represents the time lag, T_{k} is the period which has to be estimated from the fit, and τ is a characteristic time scale (the characteristic length of autocorrelation of the signal). It is worth noting here that each IMF is not purely monochromatic. This is because unlike the FFT, where each mode decomposing the time series has a unique frequency associated to it, in EMD, each mode contains different frequencies, although very close to each other. For this reason by fitting the ACF one obtains the mean period of oscillation. This is also why in Fig. 3 the fit of the ACF is not perfectly matching the ACF itself, but there are small differences.
Fig. 2 Empirical mode decomposition of the horizontal velocity of a single magnetic element into six intrinsic mode functions (h_{1} − h_{6}), and a residue (r). 

Open with DEXTER 
Fig. 3 Example of ACF fitting to estimate the main frequency of each IMF. 

Open with DEXTER 
Fig. 4 Probability density function of the period (upper panel). Position of the peaks of the PDF as a function of the peak number (lower panel). 

Open with DEXTER 
In Fig. 4 (upper panel), we plot the probability density function (PDF) of all the mean periodicities found in each IMF of the decomposition and for each magnetic element. Since the total number of frequencies is limited (one value for each of the six IMFs and for each of the 22 magnetic elements), the PDF is obtained using the kernel method, as explained in Kaiser & Schreiber (2002). The kernel method is in fact an alternative to binning a distribution and is best suited for cases in which the number of samples is limited. This approach is mainly based on the assumption that the probability density is smooth enough that one can ignore all the structures below a certain kernel width. Instead of simply counting points within a certain bin, one can consider distancedependent weights obtained by using a kernel function. In our case we made use of a Gaussian kernel function.
The PDF increases progressively at short periods and presents a series of distinct bumps corresponding to the most common periodicities found in the signals. In order to check whether these peaks have a distinct spacing, an exponential function plus a series of gaussian peaks is fitted to the data. This allows us to estimate the position of each bump. The results of this analysis are shown in Fig. 4 (lower panel). In this figure the period of each peak is plotted against the order k at which it is found. Here we start from the order k = 2 by assuming that the order k = 1 is hidden below the frequency resolution limit imposed by the length of the time series analyzed. This limit is represented by the dashed vertical line in the upper panel of Fig. 4. It is easy to see that the spacing of the most common periodicities is not random. The resulting scaling law is linear T(n) = T_{0}n where n represents the peak number and T the period.
A linear fit to the data results in a slope T_{0} corresponding to a period of 7.6 ± 0.2 min, where the error is obtained by a weighted linear regression. Here, the weights used are proportional to the reciprocal of the errors associated to the peak position as obtained by the multiGaussian nonlinear fit (LevenbergMarquardt method) of the PDF. This number is consistent with the time scale of granulation (between 5−10 min). In addition, the scaling relation of the peaks is the signature of subharmonics where T_{0} represents the fundamental period of the system.
Moreover, the EMD analysis shown in Fig. 2 also highlights the strong nonstationary nature of the signals under investigation which manifests rather clearly in the IMFs. All of them in fact display a frequency modulation and a rather common presence of wave trains with different characteristics (i.e. amplitude, modulation, and waiting time).
In the lower panel of Fig. 4, we also show the persistence of each nonstationary oscillation compared to the period of the oscillation. This piece of information is encoded in the different colors and the different widths of the points of the plot itself. This persistence is given by the characteristic decorrelation scale, τ, estimated from the fit of the ACF mentioned above. Thus the persistence represents the temporal autocorrelation of the signal, that is, the duration of a wave train. What is clear is that the shorter the period, the longer the persistence of a particular oscillation.
4. Discussion
Our analysis shows that at least the longestlived magnetic elements in our data set (22 elements in total) present several periodicities which are consistent with subharmonics. The fundamental mode is estimated to be ~7.6 ± 0.2 min, a period which is close to the evolutionary time scale expected for granulation in the solar photosphere. This means that the buffeting action that the photospheric plasma exerts on the magnetic elements induces oscillations with a period which is consistent with the granular lifetime (e.g. Title et al. 1989; Hirzberger et al. 1999; Del Moro 2004).
This is evidence for buffetinginduced transverse oscillations in photospheric magnetic elements. These oscillations are strictly nonstationary and for this reason they are investigated employing a decomposition whose modal base can be obtained directly from the signals themselves. This allows us to overcome the wellknown limitations of FFTbased techniques, when applied to nonstationary signals.
In addition, our results show the presence of a linear scaling of the periods of oscillations found the in the magnetic elements which is consistent with subharmonic oscillations.
It is worthwhile to briefly discuss here the mechanisms of excitation of subharmonics in a physical system. For a complete review of the topic we refer the reader to Han & Benaroya (2002).
Both subharmonics and superharmonics can be either excited in linear or nonlinear systems, although with some difference in the underlying mechanisms.
In linear systems, subharmonic or superharmonic response can only be expected when exact relations between the natural and the forcing frequency are satisfied. That is, the forcing frequency must be an exact fraction or multiple of the natural frequency. On the other hand, in nonlinear systems, the superharmonic or subharmonic response can be excited by a range of forcing frequencies with values not necessarily multiples or fractions of the natural frequency. It has been shown, in fact, that in nonlinear systems a response at a frequency other than the driving frequency is also observed (Linsay 1981). In particular subharmonics and superharmonics in nonlinear systems can be generated as a response to intense driving (Yen 1971), although they have a different character. While superharmonics are always generated in the presence of a nonlinear driving source, subharmonics need a certain forcing strength in order to be excited (see for example Stern 1965, for a complete treatment of the topic). Moreover, the presence of perioddoubling cascades, like that observed in our results, can be seen as the signature of chaotic excitations in nonlinear systems (Sander & Yorke 2009, 2010, to mention a few).
Matsumoto & Kitai (2010), using a local spatial crosscorrelation technique, studied the horizontal oscillations of the velocity of the photospheric plasma, which results in a broad power spectrum. For this reason, it is difficult to distinguish between a linear and nonlinear response of the flux tubes since in the linear regime the exact relation between the natural and the forcing frequency can be easily satisfied in the presence of a broad spectrum.
However, our results constitute a strong observational proof of the excitation of kink waves as a result of granular buffeting.
The flux tubes in the solar photosphere are embedded in a highβ plasma whose forcing may induce both forced and free oscillations in the flux tube itself. Stangalini et al. (2013a) found, studying a large number of photospheric magnetic concentrations, a small peak in the range 3−4 mHz of the kink oscillation spectrum, although the large error bars associated did not allowed them to identify that peak with a high confidence level.
Our results extend the results of Stangalini et al. (2013a) by analysing the power spectral density of kink oscillations in photospheric magnetic elements with the very large frequency resolution provided by the long duration of the time series associated to the elements selected (80 min). This allowed us to explore also the very low frequency part of the power spectrum (ν< 4 mHz down to 0.2 mHz) of kink oscillations.
Apart from this, our results also indicate that most of the energy is contained in higher frequency modes (P< 25 min), since the persistence of these modes is larger than the modes at lower frequency.
It is a wellknown result that the cutoff frequency of kink waves is always smaller than its acoustic counterpart (Spruit 1981). This is given by (3)where ω_{ac} represents the acoustic cutoff which is ≃5.3 mHz, and β is the gas to magnetic pressure ratio, and γ the adiabatic index. Under typical conditions the kink cutoff frequency in a flux tube is at least two times smaller than what is estimated. Of course it also depends on the inclination angle through ω_{ac} (McIntosh & Jefferies 2006). For this reason it is difficult to estimate its real value since it is not possible from our data to estimate the inclination angle of each magnetic element. However, if we assume, hypothetically, that the magnetic elements are nearly vertical, then frequencies above ~ 2.5 − 2.7 mHz will be able to propagate upwards. This is an upper limit for the cutoff frequency itself since if the magnetic element is inclined to the vertical, the cutoff will be even smaller than what is estimated, allowing the propagation of oscillations at even smaller frequencies. Regardless, the fundamental mode at ~7.6 ± 0.2 min (2.2 mHz) is outside the propagation regime and not all the magnetic elements may have a sufficient inclination to lower the cutoff frequency enough to allow this mode to propagate. However, in small magnetic elements the cutoff frequency can also be reduced by radiative losses (Khomenko et al. 2008a).
For these reasons, we believe that the subharmonic response of the flux tubes to granular buffeting does not largely contribute to cromospheric heating. However the subharmonic response observed provides useful information on the mechanisms of excitation of the flux tubes and represents a strong observational proof of the interaction between the photospheric plasma and the magnetic elements themselves.
Besides this, it is worth noting that the persistence of the wave trains is larger at higher frequencies. This means that the amount of energy eventually contained in each wave train is larger and larger as the frequency increases. This is obviously an important point that has several implications for the heating of the upper layers of the Sun’s atmosphere. Although our analysis is limited in frequency by the width of the kernel used to estimate the PDF, we can speculate that the measured increase of the persistence, seen up to ν ~ 1.6 mHz (T ~ 10 min), can in fact extend to even higher frequencies that are allowed to propagate upward. In other words, even if the lowfrequency kink oscillations were propagated they would have contributed little to the energy budget of the upper layers as their persistence is very short compared to their high frequency counterparts.
5. Conclusions
In this work we have analyzed kink oscillations of several photospheric magnetic elements by employing the EMD technique. This technique is ideal for the study of nonstationary processes. This method allows us to decompose each velocity signal associated to the magnetic element in different modes using a decomposition basis which is derived empirically from the signal itself. By studying the ACF of each mode we have estimated the main periodicities and their relative PDF.
Our findings demonstrated the presence of numerous oscillation peaks which are consistent with subharmonic oscillations of a fundamental whose period is ~7.6 ± 0.2 min. Such a period is close to the characteristic temporal scale of the photospheric convection cells, therefore we can argue that these oscillations are associated to buffeting induced oscillations.
Acknowledgments
This work has been partly supported by the PRININAF 2010 and PRINMIUR 2012 (prot. 2012P2HRCR) entitled “Il sole attivo ed i suoi effetti sul clima dello spazio e della terra” grants, funded by the Italian National Institute for Astrophysics (INAF) and Ministry of Education, Universities and Research (MIUR), respectively. 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 (UK), NASA, ESA, and NSC (Norway). We also thank the anonymous referee for the valuable comments that helped a lot to improve the manuscript.
References
 Abramenko, V. I., Carbone, V., Yurchyshyn, V., et al. 2011, ApJ, 743, 133 [NASA ADS] [CrossRef] [Google Scholar]
 Bodnárová, M., Utz, D., & Rybák, J. 2014, Sol. Phys., 289, 1543 [NASA ADS] [CrossRef] [Google Scholar]
 Bonet, J. A., Cabello, I., & Sánchez Almeida, J. 2012, A&A, 539, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Buehler, D., Lagg, A., & Solanki, S. K. 2013, A&A, 555, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chitta, L. P., van Ballegooijen, A. A., Rouppe van der Voort, L., DeLuca, E. E., & Kariyappa, R. 2012, ApJ, 752, 48 [NASA ADS] [CrossRef] [Google Scholar]
 Choudhuri, A. R., Dikpati, M., & Banerjee, D. 1993, ApJ, 413, 811 [NASA ADS] [CrossRef] [Google Scholar]
 Criscuoli, S., Del Moro, D., Giorgi, F., et al. 2012, Mem. Soc. Astron. It. Suppl., 19, 93 [Google Scholar]
 Dätig, M., & Schlurmann, T. 2004, Ocean Engineering, 31, 1783 [CrossRef] [Google Scholar]
 De Michelis, P., Consolini, G., & Tozzi, R. 2012, Nonlinear Processes in Geophysics, 19, 667 [NASA ADS] [CrossRef] [Google Scholar]
 DeForest, C. E., Hagenaar, H. J., Lamb, D. A., Parnell, C. E., & Welsch, B. T. 2007, ApJ, 666, 576 [NASA ADS] [CrossRef] [Google Scholar]
 Del Moro, D. 2004, A&A, 428, 1007 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Edwin, P., & Roberts, B. 1983, Sol. Phys., 88, 179 [NASA ADS] [CrossRef] [Google Scholar]
 Erdélyi, R., & Fedun, V. 2007, Science, 318, 1572 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Fedun, V., Shelyag, S., & Erdélyi, R. 2011, ApJ, 727, 17 [NASA ADS] [CrossRef] [Google Scholar]
 Flandrin, P., Rilling, G., & Goncalves, P. 2004, Signal Processing Letters, IEEE, 11, 112 [Google Scholar]
 Giannattasio, F., Del Moro, D., Berrilli, F., et al. 2013, ApJ, 770, L36 [NASA ADS] [CrossRef] [Google Scholar]
 Giannattasio, F., Stangalini, M., Berrilli, F., Del Moro, D., & BellotRubio, L. 2014, ApJ, 788, 137 [NASA ADS] [CrossRef] [Google Scholar]
 Han, S., & Benaroya, H. 2002, Nonlinear and Stochastic Dynamics of Compliant Offshore Structures, Solid Mechanics and Its Applications (Springer) [Google Scholar]
 Hasan, S. S., & van Ballegooijen, A. A. 2008, ApJ, 680, 1542 [NASA ADS] [CrossRef] [Google Scholar]
 Hasan, S. S., Kalkofen, W., van Ballegooijen, A. A., & Ulmschneider, P. 2003, ApJ, 585, 1138 [NASA ADS] [CrossRef] [Google Scholar]
 Hirzberger, J., Bonet, J. A., Vázquez, M., & Hanslmeier, A. 1999, ApJ, 515, 441 [NASA ADS] [CrossRef] [Google Scholar]
 Huang, N. E., Shen, Z., Long, S. R., et al. 1998, Proc. Roy. Soc. Lond. Ser. A: Mathematical, Physical and Engineering Sciences, 454, 903 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Huang, N. E., Wu, M.L. C., Long, S. R., et al. 2003, Proc. Roy. Soc. Lond. Ser. A: Mathematical, Physical and Engineering Sciences, 459, 2317 [CrossRef] [Google Scholar]
 Huang, Y., Schmitt, F. G., Lu, Z., & Liu, Y. 2008, EPL, 84, 40010 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Jafarzadeh, S., Solanki, S. K., Feller, A., et al. 2013, A&A, 549, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Jess, D. B., Mathioudakis, M., Erdélyi, R., et al. 2009, Science, 323, 1582 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Jess, D. B., Pascoe, D. J., Christian, D. J., et al. 2012a, ApJ, 744, L5 [NASA ADS] [CrossRef] [Google Scholar]
 Jess, D. B., Shelyag, S., Mathioudakis, M., et al. 2012b, ApJ, 746, 183 [NASA ADS] [CrossRef] [Google Scholar]
 Kaiser, A., & Schreiber, T. 2002, Physica D: Nonlinear Phenomena, 166, 43 [Google Scholar]
 Kalkofen, W. 1997, ApJ, 486, L145 [NASA ADS] [CrossRef] [Google Scholar]
 Keys, P. H., Mathioudakis, M., Jess, D. B., et al. 2011, ApJ, 740, L40 [NASA ADS] [CrossRef] [Google Scholar]
 Khomenko, E., Centeno, R., Collados, M., & Trujillo Bueno, J. 2008a, ApJ, 676, L85 [NASA ADS] [CrossRef] [Google Scholar]
 Khomenko, E., Collados, M., & Felipe, T. 2008b, Sol. Phys., 251, 589 [NASA ADS] [CrossRef] [Google Scholar]
 Lagg, A., Solanki, S. K., Riethmüller, T. L., et al. 2010, ApJ, 723, L164 [NASA ADS] [CrossRef] [Google Scholar]
 Laurenza, M., Vecchio, A., Storini, M., & Carbone, V. 2014, ApJ, 781, 71 [NASA ADS] [CrossRef] [Google Scholar]
 Lepreti, F., Carbone, V., Abramenko, V. I., et al. 2012, ApJ, 759, L17 [NASA ADS] [CrossRef] [Google Scholar]
 Linsay, P. S. 1981, Phys. Rev. Lett., 47, 1349 [NASA ADS] [CrossRef] [Google Scholar]
 Lites, B. W. 2011, ApJ, 737, 52 [NASA ADS] [CrossRef] [Google Scholar]
 Liu, Y., Hoeksema, J. T., Scherrer, P. H., et al. 2012, Sol. Phys., 279, 295 [NASA ADS] [CrossRef] [Google Scholar]
 Martínez González, M. J., Asensio Ramos, A., Manso Sainz, R., et al. 2011, ApJ, 730, L37 [NASA ADS] [CrossRef] [Google Scholar]
 Matsumoto, T., & Kitai, R. 2010, ApJ, 716, L19 [NASA ADS] [CrossRef] [Google Scholar]
 McIntosh, S. W., & Jefferies, S. M. 2006, ApJ, 647, L77 [NASA ADS] [CrossRef] [Google Scholar]
 Morton, R. J., Verth, G., Fedun, V., Shelyag, S., & Erdélyi, R. 2013, ApJ, 768, 17 [NASA ADS] [CrossRef] [Google Scholar]
 Morton, R. J., Verth, G., Hillier, A., & Erdélyi, R. 2014, ApJ, 784, 29 [NASA ADS] [CrossRef] [Google Scholar]
 Möstl, C., Hanslmeier, A., Sobotka, M., Puschmann, K., & Muthsam, H. J. 2006, Sol. Phys., 237, 13 [NASA ADS] [CrossRef] [Google Scholar]
 Musielak, Z. E., & Ulmschneider, P. 2003a, A&A, 400, 1057 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Musielak, Z. E., & Ulmschneider, P. 2003b, A&A, 406, 725 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Musielak, Z. E., Rosner, R., & Ulmschneider, P. 1989, ApJ, 337, 470 [NASA ADS] [CrossRef] [Google Scholar]
 Nakariakov, V. M., & Verwichte, E. 2005, Liv. Rev. Sol. Phys., 2 [Google Scholar]
 Nutto, C., Steiner, O., Schaffenberger, W., & Roth, M. 2012, A&A, 538, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Orozco Suárez, D., Katsukawa, Y., & Bellot Rubio, L. R. 2012, ApJ, 758, L38 [NASA ADS] [CrossRef] [Google Scholar]
 Rilling, G., Flandrin, P., Goncalves, P., et al. 2003, in IEEEEURASIP Workshop on Nonlinear Signal and Image Processing NSIP, 3, 8 [Google Scholar]
 Roberts, B. 1983, Sol. Phys., 87, 77 [NASA ADS] [CrossRef] [Google Scholar]
 Roberts, B., & Webb, A. R. 1978, Sol. Phys., 56, 5 [NASA ADS] [CrossRef] [Google Scholar]
 Sander, E., & Yorke, J. A. 2009, unpublished [arXiv:0910.3570] [Google Scholar]
 Sander, E., & Yorke, J. A. 2010, unpublished [arXiv:1002.3363] [Google Scholar]
 Simon, G. W., & Weiss, N. O. 1989, ApJ, 345, 1060 [NASA ADS] [CrossRef] [Google Scholar]
 Spruit, H. C. 1981, A&A, 98, 155 [NASA ADS] [Google Scholar]
 Stangalini, M. 2014, A&A, 561, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Stangalini, M., Berrilli, F., & Consolini, G. 2013a, A&A, 559, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Stangalini, M., Solanki, S. K., Cameron, R., & Martínez Pillet, V. 2013b, A&A, 554, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Steiner, O., GrossmannDoerth, U., Knoelker, M., & Schüssler, M. 1998, ApJ, 495, 468 [NASA ADS] [CrossRef] [Google Scholar]
 Stern, T. E. 1965, Theory of nonlinear networks and systems: An introduction (AddisonWesley) [Google Scholar]
 Terradas, J., Oliver, R., & Ballester, J. L. 2004, ApJ, 614, 435 [NASA ADS] [CrossRef] [Google Scholar]
 Thornton, L. M., & Parnell, C. E. 2011, Sol. Phys., 269, 13 [NASA ADS] [CrossRef] [Google Scholar]
 Title, A. M., Tarbell, T. D., Topka, K. P., et al. 1989, ApJ, 336, 475 [NASA ADS] [CrossRef] [Google Scholar]
 Tsuneta, S., Ichimoto, K., Katsukawa, Y., et al. 2008, Sol. Phys., 249, 167 [NASA ADS] [CrossRef] [Google Scholar]
 Vecchio, A., Laurenza, M., Carbone, V., & Storini, M. 2010, ApJ, 709, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Vecchio, A., Laurenza, M., Meduri, D., Carbone, V., & Storini, M. 2012, ApJ, 749, 27 [NASA ADS] [CrossRef] [Google Scholar]
 Vigeesh, G., Fedun, V., Hasan, S. S., & Erdélyi, R. 2012, ApJ, 755, 18 [NASA ADS] [CrossRef] [Google Scholar]
 Volkmer, R., Kneer, F., & Bendlin, C. 1995, A&A, 304, L1 [NASA ADS] [Google Scholar]
 Welsch, B. T., & Longcope, D. W. 2003, ApJ, 588, 620 [NASA ADS] [CrossRef] [Google Scholar]
 Yen, N.C. 1971, Subharmonic generation in acoustic systems, Tech. Rep., DTIC Document [Google Scholar]
All Figures
Fig. 1 Time averaged magnetogram obtained by averaging over the whole fourhours data set (saturated between −100 G and 100 G). The red dots indicate the mean position of the longest lived magnetic elements (≥4800 s) considered in the analysis. 

Open with DEXTER  
In the text 
Fig. 2 Empirical mode decomposition of the horizontal velocity of a single magnetic element into six intrinsic mode functions (h_{1} − h_{6}), and a residue (r). 

Open with DEXTER  
In the text 
Fig. 3 Example of ACF fitting to estimate the main frequency of each IMF. 

Open with DEXTER  
In the text 
Fig. 4 Probability density function of the period (upper panel). Position of the peaks of the PDF as a function of the peak number (lower panel). 

Open with DEXTER  
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.