Is meridional circulation important in modelling irregularities of the solar cycle?

We explore the importance of meridional circulation variations in modelling the irregularities of the solar cycle by using the flux transport dynamo model. We show that a fluctuating meridional circulation can reproduce some features of the solar cycle like the Waldmeier effect and the grand minimum. However, we get all these results only if the value of the turbulent diffusivity in the convection zone is reasonably high.

TURBULENT PUMPING OF MAGNETIC FLUX REDUCES SOLAR CYCLE MEMORY AND THUS IMPACTS PREDICTABILITY OF THE SUN'S ACTIVITY

Prediction of the Sun's magnetic activity is important because of its effect on space environment and climate. However, recent efforts to predict the amplitude of the solar cycle have resulted in diverging forecasts with no consensus. Yeates et al. have shown that the dynamical memory of the solar dynamo mechanism governs predictability, and this memory is different for advection- and diffusion-dominated solar convection zones. By utilizing stochastically forced, kinematic dynamo simulations, we demonstrate that the inclusion of downward turbulent pumping of magnetic flux reduces the memory of both advection- and diffusion-dominated solar dynamos to only one cycle; stronger pumping degrades this memory further. Thus, our results reconcile the diverging dynamo-model-based forecasts for the amplitude of solar cycle 24. We conclude that reliable predictions for the maximum of solar activity can be made only at the preceding minimum-allowing about five years of advance planning for space weather. For more accurate predictions, sequential data assimilation would be necessary in forecasting models to account for the Sun's short memory.

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.

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.

The hemispheric asymmetry of solar activity during the last century and the solar dynamo

We believe the Babcock-Leighton process of poloidal field generation to be the main source of irregularity in the solar cycle. The random nature of this process may make the poloidal field in one hemisphere stronger than that in the other hemisphere at the end of a cycle. We expect this to induce an asymmetry in the next sunspot cycle. We look for evidence of this in the observational data and then model it theoretically with our dynamo code. Since actual polar field measurements exist only from the 1970s, we use the polar faculae number data recorded by Sheeley (1991, 2008) as a proxy of the polar field and estimate the hemispheric asymmetry of the polar field in different solar minima during the major part of the twentieth century. This asymmetry is found to have a reasonable correlation with the asymmetry of the next cycle. We then run our dynamo code by feeding information about this asymmetry at the successive minima and compare the results with observational data. We find that the theoretically computed asymmetries of different cycles compare favorably with the observational data, with the correlation coefficient being 0.73. Due to the coupling between the two hemispheres, any hemispheric asymmetry tends to get attenuated with time. The hemispheric asymmetry of a cycle either from observational data or from theoretical calculations statistically tends to be less than the asymmetry in the polar field (as inferred from the faculae data) in the preceding minimum. This reduction factor turns out to be 0.43 and 0.51 respectively in observational data and theoretical simulations.

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.

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.

Stochastic fluctuations of the solar dynamo

Attempts in the past to model the irregularities of the solar cycle (such as the Maunder minimum) were based on studies of the nonlinear feedback of magnetic fields on the dynamo source terms. Since the alpha-coefficient is obtained by averaging over the turbulence, it is expected to have stochastic fluctuations, and we show that these fluctuations can explain the irregularities of the solar cycle in a more satisfactory way. We solve the dynamo equations in a slab with a single mode, taking the alpha-coefficient to be constant in space but fluctuating stochastically in time with some given amplitude and given correlation time. The same level of percentile fluctuations (about 10 %) produces no effect on an alpha-omega dynamo, but makes an alpha-2 dynamo completely chaotic. The level of irregularities in an alpha-2-omega dynamo qualitatively agrees with the solar behavior, reinforcing the conclusion of Choudhuri (1990a) that the solar dynamo is of the alpha-2-omega-type. The irregularities are found to increase on increasing either the amplitude or the correlation time of the stochastic fluctuations. The alpha-quenching mechanism tends to make the system stable against the irregularities and hence it is inferred that the alpha-quenching should not be too strong so that the irregularities are not completely suppressed. We also present a simple-minded analysis to understand why the stochastic fluctuations in the alpha-omega, alpha-2-omega and alpha-2 regimes have such different outcomes.

The solar dynamo with meridional circulation

We show that meridional circulation can have a profound influence on dynamo models for the solar cycle. Motivated by the observed tilt angles of sunspot groups we assume that the generation of the poloidal field takes place near the surface, while a shear layer of radial differential rotation produces the toroidal field at the bottom of the convection zone. Both layers are coupled by a circulation with a poleward directed flow in the upper part and an equatorward flow in the deep layers of the convection zone. The circulation forces the toroidal field belts (which are responsible for the surface activity) to move equatorward. This leads to butterfly diagrams in qualitative agreement with the observations, even if the dynamo wave would propagate poleward in the absence of circulation. This result opens the possibility to construct models for the solar cycle which are based on observational data (tilt angles, differential rotation, and meridional circulation).

Correlation Between Decay Rate and Amplitude of Solar Cycles as Revealed from Observations and Dynamo Theory

Using different proxies of solar activity, we have studied the following features of the solar cycle: i) The linear correlation between the amplitude of cycle and its decay rate, ii) the linear correlation between the amplitude of cycle and the decay rate of cycle , and iii) the anti-correlation between the amplitude of cycle and the period of cycle . Features ii) and iii) are very useful because they provide precursors for future cycles. We have reproduced these features using a flux-transport dynamo model with stochastic fluctuations in the Babcock-Leighton effect and in the meridional circulation. Only when we introduce fluctuations in meridional circulation, are we able to reproduce different observed features of the solar cycle. We discuss the possible reasons for these correlations.

Importance of meridional circulation i flux transport dynamo: the possibility of a maunder-like grand minimum

Meridional circulation is an important ingredient in flux transport dynamo models. We have studied its importance on the period, the amplitude of the solar cycle, and also in producing Maunder-like grand minima in these models. First, we model the periods of the last 23 sunspot cycles by varying the meridional circulation speed. If the dynamo is in a diffusion-dominated regime, then we find that most of the cycle amplitudes also get modeled up to some extent when we model the periods. Next, we propose that at the beginning of the Maunder minimum the amplitude of meridional circulation dropped to a low value and then after a few years it increased again. Several independent studies also favor this assumption. With this assumption, a diffusion-dominated dynamo is able to reproduce many important features of the Maunder minimum remarkably well. If the dynamo is in a diffusion-dominated regime, then a slower meridional circulation means that the poloidal field gets more time to diffuse during its transport through the convection zone, making the dynamo weaker. This consequence helps to model both the cycle amplitudes and the Maunder-like minima. We, however, fail to reproduce these results if the dynamo is in an advection-dominated regime.

On the Large-Scale Diffuse Magnetic Field of the Sun

Although the sunspots migrate towards the equator, the large-scale weak diffuse magnetic fields of the Sun migrate poleward with the solar cycle, the polar field reversing at the time of the sunspot maxima. We apply the vector model of Dikpati and Choudhuri (1994, Paper I) to fit these observations. The dynamo layer at the base of the convection zone is taken to be the source of the diffuse field, which is then evolved in the convection zone subject to meridional circulation and turbulent diffusion. We find that the longitudinally averaged observational data can be fitted reasonably well both for positive and negative values of the alpha-effect by adjusting the subsurface meridional flow suitably. The model will be extended in a future paper to include the decay of active regions as an extra source of the diffuse field, which may be necessary to explain the probable phase lag between B-tau and B-phi at lower latitudes.

Quenching of Meridional Circulation in Flux Transport Dynamo Models

Guided by the recent observational result that the meridional circulation of the Sun becomes weaker at the time of the sunspot maximum, we have included a parametric quenching of the meridional circulation in solar dynamo models such that the meridional circulation becomes weaker when the magnetic field at the base of the convection zone is stronger. We find that a flux transport solar dynamo tends to become unstable on including this quenching of meridional circulation if the diffusivity in the convection zone is less than about 2x10(11) cm(2) s(-1). The quenching of alpha, however, has a stabilizing effect and it is possible to stabilize a dynamo with low diffusivity with sufficiently strong alpha-quenching. For dynamo models with high diffusivity, the quenching of meridional circulation does not produce a large effect and the dynamo remains stable. We present a solar-like solution from a dynamo model with diffusivity 2.8x10(12) cm(2) s(-1) in which the quenching of meridional circulation makes the meridional circulation vary periodically with solar cycle as observed and does not have any other significant effect on the dynamo.

Can Superflares Occur on Our Sun?

Recent observations of Sun-like stars, similar to our Sun in their surface temperature (5600-6000 K) and slow rotation (rotational period > 10 d), using the Kepler satellite by Maehara et al. (2012, Nature, 485, 478) have revealed the existence of superflares (with energy of 10(33)-10(35) erg). From statistical analyses of these superflares, it was found that superflares with energy of 10(34) erg occur once in 800 yr, and superflares with 10(35) erg occur once in 5000 yr. In this paper, we examine whether superflares with energy of 10(33)-10(35) erg could occur on the present Sun through the use of simple order-of-magnitude estimates based on current ideas related to the mechanisms of the solar dynamo. If magnetic flux is generated by differential rotation at the base of the convection zone, as assumed in typical dynamo models, it is possible that the present Sun would generate a large sunspot with a total magnetic flux of similar to 2 x 10(23) Mx (= G cm(2)) within one solar cycle period, and lead to superflares with an energy of 10(34) erg. To store a total magnetic flux of similar to 10(24) Mx, necessary for generating 10(35) erg superflares, it would take similar to 40 yr. Hot Jupiters have often been argued to be a necessary ingredient for the generation of superflares, but we found that they do not play any essential role in the generation of magnetic flux in the star itself, if we consider only the magnetic interaction between the star and the hot Jupiter. This seems to be consistent with Maehara et al.'s finding of 148 superflare-generating solar-type stars that do not have a hot Jupiter-like companion. Altogether, our simple calculations, combined with Maehara et al.'s analysis of superflares on Sun-like stars, show that there is a possibility that superflares of 10(34) erg would occur once in 800 yr on our present Sun.