We consider to what extent the long-term dynamics of cyclic solar activity in the form of Grand Minima can be associated with random fluctuations of the parameters governing the solar dynamo. We consider fluctuations of the alpha-coefficient in the conventional Parker migratory dynamo, and also in slightly more sophisticated dynamo models, and demonstrate that they can mimic the gross features of the phenomenon of the occurrence of Grand Minima over a suitable parameter range. The temporal distribution of these Grand Minima appears chaotic, with a more or less exponential waiting time distribution, typical of Poisson processes. In contrast however, the available reconstruction of Grand Minima statistics based on cosmogenic isotope data demonstrates substantial deviations from this exponential law. We were unable to reproduce the non-Poissonic tail of the waiting time distribution either in the framework of a simple alpha-quenched Parker model, or in its straightforward generalization, nor in simple models with feedback on the differential rotation. We suggest that the disagreement may only be apparent and is plausibly related to the limited observational data, and that the observations and results of numerical modeling can be consistent and represent physically similar dynamo regimes.

A comparison between the two tracers of magnetic field mirror asymmetry in solar active regions, twist and current helicity, is presented. It is shown that for individual active regions these tracers do not possess visible similarity while averaging by time over the solar cycle, or by latitude, reveals similarities in their behaviour. The main property of the dataset is anti-symmetry over the solar equator. Considering the evolution of helical properties over the solar cycle we find signatures of a possible sign change at the beginning of the cycle, though more systematic observational data are required for a definite confirmation. We discuss the role of both tracers in the context of the solar dynamo theory.