Matérn Gaussian processes on Riemannian manifolds