The Theory and Practice of Highly Scalable Gaussian Process Regression with Nearest Neighbours