A posteriori error analysis of an augmented mixed-primal formulation for the stationary Boussinesq model