In this paper, the focus is on the reconstruction of a diffusive field and the localization of the underlying driving sources on arbitrary graphs by observing a significantly smaller subset of vertices of the graph uniformly in time. Specifically, we focus on the heat diffusion equation driven by an initial field and an external time-invariant input. When the underlying driving sources are modeled as an initial field or external input, the sources (hence the diffusive field) can be recovered from the subsampled observations without imposing any band-limiting or sparsity constraints. When the diffusion is induced by both the initial field and external input, then the field and sources can be recovered from the subsampled observations, however, by imposing band-limiting constraints on either the initial field or external input. For heat diffusion on graphs, we can compensate for the unobserved vertices with the temporal samples at the observed vertices. If the observations are noiseless, then the recovery is exact. Nonetheless, the developed least squares estimators perform reasonably well with noisy observations. We apply the developed theory for localizing and recovering hot spots on a rectangular metal plate with a cavity.