Search

When simulating full-disc helioseismic data, instrumental noise has traditionally been treated as time-independent. However, in reality, instrumental noise will often vary to some degree over time due to line of sight velocity variations and possibly degrading hardware. Here we present a new technique for simulating Birmingham Solar Oscillations Network (BiSON) helioseismic data with a more realistic analogue for instrumental noise. This is achieved by simulating the potassium solar Fraunhoffer line as observed by the BiSON instruments. Intensity measurements in the red and blue wing of the line can then be simulated and appropriate time-dependent instrumental noise can be added. The simulated time-series can then be formed in the same way as with real data. Here we present the simulation method and the first generation of a BiSON-like instrumental noise time series.

Mode fitting or "peak-bagging" is an important procedure in helioseismology allowing one to determine the various mode parameters of solar oscillations. Here we describe a way of reducing the systematic bias in the fits of certain mode parameters that are seen when using "local" fitting techniques to analyse the sun-as-a-star p-mode power spectrum. To do this we have developed a new "pseudo-global" fitting algorithm designed to gain the advantages of fitting the entire power spectrum, but without the problems involved in fitting a model incorporating many hundreds of parameters. We have performed a comparative analysis between the local and pseudo-global peak-bagging techniques by fitting the "limit" profiles of simulated helioseismic data. Results show that for asymmetric modes the traditional fitting technique returns systematically biased estimates of the central frequency parameter. This bias is significantly reduced when employing the pseudo-global routine. Similarly, we show that estimates of the background returned from the pseudo-global routine match the input values much more closely than the estimates from the local fitting method. We have also used the two fitting techniques to analyse a set of real solar data collected by the Global Oscillations at Low Frequencies (GOLF) instrument on board the ESA/NASA Solar and Heliospheric Observatory (SOHO) spacecraft. Similar differences between the estimated frequencies returned by the two techniques are seen when fitting both the real and simulated data. We show that the background fits returned by the pseudo-global routine more closely match the estimate of the background one can infer from interpolating between fits to the high and low frequency ends of the p-mode power spectrum.

This paper presents an analysis of observational data on the p-mode spectrum of the star alpha Cen B and a comparison with theoretical computations of the stochastic excitation and damping of the modes. We find that at frequencies > 4500 micro-Hz, the model damping rates appear to be too weak to explain the observed shape of the power spectral density of alpha Cen B. The conclusion rests on the assumption that most of the disagreement is due to problems modelling the damping rates, not the excitation rates, of the modes. This assumption is supported by a parallel analysis of BiSON Sun-as-a-star data, for which it is possible to use analysis of very long timeseries to place tight constraints on the assumption. The BiSON analysis shows that there is a similar high-frequency disagreement between theory and observation in the Sun. We demonstrate that by using suitable comparisons of theory and observation it is possible to make inference on the dependence of the p-mode linewidths on frequency, without directly measuring those linewidths, even though the alpha Cen B dataset is only a few nights long. Use of independent measures from a previous study of the alpha Cen B linewidths in two parts of its spectrum also allows us to calibrate our linewidth estimates for the star. The resulting calibrated linewidth curve looks similar to a frequency-scaled version of its solar cousin, with the scaling factor equal to the ratio of the respective acoustic cut-off frequencies of the two stars. The ratio of the frequencies at which the onset of high-frequency problems is seen in both stars is also given approximately by the same scaling factor., Comment: Accepted for publication in ApJ; 19 pages, 6 figures

Unresolved Doppler velocity measurements are not homogenous across the solar disc (Brookes et al. 1978). We consider one cause of the inhomogeneity that originates from the BiSON instrumentation itself: the intensity of light observed from a region on the solar disc is dependent on the distance between that region on the image of the solar disc formed in the instrument and the detector. The non-uniform weighting affects the realization of the solar noise and the amplitudes of the solar oscillations observed by a detector. An 'offset velocity', which varies with time, is observed in BiSON data and has consequences for the long-term stability of observations. We have attempted to model, in terms of the inhomogeneous weighting, the average observed offset velocity., Comment: 4 pages, 4 figures, proceedings of GONG 2008/SOHO XXI meeting

Sun-like oscillations, that is p-modes excited stochastically by convective noise, have now been observed in a number of Red Giant stars. Compared to those seen in the Sun, these modes are of large amplitude and long period, making the oscillations attractive prospects for observation. However, the low Q-factor of these modes, and issues relating to the rising background at low frequencies, present some interesting challenges for identifying modes and determining the related asteroseismic parameters. We report on the analysis procedure adopted for peak-bagging by our group at Birmingham, and the techniques used to robustly ensure these are not a product of noise. I also show results from a number of giants extracted from multi-year observations with the SMEI instrument., Comment: 6 pages, 4 figures. Contribution to the Proceedings of the Wroclaw HELAS Workshop, 2008, to be published in Communications in Asteroseismology

The Birmingham Solar-Oscillations Network (BiSON) has been collecting data for over 30yrs and so observations span nearly three 11yr solar activity cycles. This allows us to address important questions concerning the solar cycle and its effect on solar oscillations, such as: how consistent is the acoustic behaviour from one cycle to the next? We have used the p-mode frequencies observed in BiSON data from one solar activity cycle (cycle 22) to predict the mode frequencies that were observed in the next activity cycle (cycle 23). Some bias in the predicted frequencies was observed when short 108d time series were used to make the predictions. We also found that the accuracy of the predictions was dependent on which activity proxy was used to make the predictions and on the length of the relevant time series., Comment: 6 pages, 4 figures, proceedings of GONG 2008/SOHO XXI meeting

Aims: As g-mode pulsators, gamma-Doradus-class stars may naively be expected to show a large number of modes. Taking advantage of the long photometric time-series generated by the Solar Mass Ejection Imager (SMEI) instrument, we have studied the star gamma Doradus to determine whether any other modes than the three already known are present at observable amplitude. Methods: High-precision photometric data from SMEI taken between April 2003 and March 2006 were subjected to periodogram analysis with the PERIOD04 package. Results: We confidently determine three additional frequencies at 1.39, 1.87, and 2.743 cycles per day. These are above and beyond the known frequencies of 1.320, 1.364, and 1.47 cycles per day. Conclusions: Two of the new frequencies, at 1.39 and 1.87 cycles per day, are speculated to be additional modes of oscillation, with the third frequency at 2.743 cycles per day a possible combination frequency., Comment: 3 pages, 4 figures. Accepted for publication in A&A

We carried out a comparison of the signals seen in contemporaneous BiSON and GOLF data sets. Both instruments perform Doppler shift velocity measurements in integrated sunlight, although BiSON perform measurements from the two wings of potassium absorption line and GOLF from one wing of the NaD1 line. Discrepancies between the two datasets have been observed. We show,in fact, that the relative power depends on the wing in which GOLF data observes. During the blue wing period, the relative power is much higher than in BiSON datasets, while a good agreement has been observed during the red period., Comment: 7 pages, HELAS II: Helioseismology, Asteroseismology, and MHD Connections, conference proceeding

We report on the results of the latest solarFLAG hare-and-hounds exercise, which was concerned with testing methods for extraction of frequencies of low-degree solar p modes from data collected by Sun-as-a-star observations. We have used the new solarFLAG simulator, which includes the effects of correlated mode excitation and correlations with background noise, to make artificial timeseries data that mimic Doppler velocity observations of the Sun as a star. The correlations give rise to asymmetry of mode peaks in the frequency power spectrum. Ten members of the group (the hounds) applied their ``peak bagging'' codes to a 3456-day dataset, and the estimated mode frequencies were returned to the hare (who was WJC) for comparison. Analysis of the results reveals a systematic bias in the estimated frequencies of modes above approximately 1.8 mHz. The bias is negative, meaning the estimated frequencies systematically underestimate the input frequencies. We identify two sources that are the dominant contributions to the frequency bias. Both sources involve failure to model accurately subtle aspects of the observed power spectral density in the part (window) of the frequency power spectrum that is being fitted. One source of bias arises from a failure to account for the power spectral density coming from all those modes whose frequencies lie outside the fitting windows. The other source arises from a failure to account for the power spectral density of the weak l=4 and 5 modes, which are often ignored in Sun-as-a-star analysis. The Sun-as-a-star peak-bagging codes need to allow for both sources, otherwise the frequencies are likely to be biased., Comment: 11 pages; accepted for publication in MNRAS

Asteroseismology of Sun-like stars is undergoing rapid expansion with, for example, new data from the CoRoT mission and continuation of ground-based campaigns. There is also the exciting upcoming prospect of NASA's Kepler mission, which will allow the asteroseismic study of several hundred Sun-like targets, in some cases for periods lasting up to a few years. The seismic mode parameters are the input data needed for making inference on stars and their internal structures. In this paper we discuss the ease with which it will be possible to extract estimates of individual mode parameters, dependent on the mass, age, and visual brightness of the star. Our results are generally applicable; however, we look at mode detectability in the context of the upcoming Kepler observations. To inform our discussions we make predictions of various seismic parameters. To do this we use simple empirical scaling relations and detailed pulsation computations of the stochastic excitation and damping characteristics of the Sun-like p modes. The issues related to parameter extraction on individual p modes discussed here are mode detectability, the detectability and impact of stellar activity cycles, and the ability to measure properties of rotationally split components, which is dependent on the relative importance of the rotational characteristics of the star and the damping of the stochastically excited p modes., Comment: 11 pages; accepted for publication in A&A. Corrections to a few minor typos. Changes to plots to make bold symbols clearer in on-screen PDF

In a previous paper we suggested that, for a given p mode, the excitation function is the same as the component of the solar background noise that has an identical surface spherical harmonic projection (over the corresponding range of temporal frequency). An important consequence of this surmise is that the excitation of overtones of a given angular degree and azimuthal order will be correlated in time. In this note, we introduce the basic principles and a mathematical description of correlated mode excitation. We use simple, illustrative examples, involving two modes. Our treatment suggests that in the real observations, any signatures of the correlation would not appear as a correlation of the output amplitudes of overtones, but rather as subtle modifications to the power spectral density at frequencies between the central frequencies of the overtones. These modifications give a contribution to the observed peak asymmetries., Comment: 7 pages; accepted for publication in Astronomische Nachrichten

Aims: From observations of the K4III star Beta UMi we attempt to determine whether oscillations or any other form of variability is present. Methods: A high-quality photometric time series of approximately 1000 days in length obtained from the SMEI instrument on the Coriolis satellite is analysed. Various statistical tests were performed to determine the significance of features seen in the power density spectrum of the light curve. Results: Two oscillations with frequencies 2.44 and 2.92 microhertz have been identified. We interpret these oscillations as consecutive overtones of an acoustic spectrum, implying a large frequency spacing of 0.48 microhertz. Using derived asteroseismic parameters in combination with known astrophysical parameters, we estimate the mass of Beta UMi to be 1.3 +/- 0.3 solar masses. Peaks of the oscillations in the power density spectrum show width, implying that modes are stochastically excited and damped by convection. The mode lifetime is estimated at 18 +/- 9 days., Comment: 5 Pages, 3 Figures; Correction to minor typo in equation 3

We report on initial results from the first phase of Exercise 1 of the asteroFLAG hare and hounds. The asteroFLAG group is helping to prepare for the asteroseismology component of NASA's Kepler mission, and the first phase of Exercise 1 is concerned with testing extraction of estimates of the large and small frequency spacings of the low-degree p modes from Kepler-like artificial data. These seismic frequency spacings will provide key input for complementing the exoplanet search data., Comment: 11 pages, accepted for publication in Astronomische Nachrichten

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., Comment: 7 pages, 2 figures. To be published in a conference proceedings

We make predictions of the detectability of low-frequency p modes. Estimates of the powers and damping times of these low-frequency modes are found by extrapolating the observed powers and widths of higher-frequency modes with large observed signal-to-noise ratios. The extrapolations predict that the low-frequency modes will have small signal-to-noise ratios and narrow widths in a frequency-power spectrum. Monte Carlo simulations were then performed where timeseries containing mode signals and normally distributed Gaussian noise were produced. The mode signals were simulated to have the powers and damping times predicted by the extrapolations. Various statistical tests were then performed on the frequency-amplitude spectra formed from these timeseries to investigate the fraction of spectra in which the modes could be detected. The results of these simulations were then compared to the number of p-modes candidates observed in real Sun-as-a-star data at low frequencies. The fraction of simulated spectra in which modes were detected decreases rapidly as the frequency of modes decreases and so the fraction of simulations in which the low-frequency modes were detected was very small. However, increasing the signal-to-noise (S/N) ratio of the low-frequency modes by a factor of 2 above the extrapolated values led to significantly more detections. Therefore efforts should continue to further improve the quality of solar data that is currently available., Comment: 9 pages, 6 figures, accepted for publication in Astronomische Nachrichten

Frequencies, powers and damping rates of the solar p modes are all observed to vary over the 11-yr solar activity cycle. Here, we show that simultaneous variations of these parameters give rise to a subtle cross-talk effect, which we call the ``devil in the detail'', that biases p-mode frequencies estimated from analysis of long power frequency spectra. We also show that the resonant peaks observed in the power frequency spectra show small distortions due to the effect. Most of our paper is devoted to a study of the effect for Sun-as-a-star observations of the low-l p modes. We show that for these data the significance of the effect is marginal. We also touch briefly on the likely l dependence of the effect, and discuss the implications of these results for solar structure inversions., Comment: 9 pages; accepted for publication in MNRAS

We seek to address whether solar-cycle frequency shifts of the Sun's low-l p modes `distort' the underlying shapes of the mode peaks, when those peaks are observed in power frequency spectra made from data spanning large fractions, or more, of the cycle period. We present analytical descriptions of the expected profiles, and validate the predictions through use of artificial seismic timeseries data, in which temporal variations of the oscillator frequencies are introduced. Our main finding is that for the Sun-like frequency shifts the distortion of the asymmetrical Lorentzian-like profiles is very small, but also just detectible. Our analysis indicates that by fitting modes to the usual Lorentzian-like models -- which do not allow for the distortion -- rather than new models we derive, there is a bias in the mode height and linewidth parameters that is comparable in size to the observational uncertainties given by multi-year datasets. Bias in the frequency parameter gives much less cause for worry, being over an order of magnitude smaller than the corresponding frequency uncertainties. The distortion discussed in this paper may need to be considered when multi-year Sun-like asteroseismic datasets are analyzed on stars showing strong activity cycles., Comment: 7 pages; accepted for publication in MNRAS

We present new results on oscillations of the K1.5 III giant Arcturus (alpha Boo), from analysis of just over 2.5 yr of precise photometric observations made by the Solar Mass Ejection Imager (SMEI) on board the Coriolis satellite. A strong mode of oscillation is uncovered by the analysis, having frequency 3.51+/-0.03 micro-Hertz. By fitting its mode peak, we are able offer a highly constrained direct estimate of the damping time (tau = 24+/-1 days). The data also hint at the possible presence of several radial-mode overtones, and maybe some non-radial modes. We are also able to measure the properties of the granulation on the star, with the characteristic timescale for the granulation estimated to be 0.50+/-0.05 days., Comment: 6 pages, 5 figures; accepted for publication in MNRAS Letters

We use 9.5-yr of BiSON Sun-as-a-star data to search for dependence of solar-cycle parameter changes on the angular degree, l, of the data. The nature of the Sun-as-a-star observations is such that for changes measured at fixed frequency, or for changes averaged across the same range in frequency, any l dependence present carries information on the latitudinal distribution of the agent (i.e., the activity) responsible for those changes. We split the 9.5-yr timeseries into contiguous 108-d pieces, and determine mean changes in the damping of, power in, and energy supplied to the modes through the solar cycle. We also apply a careful correction to account for the deleterious effects of the ground-based BiSON window function on the results. From our full analysis we obtain a marginally significant result for the damping parameter, where the mean change is found to be weakest at l=0. The other parameters show hints of some dependence in l. Our main conclusion is that the mean fractional solar-cycle change in the l=0 damping rates is approximately 50 % smaller than was previously assumed. It had been common practice to use an average over all low-l modes; our downward revision of the radial-mode value has implications for comparisons with models of the global solar cycle changes, which are usually based on a spherically symmetric geometry., Comment: 8 pages, 3 figures, 2 tables; Submitted to A&A

We present results of fitting the 50-day time series of photometry of alpha Cen A taken by the WIRE satellite in 1999. Both power spectrum and autocovariance function (ACF) fitting techniques were used in an attempt to determine mode frequencies, rotational splittings, lifetimes and amplitudes of low-l p-modes. In all, using both techniques, we managed to fit 18 modes (seven l = 0, eight l = 1 and three l = 2) with frequencies determined to within 1 - 2 micro-Hz. These estimates are shown to be 0.6 +/- 0.3 micro-Hz lower, on average, than the frequencies determined from two other more recent studies (Bouchy & Carrier, 2002; Bedding et al. 2004) which used data gathered about 19 months after the WIRE observations. This could be indicative of an activity cycle, although due to the large uncertainty, more data would be needed to confirm this. Over a range of 1700 to 2650 micro-Hz we were also able to use the ACF fitting to determine an average lifetime of 3.9 +/- 1.4 days, and an average rotational splitting of 0.54 +/- 0.22 micro-Hz, which is the first ever reliable estimate of this parameter. In contrast to the ACF, the power spectrum fitting was shown to return significantly biased results for these parameters., Comment: 11 pages, 10 figures (some figures are made up from multiple ps files of which there are 17 in total)

We report on results from the first solar Fitting at Low-Angular degree Group (solar FLAG) hare-and-hounds exercise. The group is concerned with the development of methods for extracting the parameters of low-l solar p mode data (`peak bagging'), collected by Sun-as-a-star observations. Accurate and precise estimation of the fundamental parameters of the p modes is a vital pre-requisite of all subsequent studies. Nine members of the FLAG (the `hounds') fitted an artificial 3456-d dataset. The dataset was made by the `hare' (WJC) to simulate full-disc Doppler velocity observations of the Sun. The rotational frequency splittings of the l=1, 2 and 3 modes were the first parameter estimates chosen for scrutiny. Significant differences were uncovered at l=2 and 3 between the fitted splittings of the hounds. Evidence is presented that suggests this unwanted bias had its origins in several effects. The most important came from the different way in which the hounds modeled the visibility ratio of the different rotationally split components. Our results suggest that accurate modelling of the ratios is vital to avoid the introduction of significant bias in the estimated splittings. This is of importance not only for studies of the Sun, but also of the solar analogues that will targets for asteroseismic campaigns.

The fine-structure spacing $d_{\ell}(n) = \nu_{\ell,n} - \nu_{\ell+2,n-1}$ for low-degree solar p modes of angular degree $\ell$ and radial order n, is sensitive to conditions in the deep radiative interior of the Sun. Here, we present fine-structure spacings derived from the analysis of nearly five years of helioseismological data collected between 1991 July and 1996 February by the Birmingham Solar-Oscillations Network (BiSON). These data cover $9 \le n \le 28$ for $d_{0}(n)$, and $11 \le n \le 27$ for $d_{1}(n)$. The measured spacings are much more precise, and cover a greater range, than earlier measurements from BiSON data (Elsworth et al. 1990a). The predicted fine-structure spacings for a ``standard'' solar model are clearly excluded by the BiSON data (at $\approx 10\sigma$); models that include helium and heavy element settling provide a much better match to the observed spacings (see also Elsworth et al. 1995). Since the inclusion of core settling in solar models will tend to slightly increase the predicted neutrino flux, the BiSON fine-structure data appear to reinforce previous conclusions, i.e., an astrophysical solution to the solar neutrino problem seems unlikely., Comment: 14 pages, LaTex 2e, 1 figure (2 encapsulated .PS files); ApJ Letters, in press

Observations of the Sun with the LOWL instrument provide a homogeneous set of solar p-mode frequencies from low to intermediate degree which allows one to determine the structure of much of the solar interior avoiding systematic errors that are introduced when different data sets are combined, i.e., principally the effects of solar cycle changes on the frequencies. Unfortunately, the LOWL data set contains very few of the lowest-degree modes, which are essential for determining reliably the structure of the solar core - in addition, these lowest-degree data have fairly large associated uncertainties. However, observations made by the Birmingham Solar-Oscillations Network (BiSON) in integrated sunlight provide high-accuracy measurements of a large number of low-degree modes. In this paper we demonstrate that the low-degree mode set of the LOWL data can be successfully combined with the more accurate BiSON data, provided the observations are contemporaneous for those frequencies where the solar-cycle-induced effects are important. We show that this leads to a factor-of-two decrease in the error on the inferred sound speed in the solar core. We find that the solar sound speed is higher than in solar models for r < 0.2Rsun. The density of the solar core is, however, lower than that in solar models., Comment: Plain TeX, using MN.tex, v. 1.6. 10 pages, 10 figures. Submitted to Mon. Not. R. astr. Soc