Latin Hypercube Sampling of Gaussian Random field for Sobol’s Global Senitivity Analysis of Models with Spatial Inputs and Scalar Outputs

Nathalie Saint-Geours1, Jean-Stephane Bailly1, Frederic Grelot2, Christian Lavergne3
1. AgroParisTech, UMR TETISF-34093, Montpellier, France
2. Cemagref, UMR G-EAUF-34093, Montpellier, France
3. Mathematics and Modelling Institute of Montpellier, Montpellier 3 University, Montpellier, France
1.{saintge, bailly}@teledetection.fr; 2. frederic.grelot@cemagref.fr; 3. Christian.lavergne@univ-montp3.fr

Abstract: The variance-based Sobol' approach is one of the few global sensitivity analysis methods that is suitable for complex models with spatially distributed inputs. Yet it needs a large number of model runs to compute sensitivity indices: in the case of models where some inputs are 2D Gaussian random fields, it is of great importance to generate a relatively small set of map realizations capturing most of the variability of the spatial inputs. The purpose of this paper is to discuss the use of Latin Hypercube Sampling (LHS) of geostatistical simulations to reach better efficiency in the computation of Sobol' sensitivity indices on spatial models. Sensitivity indices are estimated on a simple analytical model with a spatial input, for increasing sample size, using either Simple Random Sampling (SRS) or LHS to generate input map realizations. Results show that using LHS rather than SRS yields sensitivity indices estimates which are slightly more precise (smaller variance), with no significant improvement of bias.

Keywords: global sensitivity analysis; Latin Hypercube Sampling; Gaussian random field; unconditionnal simulation

AttachmentSize
GeoursAccuracy2010.pdf517.96 KB