In this work we introduce and analyze a mixed virtual element method (mixed-VEM) for the two-dimensional stationary Boussinesq problem. The continuous formulation is based on the introduction of a pseudostress tensor ...