Mat\'ern Gaussian processes on Riemannian manifolds