Spectral representations uniquely define the covariance functions associated to random fields defined over spheres or spheres cross time. Covariance functions on spheres cross time are usually modeled under the assumptions of either spatial isotropy or axial symmetry, and the assumption of temporal stationarity. This paper goes beyond these assumptions. In particular, we consider the problem of spatially anisotropic covariance functions on spheres. The crux of our criterion is to escape from the addition theorem for spherical harmonics. We also challenge the problem of temporal nonstationarity in nonseparable space-time covariance functions, where space is the n-dimensional sphere. (C) 2019 Elsevier B.V. All rights reserved.