In this paper we prove the existence and regularity of a solution to a two-dimensional system of evolutionary hemivariational inequalities which describes the Boussinesq model with nonmonotone friction and heat flux. We use the time retardation and regularization technique, combined with a regularized Galerkin method, and recent results from the theory of hemivariational inequalities. (C) 2016 Elsevier Ltd. All rights reserved.