Optimal adaptive sequential spatial sampling of soil using pair-copulas

A spatial sampling design that uses pair-copulas is presented that aims to reduce prediction uncertainty by selecting additional sampling locations based on both the spatial configuration of existing locations and the values of the observations at those locations. The novelty of the approach arises in the use of pair-copulas to estimate uncertainty at unsampled locations. Spatial pair-copulas are able to more accurately capture spatial dependence compared to other types of spatial copula models. Additionally, unlike traditional kriging variance, uncertainty estimates from the pair-copula account for influence from measurement values and not just the configuration of observations. This feature is beneficial, for example, for more accurate identification of soil contamination zones where high contamination measurements are located near measurements of varying contamination. The proposed design methodology is applied to a soil contamination example from the Swiss Jura region. A partial redesign of the original sampling configuration demonstrates the potential of the proposed methodology.

[1]  Edzer Pebesma,et al.  The pair-copula construction for spatial data: a new approach to model spatial dependency , 2011 .

[2]  João Felipe Coimbra Leite Costa,et al.  Additional Samples: Where They Should Be Located , 2001 .

[3]  Ute A. Mueller,et al.  Direct minimum/maximum autocorrelation factors within the framework of a two structure linear model of coregionalisation , 2008, Comput. Geosci..

[4]  Dale L. Zimmerman,et al.  Optimal network design for spatial prediction, covariance parameter estimation, and empirical prediction , 2006 .

[5]  Timothy C. Coburn,et al.  Geostatistics for Natural Resources Evaluation , 2000, Technometrics.

[6]  Mike Rees,et al.  5. Statistics for Spatial Data , 1993 .

[7]  Pierre Goovaerts,et al.  Evaluating the probability of exceeding a site-specific soil cadmium contamination threshold , 2001 .

[8]  B. Gräler Modelling skewed spatial random fields through the spatial vine copula , 2014 .

[9]  R. Lark Optimized spatial sampling of soil for estimation of the variogram by maximum likelihood , 2002 .

[10]  João Felipe Coimbra Leite Costa,et al.  Choosing Between Two Kind of Sampling Patterns Using Geostatistical Simulation: Regularly Spaced or at High Uncertainty Locations? , 2011 .

[11]  J. W. Groenigen,et al.  Constrained optimisation of soil sampling for minimisation of the kriging variance , 1999 .

[12]  R. Lark,et al.  Optimized multi-phase sampling for soil remediation surveys , 2013 .

[13]  A. Bárdossy Copula‐based geostatistical models for groundwater quality parameters , 2006 .

[14]  A. Journel,et al.  Resampling from stochastic simulations , 1994, Environmental and Ecological Statistics.

[15]  A. Bárdossy,et al.  Geostatistical interpolation using copulas , 2008 .

[16]  Dominique Arrouays,et al.  Spatial prediction of soil properties with copulas , 2011 .

[17]  Jing Li,et al.  A Copula based observation network design approach , 2011, Environ. Model. Softw..

[18]  R. Webster,et al.  Sample adequately to estimate variograms of soil properties , 1992 .

[19]  R. Webster,et al.  Geostatistical analysis of soil contamination in the Swiss Jura. , 1994, Environmental pollution.

[20]  M. Stein,et al.  Spatial sampling design for prediction with estimated parameters , 2006 .

[21]  A. Bárdossy,et al.  Application of Copulas in Geostatistics , 2010 .

[22]  A. Shapiro Monte Carlo Sampling Methods , 2003 .

[23]  James V. Zidek,et al.  Designing environmental monitoring networks to measure extremes , 2007, Environmental and Ecological Statistics.

[24]  Peter J. Diggle,et al.  Bayesian Geostatistical Design , 2006 .