Multivariate spatio-temporal models for high-dimensional areal data with application to Longitudinal Employer-Household Dynamics