Gaussian random fields are completely characterised by their mean value and covariance function. Random fields on hyperbolic spaces have been studied to a limited extent only, namely for the case of scalar-valued fields that are not evolving over time. This paper challenges the problem of the second-order characteristics of multivariate (vector-valued) random fields that evolve temporally over hyperbolic spaces. Specifically, we characterise the continuous space–time covariance functions that are isotropic (radially symmetric) over space (the hyperbolic space) and stationary over time (the real line). Our finding is the analogue of recent findings that have been shown for the case where the space is either the n-dimensional sphere or more generally a two-point homogeneous space. Our main result can be read as a spectral representation theorem, and we also detail the main result for the subcase of covariance functions having a spectrum that is absolutely continuous with respect to the Lebesgue measure (technical details are reported below).