SMS dit: Scalable gaussian processes for characterizing multidimensional change surfaces