Clustering the Nearest Neighbor Gaussian Process

Gaussian processes are ubiquitous as the primary tool for modeling spatial data. However, the Gaussian process is limited by its 𝒪(n3) cost, making direct parameter fitting algorithms infeasible for the scale of modern data collection initiatives. The Nearest Neighbor Gaussian Process (NNGP) was introduced as a scalable approximation to dense Gaussian processes which has been successful for n∼106$$ n\sim 1{0}^6 $$ observations. This project introduces the clustered Nearest Neighbor Gaussian Process (cNNGP) which reduces the computational and storage cost of the NNGP for stationary and isotropic datasets. The accuracy of parameter estimation and reduction in computational and memory storage requirements are demonstrated with simulated data, where the cNNGP provided comparable inference to that obtained with the NNGP, in a fraction of the sampling time. An extensive simulation study is presented, with cNNGP compared with similar contemporary methods. To showcase the method's performance, we modeled biomass over the state of Maine using data collected by the Global Ecosystem Dynamics Investigation (GEDI) to generate wall‐to‐wall predictions over the state. In 20% of the time, the cNNGP produced nearly indistinguishable inference and biomass prediction maps to those obtained with the NNGP.