Neste trabalho, propomos um modelo espaço-temporal matriz-variado para ajustar duas ou mais variáveis respostas, mensuradas em diferentes pontos no tempo e locais distintos no espaço. Aqui relaxamos a hipótese de isotropia no espaço, uma vez que esta suposição não é razoável em várias situações reais. Tratamos a anisotropia do processo espacial usando o conceito de deformação espacial e lidamos com a estrutura temporal via modelos dinâmicos bayesianos. Através de estudos de simulação e também com dois exemplos, nós mostramos que a incorporação da deformação espacial ao modelo proporciona melhores interpolações em alguns cenários anisotrópicos. Apresentamos também uma solução para lidar com valores faltantes em diferentes variáveis respostas.