Dynamical behavior of strong magnetic fields in the solar convection zone

Magnetic buoyancy is thought to play an important role in the dynamical behavior of the Sun's magnetic field in the convection zone. Magnetic buoyancy is commonly thought to cause inescapable rapid loss of toroidal flux from much of the convection zone, thereby suppressing effective operation of a solar dynamo. This paper re-examines the detailed character of magnetic buoyancy, especially as it is influenced by the magnetic field's effect on heat transport and temperature gradients in the convection zone. It is suggested that suppression of convective heat transport across strong magnetic flux tubes can alter the temperature within the tubes and can subdue, or even reverse, the effect of magnetic buoyancy.     

Large-scale flow and transport of magnetic flux in the solar convection zone

Horizontal large-scale velocity field describes horizontal displacement of the photospheric magnetic flux in zonal and meridian directions. The flow systems of solar plasma, constructed according to the velocity field, create the large-scale cellular-like patterns with up-flow in the center and the down-flow on the boundaries. Distribution of the largescale horizontal eddies (with characteristic scale length from 350 to 490 Mm) was found in the broad equatorial zone, limited by 60‡ latitude circles on both hemispheres. The zonal averages of the zonal and meridian velocities, and the total horizontal velocity for each Carrington rotation during the activity cycles no. 21 and 22 varies during the 11-yr activity cycle. Plot of RMS values of total horizontal velocity is shifted about 1.6 years before the similarly shaped variation of the magnetic flux.     

Limits on the magnetic field strength at the base of the solar convection zone

Howard (1993) finds a relationship between the tilt angles of BMRs (Bipolar Magnetic Regions) and the separation between their leading and following polarities; the tilt angle increases with polarity separation. Here we present a more detailed analysis of this relationship and show that this effect constrains the strength of the magnetic field at the bottom of the convection zone to a value between 40 and 150 kG, which confirms the constraints put by D'Silva and Choudhuri (1993) based on Joy's law (the tilt-latitude relationship), through an entirely different approach.Operated by the Association of Universities for Research in Astronomy, Inc., under CooperativeOperated by the Association of Universities for Research in Astronomy, Inc., under Cooperative     

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.     

The decay of torque free toroidal magnetic fields confined beneath the solar convection zone

Calculation similar to those of Mestel and Moss (1983) are performed to investigate the decay of a toroidal field through configurations satisfying the torque free condition, imposed by the presence of a poloidal field of dipolar form confined beneath the solar convection zone. It is found that initially stable field configurations diffuse into unstable configurations on time-scales of order a few x 10⁸ yr. The results are similar to those of Tayler (1982) for a simpler field model without any dynamical constraints.     

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.     

Turbulent dynamo near tachocline and reconstruction of azimuthal magnetic field in the solar convection zone

In order to extend the abilities of the αΩ dynamo model to explain the observed regularities and anomalies of the solar magnetic activity, the negative buoyancy phenomenon and the magnetic quenching of the α effect were included in the model, as well as newest helioseismically determined inner rotation of the Sun were used. Magnetic buoyancy constrains the magnitude of toroidal field produced by the Ω effect near the bottom of the solar convection zone (SCZ). Therefore, we examined two “antibuoyancy” effects: i) macroscopic turbulent diamagnetism and ii) magnetic advection caused by vertical inhomogeneity of fluid density in the SCZ, which we call the ∇ρ effect. The Sun's rotation substantially modifies the ∇ρ effect. The reconstruction of the toroidal field was examined assuming the balance between mean‐field magnetic buoyancy, turbulent diamagnetism and the rotationally modified ∇ρ effect. It is shown that at high latitudes antibuoyancy effects block the magnetic fields in the deep layers of the SCZ, and so the most likely these deep‐rooted fields could not become apparent at the surface as sunspots. In the near‐equatorial region, however, the upward ∇ρ effect can facilitate magnetic fields of about 3000 – 4000 G to emerge through the surface at the sunspot belt. Allowance for the radial inhomogeneity of turbulent velocity in derivations of the helicity parameter resulted in a change of sign of the α effect from positive to negative in the northern hemisphere near the bottom of the SCZ. The change of sign is very important for direction of the Parker's dynamo‐waves propagation and for parity of excited magnetic fields. The period of the dynamo‐wave calculated with allowance for the magnetic quenching is about seven years, that agrees by order of magnitude with the observed mean duration of the sunspot cycles. Using the modern helioseismology data to define dynamo‐parameters, we conclude that north‐south asymmetry should exist in the meridional field. At low latitudes in deep layers of the SCZ, the αΩ dynamo excites most efficiency the dipolar mode of the meridional field. Meanwhile, in high‐latitude regions a quadrupolar mode dominates in the meridional field. The obtained configuration of the net meridional field is likely to explain the magnetic anomaly of polar fields (the apparent magnetic “monopole”) observed near the maxima of solar cycles. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On large-scale solar convection

We examine the effects of rotation about a vertical axis on thermal convection with a simple model in which an inviscid, incompressible fluid of zero thermal conductivity and electrical resistivity is contained in a thin annulus of rectangular cross-section. The initial steady state assumed is one of no motion relative to the rotating frame with constant (unstable) vertical temperature gradient and uniform toroidal magnetic field. Small periodic disturbances are then introduced and the linearized perturbation equations solved. We also determine the second-order mean circulations and magnetic fields that are forced by non-zero Reynolds and thermal stresses and magnetic field transports.The solutions have several properties which are relevant to large-scale solar phenomena if giant long-lived convection cells exist on the sun. In particular, the convective cells are tilted in latitude in the same sense as bipolar magnetic regions, and induce vertical magnetic fields with the same tilt. They transport momentum across latitude circles through Reynolds stresses and induced meridional circulations thus setting up a differential rotation. Cells which grow slowly compared to the rotation rate and have comparable dimensions in latitude and longitude transport momentum toward the equator. The cells also form a poloidal magnetic field from initial toroidal field, in a manner similar to that put forth by Parker.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

A model of the solar convection zone

A model of the convection zone is presented which matches an empirical model atmosphere (HSRA) and an interior model. A mixing length formalism containing four adjustable parameters is used. Thermodynamical considerations provide limits on two of these parameters. The average temperature-pressure relation depends on two or three combinations of the four parameters. Observational information on the structure of the outermost layers of the convection zone, and the value of the solar radius limit the range of possible parameter combinations. It is shown that in spite of the remaining freedom of choice of the parameters, the mean temperature-pressure relation is fixed well by these data.The reality of a small density inversion in the HSRA model is investigated. The discrepancy between the present model and a solar model by Mullan (1971) is discussed briefly.     

Testing turbulent convection theory in solar models – I. Structure of the solar convection zone

Turbulent convection models (TCMs) based on hydrodynamic moment equations are compared with the classical mixing-length theory (MLT) in solar models. The aim is to test the effects of some physical processes on the structure of the solar convection zone, such as the dissipation, diffusion and anisotropy of turbulence that have been ignored in the MLT. Free parameters introduced by the TCMs are also tested in order to find appropriate values for astrophysical applications. It is found that the TCMs usually give larger convective heat fluxes than the MLT does, and the heat transport efficiency is sensitively related to the dissipation parameters used in the TCMs. As a result of calibrating to the present solar values, our solar models usually have rather smaller values of the mixing length to local pressure scaleheight ratio than the standard solar model. The turbulent diffusion is found to have important effects on the structure of the solar convection zone. It leads to significantly lowered and expanded profiles for the Reynolds correlations, and a larger temperature gradient in the central part of the superadiabatic convection region but a smaller one near the boundaries of the convection zone. It is interesting to note that, due to a careful treatment of turbulence developing towards isotropic state, our non-local TCM results in radially dominated motion in the central part and horizontally dominated motion near the boundaries of the convection zone, just as what has been observed in many 3D numerical simulations. Our solar models with the TCMs give small but meaningful differences in the temperature and sound speed profiles compared with the standard solar model using the MLT.     

Large-scale convection cells in the solar atmosphere

The directions and velocities of meridional plasma motions are investigated using Doppler shifts of the magnetically non-split line Fe 557.6 nm. Possible drifts of the spectrograph were controlled by measuring nearly iodine lines from a laboratory source. The scattered light was kept low by using the Capri Coudé refractor mainly around local noon.There is a general scatter of about ±20 to ±50 m s^–1 in the yearly mean results. The scatter is up to –200 m s^–1 for the year 1985. The results are compared with published data. Although some systematic meridional plasma motions could be detected from the average of the meridional flows within the whole observing period 1982 until 1986, the final analysis suggests, that all meridional motions averaged over half a solar cycle are below ±10 m s^–1.Paper presented at the 11th European Regional Astronomical Meetings of the IAU on New Windows to the Universe, held 3–8 July, 1989, Tenerife, Canary Islands, Spain.This paper has also been presented at the IAU Colloquium No. 121Inside the Sun andAstronomy and Astrophysics 229, 224 (1990).     

Rotating convection and the solar differential rotation

We discuss the implication of a numerical experiment on rotating convection and its relevance to the construction of a model for the solar differential rotation.     

Discrete cellular scales of solar convection

The theoretical power spectrum of velocity fields and flux fluctuations at the solar photosphere is calculated using a quasi-nonlinear framework of superposition of unstable convective eigenmodes excited in the solar convection zone. It is demonstrated that this power spectrum exhibits at least three distinct peaks corresponding to granulation, mesogranulation and supergranulation. The vertical velocity and the brightness fluctuation at the solar surface are found to be correlated. The theoretical framework can be adopted for application to other types of stars in order to predict the dominant length scales in the power spectrum of convection in these stars.     

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.     

Possible solar cycle variations in the convection zone

Using data from the Global Oscillations Network Group (GONG) that covers the period from 1995 to 1998 we study the change in frequencies of solar oscillations with solar activity. From these frequencies we attempt to determine any possible variation in solar structure with solar activity. We do not find any evidence of a change in the convection zone depth or extent of overshoot below the convection zone during the solar cycle.     

Seismic study of solar convection and overshooting: results of nonlocal convection

Local mixing-length theory is incapable of describing nonlocal phenomena in stellar convection, such as overshooting. Therefore standard solar models constructed with local mixing-length theory significantly deviate from the Sun at the boundaries of the convection zone, where convection becomes less efficient and nonlocal effects are important. The differences between observed and computed frequencies mainly come from the region near the surface, while the localized difference in sound speed is just below the convective envelope. We compute a solar envelope model using Xiong’s nonlocal convection theory, and carry out helioseismic analysis. The nonlocal model has a smooth transition at the base of the convection zone, as revealed by helioseismology. It reproduces solar frequencies more accurately, and reduces the localized difference in sound speed between the Sun and standard solar models.     

Current convection in solar active regions

Current carrying magnetic fields which penetrate sunspots can be unstable to current convective modes caused by the large gradient of electrical conductivity. The linear growth rates and wavelengths of the unstable modes are found. The unstable modes produce fine-scale vortices perpendicular to the magnetic field, which overshoot well into the solar corona. The modes provide a turbulent vorticity source at the photospheric footpoints of the field. This can cause braiding and reconnection of the coronal magnetic field. The modes twist the coronal magnetic field into loops with a typical radius of 200 km, consistent with recent X-ray observations.     

Structure of the solar photospheric convection on subgranulation scales

We investigate the structure of convective flows in the solar photosphere on subgranulation scales. The solar granulation pattern is reproduced by solving the inverse problem of nonequilibrium radiation transfer on the basis of the profiles of the neutral iron line λ 523.42 nm. The wave motions are excluded by the k-ω filtration. The line-of-sight velocity has an asymmetric distribution inside the convective flows in large granules (1.5″ and larger) in the lower photosphere and at the bottom of the middle photosphere. This asymmetry is weaker in the upper photosphere. For smaller flows the distribution is more symmetric at all heights. The asymmetry of the temperature distribution is less pronounced. Large convective flows were found to have a fine structure: they are fragmentized into several smaller flows. The fine structure of large flows and spatial smearing are responsible for the observed asymmetry of the convection velocity distribution inside flows.     

Probing signatures of the alpha-effect in the solar convection zone

An attempt to extract maximum information on signatures of the alpha-effect from current helicity and twist density calculations in the solar photosphere is carried out. A possible interpretation of the results for developing the dynamo theory is discussed. The analysis shows that the surface magnetic current helicity is mainly negative/positive in the northern/southern hemispheres of the Sun. This indicates the actual alpha-effect at the photospheric level to be positive/negative, respectively. However, at the bottom of the convection zone, we may assume this effect to change the sign to negative/positive. We reveal some quantities related to the alpha-effect and discuss its spatial and temporal distribution. It is also found that there are a small number of active regions where the sign of the alpha-effect is opposite to that in most active regions. Such exceptional active regions seem to localize at certain active longitudes. We compare the determined regularities with theoretical predictions of the alpha-effect distribution in the solar convection zone.