Application and evaluation of kriging and cokriging methods on groundwater depth mapping

Groundwater and water resources management play a key role in conserving the sustainable conditions in arid and semi-arid regions. Applying some techniques that can reveal the critical and hot conditions of water resources seem necessary. In this study, kriging and cokriging methods were evaluated for mapping the groundwater depth across a plain in which has experienced different climatic conditions (dry, wet, and normal) and consequently high variations in groundwater depth in a 12 year led in maximum, minimum, and mean depths. During this period groundwater depth has considerable fluctuations. Results obtained from geostatistical analysis showed that groundwater depth varies spatially in different climatic conditions. Furthermore, the calculated RMSE showed that cokriging approach was more accurate than kriging in mapping the groundwater depth though there was not a distinct difference. As a whole, kriging underestimated the real groundwater depth for dry, wet, and normal conditions by 5.5, 2.2, and 5.3%, while cokriging underestimations were 3.3, 2, and 2.2%, respectively; which showed the unbiasedness in estimations. Results implied that in the study area farming and cultivation in dry conditions needs more attention due to higher variability in groundwater depth in short distances compared to the other climate conditions. It is believed that geostatistical approaches are reliable tools for water resources managers and water authorities to allocate groundwater resources in different environmental conditions.     

SPATIAL ANALYSIS AND EUTROPHICATION ASSESSMENT FOR CHLOROPHYLL a IN TAIHU LAKE

Spatial structure analysis and kriging analysis have been identified to be useful tools in illustrating the spatial patterns of variables. Taihu Lake is one of the largest fresh water lakes in China, and has suffered serious eutrophication in recent years due to the rapid economic development and growing environmental pollution in the Taihu Catchment. In this paper, spatial structural analysis, kriging interpolation and eutrophication assessment were carried out for chlorophyll a in the lake. Studies show that spherical model could be applied to fit all experimental variograms. Positive nuggets were observed for three directions except NE–SW direction. The variograms show some anisotropy with anisotropic ratio falling within 1.76. The spatial structural patterns of chlorophyll a in the lake were affected by factors such as distribution of pollution sources, water flow and wind. Two-dimensional ordinary block kriging was applied for interpolation process. An eutrophication assessment map was also made based on a water-quality evaluation standard. Results show that the content of chlorophyll a in Taihu Lake was quite high. The whole lake has suffered serious eutrophication. However, the eutrophic situation varied in space. Higher contents of chlorophyll a appeared mainly in the northern part of the lake.     

Groundwater levels time series sensitivity to pluviometry and air temperature: a geostatistical approach to Sfax region,Tunisia

In this paper, the pattern of groundwater level fluctuations is investigated by statistical techniques for 24 monitoring wells located in an unconfined coastal aquifer in Sfax (Tunisia) for a time period from 1997 to 2006. Firstly, a geostatistical study is performed to characterize the temporal behaviors of data sets in terms of variograms and to make predictions about the value of the groundwater level at unsampled times. Secondly, multivariate statistical methods, i.e., principal component analysis (PCA) and cluster analysis (CA) of time series of groundwater levels are used to classify groundwater hydrographs regard to identical fluctuation pattern. Three groundwater groups (A, B, and C) were identified. In group “A,” water level decreases continuously throughout the study periods with rapid annual cyclic variation, whereas in group “B,” the water level contains much less high-frequency variation. The wells of group “C” represents a steady and gradual increase of groundwater levels caused by the aquifer artificial recharge. Furthermore, a cross-correlation analysis is used to investigate the aquifer response to local rainfall and temperature records. The result revealed that the temperature is more affecting the variation of the groundwater level of group A wells than the rainfall. However, the second and the third groups are less affected by rainfall or temperature.     

Indicator and probability kriging methods for delineating Cu, Fe, and Mn contamination in groundwater of Najafgarh Block, Delhi, India

Two non-parametric kriging methods such as indicator kriging and probability kriging were compared and used to estimate the probability of concentrations of Cu, Fe, and Mn higher than a threshold value in groundwater. In indicator kriging, experimental semivariogram values were fitted well in spherical model for Fe and Mn. Exponential model was found to be best for all the metals in probability kriging and for Cu in indicator kriging. The probability maps of all the metals exhibited an increasing risk of pollution over the entire study area. Probability kriging estimator incorporates the information about order relations which the indicator kriging does not, has improved the accuracy of estimating the probability of metal concentrations in groundwater being higher than a threshold value. Evaluation of these two spatial interpolation methods through mean error (ME), mean square error (MSE), kriged reduced mean error (KRME), and kriged reduced mean square error (KRMSE) showed 3.52% better performance of probability kriging over indicator kriging. The combined result of these two kriging method indicated that on an average 26.34%, 65.36%, and 99.55% area for Cu, Fe, and Mn, respectively, are coming under the risk zone with probability of exceedance from a cutoff value is 0.6 or more. The groundwater quality map pictorially represents groundwater zones as ??desirable?? or ??undesirable?? for drinking. Thus the geostatistical approach is very much helpful for the planners and decision makers to devise policy guidelines for efficient management of the groundwater resources so as to enhance groundwater recharge and minimize the pollution level.     

Spatial assessment of water quality parameters in Jhelum city (Pakistan)

In this study, we assess the drinking water quality of Jhelum city. Two hundred and ninety-two drinking water samples were randomly collected in the study area. These samples were chemically analyzed for three key toxic (in excess) elements such as pH, total dissolved solids (TDS), and calcium. Geostatistical techniques such as variogram and kriging were used to investigate the spatial variations of these minerals across the city. The spatial structure for each element was found to be anisotropic, and thus, anisotropic variograms were used. The kriging predictions revealed significant concentrations of the above-stated elements at some locations in the study area. While comparing with the World Health Organization, United States Environmental Protection Agency, and Pakistan Environmental Protection Agency standards, the water samples were found to be unsatisfactory for drinking. We conclude that the drinking water in this region is of poor quality and needs proper treatment to make it palatable.     

Spatial evaluation of the risk of groundwater quality degradation. A comparison between disjunctive kriging and geostatistical simulation

In some previous papers a probabilistic methodology was introduced to estimate a spatial index of risk of groundwater quality degradation, defined as the conditional probability of exceeding assigned thresholds of concentration of a generic chemical sampled in the studied water system. A crucial stage of this methodology was the use of geostatistical techniques to provide an estimation of the above-mentioned probability in a number of selected points by crossing spatial and temporal information. In this work, spatial risk values were obtained using alternatively stochastic conditional simulation and disjunctive kriging. A comparison between the resulting two sets of spatial risks, based on global and local statistical tests, showed that they do not come from the same statistical population and, consequently, they cannot be viewed as equivalent in a statistical sense. At a first glance, geostatistical conditional simulation may appear to represent the spatial variability of the phenomenon more effectively, as the latter tends to be smoothed by DK. However, a close examination of real case study results suggests that disjunctive kriging is more effective than simulation in estimating the spatial risk of groundwater quality degradation. In the study case, the potentially ‘harmful event’ considered, threatening a natural ‘vulnerable groundwater system,’ is fertilizer and manure application.     

Comparison of stochastic and deterministic methods for mapping groundwater level spatial variability in sparsely monitored basins

In sparsely monitored basins, accurate mapping of the spatial variability of groundwater level requires the interpolation of scattered data. This paper presents a comparison of deterministic interpolation methods, i.e. inverse distance weight (IDW) and minimum curvature (MC), with stochastic methods, i.e. ordinary kriging (OK), universal kriging (UK) and kriging with Delaunay triangulation (DK). The study area is the Mires Basin of Mesara Valley in Crete (Greece). This sparsely sampled basin has limited groundwater resources which are vital for the island’s economy; spatial variations of the groundwater level are important for developing management and monitoring strategies. We evaluate the performance of the interpolation methods with respect to different statistical measures. The Spartan variogram family is applied for the first time to hydrological data and is shown to be optimal with respect to stochastic interpolation of this dataset. The three stochastic methods (OK, DK and UK) perform overall better than the deterministic counterparts (IDW and MC). DK, which is herein for the first time applied to hydrological data, yields the most accurate cross-validation estimate for the lowest value in the dataset. OK and UK lead to smooth isolevel contours, whilst DK and IDW generate more edges. The stochastic methods deliver estimates of prediction uncertainty which becomes highest near the southeastern border of the basin.     

Estimation of spatial distrubition of groundwater level and risky areas of seawater intrusion on the coastal region in Çarşamba Plain,Turkey, using different interpolation methods

Groundwater level plays a significant role in coastal plains. Heavy pumping and excessive use of near-coast groundwater can increase the intrusion of seawater into the aquifers. In the present study, groundwater levels were measured at 59 groundwater wells at different times during pre- and post-irrigation seasons (April and September of the year 2012) in Çar?amba Plain, Turkey. To select the best method, two deterministic interpolation methods (inverse distance weighing (IDW) with the weights of 1, 2, and 3 and radial basis function (RBF) with spline with tension (SPT) and completely regularized spline (CRS)) and two stochastic methods (ordinary kriging (OK) with spherical, exponential, and Gaussian variograms) and cokriging (COK)) were compared and then the best interpolation method was used to evaluate the spatial distribution of groundwater levels in different seasons and seasonal changes. A total of nine different techniques were tested. Also, risky areas of seawater intrusion in coastal area were determined using the best methods for two periods. The performance of these interpolation methods is evaluated by using a validation test method. Statistical indices of correlation (R ²), mean absolute error (MAE), and root-mean-square error (RMSE) were used to select and validate the best methods. Comparisons between predicted and observed values indicated RBF as the optimal method for groundwater level estimation in April and September. When the best method RBF and the worst method IDW were compared, significant differences were observed in the spatial distribution of groundwater. Results of the study also revealed that excessive groundwater withdrawals during the post-irrigation season dropped the groundwater levels up to 2.0 m in some sections. With regard to seawater intrusion, 9,103 ha of land area was determined to be highly risky and risky.     

Comparison of Groundwater Level Estimation Using Neuro-fuzzy and Ordinary Kriging

Water level in aquifer plays the main role in groundwater modeling as one of the input data. In practice, due to aspects of time and cost, data monitoring of water levels is conducted at a limited number of sites, and interpolation technique such as kriging is widely used for estimation of this variable in unsampled sites. In this study, the efficiency of the ordinary kriging (OK) and adaptive network-based fuzzy inference system (ANFIS) was investigated in interpolation of groundwater level in an unconfined aquifer in the north of Iran. The results showed that ANFIS model is more efficient in estimating the groundwater level than OK.     

GIS-based assessment and characterization of groundwater quality in a hard-rock hilly terrain of Western India

The growing population, pollution, and misuse of freshwater worldwide necessitate developing innovative methods and efficient strategies to protect vital groundwater resources. This need becomes more critical for arid/semi-arid regions of the world. The present study focuses on a GIS-based assessment and characterization of groundwater quality in a semi-arid hard-rock terrain of Rajasthan, western India using long-term and multi-site post-monsoon groundwater quality data. Spatio-temporal variations of water quality parameters in the study area were analyzed by GIS techniques. Groundwater quality was evaluated based on a GIS-based Groundwater Quality Index (GWQI). A Potential GWQI map was also generated for the study area following the Optimum Index Factor concept. The most-influential water quality parameters were identified by performing a map removal sensitivity analysis among the groundwater quality parameters. Mean annual concentration maps revealed that hardness is the only parameter that exceeds its maximum permissible limit for drinking water. GIS analysis revealed that sulfate and nitrate ions exhibit the highest (CV?>?30%) temporal variation, but groundwater pH is stable. Hardness, EC, TDS, and magnesium govern the spatial pattern of the GWQI map. The groundwater quality of the study area is generally suitable for drinking and irrigation (median GWQI?>?74). The GWQI map indicated that relatively high-quality groundwater exists in northwest and southeast portions of the study area. The groundwater quality parameter group of Ca, Cl, and pH were found to have the maximum value (6.44) of Optimum Index factor. It is concluded that Ca, Cl, and pH are three prominent parameters for cost-effective and long-term water quality monitoring in the study area. Hardness, Na, and SO₄, being the most-sensitive water quality parameters, need to be monitored regularly and more precisely.