We develop a new approach for detecting anomalies in the behavior of stochastic processes and random fields. The approach uses tensor product representations, in particular the multivariate Karhunen–Loève (KL) expansion on complex domains. From the associated eigenspaces of the covariance operator a series of nested function spaces are constructed, allowing detection of signals lying in orthogonal subspaces. In particular this can succeed even if the stochastic behavior of the signal changes either in a global or local sense. A mathematical approach is developed to locate and measure sizes of extraneous components based on construction of multilevel nested subspaces. We show examples in R and on a spherical domain S2. However, the method is flexible, allowing the detection of orthogonal signals on general topologies, including spatio-temporal domains.