Learning-based Robust Stabilization for Reduced-Order Models of 2D and 3D Boussinesq Equations