A matrix-variate dirichlet process to model earthquake hypocentre temporal patterns. (August 2022)