No Arabic abstract
Helioseismic holography is an imaging technique used to study heterogeneities and flows in the solar interior from observations of solar oscillations at the surface. Holograms contain noise due to the stochastic nature of solar oscillations. We provide a theoretical framework for modeling signal and noise in Porter-Bojarski helioseismic holography. The wave equation may be recast into a Helmholtz-like equation, so as to connect with the acoustics literature and define the holography Greens function in a meaningful way. Sources of wave excitation are assumed to be stationary, horizontally homogeneous, and spatially uncorrelated. Using the first Born approximation we calculate holograms in the presence of perturbations in sound-speed, density, flows, and source covariance, as well as the noise level as a function of position. This work is a direct extension of the methods used in time-distance helioseismology to model signal and noise. To illustrate the theory, we compute the hologram intensity numerically for a buried sound-speed perturbation at different depths in the solar interior. The reference Greens function is obtained for a spherically-symmetric solar model using a finite-element solver in the frequency domain. Below the pupil area on the surface, we find that the spatial resolution of the hologram intensity is very close to half the local wavelength. For a sound-speed perturbation of size comparable to the local spatial resolution, the signal-to-noise ratio is approximately constant with depth. Averaging the hologram intensity over a number $N$ of frequencies above 3 mHz increases the signal-to-noise ratio by a factor nearly equal to the square root of $N$. This may not be the case at lower frequencies, where large variations in the holographic signal are due to the individual contributions of the long-lived modes of oscillation.
We present an adaptation of the rotation-corrected, m-averaged spectrum technique designed to observe low signal-to-noise-ratio, low-frequency solar p modes. The frequency shift of each of the 2l+1 m spectra of a given (n,l) multiplet is chosen that maximizes the likelihood of the m-averaged spectrum. A high signal-to-noise ratio can result from combining individual low signal-to-noise-ratio, individual-m spectra, none of which would yield a strong enough peak to measure. We apply the technique to GONG and MDI data and show that it allows us to measure modes with lower frequencies than those obtained with classic peak-fitting analysis of the individual-m spectra. We measure their central frequencies, splittings, asymmetries, lifetimes, and amplitudes. The low-frequency, low- and intermediate-angular degrees rendered accessible by this new method correspond to modes that are sensitive to the deep solar interior down to the core and to the radiative interior. Moreover, the low-frequency modes have deeper upper turning points, and are thus less sensitive to the turbulence and magnetic fields of the outer layers, as well as uncertainties in the nature of the external boundary condition. As a result of their longer lifetimes (narrower linewidths) at the same signal-to-noise ratio the determination of the frequencies of lower-frequency modes is more accurate, and the resulting
Helioseismology provides important constraints for the solar dynamo problem. However, the basic properties and even the depth of the dynamo process, which operates also in other stars, are unknown. Most of the dynamo models suggest that the toroidal magnetic field that emerges on the surface and forms sunspots is generated near the bottom of the convection zone, in the tachocline. However, there is a number of theoretical and observational problems with justifying the deep-seated dynamo models. This leads to the idea that the subsurface angular velocity shear may play an important role in the solar dynamo. Using helioseismology measurements of the internal rotation and meridional circulation, we investigate a mean-field MHD model of dynamo distributed in the bulk of the convection zone but shaped in a near-surface layer. We show that if the boundary conditions at the top of the dynamo region allow the large-scale toroidal magnetic fields to penetrate into the surface, then the dynamo wave propagates along the isosurface of angular velocity in the subsurface shear layer, forming the butterfly diagram in agreement with the Parker-Yoshimura rule and solar-cycle observations. Unlike the flux-transport dynamo models, this model does not depend on the transport of magnetic field by meridional circulation at the bottom of the convection zone, and works well when the meridional circulation forms two cells in radius, as recently indicated by deep-focus time-distance helioseismology analysis of the SDO/HMI and SOHO/MDI data. We compare the new dynamo model with various characteristics if the solar magnetic cycles, including the cycle asymmetry (Waldmeiers relations) and magnetic `butterfly diagrams.
We study Doppler velocity measurements at multiple heights in the solar atmosphere using a set of six filtergrams obtained by the Helioseismic magnetic Imager on board the Solar Dynamics Observatory. There are clear and significant phase differences between core and wing Dopplergrams in the frequency range above the photospheric acoustic cutoff frequency, which indicates that these are really multi-height datasets.
We present a 3-dimensional (3D) numerical solver of the linearized compressible Euler equations (GALE -- Global Acoustic Linearized Euler), used to model acoustic oscillations throughout the solar interior. The governing equations are solved in conservation form on a fully global spherical mesh ($0 le phi le 2pi$, $0 le theta le pi$, $0 le r le R_{odot}$) over a background state generated by the standard Solar Model S. We implement an efficient pseudo-spectral computational method to calculate the contribution of the compressible material derivative dyad to internal velocity perturbations, computing oscillations over arbitrary 3D background velocity fields. This model offers a foundation for a forward-modeling approach, using helioseismology techniques to explore various regimes of internal mass flows. We demonstrate the efficacy of the numerical method presented in this paper by reproducing observed solar power spectra, showing rotational splitting due to differential rotation, and applying local helioseismology techniques to measure travel times created by a simple model of single-cell meridional circulation.
We looked for signatures of Quasi-Biennial Periodicity (QBP) over different phases of solar cycle by means of acoustic modes of oscillation. Low-degree p-mode frequencies are shown to be sensitive to changes in magnetic activity due to the global dynamo. Recently have been reported evidences in favor of two-year variations in p-mode frequencies. Long high-quality helioseismic data are provided by BiSON (Birmingham Solar Oscillation Network), GONG (Global Oscillation Network Group), GOLF (Global Oscillation at Low Frequency) and VIRGO (Variability of Solar IRradiance and Gravity Oscillation) instruments. We determined the solar cycle changes in p-mode frequencies for spherical degree l=0, 1, 2 with their azimuthal components in the frequency range 2.5 mHz < nu < 3.5 mHz. We found signatures of QBP at all levels of solar activity in the modes more sensitive to higher latitudes. The signal strength increases with latitude and the equatorial component seems also to be modulated by the 11-year envelope. The persistent nature of the seismic QBP is not observed in the surface activity indices, where mid-term variations are found only time to time and mainly over periods of high activity. This feature together with the latitudinal dependence provides more evidences in favor of a mechanism almost independent and different from the one that brings up to the surface the active regions. Therefore, these findings can be used to provide more constraints on dynamo models that consider a further cyclic component on top of the 11-year cycle.