Multi-resolution Representation of Gaussian Fields with Discontinuities