Optimal Masks for Low-Degree Solar Acoustic Modes

We suggest a solution to an important problem of observational helioseismology of the separation of lines of solar acoustic (p) modes of low angular degree in oscillation power spectra by constructing optimal masks for Doppler images of the Sun. Accurate measurements of oscillation frequencies of low-degree modes are essential for the determination of the structure and rotation of the solar core. However, these measurements for a particular mode are often affected by leakage of other p modes arising when the Doppler images are projected on to spherical-harmonics masks. The leakage results in overlaping peaks corresponding to different oscillation modes in the power spectra. In this paper we present a method for calculating optimal masks for a given (target) mode by minimizing the signals of other modes appearing in its vicinity. We apply this method to time series of 2 years obtained from Michelson Doppler Imager (MDI) instrument on board SOHO space mission and demonstrate its ability to reduce efficiently the mode leakage.Comment: to be published in Astrophys.J. Letter

The Cause of Photospheric and Helioseismic Responses to Solar Flares: High-Energy Electrons or Protons?

Analysis of the hydrodynamic and helioseismic effects in the photosphere during the solar flare of July 23, 2002, observed by Michelson Doppler Imager (MDI) on SOHO, and high-energy images from RHESSI shows that these effects are closely associated with sources of the hard X-ray emission, and that there are no such effects in the centroid region of the flare gamma-ray emission. These results demonstrate that contrary to expectations the hydrodynamic and helioseismic responses (''sunquakes") are more likely to be caused by accelerated electrons than by high-energy protons. A series of multiple impulses of high-energy electrons forms a hydrodynamic source moving in the photosphere with a supersonic speed. The moving source plays a critical role in the formation of the anisotropic wave front of sunquakes.Comment: 13 pages, 5 figures, ApJL in pres

Effect of suppressed excitation on the amplitude distribution of 5-min oscillations in sunspots

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.Comment: 12 pages, 5 figures, accepted to ApJ

Travel Time Shifts due to Amplitude Modulation in Time-Distance Helioseismology

Correct interpretation of acoustic travel times measured by time-distance helioseismology is essential to get an accurate understanding of the solar properties that are inferred from them. It has long been observed that sunspots suppress p-mode amplitude, but its implications on travel times has not been fully investigated so far. It has been found in test measurements using a 'masking' procedure, in which the solar Doppler signal in a localized quiet region of the Sun is artificially suppressed by a spatial function, and using numerical simulations that the amplitude modulations in combination with the phase-speed filtering may cause systematic shifts of acoustic travel times. To understand the properties of this procedure, we derive an analytical expression for the cross-covariance of a signal that has been modulated locally by a spatial function that has azimuthal symmetry, and then filtered by a phase speed filter typically used in time-distance helioseismology. Comparing this expression to the Gabor wavelet fitting formula without this effect, we find that there is a shift in the travel times, that is introduced by the amplitude modulation. The analytical model presented in this paper can be useful also for interpretation of travel time measurements for non-uniform distribution of oscillation amplitude due to observational effects.Comment: 17 pages, 1 figure, accepted for publication in Ap

Prediction of Sunspot Cycles by Data Assimilation Method

Despite the known general properties of the solar cycles, a reliable forecast of the 11-year sunspot number variations is still a problem. The difficulties are caused by the apparent chaotic behavior of the sunspot numbers from cycle to cycle and by the influence of various turbulent dynamo processes, which are far from understanding. For predicting the solar cycle properties we make an initial attempt to use the Ensemble Kalman Filter (EnKF), a data assimilation method, which takes into account uncertainties of a dynamo model and measurements, and allows to estimate future observational data. We present the results of forecasting of the solar cycles obtained by the EnKF method in application to a low-mode nonlinear dynamical system modeling the solar $\alpha\Omega$-dynamo process with variable magnetic helicity. Calculations of the predictions for the previous sunspot cycles show a reasonable agreement with the actual data. This forecast model predicts that the next sunspot cycle will be significantly weaker (by $\sim 30$) than the previous cycle, continuing the trend of low solar activity.Comment: 10 pages, 3 figure