Toward Waveform Heliotomography: Observing Interactions of Helioseismic Waves with a Sunspot

We investigate how helioseismic waves that originate from effective point sources interact with a sunspot. These waves are reconstructed from observed stochastic wavefields on the Sun by cross-correlating photospheric Doppler-velocity signals. We select the wave sources at different locations relative to the sunspot, and investigate the p- and f-mode waves separately. The results reveal a complicated picture of waveform perturbations caused by the wave interaction with the sunspot. In particular, it is found that for waves originating from outside of the sunspot, p-mode waves travel across the sunspot with a small amplitude reduction and slightly higher speed, and wave amplitude and phase get mostly restored to the quiet-Sun values after passing the sunspot. The f-mode wave experiences some amplitude reduction passing through the sunspot, and the reduced amplitude is not recovered after that. The wave-propagation speed does not change before encountering the sunspot and inside the sunspot, but the wavefront becomes faster than the reference wave after passing through the sunspot. For waves originating from inside the sunspot umbra, both f- and p-mode waves show significant amplitude reductions and faster speed for all propagation paths. A comparison of positive and negative time lags of cross-correlation functions shows an apparent asymmetry in the waveform changes for both the f- and p-mode waves. We suggest that the waveform variations of the helioseismic waves interacting with a sunspot found in this article can be used for developing a method of waveform heliotomography, similar to the waveform tomography of the Earth.     

Magnetic and thermal phase shifts in the local helioseismology of sunspots

Phase perturbations due to inclined surface magnetic field of active region strength are calculated numerically in quiet Sun and simple sunspot models in order to estimate and compare the direct and indirect (thermal) effects of the fields on helioseismic waves. It is found that the largest direct effects occur in highly inclined field characteristic of penumbrae, and scale roughly linearly with magnetic field strength. The combined effects of sunspot magnetic and thermal anomalies typically yield negative travel-time perturbations in penumbrae. Travel-time shifts in umbrae depend on details of how the thermal and density structure differs from the quiet Sun. The combined shifts are generally not well approximated by the sum of the thermal and magnetic effects applied separately, except at low field strengths of around 1 kG or less, or if the thermal shift is small. A useful rule-of-thumb appears to be that travel-time perturbations in umbrae are predominantly thermal, whereas in penumbrae they are mostly magnetic.     

Surface-Focused Seismic Holography of Sunspots: I. Observations

We present a comprehensive set of observations of the interaction of p-mode oscillations with sunspots using surface-focused seismic holography. Maps of travel-time shifts, relative to quiet-Sun travel times, are shown for incoming and outgoing p modes as well as their mean and difference. We compare results using phase-speed filters with results obtained with filters that isolate single p-mode ridges, and we further divide the data into multiple temporal frequency bandpasses. The f mode is removed from the data. The variations of the resulting travel-time shifts with magnetic-field strength and with the filter parameters are explored. We find that spatial averages of these shifts within sunspot umbrae, penumbrae, and surrounding plage often show strong frequency variations at fixed phase speed. In addition, we find that positive values of the mean and difference travel-time shifts appear exclusively in waves observed with phase-speed filters that are dominated by power in the low-frequency wing of the p ₁ ridge. We assess the ratio of incoming to outgoing p-mode power using the ridge filters and compare surface-focused holography measurements with the results of earlier published p-mode scattering measurements using Fourier?–?Hankel decomposition.     

Properties of High-Frequency Wave Power Halos Around Active Regions: An Analysis of Multi-height Data from HMI and AIA Onboard SDO

We study properties of waves of frequencies above the photospheric acoustic cut-off of ≈5.3 mHz, around four active regions, through spatial maps of their power estimated using data from the Helioseismic and Magnetic Imager (HMI) and Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO). The wavelength channels 1600 Å and 1700 Å from AIA are now known to capture clear oscillation signals due to helioseismic p-modes as well as waves propagating up through to the chromosphere. Here we study in detail, in comparison with HMI Doppler data, properties of the power maps, especially the so-called “acoustic halos” seen around active regions, as a function of wave frequencies, inclination, and strength of magnetic field (derived from the vector-field observations by HMI), and observation height. We infer possible signatures of (magneto)acoustic wave refraction from the observation-height-dependent changes, and hence due to changing magnetic strength and geometry, in the dependences of power maps on the photospheric magnetic quantities. We discuss the implications for theories of p-mode absorption and mode conversions by the magnetic field.     

The Effect of Abnormal Granulation on Acoustic Wave Travel Times and Mode Frequencies

Observations indicate that in plage areas (i.e. in active regions outside sunspots) acoustic waves travel faster than in the quiet Sun, leading to shortened travel times and higher p-mode frequencies. Coupled with the 11-year variation of solar activity, this may also explain the solar cycle variation of oscillation frequencies. While it is clear that the ultimate cause of any difference between the quiet Sun and plage is the presence of magnetic fields of order 100 G in the latter, the mechanism by which the magnetic field exerts its influence has not yet been conclusively identified. One possible such mechanism is suggested by the observation that granular motions in plage areas tend to be slightly “abnormal”, dampened compared to the quiet Sun. In this paper we consider the effect that abnormal granulation observed in active regions should have on the propagation of acoustic waves. Any such effect is found to be limited to a shallow surface layer where sound waves propagate nearly vertically. The magnetically suppressed turbulence implies higher sound speeds, leading to shorter travel times. This time shift Δ τ is independent of the travel distance, while it shows a characteristic dependence on the assumed plage field strength. As a consequence of the variation of the acoustic cutoff with height, Δ τ is expected to be significantly higher for higher frequency waves within the observed regime of 3 – 5 mHz. The lower group velocity near the upper reflection point further leads to an increased envelope time shift, as compared to the phase shift. p-mode frequencies in plage areas are increased by a corresponding amount, Δ ν/ν=ν Δ τ. These characteristics of the time and frequency shifts are in accordance with observations. The calculated overall amplitudes of the time and frequency shifts are comparable to, but still significantly less than (by a factor of 2 to 5), those suggested by measurements.     

Mode Conversion of Solar <Emphasis Type="Italic">p</Emphasis>-Modes in Non-Vertical Magnetic Fields

Sunspots absorb and scatter incident f- and p-modes. Until recently, the responsible absorption mechanism was uncertain. The most promising explanation appears to be mode conversion to slow magnetoacoustic-gravity waves, which carry energy down the magnetic field lines into the interior. In vertical magnetic field, mode conversion can adequately explain the observed f-mode absorption, but is too inefficient to account for the absorption of p-modes. In the first paper of the present series we calculated the efficiency of fast-to-slow magnetoacoustic-gravity wave conversion in uniform non-vertical magnetic fields. We assumed two-dimensional propagation, where the Alfvén waves decouple. In comparison to vertical field, it was found that mode conversion is significantly enhanced in moderately inclined fields, especially at higher frequencies. Using those results, Cally, Crouch, and Braun showed that the resultant p-mode absorption produced by simple sunspot models with non-vertical magnetic fields is ample to explain the observations. In this paper, we further examine mode conversion in non-vertical magnetic fields. In particular, we consider three-dimensional propagation, where the fast and slow magnetoacoustic-gravity waves and the Alfvén waves are coupled. Broadly speaking, the p-mode damping rates are not substantially different to the two-dimensional case. However, we do find that the Alfvén waves can remove similar quantities of energy to the slow MAG waves.     

Diagnostic features of sunspot regions

Five-minutes p-mode oscillations are heavily attenuated in the active sunspot region. A comparative study of wave modes and luminosity variations outside and inside the sunspot region is found to depict certain diagnostic features of sunspot regions.     

On the variations in the parameters of high-degree acoustic modes with solar cycle

We analyze the pattern of behavior of p-mode wave packets with solar cycle using TON one-day helioseismic data with a high spatial resolution. The time—distance method is used to perform this task. We make an attempt to determine the variations in the travel time of acoustic waves at maximum and minimum solar activity; at maximum activity, this time decreases by 2 s compared to that at minimum activity to a depth of 0.8R_⊙. In addition, the correlation amplitudes of acoustic wave packets from minimum to maximum solar activity were found to decrease by 10–20% for all angular distances.     

Sunspot as an isolated magnetic structure: Stability and oscillations

A simple energy model of a sunspot as a compact magnetic feature is described where the main energy contribution is provided by the coolest and most compressed part of the magnetic force tube of the spot at depths ranging from Wilson’s depression level (300–500 km) down to 2–3 thousand km. The equilibrium and stability conditions for such a system are analyzed using the variation principle, and oscillations of the system as a whole about the inferred equilibrium position are studied. The sunspot is shown to be stable in the magnetic field strength interval from 0.8–1 to 4–5 kG. The dependence of the eigenfrequency on magnetic field strength ω(B) is computed for the main oscillatory mode, where only the umbra of the sunspot takes part in oscillations, ω = ω ₁ (B). Lower subharmonics may appear in the case where penumbra too becomes involved in the oscillatory process: ω ₂ = ω ₁/2, ω ₃ = ω ₁/3. Theoretical curves agree well with the observational data obtained in Pulkovo using various independent methods: from temporal variations of sunspot magnetic field and from line-of-sight-velocity measurements. The periods of oscillations found range from 40 to 200 minutes.     

Helioseismology of Sunspots: Confronting Observations with Three-Dimensional MHD Simulations of Wave Propagation

The propagation of solar waves through the sunspot of AR?9787 is observed by using temporal cross-correlations of SOHO/MDI Dopplergrams. We then use three-dimensional MHD numerical simulations to compute the propagation of wave packets through self-similar magnetohydrostatic sunspot models. The simulations are set up in such a way as to allow a comparison with observed cross-covariances (except in the immediate vicinity of the sunspot). We find that the simulation and the f-mode observations are in good agreement when the model sunspot has a peak field strength of 3 kG at the photosphere and less so for lower field strengths. Constraining the sunspot model with helioseismology is only possible because the direct effect of the magnetic field on the waves has been fully taken into account. Our work shows that the full-waveform modeling of sunspots is feasible.     

磁场中p模/s模的转换及黑子对声波的吸收

在柱坐标下将黑子周围的环形区域（黑子除外）内的振荡分解为朝向黑子传播的（入射的）波和离开黑子传播的（出射的）波。对无黑子的环形区域内的振荡也进行了同样的分解。将黑子周围的入射波看成是被黑子磁流管磁化了的介质（介质内的磁场基本是水平的）中的波。而无黑子区的入射波看成是非磁化介质中的波。比较这两种波在固定波数下功率随频率的分布发现，在磁化介质中不同径向除ｎ的声波（ｐ模）频率系统降低，同时功率也降低，降低的功率最高达非磁化介质中波的功率的３０％。而比较在固定频率下功率随波数的分布发现，磁场中ｆ模及ｎ＝１，２，３的ｐ模的脊向高波数方向位移，功率的降低受频率调制，即声波在某些有限的频带中被吸收。这些观测表明，在磁场中ｐ模与磁声重力波（ＭＡＧ）产生了模式混合或耦合。模式混合的存在支持了模式转换作为ｐ模式被黑子吸收的机制的解释。此外，本文还分析了转换的ＭＡＧ波进入黑子磁流管（其中的磁场基本上是垂直的）后进一步被吸收，吸收的功率最高达ＭＡＧ波的２０％。在磁流管内没有进一步观测到模式的转换     

Time – Distance Modelling in a Simulated Sunspot Atmosphere

In time?–?distance helioseismology, wave travel times are measured from the cross-correlation between Doppler velocities recorded at any two locations on the solar surface. However, one of the main uncertainties associated with such measurements is how to interpret observations made in regions of strong magnetic field. Isolating the effects of the magnetic field from thermal or sound-speed perturbations has proved to be quite complex and has yet to yield reliable results when extracting travel times from the cross-correlation function. One possible way to decouple these effects is by using a 3D sunspot model based on observed surface magnetic-field profiles, with a surrounding stratified, quiet-Sun atmosphere to model the magneto-acoustic ray propagation, and analyse the resulting ray travel-time perturbations that will directly account for wave-speed variations produced by the magnetic field. These artificial travel-time perturbation profiles provide us with several related but distinct observations: i) that strong surface magnetic fields have a dual effect on helioseismic rays?–?increasing their skip distance while at the same time speeding them up considerably compared to their quiet-Sun counterparts, ii) there is a clear and significant frequency dependence of both skip-distance and travel-time perturbations across the simulated sunspot radius, iii) the negative sign and magnitude of these perturbations appears to be directly related to the sunspot magnetic-field strength and inclination, iv) by “switching off” the magnetic field inside the sunspot, we are able to completely isolate the thermal component of the travel-time perturbations observed, which is seen to be both opposite in sign and much smaller in magnitude than those measured when the magnetic field is present. These results tend to suggest that purely thermal perturbations are unlikely to be the main effect seen in travel times through sunspots, and that strong, near-surface magnetic fields may be directly and significantly altering the magnitude and lateral extent of sound-speed inversions of sunspots made by time?–?distance helioseismology.     

Microwave emission from a sunspot

From the gyroresonance brightness temperature spectrum of a sunspot, one can determine the magnetic field strength by using the property that microwave brightness is limited above a frequency given by an integer-multiple of the gyrofrequency. In this paper, we use this idea to find the radial distribution of magnetic field at the coronal base of a sunspot in the active region, NOAA 4741. The gyroresonance brightness temperature spectra of this sunspot are obtained from multi-frequency interferometric observations made at the Owens Valley Radio Observatory at 24 frequencies in the range of 4.0–12.4 GHz with spatial resolution 2.2″–6.8″. The main results of present study are summarized as follows: first, by comparison of the coronal magnetic flux deduced from our microwave observation with the photospheric magnetic flux measured by KPNO magnetograms, we show that theo-mode emission must arise predominantly from the second harmonic of the gyrofrequency, while thex-mode arises from the third harmonic. Second, the radial distribution of magnetic fieldsB(r) at the coronal base of this spot (say, 2000–4000 km above the photosphere) can be adequately fitted by $$B(r) = 1420(1 \pm 0.080)\exp \left[ { - \left( {\frac{r}{{11.05''(1 \pm 0.014)}}} \right)^2 } \right]G,$$ wherer is the radial distance from the spot center at coronal base. Third, it is found that coronal magnetic fields originate mostly from the photospheric umbral region. Fourth, although the derived vertical variation of magnetic fields can be approximated roughly by a dipole model with dipole moment 1.6 × 10³⁰ erg G^?1 buried at 11000 km below the photosphere, the radial field distribution at coronal heights is found to be more confined than predicted by the dipole model.     

Constructing Semi-Empirical Sunspot Models for Helioseismology

One goal of helioseismology is to determine the subsurface structure of sunspots. In order to do so, it is important to understand first the near-surface effects of sunspots on solar waves, which are dominant. Here we construct simplified, cylindrically-symmetric sunspot models that are designed to capture the magnetic and thermodynamics effects coming from about 500 km below the quiet-Sun τ ₅₀₀₀=1 level to the lower chromosphere. We use a combination of existing semi-empirical models of sunspot thermodynamic structure (density, temperature, pressure): the umbral model of Maltby et al. (1986, Astrophys. J. 306, 284) and the penumbral model of Ding and Fang (1989, Astron. Astrophys. 225, 204). The OPAL equation-of-state tables are used to derive the sound-speed profile. We smoothly merge the near-surface properties to the quiet-Sun values about 1 Mm below the surface. The umbral and penumbral radii are free parameters. The magnetic field is added to the thermodynamic structure, without requiring magnetostatic equilibrium. The vertical component of the magnetic field is assumed to have a Gaussian horizontal profile, with a maximum surface field strength fixed by surface observations. The full magnetic-field vector is solenoidal and determined by the on-axis vertical field, which, at the surface, is chosen such that the field inclination is 45° at the umbral – penumbral boundary. We construct a particular sunspot model based on SOHO/MDI observations of the sunspot in active region NOAA 9787. The helioseismic signature of the model sunspot is studied using numerical simulations of the propagation of f, p ₁, and p ₂ wave packets. These simulations are compared against cross-covariances of the observed wave field. We find that the sunspot model gives a helioseismic signature that is similar to the observations.     

Magnetoacoustic-gravity surface waves

The linearized theory for the parallel propagation of magnetoacoustic-gravity surface waves is developed for an interface of a horizontal magnetic field above a field-free medium. The media either side of the interface are taken to be isothermal. The dispersion relation is obtained for the case of a constant Alfvén speed. In the absence of gravity the interface may support one or two surface modes, determined by the relative temperatures and magnetism of the two media. The effect of gravity on the modes is examined and dispersion diagrams and eigenfunctions are given. In the usual - k _x diagnostic diagram, the domain of evanescence is shown to be divided into two distinct regions determining whether a given mode will have a decaying or growing vertical velocity component. In the absence of a magnetic field the transcendental dispersion relation may be rewritten as a polynomial. This polynomial possesses two acceptable solutions only one of which may exist in any given circumstances (depending on the ratio of the densities). If the gas density within the field exceeds that in the field-free medium, then the f-mode may propagate. The f-mode exists in a restricted band of horizontal wavenumber and only when the field-free medium is warmer than the magnetic atmosphere. An analytical form for the wave speed of the f-mode is obtained for small values of the Alfvén speed. It is shown that the f-mode is related to the fast magnetoacoustic surface wave, merging into that mode at short wavelengths.     

Physical Properties of Wave Motion in Inclined Magnetic Fields within Sunspot Penumbrae

At the surface of the Sun, acoustic waves appear to be affected by the presence of strong magnetic fields in active regions. We explore the possibility that the inclined magnetic field in sunspot penumbrae may convert primarily vertically-propagating acoustic waves into elliptical motion. We use helioseismic holography to measure the modulus and phase of the correlation between incoming acoustic waves and the local surface motion within two sunspots. These correlations are modeled by assuming the surface motion to be elliptical, and we explore the properties of the elliptical motion on the magnetic-field inclination. We also demonstrate that the phase shift of the outward-propagating waves is opposite to the phase shift of the inward-propagating waves in stronger, more vertical fields, but similar to the inward phase shifts in weaker, more-inclined fields.     

Enhanced p-Mode Absorption Seen Near the Sunspot Umbral – Penumbral Boundary

We investigate p-mode absorption in a sunspot using SOHO/MDI high-resolution Doppler images. The Doppler power computed from a 3.5-hour data set is used for studying the absorption in a sunspot. The result shows an enhancement in absorption near the umbral?–?penumbral boundary of the sunspot. We attempt to relate the observed absorption with the magnetic-field structure of the sunspot. The transverse component of the potential field is computed by using the observed SOHO/MDI line-of-sight magnetograms. A comparison of the power map and the computed potential field shows enhanced absorption near the umbral?–?penumbral boundary where the computed transverse field strength is higher.     

Modelling p-Mode Interaction with a Spreading Sunspot Field

Sunspots absorb and scatter incident p modes. The dominant mechanism is still uncertain. One possibility, mode conversion to slow magneto-acoustic waves, has been shown to yield results in agreement with observations for the f mode only. Absorption of p modes in simple vertical magnetic field models is too weak by an order of magnitude or more. Here we report on numerical calculations of p modes encountering a simple sunspot model with field which spreads with height. It is found that p-mode absorption is greatly enhanced by field spread, to a level consistent with observations, and it appears that it occurs preferentially in the outer regions of the spot, in line with recent results from acoustic holography.     

On Absorption and Scattering of P Modes by Small-Scale Magnetic Elements

The solar surface is characterised everywhere by the presence of small-scale magnetic structures. Their collective behaviour in the form of active regions is known to have strong influence on p-mode power. For example, sunspots and plages are strong absorbers of acoustic waves. This paper studies the effects of individual small-scale magnetic elements to understand the details of absorption of p-mode power. For this, we consider a thin magnetic flux tube and calculate the phase shifts and the absorption coefficients by numerically solving the linearised MHD equations. The phase shifts calculated from the Born Approximation are then compared for the same range of degrees. The results are discussed with a view to understanding the physical mechanism.     

Large-scale horizontal flows in the solar photosphere IV. On the vertical structure of large-scale horizontal flows

In the recent papers, we introduced a method utilised to measure the flow field. The method is based on the tracking of supergranular structures. We did not precisely know, whether its results represent the flow field in the photosphere or in some subphotospheric layers. In this paper, in combination with helioseismic data, we are able to estimate the depths in the solar convection envelope, where the detected large-scale flow field is well represented by the surface measurements. We got a clear answer to question what kind of structures we track in full-disc Dopplergrams. It seems that in the quiet Sun regions the supergranular structures are tracked, while in the regions with the magnetic field the structures of the magnetic field are dominant. This observation seems obvious, because the nature of Doppler structures is different in the magnetic regions and in the quiet Sun. We show that the large-scale flow detected by our method represents the motion of plasma in layers down to ～10 Mm. The supergranules may therefore be treated as the objects carried by the underlying large-scale velocity field.