Characterization of the groundwater flow system in the hillside and coastal aquifers of the Mersin-Tarsus region (Turkey)

In the region between Mersin and Tarsus cities, located along the Mediterranean Sea coast in southern Turkey, the demand for groundwater has increased dramatically as the available surface water supplies have already been developed. Fundamental information is required to characterize the existing groundwater system in this area in order to establish a sustainable groundwater-use policy. For this purpose, hydrochemical and environmental isotopic data were collected and integrated with available geological and hydrogeological information to develop a conceptual model of the system. Results, backed up mainly by depleted stable isotope composition and infinitesimal tritium content, suggest that most of the groundwater along the coastal zone is supplied by the neighboring mountain belt while local precipitation has also contributes to aquifer recharge. The validation of the conceptual perspective by a steady-state numerical groundwater flow model reveals that about 90% of the recharge to the aquifer system is supplied by the deep flow of karstic groundwater fed from the Taurus Mountains. Monitoring of changes in the recharge regime of the mountain sector seems to be critical in establishing future groundwater use policies.     

Hydrostratigraphic framework and hydrogeological behaviour of the Mancha Oriental System (SE Spain)

The Mancha Oriental System (MOS) has a surface area of 7,260 km², making it one of the largest carbonate aquifer systems in Spain. The system sustains about 1,000 km² of irrigated crops and supplies groundwater to 275,000 inhabitants. The economic transformation brought about by the development of extensive irrigated cropland has led to a water-balance disequilibrium of about 75 million m³/year. This input–output deficit has negative consequences in the quantity and quality of the available resources, in the river–aquifer relationship, and in the associated ecosystems as well. To understand the hydrogeological behaviour of the system, it is necessary to design a conceptual model. Further, the conceptualisation of a groundwater flow system is a requirement of the European Water Framework Directive for the characterisation of groundwater bodies. The robustness of the conceptual model depends heavily on the user capability of representing the real system. In this work, a multidisciplinary approach has been used to represent the three-dimensional geological framework and the groundwater flow conceptualisation of the MOS. Data management and three-dimensional visualisation have been carried out by means of geographical information system (GIS) tools and software for contouring and three-dimensional surface mapping.     

Intrinsic groundwater vulnerability determination at the aquifer scale: a methodology coupling travel time estimation and rating methods

Groundwater vulnerability has been subject of much research due to the valuable information it provides concerning groundwater protection and exploitation potential. Up to now, most groundwater vulnerability studies adopt subjective systems of rating the various factors and subsequently, their results are often ambiguous and contradicting. Within the present study a methodology for the estimation of intrinsic groundwater vulnerability at the aquifer scale is presented. The methodology is based on travel time estimation from specified sources of pollution to the aquifer. Besides the deterministic calculation of travel times, the methodology provides a rating system for each pollution source, based on its relative severity and the estimated threat that it poses to the aquifer. Therefore, it can be regarded as a hybrid method that couples the advantages provided by the physically based methods with those of the subjective rating systems. The methodology is applied to the Neon Sidirochorion aquifer, Northeastern Greece, an overexploited aquifer where river waters, sea waters and lake waters interact, causing groundwater quality deterioration to the aquifer. The results indicated that the proposed groundwater vulnerability assessment methodology is well capturing pollution related to saltwater intrusion and agricultural activities, while it is concluded that the conceptual model is significantly affecting the vulnerability assessment results and therefore has to be previously developed.     

Implementation and evaluation of multivariate analysis for groundwater hydrochemistry assessment in arid environments: a case study of Hajeb Elyoun–Jelma, Central Tunisia

In most arid zones, groundwater (GW) is the major source of domestic, agriculture, drinking, and industrial water. Accordingly, the monitoring of its quality by different techniques and tools is a vital issue. The purpose of this paper is the evaluation of the combination of principal components analysis (PCA) and geostatistics as a technique for (1) identifying the processes affecting the groundwater chemistry of the detrital unconfined Middle Miocene Aquifer (MMA) of the Hajeb elyoun Jelma (HJ) aquifer (Tunisia) and (2) mapping the controlling variables for groundwater quality. This work is based on a limited database recorded in 22 wells tapping the aquifer and unequally distributed in the field. The proposed approach is carried out in two steps. In the first step, the application of PCA revealed that rock–water interaction, agriculture irrigation and domestic effluents could explain 85 % of the observed variability of the chemical GW quality of the MMA. As a result, two new variables are defined: V1 (rock–water interaction influence) and V2 (irrigation and domestic effluent influence). In the second step, the spatial variability of these variables over the extent of the MMA is mapped by applying a kriging interpolation technique. The results of this study suggest that, while both natural and anthropogenic processes contribute to the GW quality of the MMA, natural impacts can be considered as the most important.     

Numerical groundwater-flow modeling to evaluate potential effects of pumping and recharge: implications for sustainable groundwater management in the Mahanadi delta region,India

Process-based groundwater models are useful to understand complex aquifer systems and make predictions about their response to hydrological changes. A conceptual model for evaluating responses to environmental changes is presented, considering the hydrogeologic framework, flow processes, aquifer hydraulic properties, boundary conditions, and sources and sinks of the groundwater system. Based on this conceptual model, a quasi-three-dimensional transient groundwater flow model was designed using MODFLOW to simulate the groundwater system of Mahanadi River delta, eastern India. The model was constructed in the context of an upper unconfined aquifer and lower confined aquifer, separated by an aquitard. Hydraulic heads of 13 shallow wells and 11 deep wells were used to calibrate transient groundwater conditions during 1997–2006, followed by validation (2007–2011). The aquifer and aquitard hydraulic properties were obtained by pumping tests and were calibrated along with the rainfall recharge. The statistical and graphical performance indicators suggested a reasonably good simulation of groundwater flow over the study area. Sensitivity analysis revealed that groundwater level is most sensitive to the hydraulic conductivities of both the aquifers, followed by vertical hydraulic conductivity of the confining layer. The calibrated model was then employed to explore groundwater-flow dynamics in response to changes in pumping and recharge conditions. The simulation results indicate that pumping has a substantial effect on the confined aquifer flow regime as compared to the unconfined aquifer. The results and insights from this study have important implications for other regional groundwater modeling studies, especially in multi-layered aquifer systems.     

地下水机器学习方法研究

随着我国地下水监测工作的高速发展,高频率高密度水位监测数据的出现催生了对其进行深入信息挖掘的需求。在传统地下水模型研究中,地下水水位监测值常位于模型构建过程的下游,当水位监测的时空密度逐渐增大时,新增信息无法有效传导至模型的规划阶段并指导概念模型的修订。文章提出了一种地下水系统补排边界的识别方法,在不建立地下水数值模型的前提下,以监测井空间位置为节点,按照德劳内原则建立三角网格。在此网格系统中,首先定义一个水力梯度变换函数gradF,以求取网格中任意位置的水力梯度;借鉴机器学习领域的优化算法,使用水力梯度场驱动含水层中随机分布质点的运行轨迹,并以此推断和识别区域内地下水补给和排泄边界。在环境地学计算平台EnviFusion-CGS中实现,并构建了详细工作流程。以山东省青岛市大沽河中下游含水层为示范区,对含水系统的补给区和排泄区的空间分布及其动态变化进行了分析,取得了良好效果。本研究为构建和修订已有含水层概念模型提供了新思路。     

Groundwater flow model, recharge estimation and sustainability in an arid region of Patagonia, Argentina

The Península Valdés, in northeastern Patagonia, Argentina, is characterised by its arid climate and the lack of perennial watercourses; thus, all economic activities depend on the groundwater resources. Water demand is mainly associated with tourism, which is centralised in Puerto Pirámides and supplied by a water desalination plant, and to sheep farming, supplied by the local aquifer. Due to the exponential growth of tourism, the government is planning to exploit groundwater and convey it by aqueduct to the abovementioned locality. The objectives of this study were to corroborate the conceptual geohydrological model, to develop a mathematical model to simulate the response of the aquifer to different scenarios, and to assess the incidence of water input into the system as a variable—a function that poses difficulties in the models for arid regions. The Visual Modflow 4.1 code was used, calibrating it in trial-and-error mode, changing the recharge and hydraulic conductivity parameters with different variants in the recharge zone and in the inclusion or exclusion of the evapotranspiration module. Results indicate the importance of the recharge analysis by treating rainfall at daily time steps. The adjusted model was exposed to four scenarios with variations in water input and in output by pumping. It can be concluded that under different input conditions, but with a controlled extraction, the system responds in a sustainable manner.     

Contribution to the hydrogeology of the Lower Cretaceous aquifer in east Central Sinai, Egypt

The present work deals with the groundwater aquifer of the Lower Cretaceous sandstone and its sustainable development in Sinai. The studied aquifer system is the most promising groundwater system in Sinai due to its wide extension, hug storage, and good quality. The objective of this paper aims to elucidate the hydrogeological characteristics of the Lower Cretaceous aquifer. The aquifer system occurs under confined conditions. The top surface of the Lower Cretaceous dips steeply towards the southwest direction with step faults. The average sand percent of the penetrated aquifer attains 54%. The main direction of groundwater flow is generally from southwest and locally is concentric to the center of study area related to the influence of the graben block. The aquifer has a hydraulic gradient generally reaches 0.0011 m/m and attains 0.0028 in central portion of study area. The aquifer parameters (effective porosity, transmissivity, and hydraulic conductivity) increase towards the northeast direction with increasing of the sand percentage. Durov diagram plot revealed that the groundwater has been a final stage evolution represented by a NaCl water type. The groundwater salinity increases towards the central of study area coinciding with groundwater flow. The groundwater salinity of the Lower Cretaceous aquifer is brackish water and varies from 2,510 to 5,256 ppm and unsuitable for drinking and domestic purposes.     

Exploring the use of alternative groundwater models to understand the hydrogeological flow processes in an alluvial context (Tiber River, Rome, Italy)

The first step of the application of a groundwater model to a specific site is the identification of the conceptual model. In this paper, a methodology to compare different conceptual models is proposed. Such a method consists in the parallel run of simple groundwater models, based on different conceptual models and in the evaluation of each model fitting. The case study proposed to illustrate this approach belongs to the alluvial series of the Tiber River (Rome, Italy).     

Evaluation of alternative conceptual models for groundwater modelling

This study evaluates the alternative conceptual models for groundwater modelling. A true model was created with a synthetic alluvial fan-plain hydrogeological framework. Various alternative conceptual models were evaluated for groundwater flow simulations. The first alternative model is a single aquifer layer model; the second alternative model is a 3-layer aquifer model; and the third model is a 5-layer model consisting of 3 aquifers separated by 2 aquitards. All models could fit very well to the observations with optimized values of hydraulic conductivities. However, the single aquifer layer model can only compute water balance components with good accuracy. The 3-layer aquifer model can be used for water balance computation and groundwater head simulation with small errors. The 5-layer model is capable of simulating water budget, groundwater head distribution and travel times with high accuracy. Multi-model analysis found only the 3rd alternative model superior.     

The Dammam aquifer in Bahrain – Hydrochemical characterization and alternatives for management of groundwater quality

 Over-abstraction of the Dammam aquifer, the principal aquifer in Bahrain, by the agricultural and domestic sectors, has led to its salinization by adjacent brackish and saline water bodies. A hydrochemical study identified the locations of the sources of aquifer salinization and delineated their areas of influence. The investigation indicates that the aquifer water quality is significantly modified as groundwater flows from the northwestern parts of Bahrain, where the aquifer receives its water by lateral underflow from eastern Saudi Arabia, to the southern and southeastern parts. Four types of salinization of the aquifer are identified: brackish-water up-flow from the underlying brackish-water zones in north-central, western, and eastern regions; seawater intrusion in the eastern region; intrusion of sabkha water in the southwestern region; and irrigation return flow in a local area in the western region. Four alternatives for the management of groundwater quality that are available to the water authorities in Bahrain are discussed and their priority areas are proposed, based on the type and extent of each salinization source, in addition to groundwater use in that area. The effectiveness of the proposed management options in controlling the degradation of water quality in the Dammam aquifer should be evaluated using simulation modeling. Received: June 1998 Revised: November 1998 Accepted: December 1998     

Fractured bedrock and saprolite hydrogeologic controls on groundwater/surface-water interaction: a conceptual model (Australia)

Hydrologic conceptual models of groundwater/surface-water interaction in a saprolite-fractured bedrock geological setting often assume that the saprolite zone is hydraulically more active than the deeper bedrock system and ignore the contribution of deeper groundwater from the fractured bedrock aquifer. A hydraulic, hydrochemical, and tracer-based study was conducted at Scott Creek, Mount Lofty Ranges, South Australia, to explore the importance of both the deeper fractured bedrock aquifer system and the shallow saprolite layer on groundwater/surface-water interaction. The results of this study suggest that groundwater flow in the deeper fractured bedrock zone is highly dynamic and is an important groundwater flow pathway along the hillslope. Deep groundwater is therefore a contributing component in streamflow generation at Scott Creek. The findings of this study suggest that hydrologic conceptual models, which treat the saprolite-fractured bedrock interface as a no-flow boundary and do not consider the deeper fractured bedrock in hydrologic analyses, may be overly simplistic and inherently misleading in some groundwater/surface-water interaction analyses. The results emphasise the need to understand the relative importance of subsurface flow activity in both of these shallow saprolite and deeper bedrock compartments as a basis for developing reliable conceptual hydrologic models of these systems.     

The typology of Irish hard-rock aquifers based on an integrated hydrogeological and geophysical approach

Groundwater flow in hard-rock aquifers is strongly controlled by the characteristics and distribution of structural heterogeneity. A methodology for catchment-scale characterisation is presented, based on the integration of complementary, multi-scale hydrogeological, geophysical and geological approaches. This was applied to three contrasting catchments underlain by metamorphic rocks in the northern parts of Ireland (Republic of Ireland and Northern Ireland, UK). Cross-validated surface and borehole geophysical investigations confirm the discontinuous overburden, lithological compartmentalisation of the bedrock and important spatial variations of the weathered bedrock profiles at macro-scale. Fracture analysis suggests that the recent (Alpine) tectonic fabric exerts strong control on the internal aquifer structure at meso-scale, which is likely to impact on the anisotropy of aquifer properties. The combination of the interpretation of depth-specific hydraulic-test data with the structural information provided by geophysical tests allows characterisation of the hydrodynamic properties of the identified aquifer units. Regionally, the distribution of hydraulic conductivities can be described by inverse power laws specific to the aquifer litho-type. Observed groundwater flow directions reflect this multi-scale structure. The proposed integrated approach applies widely available investigative tools to identify key dominant structures controlling groundwater flow, characterising the aquifer type for each catchment and resolving the spatial distribution of relevant aquifer units and associated hydrodynamic parameters.     

Groundwater mixing in a sand-island freshwater lens: density-dependent flow and stratigraphic controls

This paper focuses on a small back-barrier sand-island on the southeast coast of Queensland. The freshwater lens in the study area exhibits anomalously high short-range salinity gradients at shallow depths, which cannot be explained using a standard seawater intrusion model. The island groundwater system consists of two aquifers: a semiconfined aquifer hosting saline to hypersaline groundwater and an overlying unconfined freshwater aquifer. The deeper aquifer is semiconfined within an incised paleovalley, and groundwater flow is restricted to an east – west direction. Tidal response observations show that the tidal signal propagates far more rapidly and is of much higher magnitude in the semiconfined aquifer than the unconfined aquifer. The tidal wave-pulse amplitude is also subject to greater attenuation in the unconfined aquifer. A conceptual hydrogeological model illustrates how upwelling of hypersaline groundwater, induced by density-dependent flow and tidal pumping, has contaminated the shallow groundwater resource. Salinisation at shallow depths is restricted to an area proximal to the paleovalley aquifer. The spatial distribution of lithological heterogeneity is an initial limiting control on the movement of the upwelling saline plume. The extent of shallow groundwater contamination is also limited by the presence of a baroclinic field, resulting from lateral variations in fluid density. Hydrochemical signatures have been used to support the model hypothesis and link the salinisation of fresh groundwater with the semiconfined aquifer as opposed to the surrounding estuarine surface water. The geometry and thickness of the freshwater lens are further controlled by the presence of the largely impermeable bedrock paleosurface between 9 and 12 m depth. The combination of hypersaline groundwater and hydraulically restrictive lithology at shallow depths has produced excessive thinning of the freshwater lens, demonstrating that the application of a model such as the Dupuit – Ghyben – Herzberg relationship would grossly overestimate the available groundwater resource.     

Effects of different boundary conditions on the simulation of groundwater flow in a multi-layered coastal aquifer system (Taranto Gulf,southern Italy)

The evaluation of the accuracy or reasonableness of numerical models of groundwater flow is a complex task, due to the uncertainties in hydrodynamic properties and boundary conditions and the scarcity of good-quality field data. To assess model reliability, different calibration techniques are joined to evaluate the effects of different kinds of boundary conditions on the groundwater flow in a coastal multi-layered aquifer in southern Italy. In particular, both direct and indirect approaches for inverse modeling were joined through the calibration of one of the most uncertain parameters, namely the hydraulic conductivity of the karst deep hydrostratigraphic unit. The methodology proposed here, and applied to a real case study, confirmed that the selection of boundary conditions is among the most critical and difficult aspects of the characterization of a groundwater system for conceptual analysis or numerical simulation. The practical tests conducted in this study show that incorrect specification of boundary conditions prevents an acceptable match between the model response to the hydraulic stresses and the behavior of the natural system. Such effects have a negative impact on the applicability of numerical modeling to simulate groundwater dynamics in complex hydrogeological situations. This is particularly important for management of the aquifer system investigated in this work, which represents the only available freshwater resource of the study area, and is threatened by overexploitation and saltwater intrusion.     

Understanding groundwater processes by representing aquifer heterogeneity in the Maules Creek Catchment, Namoi Valley (New South Wales, Australia)

A FEFLOW three-dimensional (3D) groundwater model is developed to enhance the understanding of groundwater processes in the complex alluvial stratigraphy of Maules Creek Catchment (New South Wales, Australia). The aquifer vertical heterogeneity is replicated by indexing 204 lithological logs into units of high or low hydraulic conductivity, and by developing a 3D geological conceptual model with a vertical resolution based on the average lithological unit thickness for the region. The model mesh is populated with the indexed geology using nearest neighbour gridding. The calibrated model is successful in simulating the observed flow dynamics and in quantifying the important water-budget components. This indicates that the lateral groundwater flow from the mountainous region is the main inflow component of the system. Under natural conditions, the Namoi River acts as a sink of water, but groundwater abstraction increasingly removes a large amount of water each year causing dewatering of the system. The pumping condition affects the river–aquifer interaction by reversing the flow, from gaining to losing river conditions during the simulation period. The procedure is relevant for the development of groundwater models of heterogeneous systems in order to improve the understanding of the interplay between aquifer architecture and groundwater processes.     

Influence of fracture anisotropy and lithological heterogeneity on wellfield response in a fluvial sandstone aquifer of the Carboniferous Moncton Subbasin, Canada

Pump tests and geophysical logs acquired in a fluvial sandstone aquifer within the resource-rich Moncton Subbasin of New Brunswick, Canada, have been used to characterize fracture patterns and flow directions for purposes of developing a water-wellfield protection plan. Fracture patterns consist of three high-angle sets, and a low-angle set parallel to bedding. NW-trending high-angle fractures, predominantly in fluvial sandstone units, appear to be most important in controlling groundwater flow directions. This contrasts with an earlier regional hydrogeological study that attributed most flow to sub-horizontal bedding-plane fractures. Water levels monitored during a 72-h pump test revealed drawdown extension parallel to the NW-trending fracture set. Drawdown curves indicate that the aquifer is laterally constrained—likely reflecting differences in fracturing observed between the channelized sandstone and surrounding shale units. As a result, groundwater flow induced by pumping is influenced by both fracture anisotropy and by the heterogeneity of the fluvial depositional environment. Relationships observed between fracture patterns, regional geological structure and lithology provide a basis for extrapolating the conceptual model to other nearby areas in the region, where potential impacts of geological resource development on groundwater are attracting public concern.     

A proposal for an extended typology of groundwater habitats

Many governments now require the ecological condition of groundwater ecosystems to be considered when making policy decisions. However, groundwater habitats and communities occur at different spatial scales to those at which aquifers are managed, making their inclusion in policy decisions difficult. A system of classifying groundwater ecosystems is needed to aid this, yet such a typology is currently not available. Here, a hierarchical model for a typology of groundwater habitats is suggested, which considers three spatial scales: macro, meso, and local. On a macroscale, communities are influenced by biogeographic characteristics; at a mesoscale they are influenced by the hydrogeological type of the aquifer; and on local scale by the hydrological exchange with surface water and hence the oxygen and nutrient supply. Consequently, groundwater habitats are assumed to be defined by biogeographic particularities, the hydrogeological aquifer type and by localized hydrological exchange with surface water. By adopting a ‘groundwater habitat typology’, groundwater biologists may be able to work with hydrogeologists to assist in the sustainable management of groundwater by predicting groundwater biodiversity and the distribution of important or endangered species. The typology proposed here is a first step towards achieving this goal.     

Surface and subsurface conceptual model of an arid environment with respect to mid- and late Holocene climate changes

The water demand in arid regions is commonly covered by groundwater resources that date back to more humid periods of the Pleistocene and Holocene. Within the investigated arid part of SE Saudi-Arabia information about climate, groundwater levels, and pumping rates are only available for regions where groundwater extractions occur at present-day. For the prediction of the impact of long-term climate changes on groundwater resources an understanding of the hydrogeological and hydrological past and the development of the aquifers is necessary. Therefore, all available information about hydrology and hydrogeology for the past 10,000 years BP were collected and compiled to a conceptual model of the aquifer development on the Arabian Peninsula since the last Ice-Age. The climatic history was displayed by changes in precipitation, temperature and recharge during the mid-S and late Holocene. The hydrogeological development is described by groundwater ages, sea level fluctuations, movement of the coastline, and the development of sabkhas. The most sensitive parameter to describe the development of aquifer system is recharge. Present-day recharge was calculated with the hydrological model system HEC-HMS accounting for current precipitation, temperature, wind, soil types, and geomorphology. With respect to changes in precipitation and temperature over the past 10,000 years the temporal and spatial variability of groundwater recharge was calculated using empirical equations valid for semi-arid and arid settings. Further inflow into the groundwater system results from surface water infiltration in wadi beds, while natural outflow from the groundwater system occurs by discharge to the Gulf, evaporation from sabkhas, and spring discharge. Backward predictions can be verified by sedimentological observations of palaeo-river systems and lakes indicating that groundwater levels reached temporarily the surface under wetter climate conditions and ¹⁴C groundwater ages displaying groundwater residence times.     

Hydrogeological behaviour of some fault zones in a carbonate aquifer of Southern Italy: an experimentally based model

Complex flow circulation patterns are likely to be present in fault‐controlled groundwater flow systems, such as carbonate aquifers. Nevertheless, not much information is available for faults in carbonates, and their hydrogeological behaviour is often neglected in conceptual and numerical models. The understanding of this aspect of subsurface fluid flow has been improved in a carbonate aquifer, where hydrogeological investigations at site scale demonstrated the existence of fault zones that act as barriers. The hydraulic conductivity of the fault core is as low as that of siliciclastic rocks that represent the regional aquitard of the carbonate aquifer. Despite the lower permeability, the fault zones allow a significant groundwater flowthrough and a good interdependence of piezometric heads upgradient and downgradient of the faults. Because of this discontinuous heterogeneity, the aquifer looks like a basins‐in‐series system, where seasonal springs can be detected along some fault zones, as a function of groundwater level fluctuations.