Aquifer Vulnerability Assessment Based on Sequence Stratigraphic and 39Ar Transport Modeling

A large‐scale groundwater flow and transport model is developed for a deep‐seated (100 to 300 m below ground surface) sedimentary aquifer system. The model is based on a three‐dimensional (3D) hydrostratigraphic model, building on a sequence stratigraphic approach. The flow model is calibrated against observations of hydraulic head and stream discharge while the credibility of the transport model is evaluated against measurements of ³⁹Ar from deep wells using alternative parameterizations of dispersivity and effective porosity. The directly simulated 3D mean age distributions and vertical fluxes are used to visualize the two‐dimensional (2D)/3D age and flux distribution along transects and at the top plane of individual aquifers. The simulation results are used to assess the vulnerability of the aquifer system that generally has been assumed to be protected by thick overlaying clayey units and therefore proposed as future reservoirs for drinking water supply. The results indicate that on a regional scale these deep‐seated aquifers are not as protected from modern surface water contamination as expected because significant leakage to the deeper aquifers occurs. The complex distribution of local and intermediate groundwater flow systems controlled by the distribution of the river network as well as the topographical variation (Tóth 1963) provides the possibility for modern water to be found in even the deepest aquifers.     

Three-dimensional flow in the storative semiconfining layers of a leaky aquifer

An analytical solution for three-dimensional (3D) flow in the storative semiconfining layers of a leaky aquifer fully penetrated by a production well is developed in this article to provide a method from which accurate hydraulic parameters in the semiconfining layers can be derived from aquifer test data. The analysis of synthetic aquifer test data with the 3D analytical solution in the semiconfining layers provided more accurate optimal hydraulic parameters than those derived using the available quasi-two-dimensional (2D) solution. Differences between the 3D and 2D flow solutions in the semiconfining layers become larger when a no flow boundary condition is imposed at either at the top of the upper semiconfining layer or at the bottom of the lower semiconfining layer or when the hydraulic conductivity ratio of the semiconfining layer to the aquifer is larger than 0.001. In addition, differences between the 3D and 2D flow solutions in the semiconfining layers are illustrated when the thickness ratio of the semiconfining layer to the aquifer is changed. Analysis of water level data from two hypothetical and one real aquifer test showed that the 3D solution in the semiconfining layers provides lower correlation coefficients among hydraulic parameters than the 2D solution.     

A Finite Layer Formulation for Groundwater Flow to Horizontal Wells

A finite layer approach for the general problem of three‐dimensional (3D) flow to horizontal wells in multilayered aquifer systems is presented, in which the unconfined flow can be taken into account. The flow is approximated by an integration of the standard finite element method in vertical direction and the analytical techniques in the other spatial directions. Because only the vertical discretization is involved, the horizontal wells can be completely contained in one specific nodal plane without discretization. Moreover, due to the analytical eigenfunctions introduced in the formulation, the weighted residual equations can be decoupled, and the formulas for the global matrices and flow vector corresponding to horizontal wells can be obtained explicitly. Consequently, the bandwidth of the global matrices and computational cost rising from 3D analysis can be significantly reduced. Two comparisons to the existing solutions are made to verify the validity of the formulation, including transient flow to horizontal wells in confined and unconfined aquifers. Furthermore, an additional numerical application to horizontal wells in three‐layered systems is presented to demonstrate the applicability of the present method in modeling flow in more complex aquifer systems.     

A methodology for regionalizing 3‐D effective porosity at watershed scale in crystalline aquifers

An innovative approach for regionalizing the 3‐D effective porosity field is presented and applied to two large, overexploited, and deeply weathered crystalline aquifers located in southern India. The method derives from earlier work on regionalizing a 2‐D effective porosity field in that part of an aquifer where the water table fluctuates, which is now extended over the entire aquifer using a 3‐D approach. A method based on geological and geophysical surveys has also been developed for mapping the weathering profile layers (saprolite and fractured layers). The method for regionalizing 3‐D effective porosity combines water table fluctuation and groundwater budget techniques at various cell sizes with the use of satellite‐based data (for groundwater abstraction), the structure of the weathering profile, and geostatistical techniques. The approach is presented in detail for the Kudaliar watershed (983 km²) and tested on the 730 km² Anantapur watershed. At watershed scale, the effective porosity of the aquifer ranges from 0.5% to 2% in Kudaliar and between 0.3% and 1% in Anantapur, which agrees with earlier works. Results show that (a) depending on the geology and on the structure of the weathering profile, the vertical distribution of effective porosity can be very different and that the fractured layers in crystalline aquifers are not necessarily characterized by a rapid decrease in effective porosity and (b) that the lateral variations in effective porosity can be larger than the vertical ones. These variations suggest that within a same weathering profile, the density of open fractures and/or degree of weathering in the fractured zone may significantly vary from a place to another. The proposed method provides information on the spatial distribution of effective porosity that is of prime interest in terms of flux and contaminant transport in crystalline aquifers. Implications for mapping groundwater storage and scarcity are also discussed, which should help in improving groundwater resource management strategies.     

CO2 storage in the high Arctic: efficient modelling of pre‐stack depth‐migrated seismic sections for survey planning

The sequestration of CO₂ in subsurface reservoirs constitutes an immediate counter‐measure to reduce anthropogenic emissions of CO₂, now recognized by international scientific panels to be the single most critical factor driving the observed global climatic warming. To ensure and verify the safe geological containment of CO₂ underground, monitoring of the CO₂ site is critical. In the high Arctic, environmental considerations are paramount and human impact through, for instance, active seismic surveys, has to be minimized. Efficient seismic modelling is a powerful tool to test the detectability and imaging capability prior to acquisition and thus improve the characterization of CO₂ storage sites, taking both geological setting and seismic acquisition set‐up into account. The unique method presented here avoids the costly generation of large synthetic data sets by employing point spread functions to directly generate pre‐stack depth‐migrated seismic images. We test both a local‐target approach using an analytical filter assuming an average velocity and a full‐field approach accounting for the spatial variability of point spread functions. We assume a hypothetical CO₂ plume emplaced in a sloping aquifer inspired by the conditions found at the University of Svalbard CO₂ lab close to Longyearbyen, Svalbard, Norway, constituting an unconventional reservoir–cap rock system. Using the local‐target approach, we find that even the low‐to‐moderate values of porosity (5%–18%) measured in the reservoir should be sufficient to induce significant change in seismic response when CO₂ is injected. The sensitivity of the seismic response to changes in CO₂ saturation, however, is limited once a relatively low saturation threshold of 5% is exceeded. Depending on the illumination angle provided by the seismic survey, the quality of the images of five hypothetical CO₂ plumes of varying volume differs depending on the steepness of their flanks. When comparing the resolution of two orthogonal 2D surveys to a 3D survey, we discover that the images of the 2D surveys contain significant artefacts, the CO₂‐brine contact is misplaced and an additional reflector is introduced due to the projection of the point spread function of the unresolvable plane onto the imaging plane. All of these could easily lead to a misinterpretation of the behaviour of the injected CO₂. Our workflow allows for testing the influence of geological heterogeneities in the target aquifer (igneous intrusions, faults, pervasive fracture networks) by utilizing increasingly complex and more realistic geological models as input as more information on the subsurface becomes available.     

A Genetic Algorithm for Stochastic Inversion in Contaminant Subsurface Hydrology

Identifying the spatial distribution of hydrological properties of aquifers is a key problem in subsurface hydrology. The aquifer structure plays an important role in contaminant transport. Identifying the properties (primarily the hydraulic conductivity) is essentially an inversion problem that is ill-posed, non-unique and computationally intensive by definition. In this work, the non-uniqueness of the inverse problem is tackled via a novel Genetic Algorithm approach combined with a geostatistical method (Sequential Indicator Simulations) for construction of realizations of properties spatial distributions, which are modeled as random. The Genetic Algorithm cross-over operator is based on a novel concept of pilot-planes: daughter realizations adopt pilot-planes from one of their parents. In addition, each aquifer realization is conditioned on the geological hard data and is constructed by sampling the facies distribution, evaluated by indicator variograms. The approach is illustrated in two test cases: a synthetic two-dimensional (2D) case and an actual three-dimensional (3D) case. The results have shown the ability of the proposed approach to generate a set of realizations, where each individual exhibits minor deviations from the measurements. Further, a comparison between the proposed approach and direct (Monte Carlo) approach shows that the Genetic Algorithm was able to generate an ensemble of solutions with a better fitting of the measurements than the direct approach by a significantly reduced computational effort.     

Conformal mapping of groundwater flow fields with internal boundaries

An analytical approach is presented for solving problems of steady, two-dimensional groundwater flow with inhomogeneity boundaries. A common approach for such problems is to separate the problem domain into two homogeneous domains, search for solutions in each domain, and then attempt to match conditions, either exactly or approximately, along the inhomogeneity boundary. Here, we use classical solutions to problems with inhomogeneity boundaries with simple geometries, and map conformally the entire domain onto a new one. In this way, existing solutions are used to solve problems with more complex, and more practical, boundary geometries. The approach is general, but subject to some restrictions on the mapping functions that may be used.Using this approach, we develop explicit analytical solutions for two problems of practical interest. The first problem addresses aquifer interaction across a gap in an impermeable separating layer; flow regimes are defined and the interaction is quantified. The second solution represents flow in the vertical plane to a partially clogged stream bed that is partially penetrating the aquifer; the stream bed is modeled as a thin layer of low-permeability silt. Flow regimes for groundwater surface–water interaction are quantified analytically.     

A distance-drawdown aquifer test method for aquifers with areal anisotropy

A new distance-drawdown method for aquifers with anisotropy on the horizontal plane is presented. The method uses scalar transformation to convert to an equivalent, isotropic medium, thus permitting application of the Cooper-Jacob Method. The method is applicable to cases where at least one ellipse of equal drawdown can be delineated but can also be applied where no ellipse can be discerned from the data. In the latter case, a least-squares regression approach can be employed to estimate the orientation and magnitude of the anisotropy. The regression R2 value provides a quantitative assessment of the degree to which the drawdown data are indicative of a systematic areal anisotropy in the aquifer or whether the data simply reflect natural aquifer heterogeneity. In addition to confined aquifers, this methodology, like the Cooper-Jacob Method, is also applicable to unconfined aquifers either before the onset of delayed drainage or following the completion of delayed drainage provided that the u value meets the recommended criterion.     

Key features of artificial aquifers for use in modeling contaminant transport

Two large-scale (9.5 m long, 4.7 m wide, 2.6 m deep), three-dimensional artificial aquifers were constructed to investigate the influence of spatial variations in aquifer properties on contaminant transport. One aquifer was uniformly filled with coarse sand media (0.6 to 2.0 mm) and the other was constructed as a heterogeneous aquifer using blocks of fine, medium, and coarse sands. The key features of these artificial aquifers are described. An innovative deaeration tower was constructed to overcome a problem of the aquifers becoming blocked with excess air from the ground water source. A series of tracer injection experiments were conducted to test the homogeneity of the first aquifer that was purposely built as a homogeneous aquifer and to calculate values of aquifer parameters. Experimental data show that the aquifer is slightly heterogeneous, and hydraulic conductivity values are significantly higher down one side of the aquifer compared to the mean value. There was very good agreement in estimated dispersivity values between the plume area ratio methods and the curve fitting of tracer breakthrough curves. Dispersivity estimates from a full areal source injection (12.2 m2) experiment using a 1D analytical model were higher than estimates from a limited source injection (0.2 m2) experiment using a 3D model, possibly because the 1D model does not take account of the heterogeneity of hydraulic conductivity in the aquifer, thus overestimating dispersivity. Transverse and vertical dispersivity values were about five times less than the longitudinal dispersivity. There was slight sorption of Rhodamine WT onto the aquifer media.     

Coupled groundwater flow and heat transport simulation for estimating transient aquifer–stream exchange at the lowland River Spree (Germany)

Subsurface flow and heat transport near Freienbrink, NE Germany, was simulated in order to study groundwater–surface water exchange between a floodplains aquifer and a section of the lowland River Spree and an adjacent oxbow. Groundwater exfiltration was the dominant process, and only fast surface water level rises resulted in temporary infiltration into the aquifer. The main groundwater flow paths are identified based on a 3D groundwater flow model. To estimate mass fluxes across the aquifer–surface water interfaces, a 2D flow and heat transport modelling approach along a transect of 12 piezometers was performed. Results of steady‐state and transient water level simulations show an overall high accuracy with a Spearman coefficient ρ = 0.9996 and root mean square error (RMSE) = 0.008 m. Based on small groundwater flow velocities of about 10^?7 to 10^?6 ms^?1, mean groundwater exfiltration rates of 233 l m^?2 d^?1 are calculated. Short periods of surface water infiltration into the aquifer do not exceed 10 days, and the infiltration rates are in the same range. The heat transport was modelled with slightly less accuracy (ρ = 0.8359 and RMSE = 0.34 °C). In contrast to the predominant groundwater exfiltration, surface water temperatures determine the calculated temperatures in the upper aquifer below both surface water bodies down to 10 m during the whole simulation period. These findings emphasize prevailing of heat conduction over advection in the upper aquifer zones, which seems to be typical for lowland streams with sandy aquifer materials and low hydraulic gradients. Moreover, this study shows the potential of coupled numerical flow and heat transport modelling to understand groundwater–surface water exchange processes in detail. Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

Flow‐Through Stream Modeling with MODFLOW and MT3D: Certainties and Limitations

This paper aims to assess MODFLOW and MT3D capabilities for simulating the spread of contaminants from a river exhibiting an unusual relationship with an alluvial aquifer, with the groundwater head higher than the river head on one side and lower on the other (flow‐through stream). A series of simulation tests is conducted using a simple hypothetical model so as to characterize and quantify these limitations. Simulation results show that the expected contaminant spread could be achieved with a specific configuration composed of two sets of parameters: (1) modeled object parameters (hydraulic groundwater gradient, hydraulic conductivity values of aquifer and streambed), and (2) modeling parameters (vertical discretization of aquifer, horizontal refinement of stream modeled with River [RIV] package). The influence of these various parameters on simulation results is investigated, and potential complications and errors are identified. Contaminant spread from stream to aquifer is not always reproduced by MT3D due to the RIV package's inability to simulate lateral exchange fluxes between stream and aquifer. This paper identifies the need for a MODFLOW streamflow package allowing lateral stream‐aquifer interactions and streamflow routine calculations. Such developments could be of particular interest for modeling contaminated flow‐through streams.     

Modeling vertical stratification of CO2 injected into a deep layered aquifer

The vertical stratification of carbon dioxide (CO₂) injected into a deep layered aquifer made up of high-permeability and low-permeability layers, such as Utsira aquifer at Sleipner site in Norway, is investigated with a Buckley–Leverett equation including gravity effects. In a first step, we study both by theory and simulation the application of this equation to the vertical migration of a light phase (CO₂), in a denser phase (water), in 1D vertical columns filled with different types of porous media: homogeneous, piecewise homogeneous, layered periodic and finally heterogeneous. For each case, we solve the associated Riemann problems and propose semi-analytical solutions describing the spatial and temporal evolution of the light phase saturation. These solutions agree well with simulation results. We show that the flux continuity condition at interfaces between high-permeability and low-permeability layers leads to CO₂ saturation discontinuities at these interfaces and, in particular, to a saturation increase beneath low-permeability layers. In a second step, we analyze the vertical migration of a CO₂ plume injected into a 2D layered aquifer. We show that the CO₂ vertical stratification under each low-permeability layer is induced, as in 1D columns, by the flux continuity condition at interfaces. As the injection takes place at the bottom of the aquifer the velocity and the flux function decrease with elevation and this phenomenon is proposed to explain the stratification under each mudstone layer as observed at Sleipner site.     

Characterizing moldic and vuggy pore space in karst aquifers using borehole-wall,slabbed-core and thin-section images

Carbonate aquifers are prolific and important sources of potable water in many parts of the world owing to enlarged dissolution features that enhance porosity and interconnectivity. To better understand the variations of pore space in different karst aquifers, image and geospatial analyses are used to analyze pore attributes (i.e., pore area and perimeter) in images of vuggy aquifers. Pore geometry and 2D porosity derived from images of the moldic Castle Hayne and vuggy Biscayne aquifers are analyzed at three scales of observation: borehole televiewer, core and thin-section. The Castle Hayne and Biscayne aquifers are the foci of this study because the pore spaces that control the hydrologic properties in each of these aquifers are markedly different even though both of these carbonate reservoirs are prolific aquifers. Assessments of pore area, perimeter and shape index (a measure of shape complexity) indicate that pore geometries and pore complexities vary as a function pore type and scale of observation. For each aquifer type, the areas, perimeters and complexities of pores are higher at the larger scale of observation (e.g., borehole) than the smaller scale of observation (e.g., thin section). When the complexity of the moldic pores is compared to the complexity of vuggy pores, the results indicate that moldic pores are generally more complex than vuggy pores at the same scale of observation. Whereas estimates of 2D porosity from the borehole televiewer image of the vuggy aquifer are higher than those derived from the moldic aquifer, the range of 2D porosities is larger in core and thin section images for the vuggy aquifer than the moldic aquifer. A model for the development of pores is presented that suggests that the coalescence of small pores with simple shapes leads to the growth of larger pores with more complex shapes. The model suggests that the younger Biscayne aquifer is a more mature karst than the Castle Hayne aquifer.     

Sustainable groundwater management of the stressed Coastal aquifer in the Gaza region

Abstract

This study proposes an empirical approach that can lead to the sustainable management of groundwater resources. This approach enables a comprehensive understanding of an aquifer, delineates distinct hydrological scenarios, and recommends a set of operational activities for each sub-region of the aquifer. The paper focuses on the Coastal aquifer of the Gaza Strip region which has been divided into three sub-regions. The southern sub-region (WSW) is classified as scenario “+a2”, which indicates that it can be used as a multi-annual groundwater reservoir. The northern sub-region (NW-E) is designated scenario “-a2”, where the recommended operational measures include injection of freshwater in wells and cleaning of the surface environment. The third sub-region (CSE), is classified as scenario “-b2”, which requires severe management measures to correct both a negative hydrological and environmental situation. The approach also involves on-going monitoring of the aquifer, and can be considered as an empirical tool to provide preliminary guidelines for long-term groundwater management.     

Simulating pumping-induced regional land subsidence with the use of InSAR and field data in the Toluca Valley, Mexico

A multidisciplinary approach is presented here for quantifying land subsidence in a heavily pumped aquifer system with complex stratigraphy. The methodology consists in incorporating Terzaghi’s 1D instantaneous compaction principle into a 3D groundwater flow model that is then applied and calibrated to reproduce observed hydraulic heads and compaction for the Toluca Valley, Mexico. Differential Interferometric Synthetic Aperture Radar (D-InSAR), a generated 3D-geological model, extensometers, monitoring wells, and available literature are used to constrain the model. The D-InSAR measured subsidence, extensometers, and numerical simulations of subsidence agree relatively well. Simulations show that since regional subsidence began in the mid 1960s there has been up to 2 m of subsidence in the industrial corridor, where heavy pumping and thick clay layers are found. This study shows that an approach using various sources of data is useful in estimating and constraining the vertical component of the inelastic skeletal specific storage.     

A coupled vertically integrated model to describe lateral exchanges between surface and subsurface in large alluvial floodplains with a fully penetrating river

This paper presents a vertically averaged model for studying water and solute exchanges between a large river and its adjacent alluvial aquifer. The hydraulic model couples horizontal 2D Saint Venant equations for river flow and a 2D Dupuit equation for aquifer flow. The dynamic coupling between river and aquifer is provided by continuity of fluxes and water level elevation between the two domains. Equations are solved simultaneously by linking the two hydrological system matrices in a single global matrix in order to ensure the continuity conditions between river and aquifer and to accurately model two‐way coupling between these two domains. The model is applied to a large reach (about 36 km²) of the Garonne River (south‐western France) and its floodplain, including an instrumented site in a meander. Simulated hydraulic heads are compared with experimental measurements on the Garonne River and aquifer in the floodplain. Model verification includes comparisons for one point sampling date (27 piezometers, 30 March 2000) and for hydraulic heads variations measured continuously over 5 months (5 piezometers, 1 January to 1 June 2000). The model accurately reproduces the strong hydraulic connections between the Garonne River and its aquifer, which are confirmed by the simultaneous variation of the water level in the river and in piezometers located near the river bank. The simulations also confirmed that the model is able to reproduce groundwater flow dynamics during flood events. Given these results, the hydraulic model was coupled with a solute‐transport component, based on advection‐dispersion equations, to investigate the theoretical dynamics of a conservative tracer over 5 years throughout the 36 km² reach studied. Meanders were shown to favour exchanges between river and aquifer, and although the tracer was diluted in the river, the contamination moved downstream from the injection plots and affected both river banks. Copyright © 2008 John Wiley & Sons, Ltd.     

An analytical approach to capture zone delineation for a well near a stream with a leaky layer

Abstract

An analytical solution is developed to delineate the capture zone of a pumping well in an aquifer with a regional flow perpendicular to a stream, assuming a leaky layer between the stream and the aquifer. Three different scenarios are considered for different pumping rates. At low pumping rates, the capture zone boundary will be completely contained in the aquifer. At medium pumping rates, the tip of the capture zone boundary will intrude into the leaky layer. Under these two scenarios, all the pumped water is supplied from the regional groundwater flow in the aquifer. At high pumping rates, however, the capture zone boundary intersects the stream and pumped water is supplied from both the aquifer and the stream. The two critical pumping rates which separate these three scenarios, as well as the proportion of pumped water from the stream and the aquifer, are determined for different hydraulic settings.

Editor D. Koutsoyiannis; Associate editor A. Koussis

Citation Asadi-Aghbolaghi, M., Rakhshandehroo, G.R., and Kompani-Zare, M., 2013. An analytical approach to capture zone delineation for a well near a stream with a leaky layer. Hydrological Sciences Journal, 58 (8), 1813–1823.     

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.