A hierarchical multivariate spatio-temporal model for clustered climate data with annual cycles