# High-frequency Oscillations in Small Magnetic Elements Observed with Sunrise/SuFI

## Abstract

We characterize waves in small magnetic elements and investigate their propagation in the lower solar atmosphere from observations at high spatial and temporal resolution. We use the wavelet transform to analyze oscillations of both horizontal displacement and intensity in magnetic bright points found in the nm and the Ca ii H nm passbands of the filter imager on board the Sunrise balloon-borne solar observatory. Phase differences between the oscillations at the two atmospheric layers corresponding to the two passbands reveal upward propagating waves at high frequencies (up to mHz). Weak signatures of standing as well as downward propagating waves are also obtained. Both compressible and incompressible (kink) waves are found in the small-scale magnetic features. The two types of waves have different, though overlapping, period distributions. Two independent estimates give a height difference of approximately km between the two atmospheric layers sampled by the employed spectral bands. This value, together with the determined short travel times of the transverse and longitudinal waves provide us with phase speeds of km s and km s, respectively. We speculate that these phase speeds may not reflect the true propagation speeds of the waves. Thus, effects such as the refraction of fast longitudinal waves may contribute to an overestimate of the phase speed.

## 1Introduction

Wave phenomena in solar magnetic features have been proposed as a prime means of transferring energy from the solar interior and lower atmosphere to the place where it is needed to maintain the high temperatures of the mid and upper solar atmosphere as well as to power the solar wind [?]. Due to the fact that acoustic waves shock and dissipate in the chromosphere, much recent attention has focused on waves associated with the magnetic field. Most observations of propagating waves in magnetic structures are either restricted to the chromosphere and higher layers, or to large features such as sunspots. Here, we present observations of short-period waves propagating along thin, magnetic flux tubes in the lower solar atmosphere, whose transported energy is a promising source of the intensity enhancements of magnetic bright points (MBPs) seen in strong chromospheric lines such as Ca ii H [?].

Different modes of oscillation are present in the magnetized solar atmosphere, depending on their compressive/non-compressive and magnetic nature [?]. Such oscillations occur in propagating or standing states [?]. The magnetic fields act as a guide for magnetohydrodynamics (MHD) waves effectively increasing the connection between different layers of the solar atmosphere. It has been shown that longitudinal acoustic waves in the photosphere propagate along the magnetic field lines and their leakage to the upper atmospheric layers depends on the inclination of the magnetic field [?].

Magnetic waves, excited in a thin flux tube, are generally classified, for simplicity, according to their propagation speeds: (1) Alfvén waves propagate with the local Alfvén speed, e.g., inside a flux tube, (2) axisymmetric, longitudinal magneto-acoustic waves propagate with a tube speed that is smaller than both sound and Alfvén speeds (), and (3) non-axisymmetric, transverse kink waves propagating with a speed between that of the surroundings and [?].

Within thicker flux tubes, additional wave modes also propagate: the fast and the slow mode, traveling at speeds faster or slower than the speed of sound, respectively. Thus the speed of fast and slow waves depends on the ratio of and on the direction of propagation relative to the field lines. In general, both modes are compressible. The slow mode is closely related to the tube mode (longitudinal magneto-acoustic mode; [?]). While slow-mode waves dissipate in the chromosphere by forming shocks, fast waves can penetrate into the upper solar atmosphere (due to their higher phase speeds, which reduce the occurrence of shocks). Hence, the latter case, i.e., the fast waves, are of interest for understanding the heating mechanisms of the outer solar atmosphere, while the slow waves may contribute to the heating of the low-to-mid chromosphere.

While each of these waves may propagate along a flux tube, coupling between different modes may also occur [?]. Numerical simulations of wave propagation in the lower solar atmosphere have shown coupling of fast and slow magneto-acoustic waves in regions where the acoustic and Alfvén speeds nearly coincide, i.e., where the plasma- (; [?]). This level has been called the “equipartition level” [?], where equipartition between magnetic and thermal energy density is achieved. Away from the equipartition level, the waves with magnetic field-dominated and acoustic natures are decoupled. The modes of such waves in low and high media are summarized in Table of [?]. At the equipartition level, part of the energy contained in the acoustic branch can be channeled to the magnetic branch and vice versa, while, the waves are also partly transmitted through the conversion layer without changing their physical natures, but with exchanging the fast and slow labels. The former and latter cases are normally referred to as “mode conversion” and “mode transmission”, respectively [?]. The fraction of energy transmitted or converted depends on the attack angle of the wave (i.e., the angle between the wave and the magnetic field vectors at the transmission/conversion layer) itself, the wavelength, and the width of the conversion layer [?]. In addition to the fast and slow waves whose physical nature depends on the level of plasma-, (transverse) Alfvén waves may propagate in both and regions [?].

Both incompressible (kink) and compressible (e.g., sausage) waves are thought to be excited in magnetic elements through the interaction of flux tubes with surrounding granules [?]. By exploiting the high spatial resolution images provided by Sunrise/IMaX, [?] have reported the interaction between longitudinal and horizontal-velocity oscillations in small magnetic elements in the solar photosphere. Also, propagation of kink waves in small magnetic elements has been shown to be a nonlinear process [?].

[?] and [?] proposed that rapid, pulse-like motions of the flux-tube footpoints produce kink waves along the tubes that can propagate into the upper solar atmosphere. The energy that such jerky motions can potentially carry may contribute to the heating of the quiet corona [?]. [?] found that the energy of incompressible, transverse kink waves, propagating along flux tubes, are transformed into the volume filling Alfvén waves above the height where individual flux tubes merge (i.e., in the chromosphere). However, according to their model, of such Alfvénic waves are reflected at the transition region. Fast magneto-acoustic waves may also couple to Alfvén waves at the apex of their refractive path [?]. This implies a non-homogeneous field or horizontally non-homogeneous density. The propagation of magneto-acoustic and kink waves in small magnetic elements (in the lower solar atmosphere) have been studied in detail in 2D [?] and 3D [?] simulations.

We note that waves reviewed here propagate in the lower solar atmosphere. Waves in the upper layers of the atmosphere have also been observed, e.g., in coronal loops ([?]; see [?] for a review), in quiet-Sun with frequencies up to mHz using the *TRACE* spacecraft [?], or in spicules from Hinode/SOT with an average frequency of mHz ([?], and many others).

In addition, different excitation mechanisms have been proposed or speculated upon for the same wave mode observed at different heights from the solar surface. For instance, the excitation of kink waves observed in the lower solar atmosphere has been attributed to buffeting of the flux tubes by the surrounding granular flows [?], while, small-scale magnetic reconnection in the chromosphere has been proposed to drive the kink modes detected in the higher atmosphere [?]. [?] have proposed a mechanism called “magnetic pumping”, where convective downdrafts around a flux tube pump downflows inside the tube. This mechanism was shown to result in the upward propagation of (slow) magneto-acoustic waves in magnetic flux concentrations [?].

We also note that we have only focused on waves propagation in magnetic elements. In non-magnetized areas, a rich spectrum of waves such as acoustic, gravity, and surface gravity (generating resonant modes of oscillation as -modes, -modes, and -modes below the solar surface, respectively) exists [?]. Propagation of, e.g., acoustic waves in the non-magnetized atmosphere transports considerable energy flux that contributes to the heating of the chromosphere [?].

In this paper, we investigate the propagation of high-frequency waves in the lower solar atmosphere by analyzing both horizontal-displacement oscillations of MBPs and their intensity perturbations at two sampled heights observed at high spatial and temporal resolution with the Sunrise balloon-borne solar observatory. We detect high-frequency fast waves, both compressible and incompressible, propagating at the selected MBPs in the lower solar atmosphere. In Section 2 we introduce the data used in this analysis. In Section 3 we describe our analysis methods that retrieve the observational properties of the waves. We provide the results and corresponding interpretations in Section 4. The results are discussed in Section 5 and conclusions are drawn in Section 6.

## 2Observations

### 2.1Observational Data

High spatial and temporal resolution observations of the quiet-Sun disk center at 300 nm (FWHM nm) and Ca ii H nm (FWHM nm) were carried out using the Sunrise Filter Imager (SuFI; [?]) on board the Sunrise balloon-borne solar observatory [?] on 2009 June 9 (between 01:32 UT and 02:00 UT). Simultaneous full Stokes observations in the magnetically sensitive Fe i nm line recorded by the Imaging Magnetograph eXperiment (IMaX; [?]) on board Sunrise provided the photospheric magnetograms that were used to determine the magnetic properties of the features observed in the other passbands.

The seeing-free 300 nm and Ca ii H image sequences employed in this work share a common field of view of and were both recorded with a cadence of s and a spatial sampling of arcsec/pixel. The recordings at the two wavelengths are offset by s in time (300 nm images follow Ca ii H filtergrams). All SuFI images were phase diversity reconstructed with averaged wave front errors (i.e., level 3 data; [?]). Moreover, the IMaX and the two SuFI channels were aligned with sub-pixel accuracy (i.e., better than 14 km), using an approach that utilizes cross-correlation and mean squared deviation techniques applied to common sets of manually selected features (see [?] for a description of the alignment procedure).

Figure 1 shows a small part of a snapshot in each wavelength. An MBP is marked whose horizontal-displacement fluctuations as well as intensity oscillations are studied here.

### 2.2Formation Heights

We estimate the formation height by two completely independent means. We begin by estimating the heights of formation of the 300 nm band and Ca ii H line profile by computing their contribution functions (CFs) using the RH radiative transfer code [?]. The code provides the contribution to the emission as a function of height at a certain wavelength by solving both radiative transfer and statistical equilibrium equations in a given atmospheric model. The CFs for the 300 nm passband are computed in LTE (Local Thermodynamic Equilibrium) conditions. The calculations are in non-LTE for the Ca ii H line, taking partial redistribution into account [?]. A five-levels Ca ii H model atom (i.e., with Ca ii H/K and Ca ii infrared triplet lines; [?]) was used in the latter computations. Following the discussions of [?] (hereafter referred to as J13), we use the FALP model atmosphere [?] to describe the MBPs after convolving the CFs at different wavelengths by the spectral profile of the relevant SuFI filter. For more details of similar implementations, we refer the reader to section in J13. For comparison, the CFs resulting from the FALC model atmosphere (that represents an averaged quite-Sun region) are also determined.

Plotted in Figure 2 are the CFs for the two passbands and both atmospheric models. We show this figure here, although some of the results can already be found in J13, in order to highlight the difference between the CFs and the average formation heights of the 300 nm and Ca ii H passbands of Sunrise/SuFI. We note that the SuFI Ca ii H images clearly have a large contribution from the photosphere, with some contribution from the low-to-mid-chromosphere. The vertical dashed lines indicate the corresponding average formation heights. From the mean heights of formation, an average height difference of is determined between images observed in the 300 nm and the Ca ii H passbands. The sudden drop of the Ca ii H CF for the FALP at km is due to the rapid temperature increase in the upper chromosphere, leading to the ionization of Ca.

The main uncertainty in these results lies in the choice of the atmospheres. Clearly, the choice of a 1D atmosphere is a gross simplification, because the real Sun shows a rich variety of fine structure and the CFs are affected nonlinearly by changes in the atmospheres. An additional uncertainty is posed by the difference between the CFs and the response functions to individual atmospheric parameters, though response functions in non-LTE are not yet available, so that, this point is somewhat moot.

However, the results shown in Figure 2 can be understood with the help of simple physical considerations. A Wilson depression is present inside the flux tubes, so that radiation is emitted from deeper layers inside them. However, we are dealing with the height difference, rather than the absolute heights of formation. The Wilson depression occurs at both heights, making it irrelevant for the difference in formation heights. More important is the run of pressure with height. For a higher temperature inside the magnetic feature, as is required for it to be an MBP, the pressure drops more slowly with height within the magnetic feature than outside it, which increases the difference in formation height, in agreement with the CFs. The formation heights of photospheric spectral lines weakened by a higher temperature, decreases within magnetic elements [?]. However, this is not the case for the Ca ii H line core because the increased chromospheric temperature rise within magnetic elements causes the emission peaks around the core of this line to strengthen [?], so that the contribution from the chromosphere to this line increases (as indicated by the CF computations shown in Figure 2). The formation height of the 300 nm continuum also increases slightly when going from quiet-Sun to MBP, as H opacity increases with temperature. This increase is not large because MBPs are not much brighter than the average quiet-Sun at low heights. Thus we expect the difference in height of formation between the two wavelengths to increase somewhat from quiet-Sun to MBP, which supports the CF computations.

Nonetheless, we check the results obtained from the CFs by employing an entirely independent technique. To this end, we apply an analysis based on the Fourier transform of the intensity oscillations of both, the 300 nm and Ca ii H time series in a relatively large quiet-Sun FOV (i.e., ). This area was chosen because upon visual inspection it harboured only a relatively small number of MBPs (and no other larger magnetic structures) for the entire time series.

We perform the Fourier transform on the entire length of the time series, separately for each individual pixel, and compute phase-lags between the time series from their cross spectra. In the averaged quiet-Sun, acoustic waves are expected to be dominant, with km s in the photosphere. Waves propagating in a magnetic element of the chosen quiet-Sun FOV may behave differently. However, their contribution to the obtained phase-lag from the entire FOV is very small, since they only occur in a very small fraction of it. The phase-lag () can then be converted into a time lag () for each Fourier frequency () using

Figure 3 illustrates a 2D histogram of phase versus frequency. It shows the unwrapped phase-diagram, i.e., the phase-diagram that is corrected for phase jumps of . Acoustic waves with frequencies larger than the (corresponding) cut-off frequency can propagate within the atmosphere and indeed we found acoustic waves in the high-frequency range of mHz that are propagating upward. The lower boundary of the frequency range ended at the cut-off frequency of the acoustic waves, the upper boundary at the Nyquist frequency ( mHz). The diagonal ridge outlines a nearly linear trend. We presume that this ridge is due to acoustic waves. The slope of this ridge is given by the best-fit black, solid line. The yellow error bars indicate the standard deviation of the data points along the phase-lag axis. The slope reflects the travel time of the acoustic wave between the two layers and is equal to s, which corresponds to a height difference of km between the layers of the quiet solar atmosphere sampled by the 300 nm and Ca ii H passbands. The error comes from the uncertainty in determining the slope and in the sound speed of km s.

For comparison, we have over-plotted the expected phase difference () versus frequency () for linear (vertical) wave propagation between the two layers, based on three different models explained by [?]:

Adiabatic propagation in an isothermal non-stratified atmosphere (red, dashed line; Equation (Equation 1)).

In a gravitationally stratified atmosphere (white, dotted–dashed line):

where , is the acoustic cut-off frequency, m s is the gravity (assumed to be constant), and is the ratio of the specific heats for adiabatic propagation.

Non-adiabatic propagation in a gravitationally stratified atmosphere when radiative losses have been taken into account with Newton’s law of cooling (dark blue, long dashed line; originally developed by [?]):

where

is the pressure scale height and is the radiative time-scale.

We obtained the best fits of the three model curves to the observed trend using km s, s, and km, which is in agreement with the height difference of km from the best fit to the data points.

The height difference of km between the two atmospheric layers sampled by the Sunrise/SuFI Ca ii H and 300 nm passbands, obtained remarkably from analyzing the acoustic waves in the quiet atmosphere, agrees with that obtained from the CFs (i.e., km in the quiet-Sun). This gratifying agreement increases our trust in reliability of the formation height differences deduced from the CFs. For the atmosphere inside the magnetic elements, we therefore adopt the value returned by the CFs, 450 km, for this difference, but impose a conservative uncertainty of km.

## 3Analysis and a Case Study

We study wave-like phenomena in the lower solar atmosphere by analyzing oscillations in both the horizontal displacement and the intensity in small-scale magnetic features, the intensity oscillations being a tool to investigate the presence of compressive modes in the magnetic elements. Knowledge of the precise position of such elements versus time is a must in such analyses, which require stable (preferably seeing-free) observations. We focus on small magnetic bright points (i.e., MBPs), similar to those studied in J13. In short, the MBPs studied here were required to meet the following criteria: they had to (1) be located in an internetwork area, (2) be brighter than the mean intensity of the entire frame, (3) have a magnetic nature (coinciding Stokes patches with S/N), (4) live longer than five minutes (to avoid short-lived brightenings due to, e.g., high frequency oscillations or reversed granulation at the heights sampled by the Ca ii H passband; [?]), and (5) be point-like features (with a diameter smaller than 0.3 arcsec). In this way, we ensure that the selected Ca ii H MBPs differ from non-magnetic H or K grains, which are mainly due to acoustic waves [?]. The availability of polarimetric data enables us to identify magnetic fields and to distinguish between MBPs and other localized brightenings, different from other studies, such as [?]. Also, the the so-called persistent flashers (whose brightness drops below the detection limit; similar to those shown in J13) are not included in this study.

Consequently, the accurate locations as well as the intensity of such MBPs at any given time are determined using the same algorithm as described in detail in J13. We search for MBPs whose horizontal displacement and intensity show an oscillatory behavior, such as the examples presented in J13. In order to facilitate the precise localization and tracking of the MBPs, noise and extended solar brightenings were eliminated from both sets of images (i.e., granules from the 300 nm and other brigthenings due to, e.g., shock waves from the Ca ii H filtergrams). This approach returns the locations of the MBPs with an accuracy better than 0.5 pixel (). See section 3.1 of J13 for details on the image processing, detection and tracking algorithms employed here.

In order to determine whether any of the detected waves propagate within the solar atmosphere (and to measure the speed of propagation), we need to simultaneously trace horizontal displacements and intensity oscillations of the same magnetic elements at two atmospheric heights.

Thus, we found seven MBPs whose trajectories could be precisely tracked in both the 300 nm and the Ca ii H image sequences for a sufficiently long time, i.e., longer than five minutes. This latter criterion improves the frequency resolution in the observed oscillations The amplitudes of the oscillations in the seven MBPs at the two atmospheric heights, along with the periods and phase angles of the detected waves, are summarized in Table 1.

The number of MBPs investigated here is lower than the Sunrise Ca ii H BPs (with lifetimes longer than five minutes) studied in J13 because of the relatively low contrast of the 300 nm filtergrams compared to the Ca ii H images [?]. We dropped all of those MBPs that could not be precisely tracked in both filtergrams over their entire lifetimes. By restricting ourselves to the smallest MBPs, we may have discarded the longer lived ones (if larger MBPs have longer lifetimes, which, however, is not the case according to [?] if one includes splitting and merging among the causes of death). Due to the careful choice of MBPs, coupled with the seeing-free high-resolution data, we expect to get reliable results on the waves within the MBPs. Nonetheless, given the small number of MBPs in our sample, the results we obtained cannot be taken to represent the properties of quiet-Sun MBPs in general.

We note that, although uncertainties introduced by the locating algorithm as well as instrumental vibration-induced image jitter may bias the final analysis, their effects are much smaller than the motions of the MBPs under study. [?] found an average uncertainty of km s in the determined horizontal velocities of similar MBPs. This is much lower than the mean amplitudes of the horizontal velocities of both 300 nm and Ca ii H MBPs (i.e., 1.8 km s and 1.2 km s, respectively; see Table 1) in the present study. The power spectra of the residual image jitter (in both horizontal and vertical directions) measured for the Sunrise observatory show a frequency range of about [?], which is much higher than the sub-Hertz frequencies of the waves under study (see Section 4). Hence, the effect of instrumental vibration-induced jitter is negligible in our analysis, except for smearing the images. Moreover, we observe clear phase-lags between the oscillations observed in the two heights (see Section 3.2) that cannot be due to jitter or random motions.

The maximum intensity over all pixels of an MBP at any given time is measured as its intensity at that time step. [?] obtained an average photon noise of in the Ca ii H images as a mean uncertainty in determining the intensity. This noise level is much smaller (by a factor of ) than the mean amplitude of the Ca ii H intensity oscillations.

The s time difference between the two time series (i.e., 300 nm and Ca ii H image sequences) is corrected for by adding the corresponding phase angle to the final results.

### 3.1Wavelet Transform

We perform a wavelet analysis [?] in order to simultaneously localize the spectral power in both time and frequency domains. The wavelet transform is especially suitable for searching transient oscillations and for studying wave propagation within the solar atmosphere in the presence of short-lived features and waves [?].

The wavelet () is defined as the convolution of a time series with a “mother” function that is a window (envelope) whose variable width enables the analysis to capture both low/high frequencies and long/short durations simultaneously. We employ the Morlet mother function with a dimensionless frequency using a wavelet algorithm developed by [?]. This function satisfies the balance between frequency and time localization [?] and hence is suitable for investigating the propagation of waves with different ranges of frequencies. The Morlet function is a complex wavelet, consisting of a plane wave modulated by a Gaussian window, described as

where is a dimensionless time that stretches the wavelet in time by changing its scale .

The wavelet power spectrum of a time series is defined as . The cross-spectrum (or cross wavelet power spectrum) of two time series is then determined by multiplying the wavelet power spectrum of a time series with the complex conjugate of the other one. Furthermore, interaction between the two time series can be examined using a bivariate framework called wavelet coherence, which is the square of the cross-spectrum normalized by the individual power spectra [?]. The coherence level varies between zero and one representing incoherent and coherent oscillations, respectively.

[-2em] | |||||||||||||||||||

[-1.4em] | Lifetime | ||||||||||||||||||

(s) | |||||||||||||||||||

NO. | Passband | mean | peak | mean | median | mean | median | mean | peak | mean | median | mean | median | ||||||

300 nm | 528 | 3.7 | 45 | 114 | 72 | 118 | 40 | ||||||||||||

Ca ii H | 2.1 | 768 | |||||||||||||||||

300 nm | 552 | 7.8 | 59 | 67 | 99 | 140 | 79 | ||||||||||||

Ca ii H | 2.6 | 1306 | |||||||||||||||||

300 nm | 408 | 6.8 | 50 | 106 | 189 | 125 | 31 | ||||||||||||

Ca ii H | 3.4 | 882 | |||||||||||||||||

300 nm | 624 | 6.5 | 51 | 90 | 191 | 99 | 37 | ||||||||||||

Ca ii H | 7.2 | 978 | |||||||||||||||||

300 nm | 480 | 5.9 | 53 | 109 | 229 | 92 | 15 | ||||||||||||

Ca ii H | 3.3 | 705 | |||||||||||||||||

300 nm | 396 | 2.6 | 53 | 96 | 43 | 109 | 45 | ||||||||||||

Ca ii H | 3.2 | 840 | |||||||||||||||||

300 nm | 420 | 3.6 | 47 | 78 | 101 | 106 | 35 | ||||||||||||

Ca ii H | 4.3 | 1270 | |||||||||||||||||

Notes.

The cross-spectrum highlights time-frequency areas with high common power in the two time series, whereas the wavelet coherence detects regions in a time-frequency domain where the examined time series co-move, but do not necessarily possess a strong common power. Hence, while the cross-spectrum can provide sufficient information on oscillatory behaviors in a localized medium (by representing a local co-variance between two time series), the wavelet coherence is also needed for finding co-movements between perturbations at different regions (heights) in the solar atmosphere.

Since the wavelet transform of a time series with a finite length has edge artifacts (the wavelet is not totally localized in time), a “cone of influence” (COI) is introduced. Thus, the edge effects inside the COI cannot be ignored. Following [?], the COI is defined as regions where the wavelet power spectra from a discontinuity at the edge has reached e of the value at the edge.

Finally, a phase difference between a pair of time series provides information on delays in the oscillation (i.e., on wave propagation). These phase-lags are estimated from the complex and real arguments of the cross spectra.

### 3.2A Case Study: Wavelet Analysis of an MBP

In the present study, we compute wavelet power spectra of four time series: the horizontal displacement and the intensity oscillations at two sampled heights (i.e., the two wavelength bands) for a given MBP. This is done for all MBPs, but here we discuss the results in greater detail for one example. Then, we determine cross power and wavelet coherence between the power spectra of the two time series of a given type of oscillation sampled corresponding to the two atmospheric layers.

Figure 4 illustrates the wavelet power spectra of the horizontal-displacement perturbations and the intensity oscillations observed in both the Ca ii H and 300 nm passbands for the MBP highlighted in Figure 1. The background colors represent the power normalized to its maximum value. The cross-hatched areas show the COI, which defines regions that are subject to the edge effect. Hence, only the computed values outside the COI are considered for the phase analysis. The black contours indicate the 95% confidence level.

The MBP under study has a lifetime of s (for details on determining the lifetime and its uncertainty of similar MBPs we refer the reader to Section of J13). The wavelet analysis of both types of oscillations reveals periodicities with a wide range of values ( s) for this magnetic element. While the intensity perturbation of the MBP has much of the power in the period range of s, its horizontal-displacement oscillates within a higher frequency range (periods of 25-50 s). Note that the lower limit of the period range is slightly larger than the shortest detectable period of 24 s that corresponds to the Nyquist frequency of our image sequences.

Figures Figure 4(a) and Figure 4(c) show several power peaks, indicating the occurrence of horizontal-displacement oscillations (observed at the heights sampled by Ca ii H and 300 nm, respectively) with periods changing over time. By visual inspection of these power maps, we can find two high power patches in the period range s and within the time interval 20-250 s. These seem to occur in Figure 4(a) (corresponding to the oscillations observed in Ca ii H) after equivalent features in Figure 4(c) (representing the oscillations seen in the 300 nm bandpass). Another patch with similar behavior (but with weaker power) in the same period range at around time 500 s is also observed. The power patches in intensity oscillations, displayed in Figures Figure 4(b) and Figure 4(d), show that much of the power is concentrated at longer periods than that of the horizontal-velocity perturbations. We note that in all plots only small areas outside the COI correspond to high power values. However, they include patches of sufficiently large power for the oscillatory motions to be real (i.e. lying over the 95% confidence level). Moreover, we are interested in the correlations of the oscillations between the two sampled heights.

The correlations between the power spectra shown in Figures Figure 4(a) and Figure 4(c) and those in Figures Figure 4(b) and Figure 4(d) are illustrated in Figures Figure 5(a) and Figure 5(b), respectively. The cross spectra represent the regions of common high powers in the time-frequency domain. The contours mark the 95% confidence level. The COI is shown as areas with bleached colors.

The phase angles between the two atmospheric layers can be deduced from the small arrows in both panels of Figure 5. Arrows pointing to the right and to the left show in-phase and anti-phase oscillations between the two layers, respectively (i.e., representing standing waves). Arrows pointing straight down indicate that the oscillations in Ca ii H lead the oscillations in 300 nm images by , which would imply downward propagating waves.

The cross-spectrum of displacement oscillations clearly shows the two common patches of high power in Figures Figure 4(a) and Figure 4(c) mentioned earlier. However, the arrows indicate that only the leftmost patch corresponds to upwardly propagating waves. The rightmost patch in Figure 5(a) represents waves moving downward in the solar atmosphere. The intensity perturbations of common high power and long periods indicate upwardly propagating waves of a small phase shift.

The wavelet coherence of the two types of oscillations is shown in Figure 6. To ensure the reliability of the results, only values within contours of 95% confidence are included for the phase analysis. The confidence level is computed using a Monte Carlo method, because the statistical behavior of the wavelet coherence is unknown [?]. We note that relatively large values of wavelet coherence do not necessarily intimate significant co-movement between the two time series. The wavelet coherence attains large, significant values in areas which do not necessarily represent high power patches in the individual wavelet power spectra, in particular for the horizontal displacements. Therefore, the individual power spectra are not a particularly good guide to the presence of correlated oscillations at the two heights; the wavelet coherence provides better information on wave propagation.

The phase angles are also indicated on Figures Figure 6(a) and Figure 6(b), but for simplicity, only for areas with a coherence value larger than 0.5. In both panels, Figure 6a and Figure 6b, and in all regions within contours of 95% confidence that are outside the COI, the arrows tend to point upward (i.e., positive phase-angle). This indicates that the oscillations observed in the Ca ii H passband follow the ones seen in the 300 nm filtergrams, meaning that both types of waves propagate upward in the magnetic element under study.

Different angles of the upward pointing arrows represent different phase angles between the oscillations in the two layers which can be converted to time lags (i.e., the wave travel time) for specific frequencies using Equation 1.

We compute the time lags for all points of the wavelet-coherence map (e.g., Figure 6) which are located outside the COI and have a confidence level of 95% or higher.

The propagation speed of the waves () at a given frequency can then be calculated from its corresponding time lag () and the height difference () between the two atmospheric layers, i.e., .

## 4Statistics and Wave Properties

Table 1 summarizes properties of the detected oscillations and waves in the individual seven MBPs studied here. The lifetime represents the duration along which an MBP is simultaneously observed in both the 300 nm and Ca ii H images. The amplitudes of both types of oscillations, i.e., fluctuations in displacement (horizontal velocity) and in intensity of the MBPs, are also provided. In agreement with J13, the MBPs under study turned out to have high-velocity excursions (i.e., large horizontal-velocity amplitudes) over the course of their lifetimes. Such rapid “pulses” have been shown to excite kink waves in magnetic elements [?]. The MBPs detected in the 300 nm images have, on average, larger horizontal-velocity amplitudes than their Ca ii H counterparts. The mean and median values of period and of phase-angle for each oscillation observed in the individual MBPs are obtained from their computed wavelet coherence, as long as they are located outside the COI and within contours of 95% confidence (see Section 3.2).

Figures Figure 7(a) and Figure 7(b) are 2D histograms of phase-angle versus period from the wavelet coherence of all seven investigated MBPs of the horizontal velocity and the intensity oscillations, respectively. All phase-angle values obtained from all seven analyzed MBPs enter Figures Figure 7(a) and Figure 7(b), as long as they are located outside the COI and within contours of 95% confidence. It turns out that all these phase angles have a coherence value exceeding . Examples are the phase-angles obtained from the wavelet coherence spectra satisfying the above criteria in Figures Figure 6(a) and Figure 6(b).

The background colors in Figures Figure 7(a) and Figure 7(b) represent the sample density. For robust statistics, unbiased by outliers with a small probability of occurrence (due to, e.g., spurious coherence between the two oscillations), the extreme outliers, i.e., values lying outside the clusters of most of the other data points, are determined. We separate concentrated regions from the extreme outliers using the Grubbs statistical significance test [?]. This test estimates a criterion corresponding to the largest deviation from the mean value in units of the standard deviation, , for the data points in the 2D histograms. Regions with a number density smaller than this criterion are considered to host the outliers whose small probability of occurrence compared to the rich cluster of other data points may lead to biased interpretations [?]. The solid-line contours in Figures Figure 7(a) and Figure 7(b) indicate the estimated Grubbs’ criteria of and , respectively, from the mean values of the concentrated regions. The contours include 1032 and 2274 individual data points obtained from the wavelet coherence of the seven MBPs, for the horizontal displacement and intensity oscillations, respectively.

We find short-period (high-frequency) oscillations with periods (1) between 43 s and 74 s for horizontal displacements and (2) between s and s for intensity oscillations. The phase-angle distributions show almost only positive values, meaning the propagation of wave-like phenomena from the height sampled by the 300 nm passband (i.e., the lower layer) toward the height sampled by the Ca ii H filter. However, about of all occurrences seen in the intensity oscillations have a negative phase-angle, which can represent a downward propagation of the perturbations. Also, about 1% of the data points related to the intensity oscillations have a zero phase lag (standing mode).

The wave travel times corresponding to the highest peaks of the 2D histograms in Figures Figure 7(a) and Figure 7(b) were computed (using Equation (Equation 1)) as s and s for the horizontal displacement and intensity oscillations, respectively.

Both phase diagrams, particularly that of the intensity oscillations, show patches of high occurrences at distinct period ranges. We cannot, however, verify whether these separated regions represent distinct dispersion relations due to, e.g., different waves with different modes or natures.

We perform a test to inspect whether the time delays corresponding to the isolated islands in Figure 7(b) are compatible with what we expect from acoustic waves. We follow [?] and [?] and assume that the trend connecting the peaks in Figure 7(b), or a part of that, represents the acoustic dispersion relation. Consequently, the height difference between the two layers can be estimated using the wave travel time (using Equation (Equation 1)) and a sound speed of km s. In this case, we obtain km, a height difference that is small compared to the km estimated from radiative transfer (see Section 2.2). A somewhat larger , which may be found in a hot magnetic element, is not able to resolve this discrepancy. For comparison, the dispersion relation expected for acoustic waves traveling between the two atmospheric layers (with a height difference of km) is also plotted in Figure 7(b) (the dashed curve). To illustrate the effect of any uncertainty in the formation height of Ca ii H, the gray shaded area around the dashed curve shows how strongly the latter dispersion relation would change for a range of km around the 450 km height difference. The acoustic phase-lag is obviously much larger than that observed in the MBPs at all frequencies. This means that either the waves in the magnetic features propagate much faster than sound, or the comparable height difference determined using two independent approaches (i.e., from the CFs, and from the Fourier analysis of acoustic waves; see Section 2.2) are wrong. Alternatively, the short wave travel times could represent wave fronts traveling obliquely to the line of sight [?]. This would imply observations of phase speeds and not true propagation speeds of waves (see Section 5).

Since the difference in formation height between 300 nm and Ca ii H increases slightly when going to a hotter atmosphere (see Section 2.2), we conclude that the problem lies in assigning the intensity oscillations in the MBPs to acoustic or slow-mode magneto-acoustic waves.

Using a height difference of km and the wave travel times determined from the 2D histograms shown in Figure 7, we estimate an upward propagation speed of km s for the horizontal displacements and a value of km s for the intensity oscillations. If one identifies these velocities as kink or Alfvén speeds ( or , respectively), then, using

and

we can obtain a rough estimate of the field strength. In Equations (Equation 6) and (Equation 7), is the field strength in the flux tube, g/cm and g/cm are the average gas densities (averaged over the atmosphere between the two sampled heights, from FALP and FALC models) inside and outside magnetic elements, respectively. Although the Alfvén wave is incompressible, for simplicity, we use the Alfvén speed as a lower limit for the speed of the fast magnetosonic wave (they would be identical in the cold plasma limit). Applying these formulae gives kG and 1.5 kG for the kink and the Alfvén waves, respectively. The propagation speed of the kink waves resulted in an unrealistically large field strength, which suggests that either the true phase speed is lower than the value determined here, or that the wave is not a pure kink-mode wave, so that Equation (Equation 6) is not valid. We should, however, note that the Bifrost MHD simulations [?] indicate an order of magnitude lower gas densities, averaged over the atmosphere within the 500 km from the solar surface at optical depth unity (i.e., the atmosphere between the two sampled heights). Such diminished densities result in kG and 0.5 kG from the observed speeds of the kink and the Alfvén waves, respectively.

We note that the above estimations of the propagation speeds and of the field strengths are based on the time lags determined from the peaks of the phase diagrams. Patches with weaker significance in Figure 7 would result in much larger phase speeds (particularly for the intensity oscillations) of up to 200 km s. The latter velocities would be too large to be interpreted as a true propagation speed of MHD waves in a magnetic element (they would result in a much too large field strength if approximated with, e.g., Equations (Equation 6) or (Equation 7)). These seemingly very high phase speeds are likely statistically insignificant, and need not concern us too much here.

In the next section, we review a few relevant theoretical and numerical studies, based on which we will discuss our interpretations of, and speculations on, the observed travel times.

## 5Comparisons and Discussion

Theoretical models and numerical simulations of magneto-acoustic and kink waves in photospheric magnetic flux concentrations may be a better way for the interpretation of the observed propagating waves than the simple estimates that we have made in the previous section. However, we should keep in mind that such theoretical investigations have often been confined to a simplified atmosphere whose characteristics may differ, to some extent, from the actual solar atmosphere. We review a few relevant models/simulations in the following that may provide a better understanding of our observations of the propagating, high-frequency, fast transverse and longitudinal waves in small magnetic elements from Sunrise/SuFI. We note that the primary drivers of perturbations in the following papers are often motivated by observations of motions as well as the brightness of magnetic bright points.

*Wave excitation models by [?]*: They modeled the transverse velocity of magnetic flux tubes excited by their footpoint motion by analytically solving the Klein-Gordon equation for kink waves. They showed that the excitation of kink waves, due to short-duration pulses at the base of the flux tubes, can produce intermittent chromospheric emissions. Such kink waves have been shown to potentially carry enough energy to contribute to coronal heating [?]. [?] found that these transverse waves cannot represent a major contribution to chromospheric heating unless they are excited by high-frequency motions (periods of s). They speculated that such high-frequency motions could be due to intergranular turbulence (below the photospheric base) that would not be observable from the ground due to the influence of seeing.

The short periods of s that we observed in horizontal-velocity perturbations overlap with the range of periods found necessary by [?] for the transverse (kink) waves to transport enough energy to heat the chromosphere.

*Numerical simulations by [?]*: As a continuation of the work by [?], [?] extensively studied the propagation of magneto-acoustic-gravity waves within the photospheric and low chromospheric regions of a 2D magnetized atmosphere. In their numerical simulations, uniform waves are generated at a source in the photosphere, that is confined to a km region within the flux element ( km below the surface), with a driving frequency of mHz (a period of s). The propagation of horizontal and of vertical perturbations were studied separately. In addition, they investigated wave propagation in magnetic elements with different field strengths, such that the layer was located below or above the wave sources. They concluded that the slow and fast waves are decoupled in the low- and high- parts of the atmosphere but that they couple, leading to mode conversion and/or mode transmission, in areas with , i.e., where the gas and magnetic pressures are comparable. They also discussed that the intermediate Alfvén waves, that propagate in both low- and high- regions, may remain coupled with the fast mode through the whole domains. They approximated the structure of the magnetic field by a potential field, thus ignoring the current sheet that separates magnetic flux concentrations from their surroundings.

*MHD models by [?]*: As a continuation of earlier works by [?] and [?], [?] proposed that the intensity enhancement of Ca ii H MBPs are due to a heating process caused by the dissipation of high-frequency, slow magneto-acoustic waves (i.e., frequencies mHz) launched as kink waves at the base of the magnetic flux tubes. Propagation of both slow and fast waves were indeed observed in their simulations. They explained the generation of both types of waves in terms of mode conversion and mode transmission at the level where (1) horizontal (transverse) motions at the base of the photospheric flux tubes (producing slow, MHD, kink waves in regions) are partially converted to slow (longitudinal) acoustic waves propagating along magnetic field lines in regions, (2) slow, MHD, transverse waves in the medium are partially transmitted as transverse, MHD, fast waves in the regions, and (3) acoustic waves generated in the ambient medium are converted into transverse fast modes. [?] showed that while the fast transverse waves rapidly travel to the upper layers of the atmosphere, the slow, compressive acoustic waves form shocks at chromospheric heights accompanied by large temperature perturbations (around ) representing the intensity enhancement of the Ca ii H MBPs. We conjecture that the high-frequency kink as well as the longitudinal waves we observed here (or at least a fraction of them) are a result of such mode conversion and/or mode transmission processes explained by, e.g., [?] and [?].

*Numerical simulations by [?]*: they investigated the unexpectedly small time lags measured between two photospheric layers, similar to those observed by, e.g., [?]. They consider a dynamic atmosphere (i.e., convectively unstable, time-dependent and magnetic) in their numerical simulations and investigate the propagation of both spherical and plane-parallel waves, excited from different locations with respect to small magnetic elements. They found that when one of the height levels (in the solar photosphere at which the waves are observed), or even both of them, are above the level, the wave travel time significantly decreases. They explain such a decrease in time lag as the effects of: (1) convergence of formation heights inside the strong magnetic elements, (2) conversion to and detection of the fast magneto-acoustic wave above the level, and (3) the refractive propagation path of the fast magneto-acoustic wave, which can lead to very high phase speeds. The refractive propagation is produced in the presence of the inhomogeneous magnetic field, typically associated with magnetic elements. Thus a fast magneto-acoustic wave that starts out propagating upward inside a part of a flux tube, will with time and height turn into a more horizontally directed propagation path. In this case, the phase speeds would be miscalculated since the observed travel time corresponds to an obliquely propagating wave with wave fronts straddling both the upper and the lower height levels.

[?] show that the longitudinal oscillations include both slow acoustic and fast magnetic modes in the region. Unlike [?], who primarily launched the wave by displacing the base of the flux tube in their MHD model, the wave source in the simulations of [?] only has an acoustic nature. Therefore, the incident acoustic wave at the level, is partially converted to the fast magneto-acoustic wave and is partially transmitted as the slow acoustic wave propagating to the regions. [?] also noted that since the dynamic and complex level in the actual solar atmosphere may cause several levels above each other, several mode conversions may occur.

With respect to the quiet-Sun, the formation height inside magnetic elements of the 300 nm bandpass is expected to decrease due to spectral line weakening in magnetic elements, whereas it is expected to increase due to the enhanced brightness of the Ca ii H line core in magnetic elements, which strengthens the line emission [?]. Hence, we expect the formation height to diverge rather than to converge when going from quiet-Sun to a flux-tube atmosphere. The refractive wave path could be responsible for the seemingly short travel times and the resulting determined high phase speeds of the fast mode in our MBPs.

Summarizing, simulations have shown that (mainly) two types of wave modes (i.e., fast and slow) are generated at the base of small-scale magnetic flux tubes due to impacts from the side or lateral shaking while moving within intergranular lanes [?]. Bursts of perturbations can be caused by, e.g., expansion or explosion of neighboring granules as well as by intergranular turbulence, with a large range of speeds. The turbulence in the intergranular areas has been speculated to be responsible for the high-frequency range of the generated waves [?]. Pulse-like events (sometimes super-sonic; [?]) can excite transverse kink waves traveling along the flux tubes. In addition, the interaction of the magnetic elements with the convective flows generates magneto-acoustic waves that propagate along the field lines.

The following waves can propagate in the region: (1) slow, longitudinal, acoustic waves, (2) fast, transverse, magnetic (kink) waves, and (3) fast, longitudinal magnetic waves. Therefore, we may interpret our observed propagating high-frequency (fast) transverse and (fast) longitudinal waves as a result of mode conversions/transmissions.

When the magnetic field is slightly inclined, fast waves can be converted to the Alfvén mode in layers well above the level, i.e., at the fast wave reflection point [?]. The MBPs under study here do have small magnetic field inclinations with respect to the line of sight (with an average value of ), similar to those studied by [?]. If the reflection point of fast waves reaches the transition region, up to of the fast wave’s energy flux can be carried across the transition region by the Alfvén waves due to the fast-to-Alfvén mode conversion [?]. Thus, fast waves may contribute to the heating of the outer atmospheric layers. We cannot, however, verify whether the fast waves we observed in the magnetic elements can reach these atmospheric heights since our data do not sample those layers.

We note that this interpretation should be treated with caution. [?] explained characteristics of oscillations in the solar magneto-atmosphere using extensive numerical simulations. They showed that at any given location, a superposition of several distinct waves with different characteristics may be observed. These waves may propagate in different directions and may come from different locations, e.g., directly from their sources, from the equipartition level as a result of mode conversion, or locally from the interaction of various kinds of waves with -modes. Therefore, distinguishing the nature of different waves may not be straightforward from observations alone.

Similar to our findings, [?] have also observed both transverse (incompressible) and longitudinal (compressible) waves along magnetic flux tubes close to the quiet-Sun disk-center, sampled in the mid-chromosphere. They interpreted these two types of waves as () the fast MHD kink mode, measured from transverse displacements of the magnetic structures, and () the fast MHD sausage mode, determined from intensity perturbations. Our observations may represent the lower solar atmosphere’s origin of the fast waves observed in the mid to upper chromosphere by [?]. They found propagation speeds of km s for the kink oscillations and a range of km s for the longitudinal oscillations. The typical periods of s for the fast kink waves and of s for the fast longitudinal waves measured by [?] are somewhat larger than the ones we obtain here (with an overlap for the longitudinal oscillations). The size variations of our point-like magnetic elements are too small, to investigate the relationship between size and intensity variations to detect any manifestation of the sausage mode in the intensity perturbations. For somewhat larger magnetic features observed by Sunrise, however, [?] have noticed such oscillations.

The upward wave propagation in small MBPs, studied here, agrees with that of [?] who found a similar behavior in intensity oscillations of a larger number of MBPs in both observations and simulations. They investigated the wave propagation in all pixels of a relatively large FOV and concluded that much of the power is concentrated in MBPs. However, with the lower spatial resolution of their observations compared to Sunrise/SuFI, they found only low-frequency waves in the rage of mHz in intensity oscillations. The spatial resolution has been shown to have a direct correlation with the power of high-frequency oscillations, particularly on small spatial scales, with lower spatial resolution leading to a decreased sensitivity to power [?]. The factor of two longer time series of observational data used by [?] has consequently allowed them to obtain such low frequencies. However, [?] also found low-frequency intensity oscillations on the order of mHz from their MHD simulations with grid size of 25 km (i.e., spatial resolution of 50 km) and a similar frequency resolution as in our study.

[?] and [?] studied transverse oscillations and sausage-mode waves in bright, slender fibrils observed in a narrowband Ca ii H passband (FWHM nm) provided by SuFI during the second flight of the Sunrise observatory [?]. They found periods in the range of s for both types of waves, with median values of roughly 83 s and 34 s for the transverse and sausage waves, respectively. Their range of periods overlaps with those of the MBPs we observe in this study. The slender Ca ii H fibrils have been shown to outline the nearly horizontal magnetic field lines in the low solar chromosphere, slightly higher than the heights, where our Ca ii H MBPs were sampled [?]. Since such fibrils have footpoints in photospheric magnetic features, such as MBPs, we speculate that waves of the type observed here partly continue into the chromosphere and may become visible as oscillations of fibrils.

## 6Conclusions

We studied oscillations in the horizontal displacement and intensity of small-scale MBPs as observed in image sequences taken in the passbands at 300 nm and in the Ca ii H spectral line with Sunrise/SuFI. Time-series of seven MBPs reliably displaying oscillations of both types were analyzed. Although such a small number of MBPs may not be representative of quiet-Sun magnetic elements in general, the obtained properties of the oscillatory motions are free of effects introduced by seeing.

We computed a height difference of km between the two atmospheric layers sampled by the 300 nm and Ca ii H passbands by employing the RH radiative transfer code of [?]. A comparable height difference was obtained from the analysis of the intensity oscillations, likely due to acoustic waves, in a quiet-Sun area. The phase diagram of the propagating high-frequency (4-38 mHz) acoustic waves resulted in a time lag of s, implying an average height difference of km between the two atmospheric layers sampled by the Sunrise/SuFI 300 nm and Ca ii H passbands.

Our wavelet analysis of the small-scale magnetic elements yields: () consistent oscillations with high frequencies of up to mHz in intensity and up to mHz in horizontal displacement (not limited by the Nyquist frequency of mHz), () positive phase-lags between both kinds of oscillations, i.e., upward propagation of the waves, () fast waves with a phase speed of km s in the horizontal displacements (kink mode) and a propagating speed of km s in intensity oscillations (longitudinal waves), and () a relatively wide range of phase spectra, which may describe different dispersion relations and/or belong to different sources [?]. In addition, there is a slight indication of standing and downward propagating waves in the intensity oscillations.

Fast waves are of interest because they can propagate to the upper solar atmosphere and carry energy. By comparing our results with those from theoretical investigations (see Section 5), observations of the fast waves could possibly be explained as a result of mode conversion and/or mode transmission at the level (i.e., where Alfvén and sound speeds coincide; [?]). The observed high-frequency waves in the magnetic elements could be excited due to, e.g., buffeting of flux-tube footpoints by high-frequency perturbations caused by the surrounding granules and by intergranular turbulence [?].

Summarizing, we have observed high-frequency, fast, upward propagating waves using data from Sunrise unaffected by seeing. Fast waves in both, horizontal displacement and intensity oscillations appear to dominate over slow waves in our data. We speculate that the rather large propagation speeds that we deduced from our data (particularly those obtained from smaller concentrations in the phase diagram of Figure 7) could also be due to (1) refraction of the propagation path above the magnetic canopy for longitudinal waves [?], and (2) to the superposition of several wave trains, e.g., polarized kink waves for transverse waves [?].

We cannot completely rule out that we have overestimated the formation height of the Ca ii H passband, which may have led to overly large phase speeds. This is because the plane-parallel atmospheres used to compute the contribution functions may not be representative of the complex 3D solar atmosphere. Consequently, we determined the height difference between the two wavelength bands by two entirely independent means. Besides the contribution functions, we also used the time difference between the arrival of acoustic waves in the quiet-Sun at the two heights. Both methods gave very similar results in the quiet-Sun, which greatly increases our confidence in our conclusions. In order to find out more about the actual cause of these fast waves, a similar study as was done here, but using synthetic passbands from 3D MHD simulation is essential. Further studies from, e.g., simultaneous observations of multiple atmospheric layers (from the photosphere to the transitions-region/corona) might clarify to what extent the high-frequency fast waves can reach the upper solar atmosphere and how much energy they release at those heights.

The German contribution to Sunrise is funded by the Bundesministerium für Wirtschaft und Technologie through the Deutsches Zentrum für Luft- und Raumfahrt e.V. (DLR), grant No. 50 OU 0401, and by the Innovationsfond of the President of the Max Planck Society (MPG). The Spanish contribution has been funded by the Spanish MICINN under projects ESP2006-13030-C06 and AYA2009-14105-C06 (including European FEDER funds). The 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. S.J. receives support from the Research Council of Norway.