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.     

Geostatistical analysis of spatiotemporal variability of groundwater level fluctuations in Amman–Zarqa basin,Jordan: a case study

The main objective of this study was to assess the spatial and temporal variability of groundwater level fluctuations in the Amman–Zarqa basin, during the period 2001–2005. In the year 2003, as a consequence of war, there was a sudden increase in the population in this basin. Knowing that the basin is already heavily populated and witnesses most of the human and industrial activities in Jordan, this study was prompted to help make wise water resources management decisions to cope with the new situation. Data from 31 fairly distributed wells in the upper aquifer of the basin were subjected to geostatistical treatment. Kriging interpolation techniques have indicated that the groundwater flow directions remained almost constant over the years. The two main directions are SW–NE and E–W. Kriging mapped fluctuations have also showed that drop and rise events are localized in the basin. Forecasting possibilities for management purposes were tackled using autocorrelation analysis. The constructed autocorrelograms indicated, in general, the temporal dependence of seasonal water level fluctuations, and that forecasting can be carried out within a period of 3–21 months. Several suggestions were made to mitigate the drop and rise hazards in the detected sites.     

区域地下水位监测网优化设计方法

区域地下水位监测提供了定量评价含水层地下水位持续下降及其对环境影响必不可少的信息。历史上的地下水位监测网是为了评价地下水资源或监测水源地降落漏斗而设立的,目前它们已经不能适应为流域水资源综合管理提供必需的信息。本文在综述国际地下水位监测现状的基础上,介绍了区域地下水位监测网优化设计的方法。采用地理信息系统编制的地下水动态类型图为地下水位监测井位置的选择提供了坚实的水文地质基础;克里金插值法能定量评价监测网观测值绘制的地下水位等高线的精度,因而可以用来定量设计地下水位监测网;时间序列分析和统计检验提供了优化地下水位监测频率的定量标准。这些方法已被应用于北京平原、乌鲁木齐河流域和济南岩溶泉域,其成果将在本刊分期发表。     

塔里木南缘策勒绿洲地下水空间变异性与土地覆盖关系研究

基于策勒绿洲78个取样点的地下水埋深、矿化度和pH值的观测资料, 应用遥感、地理信息系统、 空间插值和地质统计学分析的方法对所取数据进行空间变异性分析.结果表明: 1)地下水埋深和pH值服从正态分布, 矿化度服从对数正态分布;2)地下水埋深、矿化度和pH值都具有强烈的空间相关性, 在步长为8 km范围之内, 地下水特征的空间变异是各向同性的, 当步长＞8 km时, 四个方向上的半变异函数发生了不同的变化;3)受盆地地形影响, 地下水埋深从南到北变浅, 东部区域埋深最浅, 最小为1.56 m;矿化度分布趋势则是从南到北逐渐增大;pH值高值区发生在绿洲东南部和中北部, pH值高的区域矿化度较低, 而pH值低的区域矿化度较高;4) 绿洲耕地主要分布于地下水位埋深为5～25 m, 且矿化度＜2.0 g·L~(-1) 的区域;林地和草地分布于绿洲边缘区, 该区域地下水埋深浅, 平均水位＜10 m, 而地下水矿化度较高, ＞2.0 g·L~(-1) .     

GIS-PModflow联合系统在松嫩平原西部潜水环境预警中的应用

针对松嫩平原西部地下水潜水位变化对环境产生影响的现状,运用GIS-PModflow联合系统确定了地下水环境预警警戒线,采用以不至于发生土壤盐碱化的潜水临界水位作为上警戒线值,以潜水开采极限深度作为下警戒线值。在对系统模型概化的基础上,进行了潜水位的数值模拟和预报。利用GIS空间分析方法对图像数据进行栅格叠加和属性判断,通过对预报水位与警戒水位对比,进行水环境预警。结果表明,研究区2015年的警情较1999年增加,且以灌溉后潜水位上升引起的土壤次生盐渍化的警情为主。     

应用Kriging方法研究格尔木河流域地下水位动态观测网的优化配置

地下水动态是水文地质研究的一个重要方面，如何获满足一定精度要求全面而合理的地下水位动态资料具有现实经济效益。本文研究格尔木河流域地下水位动态测网，通过选取理想的变差函，运有Ｋｒｉｇｉｎｇ方法进行地下水位的线性无偏最优估计和计算估计误差的标准差，结合给定的允许误差限评判地下水位动态观测网的配置是否合适，并提出调整现有观测网的方案     

Control of sea-water intrusion by salt-water pumping: Coast of Oman

A shallow alluvial coastal aquifer in the Batinah area of Oman, with sea-water intrusion that extends several kilometres inland, has been studied experimentally, analytically and numerically. The water table is proved to have a trough caused by intensive pumping from a fresh groundwater zone and evaporation from the saline phreatic surface. Resistivity traverses perpendicular to the shoreline indicated no fresh groundwater recharge into the sea. Using an analytical Dupuit-Forchheimer model, developed for the plain part of the catchment, explicit expressions for the water table, sharp interface location and stored volume of fresh water are obtained. It is shown that by the pumping of salt water from the intruded part of the aquifer, this intrusion can be mitigated. Different catchment sizes, intensities of fresh groundwater pumping, evaporation rates, water densities, sea level, incident fresh water level in the mountains and hydraulic conductivity are considered. SUTRA code is applied to a hypothetical case of a leaky aquifer with line sinks modeling fresh water withdrawal and evaporation. The numerical code also shows that pumping of saline water can pull the dispersion zone back to the shoreline.     

Distance-gradient-based variogram and Kriging to evaluate cobalt-rich crust deposits on seamounts

The spatial distribution of cobalt-rich crust thicknesses on seamounts is partly controlled by water depth and slope gradients. Conventional distance–direction-based variogram have not effectively expressed the spatial self-correlation or anisotropy of the thicknesses of cobalt-rich crusts. To estimate resources in cobalt-rich crusts on seamounts using geostatistics, we constructed a new variogram model to adapt to the spatial distribution of the thicknesses of the cobalt-rich crusts. In this model, we defined the data related to cobalt-rich crusts on seamounts as three-dimensional surface random variables, presented an experimental variogram process based on the distance–gradient or distance–“relative water depth,” and provided a theoretical variogram model that follows this process. This method was demonstrated by the spatial estimation of the thicknesses of cobalt-rich crusts on a seamount, and the results indicated that the new variogram model reflects the spatial self-correlation of the thicknesses of cobalt-rich crusts well. Substituted into the Kriging equation, the new variogram model successfully estimated the spatial thickness distribution of these cobalt-rich crusts.     

岩溶盆地地下水长观监测研究

通过覆盖型盆地岩溶盆地中取水深井1987年-2005年盆地地下水长观资料统计分析,系统探讨了岩溶地下水的水位、取水量以及抽水产生的塌陷分布规律,提出了合理开发利用地下水资源的具体措施.     

1985-2013年黑河中游流域地下水位动态变化特征

在气候变化和人类活动的影响下, 黑河流域地表水和地下水的时空分布特征发生了很大变化. 研究水系统演化及其驱动机制对流域水资源可持续管理非常关键. 基于甘肃河西黑河中游流域地下水位动态、水文气象、土地利用和灌溉统计数据, 研究了1985-2013年黑河中游流域地下水位时空变化. 结果表明: 地表水的不合理分配和耕地的扩展导致了地下水的过量开采和地下水位的剧烈变化. 1985-2004年区域地下水位以下降为主; 2005-2013年呈现下降和回升两极发展趋势, 冲洪积扇群带地下水最大下降达17.41 m, 而黑河干流沿岸地下水位最大回升了3.3 m, 地下水埋深普遍增加了1.0~3.0 m. 尽管地下水位在2005-2013年表现出回升趋势, 但干流中游盆地地下水系统处于严重负均衡状态, 制定合理的“生态分水”方案和水资源综合管理规划非常紧迫.     

Impacts of human activities on the hydrology of Baiyangdian Lake,China

Baiyangdian Lake is the biggest natural freshwater wetland in North China Plain. It provides important ecosystem services such as water regulation and supply, reed production and biodiversity protection. Baiyangdian Lake, however, was threatened by lack of freshwater in recent decades. In this paper, the hydrological changes of the lake were quantified using historical data of water level and groundwater table, and satellite images. In addition, the relationship between water level and socioeconomic development of the basin was investigated. The result revealed a significant decreasing trend of water level in the lake. Water level and groundwater table of Baiyangdian Lake decreased rapidly, caused by the great increase of water withdrawal and consumption due to socioeconomic development in the basin. In particular, population growth and the expansion of irrigated agriculture were two major contributors to the decline of water level and groundwater table. While precipitation was positively correlated with water level, it has less impact on water level and groundwater table than human activities. The diversion of water to the lake raised the water level temporarily and had significant benefits on the wetland ecosystem. The best way to solve ecological problems of Baiyangdian Lake, however, is to control the growth of population, adjust the industrial structure, control land use conversion and improve water use efficiency at the basin scale.     

Groundwater characteristics within loessic deposits: the coastal springs of Los Acantilados,Mar del Plata,Argentina

The urban growth of the southern neighborhoods of Mar del Plata City provoked significant changes in the groundwater balance of the loessic sequences. These regional loessic levels with a significant portion of volcanic ash layers were reported subject to fluoride and nitrate concentrations. Residential houses pump from sands located 70 m depth and withdraw the sewages to depths less than 5 m. These effects cause significant local and seasonal (summer) increments of the water table outcropping via springs at certain unconformities of the coastal cliffs. A mathematical model was applied to analyze the water level lowering at the productive levels, while there is a decrease in the quality of the upper levels subject to waste discharges. Much of this groundwater flow is concentrated in unconformities between different types of sediments. Human activities have affected the aquifer dynamics increasing the groundwater pumping rates and the return velocity of the sewages. This should be considered in the management strategies of coastal hydric resources.     

Spatial assessment of salinity and nitrate pollution in Amman Zarqa Basin: a case study

Intensive use of land resources in arid and semi-arid regions exert serious pressures on groundwater resources and jeopardize further socio-economical developments. The Amman-Zarqa Basin (AZB), the most vital basin in Jordan, is facing recent groundwater deterioration due to a very large increase in water demands for domestic, agricultural, and industrial uses. The objectives of this paper were to quantify the degree of contamination in the basin by evaluating the characteristics, distribution and seasonal variations of two pollution indicators (nitrate concentration and salinity) and to determine the impacts of human activities (land use) on groundwater quality. Based upon long-term data (1970–2005) of groundwater samples collected from 538 wells across the AZB, spatial analyses indicated that both indicators have a strong spatial dependence and are anisotropically distributed. Prediction maps of Ordinary Kriging and Indicator Kriging provided detailed indications of the major and minor sources of pollution in the basin. Inefficient wastewater treatment plants, industrial activities and agricultural practices were responsible for 91, 85, and 25% salinization of nearby wells, respectively. Nitrate pollution had reached 73% above threshold (50 mg/L) in some cases. The temporal analyses estimated the salinity buildup rate to be around 8 × 10⁻² μS/cm per day, while nitrate buildup rate was estimated to be around 6 × 10⁻² mg/L per day. Remote sensing and spatial analyses helped greatly in groundwater quality assessment not only in providing the environmental status of the AZB but also in delineating the potential of contamination risk zones and their correlation to human activities. Furthermore, the paper suggests some environmental protection strategies that should be adopted to protect the vital groundwater resources of the basin from further deterioration.     

岩溶塌陷易发区地下水安全开采控制——以珠三角广花盆地城市应急水源地为例

以珠三角广花盆地城市应急水源地为例,通过野外现场抽水和室内临界水力梯度试验对比,分析了地下水安全开采的控制条件。结果显示:广花盆地的岩溶塌陷是可通过控制水位降深等方法进行防控的,当地下水开采水位降深控制在9 m以内时,可大幅度降低岩溶塌陷概率。该结果可直接作为企业进行地下水安全开采的控制性指标和政府管理部门安全管控的决策性指标。由于地质条件的差异性,在实际开采地下水过程中,不能机械地依赖单一的地下水水位降深控制指标,应根据抽水影响半径以及各地不同的临界水力坡度分地段控制地下水位降深,或者改变开采方式,以达到充分利用地下水资源的目的。建议在实际抽水期间制订突发灾害预案,加强地下水监测工作,实时掌握地下水和其他影响因素变化动态。      

Assessment of alluvial aquifer heterogeneity and development of stochastic hydrofacies models for the Hat Yai Basin in Southern Thailand

In previous studies, the groundwater flow models formulated for the Hat Yai Basin were conventional and deterministic because the geologic heterogeneity of the alluvial aquifer system in the basin had not yet been assessed. This paper describes an effort to develop hydrofacies models, such that the spatial variability of the aquifer system can be represented in a systematic way. Variogram parameters that characterize the alluvial aquifer heterogeneity were determined. Based on these variogram parameters, an indicator-based geostatistical approach was used to develop hydrofacies models using sequential indicator simulation. The hydrofacies models indicate three distinct aquifer units, namely Hat Yai, Khu Tao, and Kho Hong aquifers, which is in good agreement with a conceptual model, and incorporates spatial variability as observed in field data from borehole logs. The hydrofacies models can be used in groundwater modeling and simulations.     

Development of unconfined conditions in multi-aquifer flow systems: a case study in the Rajshahi Barind,Bangladesh

Identifying flow processes in multi-aquifer flow systems is a considerable challenge, especially if substantial abstraction occurs. The Rajshahi Barind groundwater flow system in Bangladesh provides an example of the manner in which flow processes can change with time. At some locations there has been a decrease with time in groundwater heads and also in the magnitude of the seasonal fluctuations. This report describes the important stages in a detailed field and modelling study at a specific location in this groundwater flow system. To understand more about the changing conditions, piezometers were constructed in 2015 at different depths but the same location; water levels in these piezometers indicate the formation of an additional water table. Conceptual models are described which show how conditions have changed between the years 2000 and 2015. Following the formation of the additional water table, the aquifer system is conceptualised as two units. A pumping test is described with data collected during both the pumping and recovery phases. Pumping test data for the Lower Unit are analysed using a computational model with estimates of the aquifer parameters; the model also provided estimates of the quantity of water moving from the ground surface, through the Upper Unit, to provide an input to the Lower Unit. The reasons for the substantial changes in the groundwater heads are identified; monitoring of the recently formed additional water table provides a means of testing whether over-abstraction is occurring.     

Inverse modelling of aquifer parameters in basaltic rock with the help of pumping test method using MODFLOW software

The origin and movement of groundwater are the fundamental questions that address both the temporal and spatial aspects of ground water run and water supply related issues in hydrological systems. As groundwater flows through an aquifer, its composition and temperature may variation dependent on the aquifer condition through which it flows. Thus, hydrologic investigations can also provide useful information about the subsurface geology of a region. But because such studies investigate processes that follow under the Earth's shallow, obtaining the information necessary to answer these questions is not continuously easy. Springs, which discharge groundwater table directly, afford to study subsurface hydrogeological processes.The present study of estimation of aquifer factors such as transmissivity (T) and storativity (S) are vital for the evaluation of groundwater resources. There are several methods to estimate the accurate aquifer parameters (i.e. hydrograph analysis, pumping test, etc.). In initial days, these parameters are projected either by means of in-situ test or execution test on aquifer well samples carried in the laboratory. The simultaneous information on the hydraulic behavior of the well (borehole) that provides on this method, the reservoir and the reservoir boundaries, are important for efficient aquifer and well data management and analysis. The most common in-situ test is pumping test performed on wells, which involves the measurement of the fall and increase of groundwater level with respect to time. The alteration in groundwater level (drawdown/recovery) is caused due to pumping of water from the well. Theis (1935) was first to propose method to evaluate aquifer parameters from the pumping test on a bore well in a confined aquifer. It is essential to know the transmissivity (T = Kb, where b is the aquifer thickness; pumping flow rate, Q = TW (dh/dl) flow through an aquifer) and storativity (confined aquifer: S = bS_s, unconfined: S = S_y), for the characterization of the aquifer parameters in an unknown area so as to predict the rate of drawdown of the groundwater table/potentiometric surface throughout the pumping test of an aquifer. The determination of aquifer's parameters is an important basis for groundwater resources evaluation, numerical simulation, development and protection as well as scientific management. For determining aquifer's parameters, pumping test is a main method. A case study shows that these techniques have been fast speed and high correctness. The results of parameter's determination are optimized so that it has important applied value for scientific research and geology engineering preparation.     

The effects of geological heterogeneities and piezometric fluctuations on groundwater flow and chemistry in a hard-rock aquifer, southern India

Crystalline aquifers of semi-arid southern India represent a vital water resource for farming communities. A field study is described that characterizes the hydrodynamic functioning of intensively exploited crystalline aquifers at local scale based on detailed well monitoring during one hydrological year. The main results show large water-table fluctuations caused by monsoon recharge and pumping, high spatial variability in well discharges, and a decrease of well yields as the water table decreases. Groundwater chemistry is also spatially variable with the existence of aquifer compartments within which mixing occurs. The observed variability and compartmentalization is explained by geological heterogeneities which play a major role in controlling groundwater flow and connectivity in the aquifer. The position of the water table within the fracture network will determine the degree of connectivity between aquifer compartments and well discharge. The presented aquifer conceptual model suggests several consequences: (1) over-exploitation leads to a drop in well discharge, (2) intensive pumping may contribute to the hydraulic containment of contaminants, (3) groundwater quality is highly variable even at local scale, (4) geological discontinuities may be used to assist in the location of drinking-supply wells, (5) modeling should integrate threshold effects due to water-table fluctuations.     

黑河流域中游盆地地下水动态特征及其调蓄能力分析

40 多年来,我国西北地区大规模的地下水开发利用造成了部分地区地下水水位持续下降甚至泉水干涸,部分地区地下水仍能维持动态稳定,判断这类地区水资源开发利用是否具有可持续性是必须解决的重大科学问题。以黑河流域中游盆地作为研究对象,采用MK检验和连续小波分析等方法,分析长时间序列地下水水位数据的变化特点,研究区域地下水动态特征;结合区域水文地质条件,综合划分黑河流域中游盆地地下水补排平衡区与非平衡区;利用克里金插值法估算1990—2020 年盆地含水层对水资源的调节水量,并评价不同动态平衡区的调蓄能力。结果显示,黑河流域中游盆地地下水水位动态类型有:水文型、水文-开采型、开采型和蒸发-开采型4种长周期动态稳定型,过量开采型或上游过度引用地表水型2 种长周期持续下降型。黑河-梨园河倾斜平原、酒泉盆地和黑河中游下段侵蚀堆积平原的大厚度含水层是黑河流域中游盆地的地下水补排平衡区,其在1990—2001 年共输出地下水12.06×108 m3,2001—2020年共储存地下水9.06×108 m3。大厚度含水层为地下水的长周期调蓄提供了充足的空间,在合理控制开采量的前提下,该类含水层的天然调蓄能力可满足生产生活和下游生态用水需要。盆地地下水补排非平衡区,如黑河以东诸河倾斜平原、盐池盆地和榆木山山前诸小河流域等地区,目前的地下水开发利用方式和强度是不可持续的,应适当减少地下水开采量,调节盆地上游的引水量和开采量,抑制地下水资源枯竭。本研究成果可为西北干旱内陆地区水资源管理和持续开发利用提供参考。     

Assessment of aquifer vulnerability to contamination in Khanyounis Governorate, Gaza Strip—Palestine, using the DRASTIC model within GIS environment

Groundwater is a very important natural resource in Khanyounis Governorate (the study area) for water supply and development. Historically, the exploitation of aquifers in Khanyounis Governorate has been undertaken without proper concern for environmental impact. In view of the importance of quality groundwater, it might be expected that aquifer protection to prevent groundwater quality deterioration would have received due attention. In the long term, however, protection of groundwater resources is of direct practical importance because, once pollution of groundwater has been allowed to occur, the scale and persistence of such pollution makes restoration technically difficult and costly. In order to maintain basin aquifer as a source of water for the area, it is necessary to find out, whether certain locations in this groundwater basin are susceptible to receive and transmit contamination. This study aims to: (1) assess the vulnerability of the aquifer to contamination in Khanyounis governorate, (2) find out the groundwater vulnerable zones to contamination in the aquifer of the study area, and (3) provide a spatial analysis of the parameters and conditions under which groundwater may become contaminate. To achieve that, DRASTIC model within geographic information system (GIS) environment was applied. The model uses seven environmental parameters: depth of water table, net recharge, aquifer media, soil media, topography, impact of vadose zone, and hydraulic conductivity to evaluate aquifer vulnerability. Based on this model and by using ArcGIS 9.3 software, an attempt was made to create vulnerability maps for the study area. According to the DRASTIC model index, the study has shown that in the western part of the study area the vulnerability to contamination ranges between high and very high due to the relatively shallow water table with moderate to high recharge potential, and permeable soils. To the east of the previous part and in the south-eastern part, vulnerability to contamination is moderate. In the central and the eastern part, vulnerability to contamination is low due to depth of water table. Vulnerability analysis of the DRASTIC Model indicates that the highest risk of contamination of groundwater in the study area originates from the soil media. The impact of vadose zone, depth to water level, and hydraulic conductivity imply moderate risks of contamination, while net recharge, aquifer media, and topography impose a low risk of aquifer contamination. The coefficient of variation indicates that a high contribution to the variation of vulnerability index is made by the topography. Moderate contribution is made by the depth to water level, and net recharge, while impact of vadose zone, hydraulic conductivity, soil media, and Aquifer media are the least variable parameters. The low variability of the parameters implies a smaller contribution to the variation of the vulnerability index across the study area. Moreover, the “effective” weights of the DRASTIC parameters obtained in this study exhibited some deviation from that of the “theoretical” weights. Soil media and the impact of vadose zone were the most effective parameters in the vulnerability assessment because their mean “effective” weights were higher than their respective “theoretical” weights. The depth of water table showed that both “effective” and “theoretical” weights were equal. The rest of the parameters exhibit lower “effective” weights compared with the “theoretical” weights. This explains the importance of soil media and vadose layers in the DRASTIC model. Therefore, it is important to get the accurate and detailed information of these two specific parameters. The GIS technique has provided an efficient environment for analysis and high capabilities of handling large spatial data. Considering these results, DRASTIC model highlights as a useful tool that can be used by national authorities and decision makers especially in the agricultural areas applying chemicals and pesticides which are most likely to contaminate groundwater resources.