The impact of salt diapirs on the quality of carbonate karst waters,Bastak, Iran

The Gavbast karstic aquifer located in southern Iran is in direct contact with an exposed salt diapir. To assess the influence of the diapir on the quality of groundwater in the karstic aquifer, electrical conductivity, total dissolved solids, flow rate, temperature and major ion concentrations were measured at 57 sampling sites, including springs, surface waters and wells. A conceptual model of groundwater flow is proposed for the Gavbast karstic aquifer based on the geological setting, water budget, local base of erosion, and hydrochemistry of the sampling sites. The model suggests two subbasins in the Gavbast Anticline draining into two distinct discharging alluvial sections. Unexpectedly, groundwater discharging from the carbonate Gavbast aquifer is saline or brackish and water is of chloride type. The study indicates that the source of salinity of the Gavbast aquifers is infiltration of surface diapir-derived brine into the aquifer. The contribution of the diapir brine in the Gavbast karst aquifer is calculated about 4 L/s, using chloride mass balance. Construction of salt basins to evaporate brine discharging from the diapir springs is proposed to reduce the salinity of karst water. A row of strategically placed wells in the Gavbast karst aquifer would potentially exploit large volumes of fresh groundwater before it is contaminated by the salt. Such low-cost remediation should allow the agricultural exploitation of 40 km² of currently barren land.     

Implications of Kali-Hindon inter-stream aquifer water balance for groundwater management in western Uttar Pradesh

The Kali-Hindon inter-stream region extends over an area of 395 km² within the Ganga-Yamuna interfluve. It is a fertile tract for sugarcane cultivation. Groundwater is a primary resource for irrigation and industrial purposes. In recent years, over-exploitation has resulted in an adverse impact on the groundwater regime. In this study, an attempt has been made to calculate a water balance for the Kali-Hindon inter-stream region. Various inflows and outflows to and from the aquifer have been calculated. The recharge due to rainfall and other recharge parameters such as horizontal inflow, irrigation return flow and canal seepage were also evaluated. Groundwater withdrawals, evaporation from the water table, discharge from the aquifer to rivers and horizontal subsurface outflows were also estimated. The results show that total recharge into the system is 148.72 million cubic metres (Mcum), whereas the total discharge is 161.06 Mcum, leaving a deficit balance of −12.34 Mcum. Similarly, the groundwater balance was evaluated for the successive four years. The result shows that the groundwater balance is highly sensitive to variation in rainfall followed by draft through pumpage. The depths to water level are shallow in the canal-irrigated northern part of the basin and deeper in the southern part. The pre-monsoon and post-monsoon water levels range from 4.6 to 17.7 m below ground level (bgl) and from 3.5 to 16.5 m bgl respectively. It is concluded that the groundwater may be pumped in the canal-irrigated northern part, while withdrawals may be restricted to the southern portion of the basin, where intense abstraction has led to rapidly falling water table levels.     

Application of hydrochemical signatures to delineating portable groundwater resources in Ordos Basin, China

The Ordos Basin of China encompasses Shaanxi, Gansu, and Shanxi provinces, Ningxia and Inner Mongolia autonomous regions. It lacks significant surface water resources. Among the water-bearing formations, the Luohe formation, with an area of 1.32×10⁵ km², is the most prospective aquifer. Groundwater quality data collected at 211 boreholes drilled into the Luhe formation indicate a complex distribution of groundwater chemistry. The hydrochemical properties were used to study the recharge, runoff, and discharge conditions of the groundwater in Ordos Basin and to evaluate sustainable groundwater resources. In the northern part of the basin, the hydrochemistry types and the total dissolved solids (TDS) show a clear lateral transition from SEE to NWW, indicating that the groundwater gets recharge in the northwest region and discharges in the southeast region. In the southern part of the basin, maximum TDS occurs at the center of the Malian River valley, from which the TDS decreases radially. Therefore, the groundwater in the southern basin gets recharge from the southeast and southwest regions, and the Malian River valley is the discharge zone. As a result of this research, the areas with portable groundwater were delineated. They include most of the southeast region of the Sishili Ridge, east of the Ziwu Mountain, and the southwest corner of the basin. The TDS of the groundwater in these regions is less than 1 g/l, and the hydrochemistry type is either HCO₃ or HCO₃·SO₄.     

Conflicts between wetland conservation and groundwater exploitation: Two case histories in Spain

The problems in two Spanish national parks located on wetlands are analyzed. The hydrogeological and ecological characteristics of the two parks are somewhat different as are their respective degrees of deterioration. The Tablas de Daimiel National Park is located on the Central Plateau of Spain. It used to consist of a marshy area of about 20 km² around the confluence of two relatively small rivers. The area was marshy mainly because it was the natural discharge zone for a Tertiary aquifer system about 100 m thick extending over an area of some 5,000 km², composed of calcareous and detrital material of continental origin. The average annual recharge has been estimated at approximately 350 hm³/yr. Current groundwater withdrawal is around 450 hm³/yr, mainly used to irrigate a surface area of some 1,000 km². This overdevelopment has led to a continuous depletion of the regional water table and eventually to the drying out of the marshy area. Spontaneous combustion or fires caused by man have occurred in about 10 km² of the desiccated areas since the spring of 1986. The Doñana National Park is located on the estuary of the River Guadalquivir. The aquifer system of the Lower Guadalquivir estuary consists essentially of a permeable formation of unconsolidated Plioquaternary materials with an area of some 3,000 km². Under the marshy area (about 1,800 km²) the aquifer system is confined below low-permeability estuary deposits which can be over 100 m thick. Around most of the marshland the aquifer crops out and is recharged by rain. The Doñana National Park is over 700 km² in size, part of which is in the marshland and part in the recharge area where the aquifer is phreatic. In the 1970s Spain's largest irrigation project using groundwater, covering a surface area of 240 km², was planned in an area bordering on the national park. The initial project has been scaled down considerably as a result of protests by conservation groups. However, it seems likely that the water table depletion as a result of pumpage for irrigation could cause a large part of the ecotone situated at the contact-line between the marshland and the phreatic aquifer to disappear. This ecotone exists because it is a natural groundwater dischange area.     

Geochemical indicators of interbasin groundwater flow within the southern Rio Grande Valley, southwestern USA

Increased groundwater withdrawals for the growing population in the Rio Grande Valley and likely alteration of recharge to local aquifers with climate change necessitates an understanding of the groundwater connection between the Jornada del Muerto Basin and the adjoining and more heavily used aquifer in the Mesilla Basin. Separating the Jornada and Mesilla aquifers is a buried bedrock high from Tertiary intrusions. This bedrock high or divide restricts and/or retards interbasin flow from the Jornada aquifer into the Mesilla aquifer. The potentiometric surface of the southern Jornada aquifer near part of the bedrock high indicates a flow direction away from the divide because of a previously identified damming effect, but a groundwater outlet from the southern Jornada aquifer is necessary to balance inputs from the overall Jornada aquifer. Differences in geochemical constituents (major ions, δD, δ¹⁸O, δ³⁴S, and ⁸⁷Sr/⁸⁶Sr) indicate a deeper connection between the two aquifers through the Tertiary intrusions where Jornada water is geochemically altered because of a geothermal influence. Jornada groundwater likely is migrating through the bedrock high in deeper pathways formed by faults of the Jornada Fault Zone, in addition to Jornada water that overtops the bedrock high as previously identified as the only connection between the two aquifers. Increased groundwater withdrawals and lowering of the potentiometric surface of the Jornada aquifer may alter this contribution ratio with less overtopping of the bedrock high and a continued deeper flowpath contribution that could potentially increase salinity values in the Mesilla Basin near the divide.     

Modelling hydrostratigraphy and groundwater flow of a fractured and karst aquifer in a Mediterranean basin (Salento peninsula, southeastern Italy)

The control exerted by the hydrostratigraphic structure on aquifer recharge, groundwater flow and discharge along the coastal areas of a Mediterranean basin (Salento peninsula, about 5,000?km² wide, southern Italy) is assessed through the development and application of a groundwater flow model based on the reconstruction of the hydrostratigraphic architecture at the regional scale. The hydrostratigraphic model, obtained by processing surface and subsurface data, is applied to map the top of the main aquifer, which is hosted in the deep hydrostratigraphic unit corresponding to Cretaceous and Oligocene limestones with complex geometrical relationships with the sea. It is also used to estimate the aquifer recharge, which occurs by percolation through overlying younger sediments with low permeability. These data are completed with information about the soil use to estimate water abstraction for irrigation and with literature data to estimate the water abstraction for drinking and industrial purposes. The above-sketched conceptual model is the basis for a finite difference groundwater 2D pseudo-stationary flow model, which assumes the following fundamental approximations: the fractured and karst limestone hydrostratigraphic unit can be approximated, at the model scale, as a continuous medium for which the discrete Darcy??s law is valid; the transition zone between salt and fresh water is so small with respect to the grid spacing that the Ghyben?CHerzberg??s approximation for a sharp interface can be applied. Along the coastline different boundary conditions are assigned if the top of the limestone hydrostratigraphic unit lies either above the sea level (the aquifer has a free surface and fresh water is drained), or below the sea level (the aquifer is under pressure and the contact with sea occurs off-shore). The groundwater flow model correctly predicts the areas where the aquifer is fully saturated with salt water.     

Groundwater potential for irrigation in the East Oweinat area,Western Desert Egypt

The Nubia Sandstone aquifer system is one of the most extensive groundwater systems in North Africa, covering an area of about 2,000,000 km², including parts of Egypt, Libya, Sudan, and Chad. In the Western Desert of Egypt, the Nubian formation has a thermal gradient of 1.1–5°C 100 m^–1 with the exception of the East Oweinat area, located in the southern part of the Western Desert. This is the only part of this huge system where ground-water occurs under unconfmed conditions in an area where the Nubian sandstone crops out and is underlain by shallow basement rocks; in this area groundwater has no thermal characteristics. The aquifer system in the East Oweinat area attains a relatively high hydraulic conductivity. The direction of groundwater flow is generally northeastwards but is distorted at faults and fracture zones. Chemical analyses of groundwater in the area indicate a low salt content and suitability for irrigation purposes. As the estimated recharge to the area is low compared with the foreseen irrigation water requirement, the development of groundwater in the East Oweinat should be based on groundwater mining. Although the evaluation of the groundwater resources in East Oweinat has indicated that groundwater can be extracted at a rate of 4.7×10⁶ m³ d^–1, the long-term economics of extraction that can sustain large-scale development projects has to be assessed.     

Geochemical evolution of groundwater in southern Bengal Basin: The example of Rajarhat and adjoining areas,West Bengal,India

Detailed geochemical analysis of groundwater beneath 1223 km² area in southern Bengal Basin along with statistical analysis on the chemical data was attempted, to develop a better understanding of the geochemical processes that control the groundwater evolution in the deltaic aquifer of the region. Groundwater is categorized into three types: ‘excellent’, ‘good’ and ‘poor’ and seven hydrochemical facies are assigned to three broad types: ‘fresh’, ‘mixed’ and ‘brackish’ waters. The ‘fresh’ water type dominated with sodium indicates active flushing of the aquifer, whereas chloride-rich ‘brackish’ groundwater represents freshening of modified connate water. The ‘mixed’ type groundwater has possibly evolved due to hydraulic mixing of ‘fresh’ and ‘brackish’ waters. Enrichment of major ions in groundwater is due to weathering of feldspathic and ferro-magnesian minerals by percolating water. The groundwater of Rajarhat New Town (RNT) and adjacent areas in the north and southeast is contaminated with arsenic. Current-pumping may induce more arsenic to flow into the aquifers of RNT and Kolkata cities. Future large-scale pumping of groundwater beneath RNT can modify the hydrological system, which may transport arsenic and low quality water from adjacent aquifers to presently unpolluted aquifer.     

Hydrogeological framework and water balance studies in parts of Krishni–Yamuna interstream area, Western Uttar Pradesh, India

The Krishni–Yamuna interstream area is a micro-watershed in the Central Ganga Plain and a highly fertile track of Western Uttar Pradesh. The Sugarcane and wheat are the major crops of the area. Aquifers of Quaternary age form the major source of Irrigation and municipal water supplies. A detailed hydrogeological investigation was carried out in the study area with an objective to assess aquifer framework, groundwater quality and its resource potential. The hydrogeological cross section reveals occurrence of alternate layers of clay and sand. Aquifer broadly behaves as a single bodied aquifer down to the depth of 100 m bgl (metre below ground level) as the clay layers laterally pinch out. The depth to water in the area varies between 5 and 16.5 m bgl. The general groundwater flow direction is from NE to SW with few local variations. An attempt has been made to evaluate groundwater resources of the area. The water budget method focuses on the various components contributing to groundwater flow and groundwater storage changes. Changes in ground water storage can be attributed to rainfall recharge, irrigation return flow and ground water inflow to the basin minus baseflow (ground water discharge to streams or springs), evapotranspiration from ground water, pumping and ground water outflow from the basin. The recharge is obtained in the study area using Water table fluctuation and Tritium methods. The results of water balance study show that the total recharge in to the interstream region is of the order of 185.25 million m³ and discharge from the study area is of the order of 203.24 million m³, leaving a deficit balance of −17.99 million m³. Therefore, the present status of groundwater development in the present study area has acquired the declining trend. Thus, the hydrogeological analysis and water balance studies shows that the groundwater development has attained a critical state in the region.     

Estimation of groundwater recharge in a major sand and gravel aquifer in Ireland using multiple approaches

Groundwater recharge was investigated in the most extensive sand and gravel aquifer (area of approximately 200 km²) in the Republic of Ireland as part of a wider study seeking to derive recharge estimates using aquifer vulnerability mapping. The proportion of effective rainfall (total rainfall minus actual evapotranspiration) that leads to recharge is known as the recharge coefficient. The recharge investigation involved a variety of approaches, including soil moisture budgeting, well hydrograph analysis, numerical modelling and a catchment water balance. The adoption of multiple techniques provided insights on recharge and also on aquifer properties. Comparison of two soil moisture budgeting approaches (FAO Penman-Monteith with Penman-Grindley) showed how variations in the effective rainfall values from these methods influence groundwater levels simulated in a numerical groundwater model. The catchment water balance estimated the recharge coefficient to be between 81 and 85%, which is considered a reasonable range for this aquifer, where overland flow is rarely observed. The well hydrograph analysis, using a previous estimate of specific yield (0.13), gave recharge coefficients in the range of 40–80%, considered low for this aquifer: a revised specific yield of 0.19 resulted in a more reasonable range of recharge coefficients of between 70 and 100%.     

Mass balance simulation and principal components analysis applied to groundwater resources: Essaouira basin (Morocco)

The water resources of the Essaouira coastal basin (west of Marrakesh) are characteristic of a semi-arid climate and are severely impacted by the climate in terms of quantity and quality. Considering the importance of the Essaouira aquifer in the groundwater supply of a vast region (nearly 1,200 km²), a research study was conducted in order to better understand groundwater evolution in this aquifer system. It is a coastal aquifer located on the Atlantic coastline, southern Morocco, and salinization problems have been reported. Covering the Palaeozoic bedrock, the sedimentary series range from the Triassic to the Quaternary. Besides the possibility of seawater intrusion problems, the geological structures delineate a syncline bordered by the Tidzi diapir (of Triassic age) outcropping to the east and south. This is a recharge area for the aquifer, whereas the main groundwater flow direction is from SE to NW towards the Atlantic Ocean. In spite of the occurrence of calcareous and dolomitic levels, all waters in the Essaouira basin are of the Na–Cl-type. Based on a range of experimental methodologies, combined with PCA and geochemical modelling, it was possible to identify the mineralization processes occurring in the groundwater system, and the importance of the water–rock interaction in the water chemistry. Scenarios were tested using a simple mass balance model through the PHREEQC programme. The reaction path was assumed to be such that waters observed at shallow depths evolved to more mineralized waters. An important contribution of water–rock interaction in groundwater mineralization was found, corroborating the influence of preferential recharge from the Tidzi diapir in the water’s signature. Anthropogenic contamination was also identified and could lead to serious problems with groundwater degradation in the near future, in a country with scarce water resources.     

Modelling groundwater over-extraction in the southern Jordan Valley with scarce data

To deal with the challenge of groundwater over-extraction in arid and semi-arid environments, it is necessary to establish management strategies based on the knowledge of hydrogeological conditions, which can be difficult in places where hydrogeological data are dispersed, scarce or present potential misinformation. Groundwater levels in the southern Jordan Valley (Jordan) have decreased drastically in the last three decades, caused by over-extraction of groundwater for irrigation purposes. This study presents a local, two-dimensional and transient numerical groundwater model, using MODFLOW, to characterise the groundwater system and the water balance in the southern Jordan Valley. Furthermore, scenarios are simulated regarding hydrological conditions and management options, like extension of arable land and closure of illegal wells, influencing the projection of groundwater extraction. A limited dataset, literature values, field surveys, and the ‘crop water-requirement method’ are combined to determine boundary conditions, aquifer parameters, and sources and sinks. The model results show good agreement between predicted and observed values; groundwater-level contours agree with the conceptual model and expected flow direction, and, in terms of water balance, flow volumes are in accordance with literature values. Average annual water consumption for irrigation is estimated to be 29 million m³ and simulation results show that a reduction of groundwater pumping by 40% could recover groundwater heads, reducing the water taken from storage. This study presents an example of how to develop a local numerical groundwater model to support management strategies under the condition of data scarcity.     

Analysis of the exchange of groundwater and river water by using Radon-222 in the middle Heihe Basin of northwestern China

In arid regions of western China, water resources come from mountain watersheds and disappear in the desert plain. The exchange of surface water and groundwater takes place two or three times in a basin. It is essential to analyze the interaction of groundwater with surface water to use water resources effectively and predict the change in the water environment. The conventional method of analysis, however, measures only the flow of a stream and cannot determine groundwater seepage accurately. As the concentration of Radon-222 (²²²Rn) in groundwater is much higher than in surface water, the use of ²²²Rn was examined as an indicator for the analysis of the interaction between surface water and groundwater. Measurement of the ²²²Rn concentration in surface water was conducted to detect groundwater seepage into a stream in the middle Heihe Basin of northwestern China. Furthermore, the simultaneous groundwater flow into and out of a stream from the aquifers was quantified by solving the ²²²Rn mass balance equation, in which the losses of gas exchange and radioactive decay of ²²²Rn are considered. Meanwhile, river runoff was gauged to determine the exchange rates between surface water and groundwater. The result shows that ²²²Rn isotope can be used as a good environmental tracer with high sensitivity for the interaction between surface water and groundwater, especially in the fractured aquifer system, karst aquifer system and discharge basins.     

Modeling of limestone aquifer in the western part of the River Nile between Beni Suef and El Minia

A robust classification scheme for partitioning groundwater chemistry into homogeneous groups was an important tool for the characterization of Eocene limestone aquifer. The aquifer locally is composed of chalky limestone with thin clay intercalated (Samalut Fm.), the fissures, the joints, and the fractures are represented the conduits of the aquifer system. The flow patterns are conditioned by karstification processes which develop a conduit network and preserve low permeability microfractured blocks. The aquifer is mainly recharged by surrounding aquifers and agricultural wastewaters. The groundwater flows in the eastern part (due the Bahr Yossef and River Nile), which is a discharge area rather than a recharge. Twenty-eight groundwater samples was collected from the Eocene limestone aquifer and analyzed for isotopes, major, and trace elements. δD and δ¹⁸O concentrations ranged widely due to geology, infiltration of different surface waters, evaporation, and hydrogeology. The concentration of δD and δ¹⁸O isotopes is depleted in the northern zone of the northern part and western zone of the central and southern part of the study area. They are enriched due the eastern area of the central and southern part of the study area. δD vs. δ¹⁸O delineate the Pleistocene aquifer and has a strong influence than other waters on aquifer hydrogeochemistry. It is confirmed by the AquaChem outputs of the mixing proportions of different water types included in the aquifer system. Cl-δD and Cl-δ¹⁸O relationships indicate the role of evaporation especially due the eastern area of the central and southern part of the study area. This research tests the performance of the many available graphical and statistical methodologies used to classify water samples. R-mode clustering, correlation analysis, and principal component analysis were investigated. All the methods were discussed and compared as to their ability to cluster, ease of use, and ease of interpretation. Nearly most low-salinity waters are in equilibrium to supersaturate with respect to both carbonate minerals, while it is shifted to undersaturate with salinity. The inverse modeling findings clarify that the calcite, gypsum, and anhydrite dissolution increased due the northeastern area, middle zone, and southern corner of the northern, central, and southern part of the study area, respectively. The latter areas also were characterized by the lowest precipitation of the dolomite. Such areas are distinguished by much more enhancement for aquifer permeability and therefore transmissivity. The latter areas can be use as injection zone by fresh water. It can be a triple function; firstly, it recharges the saline Eocene limestone aquifer through the enhancement hydraulic conductivity and dilutes it. Secondly, it enhances much more the aquifer permeability and therefore the transmissivity. The Eocene limestone aquifer can be improved in quality and quantity by using such a model and exploits it as an alternative water resource with Quaternary aquifer and Nile water. Thirdly, it irrigates more areas to increase the income/capita. The dedolomitization represents the main hydrogeochemical process in the aquifer system. The geomedia (limestone, clay, marl, shale, and sand deposits) are in contact with water, therefore, the rock/water interaction, mixing, and ion exchange were estimated by the geochemical evolution of the groundwater systems.     

Sources of groundwater pumpage in a layered aquifer system in the Upper Gulf Coastal Plain, USA

Understanding groundwater-pumpage sources is essential for assessing impacts on water resources and sustainability. The objective of this study was to quantify pumping impacts and sources in dipping, unconfined/confined aquifers in the Gulf Coast (USA) using the Texas Carrizo-Wilcox aquifer. Potentiometric-surface and streamflow data and groundwater modeling were used to evaluate sources and impacts of pumpage. Estimated groundwater storage is much greater in the confined aquifer (2,200?km³) than in the unconfined aquifer (170?km³); however, feasibility of abstraction depends on pumpage impacts on the flow system. Simulated pre-development recharge (0.96?km³/yr) discharged through evapotranspiration (ET, ～37%), baseflow to streams (～57%), and to the confined aquifer (～6%). Transient simulations (1980–1999) show that pumpage changed three out of ten streams from gaining to losing in the semiarid south and reversed regional vertical flow gradients in ～40% of the entire aquifer area. Simulations of predictive pumpage to 2050 indicate continued storage depletion (41% from storage, 32% from local discharge, and 25% from regional discharge capture). It takes ～100?yrs to recover 40% of storage after pumpage ceases in the south. This study underscores the importance of considering capture mechanism and long-term system response in developing water-management strategies.     

神东煤炭基地1:50 000活鸡兔幅矿区水文地质图数据集

本数据集在1：50 000活鸡兔图幅内综合应用了遥感、地面调查、物探、钻探,以及水化学同位素等多种方法进行数据采集,共完成2期遥感数据解译800 km²,地面调查点377个,丰水期、枯水期地下水水位统测240点次,抽水试验10组,水体样品采集共计127组,水文地质钻探累计进尺600 m。本数据集细化了工作区内第四系孔隙含水层、白垩系洛河组孔隙-裂隙含水层富水等级,填补了侏罗系安定组、直罗组水文地质信息的空白;同时,加入了与采煤相关的次生环境地质要素,并首次将浅层含水层疏干区较为准确地勾画出来,还建立了采煤影响下的地下水的循环模式,突出了矿区水文地质调查的特色。本数据集为矿区水文地质调查工作的规划、部署与总结提供常态化支持,对保障国家能源安全、生态安全具有重要的现实意义。     

3D modeling and monitoring of karst system as a base for its evaluation and utilization: a case study from eastern Serbia

Beljanica Mountain in eastern Serbia is a part of the Carpathian Balkan arch (northern Alpine branch). It covers an area of about 300 km² and consists mostly of Jurassic and Cretaceous limestones. Numerous surface karst features, long caves and several large karstic springs located in Beljanica’s piedmont along the contact of karstic and non-karstic rocks are all indicators of an intense karstification. Currently, the large karstic water reserves of Beljanica Mountain are not properly utilized because of their distance from main consumers, the objection by national water managers that the springs lack a stable and sufficient discharge particularly during recession periods. Due to its unpolluted and high quality water, the area is a great prospect for future water supply, and provides an opportunity for artificial regulation and for the design and implementation of specific tapping structures. This paper includes an analysis of the created 3D ArcGIS model of karst interior and its correlation with historical and newly collected data of spring discharges and groundwater physico-chemical characteristics. The results of karst aquifer monitoring (both quantitative and qualitative) are linked with the results of extensive field geological and speleological survey of the upper non-saturated part of the karst (such as sinkholes, pits and caves) and with the investigation of the permanently saturated deeper part of the aquifer (including the diving methods). The model of karst interior is based on the data from the 69 caves, 15 sinks and 1,682 dolines (sinkholes) surveyed. The total length of the karst channels network, calculated using the GIS model and presented in a 3D environment, is 647 km. The catchment areas of five major springs that drain the areas are estimated to range from only 7 km² (Malo Vrelo Spring) to 124 km² (Vrelo Mlave Spring). The groundwater exploitable reserves of Beljanica karst aquifer are estimated to be over 4 m³/s. The waters are low mineralized, unpolluted and have a great potential for water supply.     

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.     

Groundwater exploitation and recharge rate estimation of a quaternary sand aquifer in Dar-es-Salaam area,Tanzania

Dar-es-Salaam City gets water supply from surface water and groundwater. The groundwater is used to supplement surface water supply and has increasingly become a major source of water supply in the city. The study area comprises three major parts: the central coastal plain with quaternary fluvial–deltaic sediments, the deltaic Mio-Pliocene clay-bound sands and gravels in the northwest and southeast and the Lower Miocene fluviatile sandstones of Pugu Hills in the west of the study area. The main objective of this study was to quantify the integrated water balance. The major source of renewable groundwater in the aquifer is rainfall. Hence, the average recharge of 256.2 mm/year (for the year 2006) to the aquifer was estimated using the balance method of Thornthwaite and Mather, which is equal to 99.4 hm³/year for the whole alluvial aquifer. This value was balanced with total groundwater abstraction of 8.59 hm³/year, baseflow to rivers of 75.7 hm³/year and discharge into the sea (15.11 hm³/year).     

Modeling of groundwater flow for Mujib aquifer, Jordan

Jordan is an arid country with very limited water resources. Groundwater is the main source for its water supply. Mujib aquifer is located in the central part of Jordan and is a major source of drinking water for Amman, Madaba and Karak cities. High abstraction rates from Mujib aquifer during the previous years lead to a major decline in water levels and deterioration in groundwater quality. Therefore, proper groundwater management of Mujib aquifer is necessary; and groundwater flow modeling is essential for proper management. For this purpose, Modflow was used to build a groundwater flow model to simulate the behavior of the flow system under different stresses. The model was calibrated for steady state condition by matching observed and simulated initial head counter lines. Drawdown data for the period 1985–1995 were used to calibrate the transient model by matching simulated drawdown with the observed one. Then, the transient model was validated by using drawdown data for the period 1996–2002. The results of the calibrated model showed that the horizontal hydraulic conductivity of the B2/A7 aquifer ranges between 0.001 and 40m/d. Calibrated specific yield ranges from 0.0001 to 0.15. The water balance for the steady state condition of Mujib aquifer indicated that the total annual direct recharge is 20.4 × 10⁶m³, the total annual inflow is 13.0 × 10⁶ m³, springs discharge is 15.3 × 10⁶ m³, and total annual outflow is 18.7 × 10⁶ m³. Different scenarios were considered to predict aquifer system response under different conditions. The results of the sensitivity analysis show that the model is highly sensitive to horizontal hydraulic conductivity and anisotropy and with lower level to the recharge rates. Also the model is sensitive to specific yield