How Gaussian Geostatistical Simulations works

How the realizations are generated

Gaussian Geostatistical Simulations works by first creating a grid of randomly assigned values drawn from a standard normal distribution (mean = 0 and variance = 1). The covariance model (from the semivariogram specified in the simple kriging layer, which is required as input for the simulation) is then applied to the raster. This ensures that raster values conform to the spatial structure found in the input dataset. The resulting raster constitutes one unconditional realization, and many more can be produced using a different raster of normally distributed values each time. Details of this method may be read in Dietrich and Newsam (1993).

If conditional simulation has been selected, the unconditional rasters are conditioned via kriging. This process uses the kriged estimate (prediction) at each location to ensure that the simulated values honor the input data values and that, on average, the kriged predictions are replicated. Details of conditioning the realizations using kriging may be read in Journel (1974).

However, if the simple kriging model includes measurement error, the input data values will not be honored (in the simple kriging layer or in the simulated realizations). In addition to this, the Gaussian Geostatistical Simulations tool has been implemented to use a continuous (or smooth) search neighborhood that avoids discontinuities in the simulated surfaces due to changes in the local neighborhood used in kriging. Details can be found in Aldworth (1998) and Gribov and Krivoruchko (2004).

To read more on geostatistical simulation concepts and for examples of conditional and unconditional simulation, refer to the help section on key concepts of geostatistical simulation.

A general workflow for Gaussian geostatistical simulation involves preparing the data, creating the realizations, back transforming the results into original units, and postprocessing the results and/or using the results as input to a transfer function (model) to assess variability in the model's output. This process is represented in the following figure:

General workflow for Gaussian Geostatistical Simulation.
General workflow for Gaussian geostatistical simulation

Simple kriging models for simulation

The Gaussian Geostatistical Simulations tool accepts any simple kriging model. However, simulation results are only valid if the input data (used to fit the semivariogram and to condition the realizations) is normally distributed. This can be accomplished by applying a normal score transformation. Declustering should be done to obtain a representative histogram from clustered data, and trends should be removed to ensure that the mean is stationary over the spatial domain. These steps are illustrated in the figures below. Figure (a) shows simple kriging with the trend removal, declustering, and normal score transformation options selected. Figure (b) shows the trend that will be removed from the dataset before declustering, normal score transformation, and variography are done. Figure (c) shows declustering by cells. Figure (d) shows the normal score transformation, which in this case uses the direct method.

Simple Kriging with the trend removal, declustering and Normal Score transform options selected.
Simple kriging with the trend removal, declustering, and normal score transformation options selected
Trend that will be removed from the data before declustering and Normal Score transform are applied.
Trend that will be removed from the data before declustering and normal score transformation are applied
Cell declustering of the data after trend removal.
Cell declustering of the data after trend removal
Normal Score transform after detrending and declustering have been performed.
Normal score transformation after detrending and declustering have been performed

Checking simulation output

Realizations should be checked to confirm the following:

Postprocessing

Once the realizations have been produced, they are usually postprocessed to obtain summaries of the results. The Gaussian Geostatistical Simulations tool allows several postprocessing options, which can be performed on the entire spatial extent of the rasters or on areas of particular interest. These areas are defined by specifying a polygon feature class in the Input statistical polygons option of the tool. Output is similar in both cases: postprocessing the entire rasters produces summary rasters, while postprocessing on polygonal areas produces a polygon output feature class that contains summary statistics for each polygon.

Postprocessing the entire raster extent

Postprocessing for areas of interest

Operators for post-processing areas of interest.
Operators for postprocessing areas of interest

Field name

Description

MIN

Minimum value of any cell in all the realizations that fall within the polygon.

MAX

Maximum value of any cell in all the realizations that fall within the polygon.

MEAN

Mean of all the cells in all the realizations that fall within the polygon.

STDDEV

Standard deviation of all the cells in all the realizations that fall within the polygon.

QUARTILE1

First quartile value of all the cells in all the realizations that fall within the polygon.

MEDIAN

Median value of all the cells in all the realizations that fall within the polygon.

QUARTILE3

Third quartile value of all the cells in all the realizations that fall within the polygon.

QUANTILE

Value corresponding to a user-specified quantile for all the cells in all the realizations that fall within the polygon.

P_THRSHLD

Percentage of cells that exceeded a user-specified threshold value, based on all the cells in all the realizations that fall within the polygon.

X_Y

The X function is applied to the values of the cells that fall within the polygon, one realization at a time. This process is equivalent to running the Zonal Statistics tool using the polygon as a zone and one realization at a time as the value grid. The Y function is applied to the values produced by the X function.

  • X can be one of the following functions: mean (MEAN), standard deviation (STD), first quartile (Q1), median (MED), third quartile (Q3), user-specified quantile (Q), or a user-specified threshold value to be exceeded (P).
  • Y can be one of the following functions: minimum (MIN), maximum (MAX), mean (MEAN), standard deviation (STD), first quartile (Q1), median (MED), or third quartile (Q3).

CELL_COUNT

The number of cells that fall within the polygon. If the cell center falls within the polygon, that cell is deemed to be in the polygon. A negative count indicates that part of the polygon falls outside the extent of the simulated raster and/or part of the polygon falls outside the clipping boundary. The negative number itself indicates the total number of cells that fall within the polygon.

SOURCE_ID

The object or feature ID of the input polygon feature class.

A description of the fields in the output polygon feature class

For both the bounding polygon and the polygon areas of interest options, raster cells are deemed to be inside the polygons if the center of the cell falls inside the polygon's boundary.

An example of conditional simulation and postprocessing output

The following figure shows the results of conditional simulation with polygon postprocessing of the output. The maps show the mean and standard deviation of 100 realizations of ozone values for each county in California. These mean and standard deviation values could be used, for example, in epidemiological studies where the occurrence of a disease needs to be compared with the average ozone value for each county.

Conditional simulation with polygon output.
Conditional simulation with polygon output

References and further reading

Aldworth, J. 1998. Spatial Prediction, Spatial Sampling, and Measurement Error. Ph.D. Thesis, Iowa State University. 186.

Chiles, J. P., and P. Delfiner. 1999. Geostatistics: Modeling Spatial Uncertainty. New York: John Wiley & Sons, 449–471.

Deutsch, C. V., and A. G. Journel. 1998. GSLIB Geostatistical Software Library and User's Guide, 2nd edition. New York: Oxford University Press, 119–141.

Dietrich, C. R., and G. N. Newsam. 1993. "A Fast and Exact Method for Multidimensional Gaussian Stochastic Simulations." Water Resources Research 29 (8): 2861–2869.

Goodchild, M. F., B. O. Parks, and L. T. Steyaert. 1993. Environmental Modeling with GIS. New York: Oxford University Press, 432–437.

Gribov, A., and K. Krivoruchko. 2004. "Geostatistical Mapping with Continuous Moving Neighborhood." Mathematical Geology 36 (2): 267–281.

Journel, A. G. 1974. "Geostatistics for Conditional Simulation of Ore Bodies." Economic Geology 69: 673–687.

Leuangthong, O., J. A. McLennan, and C. V. Deutsch. 2004. "Minimum Acceptance Criteria for Geostatistical Realizations." Natural Resources Research 13 (3): 131–141.

3/7/2014