Multivariate modeling for environmental spatio-temporal data