Estimation of Space Deformation Model for Non-stationary Random Functions

TL;DR

This paper introduces a novel method for estimating space deformation in non-stationary random functions using a single realization, combining kernel smoothing, MDS, and spline functions, outperforming stationary models in accuracy.

Contribution

It proposes a new approach for non-stationary geostatistical modeling from a single realization, advancing beyond previous multi-realization methods.

Findings

Method accurately retrieves true deformation on simulated data.

Outperforms stationary approaches on synthetic and real datasets.

Enables exploratory analysis of non-stationarity.

Abstract

Stationary Random Functions have been successfully applied in geostatistical applications for decades. In some instances, the assumption of a homogeneous spatial dependence structure across the entire domain of interest is unrealistic. A practical approach for modelling and estimating non-stationary spatial dependence structure is considered. This consists in transforming a non-stationary Random Function into a stationary and isotropic one via a bijective continuous deformation of the index space. So far, this approach has been successfully applied in the context of data from several independent realizations of a Random Function. In this work, we propose an approach for non-stationary geostatistical modelling using space deformation in the context of a single realization with possibly irregularly spaced data. The estimation method is based on a non-stationary variogram kernel estimator which serves as a dissimilarity measure between two locations in the geographical space. The proposed procedure combines aspects of kernel smoothing, weighted non-metric multi-dimensional scaling and thin-plate spline radial basis functions. On a simulated data, the method is able to retrieve the true deformation. Performances are assessed on both synthetic and real datasets. It is shown in particular that our approach outperforms the stationary approach. Beyond the prediction, the proposed method can also serve as a tool for exploratory analysis of the non-stationarity.