Summary: We study the large time behavior of solutions to a generalization of the Boussinesq system of equations in \(n \geq 2\) spatial dimensions. We establish the existence and algebraic decay of the \(L^2\)-norm of the solution.