We consider the Boussinesq PDE perturbed by a time-dependent forcing. Even though there is no smoothing effect for arbitrary smooth initial data, we are able to apply the method of self-consistent bounds to deduce the existence of smooth classical periodic solutions in the vicinity of 0. The proof is non-perturbative and relies on construction of periodic isolating segments in the Galerkin projections.