Seol mar théacs é seo: Scalable gaussian processes for characterizing multidimensional change surfaces