Can short time delays influence the variability of the solar cycle?

We present the effects of introducing results of 3D MHD simulations of buoyant magnetic fields in the solar convection zone in 2D mean-field Babcock-Leighton models. In particular, we take into account the time delay introduced by the rise time of the toroidal structures from the base of the convection zone to the solar surface. We find that the delays produce large temporal modulation of the cycle amplitude even when strong and thus rapidly rising flux tubes are considered. The study of a reduced model reveals that aperiodic modulations of the solar cycle appear after a sequence of period doubling bifurcations typical of non-linear systems. We also discuss the memory of such systems and the conclusions which may be drawn concerning the actual solar cycle variability.