Monte Carlo likelihood inference in multivariate model-based geostatistics