  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.
