We performed 3D numerical simulations of the solar surface wave field for the quiet Sun and for three models with different localized sound-speed variations in the interior with: (i) deep, (ii) shallow, and (iii) two-layer structures. We used simulated data generated by two different codes which use the same standard solar model as a background model, but utilize two different integration techniques and use different models of stochastic wave excitation. Acoustic travel times were measured from all data sets using the time-distance helioseismology technique and compared with the ray theory predictions, frequently used for helioseismic travel-time inversions. It is found that the measured travel-time shifts agree well with the ray theory in both cases with and without phase-speed filtering for the shallow and deep perturbations. This testing verifies the whole measuring-filtering-inversion procedure for sound-speed anomalies inside the Sun. It is shown, that the phase-speed filtering, frequently used to improve the signal-to-noise ratio does not introduce significant systematic errors. Results of the sound-speed inversion procedure show good agreement with the background sound-speed profiles in all cases. Due to its smoothing nature, the inversion procedure overestimates sound speed variations in areas with sharp gradients of the sound-speed profile.

The Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO/HMI) provides continuous full-disk observations of solar oscillations. We develop a data-analysis pipeline based on the time-distance helioseismology method to measure acoustic travel times using HMI Doppler-shift observations, and infer solar interior properties by inverting these measurements. The pipeline is used for routine production of near-real-time full-disk maps of subsurface wave-speed perturbations and horizontal flow velocities for depths ranging from 0 to 20 Mm, every eight hours. In addition, Carrington synoptic maps for the subsurface properties are made from these full-disk maps. The pipeline can also be used for selected target areas and time periods. We explain details of the pipeline organization and procedures, including processing of the HMI Doppler observations, measurements of the travel times, inversions, and constructions of the full-disk and synoptic maps. Some initial results from the pipeline, including full-disk flow maps, sunspot subsurface flow fields, and the interior rotation and meridional flow speeds, are presented.

Local helioseismology provides unique information about the subsurface structure and dynamics of sunspots and active regions. However, because of complexity of sunspot regions local helioseismology diagnostics require careful analysis of systematic uncertainties and physical interpretation of the inversion results. We present new results of comparison of the ring-diagram analysis and time-distance helioseismology for active region NOAA 9787, for which a previous comparison showed significant differences in the subsurface sound-speed structure, and discuss systematic uncertainties of the measurements and inversions. Our results show that both the ring-diagram and time-distance techniques give qualitatively similar results, revealing a characteristic two-layer seismic sound-speed structure consistent with the results for other active regions. However, a quantitative comparison of the inversion results is not straightforward. It must take into account differences in the sensitivity, spatial resolution and the averaging kernels. In particular, because of the acoustic power suppression, the contribution of the sunspot seismic structure to the ring-diagram signal can be substantially reduced. We show that taking into account this effect reduces the difference in the depth of transition between the negative and positive sound-speed variations inferred by these methods. Further detailed analysis of the sensitivity, resolution and averaging properties of the local helioseismology methods is necessary for consolidation of the inversion results. It seems to be important that both methods indicate that the seismic structure of sunspots is rather deep and extends to at least 20 Mm below the surface, putting constraints on theoretical models of sunspots.

We present comparison of numerical simulations of propagation of MHD waves,excited by subphotospheric perturbations, in two different ("deep" and "shallow") magnetostatic models of the sunspots. The "deep" sunspot model distorts both the shape of the wavefront and its amplitude stronger than the "shallow" model. For both sunspot models, the surface gravity waves (f-mode) are affected by the sunspots stronger than the acoustic p-modes. The wave amplitude inside the sunspot depends on the photospheric strength of the magnetic field and the distance of the source from the sunspot axis. For the source located at 9 Mm from the center of the sunspot, the wave amplitude increases when the wavefront passes through the central part of the sunspot. For the source distance of 12 Mm, the wave amplitude inside the sunspot is always smaller than outside. For the same source distance from the sunspot center but for the models with different strength of the magnetic field, the wave amplitude inside the sunspot increases with the strength of the magnetic field. The simulations show that unlike the case of the uniform inclined background magnetic field, the p- and f-mode waves are not spatially separated inside the sunspot where the magnetic field is strongly non-uniform. These properties have to be taken into account for interpretation of observations of MHD waves traveling through sunspot regions.

We use 2D numerical simulations and eikonal approximation, to study properties of MHD waves traveling below the solar surface through the magnetic structure of sunspots. We consider a series of magnetostatic models of sunspots of different magnetic field strengths, from 10 Mm below the photosphere to the low chromosphere. The purpose of these studies is to quantify the effect of the magnetic field on local helioseismology measurements by modeling waves excited by sub-photospheric sources. Time-distance propagation diagrams and wave travel times are calculated for models of various field strength and compared to the non-magnetic case. The results clearly indicate that the observed time-distance helioseismology signals in sunspot regions correspond to fast MHD waves. The slow MHD waves form a distinctly different pattern in the time-distance diagram, which has not been detected in observations. The numerical results are in good agreement with the solution in the short-wavelength (eikonal) approximation, providing its validation. The frequency dependence of the travel times is in a good qualitative agreement with observations.

Investigation of propagation, conversion, and scattering of MHD waves in the Sun is very important for understanding the mechanisms of observed oscillations and waves in sunspots and active regions. We have developed 3D linear MHD numerical model to investigate influence of the magnetic field on excitation and properties of the MHD waves. The results show that the magnetic field can substantially change the properties of the surface gravity waves (f-mode), but their influence on the acoustic-type waves (p-modes) is rather moderate. Comparison our simulations with the time-distance helioseismology results from SOHO/MDI shows that the travel time variations caused by the inclined magnetic field do not exceed 25% of the observed amplitude even for strong fields of 1400-1900 G. This can be an indication that other effects (e.g. background flows and non-uniform distribution of magnetic field) can contribute to the observed travel time variations. The travel time variations caused by the wave interaction with magnetic field are in phase with the observations for strong fields of 1400-1900 G if Doppler velocities are taken at the height of 300 km above the photosphere where plasma parameter beta<<1. The simulations show that the travel times only weakly depend on the height of velocity observation. For the photospheric level the travel times are systematically smaller on approximately 0.12 min then for the hight of 300 km above the photosphere for all studied ranges of the magnetic field strength and inclination angles. The numerical MHD wave modeling and new data from the HMI instrument of the Solar Dynamics Observatory will substantially advance our knowledge of the wave interaction with strong magnetic fields on the Sun and improve the local helioseismology diagnostics.

By analyzing numerically simulated solar oscillation data, we study the influence of non-uniform distribution of acoustic wave amplitude, acoustic source strength, and perturbations of the sound speed on the shifts of acoustic travel times measured by the time-distance helioseismology method. It is found that for short distances, the contribution to the mean travel time shift caused by non-uniform distribution of acoustic sources in sunspots may be comparable to (but smaller than) the contribution from the sound speed perturbation in sunspots, and that it has the opposite sign to the sound-speed effect. This effect may cause some underestimation of the negative sound-speed perturbations in sunspots just below the surface, that was found in previous time-distance helioseismology inferences. This effect cannot be corrected by artificially increasing the amplitude of oscillations in sunspots. For large time-distance annuli, the non-uniform distribution of wavefields does not have significant effects on the mean travel times, and thus the sound-speed inversion results. The measured travel time differences, which are used to determine the mass flows beneath sunspots, can also be systematically shifted by this effect, but only in an insignificant magnitude.

Five-minute oscillations on the Sun (acoustic and surface gravity waves) are excited by subsurface turbulent convection. However, in sunspots the excitation is suppressed because strong magnetic field inhibits convection. We use 3D simulations to investigate how the suppression of excitation sources affects the distribution of the oscillation power in sunspot regions. The amplitude of random acoustic sources was reduced in circular-shaped regions to simulate the suppression in sunspots. The simulation results show that the amplitude of the oscillations can be approximately 2-4 times lower in the sunspot regions in comparison to the quiet Sun, just because of the suppressed sources. Using SOHO/MDI data we measured the amplitude ratio for the same frequency bands outside and inside sunspots, and found that this ratio is approximately 3-4. Hence, the absence of excitation sources inside sunspots makes a significant contribution (about 50% or higher) to the observed amplitude ratio and must be taken into account in sunspot seismology.

Results of numerical 3D simulations of propagation of acoustic waves inside the Sun are presented. A linear 3D code which utilizes realistic OPAL equation of state was developed by authors. Modified convectively stable standard solar model with smoothly joined chromosphere was used as a background model. High order dispersion relation preserving numerical scheme was used to calculate spatial derivatives. The top non-reflecting boundary condition established in the chromosphere absorbs waves with frequencies greater than the acoustic cut-off frequency which pass to the chromosphere, simulating a realistic situation. The acoustic power spectra obtained from the wave field generated by sources randomly distributed below the photosphere are in good agreement with observations. The influence of the height of the top boundary on results of simulation was studied. It was shown that the energy leakage through the acoustic potential barrier damps all modes uniformly and does not change the shape of the acoustic spectrum. So the height of the top boundary can be used for controlling a damping rate without distortion of the acoustic spectrum. The developed simulations provide an important tool for testing local helioseismology.