Gaussian Process Subspace Regression for Model Reduction