Zaslat SMS: Scalable gaussian processes for characterizing multidimensional change surfaces