Integrated Geophysical Studies in the East-Indian Geothermal Province

Integrated geophysical surveys using vertical electrical sounding (VES), very low frequency (VLF) EM, radiation counting, total magnetic field and self-potential (SP) measurements are carried out to characterize the geothermal area around a hot spring in the Nayagarh district, Orissa, India that lies in the East Indian geothermal province. The study was performed to delineate the fracture pattern, contaminated groundwater movement and possible heating source. VES interpretations suggest a three- to four-layer structure in the area. Resistivity survey near the hot spring suggests that weathered and fractured formations constitute the main aquifer system and extend to 60 m depth. Current flow measured at various electrode separations normalized by the applied voltage suggests that fractures extend to a greater depth. Detailed VLF study shows that fractures extend beyond 70 m depth. VLF anomaly has also very good correlation with the total magnetic field measured along the same profiles. Study results suggest that a gridded pattern of VLF survey could map the underground conductive fracture zones that can identify the movement of contaminated groundwater flow. Therefore, precautionary measures can be taken to check further contamination by delineating subsurface conducting structures. Self potential (SP) measured over the hot spring does not show a large anomaly in favor of the presence of a sulphide mineral body. A small positive (5–15mV) SP anomaly is measured which may be streaming potential due to subsurface fluid flow. A high radiation is measured about four kilometers from the hot spring, suggesting possible radiogenic heating. However, the exact nature of the heating source and its depth is not known in the area. Deep resistivity followed by a magneto-telluric survey could reveal the deeper structures.     

Groundwater exploration using integrated geophysical techniques

The integrated approach to solving complicated geological, hydrological and environmental problems is now widely used in geophysics. Among all the geophysical methods, electrical and electromagnetic techniques are the most popular in groundwater exploration due to the close relationship between electrical conductivity and some hydrogeological properties of the aquifer (e.g. porosity, clay content, mineralization of the groundwater and degree of water saturation). Case histories presented here show that by proper combination of different techniques such as conventional direct current (DC) resistivity as well as the more advanced electromagnetic (EM) methods and the most recent nuclear magnetic resonance (NMR) tomography, the reliability of interpretation as compared to that typical for the individual methods can be significantly improved.     

Efficiency of joint use of MRS and VES to characterize coastal aquifer in Myanmar

The productivity and the water quality of coastal aquifers can be highly heterogeneous in a complex environment. The characterization of these aquifers can be improved by hydrogeological and complementary geophysical surveys. Such an integrated approach is developed in a non-consolidated coastal aquifer in Myanmar (previously named Burma).A preliminary hydrogeological survey is conducted to know better the targeted aquifers. Then, 25 sites are selected to characterize aquifers through borehole drillings and pumping tests implementation. In the same sites, magnetic resonance soundings (MRS) and vertical electrical soundings (VES) are carried out. Geophysical results are compared to hydrogeological data, and geophysical parameters are used to characterize aquifers using conversion equations. Finally, combining the analysis of technical and economical impacts of geophysics, a methodology is proposed to characterize non-consolidated coastal aquifers.Depth and thickness of saturated zone is determined by means of MRS in 68% of the sites (evaluated with 34 soundings). The average accuracy of confined storativity estimated with MRS is ± 6% (evaluated over 7 pumping tests) whereas the average accuracy of transmissivity estimation with MRS is ± 45% (evaluated using 15 pumping tests). To reduce uncertainty in VES interpretation, the aquifer geometry estimated with MRS is used as a fixed parameter in VES inversion. The accuracy of groundwater electrical conductivity evaluation from 15 VES is enough to estimate the risk of water salinity. In addition, the maximum depth of penetration of the MRS depends on the rocks' electrical resistivity and is between 20 and 80 m at the study area.     

Spatial delineation of riparian groundwater within alluvium deposit of mountainous region using Laplace equation

The interaction between surface and groundwater plays a key role in a riparian ecosystem while the size of riparian groundwater has not been typically incorporated into hydrological modelling systems. An extensive geophysical survey composed of 25 individual DC electrical resistivity profiles was conducted at the Blair–Wallis site in Wyoming. The observed resistivity images show a near‐surface aquifer interpreted as the saturated alluvium deposit along the channel, rather than the geological bedrock. Based on the electrical resistivity images, it can be inferred that only the near‐surface portion of the groundwater actively interacts with the stream flow in the mountainous and hilly watershed. This study attempted the spatial extrapolation of the measured riparian aquifer depths by means of fitting functions based on the surface topography. The analysis indicated that the boundary of the riparian aquifer well corresponds to the topographical inflexion point of the hill slope profile. It was also demonstrated that the extent of alluvium deposit, where the area of riparian aquifer is indicated, can be delineated using the slope and curvature maps in the geographic information system. Then, the parabolic and biharmonic functions were tested for the groundwater depth estimation using the developed alluvium deposit map. The proposed methodology was effective if geological diffusion processes by wind and water dominated the topography. The spatial map of the active aquifer will be useful in hydrological drought analysis because it is considered to be a main source of baseflow during dry seasons.     

Determination of aquifer geometry through geophysical methods: A case study from Quetta Valley, Pakistan

Due to increase in population and agricultural activities, the aquifer of Quetta Valley is under tremendous stress and the water table is declining at an increasing rate. This situation necessitates evaluation of the aquifer system, for which information about geometry of the aquifer is a prerequisite. However, there are no drilling-to-bedrock data available; therefore, electrical resistivity, seismic reflection and gravity methods were employed to determine geometry of the aquifer. Interpretation of vertical electrical soundings provided information about the depth-tobedrock at some specific points, whereas seismic reflection delineated bedrock topography along two lines. The depths to bedrock inferred from electrical resistivity and seismic reflection data were used as constraints in the modeling of gravity data. 2.75D gravity models were constructed along lines with a regular spacing. Map of depth-to-bedrock was prepared by contouring the depth given by the gravity models. Combination of these geophysical methods depicted the geometry of the aquifer. This example shows that in a similar geological setting proper integration of geophysical exploration methods can determine the aquifer geometry with an acceptable reliability and at an appropriate cost.     

Integration of remote sensing and electrical sounding data for hydrogeological exploration—a case study of Bakhar watershed,India / Intégration de données de télédétection et de sondages électriques pour l'exploration hydrogéologique—étude de cas du bassin versant de Bakhar,Inde

Abstract

A geographical information system (GIS) was used for the integration of hydrological data acquired using remote sensing and geoelectrical techniques to understand the groundwater condition of Bakhar watershed, Mirazpur District, UP, India. Indian remote sensing IRS-1D, LISS—III data were used to prepare a geomorphological and lineament map of the Bakhar watershed. Vertical electrical sounding (VES) was carried out in different geomorphic units, and ranges of electrical resistivity values were assigned to the different formations by calibrating electrical resistivity with borehole data. Based on these, a subsurface resistivity map and an aquifer thickness map were prepared. Several layers were superimposed using GIS techniques. Each theme was assigned a weight, depending on its influence on groundwater recharge. Each class or unit in the map was assigned a knowledge-based rank from one to four, depending on its significance in storage and transmittance of groundwater, and these were then multiplied by the layer weighting to produce a score. Based on these scores, the watershed was categorized into different groundwater potential zones. The results indicate that the eastern and northern parts of the study area have very good groundwater potential to meet the demands of water for irrigation and domestic purposes, whereas the southern region has poor groundwater potential zones. Such integrated analysis has not been attempted so far in this region for hydrogeological investigation.     

Seawater intrusion and fresh groundwater hydraulics in fjord delta aquifers inferred from ground penetrating radar and resistivity profiles—Sunndalsøra and Esebotn, western Norway

Geophysical surveys have been carried out in two fjord delta aquifers in western Norway. Geophysical methods comprised ground penetrating radar (GPR), shallow reflection seismic and resistivity measurements. The purpose of the investigations was to evaluate hydrogeological properties of the fjord delta aquifers with emphasis on the possibilities of abstracting saline groundwater for use in fish farming. At Sunndalsøra, reflection seismic was helpful in the mapping of the base of the aquifer. Resistivity profiles revealed both the spatial and seasonal variability of seawater intrusion. Penetration depths of ground penetrating radar (GPR-D) revealed the spatial variation of the boundary of the transition zone between saline and fresh groundwater. Maps of GPR-D can be used to indicate direction and magnitude of mean hydraulic gradients and thus to locate optimal sites for both fresh and saline groundwater abstraction.In fjord delta aquifers, both rivers and abandoned river channels constitute important groundwater divides. Optimal locations for saline groundwater abstraction wells are in areas separated from the rivers by abandoned river channels. In areas between the rivers and the abandoned river channels, groundwater abstraction will result in decreasing salinity of the water.     

2D and 3D resistivity inversion of Schlumberger vertical electrical soundings in Wadi El Natrun,Egypt: A case study

The Wadi El Natrun area is characterized by a very complicated geological and hydrogeological system. 45 vertical electrical soundings (Schlumberger array) were measured in the study area to elucidate the peculiarity of this unique regime, specifically the nature of waterless area. 2D and 3D resistivity inversion based on the finite element technique and regularization method were applied on the data set. 2D and 3D model resolution was investigated through the use of the Depth and Volume of Investigation Indexes. A very good matching was found between the zones of high resistivity, the waterless area, and the non-productive wells. The low resistivity zones (corresponding to Lower Pliocene clay) were also identified. The middle resistivity fresh water aquifer zones were recognized. Available results can assist in the aquifer management by selecting the most productive zone of groundwater.     

Inversion of resistivity in Magnetic Resonance Sounding

Magnetic Resonance Sounding (MRS, or Surface Nuclear Magnetic Resonance - SNMR) is used for groundwater exploration and aquifer characterization. Since this is an electromagnetic method, the excitation magnetic field depends on the resistivity of the subsurface. Therefore, the resistivity has to be taken into account in the inversion: either as a priori information or as an inversion parameter during the inversion process, as introduced in the presented paper. Studies with synthetic data show that water content and resistivity can be resolved for a low resistive aquifer even using only the amplitude of the MRS signal. However, the inversion result can be significantly improved using amplitude and phase of the MRS signal. The successful implementation of the inversion for field data shows that the resistivities derived from MRS are comparable to those from conventional geoelectric methods such as DC resistivity and transient electromagnetic. By having information about both the resistivity and the water content, MRS inversions give information about the quality of the water in the aquifer. This is of utmost interest in hydrogeological studies as this specific information cannot be determined solely by geoelectric measurements, due to the nonunique dependence of resistivity on water content and salinity.     

Using geophysical methods to define the attitude and extension of water-bearing strata in the Miocene sediments of the Pannonian Basin

In the area near the village of Jazak (southern part of Fruška Gora mountain, Serbia), hydrogeological investigations were carried out for the purpose of finding a water supply source to provide an adequate volume of water for a mineral water bottling plant. The first exploratory borehole (IBJf-1) penetrated a water-bearing layer of Miocene organogenic limestones. This aquifer has a thickness of about 30 m and a yield of only 2.2 l/s, which falls short of the required water volume (5 l/s).The objective of further exploration was to define the attitude and extension of the aquifer and thus select a more favourable site for a new exploratory borehole that would secure the required volume of water. For this purpose, geophysical exploration was carried out in 2003 through vertical electrical sounding (VES) and high-resolution 3D reflection seismic methods. The VES measurements enabled determination of aquifer depth and indicated that the water-bearing strata extend over the entire area studied. However, because of the equivalence problem, it was not possible to determine the thickness of the water-bearing stratum based solely on the VES data. Thus, the 3D seismic method was used in the second stage of investigation. A low-cost 3D seismic survey was carried out with fixed receiver lines, using a vibrator as the source of the seismic waves.From the 3D seismic data it was possible to determine the aquifer thickness. The depth of the aquifer determined by interpretation of the 3D seismic data was in accordance with the depth determined by the VES method. Based on the assumption that the hydraulic conductivity of this formation is identical or similar over the entire area, as well as the fact that the first well showed the presence of a subartesian aquifer, we proposed drilling another borehole (IBJf-2) in the zone where the data indicated that the water-bearing stratum was much thicker. The data obtained by drilling and coring were in agreement with the predicted aquifer thickness. Pumping tests showed that the water discharge in borehole IBJf-2 was 6 l/s.The results show that the objective of delineating the groundwater body by combined application of two geophysical methods (VES and 3D seismic) was successfully performed.     

Assessment of aquifer vulnerability to industrial waste water using resistivity measurements. A case study, along El-Gharbyia main drain, Nile Delta, Egypt

1D resistivity sounding and 2D resistivity imaging surveys were integrated with geological and hydrochemical data to assess the aquifer vulnerability and saltwater intrusion in the north of Nile Delta, Egypt. In the present study, the El-Gharbyia main drain was considered as a case study to map the sand bodies within the upper silt and clay aquitard. Twenty Schlumberger soundings and six 2D dipole-dipole profiles were executed along one profile close to the western side of the main drain. In addition, 14 groundwater samples and 4 surface water samples from the main drain were chemically analyzed to obtain the major and trace elements concentrations.The results from the resistivity and hydrochemical data were used to assess the protection of the groundwater aquifer and the potential risk of groundwater pollution. The inverted resistivities and thicknesses of the layers above the aquifer layer were used to estimate the integrated electrical conductivity (IEC) that can be used for quantification of aquifer vulnerability. According to the aquifer vulnerability assessment of an underlying sand aquifer, the southern part of the area is characterized by high vulnerability zone with slightly fresh to brackish groundwater and resistivity values of 11-23 Ω.m below the clay cap. The resistivity sections exhibit some sand bodies within the clay cap that lead to increase the recharging of surface waste water (650 mg/l salinity) and flushing the upper part of underlying saltwater aquifer. The region in the north has saltwater with resistivity less than 6 Ω.m and local vulnerable zones within the clay cap. The inverted 2D dipole-dipole profiles in the vulnerable zones, in combination with drilling information have allowed the identification of subsoil structure around the main drain that is highly affected by waste water.     

Electrical Resistivity Imaging and the Saline Water Interface in High-Quality Coastal Aquifers

Population growth and changing climate continue to impact on the availability of natural resources. Urbanization of vulnerable coastal margins can place serious demands on shallow groundwater. Here, groundwater management requires definition of coastal hydrogeology, particularly the seawater interface. Electrical resistivity imaging (ERI) appears to be ideally suited for this purpose. We investigate challenges and drivers for successful electrical resistivity imaging with field and synthetic experiments. Two decades of seawater intrusion monitoring provide a basis for creating a geo-electrical model suitable for demonstrating the significance of acquisition and inversion parameters on resistivity imaging outcomes. A key observation is that resistivity imaging with combinations of electrode arrays that include dipole–dipole quadrupoles can be configured to illuminate consequential elements of coastal hydrogeology. We extend our analysis of ERI to include a diverse set of hydrogeological settings along more than 100 km of the coastal margin passing the city of Perth, Western Australia. Of particular importance are settings with: (1) a classic seawater wedge in an unconfined aquifer, (2) a shallow unconfined aquifer over an impermeable substrate, and (3) a shallow multi-tiered aquifer system over a conductive impermeable substrate. We also demonstrate a systematic increase in the landward extent of the seawater wedge at sites located progressively closer to the highly urbanized center of Perth. Based on field and synthetic ERI experiments from a broad range of hydrogeological settings, we tabulate current challenges and future directions for this technology. Our research contributes to resolving the globally significant challenge of managing seawater intrusion at vulnerable coastal margins.     

Delineating coastal groundwater discharge processes in a wetland area by means of electrical resistivity imaging, 224Ra and 222Rn

Coastal groundwater discharge (CGD) plays an important role in coastal hydrogeological systems as they are a water resource that needs to be managed, particularly in wetland areas. Despite its importance, identifying and monitoring CGD often presents physical and logistical constraints, restraining the application of more traditional submarine groundwater discharge surveying techniques. Here we investigate the capability of electrical resistivity imaging (ERI) in the Peníscola wetland (Mediterranean coast, Spain). ERI surveying made it possible to identify and delineate an ascending regional groundwater flow of thermal and Ra‐enriched groundwater converging with local flows and seawater intrusion. The continuous inputs of Ra‐rich groundwater have induced high activities of Ra isotopes and ²²²Rn into the marsh area, becoming among the highest previously reported in wetlands and coastal lagoons. Geoelectrical imaging enabled inferring focused upward discharging areas, leaking from the aquifer roof through a confining unit and culminating as spring pools nourishing the wetland system. Forward modelling over idealized subsurface configurations, borehole datasets, potentiometric records from standpipe piezometers, petrophysical analysis, and four natural and independent tracers (²²⁴Ra, ²²²Rn, temperature and salinity) permitted assessing the geoelectrical model and a derived hydrogeological pattern. The research highlights the potential of ERI to improve hydrogeological characterization of subsurface processes in complex contexts, with different converging flows. Additionally, a hydrogeological conceptual model for a groundwater‐fed coastal wetland was proposed, based on the integration of surveying datasets. Copyright © 2013 John Wiley & Sons, Ltd.     

The application of ecohydrological groundwater indicators to hydrogeological conceptual models

This article reviews the application of ecohydrological indicators to hydrogeological conceptual models for earth-scientists with little or no botanical training. Ecohydrological indicators are plants whose presence or morphology can provide data about the hydrogeological setting. By examining the literature from the fields of ecohydrology, hydrogeology, geobotany, and ecology, this article summarizes what is known about groundwater indicator plants, their potential for providing information about the aquifer, and how this data can be a cost-effective addition to hydrogeological conceptual models. We conclude that the distribution and morphology of ecohydrological groundwater indicator plants can be useful to hydrogeologists in certain circumstances. They are easiest to evaluate in arid and semiarid climates. Ecohydrological groundwater indicators can provide information about the absolute depth to the water table, patterns of groundwater fluctuation, and the mineralization of the aquifer. It is shown that an understanding of the meteorological conditions of a region is often necessary to accurately interpret groundwater indicator plants and that useful data is usually obtained by observing patterns of vegetation behavior rather than interpreting individual plants. The most serious limitations to applying this source of information to hydrogeological conceptual models are the limited data in the literature and the regional nature of many indicator plants. The physical and physiological indications of the plants exist, but little effort has been made to interpret them. This article concludes by outlining several potential lines of research that could further the usefulness of ecohydrological groundwater indicators to the hydrogeological community.     

Interpretation of VLF Resistivity Data for Ground Water Contamination Surveys

Very low frequency (VLF) military communications systems provide a primary field that can be used for shallow geophysical surveys to locate ground water contamination and vertical geologic contacts. Useful properties that can be easily obtained from the interaction of the earth and the primary field are the magnitude of the vertical secondary magnetic field, the surface impedence, and the phase angle between the electrical and magnetic horizontal components. The variations in the secondary magnetic field can be related to vertical geologic contacts, such as the edges of landfill trenches. The surface impedence yields an apparent terrain conductivity, which can be used to locate low-resistivity anomalies often associated with contaminated ground water. The phase angle gives information on vertical variations in resistivity, phase angles less than 45° indicating increasing resistivity with depth. The depth of penetration of the VLF field is about one skin depth. For a frequency of 20 kHz, the skin depth in meters is approximately equal to 3.67 where p is terrain resistivity in ohmmeters.     

A comparison of helicopter‐borne electromagnetic systems for hydrogeologic studies

The increased application of airborne electromagnetic surveys to hydrogeological studies is driving a demand for data that can consistently be inverted for accurate subsurface resistivity structure from the near surface to depths of several hundred metres. We present an evaluation of three commercial airborne electromagnetic systems over two test blocks in western Nebraska, USA. The selected test blocks are representative of shallow and deep alluvial aquifer systems with low groundwater salinity and an electrically conductive base of aquifer. The aquifer units show significant lithologic heterogeneity and include both modern and ancient river systems. We compared the various data sets to one another and inverted resistivity models to borehole lithology and to ground geophysical models. We find distinct differences among the airborne electromagnetic systems as regards the spatial resolution of models, the depth of investigation, and the ability to recover near‐surface resistivity variations. We further identify systematic biases in some data sets, which we attribute to incomplete or inexact calibration or compensation procedures.     

Estimation of groundwater vulnerability to pollution based on DRASTIC in the Niipele sub-basin of the Cuvelai Etosha Basin,Namibia

Surface water is a scarce resource in Namibia with about sixty percent of Namibia's population dependent on groundwater for drinking purposes. With increasing population, the country faces water challenges and thus groundwater resources need to be managed properly. One important aspect of Integrated Water Resources Management is the protection of water resources, including protection of groundwater from contamination and over-exploitation. This study explores vulnerability mapping as a basic tool for protecting groundwater resources from pollution. It estimates groundwater vulnerability to pollution in the upper Niipele sub-basin of the Cuvelai-Etosha in Northern Namibia using the DRASTIC index. The DRASTIC index uses GIS to estimate groundwater vulnerability by overlaying different spatially referenced hydrogeological parameters that affect groundwater contamination. The study assesses the discontinuous perched aquifer (KDP) and the Ohangwena multi-layered aquifer 1 (KOH-1). For perched aquifers, point data was regionalized by a hydrotope approach whereas for KOH-1 aquifer, inverse distance weighting was used. The hydrotope approach categorized different parts of the hydrogeological system with similar properties into five hydrotopes. The result suggests that the discontinuous perched aquifers are more vulnerable than Ohangwena multi-layered aquifer 1. This implies that vulnerability increases with decreasing depth to water table because contaminants have short travel time to reach the aquifer when they are introduced on land surface. The nitrate concentration ranges between 2 and 288 mg/l in perched aquifers while in Ohangwena multi-layered aquifer 1, it ranges between 1 and 133 mg/l. It was observed that perched aquifers have high nitrate concentrations than Ohangwena 1 aquifer, which correlates well with the vulnerability results.     

Exploring Groundwater Resources and Recharge Potentialities at El-Gallaba Plain,Western Desert,Egypt

Egypt has a fast-growing population rate of 2.5%/year; consequently, there is an increase in the water demand for living and launching different development plans. Meanwhile, there is intensive construction of several dams in the upstream Nile basin countries. Thus, it is necessary to search for new water resources to overcome the expected shortages of the Nile water supply by focusing on alternative groundwater resources. El-Gallaba Plain area is one of the most promising areas in the western desert of Egypt attaining the priority for new reclamation projects; however, its hydrogeological setting is not well understood. The present work aims at identifying the recharge potential of the groundwater aquifers in El-Gallaba Plain, as well as exploring the role of geologic structures as natural conduits, and evaluating the groundwater types, origin and distribution. The integration of hydrogeophysical studies (aero and land magnetic surveys, vertical electrical sounding), hydrochemical analyses and remote sensing were successfully used for assessing the groundwater development potential. The hydrogeophysical studies show a large graben bound aquifer with thickness exceeding 220 m. The hydrochemical results indicate the presence of three major water types; Na mix, Na Cl, Na Cl HCO₃ with salinities ranging between 227 and 4324 mg/L. The aquifer receives little recharge from the western fractured calcareous plateau from past pluvial periods and scarce present flashfloods. There is no indication for recent recharge from Lake Nasser to the aquifer domain. Further modeling studies are essential for establishing sustainable abstraction levels from this aquifer.     

New insights into the hydrogeology of a basaltic shield volcano from a comparison between self-potential and electromagnetic data: Piton de la Fournaise, Indian ocean

In order to investigate aquifers, several geophysical surveys have been carried out in the Baril area of the southern flank of Piton de la Fournaise volcano on Reunion in the Indian Ocean using audiomagnetotelluric (AMT), very-low-frequency (VLF) and self-potential (SP) methods. We present the results with emphasis on a comparison between SP data and the findings of geoelectric surveys. AMT soundings have indicated, from the surface downward, three layers: (i) resistive volcanic rocks, (ii) an intermediate resistivity layer, and (iii) a conductive basement attributed to a seawater-bearing aquifer. VLF measurements allow the mapping of the first layer apparent resistivity, and therefore its bottom, when the true resistivity is supposed to be isotropic and homogenous. When this assumption does not hold, only the SP method permits the mapping of this bottom. Because of the good agreement between the SP and electromagnetic results, we propose the SP method as the first tool that should be used in studying shallow hydrogeological structures in volcanic areas.     

Geoelectrical exploration for groundwater in Al Maha Forest,Ain Jouhra,Morocco

An area of about 30 km² located in Ain Jouhra, south of Rabat, Morocco, was the subject of a geoelectric resistivity investigation. The main goal of the investigation was the assessment of the groundwater potential of the uppermost aquifer. The aquifer conditions such as depth, thickness and boundaries were also investigated. The obtained apparent resistivity curves were first analysed qualitatively and classified using simple curve shapes. Thereafter, the data were converted to resistivity and thickness pairs semi‐quantitatively by means of master curves and then quantitatively by computer modelling using ATO and Winsev software (Zohdy, 1989; Zohdy and Bisdrof, 1989). Lithological control from the available single well with a stratigraphic log aided in the correlation of the resistivity values to different rock units. Three different AB‐spacing iso‐resistivity maps, an isopach map of the main groundwater‐bearing horizon, the depth to the aquifer substratum map and five geoelectric cross‐sections were constructed. The interpretation of these soundings indicates the presence of an unconfined to semi‐confined sandy aquifer with relatively important extent and varying thickness. The maximal thickness of the aquifer is recorded in the central part of the investigated area and is thinning southwards to pinch out farther to the south. Geophysical as well as field data indicate a hydraulic connection between the upper and deeper aquifers. Indeed, the two aquifers are separated from each other by a marly substratum that is indicated throughout the area by the lowest values of the interpreted true resistivity. The value of this resistivity varies laterally, most likely due to the lateral variation in the shale‐to‐sand ratio. The altitude of the substratum decreases towards the north, and increases southwards. Regarding the availability of the groundwater in the study area, zones with high potential are theoretically expected to occur in the central part where the transversal resistance is greatest. However, sufficient water supply and high flow rates from wells intended to produce restrictively from the most upper aquifer are not likely to exist. This conclusion, which seems to be very pessimistic, is evidenced from two real field and experimental observations. The first is the rapid fall of the level of Gharnoug lake, despite the ongoing feeding by three wells. Hence, the amount of water level drop cannot be accounted for by the evaporation alone. That means that the deeper aquifer is continuously draining the upper aquifer at a high flow rate. Very low rates are recorded in all the wells that penetrated only the upper aquifer, the exception being the well that reached deeper into the lower aquifer. The flow rate in this lower aquifer measured 18 litre s^?1. Copyright © 2007 John Wiley & Sons, Ltd.