🤖 AI Summary
Traditional spatial modeling relies on stationarity and isotropy assumptions, which increasingly fail in real-world applications; meanwhile, existing nonstationary approaches often sacrifice either flexibility or computational efficiency. To address this, we propose a covariate-driven, modular covariance function—built upon the Matérn kernel—that enables parameterized, decoupled modeling of multiple sources of nonstationarity, including marginal standard deviation, geometric anisotropy, and smoothness. Our formulation ensures interpretability, computational scalability, and intuitive visualizability. By incorporating a covariate-adaptive construction mechanism and a large-scale adaptation algorithm, the method achieves superior predictive performance over state-of-the-art nonstationary models in both simulation studies and real-world analysis of Swiss precipitation data. It simultaneously attains high model fidelity and low computational overhead, demonstrating robust scalability to large spatial datasets.
📝 Abstract
The assumptions of stationarity and isotropy often stated over spatial processes have not aged well during the last two decades, partly explained by the combination of computational developments and the increasing availability of high-resolution spatial data. While a plethora of approaches have been developed to relax these assumptions, it is often a costly tradeoff between flexibility and a diversity of computational challenges. In this paper, we present a class of covariance functions that relies on fixed, observable spatial information that provides a convenient tradeoff while offering an extra layer of numerical and visual representation of the flexible spatial dependencies. This model allows for separate parametric structures for different sources of nonstationarity, such as marginal standard deviation, geometric anisotropy, and smoothness. It simplifies to a Mat'ern covariance function in its basic form and is adaptable for large datasets, enhancing flexibility and computational efficiency. We analyze the capabilities of the presented model through simulation studies and an application to Swiss precipitation data.