Geophysical and Hydrogeological Groundwater Prospectivity Mapping in the Kraaipan Granite-Greenstone Terrain , Southeast Botswana

Application of regional geophysical methods for hydrogeological purposes has increased over the last two decades especially in arid and semi-arid areas. A project to map the Kraaipan granite-greenstone terrain in southeast Botswana has recently been undertaken using regional aeromagnetic and gravity data with the aim to map the rocks at depth to understand the geology while the secondary objective was to subsequently assess the mineralization and groundwater potential in the area. An integrated analysis of the aeromagnetic and gravity data and their derived/processed products is hereby investigated for groundwater for drinking and agricultural purposes. The studies include: subsurface characterisation and delineation of structural framework suitable for groundwater exploration and determination of petrophysical relationships used to link the geophysical properties (e.g., density) to hydrological properties (e.g., porosity). The results of interpretation indicate that the rocks are under ~50 m of Kalahari cover and the study area is composed of three aquifers: the extensive hard rock aquifer (granitic and volcanic), the important (fractured) karst aquifer and the minor sedimentary aquifer. The area is dissected by an ENE-to-EW-trending dyke swarm visible on the regional aeromagnetic data and much clearer on high resolution aeromagnetic data. Minor fault and/or dyke elements of NW-SE and NE-SW trend are observed. Spectral analysis reveals three main average ensample interfaces at depths of 0.7 km, 1.99 km and 4.8 km. The linear Euler solutions maps reveal that the majority depths to top of magnetic bodies range from 40 m to 400 m throughout the survey area. The shallowest depths are the most significant one in this case as they probably relate to depth of bedrock and thickness of regolith or thickest sediments. For 2695 existing boreholes analysed, maximum borehole depth is 482 m (mean 108 m), and almost half (1263) were dry with another 972 havHow to cite this paper: Ranganai, R.T., Moidaki, M.D., King, J.G. and Bagai, Z.B. (2017) Geophysical and Hydrogeological Groundwater Prospectivity Mapping in the Kraaipan Granite-Greenstone Terrain, Southeast Botswana. Journal of Water Resource and Protection, 9, 1270-1298. https://doi.org/10.4236/jwarp.2017.911082 Received: July 5, 2017 Accepted: October 28, 2017 Published: October 31, 2017 Copyright © 2017 by authors and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/


Introduction
A project to map the Archaean Kraaipan greenstone belt (KGB) under approximately 50 m of Kalahari cover in southeast Botswana has recently been undertaken [1] using regional aeromagnetic and gravity data.The primary objectives were to map the rocks at depth to understand the 3D geology while the secondary objective was to subsequently promote mineral and groundwater exploration in the area.Groundwater represents the largest stock of accessible freshwater and accounts for about one-third of freshwater withdrawals globally ( [2] and references therein).There are increased demands for groundwater resources in the southern Africa region due to population increases and socio-economic development, for drinking, mining activities and associated settlements, small-scale irrigation, rural agro-industries, and agricultural and livestock purposes [3].Botswana is semi-arid and rainfall is unreliable such that surface water and groundwater recharge are limited.Therefore, to ensure that a secure amount of groundwater is available, systematic development and management planning should be established.Application of geophysical methods for hydrogeological purposes has become imperative especially in such arid and semi-arid areas in order to increase the probability of successful drilling of water wells.The principal reason for the growing interest in using geophysical methods for hydrogeological studies is that geophysics provides spatially distributed models of physical properties in regions that are difficult to sample using conventional hydrological methods [2].The geophysical models often reveal more detail compared with models derived from hydrogeological data, such as pump tests and observations of hydraulic heads.
Furthermore, geophysical methods are less invasive and comparatively inexpensive ( [4] and references therein).It is also possible to integrate GIS, remote sensing, and geophysical surveys to derive thematic layers of surface parameters such as resistivity, susceptibility, aquifer thickness/extent, and/or fault maps [5] [6] [7] [8] [9].Further, some researchers have used the geophysical data to determine important hydraulic properties [4] [10].
In this study, an integrated analysis of aeromagnetic and gravity data and their derived/processed products over a regional scale is investigated for groundwater for drinking and agricultural (mainly livestock) purposes in the Kraaipan granite-greenstone terrain in southeast Botswana.In hard rock areas it is not sufficient to consider only the weathered layer; deep saturated fractures in bedrock are also potential targets of groundwater exploration while dykes are known to be barriers to groundwater flow or preferential pathways [7] [9] [11].In addition, paleochannels and/or paleovalleys hidden in the subsurface also support an inherited preferential groundwater flow and recharge [9] [12] [13].The basis for the present hydrogeophysical survey is understanding the geological setting of the area, nature of stratigraphy and tectonism so as to have a better understanding of the prevailing hydrogeological conditions; thus delineating a buried framework suitable for groundwater flow and storage [7] [9] [14].Processing of the gravity and aeromagnetic data sets by 3D Euler deconvolution, shaded relief, derivatives and analytic signal techniques further map/locate the important structures.Thus the studies include: subsurface characterisation, delineation of structural frameworks suitable for groundwater exploration and analysis of existing boreholes.Determination of petrophysical relationships used to link the geophysical properties to hydrological properties is also considered [4] [6] [10].Eventually, the studies seek to quantify subsurface architecture that influence flow (such as hydrostratigraphy and preferential pathways); and estimate hydrological properties (such as porosity) and state variables (such as saturation) ( [15] [16] and references therein).

Regional Geology and Hydrogeology Ease of Use
The Kraaipan granite-greenstone terrain constitutes the western central part of the Archaean Kaapvaal craton in SE Botswana and comprises greenstones, felsic gneiss, migmatite and unfoliated granite [1] [17] [18] [19] [20].The Late Archaean Transvaal/Weterberg Supergroups as well as volcanic rocks of the Ventersdorp Supergroup and Tertiary sedimentary rocks of the Kalahari Group cover a significant portion of the area.Important regional geological units in the study area (24.3˚E to 25.8˚E and 25˚S to 25.8˚S; Figure 1 and Figure 2), which constitute the Ngwaketse and Southeast administrative districts, include the Kraaipan-Madibe greenstone belt, the Mmathethe granite, the Kanye and Lobatse Volcanic Formations, the younger Karoo sediments as well as Tertiary-Recent Kalahari sediments [1] [17] [19].The greenstones consist of three narrow NNW-striking belts (Stella, Kraaipan, Madibe; Figure 1) dominated by mafic metavolcanic rocks interlayered with ferruginous and siliceous metasedimentary rocks, mainly BIF and ferruginous chert [19] [21].The surrounding granitoid rocks include tonalitic and trondhjemitic gneisses (TTG), granodiorites and adamellites [17].The terrain is intruded by the Neoarchaean Gaborone Igneous Complex that includes A-type granitoids, anorthosites, rhyolites and subsidiary mafic rocks, represented by the Mmathethe granite [19] [22].The generalized lithostratigraphy in the area is summarised in Table 1.Kraaipan Group BIF 3410 ± 61/64 [17] There is a general paucity of geological outcrop in the region due to cover by the Kalahari sediments, and critically, the configuration and structure of rock units at depth was unknown until recently when geophysical interpretation was undertaken [1].Some minor banded iron formation (BIF) outcrops have been mapped in the Molopo River in the border area in the south with the main N-S trending outcrop constituting the ~20 km long Mosi ridge to the north [23].Kalahari sediments cover the region and ranges from 0 to 100 m thick as indicated by borehole data and geophysical modelling [1].The major fractures and lineaments have an ENE and WNW orientation, with secondary features trending N-S and S-W.The latter, which are often of tensional origin, may play an important part in groundwater circulation [24].Structural features associated with faulting are important targets for groundwater exploration, as it occurs along vertical and horizontal (bedding planes) structures in these rocks [25].
The central part largely comprises of igneous and metamorphic rocks  Basement granites/gneisses including amphibolites, quartzites and schists, with some calcrete deposits.
The extreme eastern part is under in the Barolong farms area [24].Aquifers with a fair yield are frequently encountered in fracture zones of the arenitic sediments of the Transvaal and Waterberg Supergroups but productivity decreases westwards where these rocks dip below thick Kalahari beds, obviously because of a reduction of recharge [24].The best groundwater resources occur in the Transvaal dolomites and the Karoo sandstones, with runoff from the slopes in quatzites and shales contributing to the replenishment of these aquifers.
Groundwater flow direction is generally to the SW and recharge is significant from the Molopo River which forms the border with South Africa, particularly during floods.
Almost 3000 boreholes have been drilled by government (BH) and private individuals (Z) within the study area (≤22,000 km 2 ) to date [31] (Figure 2).The large number of boreholes attests to the need for groundwater for domestic use and livestock in the area, and the government's commitment to provision of such social facilities.Clustering of boreholes appears to follow groundwater potential zones (Figure 2) and/or existing aquifers/wellfields (e.g., Ramotswa aquifer/wellfield).
Borehole clusters also partly follow settlements along main roads and/or railway line, close to communities.structure and poor siting (rather than groundwater potential) as the main cause of well failure in many developing countries, particularly in Precambrian basement aquifers.In the study area, the older surrounding gneisses and younger geological formations like the Kanye and Lobatse volcanics and sandstones (Ventersdop Supergroup) constitute important groundwater aquifers [24] [29] [33], and can be investigated using the available potential field data [7] [9] [14] [34].

Geophysical Data and Hydrogeological Problems/Investigations
Gravity and aeromagnetic data are routinely used for geological exploration to get information about subsurface geology under vegetation and regolith cover [35] [36].These types of methods can be used to identify rock units and major struc-  [11]) while in karst environments it is the developed solution channels and voids that are critical [29].A karst aquifer may be present even when there are no discernable karst landforms at the surface and gravity data are suitable for identifying/mapping them and the characteristic channels/cavities/voids.On the other hand, the thickness of sediments or the depth to basement is important for sedimentary aquifers, with both gravity and magnetic data yielding this information.In general, different hydrogeological problems to be solved that could benefit from geophysics include (e.g., Table 2): 1) The aquifer geometry (location of its bottom, top, and outlet) in alluvial and/or sedimentary deposits and in stratified volcanic layers; 2) Characteristics of fractured rocks (location of faults and/or fissures, characteristics of groundwater flow and circulation); 3) Knowledge of the hydraulic properties of the aquifer (porosity, clay content, transmissivity, permeability of the porous medium); 4) Quality of the water (fresh, salty, heavily mineralized, contaminated); 5) Monitoring of water flow (recharge) and interpretation pertaining to future effects of exploitation (discharge); 6) Geohazards as a result of withdrawal (e.g., in karst areas).
3) Linear magnetic lows (magnetite destruction due to pedogenesis and oxida-tion maghemite, hematite or limonite which are weak in spontaneous magne- 4) Narrow, linear magnetic highs  either dykes or magnetite enrichment from metamorphism.
The alluvial-illuvial pedogenic process also causes a substantial loss of magnetic minerals for the top soil and the sediments.Thus, the parent crystalline basement remains magnetically stronger than the top soil or sediments in general, and provides an ideal situation for the magnetic method to map basement buried below the soil/sediment cover [9] [41] [42].
A pseudo-lithology can also be derived from aeromagnetic data by characterizing magnetic facies of a rock using its percent magnetite and its spatial distribution [7] [9]  apparent susceptibility map [7].
The hydrogeological applications of the potential field data within the study area are summarized in Table 2.

Geophysical Data Processing Techniques
A number of commonly applied potential field transformations and image processing techniques were applied to both data sets to enhance both shallow, short wavelength features for lithological contact and structural mapping as well as medium wavelengths for the purpose of regional crustal structure.These include shaded relief imaging, vertical and horizontal derivatives, analytic signal, apparent density/susceptibility mapping, upward continuation, automatic gain control, directional cosine filter, 3D Euler deconvolution and spectral analysis [36] [37] [43] [44].Each technique helps to display the edges of source bodies and lateral contrasts in physical property, which are mainly caused by lithological and structural changes in the buried basement.Only those techniques with The original aeromagnetic data were subjected to a Reduction to the Pole (RTP) and International Geomagnetic Reference (IGRF) correction.The RTP transformation usually involves an assumption that, the total magnetizations of most rocks align parallel or anti-parallel to the earth's main field (declination = 12.58˚, inclination = 62.72˚ and IGRF total intensity value = 29,000 nT for the study area).
This results in anomalies lying directly over the causative bodies [37] leading to easier and improved interpretation.The IGRF correction removes the earth's main geomagnetic field for the particular year of survey, resulting in magnetic responses mainly due to geological features.
The initial improvement of structural elements from the aeromagnetic data was achieved by shaded relief and colour-shadow maps that rely on artificial illumination of the data mimicking topography [43] [45].The combination of the images allows accentuation of linear features perpendicular to the light source while information on the anomaly magnitude or intensity is retained in the colour component [42] [44].They also effectively display both short wavelength and medium wavelength anomalies as well as curvilinear and elliptical features (e.g., plutons and/or mafic sills) [44].
Derivatives of the potential field data are probably the most important functions for both qualitative and quantitative geological interpretation.Various combinations of vertical and horizontal derivatives were used to further remove or suppress the regional trends in the data, and assist in qualitative edge detection of the linear geological structures (dykes, faults and shear zones) [37] [38] [42] [44].The analytic signal amplitude (ASA) of potential field data is the magnitude of the squares of the derivatives, and for magnetic data it has the added advantage of independence of the magnetization direction of the source [36] [37] [46] [47].The ASA is symmetrical and peaks directly over the edges of wide bodies and the centres of narrow bodies thus providing obvious indications of dense or magnetic source geometry [42].The most useful and popular quantitative technique that uses derivatives is 3D Euler deconvolution [48] which calculates the spatial position and depth of the gravity/magnetic body, thus fully locating the source.Standard 3D Euler deconvolution is used for relatively rapid delineation of linear features (structural mapping) and contacts between the sedimentary rocks and basement rocks, and rapid depth estimation.However, the quality and/or accuracy of the results, particularly the depth estimates, depends on adequate sampling of the data, proper choice of structural index and window size The following list summarizes the important geophysical processes applied, geological structures mapped, and corresponding hydrogeological significance.
Euler deconvolution; spectral analysis  sediment thickness, depth to bedrock  maximum thickness of potential aquifer.

Aeromagnetic Data and Geostructural Correlation
Botswana and South Africa are both completely covered by regional aeromagnetic data (1 km line spacing, 150 m mean terrain clearance, ~70 m sample spacing); with higher resolution aeromagnetic data for selected areas (line spacing 150 m, flight altitude 70 m, ~20 m sample spacing).Merged data are presented in Figure 4 and Figure 5 below as colour-shadow maps, with merging processes covered in detail in [1].
To correlate magnetic anomalies with rock units, it is noteworthy that sedimentary rocks are generally non-magnetic, whereas igneous rocks rich in iron and magnesium (mafic to ultramafic) tend to be very magnetic [9] [36] [41].Granite intrusions and hornfels contact aureoles can also be magnetic [52].In this paper all aeromagnetic results are displayed as digital image maps; they are reduction to the pole map, colour-shaded total intensity map and vertical derivative map.Many more images were generated during this study, but are not presented for clarity (i.e., to reduce redundancy).The general interpretation is also aided by higher resolution on-screen images.
The map in Figure 4  thick banded iron formations (BIF) which have strong aeromagnetic signatures and which can thereby be followed along strike for considerable distances from their outcrop in the south west of Sedibeng (Figure 4).The two major magnetic anomalies (approximately 180 km in length) trend nearly N-S direction.The southeastern part (mainly in South Africa) is occupied by a series of elongated and relatively higher magnetic intensities trending in an EW-to-ENE direction

Gravity Data and Geological Correlation
The gravity data used in this study are described in detail by [1] and the station distribution and simple Bouguer gravity anomaly grid are shown in Figure 6.A total of one hundred and forty (140) new gravity data points were collected, merged with other existing gravity data, and then used to construct the gravity map of the study area, which is generally dominated by a broad high.More data are currently being acquired to improve the regional station coverage/density.The gravity maps (e.g., Figure 6) indicate that the Kraaipan greenstone rocks correlate with a major NS-trending gravity high (at Molopo, Tswaaneng and Metlobo villages), and is believed to be caused by magnetite quarzites, banded iron formations and metavolcanic rocks [1].Their peak values reach over −90 mGal, which is the strongest anomaly within the study area.Other notable highs occur over (are coincident with) the Segwagwa Complex and between Goodhope and Mmathethe in the central part, and over the Molopo farms (ultramafic) complex in the west (south of Lerolwane).The Bouguer anomaly map observations Figure 6.Simple Bouguer gravity anomaly grid/map and actual distribution of gravity stations (+ is for existing data in Botswana acquired from Geological Survey, circle represents gravity points in South Africa and triangle represents gravity data acquired in this study.Location of profile A-A' modeled for subsurface geological sections (Figure 7) is indicated.
(Figure 6) thus support the aeromagnetic anomaly maps (Figure 4 and Figure 5).There are several gravity lows with peak magnitudes of −123 m Gals also recognized in this map.The map reveals prominent gravity lows in circular and oval shapes distributed at the central and southern parts of the mapped area.It is clearly observed from several gravity maps that meta-clastic rocks at Sekhutlane village have low gravity intensities.Moreover, the results of this study show more extent and new granitoid intrusions than the previous geological mapping.
South of the study area in South Africa, three linear (and narrow) gravity highs correlate with the Madibe (M, around Mafikeng), Kraaipan (K) and Stella (S) greenstone belts (Figure 6), dominated by mafic metavolcanic rocks interlayered with ferruginous and siliceous metasedimentary rocks, mainly BIF and ferruginous chert [1].These are bounded by a series of moderate to low gravity anomalies resulting in a typical volcanic-granitic gravity character, with multiple small blocks, high gradients and quick succession of positive and negative anomalies.

Geophysical Results and Interpretation
As basement heterogeneities should affect both gravity and magnetic fields, we are preferentially looking for collocated contrasts in both of these fields.However, it should be noted that regional gravity data are commonly used to delineate the position of boundaries at deeper levels than aeromagnetic data [1] [27].The following results summarise the most significant gravity and magnetic anomalies of which at least some should represent basement faults or other significant structures within the Kraaipan granite-greenstone terrain which are of importance to groundwater exploration.

2 3 /4D Gravity Modelling and Sediment Thickness
The configurations of the main geological units in the study area were determined along one selected regional profile (see Figure 6 and Figure 7) with data sampled from the Bouguer gravity map along a detailed field traverse (~2 km spacing).The density contrast values for the various rock units and the basement were estimated from [27] [52] [53] and [54].A regional gravity structural cross-section (Figure 7) was derived/obtained using GM-Sys' 2 3 / 4 D-forward modelling technique.The gravity field was calculated iteratively for these geological models, until a good fit was reached between the observed (dots) and calculated (line) profiles.Different density values (within the standard deviation of ±30 kg/m 3 ) were tested but these did not significantly change the shapes and depth extents of the units, as noted/observed by others [54] [55].
The modelling results show that the sediment thickness and/or depth to (granitic) basement is at most 60 m, with the largest values in the west (Figure 7).
Newly mapped granitic plutons on the Bouguer anomaly map (Figure 6) are also indicated.

Depth to Susceptibility Discontinuities
The application of depth determination techniques, to determine the subsurface layer thicknesses and estimate depth to top of the causative geological bodies, is one of the most important parameters that must be outlined in order to interpret adequately the geology and structure of the area.Here, two techniques are used to calculate the depths of the causative bodies utilizing the magnetic data: spectral frequency analysis and 3D Euler deconvolution method.

Source Depth Estimates from Magnetic Spectra
Subsurface Geology: Layers & Thicknesses The Fast Fourier Transformation (FFT) was applied on the RTP aeromagnetic survey data and a two-dimensional (radially averaged) power spectrum curve obtained (Figure 8).Four main average ensample interfaces at depths of 0.3 km (S4), 0.5 km (S3), 0.86 km (S2) and 2.85 km (S1) are revealed for shallowest layer, intermediate sources and for the deepest ensemble, respectively.
The first two shallowest depths are the most significant as they map thickness of sediments and/or the weathered zone (depth to magnetic basement) important for groundwater accumulation/storage.

Depths from Euler Deconvolution of Magnetic Data Structural Mapping and Depth Estimates
The Euler deconvolution technique was carried out on the magnetic data using Several linear features indicated on the RTP colour-shadow map (Figure 5) are now clearly evident, with the boundaries of the primary elements and major domains clearly recognized in the Euler depth map (Figure 9).New structures are also evident, and their locations and depths confirmed.

Discussion
Over any area of past or present active tectonics, surface structural alignments are often visible as expressions on topographic, geologic and geophysical maps.These lineaments become less distinct with age but are still important as indicators of past stresses, displacements, and disturbance patterns significant for groundwater transport (recharge) and storage/accumulation.Because the buried sediment-basement contact is not affected by erosion, these structural patterns remain imprinted in it and could be revealed from potential fields.Alignment of anomalies in geophysical maps provides important information for structural analysis; however, it still remains an open question whether all or most magnetic and gravity lineaments correspond to basement faults, folds, or other structural patterns important for exploration.
As stated in the earlier sections, results from the enhanced aeromagnetic data and merged gravity data conform well with the previous field geological mapping, specialized geochronological, geochemical and petrological sampling by several geoscience pioneers [17] [18] [19] [20].However, these enhanced merged data together with existing aeromagnetic data yield more essential geological information, particularly those regarding contexts of tectonism and mineral/groundwater exploration.First, as suggested by [1], two major gravity and magnetic highs coincide with the western and central greenstone belts caused by magnetite and quarzites, banded iron formations and metavolcanic rocks.
Second, and more importantly, the applied processes/techniques have yielded subsurface layer thicknesses as well as location and depths of the several structures.
In this discussion, special emphases are placed on the relationship of the interpreted gravity and airborne magnetic data to geology and structure of the Kraaipan greenstone terrain and their application to groundwater exploration.

Borehole Depth in Relation to Gravity and Magnetic Basement
The borehole depths were gridded using a minimum curvature technique resulting in a grid reflecting water strike depth ranges for the study area (Figure 10), interpreted to be depth to water table.The western part/third of the area constitutes sedimentary aquifer with deep water table, followed by central area of granite-gneiss with mixed levels and then fairly shallow water table in the karst (dolomite) aquifer in the northeast.

Borehole Yield in Relation to Aquifer Types and/or Compartments/Units
The various observed and/or interpreted gravity and magnetic zones and structures correspond to different aquifer types or units (e.g., Figure 2) as alluded to above.There are high magnetic field intensities observed at both western and central zones of the study area, so the linear magnetic bodies in these zones suggest that these features are caused by intrusive rocks which correspond well with several of the known granitoid rocks comprising tonalitic and trondhjenilitic gneisses, granodiorites and adamelites with BIF horizons.Narrow and higher amplitude magnetic anomalies in western and the central parts coincided with gravity highs, particularly at Phitshane-Molopo and Mabule, suggest the existence of mafic volcanic rocks.Field investigation reveals that the rocks in these areas are made up mainly of banded iron formations (BIF), magnetite quartzites and metavolcanic rocks [19].The borehole yields (Figure 3) were also gridded using a minimum curvature technique resulting in a relatively smooth grid (Figure 11) for correlation with the gravity and magnetic zones and structures.11).This is attributed to the relatively thin (≤60 m) sedimentary cover as derived from the modelling (Figure 7).However, the extreme western study area corresponding to metalvolcanic and ultramafic rocks (Molopo farms complex) also suggests a NS-trending zone of high yields.The known karst aquifer (Ramotswa aquifer) around the Otse-Lobatse area is also reflected as a high yield zone (as expected).Overall, the map can be used to prioritise further exploration zones.

Borehole Yields in Relation to Structural Features
In this study, extraction of the geological structures and smaller-scale lineaments was based on combining multiple data attributes.Our lineament extraction mostly relied on transformations of the magnetic and gravity fields that included the vertical and horizontal gradients, analytic signal and several variants of local wavenumber mapping.Among these methods, the vertical gradient (Figure 12) and analytic signal were the most important.Thus the initial pattern of the lineaments was drawn mainly based on the vertical gradient contrasts, followed by and their refinement using other attribute maps extracted from the magnetic and gravity fields.By combining the observations from various attribute maps and reducing the effects of noise and processing artefacts, an interpretation map of the high resolution aeromagnetic data depicting major lineaments was obtained (Figure 13).
The study area is basically dominated with lineament features which can be confidently recognized from the magnetic data (Figure 4, Figure 5, Figure 9 and Figure 12).The magnetic anomalies display several trends defined by alignment of intensities, gradients and shapes of anomalies, and are best illustrated in ~EW in the south of the study area (Figure 12).The latter (~EW) regional structures discernable across the entire study area appear to be older than the others based on cross-cutting relations, and overprinting of the basement blocks of different tectonic origins and potential-field signatures (Figure 4).The aeromagnetic maps show the long continuity of lineaments of the moderate to high amplitude anomalies with the roughly north-south direction in the central zone corresponding to the greenstone (BIF) structure in geologic maps [17] [19] [20].The north-south structure is cross-cut by the east-trending lineaments.This result is in good agreement with the regional structure reported earlier by previous authors, which showed the correlation of the regional structure with the geochemical and petrological data [17] [19] [20].
We considered whether such structures could be potentially related to faulting associated with basement rocks and thus open for groundwater recharge or in-filled dykes acting as groundwater barriers.The classified existing borehole yields (Figure 3) were thus overprinted on the magnetic vertical derivative map (Figure 12) and the dominant lineaments interpreted from the HRAM data (Figure 13).On a regional scale, most of the high yield boreholes (>10 m 3 /hr) are associated with magnetic structures or features.It is evident from Figure 12 and Figure 13 that there is very little correlation between borehole positions and geophysical structures.However there is some correlation between medium to high yield boreholes and some of the major structures.

Data Synthesis
A synthesis of the various maps is presented in Table 3 based on the composite presentation in Figure 14 below, and both help define the best groundwater potential zones (A, B, C in order of priority).

Conclusions and Recommendations
Kraaipan granite-greenstone geophysical and hydrogeological data from South-East Botswana have been processed and interpreted.The results of the magnetic and gravity methods used revealed that the lithologies and tectonic framework of the Archaean basement rocks extend into Southeast Botswana from South Africa.In particular, the derivative, analytic signal and Euler deconvolution maps reveal the complex geological structure of the terrain and basement rocks suitable for groundwater potential assessment.
Figure 14.3D composite map display of the magnetic first vertical derivative, borehole yield, water strike depth and total depth (highest values in red/purple down to lowest values in deep blue) for direct correlation of these variables (Table 3).12 and Figure 14, Table 3).However, the ~E-W magnetic lineament structures appear to cut across the entire study area.were jointly evaluated and the most promising locations for positioning the hydro-wells were proposed.An analysis of various types of aquifers shows that the shallow fractured-karst aquifers can be described as very good, the shallow fractured crystalline and metamorphic aquifers as good, while the deep sedimentary aquifers are fair, in terms of quantity of groundwater.The main results and conclusions can be summarised as follows: 1) Identified dykes, faults, joints, zones of weathering and sediment thickness.In conclusion, the study shows the importance of preliminary geophysical investigations before ground borehole siting and drilling in order to improve borehole success rates and/or reduce costs.The described exploration scheme takes into account the fact that the end-users usually require that geophysical survey cost must be only a small portion of the total hydro well construction budget.
This study revealed that had the geophysical data analysis been available at the time of siting and drilling it would have guided the drilling programme and substantially reduced the number of dry boreholes and hence saved financial resources.
Further studies are planned to include determination of aquifer hydraulic properties [4] [6] [10] [57] and recharge issues-geochemistry/isotopes [4] [12] [56] which are beyond the scope of this paper but within the overall project/research goals.This will assist the sustainable management of the aquifers in this era of climate change [2].

Figure 1 .
Figure 1.Simplified geology of Kraaipan granite-greenstone terrain in SE Botswana and NW South Africa based on geological mapping and geophysical interpretation [1].

Figure 3 .
Figure 3. Borehole yield (m 3 /hr) statistics subdivided into the three main groups of poor, moderate and high or very good.Note the large number of dry boreholes.

[ 49 ]
[50].Fourier spectral analysis (or wavelength filtering) in recent years has also become a widely utilized tool for the processing and interpretation of potential field data.Energy spectral analysis provides a technique for quantitative studies of large and complex aeromagnetic or gravity data sets.The logarithm of the radial average of the energy spectrum (the square of the Fourier amplitude spectrum) is plotted versus the radial frequency.The slope of each segment provides R. T. Ranganai et al.DOI: 10.4236/jwarp.2017.9110821280 Journal of Water Resource and Protection the depth to the top of an ensemble of magnetic or gravity bodies.A typical energy spectrum for magnetic data may exhibit three parts; a deep source component, a shallow source component and a noise component.For wavenumbers in cycles/km, the depth to a statistical ensemble of sources is determined from the expression: slope = −4πh [51].
Figure 4. Regional IGRF corrected colour-shaded aeromagnetic data of the study area, showing several structural trends representing BIF horizons, dykes and faults.(Red = magnetic high; Blue = magnetic low) Blocks show high resolution aeromagnetic data coverage (e.g., Figure 5) for mineral exploration.

(Figure 5 .
Figure 5. Detailed total magnetic intensity (TMI) colour-shadow map of exploration area showing the N-S trending BIF and NE-to-ENE-trending dykes and faults.Box indicates further higher resolution mineral exploration area (magnetic, electromagnetic, radiometric data).

Figure 8 .
Figure 8. Radially averaged power spectrum of the regional RTP aeromagnetic data.The depth to a statistical ensemble of sources is determined from the expression: slope = −4πh [51].

Figure 9 .
Figure 9. Euler deconvolution solutions for high resolution aeromagnetic data for SI = 0 (W = 400 m) emphasizing shallow (≤350 m) linear features.This SI is applicable to both finite and infinite contact/fault in magnetics [49] [50].

Figure 10 .
Figure 10.Gridded borehole water strike depth reflecting water level in different aquifers, indicating relatively shallow levels in the eastern half (karst) and deep in the west (sedimentary?).Location of profile A-A' modeled for subsurface geological sections (Figure 7) is indicated.

Figure 11 .
Figure 11.Borehole yield grid in relation to known/interpreted aquifer types and/or compartments.Note the NNE-trending zone of high yields in the central part of the area/map.Location of profile A-A' modeled for subsurface geological sections (Figure 7) is indicated.

Figure 9 and
Figure 9 and Figure 12.The most prominent magnetic lineaments crossing the entire Madibe-Kraaipan terrain are at Eastern part of the survey area.These structures change their trends from NW-SE in the northeast to NE-SW and

5 m 3 /
The area is dissected by an ENE-to-EW-trending Proterozoic dyke swarm visible on the high resolution aeromagnetic data (Figure5and Figure13), and associated with the initial crustal faulting.The results of the depths interpretation indicate that the rocks are under >50 m of Kalahari cover and the study area is composed of three aquifers: the extensive hard rock aquifer (granitic and volcanic), the important (fractured) karst aquifer and the minor sedimentary aquifer.Spectral analysis reveals three main average ensample interfaces at depths of 0.7 km, 1.99 km and 4.8 km.The linear Euler solutions maps reveal that the majority depths to top of magnetic bodies range from 40 m to 400 m throughout the survey area.The shallowest depths are the most significant one for groundwater as they probably relate to depth of bedrock and thickness of regolith or thickest sediments.For 2695 existing boreholes analysed, maximum borehole depth is 482 m (mean 108 m), and almost half (1263) were dry with another 972 having low yield (1hr) and 432 yielding 6 -49 m 3 /hr and only 28 above 50 m 3 /hr (maximum ~160 m 3 /hr) and an average water strike of 64 m.Finally, all geophysical results

Table 1 .
Summary lithostratigraphic table of the study area.

Table 2 .
Hydrogeological applications of gravity/magnetic data and products.
maps/figures presented or discussed are briefly described below.

Table 3 .
Summary hydrogeophysical variables and characteristics: groundwater potential zones.Three main regions with different internal structural patterns and hydrogeological characteristics were deduced: 1) eastern zone dominated by NW-SE and NE-SW features; 2) central crystalline zone with mainly E-W structures and the major N-S BIF horizon; and 3) western sedimentary zone marked by NW-SE to WNW-ESE lineaments (Figure