Context. Dynamos in the Sun and other bodies tend to produce magnetic fields that possess magnetic helicity of opposite sign at large and small scales, respectively. The build-up of magnetic helicity at small scales provides an important saturation mechanism. Aims. In order to understand the nature of the solar dynamo we need to understand the details of the saturation mechanism in spherical geometry. In particular, we aim to understand the effects of magnetic helicity fluxes from turbulence and meridional circulation. Methods. We consider a model with only radial shear confined to a thin layer (tachocline) at the bottom of the convection zone. The kinetic α owing to helical turbulence is assumed to be localized in a region above the convection zone. The dynamical quenching formalism is used to describe the build-up of mean magnetic helicity in the model, which results in a magnetic α effect that feeds back on the kinetic α effect. In some cases we compare these results with those obtained from a model with a simple algebraic α quenching formula.Results. In agreement with earlier findings, the magnetic α effect has the opposite sign compared with the kinetic α effect and leads to a catastrophic decrease of the saturation field strength proportional to the inverse magnetic Reynolds number. At high latitudes this quenching effect can lead to secondary dynamo waves that propagate poleward because of the opposite sign of α. These secondary dynamo waves are driven by small-scale magnetic helicity instead of the small-scale kinetic helicity. Magnetic helicity fluxes both from turbulent mixing and from meridional circulation alleviate catastrophic quenching. Interestingly, supercritical diffusive helicity fluxes also give rise to secondary dynamo waves and grand minima-like episodes.