Top panel, left: radial profile of mean-field α-coefficient and the Babcock–Leighton source term due to ephemeral regions. Right: poloidal field line contours obtained from the double-ring algorithm in the northern and southern hemispheres, respectively. Middle panel: butterfly diagram generated from our simulation in the diffusion-dominated region. Bottom panel: typical variation of at the base of the solar convection zone with time.