Advances in Remote Sensing
Vol. 2  No. 2 (2013) , Article ID: 33187 , 13 pages DOI:10.4236/ars.2013.22019

Application of Remote Sensing and Geographical Information Systems in Determining the Groundwater Potential in the Crystalline Basement of Bulawayo Metropolitan Area, Zimbabwe

Constant Chuma1,2*, Oluwagbenga O. I. Orimoogunje3, Dumisani J. Hlatywayo2, Joseph O. Akinyede1

1African Regional Centre for Space Science and Technology Education in English, Obafemi Awolowo University, Ile-Ife, Nigeria

2Department of Applied Physics, Faculty of Applied Sciences, National University of Science and Technology, Bulawayo, Zimbabwe

3Department of Geography, Faculty of Social Sciences, Obafemi Awolowo University, Ile-Ife, Nigeria

Email: *chumac80@gmail.com, orimoogunje@yahoo.com, djhlatywayo@gmail.com, director@arcsstee.org

Copyright © 2013 Constant Chuma et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received February 19, 2013; revised March 21, 2013; accepted March 30, 2013

Keywords: Groundwater Potential; Lineaments; Remote Sensing; GIS; Crystalline Basement; Aquifer; Bulawayo Metropolitan

ABSTRACT

This study is concerned with the identification and delineation of aquiferous zones for potential groundwater development across Bulawayo Metropolitan from remotely sensed data and geological inference. Attempts have been made to review literatures on groundwater exploitation in the study area and the constraints to effective and sustainable management of underground water in the study area. Remote Sensing and Geographical Information Systems (GIS) is useful in the recognition and delineation of aquiferous zones for potential groundwater in crystalline basement aquifers. LANDSAT ETM+ image, SRTM data, aeromagnetic data and other ancillary data sets were utilized to extract information on the groundwater storing controlling features of this study area. Six thematic maps were produced from remote sensing data and other ancillary data—Land use/landcover, drainage density, slope map, contact density, lithology and lineament density. GIS modeling technique of the index overlay method was used to produce the groundwater potential map. The study revealed that the regional lineaments correlate with faults, fracture zones, and lithological contrasts along fold belts in the crystalline basement rocks, while the main direction of faulting and jointing is north-north-west to north with several faults oriented, to the north-north-east, parallel to the Great Dyke. Proximity to lineaments is the highest zone of increased porosity and permeability which in turn have a greater chance of accumulating groundwater. The results have shown massive spatial variability of ground water potential ranging from very good to poor. The variability closely followed variations in the structures, geology, topography/slope, drainage density and land use/land cover in the project area. This work is a reconnaissance which needs to be validated by the use of high-resolution terrain data and satellite imagery and the quantitative analysis should be done using geophysical and hydrogeological surveys.

1. Introduction

Groundwater plays a crucial role in arid and semi-arid regions in most parts of the world especially in Africa. However, in Africa most arid and semi-arid regions are underlain by crystalline basement rocks and groundwater development is complicated by the high variability of geological and hydrogeological settings of the basement aquifers [1-3]. The situation is not good for Bulawayo Metropolitan which is covered by crystalline complex rock and frequent water shortages are being faced annually due to unreliable surface water supplies [4].

There has been comparatively little research conducted on the occurrence of groundwater in the whole crystalline basement aquifer of the Bulawayo Metropolitan, despite its significance in mitigating problems of water shortages in semi-arid regions. Most of the research conducted on part of the area has focused on groundwater quality monitoring, groundwater resource evaluation and partial characterisation of the subsurface e.g. [4-11]. Consequently, there is insufficient information to explain the occurrence and distribution of groundwater in the Bulawayo basement complex aquifer. This has slowed down the exploitation and development of groundwater in the Bulawayo City as a complimentary water source to the surface water. In [2], it was also reported that the lack of understanding of the occurrence, movement and recharge of groundwater in basement terrain has frequently contributed to the unsustainable use of this resource.

Through the use remote sensing and GIS, suitable   areas for recharging basement aquifers can be delineated which have porous lithologies, maximum fractures, highly weathered region and regions of null slope. The exploration of surface features by integrating the remote sensing data and GIS techniques proved to be an essential in identifying aquiferous zones for potential groundwater exploitation in Bulawayo crystalline basement rock. This study provided preliminary information which will enhance the success of the quantitative subsurface geophysical investigation that will subsequently follow the study. The regional lineament map produced in this study is not only essential to groundwater but will be used as a reconnaissance study in mineral exploration, as well as in identifying earthquake potential zones. Zimbabwe has been facing tectonic movements in most parts of the country including the western part of Bulawayo recent that are being investigated [12,13].

The objectives of this study, therefore were to produce a regional structural lineament map of Bulawayo Metropolitan from remotely sensed data, topographic, geological inference and structural data; prepare thematic maps of the area such as lithology, lineaments, drainage density, slope and land use/land cover from remotely sensed data and other data sources; assess the hydrogeological implication of the lineaments by integrating them with the available ancillary data especially digital elevation models (DEM), and finally identify and delineate groundwater potential zones through integration of various thematic maps within a GIS environment.

2. Description of Study Area

The study area lies in the water-distressed, western region of Zimbabwe. It is located between longitude 28˚22'05"S and 28˚46'12"S and between latitude 20˚18'04"E and 20˚00'00"E. The gross command area is approximately 400 km2 in areal extent. Bulawayo is the second largest city in Zimbabwe with a population of over 1 million people, and is regarded as the industrial and business capital of Zimbabwe and is home to the National Railways of Zimbabwe because of its strategic position near Botswana and South Africa. It is the nearest large city to most tourist resort areas in Zimbabwe such as Hwange National Park, Matopo National Park and Victoria Falls. Figure 1 shows the study area and inserts of the map of Africa and the climate is semi-arid with dry winters and wet summers. The mean annual rainfall is 550 mm [14]. Rainfall is not evenly distributed throughout the year and droughts are a frequent phenomenon. Mean annual maximum and minimum temperatures for the study area is 25.8˚C and 12.7˚C, respectively.

3. Materials and Methods

3.1. Geology of the Study Area

The geology of Bulawayo is mainly constituted by the crystalline basement rock. The crystalline basement encompasses the Bulawayo greenstone belt and the surrounding granite rock. According to the geological report by [15] the greenstone belt is flanked by a belt of phyroblastic tonalitic to granodioritic gneisses, penetrated by later granitic bodies. The greenstones comprise a succession of tightly-folded metavolcanic and metasedimentary rocks with generally steep dips away from the margins of the granitic rocks. A second phase of deformation is associated with thrusting directed from the eastern part of the area over a western hinterland consisting of granitic gneiss and lower greenstone rocks. Superimposed in this phase is an episode in which a strong east-west cleavage was developed. Faulting and shearing later occurred mainly along north to north-north-west trends with some off-setting boundaries. The general lithology underlay the study area shown in Figure 2.

3.2. Preparing the Thematic Maps

This study was mainly carried out using Remote Sensing (RS) and Geographical Information Systems (GIS). In addition to that, an extensive desk study was done in collection, collation and analyzing aeromagnetic data. Remote sensing have the advantages of spatial, spectral and temporal availability of data covering large and inaccessible areas within a short time and has become a very rapid and cost effective tool in assessing, monitoring and conserving groundwater resources [16]. On the other hand, GIS is powerful environment for real time database development, especially in studies such as delineating groundwater pollution potential zones [17-19] and recharge sites [20], groundwater modeling studies [21,22] among others. The data for this study were obtained mainly from geologic maps, topographic maps, LANDSAT 7 enhanced thematic mapper plus (ETM+) image of October 2005, ASTER DEM, aeromagnetic data, geologic reports and climatological records coupled with intensive fieldwork.

Figure 1. Study area showing the roads, dams and inserts of Africa and Zimbabwe.

The methodology adopted in this study is simplified by the schematic flow chart in Figure 3. There were four correlated processing pathways followed: the LANDSAT ETM+ data, ASTER DEM data, aeromagnetic data and geological data. Remote sensing data were processed separately from the other data from which land-cover classes were determined through image classification. Additionally, lineaments in the study area were identified visually. Ancillary data were further processed to provide the layers, and classes for each layer needed for the ArcGIS Model.

The lineaments were identified by visual interpretation of both LANDSAT ETM+ and ASTER DEM processed images (Figure 4). In order to clearly identify the lineaments from LANDSAT ETM+ data, a false color composite band combination of 754 (RGB) was used which was preceded with principal component analysis and directional filtering using Environment for Visualizing Image software (ENVI 4.7). The lineaments observed on the imagery are shadowed with radiation from the sun, and were traced onto the base map. Since the area was built up there, very few lineaments were identified from

Figure 2. Map showing the geological pattern in the study area.

Figure 3. The methodology used to determine the groundwater potential zones in the study area.

Figure 4. The lineaments observed from the LANDSAT EMT (right) and ASTER DEM data (left).

the LANDSAT ETM+ image. The ASTER DEM data were converted to a hillside terrain map using ArcGIS (ArcMap). From the hillshade surface lineaments were observed by a drop in elevation for a short distance.

Drainage pattern is one of the most important indicators of hydrogeological features, because drainage pattern and density are controlled in a fundamental way by the underlying lithology [23]. The drainage map was prepared by hand digitizing the drainage lines from processed LANDSAT ETM+ and geological map. Drainage density is the length of streams per unit area (Equation (1)). It is usually measured in km/km2.

(1)

where Dd is the drainage density, Ls is the total length of streams, and As is the total area of the whole basin.

The slope amount map was prepared using contours produced from ASTER DEM 90m data and topographic map. Surfer 8 was then used to produce a three dimension surface map.

Land use and land cover (LULC) influences the occurrence of sub-surface groundwater occurrence, as a result supervised classification using the maximum likelihood classification was performed on the LANDSAT EMT+ image using ERDAS Imagine version 9.2. A false color combination of 432 (RGB) was used to identify the different themes in the image.

3.3. Groundwater Potential Modelling

In order to delineate potential groundwater sites in the project area, all the data sets were integrated using the model constructed in ArcGIS Model builder engine [24]. The final map produced by Weighted Linear Combination (WLC) where each class individual’s weight was multiplied by the map scores and then adds to the results as in Equation (2).

(2)

where S is the suitability, Wi is the weight for each map score and Xi is the individual map.

3.4. Weighting of Thematic Maps

Finally the six thematic layers were converted into a grid with related item weight generated from the Analytic Hierarchy Process (AHP) by [25,26]. The final six thematic maps were reclassified in the ArcMap according to the importance of each factor in each raster image as follows:

• Reclassification of land use/cover for analysis was done based on their character to infiltrate water into the ground and to hold water on the ground.

• Drainage density raster was reclassified based on the recharge potential.

• Pairwise comparison done and the weight calculated for slope angle was based on the fact that the flatter the topography (low slope angle) is the better are the chances for groundwater accumulation.

• Reclassified map of lineament was produced based on the weight calculated after a pair-wise comparison done based on the fact that areas closer to lineaments are the highest zone of increased porosity and permeability which in turn have greater chance of accumulating groundwater.

• Reclassification of the lithologies weight was based on the susceptibility of the rock to weathering.

• Contacts density map was based on the nature of the rocks in contact and the density of the contacts. Contacts normally act as conduit of groundwater movement in basement rocks.

In the AHP method the relative importance of each individual class within the same map and vector maps were compared to each other. As a result important matrices were produced using calculated weights from the IDRSIS module (Table 1). MatLab software was used to solve the matrices and calculate the consistence value. The grids in the integrated layer were grouped into different groundwater potential zones by a suitable logical reasoning and conditioning with emphasis on the hydrogeological condition of each grid cell. The layers were then overlaid using a weighted sum in the ArcMap module. Figure 5 illustrates the overlaying of the six factors that affects the occurrence of groundwater in the ArcSene environment.

The different polygons in the thematic layers were labeled separately and then registered. In the final thematic layer initially each of the polygons was designated on the basis of interpretation of groundwater occurrence and weights in terms of interpreting the probability level of occurrence that was assigned to each polygon.

4. Results and Discussion

Geospatial analysis of the aeromagnetic data was useful in delineating the extent of the Bulawayo crystalline basement rock as shown in Figure 6. The high magnetic susceptibility values observed indicate the presence of igneous rocks which are mainly composed of basalt and granite. These two rock formations gave what is known as the Bulawayo Crystalline Basement mentioned in [27]. Therefore, the study area is mainly underlaid by a basement aquifer. Groundwater is mainly controlled by secondary porosity, which is due to mechanical and chemical weathering.

The high magnetic values are associated with serpentinised volcanic activity whilst the sills and other intrusions that developed over the years. Crystalline basement rocks (i.e. usually of Precambrian age) have generally been subjected to multiple tectonic events under varying stress conditions that resulted in complex patterns of ductile folding and brittle fracturing in the near-surface regions of the earth’s crust [28]. According to [1] basement aquifers are essentially phreatic in character, but may respond to localized abstraction in a semi-confined fashion, if the rest water level occurs in a low permeability horizon (e.g., in a clayey regolith). The greenstone belt in Bulawayo forms the crystalline rock formation which favors the occurrence of the basement aquifer in the study area.

4.1. Regional Structural Lineaments Analysis

Lineaments being a surface manifestation of structurally controlled features such as faults, fractures and rock contacts their high density may present incipient highly connected fractures that are favorable for the accumulation of groundwater. Therefore, lineament analysis of the study area extracted from both the remotely sensed data and geological images gave important information on subsurface features that controls the movement and or storage of groundwater (Figure 7). In addition to that, lineaments like joints, fractures etc., developing generally due to tectonic stress and strain, provide important clue on surface features and are responsible for the infiltration of surface runoff into the subsurface and also for the movement and storage of groundwater [29]. Most topographic lineaments in the study area correlate with faults, fracture zones, and lithological contrasts along fold belts and in crystalline basement rocks.

A detailed structure of the Bulawayo Greenstone Belt evidently is elucidated within the overall setting and framework of deformation of the Zimbabwean craton (Figure 8). Dolerite dykes; quartz and pegmatite veins traverse the area both in concordant and discordant fashions, which represent a later phase of basic igneous activity. All the rock types except the younger intrusive have suffered several phases of deformation and have been metamorphosed. These deep-seated structures and tectonic features include, for example, buried uplifts, buried igneous plugs and salt diapirs, subsurface folds and faults (both thrust and normal faults), strata affected

Table 1. Weight applied in ArcMap weight sum scheme.

Figure 5. Integration of the thematic maps which gives the potential zones.

Figure 6. The geospatial analysis of aeromagnetic data indicating Bulawayo crystalline basement rocks.

by subtle regional warping, and stratigraphic pinch-outs [30].

From the regional structural lineament map a rose diagram was produced using Integrated Land and Water Information System 3.3 (ILWIS) software which shows the direction and magnitude of lineaments (Figure 8). Fracture zones usually occur along lineaments and often correspond to the surface drainage patterns. Fractures, dykes and zones of weathering and/or overburden deposits are targets for the exploitation of groundwater.

The main direction of the lineaments is north-westerly with more west-north-west trends in the west of the study area. The ENE-WSW to NE-SW lineaments are largely due normal faults with possibly of a small amount of transcurrent movement. As a result most of the rivers flow in the western direction of the study area as shown by the drainage map overlaid on the DEM in Figure 9. Bulawayo has four major rivers; Umatseumhlope River which runs from south to north, Khami River which runs east to west, Mpopoma River and Umguza River which run from south-east to north-west. All the small river and large rivers, which are found in the project area drain from central part to west and north direction.

4.2. Analysis of Factors Affecting the Occurrence of Groundwater

All surface and subsurface factors that affect the occurrence of groundwater were acquired, processed and analyzed. The reclassified thematic maps of all the six parameters which were incooperated in the APH system are shown in Figure 10.

Comparison of the drainage system of the area and structure has shown that the drainage system of the area is structurally controlled following lineaments directions. Dendritic and parallel drainage pattern are recognized, which are indicative of the presence of structures that act as conduits or storage for sub-surface water. Structurally controlled drainage patterns are observed in western and eastern part of the project area. Dendritic are formed on impervious rock whereas parallel drainage system form on uniform gentle slope. Overally, the drainage pattern in Bulawayo shows that the area very high drainage density is found in the west and north parts of the study area. Moderate and low drainage density concentrates in the southern and central part of the project area. The distribution is mainly attributed to the weak zones that are created due faulting, shearing and geological contacts.

Evaluation of LANDSAT ETM+ image and topographic map of year 2005 has shown that there is remarkable expansion of settlements since the independence of Zimbabwe in 1980, which negatively affect the groundwater recharge of the area. The original vegetation is more preserved in area far from the City Centre. This is due to the fact that people had made the natural forest to disappear on the flat lands near to the main road and cleared the land for construction. Generally settlements

Figure 7. Geological structures, extracted lineaments and contacts.

Figure 8. Rose diagram showing directions and lengths of lineaments in Bulawayo.

are found to be the least suitable for infiltration and after pair-wise comparison of each class weight for each class was calculated.

Based on the susceptibility of the rock to weathering, about nine classes were observed which shows have an intensive analysis of the crystalline basement rock. Results have shown that basalt rocks are easily weathered compared to the surrounding granite. Basalt rocks are associated with presence of structures such fractures, fissures, veins, joint and other structural deformations that occur in basement complex which controls the flow of ground water and also influence the rate of recharge and

Figure 9. Drainage pattern overlaid on the DEM.

discharge of the main aquiferous units. This agrees with the aeromagnetic imaging of the area which shows a contact between the granite rocks and the metabasalt rocks

Figure 10. Reclassified thematic maps of the study indicating the land use and densities of slope, drainage, lineaments and contacts.

which are penetrated intrusions in the south west part of the study area.

Figure 11 portrays a three dimensional surface map of the study area that was produced from the ASTER DEM data which was processed and analyzed in the Surfer 8 software. The surface map is showing an abrupt change

Figure 11. The three dimensional analysis of the study area.

in terrain which explains the variability of the slope in all directions. This is also an indication of the possible direction of geological structures in the study which control the movement and storage of groundwater. The slope was generated from the (Triangulated Irregular Network) TIN in the ArcMap environment. Slope map layer indicates the relative time available for runoff via land to infiltrate. The region is mainly hilly with an extensive plateau to the south and small plateaus to the north and east. The valley itself is quite narrow and elongates, with steep ridge slopes that grade into a portion of terraces that, in turn, gently slope towards the main river channel. Topography of the land surface determines the general direction of groundwater flow, and it influences groundwater recharge and discharge.

4.3. Occurrence of Groundwater in Bulawayo

The map produced from the integration of the six thematic maps was reclassified and five categories of groundwater potential zones were produced as shown in Figure 12. The amount of groundwater is mainly determined by the structural geology environment and the topographic effect. The spatial distributions of the various groundwater potential zones obtained from the model generally show regional patterns of lineaments, drainage, landform and lithology.

Spatially, very good and good categories are distributed along areas near to lineaments and less drainage density and where the lithology is affected by secondary structure and having interconnected pore spaces. This highlights importance of lineaments, geology and hydro-geomorphological parameters in the project area.  Areas with moderate groundwater prospects are attributed to contributions from combinations of the land use/cover, lithology, slope and landform. The low to poor categories of groundwater potential zones are spatially distributed mainly along ridges where slope class is very high, the lithology is compact/massive and far from lineaments.

The greenstone belt in Bulawayo forms the crystalline rock formation which favors the occurrence of the basement aquifer in the study area. Groundwater is mainly controlled by secondary porosity, which is due to mechanical and chemical weathering. The fractures and faults that are produced during the course of these processes act as conduit for groundwater movement.

Figure 13 shows the pie chart of the classified categories obtained from groundwater potential model. About 36% of the study area has high potential of groundwater and is mainly concentrated in the southern and southerneastern part of the study area. There are also high groundwater potential zone in the northern direction, however the sources area likely to be high vulnerable to sewage system contamination. A great correlation can be noticed between areas assigned as high potential areas and the location of the Matsheumhlope aquifer [4].

Areas with moderate groundwater prospects are attributed to contributions from combinations of the land use/cover, lithology, slope and landform. The low to poor categories of groundwater potential zones are spatially distributed mainly along ridges where slope class is very high, the lithology is compact/massive and far from

Figure 12. Groundwater potential zones identified from the weighted factors.

Figure 13. Pie chart of the classified categories obtained from ground-water potential model.

lineaments.

The variability in flow direction maybe attributed to variation and localization of the fractures within the basement rock. This is one of the characteristics of the basement rock which makes it difficult to exploit the groundwater available especially if the information is inadequate. Overally, the direction of flow is north and north-west as shown in Figure 14. This is mainly attributed to the drainage pattern in the study area.

5. Conclusions

The picture that emerges from this study has shown that

Figure 14. Possible flow direction of groundwater in the study area.

there is large spatial variability of ground water potential while the most promising potential zone in the area is related to volcanic rock of which is affected, by secondary structure and having interconnected pore spaces, with gentle slope and less drainage density. It can also be inferred that most of the zones with low to poor groundwater potential lie in the massive basements unit which is far from lineaments. Areas close the sewer drainage system should be avoided or buffered to prevent contaminated water.

It can also be concluded that parts of the area that are characterized with surface expression of lineament are considered to be hydrogeological insignificant due to high slope and hence constitute the poor groundwater potential zone for the area. The low drainage density areas cause more infiltration and result in good groundwater potential zones as compared to high drainage density areas. High drainage density indicates unfavourable site for groundwater occurrence especially in area the northwestern part of Bulawayo Urban.

It must, however, be added that the results of this study is qualitative, hence, it is recommended that more quantitative investigations involving geophysical technique(s), borehole drilling and pump test must be carried out to fully understand the prospect of the groundwater resources of the area and consequently design a model for the effective groundwater use and management in the area.

6. Acknowledgements

I want to extend my profound gratitude to African Regional Centre for Space Science and Technology in English (ARCSSTE-E) for according me an opportunity to study and the financial support you gave throughout this course, National University of Science and Technology (NUST) and Zimbabwe Government for the support they rendered to the success of this study. Special thanks go to Dr. O. O. Orimoogunje, the supervisor for his expert advice in remote sensing and GIS and guidance at every stage of the study. Dr. D. J. Hlatywayo, my co-supervisor, requires special mention for his untiring supervision and guidance during field investigations. All the lecturers and ARCSSTE-E staff members including the Director, Dr. J.O Akenyede H.O.D and Administrative Course Coordinator I salute you all for imparting precious knowledge in my life.

REFERENCES

  1. E. P. Wright, “The Hydrogeology of Crystalline Basement Aquifers in Africa,” In: E. P. Wright and W. G. Burgess, Eds., The Hydrogeology of Crystalline Basement Aquifers in Africa, Geological Society Special Publication, London, 1992, pp. 1-27.
  2. P. J. Chilton and S. S. D. Foster, “Hydrogeological Characterisation and Water-Supply Potential of Basement Aquifers in Tropical Africa,” Hydrologeology Journal, Vol. 3, No. 1, 1995, pp. 36-49. doi:10.1007/s100400050061
  3. D. Banks and N. S. Robins, “An Introduction to Groundwater in Crystalline Bedrock,” Geological Survey of Norway, Trondheim, 2002.
  4. F. Rusinga and A. E. Taigbenu, “Groundwater Resource Evaluation of Urban Bulawayo Aquifer,” Water SA, Vol. 31, No. 1, 2004, pp. 23-34.
  5. E. Martinelli and G. L. Hubert, “Hydrogeology. A Report to the Ministry of Energy and Water Resources and Development, National Master Plan for Rural Water Supply and Sanitation,” Interconsult/Norad, Vol. 2, 1985, pp. 109-114.
  6. J. M. C. Weaver, J. Conrad and J. F. P. Engelbrecht, “The Matsheumhlope Wellfield Project, Wellfield Evaluation after Phase 1 Drilling for Bulawayo City Water Supply,” Matabeleland Chamber of Industries, Bulawayo, 1992.
  7. E. Mangore and A. E. Taigbenu, “Land-Use Impacts on the Quality of Groundwater in Bulawayo,” Water SA, Vol. 30, No. 4, 2004, pp. 453-464 doi:10.4314/wsa.v30i4.5097
  8. E. Mangore, “Monitoring and Evaluation of Bulawayo Aquifers,” Unpublished M.Phil. Thesis, Library of National University of Science and Technology, Bulawayo, 2002.
  9. I. Muchingami, J. D. Hlatywayo and J. M. Nel, “Electrical Resistivity Survey for Groundwater Investigations and Shallow Subsurface Evaluation of the Basaltic-Greenstone Formation of the Urban Bulawayo Aquifer,” 11th WaterNet/WARFSA/GWP-SA Symposium, Victoria Falls, 11 November 2010, pp. 279-291.
  10. C. Chuma, “Investigating the Occurrence of Groundwater in the Low Yield Aquifer within Matsheumhlope Wellfield,” Unpublished M.Sc. Thesis, National University of Science and Technology Library, Bulawayo, 2011.
  11. I. Muchingami, J. D. Hlatywayo, J. M. Nel and C. Chuma, “Electrical Resistivity Survey for Groundwater Investigations and Shallow Subsurface Evaluation of the BasalticGreenstone Formation of the Urban Bulawayo Aquifer,” Journal for Physics and Chemistry of the Earth, Vol. 50- 52, 2012, pp. 44-51. doi:10.1016/j.pce.2012.08.014
  12. B. T. Shumba, D. J. Hlatywayo and V. Midzi, “Focal Mechanism Solution of The 15th March 2008, Nyamandlovu Earthquake,” South African Journal of Geology, Vol. 112, No. 3, 2009, pp. 381-386. doi:10.2113/gssajg.112.2-3.381
  13. D. J. Hlatywayo and V. Midzi, “Report on the Investigations into the Location, Causes and the Effects of the 25th June 2004 Bulawayo-Nyamandlovu Earthquake,” Goverment of Zimbabwe, Harare, 2005, 13p.
  14. Department of Meteorological Services, “Climate Data Records of Bulawayo,” Goetz Observatory Weather Station, Bulawayo, 2003.
  15. M. S. Garson and N. A. Mtsvangwa, “The Geology of Bulawayo Greenstone Belt and the Surrounding Granitic Terrain,” Zimbabwe Geological Survey, Harare, 1995.
  16. K. Sankar, “Evaluation of Groundwater Potential Zones Using Remote Sensing Data in Upper Vaigai River Basin, Tamil Nadu, India,” Journal of the Indian Society of Remote Sensing, Vol. 30, No. 3, 2002, pp. 119-129. doi:10.1007/BF02990644
  17. B. M. Evans and W. L. Myers, “A GIS-based approach to evaluating regional groundwater pollution potential with DRASTIC,” Journal of Soil and Water Conservation, Vol. 45, No. 2, 1990 pp. 242-245.
  18. J. W. Merchant, “GIS-Based Groundwater Pollution Hazard Assessment: A Critical Review of the DRASTIC Model,” Photogrammetric Engineering and Remote Sensing, Vol. 60, No. 9, 1994, pp. 1117-1128.
  19. G. P. Panagopoulos, A. K. Antonakos and N. J. Lambrakis, “Optimization of the DRASTIC Method for Groundwater Vulnerability Assessment via the Use of Simple Statistical Methods and GIS,” Hydrogeology Journal, Vol. 14, No. 6, 2006, pp. 894-911. doi:10.1007/s10040-005-0008-x
  20. T. O. Tweed, M. Leblanc, J. A. Webb and M. W. Lubczynski, “Remote Sensing and GIS for Mapping Groundwater Recharge and Discharge Areas in Salinity Prone Catchments, Southeastern Australia,” Hydrogeology Journal, Vol. 15, No. 1, 2007, pp. 75-96. doi:10.1007/s10040-006-0129-x
  21. E. Sener, A. Davraz and M. Ozcelik, “An Integration of GIS and Remote Sensing in Groundwater Investigations: A Case Study in Burdur, Turkey,” Hydrogeology Journal, Vol. 13, No. 5-6, 2005, pp. 826-834. doi:10.1007/s10040-004-0378-5
  22. S. Shahid, S. K. Nath and J. Roy, “Groundwater Potential Modelling in a Soft Rock Area Using a GIS,” International Journal of Remote Sensing, Vol. 21, No. 9, 2000, pp. 1919-1924. doi:10.1080/014311600209823
  23. J. E. Charon, “Hydrogeological Applications of ERTS Satellite Imagery,” Proceedings of UN/FAO Regional Seminar on Remote Sensing of Earth Resources and Environment, Commonwealth Science Council, Cairo, 1974, pp. 439-456.
  24. G. F. Bonham-Carter, “Geographic Information Systems for geoscientists: Modelling with GIS,” Elseviers Science Ltd., Kidlington, 1994, 398p.
  25. T. L. Saaty, “A Scaling Method for Priorities in Hierarchical Structures,” Journal of Mathematical Psychology, Vol. 15, No. 3, 1997, pp. 234-281. doi:10.1016/0022-2496(77)90033-5
  26. T. L. Saaty, “Fundamentals of Decision Making and Priority Theory with the AHP,” RWS Publications, Pittsburgh, 1994.
  27. J. E. Danielsen, T. Dahlin, R. Owen, P. Mangeya and E. Auken, “Geophysical and Hydrogeologic Investigation of Groundwater in the Karoo Stratigraphic Sequence at Sawmills in Northern Matabeleland, Zimbabwe: A Case History,” Hydrogeology Journal, Vol. 15, No. 5, 2007, pp. 945-960. doi:10.1007/s10040-007-0191-z
  28. J. W. Lloyd, “Water Resources in Hard Rock Aquifers in Arid and Semi-Arid Areas,” UNESCO Studies and Reports in Hydrogeology No. 58, United Nations Educational, Scientific and Cultural Organization (UNESCO), Paris, 1999
  29. N. S. Rao, G. K. J. Chakradhar and V. Srinivas, “Identification of Groundwater Potential Zones Using Remote Sensing Techniques in and around Gunur Town, Andhra Pradesh, India,” Journal of the Indian Society of Remote Sensing, Vol. 29, No. 1-2, 2001, pp. 69-78. doi:10.1007/BF02989916
  30. S. C. Caran, C. M. Woodruff and E. J. Thompson, “Lineament Analysis and Inference of Geologic Structure Examples from the Balcones/Ouachita Trend of Texas,” Gulf Coast Association of Geological Societies Transactions, Vol. 31, 1981, pp. 59-69.

NOTES

*Corresponding author.