The geostatistic-based spatial distribution variations of soil salts under long-term wastewater irrigation

The purpose of this study was to determine and evaluate the spatial changes in soil salinity by using geostatistical methods. The study focused on the suburb area of Beijing, where urban development led to water shortage and accelerated wastewater reuse to farm irrigation for more than 30 years. The data were then processed by GIS using three different interpolation techniques of ordinary kriging (OK), disjunctive kriging (DK), and universal kriging (UK). The normality test and overall trend analysis were applied for each interpolation technique to select the best fitted model for soil parameters. Results showed that OK was suitable for soil sodium adsorption ratio (SAR) and Na⁺ interpolation; UK was suitable for soil Cl^? and pH; DK was suitable for soil Ca²⁺. The nugget-to-sill ratio was applied to evaluate the effects of structural and stochastic factors. The maps showed that the areas of non-saline soil and slight salinity soil accounted for 6.39 and 93.61 %, respectively. The spatial distribution and accumulation of soil salt were significantly affected by the irrigation probabilities and drainage situation under long-term wastewater irrigation.     

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.     

Spatial assessment of soil salinity in the Harran Plain using multiple kriging techniques

The Harran Plain is located in the southeastern part of Turkey and has recently been developed for irrigation agriculture. It already faces soil salinity problems causing major yield losses. Management of the problem is hindered by the lack of information on the extent and geography of the salinization problem. A survey was carried out to delineate the spatial distribution of salt-affected areas by randomly selecting 140 locations that were sampled at two depths (0 to 30 and 30 to 60 cm) and analyzed for soil salinity variables: soil electrical conductivity (EC), soluble cations (Ca^2+, Mg²⁺, Na⁺, and K⁺), soluble anions (SO ₄ ^2? , Cl^?), exchangeable Na⁺ (me 100 g^?1) and exchangeable sodium percentage. Terrain attributes (slope, topographical wetness index) were extracted from the digital elevation model of the study area. Variogram analyses after log transformation and ordinary kriging (OK) were applied to map spatial patterns of soil salinity variables. Multivariate geostatistical methods—regression kriging (RK) and kriging with external drift (KED)—were used using elevation and soil electrical conductivity data as covariates. Performances of the three estimation methods (OK, RK, and KED) were compared using independent validation samples randomly selected from the main dataset. Soils were categorized into salinity classes using disjunctive kriging (DK) and ArcGIS, and classification accuracy was tested using the kappa statistic. Results showed that soil salinity variables all have skewed distribution and are poorly correlated with terrain indices but have strong correlations among each other. Up to 65 % improvement was obtained in the estimations of soil salinity variables using hybrid methods over OK with the best estimations obtained with RK using EC_0–30 as covariate. DK–ArcGIS successfully classified soil samples into different salinity groups with overall accuracy of 75 % and kappa of 0.55 (p?<?0.001).     

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.     

山区土壤重金属空间插值性能的改进

山区土壤重金属监测中,密集均匀的布点采样困难大,成本高。根据稀疏非均匀样本数据准确预测山区土壤重金属空间分布,对可视化描述检测元素的分布及趋势具有重要意义。基于云南省楚雄市南部山区表层土壤中重金属锰、钒的采样数据,对集成径向基函数神经网络(IRBFANNs)和传统插值方法:反距离(IDW)、普通克里格(OK)、径向基函数神经网络(RBFANNs),进行了不同等级采样密度下的插值性能比较。结果表明,样本点数量减少时,传统预测方法的插值精度都有所下降,而IRBFANNs算法在样本点较少情况下能够保持最好的插值精确度和稳定性,有效改进了空间插值性能。     

Groundwater depth and elevation interpolation by kriging methods in Mohr Basin of Fars province in Iran

Prediction of groundwater depth and elevation is important in quantitative water management especially in arid areas. There are several basins in southwest of Iran, in Zagross Mountain, in which the water wells are distributed along a narrow elliptic ring band around the region. To find the most applicable interpolation method, both of the groundwater depth and elevation are predicted by different kriging methods. It is found that the groundwater elevation and depth can be predicted by different methods. Furthermore, it is found that the methods in which the trend is eliminated predicted the groundwater elevation and depth in central part of the region is with less standard error. Furthermore, the methods with no trend elimination, predicted the groundwater depths with less error near the water wells. Dividing the area to hydro-geologically homogeneous sub-areas improved the interpolation precision.     

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 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.     

Comparison of ordinary and lognormal kriging on skewed data of total cadmium in forest soils of Sweden

Spatial statistical analysis of georeferenced data of total cadmium (TCd) in forest soils of Sweden was assumed to providemore advantageous maps than traditional interpolated maps. However, 264 measurements of TCd in O-horizon of forest soils displayed skewed frequency distribution. Since atypicalobservations affect badly the variogram, outliers wereidentified, different data transformations were tested andordinary (OK) and lognormal kriging (LK) scenarios werecompared based on cross-validation. Results were comparedusing overall measures of predictors, e.g. traditionalmean squared prediction error (MSPE), mean of krigingvariances, variance ratio, median of internallystandardised residuals, and assessments of classificationaccuracy, such as percentage of correctly predictedsamples and within-class MSPE.One outlier was identified based on the absolute value of skewness of value differences less or equal to one in data pairs separated at certain lag classes. Mapping categories characterised by percentage of correct classification and within-class MSPE were found to be essential in comparison of kriging results additionally to the overall measures. In comparison of kriging methods, OK predicted high values moreaccurately and LK was more effective to predict low and mediumvalues. Thus, OK was suggested for mapping high concentration of TCd and other pollutants. Percentage of correctly predictedsamples and within-class MSPE were found to be dependent on kriging method, as well as on the number and limits of categories.     

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.     

Groundwater quality assessment and its correlation with gastroenteritis using GIS: a case study of Rawal Town,Rawalpindi, Pakistan

Majority of the people of Pakistan get drinking water from groundwater source. Nearly 40 % of the total ailments reported in Pakistan are the result of dirty drinking water. Every summer, thousands of patients suffer from acute gastroenteritis in the Rawal Town. Therefore, a study was designed to generate a water quality index map of the Rawal Town and identify the relationship between bacteriological water quality and socio-economic indicators with gastroenteritis in the study area. Water quality and gastroenteritis patient data were collected by surveying the 262 tubewells and the major hospitals in the Rawal Town. The collected spatial data was analyzed by using ArcGIS spatial analyst (Moran’s I spatial autocorrelation) and geostatistical analysis tools (inverse distance weighted, radial basis function, kriging, and cokriging). The water quality index (WQI) for the study area was computed using pH, turbidity, total dissolved solids, calcium, hardness, alkalinity, and chloride values of the 262 tubewells. The results of Moran’s I spatial autocorrelation showed that the groundwater physicochemical parameters were clustered. Among IDW, radial basis function, and kriging and cokriging interpolation techniques, cokriging showed the lowest root mean square error. Cokriging was used to make the spatial distribution maps of water quality parameters. The WQI results showed that more than half of the tubewells in the Rawal Town were providing “poor” to “unfit” drinking water. The Pearson’s coefficient of correlation for gastroenteritis with fecal coliform was found significant (P?P?P?相似文献   

Improved mapping of National Atmospheric Deposition Program wet-deposition in complex terrain using PRISM-gridded data sets

High-elevation regions in the United States lack detailed atmospheric wet-deposition data. The National Atmospheric Deposition Program/National Trends Network (NADP/NTN) measures and reports precipitation amounts and chemical constituent concentration and deposition data for the United States on annual isopleth maps using inverse distance weighted (IDW) interpolation methods. This interpolation for unsampled areas does not account for topographic influences. Therefore, NADP/NTN isopleth maps lack detail and potentially underestimate wet deposition in high-elevation regions. The NADP/NTN wet-deposition maps may be improved using precipitation grids generated by other networks. The Parameter-elevation Regressions on Independent Slopes Model (PRISM) produces digital grids of precipitation estimates from many precipitation-monitoring networks and incorporates influences of topographical and geographical features. Because NADP/NTN ion concentrations do not vary with elevation as much as precipitation depths, PRISM is used with unadjusted NADP/NTN data in this paper to calculate ion wet deposition in complex terrain to yield more accurate and detailed isopleth deposition maps in complex terrain. PRISM precipitation estimates generally exceed NADP/NTN precipitation estimates for coastal and mountainous regions in the western United States. NADP/NTN precipitation estimates generally exceed PRISM precipitation estimates for leeward mountainous regions in Washington, Oregon, and Nevada, where abrupt changes in precipitation depths induced by topography are not depicted by IDW interpolation. PRISM-based deposition estimates for nitrate can exceed NADP/NTN estimates by more than 100% for mountainous regions in the western United States.     

全国近岸海域水质空间插值算法精度分析

空间插值技术对近海环境数据的分析和表征具有广泛应用价值,但目前针对各类模型适用性和精度比较的研究很少。基于2009年全国近岸海域299个站位的水质监测数据,利用交叉检验法,对3种插值算法13种模型的空间拟合结果进行分析和验证。结果表明,13种模型的模拟值与实测值显著相关(r=0.643~0.750,P0.01,n=299),说明其均适用于近岸海域水质的空间模拟;综合比较模拟结果与实测结果的相关系数、平均误差、平均绝对误差、均方差及偏离程度,排序为IDW(p2)Spline-TensionIDW(p3)IDW(p4)IDW(p1)IDW(p5)IDW(p6)Kriging-GuassianKriging-Spherical,Kriging-CircularKriging-LinearKriging-ExponentialSpline-Regularized,提出IDW(p2)是精度最优的空间插值模型。     

Estimating oil spill characteristics from oil heads in scattered monitoring wells

The results are presented of a comparison of four different methods of interpolating observed hydrocarbon depths in monitoring wells, as well as a comparison of different methods of selecting sampling points for interpolation. The results provide criteria for selecting one interpolation method over another, under different scenarios. The methods analyzed are: (1) inverse-distance weighting; (2) punctual kriging; (3) minimum surface curvature; and (4) Akima's quintic polynomial.     

Assessment of groundwater pollution in West Delhi, India using geostatistical approach

The exploration, exploitation, and unscientific management of groundwater resources in the National Capital Territory (NCT) of Delhi, India have posed a serious threat of reduction in quantity and deterioration of quality. The objective of the study is to determine the groundwater quality and to assess the risk of groundwater pollution at Najafgarh, NCT of Delhi. The groundwater quality parameters were analyzed from the existing wells of the Najafgarh and the thematic maps were generated using geostatistical concepts. Ordinary kriging and indicator kriging methods were used as geostatistical approach for preparation of thematic maps of the groundwater quality parameters such as bicarbonate, calcium, chloride, electrical conductivity (EC), magnesium, nitrate, sodium, and sulphate with concentrations equal or greater than their respective groundwater pollution cutoff value. Experimental semivariogram values were fitted well in spherical model for the water quality parameters, such as bicarbonate, chloride, EC, magnesium, sodium, and sulphate and in exponential model for calcium and nitrate. The thematic maps of all the groundwater quality parameters exhibited an increasing trend of pollution from the northern and western part of the study area towards the southern and eastern part. The concentration was highest at the southernmost part of the study area but it could not reflect correctly the groundwater pollution status. The indicator kriging method is useful to assess the risk of groundwater pollution by giving the conditional probability of concentrations of different chemical parameters exceeding their cutoff values. Thus, risk assessment of groundwater pollution is useful for proper management of groundwater resources and minimizing the pollution threat.     

Assessment of methane emission and oxidation at Air Hitam Landfill site cover soil in wet tropical climate

Methane (CH₄) emissions and oxidation were measured at the Air Hitam sanitary landfill in Malaysia and were modeled using the Intergovernmental Panel on Climate Change waste model to estimate the CH₄ generation rate constant, k. The emissions were measured at several locations using a fabricated static flux chamber. A combination of gas concentrations in soil profiles and surface CH₄ and carbon dioxide (CO₂) emissions at four monitoring locations were used to estimate the CH₄ oxidation capacity. The temporal variations in CH₄ and CO₂ emissions were also investigated in this study. Geospatial means using point kriging and inverse distance weight (IDW), as well as arithmetic and geometric means, were used to estimate total CH₄ emissions. The point kriging, IDW, and arithmetic means were almost identical and were two times higher than the geometric mean. The CH₄ emission geospatial means estimated using the kriging and IDW methods were 30.81 and 30.49 g m^?2 day^?1, respectively. The total CH₄ emissions from the studied area were 53.8 kg day^?1. The mean of the CH₄ oxidation capacity was 27.5 %. The estimated value of k is 0.138 year^?1. Special consideration must be given to the CH₄ oxidation in the wet tropical climate for enhancing CH₄ emission reduction.     

Mapping the Spatial Variability of Plant Diversity in a Tropical Forest: Comparison of Spatial Interpolation Methods

Knowledge of the spatial distribution of plant species is essential to conservation and forest managers in order to identify high priority areas such as vulnerable species and habitats, and designate areas for reserves, refuges and other protected areas. A reliable map of the diversity of plant species over the landscape is an invaluable tool for such purposes. In this study, the number of species, the exponent Shannon and the reciprocal Simpson indices, calculated from 141 quadrat sites sampled in a tropical forest were used to compare the performance of several spatial interpolation techniques used to prepare a map of plant diversity, starting from sample (point) data over the landscape. Means of mapped classes, inverse distance functions, kriging and co-kriging, both, applied over the entire studied landscape and also applied within vegetation classes, were the procedures compared. Significant differences in plant diversity indices between classes demonstrated the usefulness of boundaries between vegetation types, mapped through satellite image classification, in stratifying the variability of plant diversity over the landscape. These mapped classes, improved the accuracy of the interpolation methods when they were used as prior information for stratification of the area. Spatial interpolation by co-kriging performed among the poorest interpolators due to the poor correlation between the plant diversity variables and vegetation indices computed by remote sensing and used as covariables. This indicated that the latter are not suitable covariates of plant diversity indices. Finally, a within-class kriging interpolator yielded the most accurate estimates of plant diversity values. This interpolator not only provided the most accurate estimates by accounting for the indices' intra-class variability, but also provided additional useful interpretations of the structure of spatial variability of diversity values through the interpretation of their semi-variograms. This additional role was found very useful in aiding decisions in conservation planning.     

土壤重金属污染空间插值及其验证方法研究综述

综述了土壤重金属污染空间插值方法研究现状,分析了3类典型的插值方法(确定性插值方法、地统计学插值方法和组合插值方法),以及2种常用的精度验证方法(交叉验证法和独立的数据集验证法)各自的优缺点与适用范围,提出了加强组合插值模型优化与精度提升、利用计算机智能技术优化插值算法模型、开发针对土壤重金属领域的考虑空间插值软数据挖...     

Saltwater intrusion in coastal aquifer ??? comparison between the CIP and MOC simulation technique

Proper management and evaluation of environmental effects for coastal groundwater aquifers require accurate estimation of salt and fresh water movement. For this we compare two numerical methods to study saltwater intrusion. The two numerical methods are the CIP (constrained interpolation profile) and the MOC (method of characteristics). The two numerical methods are used to solve the salt transport equation. The two methods are applied and compared to saltwater intrusion experiments for confined and unconfined aquifer case studies. It is found that for both cases, the CIP-method is superior in terms of calculation time while keeping an acceptable numerical accuracy.     

基于反距离加权插值法评价海域水质类别空间分布

空间插值技术可以对污染物分布状况进行面状连续表征,是对点位法评价的有效补充。笔者根据现行海水水质评价方法,类比给出近岸海域水质类别空间分布评价的流程:单指标空间插值-单指标水质分类-多指标水质类别判定。研究了反距离加权插值参数的选取方法,并以点位密度不同的全国近岸海域(开放型)和渤海海域(近封闭型)为例进行了插值实验。 结果表明:对目前的点位分布情况而言,全国近岸海域参数设为p=3,N=4, C=200m较为适宜;渤海海域参数为p=2,N=5,C=200m较为适宜。反距离加权插值法在开放海域的拟合效果要优于近封闭海域。