Scalable Gaussian process models for changepoint detection and spatio-temporal predictions with large correlated data