Spatial Warped Gaussian Processes: Estimation and Efficient Field Reconstruction

A class of models for non-Gaussian spatial random fields is explored for spatial field reconstruction in environmental and sensor network monitoring. The family of models explored utilises a class of transformation functions known as Tukey g-and-h transformations to create a family of warped spatial Gaussian process models which can support various desirable features such as flexible marginal distributions, which can be skewed, leptokurtic and/or heavy-tailed. The resulting model is widely applicable in a range of spatial field reconstruction applications. To utilise the model in applications in practice, it is important to carefully characterise the statistical properties of the Tukey g-and-h random fields. In this work, we study both the properties of the resulting warped Gaussian processes as well as using the characterising statistical properties of the warped processes to obtain flexible spatial field reconstructions. In this regard we derive five different estimators for various important quantities often considered in spatial field reconstruction problems. These include the multi-point Minimum Mean Squared Error (MMSE) estimators, the multi-point Maximum A-Posteriori (MAP) estimators, an efficient class of multi-point linear estimators based on the Spatial-Best Linear Unbiased (S-BLUE) estimators, and two multi-point threshold exceedance based estimators, namely the Spatial Regional and Level Exceedance estimators. Simulation results and real data examples show the benefits of using the Tukey g-and-h transformation as opposed to standard Gaussian spatial random fields in a real data application for environmental monitoring.

Automated summary

A class of non-Gaussian spatial random field models using Tukey g-and-h transformations for spatial field reconstruction was explored. The resulting warped spatial Gaussian process models support various desirable features and have wide applicability. Statistical properties of the models were carefully characterized to obtain flexible spatial field reconstructions, deriving five different estimators for important quantities in spatial field reconstruction problems. Simulation results and real data examples demonstrated the benefits of using Tukey g-and-h transformations over standard Gaussian spatial random fields in environmental monitoring applications.