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.     

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.     

Comparison of two techniques for applying disjunctive kriging: The Gaussian anamorphosis model versus the direct statistical inference of the bivariate distributions

Indicator cokriging is an alternative to disjunctive kriging for estimation of spatial distributions. One way to determine which of these techniques is more accurate for estimation of spatial distributions is to apply each to a particular type of data. A procedure is developed for evaluation of disjunctive kriging and indicator cokriging for such an application. Application of this procedure to earthquake ground motion data found disjunctive kriging to be at least as accurate as indicator cokriging for estimation of spatial distributions of peak horizontal acceleration. Indicator cokriging was superior for all other types of earthquake ground motion data.     

Transfer function-noise modeling and spatial interpolation to evaluate the risk of extreme (shallow) water-table levels in the Brazilian Cerrados

Water regimes in the Brazilian Cerrados are sensitive to climatological disturbances and human intervention. The risk that critical water-table levels are exceeded over long periods of time can be estimated by applying stochastic methods in modeling the dynamic relationship between water levels and driving forces such as precipitation and evapotranspiration. In this study, a transfer function-noise model, the so called PIRFICT-model, is applied to estimate the dynamic relationship between water-table depth and precipitation surplus/deficit in a watershed with a groundwater monitoring scheme in the Brazilian Cerrados. Critical limits were defined for a period in the Cerrados agricultural calendar, the end of the rainy season, when extremely shallow levels (<0.5-m depth) can pose a risk to plant health and machinery before harvesting. By simulating time-series models, the risk of exceeding critical thresholds during a continuous period of time (e.g. 10 days) is described by probability levels. These simulated probabilities were interpolated spatially using universal kriging, incorporating information related to the drainage basin from a digital elevation model. The resulting map reduced model uncertainty. Three areas were defined as presenting potential risk at the end of the rainy season. These areas deserve attention with respect to water-management and land-use planning.     

On the Equivalence of the Cokriging and Kriging Systems

Simple cokriging of components of a p-dimensional second-order stationary random process is considered. Necessary and sufficient conditions under which simple cokriging is equivalent to simple kriging are given. Essentially this condition requires that it should be possible to express the cross-covariance at any lag series h using the cross-covariance at |h|=0 and the auto-covariance at lag series h. The mosaic model, multicolocated kriging and the linear model of coregionalization are examined in this context. A data analytic method to examine whether simple kriging of components of a multivariate random process is equivalent to its cokriging is given     

Comparison of approaches to spatial estimation in a bivariate context

The problem of estimating a regionalized variable in the presence of other secondary variables is encountered in spatial investigations. Given a context in which the secondary variable is known everywhere (or can be estimated with great precision), different estimation methods are compared: regression, regression with residual simple kriging, kriging, simple kriging with a mean obtained by regression, kriging with an external drift, and cokriging. The study focuses on 19 pairs of regionalized variables from five different datasets representing different domains (geochemical, environmental, geotechnical). The methods are compared by cross-validation using the mean absolute error as criterion. For correlations between the principal and secondary variable under 0.4, similar results are obtained using kriging and cokriging, and these methods are superior slightly to the other approaches in terms of minimizing estimation error. For correlations greater than 0.4, cokriging generally performs better than other methods, with a reduction in mean absolute errors that can reach 46% when there is a high degree of correlation between the variables. Kriging with an external drift or kriging the residuals of a regression (SKR) are almost as precise as cokriging.     

Criteria to Compare Estimation Methods of Regionalized Compositions

The additive logratio (alr) transformation has been used in several case studies to predict regionalized compositions using standard geostatistical estimation methods such as ordinary kriging and ordinary cokriging. It is a simple method that allows application to transformed data all the body of knowledge available for geostatistical analysis of coregionalizations without a constant sum constraint. To compare the performance of methods, it is customary to use a univariate crossvalidation approach based on the leaving-one-out technique to evaluate the performance for each attribute separately. For multivariate observations this approach is difficult to interpret in terms of overall performance. Therefore, we propose using appropriate distances in real space and in the simplex, to improve the crossvalidation approach and, going a step forward, to adapt the concept of stress from multidimensional scaling to obtain a global measure of performance for each method. The Lyons West oil field of Kansas is used to illustrate the impactof using different distances in the performance of ordinary kriging versus ordinary cokriging.     

CAD-DRASTIC: chloride application density combined with DRASTIC for assessing groundwater vulnerability to road salt application

Road salt is pervasively used throughout Canada and in other cold regions during winter. For cities relying exclusively on groundwater, it is important to plan and minimize the application of salt accordingly to mitigate the adverse effects of high chloride concentrations in water supply aquifers. The use of geospatial data (road network, land use, Quaternary and bedrock geology, average annual recharge, water-table depth, soil distribution, topography) in the DRASTIC methodology provides an efficient way of distinguishing salt-vulnerable areas associated with groundwater supply wells, to aid in the implementation of appropriate management practices for road salt application in urban areas. This research presents a GIS-based methodology to accomplish a vulnerability analysis for 12 municipal water supply wells within the City of Guelph, Ontario, Canada. The chloride application density (CAD) value at each supply well is calculated and related to the measured groundwater chloride concentrations and further combined with soil media and aquifer vadose- and saturated-zone properties used in DRASTIC. This combined approach, CAD-DRASTIC, is more accurate than existing groundwater vulnerability mapping methods and can be used by municipalities and other water managers to further improve groundwater protection related to road salt application.     

Comparative performance of indicator algorithms for modeling conditional probability distribution functions

This paper compares the performance of four algorithms (full indicator cokriging. adjacent cutoffs indicator cokriging, multiple indicator kriging, median indicator kriging) for modeling conditional cumulative distribution functions (ccdf).The latter three algorithms are approximations to the theoretically better full indicator cokriging in the sense that they disregard cross-covariances between some indicator variables or they consider that all covariances are proportional to the same function. Comparative performance is assessed using a reference soil data set that includes 2649 locations at which both topsoil copper and cobalt were measured. For all practical purposes, indicator cokriging does not perform better than the other simpler algorithms which involve less variogram modeling effort and smaller computational cost. Furthermore, the number of order relation deviations is found to be higher for cokriging algorithms, especially when constraints on the kriging weights are applied.     

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.     

Extending Multivariate Space-Time Geostatistics for Environmental Data Analysis

Environmental studies require multivariate data such as chemical concentrations with space-time coordinates. There are two general conditions related to such data: the existence of correlations among the coregionalized variables and the differences in numbers of data which occur because of insufficient data caused by measurement error or bad weather conditions. This study proposes geostatistical techniques for space-time multivariate modeling that take into consideration these correlations and data absences. These techniques consist of suitable modeling of semivariograms and cross-semivariograms for quantifying correlation structures among multivariables and of extending standardized ordinary cokriging. The tensor product cubic smoothing surface method is used for space-time semivariogram modeling. These methods are applied to the chemical component data of the Ariake Sea, a typical closed sea in southwest Japan. In order to clarify environmental changes in the Ariake Sea, the concentration data of four nutritive salts (NO₂–N, NO₃–N, NH₄–N, and PO₄–P) at 38 stations over 25 years are used as environmental indicators. For each of the kinds of data, there are spaces and times for which there is no data available. The effectiveness of the modeling of space-time semivariograms and the high estimation capability of the extended cokriging are demonstrated by cross-validation. Compared with ordinary kriging for a single variable, multivariate space-time standardized ordinary cokriging can provide a more detailed concentration map of nutritive salts and while elucidating their temporal changes over sparsely spaced data areas. In the space-time models by ordinary kriging, on the other hand, smooth trends are obvious.     

Estimating Freshwater Acidification Critical Load Exceedance Data for Great Britain Using Space-Varying Relationship Models

In this study, two distinct sets of analyses are conducted on a freshwater acidification critical load dataset, with the objective of assessing the quality of various models in estimating critical load exceedance data. Relationships between contextual catchment and critical load data are known to vary across space; as such, we cater for this in our model choice. Firstly, ordinary kriging (OK), multiple linear regression (MLR), geographically weighted regression (GWR), simple kriging with GWR-derived local means (SKlm-GWR), and kriging with an external drift (KED) are used to predict critical loads (and exceedances). Here, models that cater for space-varying relationships (GWR; SKlm-GWR; KED using local neighbourhoods) make more accurate predictions than those that do not (MLR; KED using a global neighbourhood), as well as in comparison to OK. Secondly, as the chosen predictors are not suited to providing useable estimates of critical load exceedance risk, they are replaced with indicator kriging (IK) models. Here, an IK model that is newly adapted to cater for space-varying relationships performs better than those that are not adapted in this way. However, when site misclassification rates are found using either exceedance predictions or estimates of exceedance risk, rates are intolerably high, reflecting much underlying noise in the data.     

Accuracy assessment of noise mapping on the main street

Among the sounds constituting a part of our daily lives, the undesired and disturbing ones are named noise. The noises can be defined as the acoustic energy that affects the physiological and psychological welfare of humans due to its direct or indirect effects on the human health. But, the traffic noise, which is one of the disturbing noises having wide dispersion in urban areas, affects the health of a gradually increasing number of people. The number of studies on measurement, dispersion, and mapping of traffic noise in developed and developing countries is gradually increasing. It can be seen that the studies concentrate on mapping based on the city centers and neighborhoods, where the traffic is intense. In this article, unlike the regional noise measurement studies, the transversal cross sections were established with a certain interval according to the land status along the main street, and the noise measurement points were determined using the global navigation satellite system (GNSS). The daily mean noise levels (L_den) were calculated by performing the measurements at daytime, evening, and nighttime. Using ArcGIS 10.4 software, the accuracy of noise maps created via inverse distance weighting (IDW), radial basis function (RBF), and ordinary kriging (OK) geostatistical analysis interpolation methods and the dispersion effect of traffic noise were compared. According to the visual method and error corrections, it was found that the best noise map was obtained using the RBF method. It is believed that the map provides accurate and reliable results because of the flat land surface, higher noise values measured on the road region, decrease in the noise at the points farther from the road region, and soft transitions between the noise values depending on the distance from the road. In order to take measures in order to decrease or eliminate noise on the main street, the local administrators and decision-makers can make use of these maps, the accuracy of which has been confirmed.     

Optimizing a piezometric network in the estimation of the groundwater budget: a case study from a crystalline-rock watershed in southern India

An estimate of the groundwater budget at the catchment scale is extremely important for the sustainable management of available water resources. Water resources are generally subjected to over-exploitation for agricultural and domestic purposes in agrarian economies like India. The double water-table fluctuation method is a reliable method for calculating the water budget in semi-arid crystalline rock areas. Extensive measurements of water levels from a dense network before and after the monsoon rainfall were made in a 53 km² watershed in southern India and various components of the water balance were then calculated. Later, water level data underwent geostatistical analyses to determine the priority and/or redundancy of each measurement point using a cross-validation method. An optimal network evolved from these analyses. The network was then used in re-calculation of the water-balance components. It was established that such an optimized network provides far fewer measurement points without considerably changing the conclusions regarding groundwater budget. This exercise is helpful in reducing the time and expenditure involved in exhaustive piezometric surveys and also in determining the water budget for large watersheds (watersheds greater than 50 km²).     

Estimating Regional Hydraulic Conductivity Fields—A Comparative Study of Geostatistical Methods

Geostatistical estimations of the hydraulic conductivity field (K) in the Carrizo aquifer, Texas, are performed over three regional domains of increasing extent: 1) the domain corresponding to a three-dimensional groundwater flow model previously built (model domain); 2) the area corresponding to the 10 counties encompassing the model domain (County domain), and; 3) the full extension of the Carrizo aquifer within Texas (Texas domain). Two different approaches are used: 1) an indirect approach where transmissivity (T) is estimated first and K is retrieved through division of the T estimate by the screen length of the wells, and; 2) a direct approach where K data are kriged directly. Due to preferential well screen emplacement, and scarcity of sampling in the deeper portions of the formation (> 1 km), the available data set is biased toward high values of hydraulic conductivities. Kriging combined with linear regression, simple kriging with varying local means, kriging with an external drift, and cokriging allow the incorporation of specific capacity as secondary information. Prediction performances (assessed through cross-validation) differ according to the chosen approach, the considered variable (log-transformed or back-transformed), and the scale of interest. For the indirect approach, kriging of log T with varying local means yields the best estimates for both log-transformed and back-transformed variables in the model domain. For larger regional scales (County and Texas domains), cokriging performs generally better than other kriging procedures when estimating both (log T)^∗ and T^∗. Among procedures using the direct approach, the best prediction performances are obtained using kriging of log K with an external drift. Overall, geostatistical estimation of the hydraulic conductivity field at regional scales is rendered difficult by both preferential well location and preferential emplacement of well screens in the most productive portions of the aquifer. Such bias creates unrealistic hydraulic conductivity values, in particular, in sparsely sampled areas.     

Multicriterial analysis in the investigation of favorable areas for edifications with shallow and deep foundations in the Municipality of Campos dos Goytacazes — Rio de Janeiro, Brazil

The Municipality of Campos dos Goytacazes is located on the northern coast of Rio de Janeiro, Brazil. The city developed on a floodplain and as a consequence of this occupation it has experienced problems with the water table level, drainage system, soil carrying capacity and aquifers' contamination. The purpose of this paper was to accomplish a geological and environmental study, producing interpretive and derived geotechnical maps to analyze the degree of aptitudes or restrictions of the physical environment for urban expansion. In order to do this, the methodology of multicriterial analysis was used, combined with pedological indicators resulting in the elaboration of the Favorability Maps for Buildings (Shallow and Deep Foundations). The remote sensing tools, GIS and GPS, sped up the elaboration of the Information Plans. The maps generated in this work allow us to indicate problems in the physical environment serving as a guide to urbanization. It was verified that the city expands into unsuitable areas and it is recommended as areas for expansion those considered favorable in the outcomes of this research.     

Data Configurations and the Cokriging System: Simplification by Screen Effects

Large cokriging systems arise in many situations and are difficult to handle in practice. Simplifications such as simple kriging, strictly collocated and multicollocated cokriging are often used and models under which such simplifications are, in fact, equivalent to cokriging have recently received attention. In this paper, a two-dimensional second-order stationary random process with known mean is considered and the redundancy of certain components of the data at certain locations vis-à-vis the solution to the simple cokriging system is examined. Conditions for the simple cokriging weights of these components at these locations are set to zero. The conditions generalise the notion of the autokrigeability coefficient and can, in principle, be applied to any data configuration. In specific sampling situations such as the isotopic and certain heterotropic configurations, models under which simple kriging, strictly collocated, multicollocated and dislocated cokriging are equivalent to simple cokriging are readily identified and results already available in the literature are obtained. These are readily identified and the results are already available in the literature. The advantage of the approach presented here is that it can be applied to any data configuration for analysis of permissible simplifications in simple cokriging.     

Restricted kriging for mixture of grade models

A modified type of kriging, referred to as restricted kriging (RK), is proposed in this study. The method incorporates constraints on different grade classes to restrict the influence of the samples having different likelihoods in estimation. RK is motivated by the estimation of mineral reserves when grades have highly skewed distributions. Ordinary kriging tends to produce an overly smoothed interpolated surface by underestimating high grades and overestimating low grades. The fact that ordinary kriging gives a uniform prior treatment to all samples independent of their values is a major factor associated with this smoothing effect. The new approach differentiates each grade portion by preselected cutoffs. RK is developed for a single cutoff and then extended into a general form for any finite number of cutoffs. Restricted cokriging (RCK) is also formulated to simultaneously estimate a set of random functions with restriction conditions. Methods are suggested for determination of the probabilities of occurrence of different grade portions. Finally, the new approach is demonstrated on a case study of an epithermal gold deposit.     

Preliminary flood risk assessment: the case of Athens

Flood mapping, especially in urban areas, is a demanding task requiring substantial (and usually unavailable) data. However, with the recent introduction of the EU Floods Directive (2007/60/EC), the need for reliable, but cost effective, risk mapping at the regional scale is rising in the policy agenda. Methods are therefore required to allow for efficiently undertaking what the Directive terms “preliminary flood risk assessment,” in other words a screening of areas that could potentially be at risk of flooding and that consequently merit more detailed attention and analysis. Such methods cannot rely on modeling, as this would require more data and effort that is reasonable for this high-level, screening phase. This is especially true in urban areas, where modeling requires knowledge of the detailed urban terrain, the drainage networks, and their interactions. A GIS-based multicriteria flood risk assessment methodology was therefore developed and applied for the mapping of flood risk in urban areas. This approach quantifies the spatial distribution of flood risk and is able to deal with uncertainties in criteria values and to examine their influence on the overall flood risk assessment. It can further assess the spatially variable reliability of the resulting maps on the basis of the choice of method used to develop the maps. The approach is applied to the Greater Athens area and validated for its central and most urban part. A GIS database of economic, social, and environmental criteria contributing to flood risk was created. Three different multicriteria decision rules (Analytical Hierarchy Process, Weighted Linear Combination and Ordered Weighting Averaging) were applied, to produce the overall flood risk map of the area. To implement this methodology, the IDRISI Andes GIS software was customized and used. It is concluded that the results of the analysis are a reasonable representation of actual flood risk, on the basis of their comparison with historical flood events.     

Geotechnical characterization of the upper Pleistocene–Holocene alluvial deposits of Roma (Italy) by means of multivariate geostatistics: Cross-validation results

We are presenting an attempt to evaluate the spatial variability of geotechnical parameters in the upper Pleistocene–Holocene alluvial deposits of Roma (Italy) by means of multivariate geostatistics.The upper Pleistocene–Holocene alluvial deposits of Roma are sensitive to high levels of geohazard. They occupy a sizable and significant part of the city, being the foundation for many monuments, historical neighborhoods, and archaeological areas, and the main host of the present and future subway lines. We have stored information from more than 2000 geotechnical boreholes crossing the alluvial deposits into a relational database. For the present study, only the boreholes with lithologic/textural interpretation and geotechnical information were selected. The set includes 283 boreholes and 719 samples, which have a set of geotechnical information comprising physical properties and mechanical parameters.Techniques of multivariate statistics and geostatistics were combined and compared to evaluate the estimation methods of the mechanical parameters, with special reference to the drained friction angle from direct shear test (φ′). Principal Component Analysis was applied to the dataset to highlight the relationships between the geotechnical parameters. Through cross-validation analysis, multiple linear regression, kriging, and cokriging were tested as estimators of φ′. Cross-validation demonstrates that the cokriging with granulometries as auxiliary variables is the most suitable method to estimate φ′. In addition to proving that cokriging is a good estimator of φ′, cross-validation demonstrates that input data are coherent and this allows us to use them for estimation of geotechnical parameters, although they come from different laboratories and different vintages.Nevertheless, to get the same good results of cross-validation in estimation, it is necessary for granulometries to be available at grid points. Since this information being not available at all grid points, it is expected that, in the future, textural information can be derived in an indirect way, i.e., from lithologic/textural spatial reconstructions.