Helicity of solar active regions from a dynamo model

We calculate helicities of solar active regions based on the idea that poloidal flux lines get wrapped around a toroidal flux tube rising through the convection zone, thereby giving rise to the helicity. We use our solar dynamo model based on the Babcock-Leighton α-effect to study how helicity varies with latitude and time.     

On global solar dynamo simulations

Global dynamo simulations solving the equations of magnetohydrodynamics (MHD) have been a tool of astrophysicists who try to understand the magnetism of the Sun for several decades now. During recent years many fundamental issues in dynamo theory have been studied in detail by means of local numerical simulations that simplify the problem and allow the study of physical effects in isolation. Global simulations, however, continue to suffer from the age‐old problem of too low spatial resolution, leading to much lower Reynolds numbers and scale separation than in the Sun. Reproducing the internal rotation of the Sun, which plays a crucial role in the dynamo process, has also turned out to be a very difficult problem. In the present paper the current status of global dynamo simulations of the Sun is reviewed. Emphasis is put on efforts to understand how the large‐scale magnetic fields, i.e. whose length scale is greater than the scale of turbulence, are generated in the Sun. Some lessons from mean‐field theory and local simulations are reviewed and their possible implications to the global models are discussed. Possible remedies to some current issues of solar simulations are put forward (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Towards a model for the solar dynamo

We study a mean field model of the solar dynamo, in which the non-linearity is the action of the azimuthal component of the Lorentz force of the dynamo-generated magnetic field on the angular velocity. The underlying zero-order angular velocity is consistent with recent determinations of the solar rotation law, and the form of the alpha effect is chosen so as to give a plausible butterfly diagram. For small Prandtl numbers we find regular, intermittent and apparently chaotic behaviour, depending on the size of the alpha coefficient. For certain parameters, the intermittency displays some of the characteristics believed to be associated with the Maunder minimum. We thus believe that we are capturing some features of the solar dynamo.     

Numerical simulations of three-dimensional magnetic swirls in a solar flux-tube

We aim to numerically study evolution of Alfv′en waves that accompany short-lasting swirl events in a solar magnetic flux-tube that can be a simple model of a magnetic pore or a sunspot. With the use of the FLASH code we numerically solve three-dimensional ideal magnetohydrodynamic equations to simulate twists which are implemented at the top of the photosphere in magnetic field lines of the flux-tube. Our numerical results exhibit swirl events and Alfv′en waves with associated clockwise and counterclockwise rotation of magnetic lines, with the largest values of vorticity at the bottom of the chromosphere, and a certain amount of energy flux.     

Alleviation of catastrophic quenching in solar dynamo model with nonlocal alpha‐effect

The nonlocal alpha‐effect of Babcock‐Leighton type is not prone to the catastrophic quenching due to conservation of magnetic helicity. This is shown with a dynamo model, which jointly applies the nonlocal alpha‐effect, the diamagnetic pumping, and dynamical equation for the magnetic alpha‐effect. The same model shows catastrophic quenching when the alpha‐effect is changed to its local formulation. The nonlocal model shows a preferred excitation of magnetic fields of dipolar symmetry, which oscillate with a period of about ten years and have a toroidal‐to‐polar fields ratio of about a thousand (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The solar dynamo in the light of the distribution of various sunspot magnetic classes over butterfly diagram

We present the data concerning the distribution of various sunspot magnetic classes over the solar butterfly diagram and discuss how this data can inform solar dynamo models. We use the statistics of sunspots that violate the Hale polarity law to estimate the ratio of the fluctuating and mean components of the toroidal magnetic field inside the solar convective zone. An analysis of the spatial distribution of bipolar, unipolar and complex sunspot groups in the context of simple dynamo models results in the conclusion that the mean toroidal field is relatively simple and maintains its shape during the course of the solar cycle (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Search for relationship between duration of the extended solar cycles and amplitude of sunspot cycle

Duration of the extended solar cycles is taken into the consideration. The beginning of cycles is counted from the moment of polarity reversal of large-scale magnetic field in high latitudes, occurring in the sunspot cycle n till the minimum of the cycle n + 2. The connection between cycle duration and its amplitude is established. Duration of the “latent” period of evolution of extended cycle between reversals and a minimum of the current sunspot cycle is entered. It is shown, that the latent period of cycles evolution is connected with the next sunspot cycle amplitude and can be used for the prognosis of a level and time of a sunspot maximum. The 24th activity cycle prognosis is made. The found dependences correspond to transport dynamo model of generation of solar cyclicity, it is possible with various speed of meridional circulation. Long-term behavior of extended cycle's lengths and connection with change of a climate of the Earth is considered. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A dynamo model for axisymmetric and non-axisymmetric solar magnetic fields

More and more observations are showing a relatively weak, but persistent, non-axisymmetric magnetic field co-existing with the dominant axisymmetric field on the Sun. Its existence indicates that the non-axisymmetric magnetic field plays an important role in the origin of solar activity. A linear non-axisymmetric  α²– Ω  dynamo model is derived to explore the characteristics of the axisymmetric  ( m = 0)  and the first non-axisymmetric  ( m = 1)  modes and to provide a theoretical basis with which to explain the 'active longitude', 'flip-flop' and other non-axisymmetric phenomena. The model consists of an updated solar internal differential rotation, a turbulent diffusivity varying with depth, and an α-effect working at the tachocline in a rotating spherical system. The difference between the  α²–Ω  and the  α–Ω  models and the conditions that favour the non-axisymmetric modes under solar-like parameters are also presented.     

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.     

Some recent developments in solar dynamo theory

We discuss the current status of solar dynamo theory and describe the dynamo model developed by our group. The toroidal magnetic field is generated in the tachocline by the strong differential rotation and rises to the solar surface due to magnetic buoyancy to create active regions. The decay of these active regions at the surface gives rise to the poloidal magnetic field by the Babcock-Leighton mechanism. This poloidal field is advected by the meridional circulation first to high latitudes and then down below to the tachocline. Dynamo models based on these ideas match different aspects of observational data reasonably well.     

What is a flux tube? On the magnetic field topology of buoyant flux structures

We study the topology of field lines threading buoyant magnetic flux structures. The magnetic structures, visually resembling idealized magnetic flux tubes, are generated self-consistently by numerical simulation of the interaction of magnetic buoyancy and a localized velocity shear in a stably stratified atmosphere. Depending on the parameters, the system exhibits varying degrees of symmetry. By integrating along magnetic field lines and constructing return maps, we show that, depending on the type of underlying behaviour, the stages of the evolution, and therefore the degree of symmetry, the resulting magnetic structures can have field lines with one of three distinct topologies. When the x -translational and y -reflectional symmetries remain intact, magnetic field lines lie on surfaces but individual lines do not cover the surface. When the y symmetry is broken, magnetic field lines lie on surfaces and individual lines do cover the surface. When both x and y symmetries are broken, magnetic field lines wander chaotically over a large volume of the magnetically active region. We discuss how these results impact our simple ideas of a magnetic flux tube as an object with an inside and an outside, and introduce the concept of 'leaky' tubes.     

The role of magnetic buoyancy in a Babcock-Leighton type solar dynamo

We study the effects of incorporating magnetic buoyancy in a model of the solar dynamo—which draws inspiration from the Babcock-Leighton idea of surface processes generating the poloidal field. We present our main results here.     

Equatorial magnetic helicity flux in simulations with different gauges

We use direct numerical simulations of forced MHD turbulence with a forcing function that produces two different signs of kinetic helicity in the upper and lower parts of the domain. We show that the mean flux of magnetic helicity from the small‐scale field between the two parts of the domain can be described by a Fickian diffusion law with a diffusion coefficient that is approximately independent of the magnetic Reynolds number and about one third of the estimated turbulent magnetic diffusivity. The data suggest that the turbulent diffusive magnetic helicity flux can only be expected to alleviate catastrophic quenching at Reynolds numbers of more than several thousands. We further calculate the magnetic helicity density and its flux in the domain for three different gauges. We consider the Weyl gauge, in which the electrostatic potential vanishes, the pseudo‐Lorenz gauge, where the speed of light is replaced by the sound speed, and the ‘resistive gauge’ in which the Laplacian of the magnetic vector potential acts as a resistive term. We find that, in the statistically steady state, the time‐averaged magnetic helicity density and the magnetic helicity flux are the same in all three gauges (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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)