We extend the variational multiscale proper orthogonal decomposition reduced order modeling (VMSPOD) to flows governed by double diffusive convection. We present stability and convergence analysis for it, and give results ...