Free Access
Issue
A&A
Volume 649, May 2021
Article Number A169
Number of page(s) 7
Section The Sun and the Heliosphere
DOI https://doi.org/10.1051/0004-6361/202140429
Published online 02 June 2021

© ESO 2021

1. Introduction

It is well known that large-scale solar magnetic features, such as pores and sunspots, are dominated by ∼3 mHz (∼5 min) oscillations in the photosphere and ∼5 mHz (∼3 min) oscillations in the chromosphere (Centeno et al. 2009, 2006; Nagashima et al. 2007; Khomenko et al. 2008; Felipe et al. 2010; Stangalini et al. 2011; Jess et al. 2012, 2013, 2015; Khomenko & Collados 2015, to mention but a few). The generally accepted view is that p-modes are progressively absorbed by the magnetic field of sunspots and pores (see for instance Spruit & Bogdan 1992; Spruit 1996) and converted to magneto-acoustic modes at the equipartition layer (Cally et al. 1994, 2016; Cally & Goossens 2008), which is the layer where the Alfvén and sound speeds are approximately equal. These magneto-acoustic waves can propagate upwards to the chromosphere where, due to the effect of an atmospheric cutoff frequency for the magneto-acoustic waves (Jefferies et al. 2006; Stangalini et al. 2011; Felipe 2019), the dominant period becomes ∼5 mHz. Interestingly, Schroeter & Soltau (1976) reported that the locations with high oscillatory power at ∼5 mHz were uncorrelated with those with high oscillatory power at ∼3 mHz, which Lites (1986) suggested may be a consequence of different underlying driving mechanisms. An alternative mechanism revolves around the creation of a resonance cavity in the solar atmosphere, which has recently been shown to exist by Jess et al. (2020) and independently confirmed by Felipe et al. (2020). Other localised waves can also be excited due to residual convection in the umbra (Zhugzhda & Sych 2018).

Spruit (1982), Roberts (1983), and Edwin & Roberts (1983) predicted that magnetic flux tubes, in addition to localised waves, can also support the excitation of the magneto-hydrodynamic (MHD) global resonant modes of the structure. More recently, Roberts (2019) argued that the observed velocity field in magnetic structures should therefore be regarded as the superposition of both global and local disturbances, giving rise to the complex oscillatory patterns generally seen in high resolution imaging of the solar atmosphere. This implies that filtering techniques are required to correctly identify the eigenmodes (e.g., Jess et al. 2017; Albidah et al. 2021).

However, as of now, only low-order resonant modes have been identified, mostly in small-scale magnetic structures. These include the sausage modes (Morton et al. 2011; Martínez González et al. 2011; Gafeira et al. 2017; Keys et al. 2018; Kang et al. 2019; Guevara Gómez et al. 2021), kink modes (Keys et al. 2011, 2020; Stangalini et al. 2013, 2014, 2015; Jafarzadeh et al. 2017a,b; Morton et al. 2021), and torsional Alfvén modes (Jess et al. 2009). The majority of these works were based on the analysis of intensity or velocity fluctuations (both line-of-sight [LoS] and horizontal motions). Only a few have studied the magnetic perturbations associated with, for example, area oscillations (Martínez González et al. 2011; Grant et al. 2015; Freij et al. 2016) or other spectropolarimetric diagnostics (Gilchrist-Millar et al. 2021; Keys et al. 2021) that have recently opened up the possibility of detecting true magnetic oscillations (Stangalini et al. 2018), which, in turn, would help identify clear links between such wave behaviour and the chemical abundance in the solar corona (Stangalini et al. 2021; Baker et al. 2021). However, as noted above, the velocity fields inside the umbra of pores and sunspots are generally dominated by ∼3 mHz oscillations at photospheric heights; only one exception has so far been reported to our knowledge (Stangalini et al. 2012), in which the dominance of 3 min velocity oscillations at photospheric heights was reported.

Here, we further investigate the peculiar dominance of 3 min oscillations reported by Stangalini et al. (2012) in the umbra of a magnetic pore at photospheric heights. In particular, by fully exploiting a combination of velocity and spectropolarimetric information, we are able to identify the transition of p-modes to a series of spectral features within the umbra that we interpret as the co-existence of multiple MHD resonant modes in the associated flux tubes.

The dataset (or parts thereof) we study here was previously used by Stangalini et al. (2011, 2012), by Sobotka et al. (2013, 2016), and, more recently, by Abbasvand et al. (2020). The Sobotka et al. (2016) study focused on chromospheric oscillations measured in the Ca II 854.2 nm line; the authors estimated the acoustic energy flux in the chromosphere and compared it to radiative losses in that region. They also reported variations of the power spectral density of intensity and velocity oscillations at chromospheric heights as a function of the distance from the centre of the umbra. Here we study spectropolarimetric data obtained in the photospheric Fe I 617.3 nm line, and we extend previous analyses that mostly considered Doppler velocity oscillations by including oscillations measured in circular polarisation (CP; i.e., the magnetic field). By simultaneously investigating the dependence of the power spectra of both Doppler velocity and CP on the magnetic field strength, we are able to identify global resonances of the magnetic structure. This allows us to detect, for the first time, magnetic field oscillations associated with these resonances. Compared to the power-distance plot introduced by Sobotka et al. (2013), the new B − ω diagrams presented here offer the additional benefit of being applicable to the study of magnetic structures with irregular shapes. Furthermore, being based on both velocity and magnetically derived diagnostics, they fully exploit all the information that can be extracted from 2D spectropolarimetric data.

2. Dataset and methods

The data used in this work (see Fig. 1) were acquired on 2008 October 15 at 16:30 UT with the Interferometric BIdimensional Spectrometer (IBIS; Cavallini 2006) at the Dunn Solar Telescope (New Mexico, USA). The region observed was AR11005, which appears as a small pore with a light bridge at [25.2 N, 10.0 W]. The dataset consists of 80 spectral scans in full Stokes mode, each containing 21 spectral points of the Fe I 617.3 nm line. The δλ between two consecutive spectral points was 2 pm, and the exposure time was set to 80 ms. The cadence of the data is 52 s, and the pixel scale is 0.167 arcsec. Simultaneous white light and G-band images were restored with the multi-frame blind deconvolution (MFBD; van Noort et al. 2005). Spectropolarimetric images were then co-registered and de-stretched to minimise the residual seeing aberrations uncorrected by the adaptive optics system. The calibration pipeline includes dark frame subtraction, flat fielding, and polarimetric demodulation, and it also corrects for blueshift effects (Reardon & Cavallini 2008).

thumbnail Fig. 1.

Magnetic pore investigated in the study. Top left: instantaneous Fe I 617.3 nm line core intensity. The dotted red line represents the B = 800 G contour used as a reference in Fig. 2. Top right: magnitude of CP. Bottom left: instantaneous Doppler velocity map obtained in the Fe I 617.3 nm spectral line. Bottom right: average magnetic field strength obtained with the centre-of-gravity (COG) method. The white dashed box outlines the region considered when filtering in k − ω space (see text and Fig. 3).

The Doppler velocity was estimated using the method based on the estimate of the phase of the first Fourier component, as in Schlichenmaier & Schmidt (2000). The CP is defined as the average of the absolute values of the Stokes-V profile. In this regard, only the spectral points closest to the maxima (ten points) of the lobes of the profile are considered. This is done to increase the signal-to-noise ratio and exclude spectral points close to the continuum, which do not contribute much to the polarisation signal but contribute largely to the noise.

3. Results

Here we investigate the spatial dependence of the dynamics on the magnetic flux. For this purpose we constructed a specific diagram where each column represents the average power spectrum of pixels contained within a specific magnetic field range, which was set to an interval of 80 G in our case. This plot, which we call a B − ω diagram (where B is the magnetic field strength and ω represents the oscillation frequency), is similar to the power-distance diagram investigated in Sobotka et al. (2013); it shows the modification of the power spectrum as one moves from the quiet Sun towards the inner region of the umbra (i.e., as one moves from smaller to larger values of B). However, the B − ω methodology can easily be applied to magnetic structures with irregular shapes, where the centre of the structure itself cannot be easily and unambiguously identified. This type of visualisation tool is similar to existing k − ω diagrams (where k represents the spatial wavenumber of the observations; e.g., Duvall et al. 1988; Krijger et al. 2001; Rutten & Krijger 2003; Kneer & Bello González 2011; Jess et al. 2012), only now with the emphasis placed on the magnitude of the embodied magnetic fields rather than the spatial extent of their composition, thus fully exploiting the spectropolarimetric information.

In Fig. 2 we show the B − ω diagram for both the LoS Doppler velocity (left) and the CP (middle). Here, the CP fluctuations are a measure of the perturbations of the LoS magnetic flux. In correspondence with the umbral boundary inferred from continuum intensity maps (see Fig. 1), we deduced an average magnetic field strength of ≈800 G (see also the upper-left panel of Fig. 1) that segregates the umbra from the surrounding quiet Sun, which is highlighted in the left and middle panels of Fig. 2 with vertical dashed blue lines. From Fig. 2, we observe a transition from the outside 3 mHz oscillations, corresponding to p-modes, to a series of spectral features inside the umbra (i.e., for B > 800 G). It is worth noting here that none of the umbral spectral features observed in LoS velocity share the same frequency of the p-modes observed outside the umbra (i.e., ∼3 mHz), contrary to what one would expect.

thumbnail Fig. 2.

B − ω diagram of the LoS velocity (left) and the CP (centre). The vertical dashed blue line represents the approximate position of the boundary of the umbra as inferred from intensity images. Each column represents the average power spectrum across bins equal to 80 G. The global spectra for LoS velocities and CP fluctuations, both outside and inside the magnetic structure, are shown in the right panel. These are obtained by integrating the B − ω diagram along the horizontal axis.

In the right panel of Fig. 2 we plot the average spectra inside and outside the umbra to facilitate visualisation. Here, it is clear that frequencies corresponding to the p-modes (i.e., ∼3 mHz) are observed as a reduction of power in both LoS velocity and CP signal in the umbral locations. In turn, the umbra itself is dominated by a series of peaks at approximately 2 mHz, 4–6 mHz, and > 8 mHz, which are each split into multiple subcomponents. The amplitudes of these oscillations are even larger than those of the p-modes outside the magnetic structure. It is also worth noting that the amplitude of the oscillations in the 4–6 mHz band is even larger than that of the surrounding p-modes.

Regarding the left and middle panels of Fig. 2, it is interesting to note that these panels show different patterns of spectral features. In particular, there are no evident spectral features in the B − ω (CP) diagram corresponding to the most prominent peaks in the B − ω (LoS velocity) diagram. For example, with regard to the most prominent peak in the B − ω (CP) diagram at ≈4 mHz, we note that there is no equivalent increase in power in the B − ω (LoS velocity) signal. This can be seen more clearly in the right panel of Fig. 2, where we compare the average LoS velocity and CP spectra inside the umbra. The main spectral features present in the Doppler velocity spectrum correspond to a lack of power for CP oscillations (e.g., in the 4.5–6.0 mHz band). In other words, there is no strict one-to-one correspondence between the LoS Doppler velocity and CP umbral spectra. However, this does not mean that there is absolutely no power in CP corresponding to the main spectral features in LoS velocity or vice versa. Indeed, there exist a few spectral features dominated by CP oscillations, for which a velocity signal is also detected (and vice versa), although with small amplitude.

Following the approach documented by Jess et al. (2017), we produced k − ω diagrams of LoS velocity and CP, which are shown in Fig. 3. Similar plots for this specific magnetic pore were also investigated by Sobotka et al. (2013); however, here we focus on a smaller, square field-of-view (shown in the lower-right panel of Fig. 1) to limit contamination from the surrounding quiet Sun. This helps to isolate the oscillatory signals in the umbra and reliably identify a series of horizontal power enhancements (panels a and b) in Fig. 3 – corresponding to the main spectral features already seen in the left and middle panels of Fig. 2. In the same figure we also show the root-mean-square (RMS) velocity and CP amplitude (panels c and d) and the instantaneous LoS Doppler velocities and CP signals after filtering in the k − ω band pass highlighted by the dashed contours in Fig. 3a. Here, the filter has a Gaussian profile to limit the edge effects synonymous with discontinuities in Fourier space. These maps show the presence of radial high-order oscillations, with the oscillatory pattern largely affected by the cross-sectional shape of the magnetic structure itself. Here, we also note that, surprisingly, the CP and velocity oscillations have different spatial distributions: The CP oscillations are more concentrated in the inner part of the umbra, while the velocity oscillations are dominant within an annulus close to the umbral perimeter. This can also be seen in the power maps at different frequencies that are shown in Fig. 4.

thumbnail Fig. 3.

k − ω diagrams of the LoS velocity (panel a) and CP (panel b). The dashed white line represents the acoustic theoretical fundamental mode (gk, where g is the gravitational acceleration and k the horizontal wavenumber), while the continuous black line corresponds to the Lamb line and shows the theoretically expected propagation at the sound speed (csk, where cs is the sound speed). Panels c and d: RMS amplitude of the filtered velocity and CP, respectively. Panels e and f: filtered velocities and CP signals, respectively, at one instant in time. The filter widths (in both k and ω space) are illustrated by the dashed black boxes in panels a and b, where most of the power is located.

thumbnail Fig. 4.

Power maps of vLoS (top) and CP (bottom) in different frequency bands (1 mHz width).

The intensity images of the magnetic structure indicate the presence of a light bridge dividing the structure itself into two lobes. However, despite this, the pore behaves as a single coherent structure, at least from the perspective of the wave dynamics. This fact can also be seen in Fig. 5, where we plot the time-distance diagrams of the unfiltered CP and Doppler velocity fluctuations in a slice that passes through the pore. The two lobes oscillate in phase, both for CP and Doppler velocity perturbations. Here, in order to highlight the perturbations with respect to the background, we removed the temporal average before plotting each quantity.

thumbnail Fig. 5.

Time-distance plots for the unfiltered CP and Doppler velocity fluctuations (δCP = CP − ⟨CP⟩ and δvLoS = vLoS − ⟨vLoS⟩) in a slice that passes through the entire magnetic structure (i.e., y = 6 arcsec with reference to the left panels of Fig. 3). The vertical dashed lines represent the approximate position of the light bridge.

4. Discussion

Sunspots and pores are generally dominated by 3 mHz velocity oscillations in the solar photosphere, which are normally interpreted as the result of p-mode penetration into the magnetic structure and further conversion to magneto-acoustic modes (e.g., Khomenko & Collados 2015). While, it was predicted that magnetic structures could also be affected by global resonant modes, these have only been detected in small-scale magnetic features or in a sunspot after thorough filtering (Jess et al. 2015). Roberts (2019) pointed out that the complex oscillatory behaviour inside the umbra of sunspots and pores should be regarded as the superposition of both global eigenmodes of the magnetic structure and other locally excited magneto-acoustic waves and disturbances. However, the lack of dominant frequencies other than 3 mHz in the umbra of sunspots or pores suggests that the contribution of global eigenmodes to the overall velocity field is smaller than that due to the local disturbances.

The dataset explored in this work provides an ideal opportunity to study the possible presence of global resonances within a magnetic structure. Indeed, in this case, the umbra is not dominated by 3 mHz LoS Doppler velocity fluctuations, as one would expect in the photosphere as a result of p-mode conversion, but by 5 mHz oscillations that are generally only found in the chromosphere as a consequence of atmospheric stratification and the effect of the acoustic cutoff frequency (Felipe 2019).

By combining spectropolarimetric information and Doppler velocities, we were able to characterise the wave dynamics inside the umbra and identify a series of spectral features that can be interpreted as global eigenmodes of the magnetic flux tube. These can easily be identified in a novel B − ω diagram as horizontal spectral features in LoS Doppler velocity, which arise immediately inside the perimeter boundary of the magnetic structure. This interpretation is further supported by the fact that there are no spectral features in LoS velocity at the same frequencies corresponding to p-modes, thus allowing, for the first time, a clear distinction between p-mode absorption and resonant modes. Furthermore, a B − ω diagram corresponding to CP signals allowed for the detection of magnetic oscillations within the umbral boundary, which are expected from MHD theory (Spruit 1982; Roberts 1983; Edwin & Roberts 1983). In contrast to previous approaches based on the estimation of power as a function of distance from the centre of the umbra (e.g., Sobotka et al. 2013), our approach can even be readily applied to magnetic tubes with very irregular shapes, for which the centre of the structure itself cannot be easily defined. In addition, our approach does not mix the signals from magnetic and non-magnetic pixels, thus being intrinsically insensitive to contamination from ‘quiet Sun’ pixels.

The dominant peaks in the B − ω diagram of the velocity and CP oscillations are not precisely at the same frequency. While this rules out the possibility of opacity effects and cross talk, which would show up at the same frequency by definition, this is a rather surprising result for which we cannot offer a complete explanation yet, merely some speculative ideas. One such idea is based on the detection of mixed fast- and slow-mode waves. To illustrate this here we used the standard model composed of a straight cylindrical magnetic flux tube with a purely axial magnetic field (e.g., Spruit 1981; Edwin & Roberts 1983). We assumed a uniform magnetic field, B0, pointing in the direction of the observer, which is aligned with the z-direction for convenience. Under this situation, velocity perturbations along the magnetic field, vz, correspond to velocities along the LoS, while magnetic perturbations along the equilibrium magnetic field, bz, produce changes in CP. Using the linearised ideal MHD equations, the amplitude ratio of the two perturbations for propagating magneto-acoustic linear waves can be written as

| b ¯ z | | v ¯ z | = ω s ω | ω 2 ω s 2 1 | , $$ \begin{aligned} \frac{\left|\bar{b}_z\right|}{\left|\bar{{ v}}_z\right|}=\frac{\omega _{\rm s}}{\omega }\left|\frac{\omega ^2}{\omega ^2_{\rm s}}-1\right| , \end{aligned} $$(1)

where b ¯ z = b z / B 0 $ \bar{b}_z=b_z/B_0 $ is the dimensionless magnetic field perturbation and v ¯ z = v z / c s $ \bar{v}_z=v_z/c_s $ is the dimensionless velocity fluctuation, with cs the equilibrium sound speed inside the magnetic flux tube.

In Eq. (1), ω is the frequency of the propagating wave and ωs = kzcs is the slow (acoustic) frequency for a wave with a longitudinal wavelength equal to 2π/kz. Under photospheric conditions, if the propagating wave is a slow MHD mode, then its frequency is typically slightly below the acoustic frequency (ω ≲ ωs). Therefore, from Eq. (1) we have | b ¯ z | / | v ¯ z | 1 $ \left|\bar{b}_z\right|/\left|\bar{\mathit{v}}_z\right|\ll 1 $. Conversely, for a fast MHD wave the frequency is always above the acoustic frequency (ω > ωs), meaning that the situation | b ¯ z | / | v ¯ z | > 1 $ \left|\bar{b}_z\right|/\left|\bar{\mathit{v}}_z\right| > 1 $ is most likely to occur. We note that in Eq. (1) it is assumed that the waves have a propagating nature and that the LoS is parallel to the flux tube axis. General expressions for different orientations (including standing waves) can be found in Moreels & Van Doorsselaere (2013). In other words, under these assumptions, slow and fast MHD modes should be dominated by velocity or magnetic perturbations, respectively. While a combination of fast and slow waves is an interesting possibility to explain the presence of both Doppler velocity and magnetically dominated signals, we note that the CP B − ω diagram is dominated by a 4 mHz peak, which is below the acoustic cutoff. Therefore, the above scenario, while providing an explanation for the high frequency power in CP, does not completely explain the dominant 4 mHz peak. However, we should keep in mind that CP and velocity signals are extracted from different heights in the solar atmosphere. Indeed, the polarisation signals originate from heights representative of the wings of the spectral line, while spectral points close to the spectral line core may contribute predominantly to the velocity signatures.

To illustrate this effect, in Fig. 6 we show the power spectra of the velocity oscillations computed with the bisector method by considering different spectral positions, thus spanning different geometric heights in the solar atmosphere, from the line core forming at approximately 250 km (Fleck et al. 2011) down to the base of the photosphere, corresponding to the continuum level. It is worth noting here that the line core formation height of 250 km is representative of the quiet Sun. In magnetic structures, we would expect all the relevant formation heights to be shifted down, with the height difference between velocity and CP remaining roughly the same. Although the range of heights explored by scanning the line with the bisector method is small, the power spectra of the Doppler velocity oscillations in the magnetic umbra show large variations with geometric height, an effect not observed in the corresponding quiet Sun. In particular, heights close to the continuum show a 3 mHz peak that is gradually suppressed as one moves upwards to the heights associated with the line core; however, the spectrum is still more dominated by the peaks in the 4–6 mHz band, alongside other higher frequency features. While these higher frequencies are also observed at the base of the photosphere (i.e., near the continuum level), they become the dominant sources of spectral power as one moves towards the spectral line core. It is worth stressing that in the quiet Sun there is essentially no difference in spectral power across the range of geometric heights. In addition, we also note that at intermediate heights, similar to those where the CP signal originates (i.e., the wings of the line), a 4 mHz peak is also observed. These results suggest that the method employed to measure the Doppler shifts used in this work is more sensitive to the velocities at atmospheric heights corresponding to the spectral line core. However, although at geometric heights similar to those of CP there exists a small 4 mHz peak, the change in the formation height of the signals does not explain why in CP there is a lack of power in the 4–6 mHz band, the dominant frequency band of the Doppler velocities.

thumbnail Fig. 6.

Variation of the mean periodogram of the Doppler velocity, obtained from the bisector method from the core to the continuum, for the quiet Sun (pixels below 20 G) and the magnetic pore (pixels above 700 G)

In summary, although in our opinion these two options offer some interesting discussion points, neither of them alone or in combination can completely explain the dynamic behaviour of this magnetic pore. Although we are not able to solve this puzzle at the moment, we believe it is of paramount importance to report this surprising behaviour in the hope that state-of-the-art numerical modelling or theoretical advances will help to solve this riddle in the near future. Furthermore, we stress that the different spatial distribution of the CP and velocity power maps (Fig. 4), while confirming once again the presence of a true magnetic oscillation, suggests the co-existence of two different MHD eigenmodes.

The B − ω diagrams depicted in Fig. 2 capture the transition between p-modes outside the umbra to global resonant modes inside it. This study therefore provides a unique exception where the two components of wave dynamics can be readily disentangled. However, a question remains as to whether this effect is visible due to the nearly complete suppression of the p-modes inside the umbra or if it is a result of the favourable contrast between the spectral features at 3 mHz (with respect to the others) that allows the detection of the intrinsic resonant modes of the structure. The latter possibility appears to be the case here as the amplitude of the spectral features observed inside the umbra is even larger than that of the external p-modes, hence providing significant contrast for the 3 mHz oscillations.

In our opinion, the results presented in the current study may have something to do with the excitation mechanisms of the waves and the energy associated with the external driver (e.g., turbulent convection). This suggests that the detected spectral features are associated with distinct eigenmodes, some of which show a prevalence to LoS Doppler velocity oscillations, with others associated with magnetic fluctuations. Nevertheless, the manifestation of magnetic oscillations suggests the presence of real MHD modes. If these power structures are validated as resonances, this will open up several new possibilities in magneto-helioseismology and might allow measuring magnetic fields in regions of the upper atmosphere that are currently difficult to probe.

5. Conclusions

Sunspots and pores are generally found to be dominated by 3 mHz oscillations at photospheric heights. The consensus is that this is the result of p-mode penetration and absorption at the same frequency. However, it was theoretically predicted that flux tubes should also support the excitation of MHD resonant modes, which are intrinsic global oscillations of the structure. Theory also predicts that these should have an associated Doppler velocity and magnetic signal. A unique magnetic pore observed by IBIS at high resolution shows a dominance of 5 mHz oscillations in the photosphere, instead of the typically expected 3 mHz signals. By combining Doppler velocity and magnetic information obtained from polarisation measurements, it has been possible to distinguish the transition from ambient p-modes to internal global resonances of the magnetic structure. These are seen as a series of spectral features that arise immediately inside the umbra of the magnetic structure. Some of them are also associated with magnetic oscillations, supporting the interpretation in terms of resonant MHD modes. A novel diagnostic, in the form of a B − ω diagram that combines LoS Doppler velocities and magnetic field information at the same time, clearly captures this transition, showing the progressive absorption of p-modes as one gets closer to the magnetic structure, followed by the onset of resonant modes inside the umbra.

Acknowledgments

The authors wish to acknowledge scientific discussions with the Waves in the Lower Solar Atmosphere (WaLSA; www.WaLSA.team) team, which is supported by the Research Council of Norway (project number 262622), and The Royal Society through the award of funding to host the Theo Murphy Discussion Meeting “High-resolution wave dynamics in the lower solar atmosphere” (grant Hooke18b/SCTM). DBJ is grateful to Invest NI and Randox Laboratories Ltd. for the award of a Research & Development Grant (059RDEN-1), in addition to the UK Science and Technology Facilities Council (STFC) grant ST/T00021X/1. VF, GV thank to The Royal Society, International Exchanges Scheme, collaboration with Chile (IE/170301) and Brazil (IES/R1/191114). VF and GV are grateful to the Science and Technology Facilities Council (STFC) grant ST/V000977/1 for support provided. SJ acknowledges support from the European Research Council under the European Union Horizon 2020 research and innovation program (grant agreement No. 682462) and from the Research Council of Norway through its Centres of Excellence scheme (project No. 262622). CDM would like to thank the Northern Ireland Department for the Economy for the award of a PhD studentship. DB is funded under STFC consolidated grant No. ST/S000240/1. SMJ acknowledges support under award 1829258 from the National Science Foundation. AA acknowledges the Deanship of Scientific Research (DSR), King Faisal University, Al-Hassa, KSA for the financial support under Nasher Track (grant No.186354). This research has received funding from the European Union’s Horizon 2020 Research and Innovation program under grant agreement No 82135 and 814335 (SOLARNET) and No 739500 (PRE-EST). JT and RS acknowledge the support from grant AYA2017-85465-P (MINECO/AEI/FEDER, UE). This research has made use of the IBIS-A archive.

References

  1. Abbasvand, V., Sobotka, M., Heinzel, P., et al. 2020, ApJ, 890, 22 [Google Scholar]
  2. Albidah, A. B., Brevis, W., Fedun, V., et al. 2021, Phil. Trans. R. Soc. A, 379, 20200181 [Google Scholar]
  3. Baker, D., Stangalini, M., Valori, G., et al. 2021, ApJ, 907, 16 [Google Scholar]
  4. Cally, P. S., & Goossens, M. 2008, Sol. Phys., 251, 251 [Google Scholar]
  5. Cally, P. S., Bogdan, T. J., & Zweibel, E. G. 1994, ApJ, 437, 505 [Google Scholar]
  6. Cally, P. S., Moradi, H., & Rajaguru, S. P. 2016, GMS, 216, 489 [Google Scholar]
  7. Cavallini, F. 2006, Sol. Phys., 236, 415 [NASA ADS] [CrossRef] [Google Scholar]
  8. Centeno, R., Collados, M., & Trujillo Bueno, J. 2006, ApJ, 640, 1153 [Google Scholar]
  9. Centeno, R., Collados, M., & Trujillo Bueno, J. 2009, ApJ, 692, 1211 [Google Scholar]
  10. Duvall, T. L., Jr., Harvey, J. W., Libbrecht, K. G., et al. 1988, ApJ, 324, 1158 [Google Scholar]
  11. Edwin, P. M., & Roberts, B. 1983, Sol. Phys., 88, 179 [Google Scholar]
  12. Felipe, T. 2019, A&A, 627, A169 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Felipe, T., Khomenko, E., Collados, M., et al. 2010, ApJ, 722, 131 [Google Scholar]
  14. Felipe, T., Kuckein, C., Manrique, S., et al. 2020, ApJ, 900, L29 [Google Scholar]
  15. Fleck, B., Couvidat, S., & Straus, T. 2011, Sol. Phys., 271, 27 [Google Scholar]
  16. Freij, N., Dorotovič, I., Morton, R. J., et al. 2016, ApJ, 817, 44 [Google Scholar]
  17. Gafeira, R., Jafarzadeh, S., Solanki, S. K., et al. 2017, ApJS, 229, 7 [Google Scholar]
  18. Gilchrist-Millar, C. A., Jess, D. B., Grant, S. D. T., et al. 2021, Phil. Trans. R. Soc. A, 379, 20200172 [Google Scholar]
  19. Grant, S. D. T., Jess, D. B., Moreels, M. G., et al. 2015, ApJ, 806, 132 [Google Scholar]
  20. Guevara Gómez, J. C., Jafarzadeh, S., Wedemeyer, S., et al. 2021, Phil. Trans. R. Soc. A, 379, 20200184 [Google Scholar]
  21. Jafarzadeh, S., Solanki, S. K., Gafeira, R., et al. 2017a, ApJS, 229, 9 [NASA ADS] [CrossRef] [Google Scholar]
  22. Jafarzadeh, S., Solanki, S. K., Stangalini, M., et al. 2017b, ApJS, 229, 10 [Google Scholar]
  23. Jefferies, S. M., McIntosh, S. W., Armstrong, J. D., et al. 2006, ApJ, 648, L151 [Google Scholar]
  24. Jess, D. B., Mathioudakis, M., Erdélyi, R., et al. 2009, Science, 323, 1582 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  25. Jess, D. B., Shelyag, S., Mathioudakis, M., et al. 2012, ApJ, 746, 183 [Google Scholar]
  26. Jess, D. B., Reznikova, V. E., Van Doorsselaere, T., et al. 2013, ApJ, 779, 168 [Google Scholar]
  27. Jess, D. B., Morton, R. J., Verth, G., et al. 2015, Space Sci. Rev., 190, 103 [Google Scholar]
  28. Jess, D. B., Van Doorsselaere, T., Verth, G., et al. 2017, ApJ, 842, 59 [Google Scholar]
  29. Jess, D. B., Snow, B., Houston, S. J., et al. 2020, Nat. Astron., 4, 220 [Google Scholar]
  30. Kang, J., Chae, J., Nakariakov, V. M., et al. 2019, ApJ, 877, 1 [Google Scholar]
  31. Keys, P. H., Mathioudakis, M., Jess, D. B., et al. 2011, ApJ, 740, 40 [Google Scholar]
  32. Keys, P. H., Morton, R. J., Jess, D. B., et al. 2018, ApJ, 857, 28 [Google Scholar]
  33. Keys, P. H., Reid, A., Mathioudakis, M., et al. 2020, A&A, 633, 60 [Google Scholar]
  34. Keys, P. H., Steiner, O., & Vigeesh, G. 2021, Phil. Trans. R. Soc. A, 379, 20200182 [Google Scholar]
  35. Khomenko, E., & Collados, M. 2015, Liv. Rev. Sol. Phys., 12, 6 [Google Scholar]
  36. Khomenko, E., Centeno, R., Collados, M., et al. 2008, ApJ, 676, L85 [Google Scholar]
  37. Kneer, F., & Bello González, N. 2011, A&A, 532, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Krijger, J. M., Rutten, R. J., Lites, B. W., et al. 2001, A&A, 379, 1052 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Lites, B. W. 1986, ApJ, 301, 992 [Google Scholar]
  40. Martínez González, M. J., Asensio Ramos, A., Manso Sainz, R., et al. 2011, ApJ, 730, L37 [NASA ADS] [CrossRef] [Google Scholar]
  41. Moreels, M. G., & Van Doorsselaere, T. 2013, A&A, 551, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Morton, R. J., Erdélyi, R., Jess, D. B., et al. 2011, ApJ, 729, L18 [Google Scholar]
  43. Morton, R. J., Mooroogen, K., & Henriques, V. M. J. 2021, Phil. Trans. R. Soc. A, 379, 20200183 [Google Scholar]
  44. Nagashima, K., Sekii, T., Kosovichev, A. G., et al. 2007, PASJ, 59, S631 [Google Scholar]
  45. Reardon, K. P., & Cavallini, F. 2008, A&A, 481, 897 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Roberts, B. 1983, Sol. Phys., 87, 77 [Google Scholar]
  47. Roberts, B. 2019, MHD Waves in the Solar Atmosphere (Cambridge University Press) [Google Scholar]
  48. Rutten, R. J., & Krijger, J. M. 2003, A&A, 407, 735 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Schlichenmaier, R., & Schmidt, W. 2000, A&A, 358, 1122 [NASA ADS] [Google Scholar]
  50. Schroeter, E. H., & Soltau, D. 1976, A&A, 49, 463 [Google Scholar]
  51. Sobotka, M., Švanda, M., Jurčák, J., et al. 2013, A&A, 560, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Sobotka, M., Heinzel, P., Švanda, M., et al. 2016, ApJ, 826, 49 [NASA ADS] [CrossRef] [Google Scholar]
  53. Spruit, H. C. 1981, A&A, 98, 155 [NASA ADS] [Google Scholar]
  54. Spruit, H. C. 1982, Sol. Phys., 75, 3 [Google Scholar]
  55. Spruit, H. C. 1996, Bull. Astron. Soc. India, 24, 211 [Google Scholar]
  56. Spruit, H. C., & Bogdan, T. J. 1992, ApJ, 391, L109 [Google Scholar]
  57. Stangalini, M., Del Moro, D., Berrilli, F., et al. 2011, A&A, 1, 1 [Google Scholar]
  58. Stangalini, M., Giannattasio, F., Del Moro, D., et al. 2012, A&A, 539, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. Stangalini, M., Solanki, S. K., Cameron, R., et al. 2013, A&A, 554, 115 [Google Scholar]
  60. Stangalini, M., Consolini, G., Berrilli, F., et al. 2014, A&A, 569, A102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Stangalini, M., Giannattasio, F., & Jafarzadeh, S. 2015, A&A, 577, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Stangalini, M., Jafarzadeh, S., Ermolli, I., et al. 2018, ApJ, 869, 110 [Google Scholar]
  63. Stangalini, M., Baker, D., Valori, G., et al. 2021, Phil. Trans. R. Soc. A, 379, 20200216 [Google Scholar]
  64. van Noort, M., Rouppe van der Voort, L., & Läfdahl, M. G. 2005, Sol. Phys., 228, 191 [Google Scholar]
  65. Zhugzhda, Y., & Sych, R. 2018, Res. Astron. Astrophys., 18, 105 [Google Scholar]

All Figures

thumbnail Fig. 1.

Magnetic pore investigated in the study. Top left: instantaneous Fe I 617.3 nm line core intensity. The dotted red line represents the B = 800 G contour used as a reference in Fig. 2. Top right: magnitude of CP. Bottom left: instantaneous Doppler velocity map obtained in the Fe I 617.3 nm spectral line. Bottom right: average magnetic field strength obtained with the centre-of-gravity (COG) method. The white dashed box outlines the region considered when filtering in k − ω space (see text and Fig. 3).

In the text
thumbnail Fig. 2.

B − ω diagram of the LoS velocity (left) and the CP (centre). The vertical dashed blue line represents the approximate position of the boundary of the umbra as inferred from intensity images. Each column represents the average power spectrum across bins equal to 80 G. The global spectra for LoS velocities and CP fluctuations, both outside and inside the magnetic structure, are shown in the right panel. These are obtained by integrating the B − ω diagram along the horizontal axis.

In the text
thumbnail Fig. 3.

k − ω diagrams of the LoS velocity (panel a) and CP (panel b). The dashed white line represents the acoustic theoretical fundamental mode (gk, where g is the gravitational acceleration and k the horizontal wavenumber), while the continuous black line corresponds to the Lamb line and shows the theoretically expected propagation at the sound speed (csk, where cs is the sound speed). Panels c and d: RMS amplitude of the filtered velocity and CP, respectively. Panels e and f: filtered velocities and CP signals, respectively, at one instant in time. The filter widths (in both k and ω space) are illustrated by the dashed black boxes in panels a and b, where most of the power is located.

In the text
thumbnail Fig. 4.

Power maps of vLoS (top) and CP (bottom) in different frequency bands (1 mHz width).

In the text
thumbnail Fig. 5.

Time-distance plots for the unfiltered CP and Doppler velocity fluctuations (δCP = CP − ⟨CP⟩ and δvLoS = vLoS − ⟨vLoS⟩) in a slice that passes through the entire magnetic structure (i.e., y = 6 arcsec with reference to the left panels of Fig. 3). The vertical dashed lines represent the approximate position of the light bridge.

In the text
thumbnail Fig. 6.

Variation of the mean periodogram of the Doppler velocity, obtained from the bisector method from the core to the continuum, for the quiet Sun (pixels below 20 G) and the magnetic pore (pixels above 700 G)

In the text

Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.