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.     

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.     

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

Does irrigation with reclaimed water significantly pollute shallow aquifer with nitrate and salinity? An assay in a perurban area in North Tunisia

In Tunisia, reclaimed water is increasingly used for irrigation in order to mitigate water shortage. However, few studies have addressed the effect of such practice on the environment. Thus, we attempted in this paper to assess the impact of irrigation with reclaimed water on the nitrate content and salinity in the Nabeul shallow aquifer on the basis of satellite images and data from 53 sampled wells. Ordinary and indicator kriging were used to map the spatial variability of these groundwater chemical parameters and to locate the areas where water is suitable for drinking and irrigation. The results of this study have shown that reclaimed water is not an influential factor on groundwater contamination by nitrate and salinity. Cropping density is the main factor contributing to nitrate groundwater pollution, whereas salinity pollution is affected by a conjunction of factors such as seawater interaction and lithology. The predictive maps show that nitrate content in the groundwater ranges from 9.2 to 206 mg/L while the electric conductivity ranges from 2.2 to 8.5 dS/m. The high-nitrate concentration areas underlie sites with high annual crop density, whereas salinity decreases gradually moving away from the coastline. The probability maps reveal that almost the entire study area is unsuitable for drinking with regard to nitrate and salinity levels. Appropriate measures, such as the elaboration of codes of good agricultural practices and action programs, should be undertaken in order to prevent and/or remediate the contamination of the Nabeul shallow aquifer.     

Hydrogeochemical features of groundwater of semi-confined coastal aquifer in Amol–Ghaemshahr plain, Mazandaran Province, Northern Iran

Hydrogeochemical data of groundwater from the semi-confined aquifer of a coastal two-tier aquifer in Amol–Ghaemshahr plain, Mazandaran Province, Northern Iran reveal salinization of the fresh groundwater (FGW). The saline groundwater zone is oriented at an angle to both Caspian Sea coastline and groundwater flow direction and extends inland from the coastline for more than 40 km. Spearman’s rank correlation coefficient matrices, factor analysis data, and values of C ratio, chloro-alkaline indices, and Na⁺/Cl^? molar ratio indicate that the ionic load in the FGW is derived essentially from carbonic acid-aided weathering of carbonates and aluminosilicate minerals, relict connate saline water, and ion exchange reactions. Saline groundwater samples (SGWS) (n?=?20) can be classified into two groups. SGWS of group 1 (n?=?17) represent the saline groundwater zone below the Caspian Sea level, and salinization is attributed essentially to (1) lateral intrusion of Caspian seawater as a consequence of (a) excessive withdrawal of groundwater from closely spaced bore wells located in the eastern part of the coastal zone and (b) imbalance between recharge and discharge of the two-tier aquifer and (2) upconing of paleobrine (interfaced with FGW) along deep wells. SGWS of this group contain, on average, 7.9 % of saltwater, the composition of which is similar to that of Caspian seawater. SGWS of group 2 (n?=?3) belong to the saline groundwater zone encountered above the Caspian Sea level, and salinization of the groundwater representing these samples is attributed to irrigation return flow (n?=?2) and inflow of saline river water (n?=?1).     

Geostatistical Evaluation of Spatial Variation Related to Groundwater Quality Database: Case Study for Arak Plain Aquifer,Iran

Geostatistical methods are one of the advanced techniques to interpolate groundwater quality data. Geostatistical interpolation techniques employ both the mathematical and the statistical properties of the measured points. Compiling the data distribution on spatial and temporal domain is of crucial importance in order to evaluate its quality and safety. The main purpose of this paper is to assess groundwater quality of Arak plain, Iran, by an unbiased interpolated method so called Kriging. Therefore, seven quality variables of Arak plain aquifer including TDS, SAR, EC, Na⁺, TH, Cl^?, and SO₄ ^2? have been analyzed, studied, and interpreted statistically and geostatistically. Utilized data in this study were collected from 97 water well samples in Arak plain, in 2012. After normalizing data, variogram as a geostatistical tool for defining spatial regression was calculated and experimental variograms have been plotted by GS⁺ software, then the best theoretical model was fitted to each variogram based on minimum RSS error. Cross validation was used to determine the accuracy of the estimated data. The uncertainty of the method could be well assessed via this method since the method not only gave the average error (around 0 in this study) but also gave the standard deviation of the estimations. Therefore, more than 3800 points were estimated by ordinary Kriging algorithm in places which have not been sampled. Finally, estimation maps of groundwater quality were prepared and map of estimation variance, EV, has been presented to assess the quality of estimation in each estimated point. Results showed that the Kriging method is more accurate than the traditional interpolation algorithms not honoring the spatial properties of the database.     

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.     

Hydraulics Analysis for Groundwater Flow Through Permeable Reactive Barriers

Groundwater flow modelling is an important tool in simulating and predicting hydraulic behaviour of groundwater transporting in the domain consisted of groundwater flow zone (aquifer) and permeable reactive barriers (PRBs). The aquifer regime is modelled using the Darcy equation, whereas PRBs are simulated by the Brinkman equation. By combining the above equations, the present paper is devoted to analyse the effect of permeability of aquifer (K _aq ) and reactive media (K _p ), groundwater flow velocity, and barrier size as well as assembled style of barriers on key barrier design parameter, such as the hydraulic capture zone width. The simulations indicate that the capture zone width generally decreases with a decrease in K _p at a given K _aq ; however, there is relatively little decrease in capture zone width when the K _p of freshly installed reactive media is roughly 10 times higher than the K _aq . For a given barrier system, the absolute capture zone width increases with an increase in barrier size, whereas the increment of relative capture zone width decreases dramatically. Although the capture zone width remains relatively unchanged as groundwater flow velocity increases, it results in a shorter residence time in the reactive media that may require greater barrier thickness so that contaminant levels can be reduced to regulatory limits. In addition, the effect of parallel barriers interaction on the capture zone width should be incorporated into the design. The present paper is based on the continuous reactive barrier and homogeneous barrier system. Similar methods can be used to analyse the hydrodynamics of funnel-and-gate system and heterogeneous settings by the appropriate modification.     

Predicting saltwater intrusion into aquifers in vicinity of deserts using spatio-temporal kriging

The primary objective of this study was to provide a detailed framework to use the spatio-temporal kriging to model the spatio-temporal variations of salinity data and predict saltwater intrusion into freshwater aquifers in the vicinity of deserts. EC data, measured in extraction wells in the Mahvelat plain located in the Northeastern part of Iran, available from 2007 to 2013, were used to demonstrate the developed framework. The source of data was not a well-designed measurement network. Therefore, to homogenize the data, spatial analysis was used to find EC distribution in the area in each year of study. To conduct the spatial analysis, a guideline and a systematic process were developed to select an appropriate kriging method and optimize its parameters. This process can be applied to different variables. After spatial analysis of EC data for all the years of the analysis period using empirical Bayesian kriging (EBK) method with manually optimized parameters, spatio-temporal and corresponding variogram analysis was conducted using R software. This process was based on a separable product-sum model applied to the data from 2007 to 2012. The data of 2013 and the data available for the years 1999 and 2006 were used for evaluating the performance of the spatio-temporal model. The EC distribution maps, developed for different years until 2021, show a high level of EC in the north, south, and west of the study area and growing saltwater intrusion into the central freshwater aquifer. This result can be attributed to the over-exploitation of the aquifer and hydraulic head and gradient distribution in the area. The framework provided in this study for spatio-temporal analysis of unstructured EC data is useful for groundwater managers in making proper decisions.     

Decision support model for assessing aquifer pollution hazard and prioritizing groundwater resources management in the wet Pampa plain, Argentina

This paper gives an account of the implementation of a decision support system for assessing aquifer pollution hazard and prioritizing subwatersheds for groundwater resources management in the southeastern Pampa plain of Argentina. The use of this system is demonstrated with an example from Dulce Stream Basin (1,000 km² encompassing 27 subwatersheds), which has high level of agricultural activities and extensive available data regarding aquifer geology. In the logic model, aquifer pollution hazard is assessed as a function of two primary topics: groundwater and soil conditions. This logic model shows the state of each evaluated landscape with respect to aquifer pollution hazard based mainly on the parameters of the DRASTIC and GOD models. The decision model allows prioritizing subwatersheds for groundwater resources management according to three main criteria including farming activities, agrochemical application, and irrigation use. Stakeholder participation, through interviews, in combination with expert judgment was used to select and weight each criterion. The resulting subwatershed priority map, by combining the logic and decision models, allowed identifying five subwatersheds in the upper and middle basin as the main aquifer protection areas. The results reasonably fit the natural conditions of the basin, identifying those subwatersheds with shallow water depth, loam–loam silt texture soil media and pasture land cover in the middle basin, and others with intensive agricultural activity, coinciding with the natural recharge area to the aquifer system. Major difficulties and some recommendations of applying this methodology in real-world situations are discussed.     

Groundwater resource assessment, categories, and typologies: case study, Andhra Pradesh, India

In India, groundwater assessment units are classified as overexploited areas, critical areas, semi-critical, or safe areas based on the stage of groundwater development and long-term water level trends. Intuitively, in the safe units, wells are expected to function and have good yields. Besides, in the safe units, new wells are expected to be successful. Conversely, the expectation of a successful well or wells with good yields is much lesser in the overexploited units. However, when these expectations are not met in the field, doubts are raised about the quality of assessment and its usefulness, and there is outright distrust on the agencies assessing groundwater resource by the common man as well as on the planners, administrators, and the politicians. Therefore, there is a need to present the results in a way that does not create confusion. One of the methods is to combine the assessment results with aquifer characters using geographic information system (GIS); when this is done, a whole set of newer classes emerge, which can be mapped. These classes are termed as groundwater typologies in this study. Each typology has some characteristics or traits in common, which include basic aquifer character as well as the stage of groundwater development. Thus, a class may be safe, but if the aquifer is poor, then it is separated from a class that is safe and where the aquifer is good and so on. In Andhra Pradesh, which is taken as the case study for this purpose, eight main typologies emerged, and two of these main typologies were further divided into four subtypologies each. This new way of understanding the pattern of groundwater abstraction (using GIS) has a better visual impact. Groundwater typologies are found to be much more rational and useful in developing management strategies, rather than simple listing as overexploited areas, critical areas, semi-critical areas, and safe areas as is commonly done. The typologies so delineated indicate on the map (or table) that balanced usable groundwater is in between 5 and 6 bcm/a as against the estimated balance of 20.5 bcm/a, and it is largely in poor hard rock type of aquifers, which occupy about a third of the area of the state.     

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.     

Delineation of groundwater development potential zones in parts of marginal Ganga Alluvial Plain in South Bihar, Eastern India

A part of the Gangetic Alluvial Plain covering 2,228 km², in the state of Bihar, is studied for demarcating groundwater development potential zones. The area is mainly agrarian and experiencing intensive groundwater draft to the tune of 0.12 million cubic metre per square kilometres per year from the Quaternary marginal alluvial deposits, unconformably overlain northerly sloping Precambrian bedrock. Multiparametric data on groundwater comprising water level, hydraulic gradient (pre- and post-monsoon), aquifer thickness, permeability, suitability of groundwater for drinking and irrigation and groundwater resources vs. draft are spatially analysed and integrated on a Geographical Information System platform to generate thematic layers. By integrating these layers, three zones have been delineated based on groundwater development potential. It is inferred that about 48% of the area covering northern part has high development potential, while medium and low development potential category covers 41% of the area. Further increase in groundwater extraction is not recommended for an area of 173 km², affected by over-exploitation. The replenishable groundwater resource available for further extraction has been estimated. The development potential enhances towards north with increase in thickness of sediments. Local deviations are due to variation of—(1) cummulative thickness of aquifers, (2) deeper water level resulting from localised heavy groundwater extraction and (3) aquifer permeability.     

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.     

Modelling of lindane transport in groundwater of metropolitan city Vadodara,Gujarat, India

Migration pattern of organochloro pesticide lindane has been studied in groundwater of metropolitan city Vadodara. Groundwater flow was simulated using the groundwater flow model constructed up to a depth of 60 m considering a three-layer structure with grid size of 40?×?40?×?40 m³. The general groundwater flow direction is from northeast to south and southwest. The river Vishwamitri and river Jambua form natural hydrologic boundary. The constant head in the north and south end of the study area is taken as another boundary condition in the model. The hydraulic head distribution in the multilayer aquifer has been computed from the visual MODFLOW groundwater flow model. TDS has been computed though MT3D mass transport model starting with a background concentration of 500 mg/l and using a porosity value of 0.3. Simulated TDS values from the model matches well with the observed data. Model MT3D was run for lindane pesticide with a background concentration of 0.5?μg/l. The predictions of the mass transport model for next 50 years indicate that advancement of containment of plume size in the aquifer system both spatially and depth wise as a result of increasing level of pesticide in river Vishwamitri. The restoration of the aquifer system may take a very long time as seen from slow improvement in the groundwater quality from the predicted scenarios, thereby, indicating alarming situation of groundwater quality deterioration in different layers. It is recommended that all the industries operating in the region should install efficient effluent treatment plants to abate the pollution problem.     

Groundwater quality assessment in semi-arid regions using integrated approaches: the case of Grombalia aquifer (NE Tunisia)

As many arid and semi-arid regions in the Mediterranean Basin, the Grombalia coastal aquifer (NE Tunisia) is affected by severe groundwater exploitation and contamination. Therefore, quality assessments are becoming increasingly important as the long-term protection of water resources is at stake. Multidisciplinary investigations, like the one presented in this paper, are particularly effective in identifying the different origins of mineralization within an aquifer and investigating the impact of anthropogenic activities on groundwater quality. An integrated assessment, focused on the combined use of geostatistical, geochemical and isotopic (δ¹⁸O, δ²H and ³H) tools, was performed in the Grombalia aquifer between February and March 2014. The overall goal was to study the main processes controlling aquifer salinization, with special focus to nitrate contamination. Results indicate a persisting deterioration of water quality over the whole basin except the south-eastern zone juxtaposing the recharge area of the aquifer. Nitrate contents exceed the drinking water standard (50 mg/l) in 70% of groundwater samples, mainly due to the excessive use of fertilizers and urban activities. Stable isotope measurements showed the contribution of modern rainwater to aquifer recharge and proved the presence of evaporation contributing to the salinity increase. Tritium values of groundwater samples suggested two hypotheses: the existence of mixture between old and recent water or/and the existence of two recharge periods of the aquifer, pre- and post-nuclear weapons test. Principal component analysis confirmed the geochemical interpretation, highlighting that water-rock interaction evaporation effect and intensive anthropogenic activities constitute the main processes controlling the regional groundwater mineralization.     

Geostatistical Analysis of Spatial and Temporal Variations of Groundwater Level

Groundwater and water resources management plays a key role in conserving the sustainable conditions in arid and semi-arid regions. Applying management tools which can reveal the critical and hot conditions seems necessary due to some limitations such as labor and funding. In this study, spatial and temporal analysis of monthly groundwater level fluctuations of 39 piezometric wells monitored during 12 years was carried out. Geostatistics which has been introduced as a management and decision tool by many researchers has been applied to reveal the spatial and temporal structure of groundwater level fluctuation. Results showed that a strong spatial and temporal structure existed for groundwater level fluctuations due to very low nugget effects. Spatial analysis showed a strong structure of groundwater level drop across the study area and temporal analysis showed that groundwater level fluctuations have temporal structure. On average, the range of variograms for spatial and temporal analysis was about 9.7 km and 7.2 months, respectively. Ordinary and universal kriging methods with cross-validation were applied to assess the accuracy of the chosen variograms in estimation of the groundwater level drop and groundwater level fluctuations for spatial and temporal scales, respectively. Results of ordinary and universal krigings revealed that groundwater level drop and groundwater level fluctuations were underestimated by 3% and 6% for spatial and temporal analysis, respectively, which are very low and acceptable errors and support the unbiasedness hypothesis of kriging. Although, our results demonstrated that spatial structure was a little bit stronger than temporal structure, however, estimation of groundwater level drop and groundwater level fluctuations could be performed with low uncertainty in both space and time scales. Moreover, the results showed that kriging is a beneficial and capable tool for detecting those critical regions where need more attentions for sustainable use of groundwater. Regions in which were detected as critical areas need to be much more managed for using the current water resources efficiently. Conducting water harvesting systems especially in critical and hot areas in order to recharge the groundwater, and altering the current cropping pattern to another one that need less water requirement and applying modern irrigation techniques are highly recommended; otherwise, it is most likely that in a few years no more crop would be cultivated.     

Response of groundwater to climate variation: fluctuations of groundwater level and well yields in the Halacli aquifer (Cankiri, Turkey)

The aim of this study was to investigate the response of groundwater level and well yields in the Halacli aquifer to climate variations in Central Anatolia, Turkey. The Halacli aquifer is a typical aquifer due to its vulnerability to the climate variations. The aquifer is shallow and its recharge area is small. The waters from rains and snow melts can rapidly infiltrate down to the groundwater body because the vadose zone is thin and formed from coarse material. Therefore, the groundwater system responds to the short-term recharges by raising its level. Although any exploitation did not occur, the groundwater levels have declined from 1989 to 1997. However, the groundwater levels began rising when the exploitation started in the summer of 1998. After the year 2000, although the amount and duration of yearly exploitation was constant, fluctuations of water level continued. Fluctuation of groundwater levels and well yields bewilders the water users and imperils the sustainable water management in the study area and also in arid and semi-arid regions of Turkey. In order to overcome this problem, behavior of groundwater level and discharges of the wells must be recorded and the water users must be informed about the current conditions and the possible trend in the future of the system.     

Impacts of urbanization on groundwater hydrodynamics and hydrochemistry of the Toluca Valley aquifer (Mexico)

The Toluca Valley is located on the high plains of Mexico, where there are significant industrial zones and large populations. Water needs are almost exclusively met by groundwater, which has brought about intense exploitation of the aquifer and indication of some contamination. The present study investigates the effect of urbanization, related to industrialization of the region, on groundwater in the central portion of the Toluca Valley aquifer—a zone with high population density and where the largest industrial park is located. A general decline in the groundwater level has been found over the years, at a rate of as much as 2.5 m/year. The appearance of a large drawdown cone was identified, indicating changes in the direction of groundwater flow. Also identified was the presence of several ground fissures, the location of which coincided with the drawdown cone. In hydrochemical terms, the water type is sodium-magnesium bicarbonate and this characteristic has not changed over time, although it has been possible to detect the presence of larger quantities of sulfates (up to 117 mg/L) and nitrates (up to 47 mg/L) in recent years, likely associated with contamination from industrial and urban wastewater. Factor analysis made it possible to identify ions that would characterize natural processes involving the acquisition of salts (HCO₃ ^?, Na⁺, Mg²⁺, and Si), as well as anthropic activities (SO₄ ^2?, NO₃ ^?, Cl^?, Ca²⁺, and K⁺).