Multitaper Spectral Analysis and Wavelet Denoising Applied to Helioseismic Data

Details Multitaper Spectral Analysis and Wavelet Denoising Applied to Helioseismic Data Multitaper Spectral Analysis and Wavelet Denoising Applied to Helioseismic Data

Jul 1999

adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=1999apj...519..407k&link_type=abstract

The Astrophysical Journal, Volume 519, Issue 1, pp. 407-421.

Astronomy and Astrophysics

Astronomy

18

Sun: Oscillations, Techniques: Image Processing

Scientific paper

Estimates of solar normal mode frequencies from helioseismic observations can be improved by using multitaper spectral analysis (MTSA) to estimate spectra from the time series, then using wavelet denoising of the log spectra. MTSA leads to a power spectrum estimate with reduced variance and better leakage properties than the conventional periodogram. Under the assumption of stationarity and mild regularity conditions, the log multitaper spectrum has a statistical distribution that is approximately Gaussian, so wavelet denoising is asymptotically an optimal method to reduce the noise in the estimated spectra. We find that a single m-nu spectrum benefits greatly from MTSA followed by wavelet denoising and that wavelet denoising by itself can be used to improve m-averaged spectra. We compare estimates using two different five-taper estimates (Slepian and sine tapers) and the periodogram estimate for Global Oscillation Network Group (GONG) time series at selected angular degrees l. We compare those three spectra with and without wavelet denoising, both visually and in terms of the mode parameters estimated from the preprocessed spectra using the GONG peak-fitting algorithm. The two multitaper estimates give equivalent results. The number of modes fitted well by the GONG algorithm is 20%-60% larger (depending on l and the temporal frequency) when applied to the multitaper estimates than when applied to the periodogram. The estimated mode parameters (frequency, amplitude, and width) are comparable for the three power spectrum estimates, except for modes with very small mode widths (a few frequency bins), where the multitaper spectra broaden the modes compared with the periodogram. At frequencies below 3 mHz, wavelet denoising of the log multitaper power spectra tends to increase the number of modes for which the GONG peak-fitting algorithm converges well. Close to 3 mHz, where all modes are resolved, wavelet denoising makes little difference. At higher frequencies close to the acoustic cutoff frequency, where modes are blended into ridges, wavelet denoising the multitaper spectra reduces the number of good fits. We tested the influence of the number of tapers used and found that narrow modes at low n-values are broadened to the extent that they can no longer be fitted if the number of tapers is too large. For helioseismic time series of this length and temporal resolution, the optimal number of tapers is less than 10.

Affiliated with

Also associated with

No associations

Rating

Multitaper Spectral Analysis and Wavelet Denoising Applied to Helioseismic Data does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.

If you have personal experience with Multitaper Spectral Analysis and Wavelet Denoising Applied to Helioseismic Data, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Multitaper Spectral Analysis and Wavelet Denoising Applied to Helioseismic Data will most certainly appreciate the feedback.

Rate now

Comments { 0 }

Profile ID: LFWR-SCP-O-1619373