Orthogonal representations of random fields and an application to geophysics data