Stability and optimal decay estimates for the 3D anisotropic Boussinesq equations