Anisotropy and Dissipation of Turbulence and Their Effects on Solar Models

1 INTRODUCTIONThe maing-length theory (MLT) is the most commonly used approach to calculate convective energy transport in stars and other astrophysical situations. Based on the original idea ofPrandtl (1952) that turbulent parcels trallsfer heat in a similar way as molecules of gas do inthermal conduction, the MLT assumes that convection cells, drived by buoyancy, move thlougha ~ng length 1 and release the heat they carry when they merge with their environment. Themost widely adopted f…     

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.     

Effects of errors in the solar radius on helioseismic inferences

Frequencies of intermediate-degree f modes of the Sun seem to indicate that the solar radius is smaller than what is normally used in constructing solar models. We investigate the possible consequences of an error in radius on results for solar structure obtained using helioseismic inversions. It is shown that solar sound speed will be overestimated if oscillation frequencies are inverted using reference models with a larger radius. Using solar models with a radius of 695.78 Mm and new data sets, the base of the solar convection zone is estimated to be at a radial distance of 0.7135 ± 0.0005 of the solar radius. The helium abundance in the convection zone as determined using models with an OPAL equation of state is 0.248 ± 0.001, where the errors reflect the estimated systematic errors in the calculation, the statistical errors being much smaller. Assuming that the OPAL opacities used in the construction of the solar models are correct, the surface Z / X is estimated to be 0.0245 ± 0.0006.     

Seismology of the solar convection zone

An attempt is made to infer the structure of the solar convection zone from observedp-mode frequencies of solar oscillations. The differential asymptotic inversion technique is used to find the sound speed in the solar envelope. It is found that envelope models which use the Canuto-Mazzitelli (CM) formulation for calculating the convective flux give significantly better agreement with observations than models constructed using the mixing length formalism. This inference can be drawn from both the scaled frequency differences and the sound speed difference. The sound speed in the CM envelope model is within 0.2% of that in the Sun except in the region withr > 0.99R _⊙. The envelope models are extended below the convection zone, to find some evidence for the gravitational settling of helium beneath the base of the convection zone. It turns out that for models with a steep composition gradient below the convection zone, the convection zone depth has to be increased by about 6 Mm in order to get agreement with helioseismic observations.     

Seismology of stellar envelopes: probing the outer layers of a star through the scattering of acoustic waves

The outer layers of Sun-like stars are regions of rapid spatial variation which modulate the p-mode frequencies by partially reflecting the constituent acoustic waves. With the accuracy that has been achieved by current solar observations, and that is expected from imminent stellar observations, this modulation can be observed from the spectra of the low-degree modes. We present a new and simple theoretical calculation to determine the leading terms in an asymptotic expansion of the outer phase of these modes, which is determined by the structure of the surface layers of the star. Our procedure is to compare the stellar envelope with a plane-parallel polytropic envelope, which we regard as a smooth reference background state. Then we can isolate a seismic signature of the acoustic phase and relate it to the stratification of the outer layers of the convection zone. One can thereby constrain theories of convection that are used to construct the convection zones of the Sun and Sun-like stars. The accuracy of the diagnostic is tested in the solar case by comparing the predicted outer phase with an exact numerical calculation.     

3-D non-linear evolution of a magnetic flux tube in a spherical shell: The isentropic case

We present recent 3-D MHD numerical simulations of the non-linear dynamical evolution of magnetic flux tubes in an adiabatically stratified convection zone in spherical geometry, using the anelastic spherical harmonic (ASH) code.We seek to understand the mechanism of emergence of strong toroidal fields from the base of the solar convection zone to the solar surface as active regions. We confirm the results obtained in cartesian geometry that flux tubes that are not twisted split into two counter vortices before reaching the top of the convection zone. Moreover, we find that twisted tubes undergo the poleward-slip instability due to an unbalanced magnetic curvature force which gives the tube a poleward motion both in the non-rotating and in the rotating case. This poleward drift is found to be more pronounced on tubes originally located at high latitudes. Finally, rotation is found to decrease the rise velocity of the flux tubes through the convection zone, especially when the tube is introduced at low latitudes. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Solar internal rotation rate and the latitudinal variation of the tachocline

A new set of accurately measured frequencies of solar oscillations is used to infer the rotation rate inside the Sun, as a function of radial distance as well as latitude. We have adopted a regularized least-squares technique with iterative refinement for both 1.5D inversion, using the splitting coefficients, and 2D inversion using individual m splittings. The inferred rotation rate agrees well with earlier estimates showing a shear layer just below the surface and another one around the base of the convection zone. The tachocline or the transition layer where the rotation rate changes from differential rotation in the convection zone to an almost latitudinally independent rotation rate in the radiative interior is studied in detail. No compelling evidence for any latitudinal variation in the position and width of the tachocline is found, although it appears that the tachocline probably shifts to a slightly larger radial distance at higher latitudes and possibly also becomes thicker. However, these variations are within the estimated errors and more accurate data would be needed to make a definitive statement about latitudinal variations.     

Distribution of magnetic flux on the solar surface and low-degree p-modes

The frequencies of solar p-modes are known to change over the solar cycle. There is also recent evidence that the relation between frequency shift of low-degree modes and magnetic flux or other activity indicators differs between the rising and falling phases of the solar cycle, leading to a hysteresis in such diagrams. We consider the influence of the changing large-scale surface distribution of the magnetic flux on low-degree ( l ≤3) p-mode frequencies. To that end, we use time-dependent models of the magnetic flux distribution and study the ensuing frequency shifts of modes with different order and degree as a function of time. The resulting curves are periodic functions (in simple cases just sine curves) shifted in time by different amounts for the different modes. We show how this may easily lead to hysteresis cycles comparable to those observed. Our models suggest that high-latitude fields are necessary to produce a significant difference in hysteresis between odd- and even-degree modes. Only magnetic field distributions within a small parameter range are consistent with the observations by Jiménez-Reyes et al. Observations of p-mode frequency shifts are therefore capable of providing an additional diagnostic of the magnetic field near the solar poles. The magnetic distribution that is consistent with the p-mode observations also appears reasonable compared with direct measurements of the magnetic field.     

Use of genetic algorithms to determine low-l rotational p-mode splittings at high frequencies

We present a comparative study of genetic and standard fitting routines applied to the task of extracting reliable estimates of the rotational splitting of full-disc, low-angular-degree (low- l ) solar p-mode data at high frequencies. 100 artificial proxies of a 10-yr data base of observations made by the Birmingham Solar-Oscillations Network (BiSON) were used to test the two approaches. All sets were analysed over the frequency range from 3000 to 4000 μHz.
Previous work, based on non-linear 'hill-climbing' fitting techniques, has demonstrated the unfortunate tendency for full-disc estimates of the splitting to overestimate the true, underlying values at high frequencies. Here, we show that the resulting bias is less severe when a genetic-fitting approach is adopted. This is largely the result of the number of erroneous 'null-valued' estimates of the splitting being considerably reduced: these estimates are, in effect, re-introduced into the expected normal distribution of fitted splittings. We also illustrate the diverse control one has when using a genetic algorithm as a fitting routine; this diversity is shown to allow further refinement in the estimate of the rotational splitting.
Finally, we address the issue of the reliability of the formal splitting uncertainties returned by the mode fitting, and find that complications arising from the strong anticorrelation between the splittings and their error estimates are not alleviated by the use of the genetic technique.     

Solar p modes of high degree l: coupling by differential rotation

Accurate measurements of solar p-mode frequencies and frequency splittings at high degree l require an adequate theoretical knowledge of the effects of mode coupling, induced by the variation with latitude of the angular velocity of the solar internal rotation. Earlier results for expansion coefficients of composite solutions (coupling coefficients) are due to Woodard. In this paper, the analysis is extended to allow for the dependence of the differential rotation on depth, and the result is expressed in terms of measurable quantities (the rotational splitting coefficients), which makes it convenient for diagnostic purposes. The analysis is based on the approach of quasi-degenerate perturbation theory, and is extended further to address possible effects of mode coupling in the observational line profiles. It is shown, using approximations applicable at high degree l , that the expected line profiles of composite modes in the observational power spectra are not distorted by mode coupling.     

On the accuracy of the measurements of low-ℓ solar acoustic oscillations

Using a 1154 d long measurement of solar oscillations, obtained by the Global Oscillation Network Group from 1995 June 10 to 1998 August 6, we study the dependence of the accuracy of radial p-mode parameters on the duration of the observations. It is shown that relatively rare pulses of large power lead to the decrease of the accuracy achievable for a given duration of the observations and it is usually underestimated. The corresponding correction factor to the Libbrecht formula for a frequency accuracy estimation is provided. We have also investigated the influence of the solar activity on the mode parameters soon after the solar activity minimum. There is a clearly visible increase of the radial p-mode power in the beginning of the new solar cycle while the mode frequency variations are within the corresponding error bars.     

Acoustic interferometry of the solar atmosphere: p-modes with frequencies near the 'acoustic cut-off'

High-frequency p-mode intensity data, obtained from the South Pole in 1987, 1988, 1990 and 1994, show a sharp variation in the phase-shift function and in the frequency spacings near 5.5 mHz. Using a simple theoretical model, we demonstrate that this behaviour is caused by an acoustic resonance in the atmosphere between the excitation source and the upper reflection level. We discuss the diagnostic properties of this resonance, which is sensitive to the acoustic reflectivity of the solar atmosphere and to the location and parity of the excitation source. When applied to the solar data, our model indicates that the average acoustic reflectivity increases with increasing solar activity. The model also shows that the acoustic source has composite parity and is located within one pressure scaleheight of the base of the photosphere.     

Turbulence properties in the solar convection envelope： properties of the overshooting regions

We apply the turbulent convection model (TCM) to investigate properties of tur-bulence in the solar convective envelope, especially in overshooting regions. The results show TCM gives negative turbulent heat flux uγ′T′in overshooting regions, which is sim-ilar to other nonlocal turbulent convection theories. The turbulent temperature fluctuation T′T′shows peaks in overshooting regions. Most important, we find that the downward overshooting region below the base of the solar convection zone is a thin cellular layer filled with roll-shaped convective cells. The overshooting length for the temperature gradi-ent is much shorter than that for element mixing because turbulent heat flux of downward and upward moving convective cells counteract each other in this cellular overshooting region. Comparing the models' sound speed with observations, we find that raking the convective overshooting into account helps to improve the sound speed profile of our nonlocal solar models. Comparing the p-mode oscillation frequencies with observations,we validated that increasing the diffusion parameters and decreasing the dissipation pa-rameters of TCM make the p-mode oscillation frequencies of the solar model be in betteragreement with observations.     

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.     

Determining the rotation law of the solar convection zone

The condition of minimum total dissipation is used to derive stationary rotation and azimuthal magnetic field distributions in the bulk of the solar convection zone with an upper boundary at which the relative radius is r/R=0.95. General equilibrium con figurations with symmetric and antisymmetric (about the equator) angular-velocity and field components are determined. The calculated rotation law matches the observed one in general parameters, but the decrease in angular velocity at high latitudes in theory is larger than that in observations. Besides, there are additional sharp variations in the rotation and field distributions in the theoretical curves near the generation zone of solar torsional waves. The possible cause of the latter discrepancy is discussed. The change in equilibrium distributions due to the presence of an inverse molecular-weight gradient at the base of the convection zone is also studied. This gradient is known to be produced by accelerated gravitational helium settling in the convection zone.     

Sounding stellar cycles with Kepler – I. Strategy for selecting targets

The long-term monitoring and high photometric precision of the Kepler satellite will provide a unique opportunity to sound the stellar cycles of many solar-type stars using asteroseismology. This can be achieved by studying periodic changes in the amplitudes and frequencies of the oscillation modes observed in these stars. By comparing these measurements with conventional ground-based chromospheric activity indices, we can improve our understanding of the relationship between chromospheric changes and those taking place deep in the interior throughout the stellar activity cycle. In addition, asteroseismic measurements of the convection zone depth and differential rotation may help us determine whether stellar cycles are driven at the top or at the base of the convection zone. In this paper, we analyse the precision that will be possible using Kepler to measure stellar cycles, convection zone depths and differential rotation. Based on this analysis, we describe a strategy for selecting specific targets to be observed by the Kepler Asteroseismic Investigation for the full length of the mission, to optimize their suitability for probing stellar cycles in a wide variety of solar-type stars.     

On the choice of parameters in solar-structure inversion

The observed solar p-mode frequencies provide a powerful diagnostic of the internal structure of the Sun and permit us to test in considerable detail the physics used in the theory of stellar structure. Among the most commonly used techniques for inverting such helioseismic data are two implementations of the optimally localized averages (OLA) method, namely the subtractive optimally localized averages (SOLA) and multiplicative optimally localized averages (MOLA). Both are controlled by a number of parameters, the proper choice of which is very important for a reliable inference of the solar internal structure. Here we make a detailed analysis of the influence of each parameter on the solution and indicate how to arrive at an optimal set of parameters for a given data set.     

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.     

Probing solar convection

In the solar convection zone, acoustic waves are scattered by turbulent sound speed fluctuations. In this paper the scattering of waves by convective cells is treated using Rytov's technique. Particular care is taken to include diffraction effects, which are important, especially for high-degree modes that are confined to the surface layers of the Sun. The scattering leads to damping of the waves and causes a phase shift. Damping manifests itself in the width of the spectral peak of p-mode eigenfrequencies. The contribution of scattering to the linewidths is estimated and the sensitivity of the results to the assumed spectrum of the turbulence is studied. Finally, the theoretical predictions are compared with recently measured linewidths of high-degree modes.     

A study of possible temporal and latitudinal variations in the properties of the solar tachocline

Temporal variations of the structure and the rotation rate of the solar tachocline region are studied using helioseismic data from the Global Oscillation Network Group (GONG) and the Michelson Doppler Imager (MDI) obtained during the period 1995–2000. We do not find any significant temporal variation in the depth of the convection zone, the position of the tachocline or the extent of overshoot below the convection zone. No systematic variation in any other properties of the tachocline, like width, etc., is found either. The possibility of periodic variations in these properties is also investigated. Time-averaged results show that the tachocline is prolate with a variation of about 0.02 R_⊙ in its position. Neither the depth of the convection zone nor the extent of overshoot shows any significant variation with latitude.