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.     

Flux-dominated solar dynamo model with a thin shear layer

Flux-dominated solar dynamo models have demonstrated to reproduce the main features of the large scale solar magnetic cycle, however the use of a solar like differential rotation profile implies in the the formation of strong toroidal magnetic fields at high latitudes where they are not observed. In this work, we invoke the hypothesis of a thin-width tachocline in order to confine the high-latitude toroidal magnetic fields to a small area below the overshoot layer, thus avoiding its influence on a Babcock-Leighton type dynamo process. Our results favor a dynamo operating inside the convection zone with a tachocline that essentially works as a storage region when it coincides with the overshoot layer. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Differential rotation and meridional flow in the solar convection zone and beneath

The influence of the basic rotation on anisotropic and inhomogeneous turbulence is discussed in the context of differential rotation theory. An improved representation for the original turbulence leads to a Λ‐effect which complies with the results of 3D numerical simulations. The resulting rotation law and meridional flow agree well with both the surface observations (∂Ω/∂r < 0 and meridional flow towards the poles) and with the findings of helioseismology. The computed equatorward flow at the bottom of convection zone has an amplitude of about 10 m/s and may be significant for the solar dynamo. The depth of the meridional flow penetration into the radiative zone is proportional to ν^0.5_core, where ν_core is the viscosity beneath the convection zone. The penetration is very small if the tachocline is laminar. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Self-excitation of motions in near-convection zones and the generation of solar magnetic field

In view of the recently discovered time variations in rotation velocity within the solar differentially rotating tachocline (Howe et al. 2000), we study conditions for the equilibrium and excitation of motions in nonrigidly rotating magnetized layers of the radiative zones located near the boundaries of the convection zone. The emphasis is on the possible relationship between quasi-periodic tachocline pulsations and the generation of a nonaxisymmetric magnetic field in the convection zone. This field generation is studied under the assumption that it results from a reduction in the expenditure of energy on convective heat transport. The (antisymmetric about the equator) field is shown to increase in strength if there are both a radial gradient in angular velocity and steady-state axisymmetric meridional circulation of matter. The sense of circulation is assumed to change (causing the sign of the generated field to change) after the maximum permissible field strength is reached. This is apparently attributable to the excitation of the corresponding turbulent viscosity of the medium. It is also important that the cyclic field variations under discussion are accompanied by variations in solar-type dipole magnetic field.     

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.     

Dynamo action in simulations of penetrative solar convection with an imposed tachocline

We summarize new and continuing three-dimensional spherical shell simulations of dynamo action by convection allowed to penetrate downward into a tachocline of rotational shear. The inclusion of an imposed tachocline allows us to examine several processes believed to be essential in the operation of the global solar dynamo, including differential rotation, magnetic pumping, and the stretching and organization of fields within the tachocline. In the stably stratified core, our simulations reveal that strong axisymmetric magnetic fields (of ∼ 3000 G strength) can be built, and that those fields generally exhibit a striking antisymmetric parity, with fields in the northern hemisphere largely of opposite polarity to those in the southern hemisphere. In the convection zone above, fluctuating fields dominate over weaker mean fields. New calculations indicate that the tendency toward toroidal fields of antisymmetric parity is relatively insensitive to initial magnetic field configurations; they also reveal that on decade-long timescales, the magnetic fields can briefly enter (and subsequently emerge from) states of symmetric parity.We have not yet observed any overall reversals of the field polarity, nor systematic latitudinal propagation. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dynamics of the solar tachocline – I. An incompressible study

Gough & McIntyre have suggested that the dynamics of the solar tachocline are dominated by the advection–diffusion balance between the differential rotation, a large-scale primordial field and baroclinicly driven meridional motions. This paper presents the first part of a study of the tachocline, in which a model of the rotation profile below the convection zone is constructed along the lines suggested by Gough & McIntyre and solved numerically. In this first part, a reduced model of the tachocline is derived in which the effects of compressibility and energy transport on the system are neglected; the meridional motions are driven instead by Ekman–Hartmann pumping. Through this simplification, the interaction of the fluid flow and the magnetic field can be isolated and is studied through non-linear numerical analysis for various field strengths and diffusivities. It is shown that there exists only a narrow range of magnetic field strengths for which the system can achieve a nearly uniform rotation. The results are discussed with respect to observations and to the limitations of this initial approach. A following paper combines the effects of realistic baroclinic driving and stratification with a model that closely follows the lines of work of Gough & McIntyre.     

On the possibility of a bimodal solar dynamo

A simple way to couple an interface dynamo model to a fast tachocline model is presented, under the assumption that the dynamo saturation is due to a quadratic process and that the effect of finite shear layer thickness on the dynamo wave frequency is analogous to the effect of finite water depth on surface gravity waves. The model contains one free parameter which is fixed by the requirement that a solution should reproduce the helioseismically determined thickness of the tachocline. In this case it is found that, in addition to this solution, another steady solution exists, characterized by a four times thicker tachocline and 4–5 times weaker magnetic fields. It is tempting to relate the existence of this second solution to the occurrence of grand minima in solar activity. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Latitudinal shear instability in the solar tachocline

In an attempt to explain the observed rotation profile in the solar radiative zone and the tachocline, Spiegel & Zahn proposed a model based on anisotropic turbulent angular momentum transport. Although very successful in reproducing some of the features of the solar tachocline, their model assumes without verification that the origin of the turbulence could be caused by latitudinal shear instability. This paper studies the weakly non-linear evolution of two-dimensional shear instability, in which the interaction between the global rotation profile and the Reynolds stresses can be described self-consistently. Provided that the initial rotation profile is sufficiently close to marginal stability (which is the case of the solar tachocline), the instability is shown to saturate and to relax to a marginally stable state, which differs very little from the observed rotation profile. It is therefore likely that the tachocline is in a state of marginal stability with respect to latitudinal shear instability, and shows that angular momentum transport in the tachocline is unlikely to be caused by shear-induced turbulence.     

Looking inside the butterfly diagram

The suitability of Maunder's butterfly diagram to give a realistic picture of the photospheric magnetic flux large scale distribution is discussed. The evolution of the sunspot zone in cycle 20 through 23 is described. To reduce the noise which covers any structure in the diagram, a smoothing algorithm has been applied to the sunspot data. This operation has eliminated any short period fluctuation, and given visibility to long duration phenomena. One of these phenomena is the fact that the equatorward drift of the spot zone center of mass results from the alternation of several prograde (namely, equatorward) segments with other stationary or poleward segments. The long duration of the stationary/retrograde phases as well as the similarities among the spot zone alternating paths in the cycles under examination prevent us from considering these features as meaningless fluctuations, randomly superimposed on the continuous equatorward migration. On the contrary, these features should be considered physically meaningful phenomena, requiring adequate explanations. Moreover, even the smoothed spotted area markedly oscillates. The compared examination of area and spot zone evolution allows us to infer details about the spotted area distribution inside the butterfly diagram. Links between the changing structure of the spot zone and the tachocline rotation rate oscillations are proposed. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Meridional flow profile measurements with SOHO/MDI

We present meridional flow measurements of the Sun using a novel helioseismic approach for analyzing SOHO/MDI data in order to push the current limits in radial depth. Analyzing three consecutive months of data during solar minimum, we find that the meridional flow is as expected poleward in the upper convection zone, turns equatorward at a depth of around 40 Mm (∼ 0.95 R_⊙), and possibly changes direction again in the lower convection zone. This may indicate two meridional circulation cells in each hemisphere, one beneath the other. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

Time Variability of Rotation in Solar Convection Zone From soi-mdi

The variation of rotation in the convection zone over a period of two years from mid-1996 is studied using inversions of SOI–MDI data. We confirm the existence of near-surface banded zonal flows migrating towards the equator from higher latitudes, and reveal that these banded flows extend substantially beneath the surface, possibly to depths as great as 70 Mm (10% of the solar radius). Our results also reveal apparently significant temporal variations in the rotation rate at high latitudes and in the vicinity of the tachocline over the period of study.     

A two-layer αω-dynamo model with dynamic feedback on the ω-effect

In an attempt to produce a simple representation of an interface dynamo, I examine a dynamo model composed of two one-dimensional (radially averaged) pseudo-spherical layers, one in the convection zone and possessing an α-effect, and the other in the tachocline and possessing an ω-effect. The two layers communicate by means of an analogue of Newton's law of cooling, and a dynamical back-reaction of the magnetic field on ω is provided. Extensive bifurcation diagrams are calculated for three separate values of η, the ratio of magnetic diffusivities of the two layers. I find recognizable similarities to, but also dramatic differences from, the comparable one-layer model examined by Roald &38; Thomas. In particular, the solar-like dynamo mode found previously is no longer stable in the two-layer version; in its place there is a sequence of periodic, quasi-periodic and chaotic modes probably created in a homoclinic bifurcation. These differences are important enough to provide support for the view that the solar dynamo cannot be meaningfully modelled in one dimension.     

Does the tachocline show solar cycle related changes?

The tachocline at the base of the convection zone is generally believed to be the seat of the solar dynamo. Here we investigate whether the tachocline shows any detectable change using several 72 day time-series of the Michelson Doppler Imager (MDI) Medium-l data. We do not find any clear evidence of change with time.     

Diamagnetic pumping near the base of a stellar convection zone

The property of inhomogeneous turbulence in conducting fluids to expel large‐scale magnetic fields in the direction of decreasing turbulence intensity is shown as important for the magnetic field dynamics near the base of a stellar convection zone. The downward diamagnetic pumping confines a fossil internal magnetic field in the radiative core so that the field geometry is appropriate for formation of the solar tachocline. For the stars of solar age, the diamagnetic confinement is efficient only if the ratio of turbulent magnetic diffusivity η_T of the convection zone to the (microscopic or turbulent) diffusivity η_in of the radiative interior is η_T/η_in ≥ 10⁵. Confinement in younger stars requires larger η_T/η_in. The observation of persistent magnetic structures on young solar‐type stars can thus provide evidence for the nonexistence of tachoclines in stellar interiors and on the level of turbulence in radiative cores. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The structure of the solar convective overshooting zone

Using a non-local theory of convection, we calculated the structure of the solar convection zone, paying special attention to the detailed structure of the lower overshooting zone. Our results show that an extended transition zone exists near the bottom of the convection zone, where the temperature gradient turns smoothly from adiabatic in the convection zone to radiative in solar interior. A super-radiative temperature region is found in the overshooting zone under the solar convection zone, where     ,     ,     and     . The extension of the super-radiative region (defined by     l is about 0.63  H _P (0.053 R_⊙). A careful comparison of the distribution of adiabatic sound speed and density with the local one is carried out. It is found, strikingly, that the distribution of adiabatic sound speed and density of our model is roughly consistent with the results of reversion from solar oscillation observations.     

Solar activity forecast with a dynamo model

Although systematic measurements of the Sun's polar magnetic field exist only from mid-1970s, other proxies can be used to infer the polar field at earlier times. The observational data indicate a strong correlation between the polar field at a sunspot minimum and the strength of the next cycle, although the strength of the cycle is not correlated well with the polar field produced at its end. This suggests that the Babcock–Leighton mechanism of poloidal field generation from decaying sunspots involves randomness, whereas the other aspects of the dynamo process must be reasonably ordered and deterministic. Only if the magnetic diffusivity within the convection zone is assumed to be high (of order  10¹² cm² s⁻¹  ), we can explain the correlation between the polar field at a minimum and the next cycle. We give several independent arguments that the diffusivity must be of this order. In a dynamo model with diffusivity like this, the poloidal field generated at the mid-latitudes is advected toward the poles by the meridional circulation and simultaneously diffuses towards the tachocline, where the toroidal field for the next cycle is produced. To model actual solar cycles with a dynamo model having such high diffusivity, we have to feed the observational data of the poloidal field at the minimum into the theoretical model. We develop a method of doing this in a systematic way. Our model predicts that cycle 24 will be a very weak cycle. Hemispheric asymmetry of solar activity is also calculated with our model and compared with observational data.