The solar dipole moment at activity minimum is a good predictor of the strength of the subsequent solar cycle. Through a systematic analysis using a state-of-the-art 2$\times$2 D solar dynamo model, we found that bipolar magnetic regions (BMR) with atypical characteristics can modify the strength of the next cycle via their impact on the buildup of the dipole moment as a sunspot cycle unfolds. In addition to summarizing these results, we present further effects of such "rogue" BMRs. These have the ability to generate hemispheric asymmetry in the subsequent sunspot cycle, since they modify the polar cap flux asymmetry of the ongoing cycle. We found strong correlation between the polar cap flux asymmetry of cycle i and the total pseudo sunspot number asymmetry of cycle i + 1 . Good correlation also appears in the case of the time lag of the hemispheres of cycle i + 1.

The origin of cycle-to-cycle variations in solar activity is currently the focus of much interest. It has recently been pointed out that large individual active regions with atypical properties can have a significant impact on the long term behaviour of solar activity. We investigate this possibility in more detail using a recently developed 2$\times$2D dynamo model of the solar magnetic cycle. We find that even a single "rogue" bipolar magnetic region (BMR) in the simulations can have a major effect on the further development of solar activity cycles, boosting or suppressing the amplitude of subsequent cycles. In extreme cases an individual BMR can completely halt the dynamo, triggering a grand minimum. Rogue BMRs also have the potential to induce significant hemispheric asymmetries in the solar cycle. To study the effect of rogue BMRs in a more systematic manner, a series of dynamo simulations were conducted, in which a large test BMR was manually introduced in the model at various phases of cycles of different amplitudes. BMRs emerging in the rising phase of a cycle can modify the amplitude of the ongoing cycle while BMRs emerging in later phases will only impact subsequent cycles. In this model, the strongest impact on the subsequent cycle occurs when the rogue BMR emerges around cycle maximum at low latitudes but the BMR does not need to be strictly cross-equatorial. Active regions emerging as far as $20^{\circ}$ from the equator can still have a significant impact. We demonstrate that the combined effect of the magnetic flux, tilt angle and polarity separation of the BMR on the dynamo is via their contribution to the dipole moment, $\delta D_{\mathrm{BMR}}$. Our results indicate that prediction of the amplitude, starting epoch and duration of a cycle requires an accurate accounting of a broad range of active regions emerging in the previous cycle.

In this paper we complete the presentation of a new hybrid $2\times2$D flux transport dynamo (FTD) model of the solar cycle based on the Babcock-Leighton mechanism of poloidal magnetic field regeneration via the surface decay of bipolar magnetic regions (BMRs). This hybrid model is constructed by allowing the surface flux transport (SFT) simulation described in Lemerle et al. 2015 to provide the poloidal source term to an axisymmetric FTD simulation defined in a meridional plane, which in turn generates the BMRs required by the SFT. A key aspect of this coupling is the definition of an emergence function describing the probability of BMR emergence as a function of the spatial distribution of the internal axisymmetric magnetic field. We use a genetic algorithm to calibrate this function, together with other model parameters, against observed cycle 21 emergence data. We present a reference dynamo solution reproducing many solar cycle characteristics, including good hemispheric coupling, phase relationship between the surface dipole and the BMR-generating internal field, and correlation between dipole strength at cycle maximum and peak amplitude of the next cycle. The saturation of the cycle amplitude takes place through the quenching of the BMR tilt as a function of the internal field. The observed statistical scatter about the mean BMR tilt, built into the model, acts as a source of stochasticity which dominates amplitude fluctuations. The model thus can produce Dalton-like epochs of strongly suppressed cycle amplitude lasting a few cycles and can even shut off entirely following an unfavorable sequence of emergence events.

The need for reliable predictions of the solar activity cycle motivates the development of dynamo models incorporating a representation of surface processes sufficiently detailed to allow assimilation of magnetographic data. In this series of papers we present one such dynamo model, and document its behavior and properties. This first paper focuses on one of the model's key components, namely surface magnetic flux evolution. Using a genetic algorithm, we obtain best-fit parameters of the transport model by least-squares minimization of the differences between the associated synthetic synoptic magnetogram and real magnetographic data for activity cycle 21. Our fitting procedure also returns Monte Carlo-like error estimates. We show that the range of acceptable surface meridional flow profiles is in good agreement with Doppler measurements, even though the latter are not used in the fitting process. Using a synthetic database of bipolar magnetic region (BMR) emergences reproducing the statistical properties of observed emergences, we also ascertain the sensitivity of global cycle properties, such as the strength of the dipole moment and timing of polarity reversal, to distinct realizations of BMR emergence, and on this basis argue that this stochasticity represents a primary source of uncertainty for predicting solar cycle characteristics.