Gaussian processes are efficiently approximated using GPU-accelerated Vecchia methods for geospatial data analysis.