Covariance tapering for Gaussian multivariate random fields