On the influence of kriging parameters on the cartographic output—A study in mapping subglacial topography

Geostatistics provides a suite of methods, summarized as kriging, to analyze a finite data set to describe a continuous property of the Earth. Kriging methods consist of moving window optimum estimation techniques, which are based on a least-squares principle and use a spatial structure function, usually the variogram. Applications of kriging techniques have become increasingly wide-spread, with ordinary kriging and universal kriging being the most popular ones. The dependence of the final map or model on the input, however, is not generally understood. Herein we demonstrate how changes in the kriging parameters and the neighborhood search affect the cartographic result. Principles are illustrated through a glaciological study. The objective is to map ice thickness and subglacial topography of Storglaciären, Kebnekaise Massif, northern Sweden, from several sets of radio-echo soundings and hot water drillings. New maps are presented.     

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.     

Geostatistical Mapping with Continuous Moving Neighborhood

An issue that often arises in such GIS applications as digital elevation modeling (DEM) is how to create a continuous surface using a limited number of point observations. In hydrological applications, such as estimating drainage areas, direction of water flow is easier to detect from a smooth DEM than from a grid created using standard interpolation programs. Another reason for continuous mapping is esthetic; like a picture, a map should be visually appealing, and for some GIS users this is more important than map accuracy. There are many methods for local smoothing. Spline algorithms are usually used to create a continuous map, because they minimize curvature of the surface. Geostatistical models are commonly used approaches to spatial prediction and mapping in many scientific disciplines, but classical kriging models produce noncontinuous surfaces when local neighborhood is used. This motivated us to develop a continuous version of kriging. We propose a modification of kriging that produces continuous prediction and prediction standard error surfaces. The idea is to modify kriging systems so that data outside a specified distance from the prediction location have zero weights. We discuss simple kriging and conditional geostatistical simulation, models that essentially use information about mean value or trend surface. We also discuss how to modify ordinary and universal kriging models to produce continuous predictions, and limitations using the proposed models.     

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.     

Comparison of statistical prediction methods for characterizing the spatial variability of apparent electrical conductivity in coastal salt-affected farmland

Soil salinity has been known to be problematic to land productivity and environment in the lower Yellow River Delta due to the presence of a shallow, saline water table and marine sediments. Spatial information on soil salinity has gained increasing importance for the demand of management and sustainable utilization of arable land in this area. Apparent electrical conductivity, as measured by electromagnetic induction instrument in a fairly quick manner, has succeeded in mapping soil salinity and many other soil physical and chemical properties from field to regional scales. This was done based on the correlation that existed between apparent electrical conductivity and many other soil properties. In this paper, four spatial prediction methods, i.e., local polynomial, inverse distance weighed, ordinary kriging and universal kriging, were employed to estimate field-scale apparent electrical conductivity with the aid of an electromagnetic induction instrument (type EM38). The spatial patterns estimated by the four methods using EM38 survey datasets of various sample sizes were compared with those generated by each method using the entire sample size. Spatial similarity was evaluated using difference index (DI) between the maps created using various sample sizes (i.e., target maps) and the maps generated with the entire sample size (i.e., the reference map). The results indicated that universal kriging had the best performance owing to the inclusion of residuals and spatial detrending in the kriging system. DI showed that spatial similarity between the target and reference maps of apparent electrical conductivity decreased with the reduction in sample size for each prediction method. Under the same reduction in sample size, the method retaining the most spatial similarity was universal kriging, followed by ordinary kriging, inverse distance weighed, and local polynomial. Approximately, 70 % of total survey data essentially met the need for retaining 90 % details of the reference map for universal kriging and ordinary kriging methods. This conclusion was that OK and UK were two most appropriate methods for spatial estimation of apparent electrical conductivity as they were robust with the reduction in sample size.     

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.     

Quality mapping of the Ryytimaa dolomite in western Finland

The Ryytimaa dolomite formation in western Finland is mined for fertilizing purposes in agriculture. Ordinary kriging, indicator kriging, and indicator principal components simulation were used to map the dolomite quality. The fertilizing properties depend on the relative amounts of MgO and neutralizing CaO in dolomite. Based on the chemical composition the dolomite is divided into six quality classes according to the Finnish legislation. This quality classification is only partially in accordance with the three populations inferred from the distributions of the MgO, CaO, and neutralizing CaO contents. These populations probably represent compositionally different dolomite layers, now forming complicated fold interference structures that are visible on the kriged maps. The mapping of the quality, and thus, selective mining of the dolomite deposit is possible using the quality classification based on the three populations. In contrast, selective mining of the dolomite based on the legislation is difficult—at least with techniques used in this study.     

Advantages of using the kriging interpolator to estimate the gravity surface, comparison and spatial variability of gravity data in the El Kef-Ouargha region (northern Tunisia)

A comparative study on two interpolator methods for Bouguer anomaly mapping was performed in the El Kef-Ouargha region, Tunisia. In particular, we compare the results of minimum curvature gridding and kriging. We undertake a geostatistical study using both anisotropic and isotropic analysis for studding the spatial variability and kriging the random variable: “Bouguer anomaly”. The model parameters of each method are calculated to compare the results. The anisotropic analysis most closely approximates the measured data. Consequentially, the Bouguer anomaly and vertical gradient maps using anisotropic kriging have the best correlation with the geological map, best represent density distributions at depth and most closely approximate the measured anomaly amplitudes.     

Accounting for exhaustive secondary data into the mapping of water table elevation

This paper presents the incorporation of a digital elevation model into the spatial prediction of water table elevation in Mazandaran province (Iran) using a range of interpolation techniques. The multivariate methods used are: linear regression (LR), cokriging (COK), kriging with an external drift (KED) and regression kriging (RK). The analysis is performed on 3 years (1987, 1997 and 2007) of water table elevation data from about 260 monitoring wells. Prediction performances of the different algorithms are compared with two univariate techniques, i.e. inverse distance weighting and ordinary kriging (OK), through cross validation and examination of the consistency of the generated maps with the natural phenomena. Significantly smaller prediction errors are obtained for four multivariate algorithms but, in particular, KED and RK outperform LR and COK for 3 years. The results show the potential for using elevation for a more precise mapping of water table elevation.     

Geostatistical approach for identification of transmissivity structure at Dulliu area in Taiwan

 A thorough understanding of the characteristics of transmissivity makes groundwater deterministic models more accurate. These transmissivity data characteristics occasionally possess a complicated spatial variation over an investigated site. This study presents both geostatistical estimation and conditional simulation methods to generate spatial transmissivity maps. The measured transmissivity data from the Dulliu area in Yun-Lin county, Taiwan, is used as the case study. The spatial transmissivity maps are simulated by using sequential Gaussian simulation (SGS), and estimated by using natural log ordinary kriging and ordinary kriging. Estimation and simulation results indicate that SGS can reproduce the spatial structure of the investigated data. Furthermore, displaying a low spatial variability does not allow the ordinary kriging and natural log kriging estimates to fit the spatial structure and small-scale variation for the investigated data. The maps of kriging estimates are smoother than those of other simulations. A SGS with multiple realizations has significant advantages over ordinary kriging and even natural log kriging techniques at a site with a high variation in investigated data. These results are displayed in geographic information systems (GIS) as basic information for further groundwater study. Received: 27 August 1999 · Accepted: 22 February 2000