Comparison OF Noise Properties Of GONG And MDI Time-Distance Helioseismic Data

The recently upgraded system of the ground-based Global Oscillation Network Group (GONG) network of helioseismic observatories has started to provide higher-resolution solar oscillation measurements suitable for local helioseismic studies. Selecting simultaneously observed regions on the Sun by both GONG and the space-borne Michelson Doppler Imager (MDI) instrument on board the Solar and Heliospheric Observatory (SOHO), we perform a comparative analysis of time-distance measurements focussing on the noise properties.     

Acoustic Radius Measurements from MDI and GONG

We study the temporal autocorrelation function (ACF) of global solar oscillations. It is well known that the “large frequency separation” is proportional to the solar acoustic radius. We analyze the ACF of MDI and GONG spherical-harmonic-coefficient time series for degrees ?=0?3. Acoustic radius measurements obtained from the first dominant peak locations of the ACF show a significant anticorrelation with solar cycle. This technique can be a useful tool to search for stellar activity.     

Non-adiabatic effects in the solar p-mode oscillations

Non-adiabatic effects associated with radiative loss and convective transfer add small imaginary parts to solar p-mode eigenfrequencies. An asymptotic approximation is developed to study the non-adiabatic effects. For the outermost layer where the approximation is not valid, an exact solution of the equation of non-adiabatic oscillation which results in an cigenfunction equation, is used. The inclusion of the non-adiabatic effects reduces the discrepancy between the theory and observations particularly for frequencies above 3.1 mHz.     

On the measurement precision of solar p-mode eigenfrequencies

We make use of 3456 d of observations of the low-ℓ p-mode oscillations of the Sun in order to study the evolution over time of the measurement precision of the radial eigenfrequencies. These data were collected by the ground-based Birmingham Solar-Oscillations Network (BiSON) between 1991 January and 2000 June. When the power spectrum of the complete time series is fitted, the analysis yields frequency uncertainties that are close to those expected from the returned coherence times of the modes. The slightly elevated levels compared with the prediction appear to be consistent with a degradation of the signal-to-noise ratio in the spectrum that is the result of the influence of the window function of the observations (duty cycle 71 per cent). The fractional frequency precision reaches levels of a several parts in 10⁶ for many of the modes. The corresponding errors reported from observations made by the GOLF instrument on board the ESA/NASA SOHO satellite, when extrapolated to the length of the BiSON data set, are shown to be (on average) about ∼25 per cent smaller than their BiSON counterparts owing to the uninterrupted nature of the data from which they were derived.
An analysis of the BiSON data in contiguous segments of different lengths, T , demonstrates that the frequency uncertainties scale as T ^−1/2. This is to be expected in the regime where the coherence (life) times of the modes, τ _{n ℓ}, are smaller than the observing time T (the 'oversampled' regime). We show that mode detections are only now beginning to encroach on the 'undersampled' regime (where   T < τ _{n ℓ})  .     

Intermediate degree p-mode frequency splittings near solar maximum

We report on observations of global solar Ca K-line intensity oscillations taken in May 1991 from Mees Solar Observatory, Hawaii. We measurep-mode frequency splittings for modes of spherical harmonic degrees between 20 and 129 averaged over the radial order of the modes. Our measurement of the antisymmetric component of the splittings is comparable with previous measurements and thus indicates a decrease in the latitudinal differential rotation with depth into the convection zone and the upper radiative zone. We find evidence for a 1% variation in the rotation rate of the upper convection zone roughly in phase with the solar activity cycle. Our measurement of the symmetric component of the splittings is of the same order as was reported from the previous solar maximum and is an order of magnitude larger than has been measured near solar minimum. From the degree dependence of the symmetric component of the splittings, we find evidence for an aspherical fractional sound speed perturbation located at a depth of 0.85 ± 0.05 solar radii. This perturbation has a magnitude ofc/c +9 × 10^–4 at the equator relative to the poles. Additionally, there is evidence for a near-surface aspherical sound speed perturbation of smaller magnitudec/c +4 × 10^–4 at the equator relative to the poles. If an intense global magnetic field were the dominant source of the observed symmetric component of the splittings, instead of latitudinal gradients in the sound speed, then global fields of order 10⁵ G would be required in the convection zone.     

Solar-cycle variation of the sound-speed asphericity from GONG and MDI data 1995–2000

We study the variation of the frequency splitting coefficients describing the solar asphericity in both GONG and MDI data, and use these data to investigate temporal sound-speed variations as a function of both depth and latitude during the period 1995–2000 and a little beyond. The temporal variations in even splitting coefficients are found to be correlated to the corresponding component of magnetic flux at the solar surface. We confirm that the sound-speed variations associated with the surface magnetic field are superficial. Temporally averaged results show a significant excess in sound speed around     and latitude of 60°.     

Empirical estimate of p-mode frequency shift for solar cycle 23

We have obtained empirical relations between the p-mode frequency shift and the change in solar activity indices. The empirical relations are determined on the basis of frequencies obtained from BBSO and GONG stations during solar cycle 22. These relations are applied to estimate the change in mean frequency for the cycle 21 and 23. A remarkable agreement between the calculated and observed frequency shifts for the ascending phase of cycle 23, indicates that the derived relations are independent of epoch and do not change significantly from cycle to cycle. We propose that these relations could be used to estimate the shift in p-mode frequencies for past, present and future solar activity cycles, if the solar activity index is known. The maximum frequency shift for cycle 23 is estimated to be 265±90 nHz, corresponding to a predicted maximum smoothed sunspot number 118.1±35.     

On the effect of magnetic field on the low-l solar p-mode oscillations

One of the possible magnetic field effects on the stellar pulsations is known to be a splitting in the observed frequencies. Using this knowledge in the solar convection zone, there are two aims in this work Considering the Sun as an incompressible fluid, our first objective was to investigate the variation of the physical parameters in the 30% outermost convective solar layer, during a pulsation period. The second purpose was to calculate, by means of the spherical harmonics, the shifts on the low-l p-mode frequencies which could be caused by the presence of the magnetic field in the Sun. The first order perturbation approximation was used in order to calculate analytically the resulting frequency shifts and the small perturbations on the magnetic field, as well as the physical parameters, such as density, pressure and temperature, of a Standard Solar Model excluding both rotation and magnetic field (Christensen-Dalsgaard et al., 1996) in the unperturbed equilibrium case.     

The influence of turbulence on the solar p-mode frequencies I. Free-mode approximation

The influence of turbulence on the frequencies of free acoustic modes in convection zones is considered. The frequencies are modified via the speed of sound by the turbulence-induced alterations of the effective pressure: (i) by the correlated fluctuations of temperature and density and (ii) the pressure part of the Reynolds stress. The two effects shift the frequency of low l p-modes in opposite directions. In addition, the correlation of the density fluctuations with the random velocity — the eddy-mass flow — is also relevant. It is, in a steady state, balanced by a vertical mean velocity. The balance results in a rather small net effect completely disappearing for highly nonradial oscillations. Both effects of the density fluctuations produce a redshift of the low l p-mode frequencies. The Reynolds stress, however, makes a blueshift of the frequencies relative to that computed for a laminar gas. This effect dominates for subsonic turbulences. The applied second-order correlation-approximation, however, only holds for the lowest frequencies, where the KORONAS (solar minimum) data are indicating a blueshift. Of particular importance for the present concept is the expected cycle-variations of the lineshifts, i.e. the consideration of the magnetic modification of the various contributions. Observations may show whether the suggested modifications of the solar oscillation theory are correct.     

Recent observations of high-degree solar p-mode oscillations at the Kitt Peak National Observatory

A brief summary is given of a program which is currently being carried out with the McMath telescope of the Kitt Peak National Observatory in order to study high-degree (l ≳ 150) solar p-mode oscillations. This program uses a 244 × 248 pixel CID camera and the main spectrograph of the McMath telescope to obtain velocity-time maps of the oscillations which can be converted into two-dimensional (k _h - ω) power spectra of the oscillations. Several different regions of the solar spectrum have been used in order to study the oscillations at different elevations in the solar atmosphere. The program concentrates on eastward- and westward-propagating sectoral harmonic waves so that measurements can be made of the absolute rotational velocities of the solar photospheric and shallow sub-photospheric layers. Some preliminary results from this program are now available. First, we have been unable to confirm the existence of a radial gradient in the equatorial rotational velocity as was previously suggested. Second, we have indeed been able to confirm the presence of p-mode waves in the solar chromosphere as was first suggested by Rhodes et al. (1977). Third, we have been able to demonstrate differences in photospheric and chromospheric power spectra.

     

Dynamics of the solar magnetic field from SOHO/MDI

The investigation of the dynamics of magnetic fields from small scales to the large scales is very important for the understanding of the nature of solar activity. It is also the base for producing adequate models of the solar cycle with the purpose to predict the level of solar activity. Since December 1995 the Michelson Doppler Imager (MDI) on board of the Solar and Heliospheric Observatory (SOHO) provides full disk magnetograms and synoptic maps which cover the period of solar cycle 23 and the current minimum. In this paper, I review the following important topics with a focus on the dynamics of the solar magnetic field. The synoptic structure of the solar cycle; the birth of the solar cycle (overlapping cycles 23 and 24); the relationship of the photospheric magnetic activity and the EUV solar corona, polar magnetic fields and dynamo theory (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A comparison of the use of genetic and hill-climbing algorithms to fit low-l solar p-mode spectra

We present a comparative study of low-l solar p-mode parameters extracted by genetic-algorithm and ‘standard’ hill-climbing minimisation routines. To effect this we make use of observations made in integrated sunlight by the Birmingham Solar-Oscillations Network (BiSON) and the GOLF instrument on board the ESA/NASA SOHO satellite, in addition to artificial data. We find that over the central part of the p-mode range the two fitting routines return similar results. However, at low frequencies — where the S/N in the modes is low and their resonant peaks narrow — we find that the genetic routine appears to offer more robust estimates of the underlying parameters.     

Changes in convective properties over the solar cycle: effect on p-mode damping rates

Measurements of both solar irradiance and p-mode oscillation frequencies indicate that the structure of the Sun changes with the solar cycle. Balmforth, Gough & Merryfield investigated the effect of symmetrical thermal disturbances on the solar structure and the resulting pulsation frequency changes. They concluded that thermal perturbations alone cannot account for the variations in both irradiance and p-mode frequencies, and that the presence of a magnetic field affecting acoustical propagation is the most likely explanation of the frequency change, in the manner suggested earlier by Gough & Thompson and by Goldreich et al. Numerical simulations of Boussinesq convection in a magnetic field have shown that at high Rayleigh number the magnetic field can modify the preferred horizontal length scale of the convective flow.
Here, we investigate the effect of changing the horizontal length scale of convective eddies on the linewidths of the acoustic resonant mode peaks observed in helioseismic power spectra. The turbulent fluxes in these model computations are obtained from a time-dependent, non-local generalization of the mixing-length formalism. The modelled variations are compared with p-mode linewidth changes revealed by the analysis of helioseismic data collected by the Birmingham Solar-Oscillations Network (BiSON); these low-degree (low- l ) observations cover the complete falling phase of solar activity cycle 22. The results are also discussed in the light of observations of solar-cycle variations of the horizontal size of granules and with results from 2D simulations by Steffen of convective granules.     

Microflares and thermal background of the solar corona

We investigate low-intensity microflares in the soft component of the solar X-ray radiation over the period from September through December 1995 within the framework of the Interball—Geotail project. We derived the intensity distribution of microflares and found correlations between the daily mean peak fluxes of X-ray bursts from microflares of various classes and the daily mean values of the thermal background of the solar corona.     

The subsurface radial gradient of solar angular velocity from MDI f-mode observations

We report quantitative analysis of the radial gradient of solar angular velocity at depths down to about 15 Mm below the solar surface for latitudes up to 75° using the Michelson Doppler Imager (MDI) observations of surface gravity waves (fmodes) from the Solar and Heliospheric Observatory (SOHO). A negative outward gradient of around –400 nHz/R , equivalent to a logarithmic gradient of the rotation frequency with respect to radius which is very close to –1, is found to be remarkably constant between the equator and 30° latitude. Above 30° it decreases in absolute magnitude to a very small value at around 50°. At higher latitudes the gradient may reverse its sign: if so, this reversal takes place in a thin layer extending only 5 Mm beneath the visible surface, as evidenced by the most superficial modes (with degrees l>250). The signature of the torsional oscillations is seen in this layer, but no other significant temporal variations of the gradient and value of the rotation rate there are found.     

Determination of magnetic helicity content of solar active regions from SOHO/MDI magnetograms

Chae (2001) first proposed a method of self-consistently determining the rate of change of magnetic helicity using a time series of longitudinal magnetograms only, such as taken by SOHO/MDI. Assuming that magnetic fields in the photosphere are predominantly vertical, he determined the horizontal component of velocity by tracking the displacements of magnetic flux fragments using the technique of local correlation tracking (LCT). In the present paper, after briefly reviewing the recent advance in helicity rate measurement, we argue that the LCT method can be more generally applied even to regions of inclined magnetic fields. We also present some results obtained by applying the LCT method to the active region NOAA 10365 under emergence during the observable period, which are summarized as follows. (1) Strong shearing flows were found near the polarity inversion line that were very effective in helicity injection. (2) Both the magnetic flux and helicity of the active region steadily increased during the observing period, and reached 1.2 × 10²² Mx and 8 ×10⁴² Mx², respectively, 4.5 days after the birth of the active region. (3) The corresponding ratio of the helicity to the square of the magnetic flux, 0.05, is roughly compatible with the values determined by other studies using linear-force-free modeling. (4) A series of flares took place while the rate of helicity injection was high. (5) The choice of a smaller window size or a shorter time interval in the LCT method resulted in a bigger value of the LCT velocity and a bigger value of the temporal fluctuation of the helicity rate. (6) Nevertheless when averaged over a time period of about one hour or longer, the average rate of helicity became about the same within about 10%, almost irrespective of the chosen window size and time interval, indicating that short-lived, fluctuating flows may be insignificant in transferring magnetic helicity. Our results suggest that the LCT method may be applied to 96-minute cadence full-disk MDI magnetograms or other data of similar kind, to provide a practically useful, if not perfect, way of monitoring the magnetic helicity content of active regions as a function of time.     

Determination of magnetic helicity content of solar active regions from SOHO/MDI magnetograms

Chae (2001) first proposed a method of self-consistently determining the rate of change of magnetic helicity using a time series of longitudinal magnetograms only, such as taken by SOHO/MDI. Assuming that magnetic fields in the photosphere are predominantly vertical, he determined the horizontal component of velocity by tracking the displacements of magnetic flux fragments using the technique of local correlation tracking (LCT). In the present paper, after briefly reviewing the recent advance in helicity rate measurement, we argue that the LCT method can be more generally applied even to regions of inclined magnetic fields. We also present some results obtained by applying the LCT method to the active region NOAA 10365 under emergence during the observable period, which are summarized as follows. (1) Strong shearing flows were found near the polarity inversion line that were very effective in helicity injection. (2) Both the magnetic flux and helicity of the active region steadily increased during the observing period, and reached 1.2 × 10²² Mx and 8 ×10⁴² Mx², respectively, 4.5 days after the birth of the active region. (3) The corresponding ratio of the helicity to the square of the magnetic flux, 0.05, is roughly compatible with the values determined by other studies using linear-force-free modeling. (4) A series of flares took place while the rate of helicity injection was high. (5) The choice of a smaller window size or a shorter time interval in the LCT method resulted in a bigger value of the LCT velocity and a bigger value of the temporal fluctuation of the helicity rate. (6) Nevertheless when averaged over a time period of about one hour or longer, the average rate of helicity became about the same within about 10%, almost irrespective of the chosen window size and time interval, indicating that short-lived, fluctuating flows may be insignificant in transferring magnetic helicity. Our results suggest that the LCT method may be applied to 96-minute cadence full-disk MDI magnetograms or other data of similar kind, to provide a practically useful, if not perfect, way of monitoring the magnetic helicity content of active regions as a function of time.