We introduce a regularization technique for the approximate reconstruction of spatial and time varying source terms using the observed solutions of the forward time fractional diffusion problem on a discrete set of points. The numerical method is based on computation of the derivatives of adaptive filtered versions of the noisy data by discrete mollification.