Multiple dynamo modes as a mechanism for long-term solar activity variations

Context. Solar magnetic activity shows both smooth secular changes, such as the modern Grand Maximum, and quite abrupt drops that are denoted as grand minima, such as the Maunder Minimum. Direct numerical simulations (DNS) of convection-driven dynamos offer one way of examining the mechanisms behind these events. Aims. In this work, we analyze a solution of a solar-like DNS that was evolved for roughly 80 magnetic cycles of 4.9 years and where epochs of irregular behavior are detected. The emphasis of our analysis is to find physical causes for such behavior. Methods. The DNS employed is a semi-global (wedge-shaped) magnetoconvection model. For the data analysis tasks we use Ensemble Empirical Mode Decomposition and phase dispersion methods, as they are well suited for analyzing cyclic (non-periodic) signals. Results. A special property of the DNS is the existence of multiple dynamo modes at different depths and latitudes. The dominant mode is solar-like (equatorward migration at low latitudes and poleward at high latitudes). This mode is accompanied by a higher frequency mode near the surface and at low latitudes, showing poleward migration, and a low-frequency mode at the bottom of the convection zone. The low-frequency mode is almost purely antisymmetric with respect to the equator, while the dominant mode has strongly fluctuating mixed parity. The overall behavior of the dynamo solution is extremely complex, exhibiting variable cycle lengths, epochs of disturbed and even ceased surface activity, and strong short-term hemispherical asymmetries. Surprisingly, the most prominent suppressed surface activity epoch is actually a global magnetic energy maximum; during this epoch the bottom toroidal magnetic field obtains a maximum, demonstrating that the interpretation of grand minima-type events is non-trivial. The hemispherical asymmetries are seen only in the magnetic field, while the velocity field exhibits considerably weaker asymmetry. Conclusions. We interpret the overall irregular behavior as being due to the interplay of the different dynamo modes showing different equatorial symmetries, especially the smoother part of the irregular variations being related to the variations of the mode strengths, evolving with different and variable cycle lengths. The abrupt low-activity epoch in the dominant dynamo mode near the surface is related to a strong maximum of the bottom toroidal field strength, which causes abrupt disturbances especially in the different