A transdisciplinary analysis of water problems in the mountainous karst areas of Morocco

Water management and engineering in the karstic High Atlas of Morocco are difficult tasks under the prevailing geological, hydrogeological, geomorphological, vegetational and climatic conditions. It is important to be able to understand and predict the characteristics and availability of water for future water planning in the region under changing climatic and agricultural conditions. An interdisciplinary analysis of problems and adequate hydrological modelling tools developed by geologists, hydrologists and biologists are necessary. The karst areas of the High Atlas Mountains are characterised by impermeable triassic basalt underlying substantial subsurface reservoirs with high potential discharge rates. The karst groundwater aquifers are extensive but largely unknown in dimension, probably with a hierarchical network of groundwater flow paths. It is estimated that approximately 70% of the surface water is directly lost to groundwater. Steep landslide- and debris flow prone slopes exist next to coarse-grained, highly porous river beds. Infrequent, high intensity rainfall or snowmelt causes a particularly high flood risk to these karst areas. In addition, agriculture and land use changes have degraded the karst areas. The most important driving forces for degradation include permanent overgrazing even during droughts and the use of firewood by a continually growing population. Large scale degradation of vegetation has occurred in the oro-mediterranean (mountainous Mediterranean) zone, between 2600 and 3400 m which coincides with the most important zone for karstic groundwater creation. The combination of high amounts of groundwater flow and rapid surface flow due to sparse vegetation has increased the problems of flood flow.     

Hydrogeology of karst aquifer systems in SW Turkey and an assessment of water quality and contamination problems

The aim of this research is to determine the relationship between groundwater flow and water quality of different ground and surface water basins in the southwest Turkey. In addition, groundwater vulnerability is assessed taking into consideration groundwater flow and quality. The autochthonous Beydaglari limestone is the major karstic aquifer in the region. According to the groundwater level map of alluvium aquifers in the basins, groundwater discharge toward the carbonate aquifer is direct and indirect. The hydrogeological connection between ground and surface water basins occurs via the karstic aquifer located at the bottom of the alluvium bottom. In Egirdir lake, water also discharges in the karstic aquifer via karstic sinkholes at the western border of the lake. In the research area, general groundwater discharge is toward the Mediterranean Sea by means of autochthonous carbonate system, according to hydrogeological investigations, research of lineament and hydraulic conductivities. This result is supported by the locations of lineaments and shore springs discharging from the limestone. In addition, spreading of contaminants via karstic aquifer to great distance has been clearly identified.     

Numerical modelling of the saline interface in coastal karstic aquifers within a conceptual model uncertainty framework

Numerical modelling is increasingly used as a tool for improving management strategies in aquifers and to support the design of comprehensive projects considering natural and anthropogenic processes. Overall, numerical simulation in karstic aquifers poses a major scientific challenge due to the non-Darcian groundwater flow dynamics. In specific cases, the equivalent porous medium approach has shown acceptable results, particularly in poorly karstified aquifers with regional/subregional scales such as this case. The Yucatan coastal karstic aquifer (Mexico) has been defined as a complex regional heterogeneous system, partially confined, thus allowing the discussion of multiple conceptual models. In this research, a two-dimensional numerical model of flow and transport was implemented using SEAWAT for the NW Yucatan aquifer. Four likely conceptual models were audited, calibrated and verified using hydrogeological field data, to select the best one, considering their fit and complexity. The numerical model accuracy was evaluated using the root-mean-square error, Nash Sutcliffe efficiency and the Pearson coefficient. The Akaike information criterion and Bayesian information criterion were included for evaluating the complexity of the numerical models. In addition, the signal of tide propagation into the aquifer was assessed as a proxy to improve the numerical calibration process. Results show that the most complex numerical model has a better calibration than the simpler models, but the model accuracy is worse when compared to less complex numerical models in the verification exercise. This research offers enhancement in the knowledge of numerical modelling in heterogeneous coastal aquifers within a conceptual-model uncertainty setting.

     

黑河流域河道径流和人类活动对地下水动态的影响

文章概述了西北内陆干旱区黑河流域地下水系统的构成及第四纪含水层分布,在系统收集整理1986~2009年间黑河流域河道径流、地下水位、灌溉引水资料及补充调查成果的基础上,重点分析了径流变化和人类活动对流域地下水动态的影响。根据高程梯度、地貌类型和径流对平原地下水系统的补给作用,流域地下水动态类型可分为3个典型区：1）山前潜水含水层强入渗补给区,河道过水量控制着年地下水位变幅;2）冲积扇边缘细土平原承压含水层区,人类活动干预下的地表-地下水复杂交互过程使地下水位多年变化呈动态平衡状态;3）下游荒漠平原入渗生态补给区,政策性分水计划强烈影响地下水动态变化和空间分布。地表-地下水交互过程强烈影响地下水系统的动态平衡,且在不同区域呈现的时间尺度不同。中游绿洲农业区的人类活动通过政策性的水量再分配及土地利用空间格局变化影响地下水的循环路径及补径排条件。通过建立地表-地下水-生态水文过程耦合模型是量化动态水循环过程、合理管理地表-地下水资源的迫切需求,是未来干旱区水资源可持续利用研究的方向。     

Numerical simulation of groundwater flow and contamination transport in Nahavand plain aquifer,west of Iran

Numerical simulation of groundwater flow used for the estimation of hydraulic and hydrologic parameters which is an important tool for the management of aquifers. This study presents the results of a mathematical model developed for the simulation of groundwater flow in Nahavand plain aquifer in the southwest Hamadan province. For this purpose Groundwater Modeling Software (GMS) was used which supports the MODFLOW-2000 code. After gathering required data such as the hydrological, hydrogeological and topography maps, a 3D hydrogeological model of plain was constructed with borehole and surface elevation data. Then MODFLOW was used for simulation of flow. After initial simulation of the flow, the model was calibrated in steady state with trial-and-error and parameter estimation methods the observed head of groundwater table monitoring data of 1997. Results of calibration show that error between observed head and computed head is in allowable range. Also results of computed head with model show that groundwater flow is in the direction of the dominate slope (southeast to northwest). Finally MODPATH code which simulates advective transport of particles was used for estimation of flow path and source of contaminants.     

Hydrological and hydrogeological characterization of groundwater and river water in the North Hanoi industrial area,Vietnam

Degradation of groundwater quality by human activities is a widespread environmental problem in Vietnam. Groundwater there is a major source of water for domestic and industrial purposes. This paper reviews, compiles, and comprehensively analyzes spatiotemporal variations of hydrological and hydrogeological characteristics of shallow and deep groundwater aquifers in northern Hanoi industrial zones and in nearby Red River water. Groundwater level, electrical conductivity, and water temperature were measured in six monitoring wells, complemented by anion, cation, and stable isotope analyses of ground and surface water. The results show that the groundwater in both shallow and deep aquifers was fresh, but mainly calcium-bicarbonate type contaminants and human activities affect groundwater and surface water composition. With the goal of devising sustainable water use regulations, more research must be directed toward long-term monitoring of groundwater and surface water quality, as well as toward detailed investigation of the hydraulic characteristics of local aquifers in the study area.     

Development of a 3D groundwater model based on scarce data: the Wadi Kafrein catchment/Jordan

In this work, a 3D groundwater flow model integrating all important geological features of the hydrogeological system is developed to investigate hydrological processes in the Wadi Kafrein area of Jordan. A large amount of available geological and hydrological data is integrated to construct a 3D groundwater flow model for the Wadi Kafrein area. Using the newly developed mapping approach, the translation of the highly detailed geological formations to an unstructured finite element grids, can be accomplished with high precision. The existing data set for model calibration is scarce, which is a typical situation for many hydrogeological case studies. At first, the steady state calibration of the groundwater model is carried out based on the observation wells. Then, the time and space-dependent recharge from precipitation are applied at the top surface of the finite element model. The transient simulation is conducted during the period of 1996–2008 considering the abstraction rates of the production wells and discharge of the springs. The calculated water levels are close to the observed values. The difference is partly caused by return flows from irrigation and the groundwater inflow from the adjacent aquifers which are not taken into consideration so far. Since the Wadi Kafrein area is an important agricultural area in the semiarid region of the Lower Jordan Valley, the model developed in this study can be regarded as a useful tool for analyzing the hydrological processes and improving groundwater management practices elsewhere affected by similar geological and hydrogeological conditions.     

Application of MODFLOW for simulating groundwater flow in the Trifilia karst aquifer, Greece

Karst aquifers can have a complex flow as a result of the formation of large conduits from dissolution features. As a result, a three-dimensional finite-difference groundwater flow model (equivalent porous media) may not apply as the dual porosity nature of karst features and the effects of turbulent flow cannot be directly simulated. Statistical analysis of karst hydrographs of the Trifilia aquifer in Greece showed the existence of a slightly karstified mass with high primary porosity that regulates the flow. An equivalent porous media model was developed to simulate the Trifilia karst aquifer using MODFLOW. Steady state and transient state calibration gave encouraging results for the equivalent porous media approach, which does not consider pipe flow or turbulence. Detailed hydrogeological research conducted in the area helped define the aquifer hydraulic conductivity zones and extent; and flux to/from the aquifer. Only hydraulic conductivity and specific yield were adjusted during calibration, as the flux to/from the system was considered known and applied as boundary conditions. Small mean absolute and RMS piezometric head error of the model under both steady and transient state conditions were achieved.     

Management of karst aquifers in Serbia for water supply

Karst aquifers represent one of the main sources of water supply in Serbia. Wide distributions of karstic areas, abundant reserves, and excellent quality of karst groundwater have been the reasons for its extensive use in water supply systems throughout the country. In Serbia, 70 karstic sources have been tapped for centralized water supply, with the assessed minimal capacity exceeding 4.5 m³/s. Most of the large cities of Eastern and Western Serbia currently use karst groundwater for water supply; however, due to unstable flow regime when only natural springflow is tapped, numerous problems arise during the recession period (summer–autumn). During the last two decades, after favorable conditions had been confirmed through hydrogeological survey and feasibility studies, several newly constructed systems for artificial control of karst aquifers (mostly in Eastern Serbia) resulted in significantly improved water supply.     

A multi-analysis remote-sensing approach for mapping groundwater resources in the karstic Meo Vac Valley, Vietnam

Remote sensing data can be integrated with analyses of topography, structural geology, hydrogeology and geophysics. The integration gives premises for the delineation of zones of potential groundwater resources in strongly fractured and karstified deep aquifers in the uplifted Meo Vac Highland, northern Vietnam. Remote sensing analysis outlines geological faults with hydrogeological significance. These faults are combined with a derived lineament density map, interpreted analysis of surface flow direction and existing hydrogeological data, resulting in indications of groundwater flow direction. An analysis of the SPOT 5 band ratio 4/1, together with indications of surface-flow direction in low terrains, results in a determination of underground cavern passages. The delineated zones of potential groundwater resources are verified by detailed hydrogeological field surveys and geophysical measurements. Remote sensing analysis is shown to effectively contribute to the investigation of groundwater resources for a hydrogeologically complex area.     

Hydrogeological characteristics of some deep siphonal springs in Serbia and Montenegro karst

In terms of hydrogeological, engineering-geological, and hydrotechnical tapping in karst in relation to ground waters, karst channels, springs and ponors, speleodiving is the only research method which enables direct observation, studying and exact geological mapping of karst channels and caverns. Data collected during speleodiving research contribute considerably to the analysis of karst evolution process in the given region, which is very important in evaluating the depth of karstification and determining the main direction of the groundwater flow. In the past 30 years in Serbia and Montenegro, speleodivers have investigated over 40 siphonal springs, active cave channels and ponors, of which more than 20 are proof of deep siphonal circulation in karstic aquifers. The karstic springs are the most interesting phenomenon from a hydrogeological view point, and their investigations need particular attention. Most of significant karstic springs are on the rims of erosion basins—perimeters of karst poljes, river valleys, sea coasts and contact areas between karst aquifers and hydrogeological barriers. General characteristics of the spring regime are the direct correlation between precipitation and spring discharge. Moreover, the hydrogeological regime of these springs also depends on the size of the catchment area, karstic aquifer retardation capacity, total porosity, as well as lithological and structural characteristics.     

Compilation simulation of surface water and groundwater resources using the SWAT-MODFLOW model for a karstic basin in Iran

Modeling of karstic basins can provide a better understanding of the interactions between surface water and groundwater, a more accurate estimation of infiltrated water amount, and a more reliable water balance calculation. In this study, the hydrological simulation of a karstic basin in a semiarid region in Iran was performed in three different stages. In the first stage, the original SWAT model was used to simulate surface-water flow. Then, the SWAT-MODFLOW conjunctive model was implemented according to the groundwater characteristics of the study area. Finally, due to the karstic characteristics of the region and using the CrackFlow (CF) package, the SWAT-MODFLOW-CF conjunctive model was developed to improve the simulation results. The coefficient of determination (R²) and the Nash-Sutcliffe efficiency coefficient (NSE) as error evaluation criteria were calculated for the models, and their average values were 0.63 and 0.57 for SWAT, 0.68 and 0.61 for SWAT-MODFLOW, 0.73 and 0.7 for SWAT-MODFLOW-CF, respectively. Moreover, the mean absolute error (MAE) and root mean squared error (RMSE) of the calibration for groundwater simulation using the SWAT-MODFLOW model were 1.23 and 1.77 m, respectively. These values were 1.01 and 1.33 m after the calibration of the SWAT-MODFLOW-CF model. After modifying the CF code and keeping the seams and cracks open in both dry and wet conditions, the amount of infiltrated water increased and the aquifer water level rose. Therefore, the SWAT-MODFLOW-CF conjunctive model can be proposed for use in karstic areas containing a considerable amount of both surface water and groundwater resources.

     

The hydrochemical identification of groundwater flowing to the Bet She’an-Harod multiaquifer system (Lower Jordan Valley) by rare earth elements,yttrium, stable isotopes (H,O) and Tritium

The Bet She’an and Harod Valleys in Israel are regional recipients and mixing zones for groundwater draining from a multiple aquifer system, which includes carbonate and basalt aquifers and deep-seated pressurized brines. The aquifers drain through two types of outlets, distinct and mixed. The latter type is mainly conditioned by the occurrence of fault-blocks related to the Jordan Rift system, which act as connecting media between the aquifers and facilitate interaquifer flow. Conjoint application of rare earth element distribution and water isotopes enables detection of the local areas replenishment by rainfall infiltration and, in connection with the position of wells or springs, the identification of groundwater flow paths. Once stationary equilibria are established changes of REY composition between REY in groundwater and their surface adsorption, are negligible. In areas with little soil coverage and vegetation even recharge over young Tertiary and diagenetic Cretaceous limestones is distinguishable by their REY distribution patterns. Groundwater recharged over Tertiary limestones show higher REY abundance and more significant Ce anomalies than those derived from the Cretaceous limestones. Weathering of alkali olivine basalts leads to REY patterns in groundwater depleted in the middle REE. The improved knowledge of the hydrological systems is thought to be useful for regional hydrogeological modeling and for designing rational water management schemes.     

Physically based modelling of shallow landslide sediment yield at a catchment scale

 A shallow landslide erosion and sediment yield component, applicable at the basin scale, has been incorporated into the physically based, spatially distributed, hydrological and sediment transport modelling system, SHETRAN. The component determines when and where landslides occur in a basin in response to time-varying rainfall and snowmelt, the volume of material eroded and released for onward transport, and the impact on basin sediment yield. Derived relationships are used to link the SHETRAN grid resolution (up to 1 km), at which the basin hydrology and final sediment yield is modelled, to a subgrid resolution (typically around 10–100 m) at which landslide occurrence and erosion is modelled. The subgrid discretization, landslide susceptibility and potential landslide impact are determined in advance using a geographic information system (GIS), with SHETRAN then providing information on temporal variation in the factors controlling landsliding. The ability to simulate landslide sediment yield is demonstrated by a hypothetical application based on a catchment in Scotland. Received: 30 October 1996 · Accepted: 25 June 1997     

Monthly karstic spring flow forecasting using a sequential gaussian simulation technique

Analyzing groundwater hydrologic equations related to karstic aquifers and spring hydrograph simulation have become the focus of many researches. Having double or triple porosity structure, mixed flow nature, and varying conduit permeability have made these formations become complex heterogenic systems with great temporal and spatial hydrodynamic variability. In this paper, a conditional sequential gaussian simulation (SGS) is used to simulate monthly flow data of five karstic springs with different hydrogeological properties, located in Zagros Mountain Chain, in western Iran. To evaluate the performance of the SGS algorithm, the results are compared with those of an autoregressive integrated moving average (ARIMA) model. The results demonstrate the efficiency of the SGS model in simulation of monthly flows compared to the ARIMA model. They also show the suitability of this model for handling uncertainty associated with karstic spring flows through generation of several equally probable stochastic realizations.     

Water vulnerability assessment in karst environments: a new method of defining protection areas using a multi-attribute approach and GIS tools (EPIK method)

 Groundwater resources from karst aquifers play a major role in the water supply in karst areas in the world, such as in Switzerland. Defining groundwater protection zones in karst environment is frequently not founded on a solid hydrogeological basis. Protection zones are often inadequate and as a result they may be ineffective. In order to improve this situation, the Federal Office for Environment, Forests and Landscape with the Swiss National Hydrological and Geological Survey contracted the Centre of Hydrogeology of the Neuchatel University to develop a new groundwater protection-zones strategy in karst environment. This approach is based on the vulnerability mapping of the catchment areas of water supplies provided by springs or boreholes. Vulnerability is here defined as the intrinsic geological and hydrogeological characteristics which determine the sensitivity of groundwater to contamination by human activities. The EPIK method is a multi-attribute method for vulnerability mapping which takes into consideration the specific hydrogeological behaviour of karst aquifers. EPIK is based on a conceptual model of karst hydrological systems, which suggests considering four karst aquifer attributes: (1) Epikarst, (2) Protective cover, (3) Infiltration conditions and (4) Karst network development. Each of these four attributes is subdivided into classes which are mapped over the whole water catchment. The attributes and their classes are then weighted. Attribute maps are overlain in order to obtain a final vulnerability map. From the vulnerability map, the groundwater protection zones are defined precisely. This method was applied at several sites in Switzerland where agriculture contamination problems have frequently occurred. These applications resulted in recommend new boundaries for the karst water supplies protection-zones. Received: 27 October 1997 · Accepted: 4 July 1998     

Tracer tests and the structure of permeability in the Corallian limestone aquifer of northern England, UK

The Corallian limestone of northern England (UK) is widely exploited for water supplies and exhibits the karstic phenomena of sinking rivers, conduit development and groundwater velocities of several kilometres per day. To test a number of model-derived source protection zones and elucidate contaminant transport mechanisms in the aquifer, three tracer tests were conducted from a set of swallow-holes draining the River Derwent toward public water supply wells in the eastern part of the aquifer. Tracers used included: Enterobacter cloacae (bacteriophage), Photine C (optical brightener), sodium fluorescein (fluorescent dye) and sulphur hexafluoride (dissolved gas), the varying properties of which make them suitable analogues for different types of potential contaminant. Observed tracer transport times and arrival patterns indicate that tracer transport occurs through karstic channels embedded in a network of primary fissures which exert control over tracer concentrations once initial tracer plumes have passed. A dipole flow system is observed between the swallow-holes and the closest abstraction well, whilst previously modelled source protection zones do not accurately reflect either groundwater velocity or those areas of the aquifer supplying the wells. These findings imply that managing such aquifers for potential contamination should rely upon empirical tracer evidence for source-protection zone modelling.     

Using non-conservative tracers to characterise karstification processes in the Merinos-Colorado-Carrasco carbonate aquifer system (southern Spain)

The systematic sampling of the chemical composition of the groundwater from five karst springs (including an overflow spring) and one outflowing borehole have permitted to determine distinctive chemical changes in the waters that reflect the geochemical processes occurring in a carbonate aquifer system from southern Spain. The analysis of the dissolution parameters revealed that geochemical evolution of the karst waters basically depends on the availability of the minerals forming aquifer rocks and the residence time within the aquifers. In the three proposed scenarios in the aquifers, which include the preferential flow routines, the more important geochemical processes taking place during the groundwater flow from the recharge to the discharge zones are: CO₂ dissolution and exsolution (outgassing), calcite net dissolution, calcite and dolomite sequential dissolution, gypsum/anhydrite and halite dissolution, de-dolomitization and calcite precipitation. A detailed analysis of the hydrochemical data set, saturation indices of the minerals and partial pressure of CO₂ in the waters joined to the application of geochemical modelling methods allowed the elaboration of a hydrogeochemical model of the studied aquifers. The developed approach contributes to a better understanding of the karstification processes and the hydrogeological functioning of carbonate aquifers, the latter being a crucial aspect for the suitable management of the water resources.     

Modelling an inland Delta aquifer system to evolve pre-development management schemes: a case study in Upper Thamalakane River valley, Botswana, southern Africa

 Groundwater modelling studies have been found to be a potential tool in planning the pre-development management of groundwater resources in newly developing aquifer systems. One such study was attempted in Upper Thamalakane River valley, Okavango Delta, Botswana (southern Africa). There are three major aquifers separated by two aquitards in the valley portion. The top two aquifers are freshwater bearing zones and the bottom one is saline. The hydrological set-up of the basin is complex, as the groundwater flow directions are opposite in the upper-unconfined and in the lower-confined aquifers. A preliminary multilayer model was developed for this aquifer system by making use of only available data. The hydrodynamic behavior was then studied under two prediction scenarios to evolve appropriate management decisions for locating the well field (large diameter wells) in the upper aquifer by making use of induced river infiltration during the flood season. The aquifer response for variable river-flow conditions was studied and the induced river infiltration was quantified. Received: 27 August 1998 · Accepted: 8 March 1999