The influence of the Coriolis force on flux tubes rising through the solar convection zone

In order to study the effect of the Coriolis force due to solar rotation on rising magnetic flux, the authors consider a flux ring, azimuthally symmetric around the rotation axis, starting from rest at the bottom of the convection zone, and then follow the trajectory of the flux ring as it rises. If it is assumed that the flux ring remains azimuthally symmetric during its ascent, then the problem can be described essentially in terms of two parameters: the value of the initial magnetic field in the ring when it starts, and the effective drag experienced by it. For field strengths at the bottom of the convection zone of order 10,000 G or less, it is found that the Coriolis force plays a dominant role and flux rings starting from low latitudes at the bottom are deflected and emerge at latitudes significantly poleward of sunspot zones.

The evolution of loop structures in flux rings within the solar convection zone

Choudhuri and Gilman (1987) considered certain implications of the hypothesis that the magnetic flux within the Sun is generated at the bottom of the convection zone and then rises through it. Taking flux rings symmetric around the rotation axis and using reasonable values of different parameters, they found that the Coriolis force deflects these flux rings into trajectories parallel to the rotation axis so that they emerge at rather high latitudes. This paper looks into the question of whether the action of the Coriolis force is subdued when the initial configuration of the flux ring has non-axisymmetries in the form of loop structures. The results depend dramatically on whether the flux ring with the loops lies completely within the convection zone or whether the lower parts of it are embedded in the stable layers underneath the convection zone. In the first case, the Coriolis force supresses the non-axisymmetric perturbations so that the flux ring tends to remain symmetric and the trajectories are very similar to those of Choudhuri and Gilman (1987). In the second case, however, the lower parts of the flux ring may remain anchored underneath the bottom of the convection zone, but the upper parts of the loops still tend to move parallel to the rotation axis and emerge at high latitudes. Thus the problem of the magnetic flux not being able to come out at the sunspot latitudes still persists after the non-axisymmetries in the flux rings are taken into account.

Effect of Turbulence on Emerging Magnetic Flux Tubes in the Convection Zone

Working under the hypothesis that magnetic flux in the sun is generated at the bottom of the convection zone, Choudhuri and Gilman (1987; Astrophys. J. 316, 788) found that a magnetic flux tube symmetric around the rotation axis, when released at the bottom of the convection zone, gets deflected by the Coriolis force and tends to move parallel to the rotation axis as it rises in the convection zone. As a result, all the flux emerges at rather high latitudes and the flux observed at the typical sunspot latitudes remains unexplained. Choudhuri (1989; Solar Physics, in press) finds that non-axisymmetric perturbations too cannot subdue the Coriolis force. In this paper, we no longer treat the convection zone to be passive as in the previous papers, but we consider the role of turbulence in the convection zone in inhibiting the Coriolis force. The interaction of the flux tubes with the turbulence is treated in a phenomenological way as follows: (1) Large scale turbulence on the scale of giant cells can physically drag the tubes outwards, thus pulling the flux towards lower latitudes by dominating over the Coriolis force. (2) Small scale turbulence of the size of the tubes can exchange angular momentum with the tube, thus suppressing the growth of the Coriolis force and making the tubes emerge at lower latitudes. Numerical simulations show that the giant cells can drag the tubes and make them emerge at lower latitufes only if the velocities within the giant cells are unrealistically large of if the radii of the flux tubes are as small as 10 km. However, small scale turbulence can successfully suppress the growth of the Coriolis force if the tubes have radii smaller than about 300 km which may not be unreasonable. Such flux tubes can then emerge at low latitudes where sunspots are seen.

An elementary introduction to solar dynamo theory

The cyclically varying magnetic field of the Sun is believed to be produced by the hydromagnetic dynamo process. We first summarize the relevant observational data pertaining to sunspots and solar cycle. Then we review the basic principles of MHD needed to develop the dynamo theory. This is followed by a discussion how bipolar sunspots form due to magnetic buoyancy of flux tubes formed at the base of the solar convection zone. Following this, we come to the heart of dynamo theory. After summarizing the basic ideas of a turbulent dynamo and the basic principles of its mean field formulation, we present the famous dynamo wave solution, which was supposed to provide a model for the solar cycle. Finally we point out how a flux transport dynamo can circumvent some of the difficulties associated with the older dynamo models.

Reynolds stress and heat flux in spherical shell convection

Context. Turbulent fluxes of angular momentum and heat due to rotationally affected convection play a key role in determining differential rotation of stars. Aims. We compute turbulent angular momentum and heat transport as functions of the rotation rate from stratified convection. We compare results from spherical and Cartesian models in the same parameter regime in order to study whether restricted geometry introduces artefacts into the results. Methods. We employ direct numerical simulations of turbulent convection in spherical and Cartesian geometries. In order to alleviate the computational cost in the spherical runs and to reach as high spatial resolution as possible, we model only parts of the latitude and longitude. The rotational influence, measured by the Coriolis number or inverse Rossby number, is varied from zero to roughly seven, which is the regime that is likely to be realised in the solar convection zone. Cartesian simulations are performed in overlapping parameter regimes. Results. For slow rotation we find that the radial and latitudinal turbulent angular momentum fluxes are directed inward and equatorward, respectively. In the rapid rotation regime the radial flux changes sign in accordance with earlier numerical results, but in contradiction with theory. The latitudinal flux remains mostly equatorward and develops a maximum close to the equator. In Cartesian simulations this peak can be explained by the strong 'banana cells'. Their effect in the spherical case does not appear to be as large. The latitudinal heat flux is mostly equatorward for slow rotation but changes sign for rapid rotation. Longitudinal heat flux is always in the retrograde direction. The rotation profiles vary from anti-solar (slow equator) for slow and intermediate rotation to solar-like (fast equator) for rapid rotation. The solar-like profiles are dominated by the Taylor-Proudman balance.

The origin of the solar magnetic cycle

After summarizing the relevant observational data, we discuss how a study of flux tube dynamics in the solar convection zone helps us to understand the formation of sunspots. Then we introduce the flux transport dynamo model and assess its success in modelling both the solar cycle and its departures from strictly periodic behaviour.

POLAR NETWORK INDEX AS A MAGNETIC PROXY FOR THE SOLAR CYCLE STUDIES

The Sun has a polar magnetic field which oscillates with the 11 yr sunspot cycle. This polar magnetic field is an important component of the dynamo process which operates in the solar convection zone and produces the sunspot cycle. We have direct systematic measurements of the Sun's polar magnetic field only from about the mid-1970s. There are, however, indirect proxies which give us information about this field at earlier times. The Ca-K spectroheliograms taken at the Kodaikanal Solar Observatory during 1904-2007 have now been digitized with 4k x 4k CCD and have higher resolution (similar to 0.86 arcsec) than the other available historical data sets. From these Ca-K spectroheliograms, we have developed a completely new proxy (polar network index, hereafter PNI) for the Sun's polar magnetic field. We calculate PNI from the digitized images using an automated algorithm and calibrate our measured PNI against the polar field as measured by the Wilcox Solar Observatory for the period 1976-1990. This calibration allows us to estimate the polar fields for the earlier period up to 1904. The dynamo calculations performed with this proxy as input data reproduce reasonably well the Sun's magnetic behavior for the past century.

Mechanisms for mhd poynting flux generation in simulations of solar photospheric magnetoconvection

We investigate the generation mechanisms of MHD Poynting flux in the magnetized solar photosphere. Using radiative MHD modeling of the solar photosphere with initial magnetic configurations that differ in their field strength and geometry, we show the presence of two different mechanisms for MHD Poynting flux generation in simulations of solar photospheric magnetoconvection. The weaker mechanism is connected to vertical transport of weak horizontal magnetic fields in the convectively stable layers of the upper photosphere, while the stronger is the production of Poynting flux in strongly magnetized intergranular lanes experiencing horizontal vortex motions. These mechanisms may be responsible for the energy transport from the solar convection zone to the higher layers of the solar atmosphere.

Reconstruction des structures magnéto-convectives solaires sous une région active, par l’utilisation conjointe d’un modèle de convection anélastique et d’une méthode d’assimilation de données.

L’utilisation d’une méthode d’assimilation de données, associée à un modèle de convection anélastique, nous permet la reconstruction des structures physiques d’une partie de la zone convective située en dessous d’une région solaire active. Les résultats obtenus nous informent sur les processus d’émergence des tubes de champ magnétique au travers de la zone convective ainsi que sur les mécanismes de formation des régions actives. Les données solaires utilisées proviennent de l’instrument MDI à bord de l’observatoire spatial SOHO et concernent principalement la région active AR9077 lors de l’ ́évènement du “jour de la Bastille”, le 14 juillet 2000. Cet évènement a conduit à l’avènement d’une éruption solaire, suivie par une importante éjection de masse coronale. Les données assimilées (magnétogrammes, cartes de températures et de vitesses verticales) couvrent une surface de 175 méga-mètres de coté acquises au niveau photosphérique. La méthode d’assimilation de données employée est le “coup de coude direct et rétrograde”, une méthode de relaxation Newtonienne similaire à la méthode “quasi-linéaire inverse 3D”. Elle présente l’originalité de ne pas nécessiter le calcul des équations adjointes au modèle physique. Aussi, la simplicité de la méthode est un avantage numérique conséquent. Notre étude montre au travers d’un test simple l’applicabilité de cette méthode à un modèle de convection utilisé dans le cadre de l’approximation anélastique. Nous montrons ainsi l’efficacité de cette méthode et révélons son potentiel pour l’assimilation de données solaires. Afin d’assurer l’unicité mathématique de la solution obtenue nous imposons une régularisation dans tout le domaine simulé. Nous montrons enfin que l’intérêt de la méthode employée ne se limite pas à la reconstruction des structures convectives, mais qu’elle permet également l’interpolation optimale des magnétogrammes photosphériques, voir même la prédiction de leur évolution temporelle.

The effect of closed boundary conditions on a stationary dynamo

One of two boundary conditions generally assumed in solutions of the dynamo equation is related to the disappearance of the azimuthal field at the boundary. Parker (1984) points out that for the realization of this condition the field must escape freely through the surface. Escape requires that the field be detached from the gas in which it is embedded. In the case of the sun, this can be accomplished only through reconnection in the tenuous gas above the visible surface. Parker concludes that the observed magnetic activity on the solar surface permits at most three percent of the emerging flux to escape. He arrives at the conclusion that, instead of B(phi) = 0, the partial derivative of B(phi) to r is equal to zero. The present investigation is concerned with the effect of changing the boundary condition according to Parker's conclusion. Implications for the solar convection zone are discussed.

Les oscillations torsionnelles dans la zone de convection solaire

Nous analysons les oscillations torsionnelles se développant dans une simulation magnétohydrodynamique de la zone de convection solaire produisant des champs magnétiques de type solaire (champs axisymétriques subissant des inversions de polarités régulières sur des échelles temporelles décadaires). Puisque ces oscillations sont également similaires à celles observées dans le Soleil, nous analysons les dynamiques zonales aux grandes échelles. Nous séparons donc les termes aux grandes échelles (force de Coriolis exercée sur la circulation méridienne et les champs magnétiques aux grandes échelles) de ceux aux petites échelles (les stress de Reynolds et de Maxwell). En comparant les flux de moments cinétiques entre chacune des composantes, nous nous apercevons que les oscillations torsionnelles sont maintenues par l’écoulement méridien aux grandes échelles, lui même modulé par les champs magnétiques. Une analyse d’échange d’énergie confirme ce résultat, puisqu’elle montre que seul le terme comprenant la force de Coriolis injecte de l’énergie dans l’écoulement. Une analyse de la dynamique rotationnelle ayant lieu à la limite de la zone stable et de la zone de convection démontre que celle-ci est fortement modifiée lors du passage de la base des couches convectives à la base de la fine tachocline s’y formant juste en-dessous. Nous concluons par une discussion au niveau du mécanisme de saturation en amplitude dans la dynamo s’opérant dans la simulation ainsi que de la possibilité d’utiliser les oscillations torsionnelles comme précurseurs aux cycles solaires à venir.

Locating the seat of the solar dynamo

The hypothesis that the solar dynamo operates in a thin layer at the bottom of the convection zone is addressed. Recent work on the question whether the magnetic flux can be made to emerge at sunspot latitudes is reviewed. It is concluded that this hypothesis can fit the observational facts only if there is turbulence with a length scale of a few hundred kilometers in and around the dynamo region.

A theoretical model for the magnetic helicity of solar active regions

Active regions on the solar surface are known to possess magnetic helicity, which is predominantly negative in the northern hemisphere and positive in the southern hemisphere. Choudhuri et al. [Choudhuri, A.R. On the connection between mean field dynamo theory and flux tubes. Solar Phys. 215, 31–55, 2003] proposed that the magnetic helicity arises due to the wrapping up of the poloidal field of the convection zone around rising flux tubes which form active regions. Choudhuri [Choudhuri, A.R., Chatterjee, P., Nandy, D. Helicity of solar active regions from a dynamo model. ApJ 615, L57–L60, 2004] used this idea to calculate magnetic helicity from their solar dynamo model. Apart from getting broad agreements with observational data, they also predict that the hemispheric helicity rule may be violated at the beginning of a solar cycle. Chatterjee et al. [Chatterjee, P., Choudhuri, A.R., Petrovay, K. Development of twist in an emerging magnetic flux tube by poloidal field accretion. A&A 449, 781–789, 2006] study the penetration of the wrapped poloidal field into the rising flux tube due to turbulent diffusion using a simple 1-d model. They find that the extent of penetration of the wrapped field will depend on how weak the magnetic field inside the rising flux tube becomes before its emergence. They conclude that more detailed observational data will throw light on the physical conditions of flux tubes just before their emergence to the photosphere.

Locating the seat of the solar dynamo

The hypothesis that the solar dynamo operates in a thin layer at the bottom of the convection zone is addressed. Recent work on the question whether the magnetic flux can be made to emerge at sunspot latitudes is reviewed. It is concluded that this hypothesis can fit the observational facts only if there is turbulence with a length scale of a few hundred kilometers in and around the dynamo region.