This article deals with Neumann boundary optimal control problems associated with the Boussinesq equations including solid media. These problems are rst put into an appropriate mathematical formulation. Then the existence of optimal solutions is proved. The use of Lagrange multiplier techniques is justified and an optimality system of equations is derived.