Optimal Masks for Low-Degree Solar Acoustic Modes

We suggest a solution to an important problem in 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 harmonic masks. The leakage results in overlapping peaks corresponding to different oscillation modes in the power spectra. In this Letter, 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 yr obtained from the Michelson Doppler Imager instrument on board the Solar and Heliospheric Observatory space mission and demonstrate its ability to reduce efficiently the mode leakage.     

Global-Oscillation Eigenfunction Measurements of Solar Meridional Flow

We describe and apply a new helioseismic method for measuring solar subsurface axisymmetric meridional and zonal flow. The method is based on a theoretical model of the response of global-oscillation eigenfunctions to the flow velocity and uses cross spectra of the time-varying coefficients in the spherical-harmonic expansion of the photospheric Doppler-velocity field. Eigenfunction changes modify the leakage matrix, which describes the sensitivity of the spherical-harmonic coefficients to the global-oscillation modes. The form of the leakage matrix in turn affects the theoretically expected spherical-harmonic cross spectra. Estimates of internal meridional and zonal flow were obtained by fitting the theoretical flow-dependent cross spectra to spherical-harmonic cross spectra computed from approximately 500 days of full-disk Dopplergrams from the Helioseismic and Magnetic Imager (HMI) on the SDO spacecraft. The zonal-flow measurements, parameterized in the form of “a” coefficients, substantially agree with measurements obtained from conventional global-mode-frequency analysis. The meridional-flow estimates, in the form of depth-weighted averages of the flow velocity, are similar to estimates obtained from earlier analyses, for oscillation modes that penetrate the outermost one-third of the convection zone. For more deeply penetrating modes, the inferred flow velocity increases significantly with penetration depth, indicating the need for either a modification of the simple conveyor-belt picture of meridional flow or improvement in the cross-spectral model.     

On model predictions of the power spectral density of radial solar p modes

We investigate the frequency dependence of the power spectral density of low-degree solar p modes by comparing measurements with the results of a stochastic-excitation model. In the past it was common practice to use the total power in such investigations. Using the maximum of the power spectral density instead provides a direct comparison with the measured mode heights in the observed power spectrum. This method permits a more careful calibration of the adjustable parameters in the excitation model, a model which we present here, for the first time, in a format that precisely and unambiguously relates the amplitudes of the modes of oscillation to the Reynolds stress in the equilibrium model. We find that errors in the theory of the linear mode damping rates, particularly at low frequency, have a dramatic impact on the predictions of the mode heights in the spectral density, whereas parameter changes in the stochastic excitation model, within a plausible domain of parameter space, have a comparatively small effect.     

On measuring low-degree p-mode frequency splitting with full-disc integrated data

The standard method of measuring rotational splitting from solar full-disc oscillation data, based on maximum-likelihood fitting of multi-Lorentzian profiles to oscillation power spectra, systematically overestimates the splitting. One of the reasons is that the maximum likelihood estimators (MLE) become unbiased only asymptotically as the number of data tends to infinity; for a finite data set they are often biased, inducing a systematic error. In this paper we assess by Monte Carlo simulations the amount of systematic error in the splitting measurement, using artificially generated power spectra. The simulations are carried out for multiplets of degree     2 and 3 with various signal-to-noise ratios, linewidths and observing times. We address the possible use of non-MLE estimators that could provide a smaller or negligible systematic error. The implication for asteroseismology is also discussed.     

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.     

Variability of Solar Five-Minute Oscillations in the Corona as Observed by the Extreme Ultraviolet Spectrophotometer (ESP) on the Solar Dynamics Observatory/Extreme Ultraviolet Variability Experiment (SDO/EVE)

Solar five-minute oscillations have been detected in the power spectra of two six-day time intervals from soft X-ray measurements of the Sun observed as a star using the Extreme Ultraviolet Spectrophotometer (ESP) onboard the Solar Dynamics Observatory (SDO)/Extreme Ultraviolet Variability Experiment (EVE). The frequencies of the largest amplitude peaks were found to match the known low-degree (?=0?–?3) modes of global acoustic oscillations within 3.7 μHz and can be explained by a leakage of the global modes into the corona. Due to the strong variability of the solar atmosphere between the photosphere and the corona, the frequencies and amplitudes of the coronal oscillations are likely to vary with time. We investigated the variations in the power spectra for individual days and their association with changes of solar activity, e.g. with the mean level of the EUV irradiance, and its short-term variations caused by evolving active regions. Our analysis of samples of one-day oscillation power spectra for a 49-day period of low and intermediate solar activity showed little correlation with the mean EUV irradiance and the short-term variability of the irradiance. We suggest that some other changes in the solar atmosphere, e.g., magnetic fields and/or inter-network configuration may affect the mode leakage to the corona.     

Self-similar pressure oscillations in neutron star envelopes as probes of neutron star structure

We study eigenmodes of acoustic oscillations of high multipolarity l ∼ 100–1000 and high frequency (∼100 kHz), localized in neutron star envelopes. We show that the oscillation problem is self-similar. Once the oscillation spectrum is calculated for a given equation of state (EOS) in the envelope and given stellar mass M and radius R , it can be rescaled to a star with any M and R (but the same EOS in the envelope). For l ≳ 300, the modes can be subdivided into the outer and inner ones. The outer modes are mainly localized in the outer envelope. The inner modes are mostly localized near the neutron drip point, being associated with the softening of the EOS after the neutron drip. We calculate oscillation spectra for the EOSs of cold-catalyzed and accreted matter and show that the spectra of the inner modes are essentially different. A detection and identification of high-frequency pressure modes would allow one to infer M and R and determine also the EOS in the envelope (accreted or ground state) providing a new, potentially powerful method to explore the main parameters and internal structure of neutron stars.     

Gapfilling interrupted helioseismic data with the EM algorithm

Helioseismic data are often interrupted by gaps, which diminish the quality of the data. In the frequency domain, these gaps lead to systematical effects with misleading interpretation of the power spectra. We propose a gap filling method that is based on modeling solar oscillation data with a statistical process, i.e., the stochastic nature of a single oscillation is taken into account by regarding it as realization of an autoregressive (AR) processes of second order. From the whole oscillation time series given as the superposition of the realization of many excited modes, the process parameters are estimated via the expectation maximization (EM) algorithm. Then the estimated model is used to predict the further course of the oscillatory process during occurring gaps. We demonstrate the applicability of this procedure on the basis of both simulations and data obtained with the DIFOS satellite experiment suffering from gaps of 30 min duration occurring regularly every 90 min due to the orbit around the Earth.     

Fourier Analysis of Gapped Time Series: Improved Estimates of Solar and Stellar Oscillation Parameters

Quantitative helioseismology and asteroseismology require very precise measurements of the frequencies, amplitudes, and lifetimes of the global modes of stellar oscillation. The precision of these measurements depends on the total length (T), quality, and completeness of the observations. Except in a few simple cases, the effect of gaps in the data on measurement precision is poorly understood, in particular in Fourier space where the convolution of the observable with the observation window introduces correlations between different frequencies. Here we describe and implement a rather general method to retrieve maximum likelihood estimates of the oscillation parameters, taking into account the proper statistics of the observations. Our fitting method applies in complex Fourier space and exploits the phase information. We consider both solar-like stochastic oscillations and long-lived harmonic oscillations, plus random noise. Using numerical simulations, we demonstrate the existence of cases for which our improved fitting method is less biased and has a greater precision than when the frequency correlations are ignored. This is especially true of low signal-to-noise solar-like oscillations. For example, we discuss a case where the precision of the mode frequency estimate is increased by a factor of five, for a duty cycle of 15%. In the case of long-lived sinusoidal oscillations, a proper treatment of the frequency correlations does not provide any significant improvement; nevertheless, we confirm that the mode frequency can be measured from gapped data with a much better precision than the 1/T Rayleigh resolution.     

Unstable non-radial modes in radial pulsators: theory and an example

We study the possibility of the excitation of non-radial oscillations in classical pulsating stars. The stability of an RR Lyrae model is examined through non-adiabatic non-radial calculations. We also explore stability in the presence of non-linear coupling between radial and non-radial modes of nearly identical frequency.   In our model, a large number of unstable low-degree (ℓ = 1,2) modes have frequencies in the vicinity of unstable radial mode frequencies. The growth rates of such modes, however, are considerably smaller than those of the radial modes. We also recover an earlier result that at higher degrees (ℓ = 6–12) there are modes trapped in the envelope with growth rates similar to those of radial modes.   Subsequently, monomode radial pulsation of this model is considered. The destabilizing effect of the 1:1 resonance between the radial mode and nearby non-radial modes of low degrees is studied, with the assumption that the excited radial mode saturates the linear instability of all other modes. The instability depends on the radial mode amplitude, the frequency difference, the damping rate of the non-radial mode, and the strength of the non-linear coupling between the modes considered. At the pulsation amplitudes typical for RR Lyrae stars, the instability of the monomode radial pulsation and the concomitant resonant excitation of some non-radial oscillation modes is found to be very likely.     

Spatial distributions of sunspot oscillation modes at different temperatures

Three-and five-minute sunspot oscillations have different spatial distributions in the solar atmospheric layers. The spatial distributions are crucial for revealing the physical origin of sunspot oscillations and to investigate their propagation. In this study, six sunspots observed by Solar Dynamics Observatory/Atmospheric Imaging Assembly were used to obtain the spatial distributions of three-and five-minute oscillations. The fast Fourier transform method is applied to represent the power spectra of oscillation modes. We find that, from the temperature minimum to the lower corona, the powers of the fiveminute oscillation exhibit a circle-shape distribution around its umbra, and the shapes gradually expand with temperature increase. However, the circle-shape disappears and the powers of the oscillations appear to be very disordered in the higher corona. This indicates that the five-minute oscillation can be suppressed in the high-temperature region. For the three-minute oscillations, from the temperature minimum to the high corona, their powers mostly distribute within an umbra, and part of them are located at the coronal fan loop structures. Moreover, those relative higher powers are mostly concentrated in the position of coronal loop footpoints.     

Gaussian process modelling of asteroseismic data

The measured properties of stellar oscillations can provide powerful constraints on the internal structure and composition of stars. To begin this process, oscillation frequencies must be extracted from the observational data, typically time series of the star's brightness or radial velocity. In this paper, a probabilistic model is introduced for inferring the frequencies and amplitudes of stellar oscillation modes from data, assuming that there is some periodic character to the oscillations, but that they may not be exactly sinusoidal. Effectively, we fit damped oscillations to the time series, and hence the mode lifetime is also recovered. While this approach is computationally demanding for large time series (>1500 points), it should at least allow improved analysis of observations of solar-like oscillations in subgiant and red giant stars, as well as sparse observations of semiregular stars, where the number of points in the time series is often low. The method is demonstrated on simulated data and then applied to radial velocity measurements of the red giant star  ξ Hydrae  , yielding a mode lifetime between 0.41 and 2.65 d with 95 per cent posterior probability. The large frequency separation between modes is ambiguous, however we argue that the most plausible value is 6.3 μHz, based on the radial velocity data and the star's position in the Hertzsprung–Russell diagram.     

Constraints on black hole spins with a general relativistic accretion disk corona model

The peaks in the spectra of the accretion disks surrounding massive black holes in quasars are in the far-UV or soft X-ray band, which are usually not observed. However, in the disk corona model, soft photons from the disk are Comptonized to high energy in the hot corona, and the hard X-ray spectra(luminosity and spectral shape) contain information on the incident spectra from the disk. The values of black hole spin parameter a*are inferred from the spectral fitting, which are spread over a large range, ~-0.94 to 0.998. We find that the inclination angles and mass accretion rates are well determined by the spectral fitting, but the results are sensitive to the accuracy of black hole mass estimates. No tight constraints on the black hole spins are achieved, if the uncertainties in black hole mass measurements are a factor of four,which are typical for the single-epoch reverberation mapping method. Recently, the accuracy of black hole mass measurement has been significantly improved to 0.2- 0.4 dex with the velocity resolved reverberation mapping method. The black hole spin can be well constrained if the mass measurement accuracy is50%. In the accretion disk corona scenario, a fraction of power dissipated in the disk is transported into the corona, and therefore the accretion disk is thinner than a bare disk for the same mass accretion rate,because the radiation pressure in the disk is reduced. We find that the thin disk approximation, H/R0.1,is still valid if 0.3 m 0.5, provided half of the dissipated power is radiated in the corona above the disk.     

STRUCTURE AND ROTATION OF THE SOLAR INTERIOR: INITIAL RESULTS FROM THE MDI MEDIUM-L PROGRAM

The medium-l program of the Michelson Doppler Imager instrument on board SOHO provides continuous observations of oscillation modes of angular degree, l, from 0 to 300. The data for the program are partly processed on board because only about 3% of MDI observations can be transmitted continuously to the ground. The on-board data processing, the main component of which is Gaussian-weighted binning, has been optimized to reduce the negative influence of spatial aliasing of the high-degree oscillation modes. The data processing is completed in a data analysis pipeline at the SOI Stanford Support Center to determine the mean multiplet frequencies and splitting coefficients. The initial results show that the noise in the medium-l oscillation power spectrum is substantially lower than in ground-based measurements. This enables us to detect lower amplitude modes and, thus, to extend the range of measured mode frequencies. This is important for inferring the Sun's internal structure and rotation. The MDI observations also reveal the asymmetry of oscillation spectral lines. The line asymmetries agree with the theory of mode excitation by acoustic sources localized in the upper convective boundary layer. The sound-speed profile inferred from the mean frequencies gives evidence for a sharp variation at the edge of the energy-generating core. The results also confirm the previous finding by the GONG (Gough et al., 1996) that, in a thin layer just beneath the convection zone, helium appears to be less abundant than predicted by theory. Inverting the multiplet frequency splittings from MDI, we detect significant rotational shear in this thin layer. This layer is likely to be the place where the solar dynamo operates. In order to understand how the Sun works, it is extremely important to observe the evolution of this transition layer throughout the 11-year activity cycle.     

Cosmological constraints on agegraphic dark energy in DGP braneworld gravity

A proposal to study the original and new agegraphic dark energy in DGP braneworld cosmology is presented in this work. To verify our model with the observational data, the model is constrained by a variety of independent measurements such as Hubble parameter, cosmic microwave background anisotropies, and baryon acoustic oscillation peaks. The best fitting procedure shows the effectiveness of agegraphic parameter n in distinguishing between the original and new agegraphic dark energy scenarios and subsequent cosmological findings. In particular, the result shows that in both scenarios, our universe enters an agegraphic dark energy dominated phase.     

A modified peak‐bagging technique for fitting low‐𝓁 solar p‐modes

We introduce a modified version of a standard power spectrum ‘peak‐bagging’ technique which is designed to gain some of the advantages that fitting the entire low‐degree p‐mode power spectrum simultaneously would bring, but without the problems involved in fitting a model incorporating many hundreds of parameters. Employing Monte‐Carlo simulations we show that by using this modified fitting code it is possible to determine the true background level in the vicinity of the p‐mode peaks. In addition to this we show how small biases in other mode parameters, which are related to inaccurate estimates of the true background, are also consequently removed. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Variation of the low-degree solar acoustic mode parameters over the solar cycle

VIRGO/SPM is a helioseismic sunphotometer on board SOHO that observes the disk-integrated sunlight irradiance at three different colors (red, green, and blue). The data obtained for SPM since the beginning of the SOHO mission, April 1996, to March 2001 have been used to study the differences of the p-mode parameters during the solar activity cycle. These time series have been divided in sub-series of 100 days, transformed to power spectra and averaged in sets of three to yield a total number of six averaged power spectra (around one per year). A new way of analyzing the power spectrum has been applied to the six power spectra of each color; it consists of fitting the whole p-mode spectrum at once with a unique background. The results for the frequencies, line widths, power, mode energy, energy rate fed in the mode and splittings along the activity cycle are found, compared and discussed.     

Helioseismic Holography of an Artificial Submerged Sound Speed Perturbation and Implications for the Detection of Pre-emergence Signatures of Active Regions

We use a publicly available numerical wave-propagation simulation of Hartlep et al. (Solar Phys. 268, 321, 2011) to test the ability of helioseismic holography to detect signatures of a compact, fully submerged, 5 % sound-speed perturbation placed at a depth of 50 Mm within a solar model. We find that helioseismic holography employed in a nominal “lateral-vantage” or “deep-focus” geometry employing quadrants of an annular pupil can detect and characterize the perturbation. A number of tests of the methodology, including the use of a plane-parallel approximation, the definition of travel-time shifts, the use of different phase-speed filters, and changes to the pupils, are also performed. It is found that travel-time shifts made using Gabor-wavelet fitting are essentially identical to those derived from the phase of the Fourier transform of the cross-covariance functions. The errors in travel-time shifts caused by the plane-parallel approximation can be minimized to less than a second for the depths and fields of view considered here. Based on the measured strength of the mean travel-time signal of the perturbation, no substantial improvement in sensitivity is produced by varying the analysis procedure from the nominal methodology in conformance with expectations. The measured travel-time shifts are essentially unchanged by varying the profile of the phase-speed filter or omitting the filter entirely. The method remains maximally sensitive when applied with pupils that are wide quadrants, as opposed to narrower quadrants or with pupils composed of smaller arcs. We discuss the significance of these results for the recent controversy regarding suspected pre-emergence signatures of active regions.     

The Photospheric Convection Spectrum

Spectra of the cellular photospheric flows are determined from observations acquired by the MDI instrument on the SOHO spacecraft. Spherical harmonic spectra are obtained from the full-disk observations. Fourier spectra are obtained from the high-resolution observations. The p-mode oscillation signal and instrumental artifacts are reduced by temporal filtering of the Doppler data. The resulting spectra give power (kinetic energy) per wave number for effective spherical harmonic degrees from 1 to over 3000. Significant power is found at all wavenumbers, including the small wavenumbers representative of giant cells. The time evolution of the spectral coefficients indicates that these small wavenumber components rotate at the solar rotation rate and thus represent a component of the photospheric cellular flows. The spectra show distinct peaks representing granules and supergranules but no distinct features at wavenumbers representative of mesogranules or giant cells. The observed cellular patterns and spectra are well represented by a model that includes two distinct modes – granules and supergranules.     

Measurements of Frequencies of Solar Oscillations from the Mdi Medium-l Program

Inversions of solar internal structure employ both the frequencies and the associated uncertainties of the solar oscillation modes as input parameters. In this paper we investigate how systematic errors in these input parameters may affect the resulting inferences of the sun's internal structure. Such systematic errors are likely to arise from inaccuracies in the theoretical models which are used to represent the spectral lines in the observational power spectra, from line blending, from asymmetries in the profiles of these lines, and from other factors. In order to study such systematic effects we have employed two different duration observing runs (one of 60 days and the second of 144 days) obtained with the Medium-l Program of the Michelson Doppler Imager experiment onboard the SOHO spacecraft. This observing program provides continuous observations of solar oscillation modes having angular degrees, l, ranging from 0 to ∼ 300. For this study intermediate- and high-degree p-mode oscillations having degrees less than 251 were employed. In the first of our tests we employed two different methods of estimating the modal frequencies and their associated uncertainties from the 144-day observational power spectra. In our second test we also repeated both methods of frequency estimation on the 60-day time series in order to assess the influence of the duration of the observed time series on the computed frequencies and uncertainties. In a third test we investigated the sensitivity of the computed frequencies to the choice of initial-guess, or ‘seed’ frequencies that are used in the frequency estimation codes. In a fourth test we attempted to investigate the possible systematic frequency errors which are introduced when the observational asymmetry in the p-mode peaks is ignored. We carried out this particular test by fitting simple models of asymmetric line profiles to the peaks in the observational power spectra. We were then able to compute the differences between those frequencies and our previous frequencies which had been obtained using the assumption that all of the observational peaks were symmetric in shape. In order to study the possible influence of the two different frequency estimation methods upon the radial profile of the internal sound speed, we carried out four parallel structural inversions using the different sets and subsets of frequency estimates and uncertainties as computed from the 144-day observing run as inputs. The results of these four inversions confirm the previous finding by the GONG project (Gough et al., 1996) and by the MDI Medium-l Program (Kosovichev et al., 1997) that, in a thin layer just beneath the convection zone, helium appears to be less abundant than predicted by theory. However, differences in our four inverted radial sound speed profiles demonstrate that the currently-available techniques for determining the frequencies of the Medium-l oscillation peaks introduce systematic errors which are large enough to affect the results of the structural inversions. Moreover, based upon the differences in these four inverted sound speed profiles, it appears that the choice of which subset of modes is included in a particular inversion and which modes are not included may also be introducing systematic errors into our current understanding of solar internal structure. Hence, it appears to be very important that consistent sets of modal selection criteria be employed. Finally, at least one of the two frequency estimation codes which we used was not sensitive to changes in the input ‘seed’ frequencies which were employed as initial guesses for that code. This result allays fears that the difference in the helium abundance between the sun and the reference solar model in the thin layer beneath the convection zone which was mentioned above might have been due to the particular seed frequencies which were employed in the earlier inversions. Since this thin layer may likely be the place where the solar dynamo operates, it will be extremely important to observe any possible evolution of this transition layer throughout the upcoming 11-year activity cycle. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1004963425123