Probabilistic inference of subsurface heterogeneity and interface geometry using geophysical data