In this paper, we study two inverse problems for the nonlinear Boussinesq system for incompressible viscoelastic non-isothermal Kelvin-Voigt fluids. The studying inverse problems consist of determining an intensities of density of external forces and heat source under given integral overdetermination conditions. Two types of boundary conditions for the velocity v are considered: sticking and sliding conditions on boundary. In both cases of these boundary conditions, the local and global in time existence and uniqueness of weak and strong solutions are established under suitable assumptions on the data. The large time behavior of weak solutions is also studied.