A project to map the Archaean Kraaipan greenstone belt (KGB) under approximately 50 m of Kalahari cover in southeast Botswana has recently been undertaken  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 (  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  . 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  . 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 (  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      . Further, some researchers have used the geophysical data to determine important hydraulic properties   .
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    . In addition, paleochannels and/or paleovalleys hidden in the subsurface also support an inherited preferential groundwater flow and recharge    . 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    . 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    . 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) (   and references therein).
2. 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      . 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    . 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   . The surrounding granitoid rocks include tonalitic and trondhjemitic gneisses (TTG), granodiorites and adamellites  . 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   . The generalized lithostratigraphy in the area is summarised in Table 1.
Figure 1. Simplified geology of Kraaipan granite-greenstone terrain in SE Botswana and NW South Africa based on geological mapping and geophysical interpretation  .
Table 1. Summary lithostratigraphic table of the study area.
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  . 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  . Kalahari sediments cover the region and ranges from 0 to 100 m thick as indicated by borehole data and geophysical modelling  . 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  . Structural features associated with faulting are important targets for groundwater exploration, as it occurs along vertical and horizontal (bedding planes) structures in these rocks  .
The Botswana National Master Plan Study   recognises 4 aquifer types based on the type of water transmissive voids:
1) Fractured à Ghanzi Group, patches of Waterberg, Transvaal, Ventersdorp, and Archaean rocks;
2) Fractured porous à most common (e.g. Ntane Sandstone, Ecca & Dwyka groups);
3) Porous à Kalahari sands and alluvium (perched aquifers);
4) Karstified fractured à least common (dolomite and limestone) (SE Botswana).
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 Transvaal and Waterberg Supergroups with minor/some igneous and metamorphic rocks à Dolomite formation: dolomite and subordinate chert, shale and quartzite; quartzite, shale, sandstone and conglomerate. Ventersdop Supergroup à Lobatse/Kanye volcanics: tuff, felsite, greywacke, felsic porphyry and agglomerate.
Within the study area, the western part is mainly under Kalahari beds and Karoo Supergroup rocks with minor Waterberg Supergroup à sand, clay, calcrete and pan sediments (up to 100 m thick); mudstone, sandstone, carbonaceous shale; quartzite, shale, sandstone and conglomerate.
With regards groundwater potential in the study area (Figure 2), the granite complex and the volcanics in general show poor aquifer characteristics (low primary porosity) because of a lack of well-connected open fracture systems although fair yields are reported from boreholes in volcanics west of Lobatse. The gneisses in general form slightly better water-bearing layers with several fair-yielding boreholes in the Barolong farms area  . 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  . The best groundwater resources occur in the Transvaal dolomites and the Karoo sandstones,
Figure 2. Distribution of existing boreholes in relation to groundwater potential zones(yellow = poor,but locally fair,light green=fair,but variable and dark green= high but bariable) in part of the study area(source: Zones/areas with white dots indicate occasional or no recharge with only the eastern and norteast parts indicated as receiving regular recharge).
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 km2) to date  (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.
There were a few data base problems noted, mainly data not consistently reliable; thus borehole screening included discarding/rejecting those with water strike depth and/or yield but with zero total drill depth (89), with the remainder sorted by yield. A few others have doubtful location, plotting in South Africa. For 2695 existing boreholes retained and analysed, maximum borehole depth is 482 (mean 108), and almost half (1263) were dry with another 972 having low yield (1 - 5 m3/hr) and 432 yielding 6 - 49 m3/hr and only 28 above 50 m3/hr (maximum ~160 m3/hr) (Figure 3). Average water strike depth is 64 m and the maximum is 310 m (maximum drill depth of 482 m). There is multiple water strike for some boreholes (710), suggesting presence of confined aquifers.
BHs with WStrk ≥ 1 m à 1719
BHs with Yield ≥ 1 m3/hr à 1432
BHs with Yield ≥ 6 m3/hr à 460
References    note poor understanding of the hydrogeological
Figure 3. Borehole yield (m3/hr) statistics subdivided into the three main groups of poor, moderate and high or very good. Note the large number of dry boreholes.
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    , and can be investigated using the available potential field data     .
3. 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   . These types of methods can be used to identify rock units and major structures by applying filter enhancements, image processing and interpretation methods   . In addition, these products are useful for identifying target areas for mineral and groundwater exploration       .
Groundwater in crystalline/metamorphic areas is mainly controlled by two factors: tectonic fracturing (intensity & intersections) and weathering development (deep/thick desirable    ) while in karst environments it is the developed solution channels and voids that are critical  . 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).
Faults are usually interpreted from aeromagnetic (AM) data as:
1) Sharp gradients, boundaries or contacts between areas of different magnetic relief, level (magnitude) or texture à morphology    .
2) Disruptions and/or deflections of magnetic zones/features/trends à displacement (wrench faults/shears) with lateral movements    .
3) Linear magnetic lows (magnetite destruction due to pedogenesis and oxidationà maghemite, hematite or limonite which are weak in spontaneous magnetization at an ambient temperature   ).
4) Narrow, linear magnetic highs à either dykes or magnetite enrichment from metamorphism.
5) Weathering à areas of less quartz-rich rocks/facies   .
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    .
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   à apparent susceptibility map  .
The hydrogeological applications of the potential field data within the study area are summarized in Table 2.
4. 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     . 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
Table 2. Hydrogeological applications of gravity/magnetic data and products.
maps/figures presented or discussed are briefly described below.
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  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   . 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   . 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)  .
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)     . 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     . 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  . The most useful and popular quantitative technique that uses derivatives is 3D Euler deconvolution  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   .
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 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  .
The following list summarizes the important geophysical processes applied, geological structures mapped, and corresponding hydrogeological significance.
Pseudoshading and derivatives à structural mapping (faults, dykes, BIFs) à either groundwater conduits/pathways or barriers.
Susceptibility; Euler deconvolution; spectral analysis à saprolite/regolith mapping (weathered zone thickness; depth to bedrock) à approximate thickness of potential aquifer.
Euler deconvolution; spectral analysis à sediment thickness, depth to bedrock à maximum thickness of potential aquifer.
5. Geophysical Data
5.1. 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  .
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    . Granite intrusions and hornfels contact aureoles can also be magnetic  . 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 has the RTP applied and International Geomagnetic Reference Field (IGRF) already applied and thus reflects residual magnetic intensity level ranges from around −330.0 nT to over 280.0 nT with the shapes of the anomalies approximately centered over their respective causative bodies. The map shows rather complex crustal magnetization pattern with different locations of high and low magnetic intensities. It exhibits some different types of positive (magenta and red colours) and negative (light and dark blue colours) anomalies. Accordingly, the western and the northern parts of the map area are occupied by two major positive magnetic belts. This map elucidates two major positive magnetic belts trending E-W in the north and NW-SE occupied in the northwestern part (Mabule to Lerolwane). The Segwagwa Complex is represented by a large roughly oval (NNW-trending) highly magnetic feature with average width of about 10 km. Another feature of probably similar composition is seen around Gasita, though with irregular shape (Figure 4). The map also reveals prominent wide negative zones in elongated and irregular shapes distributed at the central parts of the mapped area. The major negative anomaly occurs at the central part (see Mmakgori, Metlobo and Mokgoriapitse villages). Based on the results of enhanced maps, it is quite likely that magnetic quiet areas are widely distributed in the central part of the study area and usually exhibit very low magnetic relief. These areas (Mmakgori to Metlobo) correspond to areas covered by a blanket of Tertiary-Recent Kalahari sediments  . These areas are faulted and folded as recognized both in the field and the geophysical maps. The magnetically moderate areas, such as those in the Papatlo-Ramatlabama and Metlojane villages have accentuated magnetic relief with lineaments and anomalies having amplitudes of up to 200 nT. Most of the magnetic anomalies in these areas are observed over plutonic rock exposures.
Structurally, the colour-shading accentuates several linear and curvilinear features that could be dykes, faults and/or intrusions clearly observed, even without any major wavelength filtering of the data. The Kraaipan greenstone belt contains
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.
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 (good examples around Mafikeng), cross-cut by N-S lineaments of the Kraaipan BIF in the centre of the map (Molopo). The northeastern corner is marked by three nearly NW-SE trending lineaments/features of moderate to high positive magnetic intensities centered around Kanye, Ranaka and Lekgolobotlo but converging around Digawana.
The high resolution aeromagnetic (HRAM) data (Figure 5) show that the southeastern area is marked by NE-SW-trending high-frequency, short wavelength positive magnetic intensities, corresponding to a Proterozoic dyke swarm and faults or other lineaments  . Based on this map, the Kraaipan BIF is clear
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).
for over 30 km from west of Molopo. All these structures are further elucidated in the processed maps and their hydrogeological significance discussed later.
5.2. Gravity Data and Geological Correlation
The gravity data used in this study are described in detail by  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  . 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  . 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.
6. 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   . 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.
6.1. 23/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    and  . A regional gravity structural cross-section (Figure 7) was derived/obtained using GM-Sys’ 23/4D-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/m3) were tested but these did not significantly change the shapes and depth extents of the units, as noted/observed by others   .
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.
Figure 7. Kraaipan greenstone belt 23/4D gravity model along profile A-A’ (see Figure 6). The numbers inside the model bodies are densities in kg m−3. Densities not shown are for sediments 2500 kg m−3, granitic pluton 2620 kg m−3 and granitic gneiss 2750 kg m−3. Vertical Exaggeration (V.E) = 11.0.
6.2. 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.
6.2.1. 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.
6.2.2. Depths from Euler Deconvolution of Magnetic Data
Structural Mapping and Depth Estimates
The Euler deconvolution technique was carried out on the magnetic data using
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  .
Geosoft program (Oasis montaj 8.2). Structural indices SI = 0, 1 and 2 and window size of 12 grid cells (3000 m) were used for interpreting contacts, faults and other curvilinear source types   . The Euler solution at different depths for SI = 0 is shown in Figure 9. The depths to the detected structures range between 0.0 m (outcropping) and just below 350.0 m, which is also expected for a window size of 400 m.
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.
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     . However, these enhanced merged data together with existing aeromagnetic data yield more essential geological information, particularly those regarding contexts of tectonism and
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   .
mineral/groundwater exploration. First, as suggested by  , 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.
7.1. 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
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.
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.
7.2. 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  . 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.
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.
The extreme northern area of partly exposed granitic-gneiss terrain (hard rock aquifer) corresponds to high yields. The central area depicts two conspicuous NS-trending linear belts of high yields over the main BIF while the southeast area is characterised by moderate to high yields. The west-central part mainly covered by the Kalahari sands (sedimentary aquifer) surprisingly shows the lowest yields (Figure 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.
7.3. 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)
Figure 12. Existing borehole yield classifications in relation to aeromagnetic first vertical derivative structures (D = low yield range: 1 - 5 m3/hr; o = medium/average yield range: 6 - 9 m3/hr, and à = high yield range: 10 - 160 m3/hr).
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 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 ~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).
Figure 13. Existing borehole yield classifications in relation to interpreted aeromagnetic structures for HRAM area (D = low yield range: 1 - 5 m3/hr; o = medium/average yield range: 6 - 9 m3/hr, and à = high yield range: 10 - 160 m3/hr). Clearly, boreholes with low yield range are the majority.
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    . 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    .
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 m3/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.
7.4. 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).
8. 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).
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 12 and Figure 14, Table 3). However, the ~E-W magnetic lineament structures appear to cut across the entire study area. The area is dissected by an ENE-to-EW-trending Proterozoic dyke swarm visible on the high resolution aeromagnetic data (Figure 5 and Figure 13), 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 (1 - 5 m3/hr) and 432 yielding 6 - 49 m3/hr and only 28 above 50 m3/hr (maximum ~160 m3/hr) and an average water strike of 64 m. Finally, all geophysical results 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.
2) Correlation of structural map with hydrogeological and drill hole information.
3) Very little correlation between interpreted hydrogeophysical features and the existing borehole locations.
4) Shallow fractured-karst aquifers à very good; but prone to pollution (e.g., pit latrines).
5) Fractured crystalline and metamorphic aquifers à good (moderate depth and cost).
6) Deep sedimentary aquifers à fair (expensive/costly; possibility of brackish water due to poor circulation)    .
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     and recharge issues-geochemistry/isotopes    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  .
Sincere thanks to Department of Geological Survey of Botswana and Department of Water Affaires for providing the existing geophysical data and borehole data, respectively. The work was financially supported by Office of Research and Development (ORD) of University of Botswana―Grant R0699. The reviewers’ critical comments and the editor’s suggestions are greatly appreciated as they improved the paper.
 Ramotoroko, C.D., Ranganai, R.T. and Nyabeze, P. (2016) Extension of the Archaean Madibe-Kraaipan Granite-Greenstone Terrane in Southeast Botswana: Constraints from Gravity and Magnetic Data. Journal of African Earth Sciences, 123, 39-56.
 Ranganai, R.T., King, J.G., Koosimile, D.I. and Ngwisanyi, T.H. (2003) Geophysical Methods for Groundwater Resources Assessment, Development and Pollution Mapping and Monitoring: Examples from Botswana and Zimbabwe. 4th WATERNET/WARFSA Annual Symposium, Gaborone, 15-17 October 2003.
 Araffa, S.A.S. (2013) Delineation of Groundwater Aquifer and Subsurface Structures on North Cairo, Egypt, using Integrated Interpretation of Magnetic, Gravity, Geoelectrical and Geochemical Data. Geophysical Journal International, 192, 94-112.
 Murty, B.V.S. and Raghavan, V.K. (2002) The Gravity Method in Groundwater Exploration in Crystalline Rocks: A Study in the Peninsular Region of Hyderabad, India. Hydrogeology Journal, 10, 307-321. https://doi.org/10.1007/s10040-001-0184-2
 Ranganai, R.T. and Ebinger, C.J. (2008) Aeromagnetic and LANDSAT TM Structural Interpretation for Identifying Regional Groundwater Exploration Targets, South-Central Zimbabwe Craton. Journal of Applied Geophysics, 65, 73-83.
 Monteiro Santos, F.A., Sultan, S.A., Represas, P. and El Sorady, A.L. (2006) Joint Inversion of Gravity and Geoelectrical Data for Groundwater and Structural Investigation: Application to the Northwestern Part of Sinai, Egypt. Geophysical Journal International, 165, 705-718.
 Binley, A., Hubbard, S.S., Huisman, J.A., Revil, A., Robinson, D.A., Singha, K. and Slater, L.D. (2015) The Emergence of Hydrogeophysics for Improved Understanding of Subsurface Processes over Multiple Scales. Water Resources Research, 51, 3837-3866.
 Anhaeusser, C.R. and Walraven, F. (1999) Episodic Granitoid Emplacement in the Western Kaapvaal Craton: Evidence from the Archaean Kraaipan Granite-Greenstone Terrain, South Africa. Journal of African Earth Sciences, 8, 289-309.
 Mapeo, R.B.M., Armstrong, R.A., Kampunzu, A.B. and Ramokate, L.V. (2004) SHRIMP U-Pb Zircon Ages of Granitoids from the Western Domain of the Kaapvaal Craton, Southeastern Botswana: Implications for Crustal Evolution. South African Journal of Geology, 107, 159-172.
 Poujol, M., Hirner, A.J., Armstrong, R.A. and Anhaeusser, C.A. (2008) U-Pb SHRIMP Data for the Madibe Greenstone Belt: Implications for Crustal Growth on the Western Margin of the Kaapvaal Craton, South Africa. South African Journal of Geology, 111, 67-78.
 Hammond, N.Q. and Moore, J.M. (2006) Archaean Lode Gold Mineralisation in Banded Iron Formation at the Kalahari Goldridge Deposit, Kraaipan Greenstone Belt, South Africa. Mineralium Deposita, 41, 483-503.
 Bagai, Z. and Odumetse, M.J. (2011) Mutiple Phases in the Mmathethe Granite in SE Botswana; Evidence for a Succession of Intrusive Episodes. 23rd Colloquium of African Geology CAG23, Johannesburg, 10-14 January 2011, Abstracts Volume, 22.
 Cobbing, J.E., Hobbs, P.J., Meyer, R. and Davies, J. (2008) A Critical Overview of Transboundary Aquifers Shared by South Africa. Hydrogeology Journal, 16, 1207-1214.
 Mapeo, R.B.M. and Wingate, M.T.D. (2009) SHRIMP U-Pb Zircon Age for the Segwagwa-Masoke Igneous Complex of Southeastern Botswana and Implications for the Deformation History of the Achaean to Proterozoic Transvaal Supergroup of Southern Africa. Journal of African Earth Sciences, 54, 97-106.
 Stettler, E.H., de Beer, J.H., Eberle, D., Ludden, J. And Mareschal, M. (1997) Geophysics and Deep Structures. In: de Wit, M.J. and Ashwal, L.D., Eds., Greenstone Belts, Oxford Monograph on Geology and Geophysics, Clarendon Press, Oxford, 567-580.
 Poujol, M., Anhaeusser, C.R. and Armstrong, R.A. (2002) Episodic Granitoid Emplacement in the Archaean Amalia-Kraaipan Terrane, South Africa: Confirmation from Single Zircon U-Pb Geochronology. Journal of African Earth Sciences, 35, 147-161.
 DWA (1991) Botswana National Water Master Plan Study, Vol. 5 Hydrogeology, Ministry of Minerals, Energy and Water Resources (MMWER), Gaborone DWA, 2006. Botswana National Water Master Plan. Final Report Volume 5 Hydrogeology, by SMEC in Association with EHES, Department of Water Affairs, Gaborone.
 Bakundukize, C., Mtoni, Y., Martens, K., Van Camp, M. and Walraevens, K. (2016) Poor Understanding of the Hydrogeological Structure Is a Main Cause of Hand-Dug Wells Failure in Developing Countries: A Case Study of a Precambrian Basement Aquifer in Bugesera Region (Burundi). Journal of African Earth Sciences, 121, 180-199.
 Dietvorst, E.J.L., De Vries, J.J. and Gieske, A. (1991) Coincidence of Well Fields and Tectonic Basins in the Precambrian Shield Area of Southeast Botswana. Groundwater, 29, 869-877.
 Terblanche, M. and Stroebel, L. (2013) Case Study: Groundwater Exploration in the Basement Granites of the Kraaipan Group using Aeromagnetic Surveys. Groundwater: A New Paradigm. 13th Biennial Groundwater Division Conference & Exhibition, Durban, 17-19 September 2013. http://gwd.org.za/sites/gwd.org.za/files/
 Gibson, R.I. and Millegan, P.S. (1998) Geologic Applications of Gravity and Magnetics: Case Histories. Geophysical Reference Series 8, 162, Society of Exploration Geophysicists, Tulsa.
 Hinze, W.J., von Frese, R.R.B. and Saad, A.H. (2013) Gravity and Magnetic Exploration: Principles, Practices, and Applications. Cambridge University Press, 512 p.
 Bierlein, F.P., Murphy, F.C., Weinberg, R.F. and Lees, T. (2006) Distribution of Orogenic Gold Deposits in Relation to Fault Zones and Gravity Gradients: Targeting Tools Applied to the Eastern Goldfields, Yilgarn Craton, Western Australia. Mineralium Deposita, 41, 107-126.
 Reeves, C.V., Reford, S.W. and Milligan, P.R. (1997) Airborne Geophysics-Old Methods, New Images. In: Gubbins, A.G., Ed., Proceedings of Exploration: 4th Decennial International Conference on Mineral Exploration, 13-30.
 Broome, H.J. (1990) Generation and Interpretation of Geophysical Images with Examples from the Rae Province, Northwestern Canada Shield. Geophysics, 55, 977-997.
 MacLeod, I.N., Jones, K. and Dai, T.F. (1993) 3-D Analytic Signal in the Interpretation of Total Magnetic Field Data at Low Magnetic Latitudes. Exploration Geophysics, 24, 679-688.
 Reid, A.B., Allsop, J.M., Granser, H., Millet, A.J. and Somerton, I.W. (1990) Magnetic Interpretation in Three Dimensions using Euler Deconvolution. Geophysics, 55, 80-91.
 Reid, A.B., Ebbing, J. and Webb, S.J. (2014) Avoidable Euler Errors—The Use and Abuse of Euler Deconvolution Applied to Potential Fields. Geophysical Prospecting, 62, 1162-1168.
 Gwavava, O. and Ranganai, R.T. (2009) The Geology and Structure of the Masvingo Greenstone Belt and Adjacent Granite Plutons from Geophysical Data, Zimbabwe Craton. South African Journal of Geology, 112, 119-132.
 Levin, M., Gat, J.R. and Issar, A. (1980) Precipitation, Flood and Groundwaters of the Negev Highlands: An Isotopic Study of Desert Hydrology. In: Proceedings Advisory Group Meeting, Arid Zone Hydrology: Investigations with Isotope Techniques, International Atomic Energy Agency, IAEA-AG-158/1, 3-22, Vienna.
 Soupios, P.M., Kouli, M., Vallianatos, F., Vafidis, A. and Stavroulakis, G. (2007) Estimation of Aquifer Hydraulic Parameters from Surficial Geophysical Methods: A Case Study of Keritis Basin in Chania (Crete-Greece). Journal of Hydrology, 338, 122-131.