Application of WetSpass model to estimate groundwater recharge variability in the Nile Delta aquifer

Long-term groundwater recharge from rainfall in the Nile Delta is needed as an input for integrated groundwater modelling in the Nile Delta aquifer for more accurate simulation. The main objective is to estimate the spatial and temporal variation of groundwater recharge from rainfall in the Nile Delta aquifer. Water and Energy Transfer between Soil, Plants and Atmosphere under quasi-Steady State (WetSpass) model parameters were identified for the Nile Delta based on the available meteorological data for the area collected in 1991 and 2000. The collected data were rainfall, temperature, wind speed and evapotranspiration. Geomorphological characteristics, such as soil type, topography, groundwater depth and slope, were also collected as input data for the WetSpass model. ENVI software was used to come up with land use classification based on available land cover images of the Nile Delta for 1972, 1984, 1990, 2000 and 2009. The WetSpass model was calibrated by comparing the simulated groundwater recharge with the calculated one by using the water balance equation model. The results indicated close agreement in groundwater recharge between the two model outputs with R ² of 0.99 and 0.94, while the root-mean-square errors (RMSEs) were 4.86 and 9.39 mm for 1991 and 2000, respectively. The WetSpass model was then applied in respect of 1970, 1980, 1990 and 2010 for the purpose of validation. The overall RMSE and R ² for the 6 years were 8.83 mm and 0.88, respectively. The results of the WetSpass calibrated model provide information to support integrated groundwater modelling. The results reveal that WetSpass works well in simulating the components of the hydrological balance in the Nile Delta.     

Estimating groundwater recharge using GIS-based distributed water balance model in an environmental protection area in the city of Sete Lagoas (MG), Brazil

Improvement in modern water resource management has become increasingly reliant on better characterizing of the spatial variability of groundwater recharge mechanisms. Due to the flexibility and reliability of GIS-based index models, they have become an alternative for mapping and interpreting recharge systems. For this reason, an index model by integrating water balance parameters (surface runoff, actual evapotranspiration, and percolation) calculated by Thornthwaite and Mather’s method, with maps of soil texture, land cover, and terrain slope, was developed for a sustainable use of the groundwater resources. The Serra de Santa Helena Environmental Protection Area, next to the urbanized area of Sete Lagoas (MG), Brazil, was selected as the study area. Rapid economic growth has led to the subsequent expansion of the nearby urban area. Large variability in soil type, land use, and slope in this region resulted in spatially complex relationships between recharge areas. Due to these conditions, the study area was divided into four zones, according to the amount of recharge: high (>?100 mm/year), moderate (50–100 mm/year), low (25–50 mm/year), and incipient (>?25 mm/year). The technique proved to be a viable method to estimate the spatial variability of recharge, especially in areas with little to no in situ data. The success of the tool indicates it can be used for a variety of groundwater resource management applications.     

Simulation of spatial and temporal trends in nitrate concentrations at the regional scale in the Upper Dyle basin,Belgium

Models are the only tools capable of predicting the evolution of groundwater systems at a regional scale, by taking into account a large amount of information. This study presents the association of a water balance model (WetSpass) with a groundwater flow and solute transport model (SUFT3D, saturated and unsaturated flow and transport in 3D) in order to simulate the present and future groundwater quality in terms of nitrate in the Upper Dyle basin (439 km²) Belgium. The HFEMC (hybrid finite element mixing cell) method implemented in the SUFT3D code is used to model groundwater flow and nitrate transport. Spatially distributed recharge, modelled with WetSpass, is considered for prescribing the recharge to the groundwater flow model. The feasibility of linking the WetSpass model with the finite-elements SUFT3D code is demonstrated. Time evolution and distribution of nitrate concentration are then simulated using the calibrated model. Nitrate inputs are spatially distributed according to land use. The spatial simulations and temporal trends are compared with previously published data on this aquifer and show good results.     

Holistic approach of GIS based Multi-Criteria Decision Analysis (MCDA) and WetSpass models to evaluate groundwater potential in Gelana watershed of Ethiopia

Appropriate quantification and identification of the groundwater distribution in a hydrological basin may provide necessary information for effective management, planning and development of groundwater resources. Groundwater potential assessment and delineation in a highly heterogeneous environment with limited Spatiotemporal data derived from Gelana watershed of Abaya Chamo lake basin is performed, using integrated multi-criteria decision analysis (MCDA), water and energy transfer between soil and plant and atmosphere under quasi-steady state (WetSpass) models. The outputs of the WetSpass model reveal a favorable structure of water balance in the basin studied, mainly using surface runoff. The simulated total flow and groundwater recharge are validated using river measurements and estimated baseflow at two gauging stations located in the study area, which yields a good agreement. The WetSpass model effectively integrates a water balance assessment in a geographical information system (GIS) environment. The WetSpass model is shown to be computationally reputable for such a remote complex setting as the African rift, with a correlation coefficient of 0.99 and 0.99 for total flow and baseflow at a significant level of p-value<0.05, respectively. The simulated annual water budget reveals that 77.22% of annual precipitation loses through evapotranspiration, of which 16.54% is lost via surface runoff while 6.24% is recharged to the groundwater. The calibrated groundwater recharge from the WetSpass model is then considered when determining the controlling factors of groundwater occurrence and formation, together with other multi-thematic layers such as lithology, geomorphology, lineament density and drainage density. The selected five thematic layers through MCDA are incorporated by employing the analytical hierarchy process (AHP) method to identify the relative dominance in groundwater potential zoning. The weighted factors in the AHP are procedurally aggregated, based on weighted linear combinations to provide the groundwater potential index. Based on the potential indexes, the area then is demarcated into low, moderate, and high groundwater potential zones (GWPZ). The identified GWPZs are finally examined using the existing groundwater inventory data (static water level and springs) in the region. About 70.7% of groundwater inventory points are coinciding with the delineated GWPZs. The weighting comparison shows that lithology, geomorphology, and groundwater recharge appear to be the dominant factors influence on the resources potential. The assessment of groundwater potential index values identify 45.88% as high, 39.38% moderate, and 14.73% as low groundwater potential zones. WetSpass model analysis is more preferable in the area like Gelana watershed when the topography is rugged, inaccessible and having limited gauging stations.     

Comparing approaches for modeling spatially distributed direct recharge in a semi-arid region (Okanagan Basin, Canada)

Spatially distributed recharge is compared at two different scales using three different modeling approaches within the semi-arid Okanagan Basin, British Columbia, Canada. Regional recharge was modeled by mapping results for one-dimensional soil columns from the water-balance code HELP (Hydrologic Evaluation of Landfill Performance, V3.80D). The regional model was then compared to two, independently derived, local-scale models to ensure local trends were captured in the regional model, and to compare modeling methods. Average annual recharge, predicted by the regional model, varied from no recharge to 186 mm/yr. For the north Okanagan (Vernon area), regional estimates were compared to Richards’ equation-based MIKE-SHE (V2007) estimates, which showed a significant difference in average annual recharge: 7 mm/yr (MIKE-SHE) and 109 mm/yr (HELP). In the south Okanagan (Oliver area), regional estimates were compared to high-resolution, local HELP estimates. Similar values of average annual recharge were obtained: 34 mm/yr (local) and 42 mm/yr (regional). A comparison with measured actual evapotranspiration data in the north Okanagan, showed HELP over-predicted recharge compared to MIKE-SHE by under-predicting evapotranspiration during summer months. Thus, the use of HELP in semi-arid areas may be limited if accurate estimates of recharge are needed. However, results may give satisfactory groundwater model calibrations results because of high uncertainty in hydraulic properties.     

Effects of climate change on the groundwater system in the Grote-Nete catchment,Belgium

The effects of climate change on the groundwater systems in the Grote-Nete catchment, Belgium, covering an area of 525 km², is modeled using wet (greenhouse), cold or NATCC (North Atlantic Thermohaline Circulation Change) and dry climate scenarios. Low, central and high estimates of temperature changes are adopted for wet scenarios. Seasonal and annual water balance components including groundwater recharge are simulated using the WetSpass model, while mean annual groundwater elevations and discharge are simulated with a steady-state MODFLOW groundwater model. WetSpass results for the wet scenarios show that wet winters and drier summers are expected relative to the present situation. MODFLOW results for wet high scenario show groundwater levels increase by as much as 79 cm, which could affect the distribution and species richness of meadows. Results obtained for cold scenarios depict drier winters and wetter summers relative to the present. The dry scenarios predict dry conditions for the whole year. There is no recharge during the summer, which is mainly attributed to high evapotranspiration rates by forests and low precipitation. Average annual groundwater levels drop by 0.5 m, with maximum of 3.1 m on the eastern part of the Campine Plateau. This could endanger aquatic ecosystem, shrubs, and crop production.     

Estimation of groundwater recharge using water balance coupled with base-flow-record estimation and stable-base-flow analysis

In this paper, the long-term mean annual groundwater recharge of Taiwan is estimated with the help of a water-balance approach coupled with the base-flow-record estimation and stable-base-flow analysis. Long-term mean annual groundwater recharge was derived by determining the product of estimated long-term mean annual runoff (the difference between precipitation and evapotranspiration) and the base-flow index (BFI). The BFI was calculated from daily streamflow data obtained from streamflow gauging stations in Taiwan. Mapping was achieved by using geographic information systems (GIS) and geostatistics. The presented approach does not require complex hydrogeologic modeling or detailed knowledge of soil characteristics, vegetation cover, or land-use practices. Contours of the resulting long-term mean annual P, BFI, runoff, groundwater recharge, and recharge rates fields are well matched with the topographical distribution of Taiwan, which extends from mountain range toward the alluvial plains of the island. The total groundwater recharge of Taiwan obtained by the employed method is about 18 billion tons per year.An erratum to this article can be found at     

Assessment of groundwater vulnerability to pollution in Dire Dawa,Ethiopia using DRASTIC

Assessment of groundwater vulnerability to pollution is an effective tool for the delineation of groundwater protection zones. DRASTIC approach was used to determine vulnerability zones in Dire Dawa groundwater basin, a semiarid region of Ethiopia. Maps of the seven DRASTIC parameters were prepared. GIS-ArcView was used for mapping and performing weighted-overlay analysis. The result of the analysis indicated that eastern part of the study area, in which Dire Dawa town is located, is highly vulnerable. A low aquifer vulnerability class was determined for the western portion of the study area as a result of greater groundwater depths, higher relative soil-clay content, and relatively low recharge rates for this area. The area between the two zones is of medium vulnerability. Observed nitrate concentrations in boreholes are in accordance with the vulnerability map. Some of the boreholes in the Sabian well field (Dire Dawa area) already deliver groundwater with nitrate levels significantly exceeding health standards set by the World Health Organization, while boreholes in the western part (low vulnerability zone) contain almost no nitrate. The result of this study is useful for risk assessments and for the development of effective groundwater management strategies for this region and others like it.     

Application of the water balance model J2000 to estimate groundwater recharge in a semi-arid environment: a case study in the Zarqa River catchment, NW-Jordan

Pollution and overexploitation of scarce groundwater resources is a serious problem in the Zarqa River catchment, Jordan. To estimate this resource’s potential, the amount and spatial distribution of groundwater recharge was calculated by applying the hydrological model J2000. The simulation period is composed of daily values gathered over a 30-year period (July 1977 to June 2007). The figure finally obtained for estimated groundwater recharge of the Zarqa River catchment is 105 × 10⁶ m³ per year (21 mm a^?1). This is 19 % higher than the value previously assumed to be correct by most Jordanian authorities. The average ratio of precipitation to groundwater recharge is 9.5 %. To directly validate modelled groundwater recharge, two independent methods were applied in spring catchments: (1) alteration of stable isotope signatures (δ¹⁸O, δ²H) between precipitation and groundwater and (2) the chloride mass balance method. Recharge rates determined by isotopic investigations are 25 % higher, and recharge rates determined by chloride mass balance are 9 % higher than the modelled results for the corresponding headwater catchments. This suggests a reasonably modelled safe yield estimation of groundwater resources.     

多年冻土区下垫面条件对坡面关键水循环过程的影响分析

为从整体上认识多年冻土流域水循环过程基本规律及其对下垫面条件变化响应,以长江源区风火山小流域为例,基于2016—2019年的水文气象要素的野外观测与计算,分析了坡面尺度上水分入渗、蒸散发、活动层内部水热条件以及冻结层上地下水等关键水循环过程的变化特征及其对下垫面条件变化的响应。研究结果表明:①风火山小流域生长季实际蒸散发的多年平均值为472.1±42.9 mm,实际蒸散发的气象影响因子排序为:净辐射(敏感系数S_Rn =1.22,相关系数R=0.93)＞气温(S_Ta =0.33,R=0.84)＞相对湿度(S_HR=0.32,R=0.46)＞风速(S_U =-0.25,R=-0.27),坡面尺度上实际蒸散发与植被覆盖度以及海拔高度正相关;②初始和稳定入渗速率均随坡位的升高而增大,对于稳定入渗率,初始融化期、完全融化期和初始冻结期,坡顶(1.07 mm/min、0.63 mm/min、0.88 mm/min)>坡中(0.29 mm/min、0.45 mm/min、0.21 mm/min)>坡底(0.11 mm/min、0.30 mm/min、0.10 mm/min),而植被覆盖度对入渗速率的影响在不同冻融阶段表现出差异;③阴、阳坡面和不同坡位冻融起始时间差异体现在土壤水分、地下水位变化上,解释了多年冻土流域产流区从阳坡发展至阴坡,从低海拔发展至高海拔的时空差异原因;④ 50 cm深度以下的活动层冻融状况显著影响着地下水位动态,地形、植被和土壤质地差异对热量传输的影响共同导致了坡面尺度上冻结层上地下水动态的空间差异性。     

Watershed-scale response of groundwater recharge to inter-annual and inter-decadal variability in precipitation (Alberta,Canada)

Groundwater recharge sets a constraint on aquifer water balance in the context of water management. Historical data on groundwater and other relevant hydrological processes can be used to understand the effects of climatic variability on recharge, but such data sets are rare. The climate of the Canadian prairies is characterized by large inter-annual and inter-decadal variability in precipitation, which provides opportunities to examine the response of groundwater recharge to changes in meteorological conditions. A decadal study was conducted in a small (250 km²) prairie watershed in Alberta, Canada. Relative magnitude of annual recharge, indicated by water-level rise, was significantly correlated with a combination of growing-season precipitation and snowmelt runoff, which drives depression-focussed infiltration of meltwater. Annual precipitation was greater than vapour flux at an experimental site in some years and smaller in other years. On average precipitation minus vapour flux was 10 mm y^?1, which was comparable to the magnitude of watershed-scale groundwater recharge estimated from creek baseflow. Average baseflow showed a distinct shift from a low value (4 mm y^?1) in 1982–1995 to a high value (15 mm y^?1) in 2003–2013, indicating the sensitivity of groundwater recharge to a decadal-scale variability of meteorological conditions.     

Groundwater flood of a river terrace in southwest Wisconsin,USA

Intense rainstorms in 2008 resulted in wide-spread flooding across the Midwestern United States. In Wisconsin, floodwater inundated a 17.7-km² area on an outwash terrace, 7.5 m above the mapped floodplain of the Wisconsin River. Surface-water runoff initiated the flooding, but results of field investigation and modeling indicate that rapid water-table rise and groundwater inundation caused the long-lasting flood far from the riparian floodplain. Local geologic and geomorphic features of the landscape lead to spatial variability in runoff and recharge to the unconfined sand and gravel aquifer, and regional hydrogeologic conditions increased groundwater discharge from the deep bedrock aquifer to the river valley. Although reports of extreme cases of groundwater flooding are uncommon, this occurrence had significant economic and social costs. Local, state and federal officials required hydrologic analysis to support emergency management and long-term flood mitigation strategies. Rapid, sustained water-table rise and the resultant flooding of this high-permeability aquifer illustrate a significant aspect of groundwater system response to an extreme precipitation event. Comprehensive land-use planning should encompass the potential for water-table rise and groundwater flooding in a variety of hydrogeologic settings, as future changes in climate may impact recharge and the water-table elevation.     

Estimation of groundwater recharge in arid, data scarce regions; an approach as applied in the El Hawashyia basin and Ghazala sub-basin (Gulf of Suez, Egypt)

In this study, an approach for runoff and recharge estimations that can be applied in arid regions which suffer from lack of data is presented. Estimating groundwater recharge in arid regions is an extremely important but difficult task, the main reason is the scarcity of data in arid regions. This is true for the Eastern Egyptian Desert where groundwater is used for irrigation purposes in agricultural reclamation along the Red Sea coast line. As a result of the scarcity of hydrologic information, the relation between rainfall and runoff was calculated depending on the paleo-flood hydrology information. Two models were used to calculate the rainfall–runoff relationships for El Hawashyia basin and Ghazala sub-basin. Two computer programs known as Gerinne (meaning channel in German) and SMADA6 (Stormwater Management and Design Aid, version 6) were conjunctively used for this purpose. As a result of the model applied to El Hawashyia basin, a rainfall event of a total of 18.3 mm with duration 3 h at the station of Hurghada, which has an exceedance probability of 5–10 %, produces a discharge volume of 10.2 × 10⁶ m³ at the delta, outlet of the basin, as 4.7 mm of the rainfall infiltrates (recharge). For the Ghazala sub-basin, the model yields a runoff volume of 3.16 × 10⁶ m³ transferred from a total rainfall of 25 mm over a period of 3 h, as 3.2 mm of it was lost as infiltration.     

The effects of a dry sand layer on groundwater recharge in extremely arid areas: field study in the western Hexi Corridor of northwestern China

Evaporation capacity is an important factor that cannot be ignored when judging whether extreme precipitation events will produce groundwater recharge. The evaporation layer’s role in groundwater recharge was evaluated using a lysimeter simulation experiment in the desert area of Dunhuang, in the western part of the Hexi Corridor in northwestern China’s Gansu Province. The annual precipitation in the study area is extremely low, averaging 38.87 mm during the 60-year study period, and daily pan evaporation amounts to 2,486 mm. Three simulated precipitation regimes (normal, 10 mm; ordinary annual maximum, 21 mm; and extreme, 31 mm) were used in the lysimeter simulation to allow monitoring of water movement and weighing to detect evaporative losses. The differences in soil-water content to a depth of 50 cm in the soil profile significantly affected rainfall infiltration during the initial stages of rainfall events. It was found that the presence of a dry 50-cm-deep sand layer was the key factor for “potential recharge” after the three rainfall events. Daily precipitation events less than 20 mm did not produce groundwater recharge because of the barrier effect created by the dry sand. Infiltration totaled 0.68 mm and penetrated to a depth below 50 cm with 31 mm of rainfall, representing potential recharge equivalent to 1.7 % of the rainfall. This suggests that only extreme precipitation events offer the possibility of recharge of groundwater in this extremely arid area.     

Delineation of potential recharge area using a hybrid model,case of Djelfa Hadjia watershed

Groundwater potential map is important for environmental assessment and water resources management. In this work, a groundwater recharge potential map was established for the watershed of Oued Djelfa Hadjia in Algeria, based on new multiparameters hybrid model. The model has hydroclimatic parameters, geological settings, slope factor, and stream network density factor as inputs. The groundwater recharge estimated by the model range from 0.71 to 14 mm. The model allows delineation of potential area of recharge. The total water abstraction in Djelfa city is about of 14 hm³; however, the calculated groundwater recharge is about 3 mm/year (min 0.71 mm and max 14 mm), which correspond to an average recharge volume of 3.9 hm³ which mean that the aquifer is under over exploitation.     

Groundwater recharge rates for regional groundwater modelling: a case study using GROWA in the Lower Rhine lignite mining area, Germany

Groundwater recharge rates calculated with the GROWA model have been applied as the recharge boundary condition for the regional groundwater model Rurscholle. This model simulates groundwater dynamics in the Pleistocene aquifers of the Lower Rhine lignite mining area (Germany). GROWA uses an area-differentiated approach to calculate recharge rates depending on runoff-relevant site characteristics, which are represented by a set of baseflow indices. The regional accuracy of the coupled groundwater and GROWA models has been checked using groundwater hydrographs as validation criteria. The results suggest that the current (unadjusted) version of GROWA underestimates the regional groundwater recharge rate by 10–20 mm/yr. The comparative analysis identified areas where recharge calculations could be improved by adjusting the baseflow indices for areas where runoff is dominated by slope, low water-logging and a low degree of sealing. Using the adjusted set of baseflow indices, the mean groundwater recharge rate of the Rurscholle region was modelled as approx. 170 mm/yr. This study highlights the benefit of using a coupled approach and being able to independently calibrate and validate groundwater recharge boundary conditions in regional groundwater models.     

Spatio-temporal variation of groundwater recharge in response to variability in precipitation, land use and soil in Yanqing Basin, Beijing, China

It is important to understand how groundwater recharge responds to precipitation variability in space and time, especially in those areas such as Yanqing Basin (China), where groundwater represents the sole water resource. A simple soil-water balance method is applied for spatio-temporal simulation of groundwater recharge in Yanqing Basin from 1981 to 2000. It was implemented on a monthly time step considering the effects of land use and soil texture. The area-average recharge associated with various land uses and soil textures was then compared with zonal analysis using a geographic information system (GIS). The main findings include: (1) the variation of groundwater recharge follows precipitation changes, either at yearly or seasonal intervals, (2) land use plays a more influential role in groundwater recharge than soil texture in this area, and (3) the water table quickly rises in response to recharge in the shallow parts of the aquifer, while there is a delay of 0.5–1.0?years where the groundwater level is at depth 4–10?m. The application demonstrates how spatio-temporal analysis can be utilized for groundwater-recharge estimation through distributed modeling and GIS.     

太行山山前平原地下水补给规律分析:以河北栾城为例

采用水均衡法与通量法对太行山山前平原(河北栾城)进行了地下水补给评价,结果表明:近20 a地下水入渗补给量为8～172mm/a(水均衡法),平均值49.3 mm/a,变化幅度大.其中2004年入渗补给量(63mm/a)比采用通量法计算结果(28.2mm/a)大一倍,表明土壤优先流对入渗补给作用显著.年入渗补给量总体上随有效降雨量增加、地下水位下降速率减少、作物实际蒸发蒸腾量减少而增加,但与有效灌溉量和灌溉量之和的相关性不显著.通过长序列地下水位与灌溉量等数据估算得到研究区侧向净补给量为125.2 mm/a,约占地下水总补给量的71.7％(175.4 mm/a),说明地下水主要以侧向补给为主.同时计算了入渗补给与侧向补给地下水的延迟时间,分别为60 a和1个月.     

Use of the subsurface thermal regime as a groundwater-flow tracer in the semi-arid western Nile Delta,Egypt

Temperature profiles from 25 boreholes were used to understand the spatial and vertical groundwater flow systems in the Western Nile Delta region of Egypt, as a case study of a semi-arid region. The study area is located between the Nile River and Wadi El Natrun. The recharge areas, which are located in the northeastern and the northwestern parts of the study area, have low subsurface temperatures. The discharge areas, which are located in the western (Wadi El Natrun) and southern (Moghra aquifer) parts of the study area, have higher subsurface temperatures. In the deeper zones, the effects of faults and the recharge area in the northeastern direction disappear at 80 m below sea level. For that depth, one main recharge and one main discharge area are recognized. The recharge area is located to the north in the Quaternary aquifer, and the discharge area is located to the south in the Miocene aquifer. Two-dimensional groundwater-flow and heat-transport models reveal that the sealing faults are the major factor disturbing the regional subsurface thermal regime in the study area. Besides the main recharge and discharge areas, the low permeability of the faults creates local discharge areas in its up-throw side and local recharge areas in its down-throw side. The estimated average linear groundwater velocity in the recharge area is 0.9 mm/day to the eastern direction and 14 mm/day to the northwest. The average linear groundwater discharge velocities range from 0.4 to 0.9 mm/day in the southern part.     

Environmental assessment of soil,groundwater, and surface water quality in the south of the Riyadh,Saudi Arabia

To understand the impact of the influence of treated wastewater, a study was undertaken in the downstream side of Wadi Hanifa in the southern part of Riyadh City in Saudi Arabia. Chemical elements from 17 groundwater samples, 9 surface water samples (treated wastewater), and 14 soil samples were analyzed. Water facies analyses showed that both groundwater and surface water belong to the SO₄-Cl class. However, the groundwater is characterized by high salinity (average value of 3547 mg/l), which may be result of the greater rock-water interaction and limited rainfall recharge. The NO₃ values are also high in the groundwater samples (average value of 40 mg/l) and are mainly attributed to the agricultural practices in the study area. The surface water samples (treated wastewater from the plant) shows an average salinity value of 1232 mg/l and is better suited for irrigation purposes. Heavy elements analyzed in the soil samples show high concentrations of all the elements except Mn and Ni as compared to their background concentration. Enrichment factor and index of geoaccumulation calculated from the soil samples indicate that they are severely enriched with respect to Cd and Se. The spatial distribution maps were prepared based on kriging interpolation technique to estimate the concentrations of the analyzed elements at unknown locations. The treated wastewater in the study area is better suited for agricultural and domestic activities as compared to groundwater.