A latent variable approach to modelling multivariate geostatistical skew-normal data