Scalable Bayesian reduced-order models for high-dimensional multiscale dynamical systems