Groundwater constitutes a very important natural resource for mankind. It occurs almost everywhere beneath the surface of the earth within saturated zones where the hydro-static pressure is equal to or greater than atmospheric pressure. Most of groundwater comes from precipitation and is stored in, and moves slowly through layers of soil, sand and rocks called aquifers, which may either be igneous, sedimentary or metamorphic. Despite the fact that Cameroon is endowed with abundant water resources, the quest for suitable and viable portable water sources remains a challenge for different water stakeholders. Thus as the need for portable water increases, great care should be taken to ensure sustainability of the water resources with respect to quantity and quality. Challenges of water resources exploitation combined with an increasing population, therefore, presents greater challenges to stakeholders at all levels   .
Mount Cameroon whose geographic occurrence is in Buea, in the South West region of Cameroon, is a 4095 m high and 50 km wide volcano along the Cameroon Volcanic Line. This mountain is predominantly made up of volcanic materials (basalts, lavas and pyroclastics) and containing several faults and fractures  . The coarse nature of the pyroclastics and the presence of fractured lavas make Mount Cameroon  potential volcanic aquifer system. The potential of volcanic rocks as aquifers strongly depends on their hydraulic characteristics, including the hydraulic conductivity, transmissivity and storativity. Hydraulic conductivity values range from orders of 10−8 (unfractured igneous and metamorphic rocks) to 104 (lava flows and gravel). Transmissivity is the capacity of an aquifer to transmit water of the prevailing kinematic viscosity. It is the product of the hydraulic conductivity and the aquifer’s thickness  . Determination of hydraulic characteristics of an aquifer could be done by carrying out geo-electric and borehole investigations. A combination of both geo-electrical and borehole investigations provide optimum information on hydraulic characteristics of an aquifer  . Pumping tests are the only borehole investigations that provide simultaneous information on the hydraulic behaviour of the well (borehole), the reservoir and the reservoir boundaries, which are essential for efficient aquifer and well field management  . Pumping tests areused for two purposes: 1) to determine the performance characteristics of a well and 2) to determine the hydraulic characteristics of an aquifer. A geo-electric investigation could involve both vertical electric sounding (VES) and constant separation traversing (CST). VES is used for depth sounding to determine the vertical variation of resistivity, while CST is used for horizontal variation of resistivity  . VES and CST allow an electrical image of the subsurface to be created that can be used to evaluate groundwater distribution in the aquifer  . These methods function as complementary tools, reduce uncertainties in subsurface exploration and assist in determining the most reliable location for groundwater exploitation  .
Shortage of potable water for human consumption as well as industrial needs has grown and continues to grow across the country. The supply by the national water company in metropolis is grossly inadequate and most people depend on surface water and water from hand dug wells, while the most privileged are supplied by water from shallow or deep boreholes. While ground water is easily obtainable in sedimentary areas, the situation is more complex in bedrock areas. The location of possible aquifers in such areas could be possible through geo-electrical and borehole investigations. Currently in Cameroon, there is an observable increase in demand for bottled water as people believe such water may have a better level of purification than that provided by the local portable water providers. More so, bottled water can get to remote areas easily than establishing portable viable portable water schemes in such areas. The objective of this paper is to determine the potential of basalts in the study area as a source of water for a water bottling and soft drink plant, using geo-electrical and borehole investigations.
2. Geologic Setting
The study area is located at the south western flank of Mount Cameroon (Figure 1). The mountain consists of more than a hundred cones which originate from moderate explosive activity  . It is the only active volcano  of a 1600 km long series of volcanoes known as the Cameroon Volcanic Line (CVL).
The CVL is an alignment of Tertiary continental and oceanic volcanic massifs trending 30˚N from Pagulu Island to Lake Chad  . It is made up of oceanic horsts, namely, Pagalu, Bioko, Sao Tome and Principe, and continental horsts, namely, Mounts Cameroon, Rumpi, Manenguba, Bamboutos, Mbam, and Oku  . The different rock types along the CVL are dominated by granites or syenites, and subordinate intermediate, basic and ultrabasic rocks  . The Mount Cameroon’s eruptions are generally made up of both mafic and silisic lavas. Mafic lavas are associatedwith strombolian-styleexplosive activity along fissures   , which comprise basanites, alkali basalts, hawaiites and mugearites, and occasionally pahoehoe lavas; while the silicic lavas comprise benmoreites, trachytes, phonolites and/or rhyolites  .
Structurally, three major tectonic axes trending northeast-southwest control the volcanic activity of the mountain: the Batoke axis (030˚ - 040˚), the Debundscha axis (060˚ - 070˚) and the Limbe axis (140˚ - 150˚)  and  . The 060˚ - 070˚ corresponds to a major Precambrian shear zone considered as the continuation of the Pernambuco lineament in northeast of Brazil, prior to the continental separation  . Two main faults oriented SSW-NNE are found around Mount: the Tiko fault to the East and the Boa fault to the West.  also identified the Likomba, Esuke, Ombéand Tole-Ekona faults on the East flank. According to  , three factors control the availability of groundwater in the Mount Cameroon region: 1) the degree and depth of weathering of the volcanic rocks, 2) the
Figure 1. Geologic and locationmap of the study area.
absence or presence of permeable surficial materials to allow infiltration into the subsurface, and 3) the presence or absence of faults, joints, fractures and other lineaments.
3. Methods of Study
3.1. Resistivity Method
The geophysical method adopted for this study is the resitivity survey in particular the Schlumberger configuration. The measurements are obtained by introducing a direct current into the ground through a pair of electric current electrodes (A and B) and measuring the resulting potential difference between the electrodes M and N. For a homogeneous ground and an arbitrary electrode arrangement the resistivity can be calculated knowing the current I and the potential difference dV between point M and N by the following equations as elaborated in  .
3.2. Resistivity from VES
At the locations of the four gradient zones as per CST curve, four successive VES curves were recorded. Since the sites with low gradient are expected to be favorable for borehole construction, VES is supposed to reveal the variation of electrical resistivity with depth and in doing so, suggest the possible aquiferous layers and therefore targets for subsequent drilling operations. The apparent resistively values (ρa) were plotted against the current electrodes separation (AB/2) on log-log graph papers.
The vertical electrical soundings were later inverted i.e. quantitatively interpreted using IP2WINautomated inversion software  . This program was used to perform quantitative analysis and interpretation of the field curves. The models resulting from this inversion are illustrated on Figure 3 and Figure 4.
3.2.1. Constant Separation Traversing (CST)
In constant separation traversing, a fixed electrode spacing is chosen and the whole electrode array is moved along a line of traverse or grid after each measurement is made. The value of apparent resistivity is plotted, generally, at the geometric center of the electrode array. Maximum apparent resistivity anomalies are obtained by orienting the profiles perpendicular to the strike of the geologic structure. The results are presented as apparent resistivity profiles or apparent resistivity maps or both  . A 250 meters long profile was covered by 50 measurement points taken at 5 meter interval with a Schlumberger configuration array of AB = 240 m and MN = 10 m the apparent resistivity curve was plotted against the measurement points (Figure 2).
3.2.2. The Vertical Electrical Sounding (VES)
The vertical electrical sounding (VES) consists of a succession of apparent resistivity measurements made with an increasing electrode separation, while keeping the center of the configuration and its orientation fixed  . The variation in
Figure 2. CST curve recorded along the profile presenting fractured zones identified as potential drilling sites.
apparent resistivity will essentially be due to the increasing depth of investigation. Practically, as the spacing between the current electrodes is increased about a center, the total volume of earth included in the measurement also increases both vertically and horizontally   . The field procedure consists in measuring the apparent resistivity as the midpoint of the array is kept fixed while the distance between the current is progressively increased and plotting the apparent resistivity values against half the current electrode spacing on bilogarithmic graph paper. At the locations of the four gradient zones as per CST curve, four successive VES were recorded and plotted as VES 1-4, Figures 4(a)-(d). These data were acquired with ABEM Terrameter 1000 employing the VES and CST techniques of the Schlumberger configuration.
At the positions corresponding to the high gradient zones as per CST curve, three boreholes (F1, F2 and F3) were drilled using rotary drilling techniques to depths indicated by the VES.
3.4. Aquifer Tests
Step-drawdown pumping tests and constant rate pumping tests  were carried out in the boreholes after the holes were equipped and developed. The step-drawdown pumping test was used to evaluate the performance of the well with regard to yield. This test (step-drawdown) was done in three boreholes (F1, F2 and F3) in three steps, with an interval of one hour and increasing pumping rates (Table 1).
The constant rate pumping test was done at pumping rates of 4.4 l/s for 24 hours at F1, 6.9 l/s for 48 hours at F2 and 11.1 l/s for 48 hours at F3. Drawdowns were taken manually using a water level indicator. These values were then used to determine the hydrodynamic parameters of the aquifer  . The goal of this test was to establish the reaction of the aquifer in case of high output and to determine the maximum discharge allowable. The water level was measured using a water level indicator. The results were analysed using the is-Jacob’s method.
3.5. Water Quality
During the constant rate pumping test, water samples were collected and tested in-situ and the results compared with other bottled water brands. The pH, conductivity, iron content and TDS were measured using HANNA Instruments.
4.1. Resistivity from CST
The obtained values of apparent resistivity at each of the sampling points were plotted to generate the apparent resistivity curve along this profile as shown in Figure 3.
K = geometrical factor for a symmetrical array with O middle of AB.
4.2. Resistivity from VES
At the locations of the four gradient zones as per CST curve, four successive VES were recorded. Since the sites with low gradient are expected to be favorable for borehole construction, VES is supposed to reveal the variation of electrical resistivity with depth and in doing so, suggest the possible aquiferous layers and therefore targets for subsequent drilling operations. The apparent resistively
Table 1. Duration and pumping rates of the step-drawdown pumping test carried out in the three boreholes.
values (ρa) were plotted against the current electrodes separation (AB/2) on log-log graph papers (Figure 3 and Figure 4). The vertical electrical soundings were later inverted using IP2WIN automated inversion software  . This program was used to perform quantitative analysis and interpretation of the field curves. The models resulting from these inversions are illustrated on Figure 3 and Figure 4.
Three boreholes were drilled, including F1 (40 m), F2 (80 m) and F3 (110 m). A lithology of the drilled boreholes is shown in Figure 5. A fourth borehole F4 is planned to be drilled on spot of the last major anomalous resistivity gradient of the CST curve. The proposed target is 130 m as per qualitative and quantitative interpretations and the expected lithologies from extrapolation from the previous borehole are summarized on Figure 5.
4.4. Hydraulic Properties
4.4.1. Discharge and Specific Capacity
The results of the step drawdown pumping tests carried out at the three boreholes are found in Table 2. The productivity of F1 characterized by a specific discharge of 3.8 m3/h/m for 14 m3/h. This productivity was observed for one
Figure 3. VES curves recorded at the locations of the first two gradient zones as per CST curve with the corresponding IP2W in models.
Figure 4. VES curves recorded at the locations of the last two gradient zones as per CST curve with the corresponding IP2Win models.
hour pumping duration. F2 has a specific discharge of 1.9 m3/h/m while F3 has the highest specific capacity (5.1 m3/h/m) after one hour of pumping.
The plots of the drawdown versus the pumping rates for the 3 boreholes are shown in Figures 6-8.
The results of the three steps carried out on the three boreholes (F1, F2 and F3) fall on parabolas. The characteristic curve of F1 can be extrapolated to 10 l/s to obtain a drawdown of 10 m (Figure 6); the characteristic curve of F2 can be extrapolated to 10 l/s to obtain a drawdown of about 20 m (Figure 7), and the characteristic curve of F3 can be extrapolated to 20 l/s to obtain a drawdown of 27 m (Figure 8).
4.4.2. Transmissivity and Hydraulic Conductivity from Constant Rate Pumping
During the constant rate pumping test, the specific discharge equals 4.41 m3/h/m after 24 hours pumping at a rate of 16 m3/h. The drawdown stabilized at 4 m. This productivity is also very good. The results of the constant rate pumping test carried in the three boreholes are shown in Table 3.
The recession and recovery curves of all three boreholes are shown in Figures 9-11.
Figure 5. Lithologs of Ombe area.
Table 2. Results of the step-drawdown pumping test.
Table 3. Transmissivity and hydraulic conductivity from constant rate pumping test.
Figure 6. Characteristic curve of F1 (Drawdown vs pumping rate).
Figure 7. Characteristic curve of F2 (Drawdown vs pumping rate).
Figure 8. Characteristic curve of F3 (Drawdown vs pumping rate).
Figure 9. Recession and recovery curves of F1.
The maximum exploitation rates of the different boreholes obtained from a software simulation extrapolation after pumping the boreholes for a month, a year and ten years are found in Table 4. The lowest tolerable drawdown at F1 is 10 m. This drawdown will not be attained after 10 years of exploitation at an average pumping rate of 4.0 l/s. F2 can be exploited at an average pumping rate of
Figure 10. Recession and recovery curves of F2.
5.9 l/s, without the lowest tolerable drawdown at 22 m being attained after 10 years of exploitation. Similarly, the lowest tolerable drawdown of 43 m at F3 will not be attained after 10 years of exploitation at a pumping rate of 18 l/s.
4.4.3. Water Quality
Chemical analysis of the water samples from the boreholes showed that bicarbonate is the most dominant anion (231.8 - 268.4 mg/l), whereas calcium is the
Figure 11. Recession and recovery curves of F3.
most dominant cation (40 - 72 mg/l) (Table 5). concentrations range between 0 (F2) and 1.98 mg/l (F1). pH ranges between 7.6 (F2) and 7.83 (F3), which is indicative of neutral waters; and EC values ranges 323 (F2) - 361 µS/cm (F1).
Table 4. Maximum exploitation rates of the different boreholes.
Table 5. Selected chemical and physico-chemical parameters of the boreholes and WHO permissible limits.
A comparison was made between these samples and available water brands on the local market (Figure 12). The samples from the boreholes had better physico-chemical characteristics that those of local water brands (Alfa, Beta, Gamma and Delta). Water samples from these boreholes contain more than 60 mg/lof Ca2+, while the available water brands contain less than 50 mg/l; the borehole water sample also contain more Mg2+ than all water brands. pH values are similar in all water samples; whereas EC varies considerably.
5. Interpretation and Discussion
5.1. Structure and Aquiferous Zones
The CST profile as shown in Figure 3 is characterized by a series of highs and lows with the highest peak value of 70 Ohm・m and lowest value of 27 Ohm・m. The curve presents four major horizontal gradient zones with significant drop in the values of apparent resistivity. These gradient zones suggest an underlying fractured and/or weathered basement and these fractures are potential aquifers  . These results are good and in agreement with those obtained by  who used the same technique to decipher the existing subsurface stratification in the River Segen Dam Site in northern Ethiopia. Similar conclusions were reach by  . On the basis of the results obtained, points F1 (5, 37), F2 (21, 31), F3 (30, 26) and F4 (35, 27) were selected as more probable site for further investigation by VES and potential sites for drilling.
5.2. Resistivity of the Formations
The four VES curves recorded at F1, F2, F3 and F4 are very similar in shape as well as in amplitude and all show a four-layered earth profile type KHKH, corresponding to a typical conductive/resistive/conductive/resistive geoelectrical
Figure 12. Comparison of the physico-chemical characteristics of the water samples in the boreholes and those of local brands.
section and implying that in addition to the presence of a top unconfined aquifer, there are two confined fractured zone within the fresh basement. The results of the four VES curves are summarized in Table 6.
The general shape of the different curves is indicative of a geological section made of resistive top soil, which becomes more conductive at around 5 - 6 m depth, with an apparent resistivity varying from 135 to 77 Ohm・m approximately. This depth corresponds to the top of the clayey layer with basaltic boulders. This layer overlies a layer of relatively lower resistivity that extends to the depth of about 40 m and corresponds to a porous aquifer made up of fractured/altered basalts.
Beneath this layer is found a resistive bed of unfractured/unaltered basalt, extending from 40 to 55/60 m. Another low resistivity layer, possible aquifer is detected around 80 m depth and extends to 120 m and corresponds to much fractured basalts. Below the depth of 120 m, the bedrock is made of fresh unfractured basalt. It is worth noticing that all the VES curves reveal low variations of apparent resistivities (80 to 175 Ohm・m) up to 80 m depth. This corresponds to the interval of deposition of altered volcanic materials.
5.3. Lithofacies and Aquiferous Formations
The geological sections obtained from the processing of drilling data have been compared with the results of the vertical electrical sounding. From these comparisons, it appears that the interpretation of sounding curves corroborate the results obtained from the lithological sections made up of a top soil, sandy clay, fresh basalt, fractured basalt, fresh basalt and a last layer of fractured basalt. This means that the number and the order of occurrence of the aquifers match the drilling results. The only discrepancy is in the depths prediction from the VES data that do not match with the actual drilling depths. These discrepancies could
Table 6. Results of the interpretation of the VES curves.
be explained by the suggestion that the weathered layer thickness includes partly weathered the semi weathered layer as well or the problem of equivalence. The identified aquifers from the different boreholes are shown in Table 7.
5.4. Hydraulic Properties
The different transmissivities obtained are less than the Deccan volcanic aquifer in India  . The same  obtained mean transmissivities of 3 × 10−3 m2/s in the sheet joints of vesicular amygdaloidal basalt unit; 5 × 10−4 m2/s in the vertical joints in compact basalt unit, and 2.2 × 10−3 m2/s from boreholes crossing both sheet joints in vesicular amygdaloidal basalt and vertical joints in compact basalt. However, the hydraulic conductivities obtained fall in the ranges obtained by  in the Snake River Plain Aquifer, in Idaho (USA). Results obtained from  showed that most of the hydraulic conductivities of wells found in thin, tube-fedpahoehoe flows and near-vent volcanic deposits (shelly pahoehoe and slab pahoehoe flows and bedded scoria, spatter, and ash) range between 3.5 × 10−6 m/s and 3.5 × 10−4 m/s.
5.5. Groundwater Quality
To an important extent, the chemical composition of mineral waters is determined by the composition of the rock it is abstracted from. Interactions between water and rock influence the mineral character of groundwater  by affecting its quality. The abundance of bicarbonate ions might be due to the fact that calcium carbonate occurs widely as secondary minerals in igneous rocks  . It is also possible that the high concentrations are related to CO2 degassing of Mount Cameroon  . There might have been a reaction between this CO2 and groundwater; leading to the formation of HCO3 or between biogenic and atmospheric CO2 and groundwater    .
Potassium concentrations in volcanic rocks are usually low due to the resistance to weathering of potassium-bearing minerals and its fixation in the formation of clay minerals  . The study area also has distinct depletion in potassium  . Moreover, the lavas in this area form a basanite series with minerals such as clinopyroxenes (magnesium and/or iron rich), olivines (forsterite) and plagioclase (calcium rich). Therefore the principal cations in the area are calcium
Table 7. Depths of different aquifers in the area.
and magnesium ions. These cations are also found dominant in groundwater samples. Analysis of water samples from the three boreholes showed that their physico-chemical characteristics are within the  norms.
Geoelectrical survey was conducted on the southwestern flank of Mt. Cameroon, area underlain by basalts. The occurrence of groundwater in these rocks was order to determine their potentials as a source of water for a water bottling and soft drink plant. To achieve this, the constant separation traversing (CST) of the resistivity method was used and the data qualitatively analyzed to determine the occurrence of fractures as most probable borehole sites. As a result of this procedure, four sites were selected and at each of the preselected sites, a vertical electrical sounding (VES) was conducted and interpreted. The results of the interpretation of the VES data in the area showed a four layered earth profile model type KHKH, corresponding to a typical weathered/fractured confined aquifer type curve. This profile was subsequently confirmed from borehole information after drilling at three of the four selected sites. A fourth borehole still to be drilled is expected to confirm the already obtained results. The comparison of results of VES curve interpretations with the geological sections obtained from boreholes curves is in good agreement. It shows that the number and the order of occurrence of the aquifers estimated from VES data match the drilling results. The only discrepancy is in the depths prediction from the VES data that do not match with the actual drilling depths, discrepancies that could be explained by the suggestion that the weathered layer from VES could include partly weathered the semi weathered layer as well. A Step-drawdown pumping test was done in three boreholes (F1, F2 and F3) in three steps, with an interval of an hour and increasing pumping rates. A constant rate pumping test was done at pumping rates of 4.4 l/s for 24 hours at F1, 6.9 l/s for 48 hours at F2 and 11.1 l/s for 48 hours at F3. The results of these tests were analysed using Theis-Jacob’s method, and the transmissivity, hydraulic conductivity, the storage coefficient and the radius of influence were computed. These results reveal the existence of a high permeability zone around the three boreholes, with very good transmissivity that can safely produce water at a rate of 4.0 l/s, 5.9 l/s and 18 l/s for F1, F2 and F3, respectively, for a minimum of 10 years. Chemical analysis done in-situ during the pumping test and on samples at the laboratory reveal that this water has a lower iron content and better mineralisation compared with other bottled water brands sold on the Cameroon market.
Vertical Electrical Sounding
 Teikeu, W.A., Njandjock, P.N., Ndougsa-Mbarga, T. and Tabod, T.C. (2012) Geoelectric Investigation for Groundwater Exploration in Yaounde Area, Cameroon. International Journal of Geosciences, 3, 640-649.
 Ateba, B., Dorbath, C., Dorbath, L., Ntepe, N., Frogneux, M., Aka, F.T., Hell, J.V., Delmond, J.C. and Manguelle, D. (2009) Eruptive and Earthquake Activities Related to the 2000 Eruption of Mount Cameroon Volcano (West Africa). Journal of Volcanology and Geothermal Research, 179, 206-216.
 Ako, A.A., Jun, S., Takahiro, H., Makoto, K., Akoachere, R.A., George, E.N., Gloria, E.T.E. and Alain, L.F.T. (2012) Spring Water Quality and Usability in the Mount Cameroonarea Revealed by Hydrogeochemistry. Environmental Geochemistry and Health, 34, 615-639.
 Yeh, H.F., Lin, H.I., Wu, C.S, Hsu, K.C., Lee, J.W. and Lee, C.H. (2015) Electrical Resistivity Tomography Applied to Groundwater Aquifer at Downstream of Chih-Ben Creek Basin, Taiwan. Environmental Earth Sciences, 73, 4681-4687.
 Fields, J. and Halihan, T. (2015) Preliminary Electrical Resistivity Surveys of Mount Judea Alluvial Sites. 2nd Quarter 2015 Report. Boone Pickens School of Geology, Oklahoma State University, 24 p.
 Akinrinade, O.J. and Adesina, R.B. (2016) Hydrogeophysical Investigation of Groundwater Potential and Aquifer Vulnerability Prediction in Basement Complex Terrain—A Case Study from Akure, Southwestern Nigeria. De Gruyter Open.
 Wandji, P., Tsafack, J.P.F., Bardintzeff, J.M., Nkouathio, D.G., Dongmo, A.K., Bellon, H. and Guillou, H. (2009) Xenoliths of Dunites, Wehrlites and Clinopyroxenites in the Basanites from Batoke Volcanic Cone (Mount Cameroon, Central Africa) Petrogenetic Implications. Mineralogy and Petrology.
 Njome, M.S., Suh, C.E. and de Wit, M.J. (2009) The Mount Cameroon Volcano, West Africa: An Active Link between Recent Eruptive and Mantle Signatures of the Deep Past beneath the Margins of Africa. 11th SAGA Biennial Technical Meeting and Exhibition, 533-539.
 Kamgang, P., Chazot, G., Njonfang, E., Ngongang, N.B.T. and Tchoua, F.M. (2013) Mantle Sources and Magma Evolution beneath the Cameroon Volcanic Line: Geochemistry of Mafic Rocks from Bamenda Mountains (NW Cameroon). Gondwana Research, 24, 727-741.
 Ngako, V., Njonfang, E., Aka, F.T., Affaton, P. and Nnange, J.M. (2006) The North-South Paleozoic to Quaternary Trend if Alkaline Magmatism from Niger-Nigeria to Cameroon: Complex Interaction between Hotspots and Precambrian Faults. Journal of African Earth Sciences, 45, 241-256.
 Suh, C.E., Sparks, R.S.J., Fitton, J.G., Ayonghe, S.N., Annen, C., Nana, R. and Luckman, A. (2003) The 1999 and 2000 Eruptions of Mount Cameroon: Eruption Behaviour and Petrochemistry of Lava. Bulletin of Volcanology, 65, 267-281.
 Ayonghe, S.N., Ntasin, E.B., Samalang, P. and Suh, C.E. (2004) The June 27, 2001 Landslide on Volcanic Cones in Limbe, Mount Cameroon, West Africa. Journal of African Earth Sciences, 39, 435-439.
 Ngako, V., Affaton, P., Nnange, J.M. and Njako, T. (2003) Pan African Tectonic Evolution in Central and Southern Cameroon: Transpression and Transtension during Sinistral Shear Movements. Journal of AfricanEarth Science, 36, 207-214.
 Akoachere, R.A., Mbotake, I.T., Yaya, O.O. and Oteze, G.E. (2007) Using Morphometric Analysis to Evaluate Groundwater Hydraulic Characteristics in the Piedmont Regions of Mount Cameroon. 19th Annual Conference of Nigerian International Association of Hydrogeologists.
 Raimi, J., Abdulkarim, M.S., Hamidu, I. and Arabi, A.S. (2011) Application of Schlumberger Vertical Electrical Sounding for Determination of Suitable Sites for Construction of Boreholes for Irrigation Scheme within a Basement Complex. International Journal of Multidisciplinary Science and Engineering, 2, 81-84.
 Bairu, A. (2013) Application of Vertical Electrical Sounding and Horizontal Profiling Methods to Decipher the Existing Subsurface Stratification in River Segen Dam Site, Tigray, Northern Ethiopia. Journal of Environment and Earth Science, 3, 42-54.
 Abraham, J.D. and Lucius, J.E. (2004) Direct Current Resistivity Profiling to Study Distribution of Water in the Unsaturated Zonenear the Amargosa Desert Research Site, Nevada. U.S. Geological Survey, Reston, Virginia 2004 Revised and Reprinted.
 Kulkrani, H., Deolankar, S.B., Lalwani, A., Joseph, B. and Pawar, S. (2000) Hydrogeological Framework of the Deccan Basalt Groundwater Systems, West-Central India. Hydrogeology Journal, 8, 368-378.
 Anderson, S.R., Kuntz, M.A. and Davis, L.C. (1999) Geologic Controls of Hydraulic Conductivity in the Snake River Plain Aquifer at and near the Idaho National Engineering and Environmental Laboratory, Idaho. Water-Resources Investigations Report 99-4033, US Geological Survey, Idaho, 44 p.
 Edmunds, W.M., Shand, P., Hart, P. and Ward, R.S. (2003) The Natural (Baseline) Quality of Groundwater: A UK Pilot Study. The Science of the Total Environment, 310, 25-35.
 Suh, C.E., Luhr, J.F. and Njome, M.S. (2008) Olivine-Hosted Glass Inclusions from Scoriae Erupted in 1954-2000 at Mount Cameroon Volcano, West Africa. Journal of Volcanology and Geothermal Research, 169, 1-33.
 Jeong, C.H. (2001) Effects of Land Use and Urbanisation on Hydrochemistry and Contamination of Groundwater from Taejon Area, Korea. Journal of Hydrology, 253, 194-210.
 Olobaniyi, S.B. and Owoyemi, F.B. (2006) Characterization by Factor Analysis of the Chemical Facies of Groundwater in the Deltaic Plain Sands Aquifer of Warri, Western Niger Delta, Nigeria. African Journal of Science and Technology-Science and Engineering Series, 7, 73-81.
 Rajmohan, N. and Elango, L. (2006) Hydrogeochemistry and Its Relation to Groundwater Level Fluctuation in the Palar and Cheyyar River Basins, Southern India. Hydrological Processes, 20, 2415-2427.