Hydrological Evaluation of the Groundwater Potential in the Fractured Karoo Aquifer Using Magnetic and Electrical Resistivity Methods: Case Study of the Balfour Formation, Alice, South Africa

,


Introduction
Groundwater resource is becoming essential and accounts for a significant proportion of the water used for various purposes [1]. Many rural communities rely on dams and rivers' surface water supplies, which might be completely or almost dry during dry seasons [2]. Groundwater is mostly odorless and contains low dissolved solids [3]. The benefits of groundwater, being the source of potable water, cannot be underestimated, mainly in places whereby the population is predominantly rural, and demand is distributed across wide areas [4]. Finding sustainable, hygienic, and potable water is a big challenge that never ends because it promotes the development of every community [5]. Surface water is becoming increasingly scarce due to the strain of industrial growth in several areas of the globe. Due to this, water scarcity is now the third most significant global problem since freshwater makes up only 3% of all water on Earth [6]. Groundwater is the primary source of drinking water for over 60% of people worldwide [7]. It makes up roughly 40% of the total water supply by the government in the UK, additionally about 99% in Denmark, 60% in the USA [7,8], and approximately 80% in Germany [7]. In South Africa, 22% of communities rely solely on groundwater; another 34% combine it with surface water [9].
The limited and uneven water distribution in South Africa is severely impacted by climate change and the development of invasive alien plant species [10]. Owing to population growth and migration, the development of rural areas puts a strain on South Africa's water system. According to [11], nine million South Africans are not provided with drinkable water, and one-third of the country's population lacks access to water infrastructure. In addition, about five provinces in South Africa were announced as disaster zones in 2015 as a result of the severe drought that unfortunately hit the nation. Several countries in semiarid and arid areas, including South Africa, are classified as water-stressed [11] and rely on groundwater as their primary water source for numerous purposes. Studying the properties of aquifers and the pattern of groundwater movement is vital for the country's water resource development.
Groundwater investigation is often based on assessing the hydrological indicators and validation using a geological and geophysical application in hydrology [12]. Fault zones, lineaments, and fractures change the characteristics of a hydrogeological medium [13]. Several studies have established that drilling wells through faults or cracks increase borehole yield [14][15][16]. [17] notes that the magnetic technique is the fast and most affordable method for mapping the weathered/fractured basement and investigating surface-to-subsurface geological formations. Numerous authors have successfully applied the magnetic technique to map the subsurface features [17][18][19]. Electromagnetic, electrical resistivity, or seismic surveys are usually conducted after a magnetic survey to obtain accurate subsurface information. The electrical resistivity technique is the most often utilized geophysical method for groundwater studies, among others, as it could give a strong response under the existing subsurface conditions [20,21].
The integrated applications of magnetic and electrical resistivity techniques provide a better understanding of subsurface geology [22]. Many researchers have used magnetic and resistivity methods as a combined geophysical tool to efficiently investigate groundwater [23][24][25]. Ground magnetic survey and vertical electrical sounding (VES) were combined to understand the subsurface geology and ascertain groundwater availability in Alice. Furthermore, this was done to show the benefits of magnetic and resistivity methods as proficient geophysical tools for investigating deep and shallow aquifers.
One of the provinces struggling to provide clean water to its residents is the Eastern Cape, where Alice is situated (Figure 1). There is essentially no groundwater in most aquifers [26]. The municipality supplies Alice with piped water every day. There was no water supply scarcity from the reservoir in the municipality until the reservoir turned out to be the major source of water supply to recently developed neighborhoods near Alice. Previous studies conducted in Alice focused on the contaminant delineation of a landfill site [27] using combined induced polarization and electrical resistivity methods and the characterization of the weathered/fractured aquifer using geological mapping, data collection, and borehole drilling and logging [28]. To date, no studies have been conducted to evaluate Alice's groundwater resources. This knowledge gap could hinder efficient groundwater production. To better understand Alice's groundwater resources, a thorough investigation is required to know the overburden thickness and depth to bedrock and vertical and horizontal changes in electrical resistivity with depth and select subsurface aquifer zones for borehole drilling.

Geology
The Karoo Basin is a sedimentary basin that covers about 600,000 km 2 of central and southern South Africa. Its sedimentation occurred about 300 million years ago during the late Carboniferous, up until the breakup of Gondwana in the Middle Jurassic. The breakup was mainly because of the deposition of a large igneous province [29], whereby dolerite sills and dykes intruded the surrounding sediments and induced the formation of fractures. The dolerite sills and dykes are now responsible for the recharge and drainage patterns, as well as controlling the morphology and influencing the emergence of many seepages and springs [28]. [30] reported that this basin is a retro-arc foreland basin, which occurred when the Paleo-pacific plate underwent subduction at a shallow angle underneath Gondwana. The Karoo Supergroup is stratigraphically divided into Dwyka, Ecca, Beaufort, Stormberg, and Drakensberg Groups [31]. The groups are mainly composed of sedimentary sequences, apart from the Drakensberg Group, comprising volcanic rocks [31]. [32] reported that the Karoo Basin protects a stratum of glaciomarine to terrestrial sequences that are formed in a diversity of tectonically influenced depositories under increasingly arid climatic settings. The southern part of the depositional environment of the Karoo varied from glacial (Dwyka) to fluvial (Ecca) over time.
The Dwyka Group is formed at the bottom of the Karoo Basin (Table 1). In the southern part of the Karoo, the group unconformably overlies the Cape Supergroup [33]. The group is 600 m-750 m thick, comprising diamictites and shales [34]. It is roughly 2340 m thick [34]. It comprises six formations, which are Prince Albert, Whitehill, Collingham, Ripon, Fort Brown, and Waterford Formations [28]. The Beaufort Group, consisting of the lower Adelaide Subgroup and the upper Tarkastad Subgroup [28], comprises alternating sandstones and mudrocks and has attained a thickness above 500 m and a surface area of about 200,000 km 2 [34]. The group's age ranges from the Middle Permian to the Middle Triassic period [35].
The study area is situated in the Balfour Formation within the Beaufort Group. It is made up of a succession of sedimentary rocks, including a succession of mudstones with subordinate interbedded sandstones, siltstone, and shale. These rocks were eventually intruded by the dolerites in the Karoo, forming sills and dykes [36]. The variation in its lithology is characterized by alternating sandstone-rich and mudstone-rich members [37]. This formation is divided into five stratigraphic members, namely, Oudeberg, Daggaboersnek, Barberskrans, Elandsberg, and Palingkloof [35]. Alice is geologically located in the Daggaboersnek Member ( Figure 2). The Daggaboersnek Member is made up of sandstones with mudrock intercalations. According to [38], shale, siltstone, sandstone, and mudstones make up the majority of the low-permeability Karoo aquifers. Many boreholes drilled in the Karoo Formation are recorded to be of low yields [26].

Materials and Methods
3.1. Magnetic Study. Magnetic data were acquired from 8458 points in Alice using the G-857 magnetometer. Two magnetometers were used for the ground magnetic survey. The magnetometers were checked to know whether they both measured similar values at a given time in the same place. One magnetometer was fixed at the base station to account for diurnal variation. There were no magnetic objects in the area where this base station was located. The base station readings were taken at an interval of 30 s. The magnetic intensity of the research area was measured using a second magnetometer along the traverse. The intertraverse spacing was 20 m. Readings along the traverses were recorded at an interval of 5 s. These traverses were located far from significant metallic objects. The magnetic readings and time of the recording were taken and stored automatically by the magnetometer. The time of recording is needed for the correction of diurnal variations. The Garmin GPS was mounted on the G-857 magnetometer nonmetallic sensor pole, and it was configured to record locations when taking magnetic readings along the traverse.
3.1.1. Magnetic Data Analysis. The ground magnetic data were processed using Geosoft software, Oasis Montaj 9.10, for filtering and enhancement to improve the interpretive power of the magnetic data. Data enhancements (i.e., reduced magnetic pole, total horizontal derivative, vertical derivative, and analytical signal) were applied to the residual magnetic data to identify the edges of geological features and enhance deeper magnetic sources [41,42]. Magnetic depth slicing was successfully conducted with the aid of GETECH GETGRID software. The magnetic results are displayed using spatial maps ("Reduced to the pole," "Total horizontal derivative," "Vertical derivative," and "Analytical signal" sections). The enhancement techniques are summarized below.
(1) Reduced to the Pole (RTP). The analysis of the residual magnetic field intensity map is generally inaccurate due to the movement of magnetic anomalies away from their magnetic sources [41]. The magnetic inclination of 63.63°and declination of -28.02°were applied to the total magnetic field map after removing the diurnal variation to produce the RTP map [43]. The RTP map allows the elimination of anomalous skewness while maintaining the data configuration [44]. This computation process aligns the boundaries of the magnetic anomalies with their buried sources [45]. This was accomplished by using Equation (1) to determine the differentiation of the source components in the vertical direction [45].
where T y ðrÞ = RTP magnetic anomaly depends on the vertical magnetization (y), r is the radius of the magnetization field from the magnetic source to the point of observation, ρ is the direction of the anomalous magnetization, and b is the direction of the magnetic field of the earth magnetic field.
(2) Analytical Signal (AS) Mapping. The analytical signal (AS) is produced by combining the vertical and horizontal gradients of magnetic anomalies. The analytical signal has both direction and amplitude. However, the analytical signal is generally defined by its "amplitude function" (i.e., the 3 International Journal of Geophysics square root of the sum of the squares of the horizontal and vertical derivatives of the Earth's magnetic field) [46]. The data in the wavenumber domain was processed and filtered using an inverse fast Fourier transform (FFT) algorithm. The filter's capacity to provide a total value precisely over the magnetic source, as well as a depth estimate, makes it a suitable method for interpreting magnetic data [47]. [48] states that the benefit of this magnetic data enhancement technique is that it does not depend on specific preconceptions regarding the body's magnetization direction, and it usually has a positive amplitude function. The analytical signal was computed using the following equation: where Aðx, yÞ is the amplitude of the AS at ðx, yÞ, f is the observed magnetic field at ðx, yÞ, and x, y, and z are the edges of the magnetic structure.
(3) First Vertical Derivative (FVD). The vertical gradients of potential fields are potentially important in interpreting magnetic data because they map out linear structures such as dykes or faults [39]. [49] noted that the FVD technique is useful in enhancing magnetic anomalies caused by shallow sources, narrowing the size of the anomalies, and identifying the causative bodies precisely. The vertical derivative map responds very well to local effects than regional impacts, producing a sharper image than the total field magnetic map [49].
(4) Total Horizontal Derivative (THD). Horizontal derivatives (THD) are magnetic vector data that provide more information about the directional variations of the total magnetic field [50]. They are crucial in mapping linear features like dykes or fault zones from magnetic data [51]. Horizontal derivative maps are derivative products that point out discontinuities in anomalous patterns and textures. The Pythagorean sum of the gradients in orthogonal directions is used to compute an anomaly field's horizontal derivative. The horizontal gradient magnitude THDðx, yÞ for the magnetic field Fðx, yÞ is given by [52] as follows: where THDðx, yÞ is the amplitude of the THD at ðx, yÞ, f is the observed magnetic field at ðx, yÞ, and x and y are the edges of the magnetic structure.
(5) Depth Slicing Technique. This technique uses the linear filter to estimate the depth of the magnetic anomaly source as a plot of the gradient of the Earth's magnetic field all through the subsidence of the magnetic anomaly [50]. [39] noted that the magnetic signal is considered the outcome of several uncorrelated random processes, according to the Wiener filtering principle, which is the basis for depth slicing. Outlining the impacts of shallow magnetic sources from deeper magnetic sources enables the differentiation of many data components. The majority of magnetic anomalies seen at the surface typically come from shallow depths, whereas unclear magnetic anomalies generally come from deeper depths. The depth estimation is calculated using the following equation: where h represents the depth of the anomalous body and b represents the Nyquist wavenumber component's slope.   (Figure 3). The vertical electrical sounding (VES) method using the Schlumberger electrode configuration (Figure 3), which is effective in the investigation of groundwater, was used for the resistivity measurements to determine the apparent resistivity [53]. The Schlumberger array is successfully utilized for groundwater studies because of its flexibility, easy interpretation, and durability of the instrument used for the survey [54]. The survey was completed at twenty-five (25) VES stations within the study area ( Figure 4). Four electrodes were arranged collinearly, and two current electrodes, A and B, were used for deep current penetration ( Figure 3) while measuring the potential difference using another pair of electrodes (M and N) [55]. The current electrodes (AB) were gradually moved outwardly, leaving the potential difference electrodes (MN) fixed to increase the depth range of the particular station being investigated. At some points, the ratio AB/MN became too large, which led to a drop in the potential difference values, making the readings inaccurate. Therefore, it was necessary to increase the potential difference electrodes (MN), and the distance between the electrodes was kept at MN ≤ AB/5 in order to provide accurate resistivity data. The potential difference readings from the instrument were maintained at a value greater than 5 mV. The current electrodes (AB/2) reached a distance of 250 m, while the potential electrodes (MN) started from 0.5 m and subsequently extended to a suitable distance as the current electrodes (AB) were moved symmetrically.

Vertical Electrical
For the Schlumberger electrode configuration, the instrument measured the apparent resistivity (ρ a ) using the following equation: where k is the geometric factor, ΔV is the potential difference, and I is the electric current. The geometric factor (k) for the Schlumberger array is calculated using the following equation: where s is the distance between two current electrodes AB and a is the distance between two potential electrodes MN. The results of the VES are presented as electrical sounding curves and geoelectric cross-sections.

Data Analysis and Resistivity
Interpretation. The obtained VES data were analyzed with the aid of the partial curve matching technique and iteration modeling method using WINRESIST software to produce the geoelectric models and SURFER-10 software to generate the geoelectric sections [56]. The partial curve matching technique is based on a preliminary evaluation that provides the thickness and resistivity of different geoelectric layers [57]. The thickness, resistivity, and depth of a layer are taken into account when deciding if a layer is fractured or weathered, which makes the layer suitable for groundwater storage. The decision is accomplished by comparing each layer's resistivity measurements with the standard resistivity values of groundwater.
The resistivity values of groundwater vary between 10 and 100 Ωm since the quantity of dissolved salts in groundwater lowers its resistivity [58]. The VES points with a layer having resistivity values between 10 and 200 Ωm, large thickness, and a depth greater than 30 m are therefore chosen, considering the water table and depths of the existing borehole from previous studies [26,28].  International Journal of Geophysics

Results
The geophysical techniques of magnetic and electrical resistivity surveys were combined to examine the subsurface geological structures which influence the distribution of groundwater using magnetic data and evaluate the thickness of water-bearing areas and depth using the electrical resistivity data, hence identifying potential groundwater zones.

Magnetic.
The results of the magnetic survey are presented in maps. The residual magnetic map is presented in Figure 5. The amplitude of the total magnetic map ranges from 21756.8 to 28514.9 nT. However, this map could not be accurately interpreted owing to the effect of the nonzero magnetic inclination. Therefore, it was reduced to the pole (RTP) and displayed as an RTP map in Figure 6.      (Figure 10(e)) does not show the narrow anomalies caused by shallow magnetic sources that can be seen on the depth slices 1 to 4 maps (Figures 10(a)-10(d)). This is due to the shallow sources' disappearance at a depth of about 31 m.

Vertical Electrical Sounding (VES).
The VES interpretation showed four geoelectric layers, and the curves consist of HA and HK types. About 64% of the overall soundings are of the HA curve type, whereas 36% are of the HK type ( Figure 11).   Table 2. In this article, one representative VES curve from each category is presented (Figure 12). The HA curve type is denoted by VES modeled curve site 17 acquired along the University of Fort Hare Farm (Figure 12(a)). The curve begins with a top layer having a resistivity value of 31 Ωm and a thickness of 1.1 m. The resistivity of the second layer is 11 Ωm, and the thickness is 2.7 m, which is inferred to be mudstone. The third layer's resistivity is 136 Ωm, and its thickness is 30.2 m. This conductive layer is likely to be weathered/fractured sandstone. The fourth layer's resistivity is 5512 Ωm. The resistive fourth layer may be due to the presence of bedrock that has been unweathered. The thickness of the fourth layer is undefined, given that it is the last layer. The HK curve type is represented by VES modeled curve site 4 obtained along the University Cricket Field (Figure 12(b)). The surface terrain is composed of resistive topsoil (127 Ωm), and its thickness is 1.8 m. The second layer is 7.4 m thick and has a low resistivity value of 29 Ωm, which is inferred to be mudstone. The third layer is 72.6 m thick and has a resistivity value of 195 Ωm, which corresponds to sandstone. The resistivity value of the fourth layer is 112 Ωm, which corresponds to mudstone.  Figures 13(a) and 13(b).
The geoelectric section acquired along a line in East Camp is presented in Figure 13

Discussion
5.1. Magnetic. The maximum and minimum anomaly frequencies, which are connected to the local distribution of magnetization, are displayed by the magnetic maps. A change in magnetic susceptibility generally implies that the rock's composition has changed [60]. Magnetic anomalies are shown on magnetic maps based on their characteristics and structure. [61] states that the shapes of the anomalies enable the identification of the fundamental geological features of the buried source. The RTP magnetic map of the study area shows different regions of low and high intensity. The low magnetic regions consist of rocks with low magnetic susceptibility, situated in Alice's central and southwest parts. The high magnetic zones, mainly the northeastern and southwestern parts, occur owing to magnetic sources at depth [59] and are inferred to be Karoo dolerite intrusions. According to [62], the RTP map's high-gradient regions may suggest the occurrence of high-gradient zones. [63]   11 International Journal of Geophysics documents that knowing the geometry and distribution of the Karoo dolerite will assist groundwater exploration companies in drilling wells to intersect the dolerite-dominated areas. Hence, the THD map reveals several dolerite intrusions (dykes and sills), which may be targeted for potential groundwater exploration. Several researchers [26,64,65] have documented that most boreholes drilled alongside linear geological structures, mainly dolerite dykes, are more productive than those drilled away from geological structures.
Meanwhile, there is a positive correlation between the low magnetic areas seen on the RTP map and the location of the dolerite intrusions on the THD map, which implies that the low magnetic intensity of RTP possibly represents the area of occurrence of dolerite dykes and sills on the THD map. This finding aligns with [66], who stated that the dolerites are embedded in neotectonic fractures. However, the significant aquifer in the Karoo is anticipated to have low yield and permeability at the contact area of host rocks and dolerites, which aligns with the reports of [16,65]. The aquifer system in Alice is possibly hosted in the sedimentary rocks within the Daggaboersnek Member of the Karoo Supergroup [26,37]. [28] conducted research that focuses on characterizing, borehole drilling, and measuring the flow parameters of the rock aquifer system within the current study area. Superimposed on the AS map are the two (2) drilled borehole points in the study area ( Figure 14). The boreholes are sited in the NE region. This was done to evaluate the correlation between geological features (sills and dykes) and the borehole locations to examine the influence of these geological features on groundwater flow using the borehole data. holes. It appears that the formation may be extensively fractured because it is extremely loose between depths of 54 and 65 m. Figure 14 shows that the boreholes are drilled on a silllike structure in the NE part. Meanwhile, [64] noted that fracturing might be improved at the bottom of curving sills or dyke-sill intersection points, thus influencing the migration and distribution of groundwater. Therefore, the loose formation and the high blow yield noted at depths of 54 and 65 m indicated that the boreholes might have been drilled through a fracture caused by the dolerite intrusion. The finding aligns with [14], who posited that the dolerite intrusion in the study area caused the area's formation to   Figure 14: Overlay of the two boreholes on the analytical signal map of Alice for correlation purposes. 13 International Journal of Geophysics be extensively fractured, which made them suitable aquifers due to secondary permeability. In addition, it was noted that the Permian mudstones and sandstones serve as the primary aquifers in the Karoo Basin, as well as the dolerite dykes and sills. Furthermore, the finding on the high blow yield due to boreholes drilled through the fracture aligns with the report of [67] on the positive yield analysis of boreholes drilled into a fractured structure. According to the classification of several well-known rocks using resistivity values given by [68], layers having low resistivity (1-20 Ωm) in the study area could probably be associated with clay, and the intermediate resistivity layers (20-5000 Ωm) could be related to mudstone and sandstone. High resistivity values (>5000 Ωm) of the bedrock observed are associated with sediments found in the Balfour Formation, which were intruded by dolerites [27].
In addition, the findings on outcrops and data from past studies [26,28,34,36,37,40] aided in designing the lithological 1D models associated with each sounding curve using the 1D geoelectrical models generated. Therefore, it is acknowledged that the study area shows various lithological models composed of different types of sedimentary rocks   14 International Journal of Geophysics (sandstone, mudstone, and clays, which are intruded by dolerites) [40]. The findings align with [34], who stated that the Balfour Formation comprises sedimentary rocks, like sandstone, shale, mudstone, and alternating siltstone. The lithological 1D models from the HA-and HK-type curves from the study area show that a viable aquifer could be found in the third layer, which has resistivity values between 136 and 352 Ωm, 9.9 to 172 m thick, and depths ranging from 16.1 to 163 m. This third layer comprises sandstone and mudstone, in which the major substantial aquifers located within the Karoo Supergroup are found [26,28,34,37]. In the geoelectric section (Figure 13(a)), the sandstone layer is more pronounced and, in some places, intruded by dolerites. The noticeable deformation brought on by the dolerites present in the Karoo enhances the permeability of its aquifer system in the contact areas of the dyke. In VES 24 and 25, the dolerite intruded almost to the surface of this layer. The two VES stations are situated near Somgxada Mountain in the study area.
The geoelectric section shows an electrical discontinuity between VES stations 21 and 23, which could be caused by the passage of a fault that influences the aquifer's groundwater flow. The low gaps between high resistive bodies indicate fault regions that may serve as paths for groundwater passage, thus increasing groundwater localization [68]. Similarly, [69] posited that fault zones have excellent storage capacity and are highly productive. The aquifer around VES 22 lies within the weathered/fractured sandstone, about 59.7 m thick, with a depth of 62 m. The deepest aquifer at approximately 110 m depth could be seen at VES station 14 ( Figure 13(b)). This VES 14 station aquifer seems preferable for groundwater exploitation because of its greater depth and low resistivity value. The geoelectric sections (Figures 13(a) and 13(b)) show that the aquifer zone could be found in the third layer (i.e., sandstone). The finding corroborates the results of [16], who suggested that the most important type of aquifers found in the Karoo are developed in weathered sedimentary rocks (sandstone, mudstone, siltstone, and shale) characterized by low permeability.
To enable the lithocorrelation, two vertical electrical soundings (VES 6 and VES 8) were located beside two boreholes (BH1 and BH2). The geoelectric models and sections were created using other VES stations. The resistivity records of VES 6 and 8 show a high correspondence with the borehole data of BH1 and BH 2 (Figures 15(a)  and 15(b)), where the lithology at a depth of 28 m at the borehole BH1 is sandstone and at VES 6 is sandstone. The depth of the sandstone aquifer at VES 8 is above 48 m, while borehole BH2 strikes the water at depths of 54 and 66 m. Therefore, the results of the geoelectrical data correlated well with the borehole log data from the study area's boreholes drilled by [28].

Conclusions
This research shows the vital contribution of the magnetic and resistivity methods as a practical tool for investigating potential groundwater zones in Alice, Eastern Cape, South Africa, and other areas where applicable. The ground mag-netic results show that the anomaly signature amplitudes vary, suggesting that the magnetic body source is not uniformly distributed along each profile across the study area. The variation between positive and negative magnetic anomalies indicates the presence of subsurface lineaments, which may influence the groundwater movement in the study area. As observed on the magnetic depth slice maps, the magnetic anomalies become broader up to the depth of 18.8 m and later fade away at a depth of around 31 m. This finding suggests that Karoo dolerite intrusions were intruded close to the surface. According to the findings of the VES, the subsurface is composed of four layers of HK and HA types that characterize the types of sounding curves and the vertical changes in Alice. From the results of the VES curves, 15 VES stations (VES 1, 2, 3, 4, 5, 6, 7, 8, 9, 12, 13, 14, 15, 17, and 22) are selected as potential groundwater zones, with aquifer layer resistivity values between 136 and 352 Ωm. The thickness of the aquifer layer varies between 9.9 and 143.9 m, with depths between 16.1 and 163 m. The aquifer zone is inferred to be the third layer, and the depth of the aquifer zone is considered between 30 and 160 m. However, it is recommended that water should be collected from each successful borehole for bacteriological and physicochemical tests to ensure that the water is safe and suitable for various purposes. The results from the geoelectric sections further revealed areas of weakness, interpreted as weathered/fractured zones accountable for groundwater accumulation in the study area. Low resistivity zones indicate highly faulted and weathered/fractured zones; therefore, they are suitable for groundwater storage. Meanwhile, high resistivities are linked to dolerite intrusions. The borehole lithologic log correlated well with the magnetic and electrical resistivity survey results. Thus, these geophysical methods are efficient in locating potential groundwater zones. The results will guide the drilling of boreholes aiming to harness the groundwater resources in Alice. International Journal of Geophysics