In stratified atmospheres, acoustic waves can only propagate if their frequency is above the cutoff value. Different theories provide different cutoff values. We developed an alternative method to derive the cutoff frequency in several standard solar models, including various quiet-Sun and umbral atmospheres. We performed numerical simulations of wave propagation in the solar atmosphere. The cutoff frequency is determined from the inspection of phase difference spectra computed between the velocity signal at two atmospheric heights. The process is performed by choosing pairs of heights across all the layers between the photosphere and the chromosphere, to derive the vertical stratification of the cutoff in the solar models. The cutoff frequency predicted by the theoretical calculations departs significantly from our measurements. In quiet-Sun atmospheres, the cutoff shows a strong dependence on the magnetic field for adiabatic wave propagation. When radiative losses are taken into account, the cutoff frequency is greatly reduced and the variation of the cutoff with the strength of the magnetic field is lower. The effect of the radiative losses in the cutoff is necessary to understand recent quiet-Sun and sunspot observations. In the presence of inclined magnetic fields, our numerical calculations confirm the reduction of the cutoff frequency due to the reduced gravity experienced by waves propagating along field lines. An additional reduction is also found in regions with significant changes in the temperature, due to the lower temperature gradient along the path of field-guided waves. Our results show that the cutoff values are not correctly captured by theoretical estimates. In addition, most of the widely-used analytical cutoff formulae neglect the impact of magnetic fields and radiative losses, whose role is critical to determine the evanescent or propagating nature of the waves.

The analysis of waves in the visible side of the Sun allows the detection of active regions in the farside through local helioseismology techniques. The knowledge of the magnetism in the whole Sun, including the non-visible hemisphere, is fundamental for several space weather forecasting applications. Seismic identification of farside active regions is challenged by the reduced signal-to-noise, and only large and strong active regions can be reliable detected. Here we develop a new methodology to improve the identification of active region signatures in farside seismic maps. We have constructed a deep neural network that associates the farside seismic maps obtained from helioseismic holography with the probability of presence of active regions in the farside. The network has been trained with pairs of helioseismic phase shift maps and Helioseismic and Magnetic Imager magnetograms acquired half a solar rotation later, which were used as a proxy for the presence of active regions in the farside. The method has been validated using a set of artificial data, and it has also been applied to actual solar observations during the period of minimum activity of the solar cycle 24. Our approach shows a higher sensitivity to the presence of farside active regions than standard methods applied up to date. The neural network can significantly increase the number of detected farside active regions, and will potentially improve the application of farside seismology to space weather forecasting.

Imaging spectrographs are popular instruments used to obtain solar data. They record quasi-monochromatic images at selected wavelength positions. By scanning the spectral range of the line, it is possible to obtain bidimensional maps of the FoV with a moderate spectral resolution. In this work, we evaluate the quality of spectropolarimetric inversions obtained from various wavelength samplings during umbral flashes. We computed numerical simulations of nonlinear wave propagation in a sunspot and constructed synthetic Stokes profiles in the Ca II 8542 \AA\ line during an umbral flash using the NLTE code NICOLE. The spectral resolution of the Stokes profiles was downgraded to various cases with differences in the wavelength coverage. A large set of wavelength samplings was analyzed and the performance of the inversions was evaluated by comparing the inferred chromospheric temperature, velocity, and magnetic field with the actual values at the chromosphere of the numerical simulation. The errors in the inverted results depend to a large extent on the location of the wavelength points across the profile of the line. The inferred magnetic field improves with the increase of the spectral resolution. In the case of velocity and temperature, low spectral resolution data produce a match of the inverted atmospheres with the actual values comparable to wavelength samplings with finer resolution, while providing a higher temporal cadence in the data acquisition. We validated the NLTE inversions of spectropolarimetric data from the Ca II 8542 \AA\ during umbral flashes, during which the atmosphere undergoes sudden dramatic changes due to the propagation of a shock wave. Our results favor the use of fine spectral resolution for analyses that focus on the inference of the magnetic field, whereas the estimation of temperature and velocity fluctuations can be performed with lower spectral resolution.

Sunspot umbrae show a change in the dominant period of their oscillations from five minutes in the photosphere to three minutes in the chromosphere. In this paper, we explore the two most popular models proposed to explain the three-minute oscillations: the chromospheric acoustic resonator and the propagation of waves with frequency above the cutoff value directly from lower layers. We employ numerical simulations of wave propagation from the solar interior to the corona. Waves are driven by a piston at the bottom boundary. We have performed a parametric study of the measured chromospheric power spectra in a large number of numerical simulations with differences in the driving method, the height of the transition region (or absence of transition region), the strength of the vertical magnetic field, and the value of the radiative cooling time. We find that both mechanisms require the presence of waves with period in the three-minute band at the photosphere. These waves propagate upward and their amplitude increases due to the drop of the density. Their amplification is stronger than that of evanescent low-frequency waves. This effect is enough to explain the dominant period observed in chromospheric spectral lines. However, waves are partially trapped between the photosphere and the transition region, forming an acoustic resonator. This chromospheric resonant cavity strongly enhances the power in the three-minute band. The chromospheric acoustic resonator model and the propagation of waves in the three-minute band directly from the photosphere can explain the observed chromospheric three-minute oscillations. They are both important in different scenarios. Resonances are produced by waves trapped between the temperature minimum and the transition region. Strong magnetic fields and radiative losses remove energy from the waves inside the cavity, resulting in weaker amplitude resonances.

Solar active regions show a wide variety of oscillatory phenomena. The presence of the magnetic field leads to the appearance of several wave modes, whose behavior is determined by the sunspot thermal and magnetic structure. We aim to study the relation between the umbral and penumbral waves observed at the high photosphere and the magnetic field topology of the sunspot. Observations of the sunspot in active region NOAA 12662 obtained with the GREGOR telescope (Observatorio del Teide, Spain) were acquired on 2017 June 17. The data set includes a temporal series in the Fe I 5435 \AA\ line obtained with the imaging spectrograph GREGOR Fabry-P\'erot Interferometer (GFPI) and a spectropolarimetric raster map acquired with the GREGOR Infrared Spectrograph (GRIS) in the 10830 \AA\ spectral region. The Doppler velocity deduced from the restored Fe I 5435 \AA\ line has been determined, and the magnetic field vector of the sunspot has been inferred from spectropolarimetric inversions of the Ca I 10839 \AA\ and the Si I 10827 \AA\ lines. A two-armed spiral wavefront has been identified in the evolution of the two-dimensional velocity maps from the Fe I 5435 \AA\ line. The wavefronts initially move counterclockwise in the interior of the umbra, and develop into radially outward propagating running penumbral waves when they reach the umbra-penumbra boundary. The horizontal propagation of the wavefronts approximately follows the direction of the magnetic field, which shows changes in the magnetic twist with height and horizontal position. The spiral wavefronts are interpreted as the visual pattern of slow magnetoacoustic waves which propagate upward along magnetic field lines. Their apparent horizontal propagation is due to their sequential arrival to different horizontal positions at the formation height of the Fe I 5435 \AA\ line, as given by the inclination and orientation of the magnetic field.

In the solar atmosphere, the acoustic cutoff frequency is a local quantity which depends on the atmospheric height. It separates the low-frequency evanescent waves from the high-frequency propagating waves. We measure the cutoff frequency of slow magnetoacoustic waves at various heights of a sunspot umbra and compare the results with the estimations from several analytical formulae. We analyzed the oscillations in the umbra of a sunspot belonging to active region NOAA 12662 observed in the 10830 \AA\ spectral region with the GREGOR Infrared Spectrograph and in the Fe I 5435 \AA\ line with the GREGOR Fabry-P\'erot Interferometer. Both instrumets are attached to the GREGOR telescope at the Observatorio del Teide, Tenerife, Spain. We have computed the phase and amplification spectra between the velocity measured from different pairs of lines that sample various heights of the solar atmosphere. The cutoff frequency and its height variation have been estimated from the inspection of the spectra. At the deep umbral photosphere the cutoff frequency is around 5 mHz and it increases to 6 mHz at higher photospheric layers. At the chromosphere the cutoff is $\sim 3.1$ mHz. The comparison of the observationally determined cutoff with the theoretically predicted values reveals an agreement in the general trend and a reasonable match at the chromosphere, but also significant quantitative differences at the photosphere. Our analyses show strong evidence of the variation of the cutoff frequency with height in a sunspot umbra, which is not fully accounted for by current analytical estimations. This result has implications for our understanding of wave propagation, the seismology of active regions, and the evaluation of heating mechanisms based on compressible waves.

The use of instruments that record narrow band images at selected wavelengths is a common approach in solar observations. They allow the scanning of a spectral line by sampling the Stokes profiles with 2D images at each line position, but require a compromise between spectral resolution and temporal cadence. We evaluate the impact of the time-dependent acquisition of different wavelengths on the inversion of spectropolarimetric profiles from chromospheric lines during umbral flashes. Simulations of non-linear wave propagation in a sunspot were performed with the code MANCHA. Synthetic Stokes parameters in the Ca II 8542 A line in NLTE were computed for an umbral flash using the code NICOLE. Artificial profiles with the same wavelength coverage and temporal cadence from reported observations were constructed and inverted. The inferred atmospheric stratifications were compared with the original models. The inferred atmospheres provide a reasonable characterization of the thermodynamic properties of the atmosphere during most of the phases of the umbral flash. Only at the early stages of the flash, when the shock wave reaches the formation height of the line, the Stokes profiles present apparent wavelength shifts and other spurious deformations. These features are misinterpreted by the inversion code, which can return unrealistic atmospheric models from a good fit of the Stokes profiles. The misguided results include flashed atmospheres with strong downflows, even though the simulation exhibits upflows during the umbral flash, and large variations in the magnetic field strength. Our analyses validate the inversion of Stokes profiles acquired by sequentially scanning certain selected wavelengths of a line profile, even in the case of rapidly-changing events such as umbral flashes. However, the inversions are unreliable during a short period at the development phase of the flash.

We investigate the properties of a sunspot light-bridge, focusing on the changes produced by the impact of a plasma blob ejected from a C-class flare. We observed a sunspot in active region NOAA 12544 using spectropolarimetric raster maps of the four Fe I lines around 15655 \AA\ with the GREGOR Infrared Spectrograph (GRIS), narrow-band intensity images sampling the Fe I 6173 \AA\ line with the GREGOR Fabry-P\'erot Interferometer (GFPI), and intensity broad band images in G-band and Ca II H band with the High-resolution Fast Imager (HiFI). All these instruments are located at the GREGOR telescope at the Observatorio del Teide, Tenerife, Spain. The data cover the time before, during, and after the flare event. The analysis is complemented with Atmospheric Imaging Assembly (AIA) and Helioseismic and Magnetic Imager (HMI) data from the Solar Dynamics Observatory (SDO). The physical parameters of the atmosphere at differents heights were inferred using spectral-line inversion techniques. We identify photospheric and chromospheric brightenings, heating events, and changes in the Stokes profiles associated to the flare eruption and the subsequent arrival of the plasma blob to the light bridge, after traveling along an active region loop. The measurements suggest that these phenomena are the result of reconnection events driven by the interaction of the plasma blob with the magnetic field topology of the light bridge.

Improving methods for determining the subsurface structure of sunspots from their seismic signature requires a better understanding of the interaction of waves with magnetic field concentrations. We aim to quantify the impact of changes in the internal structure of sunspots on local helioseismic signals. We have numerically simulated the propagation of a stochastic wave field through sunspot models with different properties, accounting for changes in the Wilson depression between 250 and 550 km and in the photospheric umbral magnetic field between 1500 and 3500 G. The results show that travel-time shifts at frequencies above approximately 3.50 mHz (depending on the phase-speed filter) are insensitive to the magnetic field strength. The travel time of these waves is determined exclusively by the Wilson depression and sound-speed perturbation. The travel time of waves with lower frequencies is affected by the direct effect of the magnetic field, although photospheric field strengths below 1500 G do not leave a significant trace on the travel-time measurements. These results could potentially be used to develop simplified travel-time inversion methods.

Photospheric waves in sunspots moving radially outwards at speeds faster than the characteristic wave velocities have been recently detected. It has been suggested that they are the visual pattern of $p-$modes excited around 5 Mm beneath the sunspot. Using numerical simulations, we have performed a parametric study of the waves observed at the photosphere and higher layers produced by sources located at different depths. The observational measurements are consistent with waves driven between approximately 1 Mm and 5 Mm below the sunspot surface.

Active regions are the most prominent manifestations of solar magnetic fields; their generation and dissipation are fundamental problems in solar physics. Light bridges are commonly present during sunspot decay, but a comprehensive picture of their role in the removal of photospheric magnetic field is still missing. We study the three dimensional configuration of a sunspot and in particular its light bridge during one of the last stages of its decay. We present the magnetic and thermodynamical stratification inferred from full Stokes inversions of the photospheric Si I 10827 \AA\ and Ca I 10839 \AA\ lines obtained with the GREGOR Infrared Spectrograph of the GREGOR telescope at Observatorio del Teide, Tenerife, Spain. The analysis is complemented by a study of continuum images covering the disk passage of the active region, which are provided by the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory. The sunspot shows a light bridge with penumbral continuum intensity that separates the central umbra from a smaller umbra. We find that in this region the magnetic field lines form a canopy with lower magnetic field strength in the inner part. The photospheric light bridge is dominated by gas pressure (high-$\beta$), as opposed to the surrounding umbra where the magnetic pressure is higher. A convective flow is observed in the light bridge. This flow is able to bend the magnetic field lines and to produce field reversals. The field lines close above the light bridge and become as vertical and strong as in the surrounding umbra. We conclude that it develops because of two highly magnetized regions which come closer during the sunspot evolution.

Stellar prominences detected in rapidly rotating stars serve as probes of the magnetism in the corona of cool stars. We have synthesized the temporal evolution of the Stokes profiles generated in the He I 10830 and 5876 A triplets during the rotation of a prominence around a star. The synthesis was performed with the HAZEL code using a cloud model in which the prominence is characterized by a slab located at a fixed latitude and height. It accounts for the scattering polarization and Zeeman and Hanle effects. Several cases with different prominence magnetic field strengths and orientations have been analyzed. The results show an emission feature that drifts across the profile while the prominence is out of the stellar disk. When the prominence eclipses the star, the intensity profile shows an absorption. The scattering induced by the prominence generates linear polarization signals in Stokes Q and U profiles, which are modified by the Hanle effect when a magnetic field is present. Due to the Zeeman effect, Stokes V profiles show a signal with very low amplitude when the magnetic field along the line-of-sight is different from zero. The estimated linear polarization signals could potentially be detected with the future spectropolarimeter MIRADAS, to be attached to GTC telescope.

Wave propagation through sunspots involves conversion between waves of acoustic and magnetic character. In addition, the thermal structure of sunspots is very different than that of the quiet Sun. As a consequence, the interpretation of local helioseismic measurements of sunspots has long been a challenge. With the aim of understanding these measurements, we carry out numerical simulations of wave propagation through sunspots. Helioseismic holography measurements made from the resulting simulated wavefields show qualitative agreement with observations of real sunspots. We use additional numerical experiments to determine, separately, the influence of the thermal structure of the sunspot and the direct effect of the sunspot magnetic field. We use the ray approximation to show that the travel-time shifts in the thermal (non-magnetic) sunspot model are primarily produced by changes in the wave path due to the Wilson depression rather than variations in the wave speed. This shows that inversions for the subsurface structure of sunspots must account for local changes in the density. In some ranges of horizontal phase speed and frequency there is agreement (within the noise level in the simulations) between the travel times measured in the full magnetic sunspot model and the thermal model. If this conclusion proves to be robust for a wide range of models, it would suggest a path towards inversions for sunspot structure.

Spectropolarimetric temporal series from Fe I $\lambda$ 6301.5 \AA\ and Ca II infrared triplet lines are obtained by applying the Stokes synthesis code NICOLE to a numerical simulation of wave propagation in a sunspot umbra from MANCHA code. The analysis of the phase difference between Doppler velocity and intensity core oscillations of the Fe I $\lambda$ 6301.5 \AA\ line reveals that variations in the intensity are produced by opacity fluctuations rather than intrinsic temperature oscillations, except for frequencies between 5 and 6.5 mHz. On the other hand, the photospheric magnetic field retrieved from the weak field approximation provides the intrinsic magnetic field oscillations associated to wave propagation. Our results suggest that this is due to the low magnetic field gradient of our sunspot model. The Stokes parameters of the chromospheric Ca II infrared triplet lines show striking variations as shock waves travel through the formation height of the lines, including emission self-reversals in the line core and highly abnormal Stokes V profiles. Magnetic field oscillations inferred from the Ca II infrared lines using the weak field approximation appear to be related with the magnetic field strength variation between the photosphere and the chromosphere.

The helioseismic properties of the wave scattering generated by monolithic and spaghetti sunspots are analyzed by means of numerical simulations. In these computations, an incident f or p1 mode travels through the sunspot model, which produces absorption and phase shift of the waves. The scattering is studied by inspecting the wavefield, computing travel-time shifts, and performing Fourier-Hankel analysis. The comparison between the results obtained for both sunspot models reveals that the differences in the absorption coefficient can be detected above noise level. The spaghetti model produces an steep increase of the phase shift with the degree of the mode at short wavelengths, while mode-mixing is more efficient for the monolithic model. These results provide a clue for what to look for in solar observations to discern the constitution of sunspots between the proposed monolithic and spaghetti models.

We use numerial simulations to study the absorption and phase shift of surface-gravity waves caused by groups of magnetic flux tubes. The dependence of the scattering coefficients with the distance between the tubes and their positions is analyzed for several cases with two or three flux tubes embedded in a quiet Sun atmosphere. The results are compared with those obtained neglecting completely or partially multiple scattering effects. We show that multiple scattering has a significant impact on the absorption measurements and tends to reduce the phase shift. We also consider more general cases of ensembles of randomly distributed flux tubes, and we have evaluated the effects on the scattering measurements of changing the number of tubes included in the bundle and the average distance between flux tubes. We find that for the longest wavelength incoming waves multiple scattering enhances the absorption, and its efficiency increases with the number of flux tubes and the reduction of the distance between them.

The conversion of fast waves to the Alfven mode in a realistic sunspot atmosphere is studied through three-dimensional numerical simulations. An upward propagating fast acoustic wave is excited in the high-beta region of the model. The new wave modes generated at the conversion layer are analyzed from the projections of the velocity and magnetic field in their characteristic directions, and the computation of their wave energy and fluxes. The analysis reveals that the maximum efficiency of the conversion to the slow mode is obtained for inclinations of 25 degrees and low azimuths, while the Alfven wave conversions peaks at high inclinations and azimuths between 50 and 120 degrees. Downward propagating Alfven waves appear at the regions of the sunspot where the orientation of the magnetic field is in the direction opposite to the wave propagation, since at these locations the Alfven wave couples better with the downgoing fast magnetic wave which are reflected due to the gradients of the Alfven speed. The simulations shows that the Alfven energy at the chromosphere is comparable to the acoustic energy of the slow mode, being even higher at high inclined magnetic fields.

The scattering of f-modes by magnetic tubes is analyzed using three-dimensional numerical simulations. An f-mode wave packet is propagated through a solar atmosphere embedded with three different flux tube models which differ in radius and total magnetic flux. A quiet Sun simulation without a tube present is also performed as a reference. Waves are excited inside the flux tube and propagate along the field lines, and jacket modes are generated in the surroundings of the flux tube, carrying 40% as much energy as the tube modes. The resulting scattered wave is mainly an f-mode composed of a mixture of m=0 and m=+/-1 modes. The amplitude of the scattered wave approximately scales with the magnetic flux. A small amount of power is scattered into the p_1-mode. We have evaluated the absorption and phase shift from a Fourier-Hankel decomposition of the photospheric vertical velocities. They are compared with the results obtained from the emsemble average of 3400 small magnetic elements observed in high-resolution MDI Doppler datacubes. The comparison shows that the observed dependence of the phase shift with wavenumber can be matched reasonably well with the simulated flux tube model. The observed variation of the phase-shifts with the azimuthal order $m$ appears to depend on details of the ensemble averaging, including possible motions of the magnetic elements and asymmetrically shaped elements.

We aim at reproducing the height dependence of sunspot wave signatures obtained from spectropolarimetric observations through 3D MHD numerical simulations. A magneto-static sunspot model based on the properties of the observed sunspot is constructed and perturbed at the photosphere introducing the fluctuations measured with the \SiI\ $\lambda$ 10827 \AA\ line. The results of the simulations are compared with the oscillations observed simultaneously at different heights from the \HeI\ $\lambda$ 10830 \AA\ line, the \CaIIH\ core and the \FeI\ blends in the wings of the \CaIIH\ line. The simulations show a remarkable agreement with the observations. They reproduce the velocity maps and power spectra at the formation heights of the observed lines, as well as the phase and amplification spectra between several pair of lines. We find that the stronger shocks at the chromosphere are accompanied with a delay between the observed signal and the simulated one at the corresponding height, indicating that shocks shift the formation height of the chromospheric lines to higher layers. Since the simulated wave propagation matches very well the properties of the observed one, we are able to use the numerical calculations to quantify the energy contribution of the magneto-acoustic waves to the chromospheric heating in sunspots. Our findings indicate that the energy supplied by these waves is too low to balance the chromospheric radiative losses. The energy contained at the formation height of the lowermost \SiI\ $\lambda$ 10827 \AA\ line in the form of slow magneto-acoustic waves is already insufficient to heat the higher layers, and the acoustic energy which reaches the chromosphere is around 3-9 times lower than the required amount of energy. The contribution of the magnetic energy is even lower.

We study the propagation of waves from the photosphere to the chromosphere of sunspots. From time series of cospatial Ca II H (including its line blends) intensity spectra and polarimetric spectra of Si I 1082.7 nm and He I 1083.0 nm we retrieve the line-of-sight velocity at several heights. The analysis of the phase difference and amplification spectra shows standing waves for frequencies below 4 mHz and propagating waves for higher frequencies, and allows us to infer the temperature and height where the lines are formed. Using these observational data, we have constructed a model of sunspot, and we have introduced the velocity measured with the photospheric Si I 1082.7 nm line as a driver. The numerically propagated wave pattern fits reasonably well with the observed using the lines formed at higher layers, and the simulations reproduce many of the observed features. The observed waves are slow MHD waves propagating longitudinally along field lines.

We analyze the propagation of waves in sunspots from the photosphere to the chromosphere using time series of co-spatial Ca II H intensity spectra (including its line blends) and polarimetric spectra of Si I 10827 and the He I 10830 multiplet. From the Doppler shifts of these lines we retrieve the variation of the velocity along the line-of-sight at several heights. Phase spectra are used to obtain the relation between the oscillatory signals. Our analysis reveals standing waves at frequencies lower than 4 mHz and a continuous propagation of waves at higher frequencies, which steepen into shocks in the chromosphere when approaching the formation height of the Ca II H core. The observed non-linearities are weaker in Ca II H than in He I lines. Our analysis suggests that the Ca II H core forms at a lower height than the He I 10830 line: a time delay of about 20 s is measured between the Doppler signal detected at both wavelengths. We fit a model of linear slow magnetoacoustic wave propagation in a stratified atmosphere with radiative losses according to Newton's cooling law to the phase spectra and derive the difference in the formation height of the spectral lines. We show that the linear model describes well the wave propagation up to the formation height of Ca II H, where non-linearities start to become very important.

Waves observed in the photosphere and chromosphere of sunspots show complex dynamics and spatial patterns. The interpretation of high-resolution sunspot wave observations requires modeling of three-dimensional non-linear wave propagation and mode transformation in the sunspot upper layers in realistic spot model atmospheres. Here we present the first results of such modeling. We have developed a 3D non-linear numerical code specially designed to calculate the response of magnetic structures in equilibrium to an arbitrary perturbation. The code solves the 3D nonlinear MHD equations for perturbations; it is stabilized by hyper-diffusivity terms and is fully parallelized. The robustness of the code is demonstrated by a number of standard tests. We analyze several simulations of a sunspot perturbed by pulses of different periods at subphotospheric level, from short periods, introduced for academic purposes, to longer and realistic periods of three and five minutes. We present a detailed description of the three-dimensional mode transformation in a non-trivial sunspot-like magnetic field configuration, including the conversion between fast and slow magneto-acoustic waves and the Alfv\'en wave, by calculation of the wave energy fluxes. Our main findings are the following: (1) the conversion from acoustic to the Alfv\'en mode is only observed if the the driving pulse is located out of the sunspot axis, but this conversion is energetically inefficient; (2) as a consequence of the cut-off effects and refraction of the fast magneto-acoustic mode, the energy of the evanescent waves with periods around 5 minutes remains almost completely below the level beta=1; (3) waves with frequencies above the cut-off propagate field-aligned to the chromosphere and their power becomes dominating over that of evanescent 5-minute oscillations, in agreement with observations.

Three-dimensional numerical simulations of magnetoacoustic wave propagation are performed in a sunspot atmosphere with a computational domain covering from the photosphere to the chromosphere. The wave source, with properties resembling the solar spectrum, is located at different distances from the axis of the sunspot for each simulation. These results are compared with the theory of mode transformation and also with observational features. Simulations show that the dominant oscillation frequency in the chromosphere decreases with the radial distance from the sunspot axis. The energy flux of the different wave modes involved, including de Alfv\'en mode, is evaluated and discussed.

We present some results obtained from the synthesis of Stokes profiles in small-scale flux tubes with propagating MHD waves. To that aim, realistic flux tubes showing internal structure have been excited with 5 min period drivers, allowing non-linear waves to propagate inside the magnetic structure. The observational signatures of these waves in Stokes profiles of several spectral lines that are commonly used in spectropolarimetric measurements are discussed.

We present results of non-linear, 2D, numerical simulations of magneto-acoustic wave propagation in the photosphere and chromosphere of small-scale flux tubes with internal structure. Waves with realistic periods of three to five minutes are studied, after applying horizontal and vertical oscillatory perturbations to the equilibrium model. Spurious reflections of shock waves from the upper boundary are minimized thanks to a special boundary condition. This has allowed us to increase the duration of the simulations and to make it long enough to perform a statistical analysis of oscillations. The simulations show that deep horizontal motions of the flux tube generate a slow (magnetic) mode and a surface mode. These modes are efficiently transformed into a slow (acoustic) mode in the vA < cS atmosphere. The slow (acoustic) mode propagates vertically along the field lines, forms shocks and remains always within the flux tube. It might deposit effectively the energy of the driver into the chromosphere. When the driver oscillates with a high frequency, above the cut-off, non-linear wave propagation occurs with the same dominant driver period at all heights. At low frequencies, below the cut-off, the dominant period of oscillations changes with height from that of the driver in the photosphere to its first harmonic (half period) in the chromosphere. Depending on the period and on the type of the driver, different shock patterns are observed.