Bayesian Quadrature on Riemannian Data Manifolds