Multivariate isotropic random fields on spheres: Nonparametric Bayesian modeling and Lp fast approximations