Extreme solar activity fluctuations and the occurrence of solar grand minima and maxima episodes, are well established, observed features of the solar cycle. Nevertheless, such extreme activity fluctuations and the dynamics of the solar cycle during Maunder minima-like episodes remain ill-understood. We explore the origin of such extreme solar activity fluctuations and the role of dual poloidal field sources, namely the Babcock-Leighton mechanism and the mean-field alpha effect in the dynamics of the solar cycle. We mainly concentrate on entry and recovery from grand minima episodes such as the Maunder minimum and the dynamics of the solar cycle. We use a kinematic solar dynamo model with a novel set-up in which stochastic perturbations force two distinct poloidal field alpha effects. We explore different regimes of operation of these poloidal sources with distinct operating thresholds, to identify the importance of each. The perturbations are implemented independently in both hemispheres which allows one to study the level of hemispheric coupling and hemispheric asymmetry in the emergence of sunspots. From the simulations performed we identify a few different ways in which the dynamo can enter a grand minima episode. While fluctuations in any of the $\alpha$ effects can trigger intermittency we find that the mean-field alpha effect is crucial for the recovery of the solar cycle from a grand minima episode which a Babcock-Leighton source alone, fails to achieve. Our simulations also demonstrate other cycle dynamics. We conclude that stochastic fluctuations in two interacting poloidal field sources working with distinct operating thresholds is a viable candidate for triggering episodes of extreme solar activity and that the mean-field alpha effect capable of working on weak, sub-equipartition fields is critical to the recovery of the solar cycle following an extended solar minimum.

Surface observations indicate that the speed of the solar meridional circulation in the photosphere varies in anti-phase with the solar cycle. The current explanation for the source of this variation is that inflows into active regions alter the global surface pattern of the meridional circulation. When these localized inflows are integrated over a full hemisphere, they contribute to the slow down of the axisymmetric poleward horizontal component. The behavior of this large scale flow deep inside the convection zone remains largely unknown. Present helioseismic techniques are not sensitive enough to capture the dynamics of this weak large scale flow. Moreover, the large time of integration needed to map the meridional circulation inside the convection zone, also masks some of the possible dynamics on shorter timescales. In this work we examine the dynamics of the meridional circulation that emerges from a 3D MHD global simulation of the solar convection zone. Our aim is to assess and quantify the behavior of meridional circulation deep inside the convection zone, where the cyclic large-scale magnetic field can reach considerable strength. Our analyses indicate that the meridional circulation morphology and amplitude are both highly influenced by the magnetic field, via the impact of magnetic torques on the global angular momentum distribution. A dynamic feature induced by these magnetic torques is the development of a prominent upward flow at mid latitudes in the lower convection zone that occurs near the equatorward edge of the toroidal bands and that peaks during cycle maximum. Globally, the dynamo-generated large-scale magnetic field drives variations in the meridional flow, in stark contrast to the conventional kinematic flux transport view of the magnetic field being advected passively by the flow.

The form of the solar meridional circulation is a very important ingredient for mean field flux transport dynamo models. Yet a shroud of mystery still surrounds this large-scale flow, given that its measurement using current helioseismic techniques is challenging. In this work we use results from 3D global simulations of solar convection to infer the dynamical behavior of the established meridional circulation. We make a direct comparison between the meridional circulation that arises in these simulations and the latest observations. Based on our results we argue that there should be an equatorward flow at the base of the convection zone at mid latitudes, below the current maximum depth helioseismic measures can probe (0.75 R). We also provide physical arguments to justify this behaviour. The simulations indicate that the meridional circulation undergoes substantial changes in morphology as the magnetic cycle unfolds. We close by discussing the importance of these dynamical changes for current methods of observation that involve long averaging periods of helioseismic data. Also noteworthy is the fact that these topological changes indicate a rich interaction between magnetic fields and plasma flows, which challenges the ubiquitous kinematic approach used in the vast majority of mean field dynamo simulations.

Numerical simulations that reproduce solar-like magnetic cycles can be used to generate long-term statistics. The variations in N-S hemispheric cycle synchronicity and amplitude produced in simulations has not been widely compared to observations. The observed limits on asymmetry show that hemispheric sunspot area production is no more than 20% asymmetric for cycles 12-23 and phase lags do not exceed 20% (2 yrs) of the total cycle period. Independent studies have found a long-term trend in phase values as one hemisphere leads the other for ~four cycles. Such persistence in phase is not indicative of a stochastic phenomenon. We compare the findings to results from a numerical simulation of solar convection recently produced with the EULAG-MHD model. This simulation spans 1600 yrs and generated 40 regular, sunspot-like cycles. While the simulated cycle length is too long and the toroidal bands remain at too high of latitudes, some solar-like aspects of hemispheric asymmetry are reproduced. The model reproduces the synchrony of polarity inversions and onset of cycle as the simulated phase lags do not exceed 20% of the cycle period. Simulated amplitude variations between the N and S hemispheres are larger than observed in the Sun. The simulations show one hemisphere persistently leads the other for several successive cycles, placing an upper bound on the efficiency of transequatorial magnetic coupling mechanisms. These include magnetic diffusion, cross-equatorial mixing within elongated convective rolls and transequatorial meridional flow cells. One or more of these processes may lead to magnetic flux cancellation whereby the oppositely directed fields come in close proximity and cancel each other across the magnetic equator late in the solar cycle. We discuss the discrepancies between model and observations and the constraints they pose on possible mechanisms of hemispheric coupling.

This article reviews some of the leading results obtained in solar dynamo physics by using temporal oscillator models as a tool to interpret observational data and dynamo model predictions. We discuss how solar observational data such as the sunspot number is used to infer the leading quantities responsible for the solar variability during the last few centuries. Moreover, we discuss the advantages and difficulties of using inversion methods (or backward methods) over forward methods to interpret the solar dynamo data. We argue that this approach could help us to have a better insight about the leading physical processes responsible for solar dynamo, in a similar manner as helioseismology has helped to achieve a better insight on the thermodynamic structure and flow dynamics in the Sun's interior.

Fluctuations in the Sun's magnetic activity, including episodes of grand minima such as the Maunder minimum have important consequences for space and planetary environments. However, the underlying dynamics of such extreme fluctuations remain ill-understood. Here we use a novel mathematical model based on stochastically forced, non-linear delay differential equations to study solar cycle fluctuations, in which, time delays capture the physics of magnetic flux transport between spatially segregated dynamo source regions in the solar interior. Using this model we explicitly demonstrate that the Babcock-Leighton poloidal field source based on dispersal of tilted bipolar sunspot flux, alone, can not recover the sunspot cycle from a grand minimum. We find that an additional poloidal field source effective on weak fields--the mean-field alpha-effect driven by helical turbulence--is necessary for self-consistent recovery of the sunspot cycle from grand minima episodes.

We develop and subsequently explore the solution space of a simple flux transport dynamo model that incorporates a time dependent large scale meridional circulation. Based on recent observations we prescribed an analytical form for the amplitude of this circulation and study its impact in the evolution of the magnetic field. We find that cyclic variations in the amplitude and frequency of the meridional flow affect the strength of the solar cycle. Variations in the amplitude of the fluctuations influence the shape of the solar cycle but are only relevant to the cycle's strength variations when they occur at a frequency different from or out of phase of the solar cycle's.

In this work we use a low order dynamo model and study under which conditions can it reproduce solar grand minima. We begin by building the phase space of a proxy for the toroidal component of the solar magnetic field and we develop a model, derived from mean field dynamo theory, that gives the time evolution of the toroidal field. This model is characterized by a non-linear oscillator whose coefficients retain most of the physics behind dynamo theory. In the derivation of the model we also include stochastic oscillations in the $\alpha$ effect. We found no evidences that stochastic fluctuations in a linear $\alpha$ effect can trigger grand minima episodes in this model. In contrast, the model used points out that possible mechanism that can trigger grand minima should involve the meridional circulation, magnetic diffusivity or field intensification by buoyancy driven instabilities.