Gaussian process product models for nonparametric nonstationarity