Some peculiarities of the power spectrum of the 5-min solar brightness oscillations

Several results of the solar irradiance observations carried out with the IPHIR experiment on board the Soviet spacecraft within the PHOBOS Project are considered. The time variations of the power spectrum of 5-min solar oscillations is discussed. It is found that the lines are generated inhomogeneously in the spectrum as a whole. Each of the frequencies is excited separately at a definite time interval. The duration of active periods does not exceed 2 or 3 hours.     

Influence of Low-Degree High-Order p-Mode Splittings on the Solar Rotation Profile

The solar rotation profile is well constrained down to about 0.25R _⊙ thanks to the study of acoustic modes. Since the radius of the inner turning point of a resonant acoustic mode is inversely proportional to the ratio of its frequency to its degree, only the low-degree p modes reach the core. The higher the order of these modes, the deeper they penetrate into the Sun and thus they carry more diagnostic information on the inner regions. Unfortunately, the estimates of frequency splittings at high frequency from Sun-as-a-star measurements have higher observational errors because of mode blending, resulting in weaker constraints on the rotation profile in the inner core. Therefore inversions for the solar internal rotation use only modes below 2.4 mHz for ?≤3. In the work presented here, we used an 11.5-year-long time series to compute the rotational frequency splittings for modes ?≤3 using velocities measured with the GOLF instrument. We carried out a theoretical study of the influence of the low-degree modes in the region from 2 to 3.5 mHz on the inferred rotation profile as a function of their error bars.     

First View of the Solar Core from Golf Acoustic Modes

After 8 months of nearly continuous measurements the GOLF instrument, aboard SOHO, has detected acoustic mode frequencies of more than 100 modes, extending from 1.4 mHz to 4.9 mHz. In this paper, we compare these results with the best available predictions coming from solar models. To verify the quality of the data, we examine the asymptotic seismic parameters; this confirms the improvements achieved in solar models during the last decade. Using the GOLF set of frequencies for l=0, 1, 2, 3 combined with the LOWL second year data set for l > 3 we then carry out inversions to infer properties of the solar core. This largely confirms the previous results down to around 0.1 R⊙, while there remain differences, even closer to the centre, where the present study shows an extreme sensitivity of the inversion results to the values of the frequencies. We finally consider physical processes which may influence directly or indirectly the solar core structure.     

Properties of Flares-Generated Seismic Waves on the Sun

The purpose of this article is to carry out a power-spectrum analysis of the Super-Kamiokande five-day dataset that takes account of the asymmetry in the error estimates. Whereas for symmetrical error estimates the likelihood analysis involves a linear optimization procedure, for asymmetrical error estimates it involves a nonlinear optimization procedure. For most frequencies there is little difference between the power spectra derived from analyses of symmetrized error estimates and from asymmetrical error estimates, but this is not the case for the principal peak in the power spectrum at 9.43 yr ⁻¹. A likelihood analysis that takes account of the error asymmetry leads to a peak with power 13.24 at that frequency, and a Monte Carlo analysis shows that there is a chance of only 0.1% of finding a peak this big or bigger in the search band 1 – 36 yr ⁻¹. From this perspective, power-spectrum analysis that takes account of asymmetry of the error estimates gives evidence for variability that is significant at the 99.9% level. We comment briefly on an apparent discrepancy between power-spectrum analyses of the Super-Kamiokande and SNO solar neutrino experiments.     

measurement of noise power spectrum in photon counting PMT's and the effect of the observed (1/f) noise on astronomical photometry

In astronomical photometry, the sensitivity of observations is limited by the dark counts of the photomultiplier tube. In the present work, the effect of dark count noise in photon counting systems is investigated by theory and experimental measurements. Dark counts are considered to be originating from two sources, namely: dc fluctuations and random pulses.Experimental measurements were carried out to determine noise effects in different operating regions of noise dominance. The results provide strong evidence that: in normal operating mode, where the effect of random pulses is dominant, dark counts do not follow Poisson statistics. The observed noise shows strong (1/f) power spectrum, where the observed noise power is found to increase with time of observation.The results are important in photon counting systems operating under dark count limited mode. The conclusions drawn can be useful in obtaining more accurate error estimates and in assessing astronomical photometric observations and data reduction techniques.     

The visibility of low‐frequency solar acoustic modes

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. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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     

Periodicity of North–South asymmetry of sunspot area revisited: Cepstrum analysis

The North–South asymmetry in solar activity is an important solar property. During recent years the presence of short and intermediate-term periodicities of the North–South asymmetry has been reported, implying it could be an observational constraint that should be explained by the solar dynamo action. We attempt to confirm reported periodicities by employing an independent method so that observational constraints are established. We adopt the cepstrum technique to deconvolve the signal and the noise so that spurious peaks due to random noise are eliminated in the Fourier power spectrum. We also demonstrate how effectively it removes random noise using artificial data generated by assuming that the asymmetry of the sunspot area is characterized by random noise superposed on a slowly varying sinusoidal background function. We find that (i) the main periodicity of the North–South difference corresponding to ～9 years is present, (ii) other periodicities recently claimed, such as, of ～1.4, ～3.8, ～43 years, can not survive the deconvolution process, suggesting they seem due to stochastic random noise. We conclude by pointing out a possible implication of our finding on the nature of the North–South asymmetry of solar activity.     

Variations of the low‐degree solar p‐mode frequencies for 10 years of GOLF observations

We experiment with a method of measuring the frequency of solar p modes, intended to extend the passband for the variations of the frequency spectrum as high as possible. So far this passband is limited to a fraction of μ Hz for the classical analysis based on numerical fits of a theoretical line profile to a power spectrum averaged over periods lasting at least several weeks. This limit for the present analysis can be shifted to the mHz range, corresponding to some of the “5 min” oscillations, but in this range we use a lower resolution which allows us to separate odd and even p modes. We show an example of the results for long term variations and apply this analysis to search for a modulation of the p‐mode frequency spectrum by asymptotic series of solar g modes. A faint signal is found in the analysis of 10 years of GOLF data. This very preliminary result possibly indicates the detection of a small number of g modes of degree l = 1. A tentative determination of an observational value of the parameter P₀ follows. P₀ is the scaling factor of the asymptotic series of g modes and is a key data for solar core physics. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Angular power spectrum of CMB anisotropy from WMAP

The remarkable improvement in the estimates of different cosmological parameters in recent years has been largely spearheaded by accurate measurements of the angular power spectrum of cosmic microwave background (CMB) radiation. This has required removal of foreground contamination as well as detector noise bias with reliability and precision. Recently, a novel model-independent method for the estimation of CMB angular power spectrum from multi-frequency observations has been proposed and implemented on the first year WMAP (WMAP-1) data by Saha et al. [Saha, R., Jain, P., Souradeep, T., 2006. ApJL, 645, L89]. We review the results from WMAP-1 and also present the new angular power spectrum based on three years of the WMAP data (WMAP-3). Previous estimates have depended on foreground templates built using extraneous observational input to remove foreground contamination. This is the first demonstration that the CMB angular spectrum can be reliably estimated with precision from a self contained analysis of the WMAP data. The primary product of WMAP are the observations of CMB in 10 independent difference assemblies (DA) distributed over five frequency bands that have uncorrelated noise. Our method utilizes maximum information available within WMAP data by linearly combining DA maps from different frequencies to remove foregrounds and estimating the power spectrum from the 24 cross-power spectra of clean maps that have independent noise. An important merit of the method is that the expected residual power from unresolved point sources is significantly tempered to a constant offset at large multipoles (in contrast to the l² contribution expected from a Poisson distribution) leading to a small correction at large multipoles. Hence, the power spectrum estimates are less susceptible to uncertainties in the model of point sources.     

Effects of destriping errors on estimates of the CMB power spectrum

Destriping methods for constructing maps of the cosmic microwave background (CMB) anisotropies have been investigated extensively in the literature. However, their error properties have been studied in less detail. Here we present an analysis of the effects of destriping errors on CMB power spectrum estimates for Planck -like scanning strategies. Analytic formulae are derived for certain simple scanning geometries that can be rescaled to account for different detector noise. Assuming Planck -like low-frequency noise, the noise power spectrum is accurately white at high multipoles  (ℓ≳ 50)  . Destriping errors, though dominant at lower multipoles, are small in comparison to the cosmic variance. These results show that simple destriping map-making methods should be perfectly adequate for the analysis of Planck data and support the arguments given in an earlier paper in favour of applying a fast hybrid power spectrum estimator to CMB data with realistic '1/ f ' noise.     

基于盲退卷积的AIA图像增强方法

太阳图像中存在各种不同尺度、亮度和结构的物理活动现象,由于太阳日冕高动态活动和传感器设备等因素的影响,太阳图像成像质量不佳。根据太阳动力学天文台(Solar Dynamic Observatory,SDO)的大气成像仪(Atmospheric Imaging Assenbly,AIA)拍摄不同波段数据结构的动态范围大、噪声大、结构相对模糊等特点,提出一种基于盲退卷积的图像增强方法。首先对图像进行去噪和降低动态范围的处理,基于图像功率谱的分布假设,从原图中估计点扩散函数(Point Spread Function,PSF)的功率谱;然后使用相位提取算法恢复点扩散函数的相位,再退卷积得出较高质量的目标图像;最后通过轮廓切片分析、功率谱分析以及点扩散函数分析对增强结果进行定量和定性评价。实验结果表明,相比现有的图像增强方法,该方法在有效增强太阳日冕图像细节结构的同时,能够复原原图中因模糊无法识别的结构。     

Inverting the angular correlation function

The two point angular correlation function is an excellent measure of structure in the Universe. To extract from it the three-dimensional power spectrum, one must invert Limber's equation. Here we perform this inversion using a Bayesian prior constraining the smoothness of the power spectrum. Among other virtues, this technique allows for the possibility that the estimates of the angular correlation function are correlated from bin to bin. The outputs of this technique are estimators for the binned power spectrum and a full covariance matrix. Angular correlations mix small and large scales but after the inversion, small-scale data can be trivially eliminated, thereby allowing for realistic constraints on theories of large-scale structure. We analyse the automated plate measurement (APM) catalogue as an example, comparing our results with previous results. As a by-product of these tests, we find – in rough agreement with previous work – that APM places stringent constraints on cold dark matter inspired models, with the shape parameter constrained to be 0.25±0.04 (using data with wavenumber k ≤0.1  h  Mpc⁻¹). This range of allowed values uses the full power spectrum covariance matrix, but assumes negligible covariance in the off-diagonal angular correlation error matrix, which is estimated with a large angular resolution of 0.5° (in the range 0.5° and 20°).     

Uncovering the Bias in Low-Degree p-Mode Linewidth Fitting

Obtaining reliable estimates of linewidths in the power spectra of low-degree p modes is problematic at low frequency. In this regime, the mode coherence time increases with decreasing frequency, often causing the modes to be unresolved in relatively long duration spectra. The signal-to-noise ratio is also less favourable at low frequency, resulting in fits to power spectra underestimating the true linewidth of the p modes owing to the tails of the Lorentzian peaks becoming dominated by the background noise. We use a numerical simulation approach to assess the effect of this bias on the fitted widths of p-mode peaks and calculate observational duration limits required to obtain an unbiased estimate of the p-mode linewidth as a function of frequency. This is done in four different cases, where the precision of the artificial data is set at 0.25, 0.50, 0.75, and 1.00 m?s^?1 by adding random?scatter to increase the sample standard deviation per 40-second measurement. In all cases, the observational duration required to accurately obtain width estimates increases beyond that required for sufficient spectral resolution below a certain threshold frequency. For modes at ≈?1500 μHz, with an amplitude of approximately ten times the background, observations of up to 972 days are required to obtain an unbiased estimate of the linewidth. This is equivalent to ≈?18 times the coherence time of the corresponding p modes.     

Helioseismic Constraints on Rotation and Magnetic Fields in the Solar Core

In recent years, more and more precise measurements have been made of solar oscillation frequencies and line widths. From space, the Solar and Heliospheric Observatory/Michelson Doppler Imager (MDI) data has led to much progress. From the ground, networks, like Global Oscillation Network Group (GONG), Taiwanese Oscillation Network (TON), and Birmingham Solar Oscillations Network (BiSON) have also led to much progress. The sharpened and enriched oscillation spectrum of data have been critically complemented by advances in the treatments of the opacities and the equation of state. All of this has led to a significantly more precise probing of the solar core. Here we discuss the progress made and suggest how the core may be better probed with seismic data on-hand. In particular, we review our knowledge of the rotation and structure of the core. We further argue that much may be learned about the core by exploiting the line width data from the aforementioned sources. Line-width data can be used to place sharper constraints on core properties, like the degree to which the Sun rotates on a single axis and the upper limit on magnetic fields that may be buried in the core.     

Combined Analysis of Solar Neutrino and Solar Irradiance Data: Further Evidence for Variability of the Solar Neutrino Flux and Its Implications Concerning the Solar Core

A search for any particular feature in any single solar neutrino dataset is unlikely to establish variability of the solar neutrino flux since the count rates are very low. It helps to combine datasets, and in this article we examine data from both the Homestake and GALLEX experiments. These show evidence of modulation with a frequency of 11.85 year⁻¹, which could be indicative of rotational modulation originating in the solar core. We find that precisely the same frequency is prominent in power spectrum analyses of the ACRIM irradiance data for both the Homestake and GALLEX time intervals. These results suggest that the solar core is inhomogeneous and rotates with a sidereal frequency of 12.85 year⁻¹. From Monte Carlo calculations, it is found that the probability that the neutrino data would by chance match the irradiance data in this way is only 2 parts in 10 000. This rotation rate is significantly lower than that of the inner radiative zone (13.97 year⁻¹) as recently inferred from analysis of Super-Kamiokande data, suggesting that there may be a second, inner tachocline separating the core from the radiative zone. This opens up the possibility that there may be an inner dynamo that could produce a strong internal magnetic field and a second solar cycle.     

General presentation of a single IRIS site raw data analysis problem

A complete software package has been built for the calibration in m s ^–1 of the velocity residuals due to solar oscillations in the raw IRIS (International Research on the Interior of the Sun) data. It takes into account all known astronomical components contributing to the line-of-sight velocity between the instrument and the solar surface, and also the apparent velocity due to the non-uniform integration of the solar rotation as seen through an inhomogeneous Earth atmosphere. The IRIS data itself is used for the estimation of the nonlinear instrumental response to the velocity, and the residual can be directly obtained in velocity units, without low frequency filtering. On a day of typical photometric sky quality, the power spectrum obtained appears to be solar noise limited.     

Inhomogeneity in the solar core

In recent years, normal-mode helioseismology has shown that the spherically averaged sound-speed distribution throughout the solar interior is in remarkable agreement with suitable standard solar models. This implies that any deviation of the theoretical models from the Sun has only a very small influence on the oscillation frequency spectrum (excluding the contributions from the uncertain near-surface layers). Nevertheless, it is important to determine whether the Sun really is very similar to a standard model, or whether there are substantial differences. This is especially important of the energy-generating core, particularly because it is likely to be necessary to understand the conditions under which the nuclear reactions are taking place in order to utilize neutrino detectors to the full to measure the properties of neutrino transitions.     

Cluster formation and the Sunyaev–Zel'dovich power spectrum in modified gravity: the case of a phenomenologically extended DGP model

We investigate the effect of modified gravity on cluster abundance and the Sunyaev–Zel'dovich (SZ) angular power spectrum. Our modified gravity is based on a phenomenological extension of the Dvali–Gabadadze–Porrati model which includes two free parameters characterizing deviation from Λ cold dark matter cosmology. Assuming that Birkhoff's theorem gives a reasonable approximation, we study the spherical collapse model of structure formation and show that while the growth function changes to some extent, modified gravity gives rise to no significant change in the linear density contrast at collapse time. The growth function is enhanced in the so called normal branch, while in the 'self-accelerating' branch it is suppressed. The SZ angular power spectrum is computed in the normal branch, which allows us to put observational constraints on the parameters of the modified gravity model using small scale cosmic microwave background observation data.     

Phase differences between irradiance and velocity low degree solar acoustic modes revisited

Since 1984, simultaneous observations of irradiance and velocity solar acoustic modes, have been carried out by several authors in order to measure the phase difference between irradiance and velocity modes. Following the earliest observations with stratospheric balloon (Frolich and van Der Raay, 1984), a two ground-based stations (Tenerife and Baja California) were established (Jimenez et al, 1990) obtaining coherence results in the frequency range from 2.5 mHz to 4.3 mHz. These phase differences between irradiance and velocity solar acoustic modes are interpreted in terms of the non-adiabatic behaviour of the solar atmosphere. In 1988 the IPHIR (Frolich et al, 1988) instrument flown on the PHOBOS-2 mission to Mars and measured the solar irradiance during 150 consecutive days. The best velocity observations obtained in Tenerife for this period were compared with IPHIR data to compute the phase differences (Schrijver et al, 1991). The final conclusion is that good agreement is attained between space quadsi-space and ground observations which yield a phase diffenrece of about -125 degrees in the frequency range 2.5 mHz to 4.2 mHz, with a slight increase suggested by the data running up to 4.6 mHz.