Spatial processes exhibit nonstationarity in many climate and environmental applications. Convolution-based approaches are often used to construct nonstationary covariance functions in Gaussian processes. Although convolution-based models are flexible, their computation is extremely expensive when the data set is large. Most existing methods rely on fitting an anisotropic, but stationary model locally, and then reconstructing the spatially varying parameters. In this study, we propose a new estimation procedure to approximate a class of nonstationary Matérn covariance functions by local-polynomial fitting the covariance parameters. The proposed method allows for efficient estimation of a richer class of nonstationary covariance functions, with the local stationary model as a special case. We also develop an approach for a fast high-resolution simulation with nonstationary features on a small scale and apply it to precipitation data in climate model outputs.
|Original language||English (US)|
|Number of pages||23|
|State||Published - Oct 8 2018|
Bibliographical noteKAUST Repository Item: Exported on 2020-10-01
Acknowledgements: This research was supported by funding from King Abdullah University of Science and Technology (KAUST). We would like to thank the editor, associate editor, and two anonymous reviewers for their valuable comments.