Finite element approximation of the modified Boussinesq equations using a stabilized formulation