Kriging with imprecise (fuzzy) variograms. II: Application

The geostatistical analysis of soil liner permeability is based on 20 measurements and imprecise prior information on nugget effect, sill, and range of the unknown variogram. Using this information, membership functions for variogram parameters are assessed and the fuzzy variogram is constructed. Both kriging estimates and estimation variances are calculated as fuzzy numbers from the fuzzy variogram and data points. Contour maps are presented, indicating values of the kriged permeability and the estimation variance corresponding to selected membership values called levels.     

Variogram models for regional groundwater geochemical data

Four variogram models for regional groundwater geochemical data are presented. These models were developed from an empirical study of the sample variograms for more than 10 elements in groundwaters from two geologic regions in the Plainview quandrangle, Texas. A procedure is given for the estimation of the variogram in the isotropic and anisotropic case. The variograms were found useful for quantifying the differences in spatial variability for elements within a geologic unit and for elements in different geologic units. Additionally, the variogram analysis enables assessment of the assumption of statistical independence of regional samples which is commonly used in many statistical procedures. The estimated variograms are used in computation of kriged estimates for the Plainview quadrangle data. The results indicate that an inverse distance weighting model was superior for prediction than simple kriging with the particular variograms used.     

Kriging of water levels in the Souss aquifer,Morocco

Universal kriging is applied to water table data from the Souss aquifer in central Morocco. The procedure accounts for the spatial variability of the phenomenon to be mapped. With the use of measured elevations of the water table, an experimental variogram is constructed that characterizes the spatial variability of the measured water levels. Spherical and Gaussian variogram models are alternatively used to fit the experimental variogram. The models are used to develop contour maps of water table elevations and corresponding estimation variances. The estimation variances express the reliability of the kriged water table elevation maps. Universal kriging also provides a contour map of the expected elevation of the water table (drift). The differences between the expected and measured water table elevations are called residuals from the drift. Residuals from the drift are compared with residuals obtained by more traditional least-squares analysis.     

On Visualization for Assessing Kriging Outcomes

Extant opinion about kriging is that all weights should be positive. Visualizations rendered by converting kriged grids to digital images are presented to show that negative weights may be beneficial to some spatial problems. In particular, variogram models with zero-valued nuggets, already well known to minimize smoothing through kriging, result in a visual resolution substantially superior to that from kriging with a variogram model having a nonzero nugget value in application to satellite acquired data. Negative weights are more likely when using variogram models with zero-valued nuggets, but resultant visualizations often show a smoother transition between extreme data values. This is true even when a variogram model having a nugget value of zero is not optimum with respect to mean square error, as is demonstrated using a nitrate data set. An analogy to digital image processing is used to suggest that the influence of negative weights in kriging is similar to a high-boost kernel.     

Kriging with imprecise (fuzzy) variograms. I: Theory

Imprecise variogram parameters are modeled with fuzzy set theory. The fit of a variogram model to experimental variograms is often subjective. The accuracy of the fit is modeled with imprecise variogram parameters. Measurement data often are insufficient to create good experimental variograms. In this case, prior knowledge and experience can contribute to determination of the variogram model parameters. A methodology for kriging with imprecise variogram parameters is developed. Both kriged values and estimation variances are calculated as fuzzy numbers and characterized by their membership functions. Besides estimation variance, the membership functions are used to create another uncertainty measure. This measure depends on both homogeneity and configuration of the data.     

FuzzyKrig: a comprehensive matlab toolbox for geostatistical estimation of imprecise information

Using kriging has been accepted today as the most common method of estimating spatial data in such different fields as the geosciences. To be able to apply kriging methods, it is necessary that the data and variogram model parameters be precise. To utilize the imprecise (fuzzy) data and parameters, use is made of fuzzy kriging methods. Although it has been 30 years since different fuzzy kriging algorithms were proposed, its use has not become as common as other kriging methods (ordinary, simple, log, universal, etc.); lack of a comprehensive software that can perform, based on different fuzzy kriging algorithms, the related calculations in a 3D space can be the main reason. This paper describes an open-source software toolbox (developed in Matlab) for running different algorithms proposed for fuzzy kriging. It also presents, besides a short presentation of the fuzzy kriging method and introduction of the functions provided by the FuzzyKrig toolbox, 3 cases of the software application under the conditions where: 1) data are hard and variogram model parameters are fuzzy, 2) data are fuzzy and variogram model parameters are hard, and 3) both data and variogram model parameters are fuzzy.     

Orthonormal residuals in geostatistics: Model criticism and parameter estimation

In linear geostatistics, models for the mean function (drift) and the variogram or generalized covariance function are selected on the basis of the modeler's understanding of the phenomenon studied as well as data. One can seldom be assured that the most appropriate model has been selected; however, analysis of residuals is helpful in diagnosing whether some important characteristic of the data has been neglected and, ultimately, in providing a reasonable degree of assurance that the selected model is consistent with the available information. The orthonormal residuals presented in this work are kriging errors constructed so that, when the correct model is used, they are uncorrelated and have zero mean and unit variance. It is suggested that testing of orthonormal residuals is a practical way for evaluating the agreement of the model with the data and for diagnosing model deficiencies. Their advantages over the usually employed standardized residuals are discussed. A set of tests are presented. Orthonormal residuals can also be useful in the estimation of the covariance (or variogram) parameters for a model that is considered correct.     

New variogram modeling method using MGGP and SVR

A critical step for kriging in geostatistics is estimation of the variogram. Traditional variogram modeling comprise of the experimental variogram calculation, appropriate variogram model selection and model parameter determination. Selecting of the variogram model and fitting of model parameters is the most controversial aspect of geostatistics. Shapes of valid variogram models are finite, and sometimes, the optimal shape of the model can not be fitted, leading to reduced estimation accuracy. In this paper, a new method is presented to automatically construct a model shape and fit model parameters to experimental variograms using Support Vector Regression (SVR) and Multi-Gene Genetic Programming (MGGP). The proposed method does not require the selection of a variogram model and can directly provide the model shape and parameters of the optimal variogram. The validity of the proposed method is demonstrated in a number of cases.     

Robustness of variograms and conditioning of kriging matrices

Current ideas of robustness in geostatistics concentrate upon estimation of the experimental variogram. However, predictive algorithms can be very sensitive to small perturbations in data or in the variogram model as well. To quantify this notion of robustness, nearness of variogram models is defined. Closeness of two variogram models is reflected in the sensitivity of their corresponding kriging estimators. The condition number of kriging matrices is shown to play a central role. Various examples are given. The ideas are used to analyze more complex universal kriging systems.Research performed while on leave at Centre de Geóstatistique et de Morphologie Mathématique, Fontainebleau.     

Robust kriging—A proposal

Geological data frequently have a heavy-tailed normal-in-the-middle distribution, which gives rise to grade distributions that appear to be normal except for the occurrence of a few outliers. This same situation also applies to log-transformed data to which lognormal kriging is to be applied. For such data, linear kriging is nonrobust in that (1)kriged estimates tend to infinity as the outliers do, and (2)it is also not minimum mean squared error. The more general nonlinear method of disjunctive kriging is even more nonrobust, computationally more laborious, and in the end need not produce better practical answers. We propose a robust kriging method for such nearly normal data based on linear kriging of an editing of the data. It is little more laborious than conventional linear kriging and, used in conjunction with a robust estimator of the variogram, provides good protection against the effects of data outliers. The method is also applicable to time series analysis.     

Evaluating Thickness of Bauxite Deposit Using Indicator Geostatistics and Fuzzy Estimation

Evaluating the geological properties of a mineral deposit is a fundamental task for mine planning and it requires an assessment of reserve parameters such as thickness and grade. This paper presents a linguistic model for estimating bauxite thickness within a deposit in a fuzzy environment using indicator geostatistics and fuzzy modeling. The proposed model consists of two main stages: determining the orebody boundary and estimating the thickness. In order to estimate the thickness, a rule‐based fuzzy inference mechanism has been developed based on data statistics. Results and performance of the model have been compared with that of a well‐known conventional technique, geostatistics (kriging), and it is shown that the proposed model has bigger estimation power. In addition, the fuzzy approach is more flexible than the kriging approach. The fuzzy methodology used in the present paper is convenient for modeling reserve parameters.     

Bayesian Maximum Entropy Analysis and Mapping: A Farewell to Kriging Estimators?

The Bayesian Maximum Entropy (BME) method of spatial analysis and mapping provides definite rules for incorporating prior information, hard and soft data into the mapping process. It has certain unique features that make it a loyal guardian of plausible reasoning under conditions of uncertainty. BME is a general approach that does not make any assumptions regarding the linearity of the estimator, the normality of the underlying probability laws, or the homogeneity of the spatial distribution. By capitalizing on various sources of information and data, BME introduces an epistemological framework that produces predictive maps that are more accurate and in many cases computationally more efficient than those derived by traditional techniques. In fact, kriging techniques can be derived as special cases of the BME approach, under restrictive assumptions regarding the prior information and the data available. BME is a more rigorous approach than indicator kriging for incorporating soft data. The BME formulation, in fact, applies in a spatial or a spatiotemporal domain and its extension to the case of block and vector random fields is straightforward. New theoretical results are presented and numerical examples are discussed, which use the BME approach to account for important sources of knowledge in a systematic manner. BME can be useful in practical situations in which prior information can be used to compensate for the limited amount of measurements available (e.g., preliminary or feasibility study levels) or soft data are available that can be combined with hard data to improve mapping significantly. BME may be then viewed as an effort towards the development of a more general framework of spatial/temporal analysis and mapping, which includes traditional geostatistics as its limiting case, and it also provides the means to derive novel results that could not be obtained by traditional geostatistics.     

Indicator-Based Geostatistical Models For Mapping Fish Survey Data

Marine research survey data on fish stocks often show a small proportion of very high-density values, as for many environmental data. This makes the estimation of second-order statistics, such as the variance and the variogram, non-robust. The high fish density values are generated by fish aggregative behaviour, which may vary greatly at small scale in time and space. The high values are thus imprecisely known, both in their spatial occurrence and order of magnitude. To map such data, three indicator-based geostatistical methods were considered, the top-cut model, min–max autocorrelation factors (MAF) of indicators, and multiple indicator kriging. In the top-cut and MAF approaches, the variable is decomposed into components and the most continuous ones (those corresponding to the low and medium values) are used to guide the mapping. The methods are proposed as alternatives to ordinary kriging when the variogram is difficult to estimate. The methods are detailed and applied on a spatial data set of anchovy densities derived from a typical fish stock acoustic survey performed in the Bay of Biscay, which show a few high-density values distributed in small spatial patches and also as solitary events. The model performances are analyzed by cross-validating the data and comparing the kriged maps. Results are compared to ordinary kriging as a base case. The top-cut model had the best cross-validation performance. The indicator-based models allowed mapping high-value areas with small spatial extent, in contrast to ordinary kriging. Practical guidelines for implementing the indicator-based methods are provided.     

Kriging in the Presence of Locally Varying Anisotropy Using Non-Euclidean Distances

A stationary specification of anisotropy does not always capture the complexities of a geologic site. In this situation, the anisotropy can be varied locally. Directions of continuity and the range of the variogram can change depending on location within the domain being modeled. Kriging equations have been developed to use a local anisotropy specification within kriging neighborhoods; however, this approach does not account for variation in anisotropy within the kriging neighborhood. This paper presents an algorithm to determine the optimum path between points that results in the highest covariance in the presence of locally varying anisotropy. Using optimum paths increases covariance, results in lower estimation variance and leads to results that reflect important curvilinear structures. Although CPU intensive, the complex curvilinear structures of the kriged maps are important for process evaluation. Examples highlight the ability of this methodology to reproduce complex features that could not be generated with traditional kriging.     

Mapping Soil Particle Size Fractions Using Compositional Kriging,Cokriging and Additive Log-ratio Cokriging in Two Case Studies

Information on the spatial distribution of soil particle-size fractions (psf) is required for a wide range of applications. Geostatistics is often used to map spatial distribution from point observations; however, for compositional data such as soil psf, conventional multivariate geostatistics are not optimal. Several solutions have been proposed, including compositional kriging and transformation to a composition followed by cokriging. These have been shown to perform differently in different situations, so that there is no procedure to choose an optimal method. To address this, two case studies of soil psf mapping were carried out using compositional kriging, log-ratio cokriging, cokriging, and additive log-ratio cokriging; and the performance of Mahalanobis distance as a criterion for choosing an optimal mapping method was tested. All methods generated very similar results. However, the compositional kriging and cokriging results were slightly more similar to each other than to the other pair, as were log-ratio cokriging and additive log-ratio cokriging. The similar results of the two methods within each pair were due to similarities of the methods themselves, for example, the same variogram models and prediction techniques, and the similar results between the two pairs were due to the mathematical relationship between original and log-ratio transformed data. Mahalanobis distance did not prove to be a good indicator for selecting an optimal method to map soil psf.     

Disjunctive Kriging with Hard and Imprecise Data

This paper presents a methodology for assessing local probability distributions by disjunctive kriging when the available data set contains some imprecise measurements, like noisy or soft information or interval constraints. The basic idea consists in replacing the set of imprecise data by a set of pseudohard data simulated from their posterior distribution; an iterative algorithm based on the Gibbs sampler is proposed to achieve such a simulation step. The whole procedure is repeated many times and the final result is the average of the disjunctive kriging estimates computed from each simulated data set. Being data-independent, the kriging weights need to be calculated only once, which enables fast computing. The simulation procedure requires encoding each datum as a pre-posterior distribution and assuming a Markov property to allow the updating of pre-posterior distributions into posterior ones. Although it suffers some imperfections, disjunctive kriging turns out to be a much more flexible approach than conditional expectation, because of the vast class of models that allows its computation, namely isofactorial models.     

Comparison of kriging techniques in a space-time context

Space-time processes constitute a particular class, requiring suitable tools in order to predict values in time and space, such as a space-time variogram or covariance function. The space-time co-variance function is defined and linked to the Linear Model of Coregionalization under second-order space-time stationarity. Simple and ordinary space-time kriging systems are compared to simple and ordinary cokriging and their differences for unbiasedness conditions are underlined. The ordinary space-time kriging estimation then is applied to simulated data. Prediction variances and prediction errors are compared with those for ordinary kriging and cokriging under different unbiasedness conditions using a cross-validation. The results show that space-time kriging tend to produce lower prediction variances and prediction errors that kriging and cokriging.     

Variogram Model Selection via Nonparametric Derivative Estimation

Before optimal linear prediction can be performed on spatial data sets, the variogram is usually estimated at various lags and a parametric model is fitted to those estimates. Apart from possible a priori knowledge about the process and the user's subjectivity, there is no standard methodology for choosing among valid variogram models like the spherical or the exponential ones. This paper discusses the nonparametric estimation of the variogram and its derivative, based on the spectral representation of positive definite functions. The use of the estimated derivative to help choose among valid parametric variogram models is presented. Once a model is selected, its parameters can be estimated—for example, by generalized least squares. A small simulation study is performed that demonstrates the usefulness of estimating the derivative to help model selection and illustrates the issue of aliasing. MATLAB software for nonparametric variogram derivative estimation is available at http://www-math.mit.edu/~gorsich/derivative.html. An application to the Walker Lake data set is also presented.     

Application of geostatistics to identify gold-rich areas in the Finisterre-Fervenza region,NW Spain

Three univariate geostatistical methods of estimation are applied to a geochemical data set. The studied methods are: ordinary kriging (cross-validation), factorial kriging, and indicator kriging. These techniques use the probabilistic and spatial behaviour of geochemical variables, giving a tool for identifying potential anomalous areas to locate mineralization. Ordinary kriging is easy to apply and to interpret the results. It has the advantage of using the same experimental grid points for its estimates, and no additional grid points are needed. Factorial kriging decomposes the raw variable into as many components as there are identified structures in the variogram. This, however, is a complex method and its application is more difficult than that of ordinary or indicator kriging. The main advantages of indicator kriging are that data are used by their rank order, being more robust about outlier values, and that the presentation of results is simple. Nevertheless, indicator kriging is incapable of separating anomalous values and the high values from the background, which have a behaviour different to the anomaly. In this work, the results of the application of these 3 kriging methods to a set of mineral exploration data obtained from a geochemical survey carried out in NW Spain are presented. This area is characterised by the presence of Au mineral occurrences. The kriging methods were applied to As, considered as a pathfinder of Au in this area. Numerical treatment of Au is not applicable, because it presents most values equal to the detection limit, and a series of extreme values. The results of the application of ordinary kriging, factorial kriging and indicator kriging to As make possible the location of a series of rich values, sited along a N–S shear zone, considered a structure related to the presence of Au.