The large-scale magnetic field in the Sun varies with a period of approximately 22 years, although the amplitude of the cycle is subject to long-term modulation with recurrent phases of significantly reduced magnetic activity. It is believed that a hydromagnetic dynamo is responsible for producing this large-scale field, although this dynamo process is not well understood. Within the framework of mean-field dynamo theory, our aim is to investigate how competing mechanisms for poloidal field regeneration (namely a time delayed Babcock-Leighton surface alpha-effect and an interface-type alpha-effect), can lead to the modulation of magnetic activity in a deep-seated solar dynamo model. We solve the standard alpha-omega dynamo equations in one spatial dimension, including source terms corresponding to both of the the competing alpha-effects in the evolution equation for the poloidal field. This system is solved using two different methods. In addition to solving the one-dimensional partial differential equations directly, using numerical techniques, we also use a local approximation to reduce the governing equations to a set of coupled ordinary differential equations (ODEs), which are studied using a combination of analytical and numerical methods. In the ODE model, it is straightforward to find parameters such that a series of bifurcations can be identified as the time delay is increased, with the dynamo transitioning from periodic states to chaotic states via multiply periodic solutions. Similar transitions can be observed in the full model, with the chaotically modulated solutions exhibiting solar-like behaviour.

We discuss asymptotic solutions of the kinematic $\alpha\omega$-dynamo in a thin disc (slab). Focusing upon the strong dynamo regime, in which the dynamo number $D$ satisfies $|D|\gg1$, we resolve uncertainties in the earlier treatments and conclude that some of the simplifications that have been made in previous studies are questionable. Comparing numerical solutions with asymptotic results obtained for $|D|\gg1$ and $|D|\ll1$ we find that the asymptotic solutions give a reasonably accurate description of the dynamo even far beyond their formal ranges of applicability. Indeed, our results suggest a simple analytical expression for the growth rate of the mean magnetic field that remains accurate in the range $-200< D< -10$ (which is appropriate for dynamos in spiral galaxies and accretion discs). Finally, we analyse the role of various terms in the dynamo equations to clarify the fine details of the dynamo process.