Towards Scalable Gaussian Process Modeling