Scalable Gaussian-process regression and variable selection using Vecchia approximations