Stav dette: Scalable gaussian processes for characterizing multidimensional change surfaces