Spatially Aggregated Gaussian Processes with Multivariate Areal Outputs