Determination of Groundwater Potential Zones on the Eastern Slope of Mount Cameroon using Geospatial Techniques and Seismoelectric Method

doi:10.21203/rs.3.rs-2990795/v1

Download PDF

Research Article

Determination of Groundwater Potential Zones on the Eastern Slope of Mount Cameroon using Geospatial Techniques and Seismoelectric Method

https://doi.org/10.21203/rs.3.rs-2990795/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

In complex hard-rock terrain where the population suffers from water scarcity, the use of site-specific single-method surface geophysical survey to prospect for groundwater is a common practice. However, this has not completely solved the problem of drilling dry wells. To combat this failure, geospatial data and the analytical hierarchy process (AHP) are being used as supporting techniques to increase the chance of success. This study had the aim of generating groundwater potential zones (GWPZ) map of the complex strato-volcanic terrain of Buea by integrating geospatial techniques, AHP and the seismoelectric method. Six factors influencing groundwater recharge were used to define the groundwater potential zones. Via the AHP and weighted overlay methods, five classes of groundwater potential zones were delineated: very poor, poor, moderate, good, and very good. More than half of the study area has good to very good groundwater potential. Lastly, the GWPZ model is validated with an iso-conductivity map from a seismoelectric survey and existing static water level data. The iso-conductivity map revealed four main conductivity zones that correlate positively with the GWPZ map. The northwestern part of the study area is characterised by moderate groundwater potential, poor formation conductivity and dry boreholes. There is an uneven distribution of groundwater and variable water table depth in the area. The results of this study are very encouraging, and the integrated approach used has proven to be efficient in determining groundwater potential zones in complex volcanic terrain.

Analytical hierarchy process

Geospatial techniques

Seismoelectric method

Groundwater potential

volcanic terrain

More than 65% of the study area has good to very good groundwater potential and these zones include the older basanite and lowland regions.
The conductivities of the younger and recent basanites are very low coinciding with moderate to poor groundwater potential and occurrence.
The groundwater table of Buea varies immensely, and it is directly proportional to the elevation.

Inadequate management practices, lack of national water policies, poor maintenance and deterioration of infrastructures, corruption, and the lack of adequate funding are some of the problems in developing countries that have resulted in water scarcity [1]–[4]. As a result, the population is now dependent on groundwater as an alternative and most reliable source of potable water. In the absence of springs, the only other access to groundwater is through wells and boreholes. Nevertheless, the global distribution of groundwater is very uneven [5], [6], especially in active volcanic terrains where there is constant deformation and change in landforms. Additionally, the global groundwater table depth varies extensively [7] even within small regions like Buea. These are some of the causes of unproductive and deeper water level boreholes in Buea.

Buea, the study area, is situated on the active Mount Cameroon volcano and is noted for water scarcity especially in the upland villages of Bova, Bonakanda, Bokova, Woteva, Ewonda, Bokwango and Bwassa with ever expanding populations. To solve this problem, individuals have been extracting groundwater through boreholes. However, even with the practice of site-specific geophysical surveys, there have been drilling failures with many existing dry boreholes. These frequent drilling failures are due to the geological complexity of the volcanic terrain, the lack of detailed geological and hydrogeological studies, and insufficient geophysical data. Geophysical surveys are more effective when used with a good understanding of the geology and hydrogeological conditions of the area to be exploited. Hence, the assessment and zoning of the groundwater potential and the planning and management of groundwater resource is fundamental.

In addition to geophysical surveys, geospatial techniques can be used for groundwater exploration [8]–[10]. Mapping of groundwater potential zones (GWPZ) is the prerequisite to detailed geophysical surveys for borehole placement. To evaluate the groundwater potential of a region and enhance drilling success, researchers around the world use remote sensing and GIS techniques which are time and cost effective [8], [9], [11]–[14]. Groundwater potential can be inferred from many groundwater recharge influencing factors such as lineament density, land cover, land use, drainage density, vegetation, aspect, rainfall, soil type, e.t.c., depending on regional peculiarities [12], [15]–[18]. These factors are classified and weighted using multi-influencing factor (MIF) and Analytic Hierarchy Process (AHP) methods [8], [12], [19].

To help solve the water scarcity problem in Buea, this study is the first to integrate remote sensing, geospatial techniques and geophysical survey to demarcate the different groundwater potential zones. This was achieved by using factors influencing groundwater recharge potential in the area, which are lithology, slope, drainage density, lineament density, elevation, and land use. Each thematic map, including lithology, elevation, slope, drainage density, lineament density and land use provides a firm indication of the occurrence of groundwater and its potential. Further validation of the groundwater potential zones was done using existing springs and boreholes with the static groundwater level data. The results obtained from this study will serve as a guide and tool for siting boreholes, to aid decision-makers and local authorities for future planning projects such as the zoning of new settlements and groundwater resource management.

Buea is situated on the eastern slope of Mount Cameroon (Fig. 1), which is the highest mountain in West Africa standing at 4095 m.a.s.l. [20]. It is the capital of Fako Division in the Southwest Region of Cameroon with a coastal equatorial climate [21], having two major seasons: the dry season and the rainy season. The vegetation type in this region is the tropical evergreen forest that ranges from primary thick forest, secondary forest, shrubs to savannah towards the peak of the mountain [22].

The geology of Buea is indicative of that of Mount Cameroon, an active volcano that last erupted in the year 2000 [23], [24]. Mount Cameroon is one of the volcanic centres of the Cameroon Volcanic Line (CVL) that runs from Pagalu Island in the Atlantic Ocean to Lake Chad [25]–[29]. Mount Cameroon is a strato-volcano that erupts both effusions of mafic lava and strombolian-style explosive lava [23]. The lithology of the Buea is composed of basanites, trachyte and trachybasalt [20], [23], [24], [30] (Fig. 1). These rocks show extensive fracturing. The fractures have a more or less N-S orientation.

To achieve the main aim of the current study, we followed the procedure indicated in Fig. 2. Remotely sensed data were used to generate thematic maps, which were integrated to produce the GWPZ map. An area of 12 km x 11 km of Buea was selected for the case study. Field verification and validation of the groundwater potential zones were conducted using the seismoelectric surface geophysical method and existing springs and static groundwater level data. The different stages are discussed in the following sections.

3.1. Geospatial techniques

Geospatial techniques were used to delineate the different groundwater potential zones. The steps used in achieving this were: (1) identifying and mapping the different surface factors that are believed to control the recharge, runoff, occurrence and storage of groundwater; (2) reclassification of each factor with respect to its influence on recharge; (3) assigning each thematic layer a rank and weight based on its relative importance to groundwater recharge; and finally, (4) integrating all the thematic maps in ArcGIS Pro using the weighted overlay method to generate the GWPZ map [10].

Six factors that include lithology, elevation, slope, lineament density, drainage density and land use were used in this study. The elevation, slope and drainage density thematic maps were produced from SRTM DEM data with 30 m spatial resolution (SRTM1N04E009V3) using ArcGIS Pro. Lineament density and land use thematic maps were generated from band 7 of Landsat 8 OLI data (LC08_L1TP_187057_20170507_20170515_01_T1) that covers the study area [31], using ArcGIS Pro. The slope, elevation, drainage and lineaments densities were generated using the Spatial Analyst Tools while the land use map was created using the Image Classification Tools of ArcGIS Pro. The lithology map was reproduced from the geological map of Mount Cameroon after Ateba et al. (2009).

Lineament extraction to produce the lineament density map is a rigorous process that involves pre-processing, image enhancement and image transformation of the Landsat 8 multispectral images using ENVI and PCI Geomatica software [32], [33]. Pre-processing involves radiometric correction (i.e., Radiometric Calibration and FLAASH Atmospheric Correction) and data resizing to the area of interest. Image enhancement involves Spectral Enhancement and Spatial Filtering using Directional Convolution Filter. Image transformation involves Image Sharpening using Gram_Schmidt Pan sharpening, Principal Component Analysis, Minimum Noise Fraction (MNF), and Independent Component Analysis [33], [34]. The lineaments were then extracted from the processed band 7 images using an automatic method[33], [35] with PCI Geomatica Banff. The Landsat image is processed to increase the qualities of tones and hues, image textures, lineaments and their trends [32]. Lineaments here represent fractures, joints, and faults that act as pathways for water infiltration and groundwater flow.

3.2. Analytic Hierarchy Process (AHP)

To avoid errors and to preserve the geometric properties of the data imagery, all the thematic layers were harmonised and projected onto a common Universal Transverse Mercator (UTM) coordinate system (WGS84 Zone 32 N) at 30 m resolution. The six thematic maps produced from the above-mentioned processes were used as the six factors to generate the GWPZ map using the AHP. The AHP is a multi-criteria decision making technique originally developed by Prof. Thomas Saaty [36]. It relies on weights calculation by finding the dominant right eigenvector of a positive reciprocal decision matrix based on pairwise comparison inputs and using the fundamental scale of values to represent the intensities of importance [36], [37]. Using the Analytic Hierarchy Process Online Software (AHP-OS) developed by Klaus Goepel the six criteria (factors) were used for the pairwise comparison (pwc), based on the equations below [37].

$$pwc=\left({a}_{1},{a}_{2}\dots {a}_{npc}\right),\left({x}_{1},{x}_{2}, \dots {x}_{npc}\right)$$

With integers 𝑎_𝑖 ∈ [0, 1], 𝑥_𝑖 ∈ [1, 𝑀], M = 9 and 𝑖=1 …𝑛𝑝𝑐, where 𝑛𝑝𝑐 is the number of pairwise comparisons and M is the decision matrix. For n criteria, the n x n decision matrix is then filled from pwc.

$$npc=\frac{{n}^{2}-n}{2}$$

The six influencing factors with their weighted percentages were geoprocessed using Weighted Overlay under the Spatial Analyst Tools of ArcGIS Pro to generate the GWPZ map. All the input criteria were reclassified to integer rasters before they were used in Weighted Overlay. An output scale value of 1–5 selected from the Scale setting was used to assign scale values for the criteria after setting the Evaluation Scale at 1–5, where 1 represents the lowest and 5 is the most favourable for groundwater recharge and occurrence.

3.3. Verification and Validation

Validation and field verification of the GWPZ map was carried out via surface geophysical surveys using the seismoelectric method and investigating the existing springs and static water level.

The seismoelectric method involves the propagation of seismic wave in a medium that induces an electrical field or causes radiation of an electromagnetic wave [38]. Hence, the seismoelectric method involves the use of a seismic source and electric or magnetic field receivers. The ATS GeoSuite App was used as the receiver, while a sledgehammer generated the seismic wave. Sounding was done on 253 points spread over the study area. The data were collected during the dry season to avoid seasonal aquifers and raised water table due to infiltration from rainwater. It was entered into an Excel sheet, which was added to ArcGIS Pro and then inserted into the geostatistical kriging method for interpolation to produce an iso-conductivity map of the study area that was used as field verification and validation of the GWPZ.

Other means of validation used are existing springs and the static water level in boreholes. Data from 241 producing boreholes, 9 dry boreholes and 19 springs were used to validate the GWPZ.

The results obtained on various aspects namely lithology, elevation, land use, lineament density, drainage density, slope, seismoelectric survey, existing springs and static water level are discussed below.

4.1. Factors Influencing Groundwater Recharge

4.1.1. Lithology

Lithology is the most important factor controlling groundwater recharge and occurrence in the study area (Table 1). For groundwater to be stored in aquifers, there must exist suitable lithological conduits for infiltration and groundwater flow. These conduits are provided by the presence of fractures and the high primary porosity of pyroclastic rocks that increase with age. The lithology of Buea is made of basanites, trachytes and trachybasalts of different ages (Fig. 3a). These rocks occur as alternate layers of lava flow and pyroclastic rocks. The oldest basanites and trachytes/trachybasalts occur in the central and southeastern part of the study area and divide the younger rocks of the north and southwest. The recent rocks can be found in the northwest, ascending Mount Cameroon. Based mainly on permeability, the rocks were reclassified and then ranked according to their influence on groundwater recharge (Table 1).

Table 1

Class ranking and percentage of thematic weight.
Factor/Criterion (units)	Class/Value range	Rating	Reclassified value	Scale/Rank	Weight/Priority
Lithology	Younger basalts1 Recent basalts Younger basalts2 Older basalts	1	1 2 3 4	3 1 3 5	37.2
Elevation (m)	420.001–720 720.001–1063 1063.01–1500 1500.01–2112 2112.01–3007	2	1 2 3 4 5	5 4 3 2 1	24.5
Slope (°)	0.00000728–3.73 3.74–5.11 5.12–8.84 8.85° − 18.9 19° − 45.8	3	1 2 3 4 5	5 5 4 2 1	16.9
Lineament density (km/km²)	Very low Low Moderate High Very high	4	1 2 3 4 5	3 3 4 5 5	12.0
Drainage density (km/km²)	Very low Low Moderate High Very high	5	1 2 3 4 5	5 5 4 3 2	5.1
Land use	Build up / Developed Forest Herbaceous Planted / Cultivated	6	1 2 3 4	3 1 3 4	4.2

4.1.2. Elevation

Elevation is the second most important factor that greatly influences groundwater recharge in the study area and has a weight value of 24.5 (Table 1). Lower elevations offer increased residence time for surface water, thus greater infiltration and lesser runoff. Lesser amount of infiltration and larger amount of runoff are associated with areas of higher elevation [12]. The elevation of the study area ranges from 420 m to 3007 m from east to west (Fig. 3b) above sea level. The elevation is divided into five subclasses, with the 420 m to 720 m subclass occupying the greatest area. The five subclasses were further reclassified and assigned ranks. The lowest elevation class was given a scale of 5.

4.1.3. Slope

Slope is one of the main factors that influence the recharge of aquifers in terms of surface water infiltration and runoff. It is measured in degrees. The degree of slope determines whether a zone favours infiltration due to high residence time or favours runoff [11]. The steeper the slope, the higher the runoff and the lesser the infiltration of surface water. The steepest slopes occur in the western part at the highest elevation. These slopes of greater than 19° have very low infiltration and cover about 1/5 of the total surface area of the study area (Fig. 4a). More than half of the study area is characterised by gentle slopes of less than 9°. The slope map is subdivided into five slope classes; 0.00–3.73°, 3.74–5.11°, 5.12–8.84°, 8.85° − 18.9°, 19° − 45.8° and further reclassified into very low, low, moderate, high, and very high (Table 1).

4.1.4. Land use

Land use is the least factor influencing recharge in this study (Table 1). Nonetheless, it is a significant factor affecting groundwater recharge process [39]. The land use includes build up or developed areas, forest, herbaceous and cultivated areas (Fig. 4b), classified using supervised classification on the composite image of band 423. The reclassified and ranked land use is shown in Table 1. More than 50% of the study area is dominated by cultivated and herbaceous land. The cultivated and herbaceous land have the most effect on surface water infiltration into aquifers by increasing infiltration rate. The forest zone has a thick canopy and occurs on steep slopes, thus reducing surface water percolation. Build up and developed areas create a compacted terrain that seals the ground surface, thus preventing or reducing groundwater recharge.

4.1.5. Lineament density

Lineament density is the total length of linear features per unit area [40]. The greater the amount of lineament concentration and extension in an area, the greater the lineament density, and the greater the lineament density, the higher the probability of groundwater recharge. Lineaments here may represent surface expressions of faults, folds, fractures, lithological boundaries, or rock contacts, through which surface water could percolates or groundwater could flows. The study area has a high lineament density (Fig. 5a). The highest lineament density occurs in Mile 16, Mile 17, Molyko, Small Soppo, Buea Town, Lissoka, and between Bokova and Bonakanda. The northwest of Ewonda, Bova and Bonakanda show the lack of lineament. High lineament density zones are areas of high infiltration and recharge. The high lineament density zones were ranked the highest on a scale of 5 while the low lineament density zones were given ranks of 3 (Table 1).

4.1.6. Drainage density

Drainage density is the ratio of the overall length of channels per unit of area [41]. The drainage pattern in the study area is dendritic, with most of the drainage originating from the mountain and highlands in the west (Fig. 5b). The drainage density is measured in km/km² and varies from 11.54 to 96.81 km/km². Drainage reduces infiltration due to runoff. The low drainage density areas are positive for high infiltration rates, whereas the high drainage density areas are characterised by low infiltration rates. Low drainage density zones are found in the eastern lowland of the study area (Fig. 5b and Fig. 3b). The drainage density map was reclassified into five categories of very low, low, moderate, high and very high (Table 1).

4.2. Groundwater potential zonation

The groundwater potential zones of the study area were delineated by integrating lithology, elevation, drainage density, slope, lineament density and land use thematic maps using geospatial techniques. The resulting weights for the influencing factors based on 15 pairwise comparisons by the AHP-OS are presented in Table 2. Figure 6 shows the bar chart of the resulting weights and their percentage error. The resulting weights are based on the principal eigenvector of the decision matrix below (Table 3). Through the weighted overlay technique five groundwater potential zones were delineated, namely very good, good, moderate, poor, and very poor (Fig. 7). From the groundwater potential zones map, the best (very good) groundwater potential zone is concentrated in the southeastern part of the study area, whereas the northwestern region has the poorest groundwater potential. Large parts of the study area fall within the good and very good groundwater potential zones.

Table 2

Resulting weights for the six factors based on 15 pairwise comparisons. Consistency Ratio CR = 4.5%
Factor	Weight	Rank
Lithology	32.9%	1
Elevation	22.7%	2
Slope	20.5%	3
Lineament density	12.2%	4
Drainage density	7.1%	5
Land use	4.5%	6

Table 3

Decision Matrix. Principal eigen value = 6.284, Eigenvector solution: 6 iterations, delta = 8.2E-8
	Lithology	Elevation	Slope	Lineament density	Drainage density	Land use
Lithology	1	2.00	2.00	3.00	4.00	5.00
Elevation	0.50	1	1.00	3.00	4.00	4.00
Slope	0.50	1.00	1	3.00	2.00	4.00
Lineament density	0.33	0.33	0.33	1	3.00	4.00
Drainage density	0.25	0.25	0.50	0.33	1	2.00
Land use	0.20	0.25	0.25	0.25	0.50	1

4.3 Model validation

4.3.1. Seismoelectric survey

The results of the 1D seismoelectric grid survey of the study area are presented in Fig. 8. There is great variability in the conductivity of the formations between 40 to 60 m below the surface. Bonakanda, Bova and the northwestern zone of the map show very low conductivity. Most of the study area is composed of moderate to high formation conductivity. High conductivity implies high electrical conductance of the medium due to the presence of water. The highest conductivities were registered in Mile 16, Mile 17, Muea and Molyko in the southeast.

4.5.2. Existing springs and static water level

The existing boreholes in the study area were mapped and their static water levels were used in validating the groundwater potential zones map. The study area has a lot of producing boreholes, the majority of which occur in the very good groundwater potential zones (Fig. 9). The static water level in the boreholes varies considerably, from 10 m in the very good groundwater potential zone to 100 m in the good to moderate zones. Springs only occur in the very good groundwater zone, particularly where the groundwater table is low (10–40 m below the surface). Dry boreholes occur generally in the moderate groundwater potential zone at Bova, Bonakanda and Bokova. These holes range from 60–110 m deep.

5.1. Groundwater potential zones

The six thematic maps that were integrated to generate the GWPZ map have different degrees of influence on aquifer recharge and thus groundwater occurrence and distribution potential as demonstrated in similar studies [10], [11], [15], [39]. The weighted thematic maps have a consistency ratio of 4.5% (0.045) (see Table 2) implying the reasoning of the AHP is consistent [15], [18] and the GWPZ map is reliable. The influencing attributes under lithology, elevation, slope, lineament density, drainage density and land use, have varied groundwater potentials classified as very poor to very good (Table 4). The importance of each thematic layer to groundwater recharge and occurrence are further discussed.

Low elevation, low slope angles, high lineament densities, low drainage densities, herbaceous and cultivated land, and the older permeable basanites epitomise good to very good groundwater potential zones in the study area (Table 4). On the other hand, poor groundwater potential zones are typified by high elevation, high slope, low lineament density (little or no surface lineaments), high drainage density, forest and recent basanites. Large parts of the study area fall within the good and very good groundwater potential zones. This is mainly due to the geological nature of the terrain that consists of old weathered and fractured basanites of high permeability, a similar situation noted by Serele et al. (2020) in Madagascar. Geology (lithology and fractures) and topography also play major roles in infiltration and runoff [42]. Hence, the influence of land use on groundwater recharge is dependent on the underlying geologic conditions and the topography [43] other than manmade structures. This explains why we have high groundwater potential in build-up and developed areas, even though these areas are characterised by low infiltration and high runoff.

The ability of formations to host groundwater depends on some geologic factors like rock age, fracturing and unconformity, all of which influence infiltration, percolation and groundwater flow. Age increases the porosity of hard-rocks, as they become more porous and permeably due to weathering and deformation. These elucidate why the older basanites zone has the greatest groundwater potential (Fig. 7 and Table 4). Fracturing results in secondary porosity and permeability. Thus, high lineament density in the study area might be revealing of secondary porosity and permeability represented by fracture systems serving as conduits for the movement and storage of groundwater [12].

Elevation and slope have an inverse relationship with infiltration but a direct proportionality with runoff because surface runoff is low on low slopes and plains allowing more time for surface water infiltration into the subsurface and vice versa [8], [16], [43]. Low elevations and gentle slopes in the eastern part of the study area have good to very good groundwater potential (Fig. 7 and Table 4). Similarly, low drainage density, which is associated with low elevations and gentle slopes favours reduced runoff, consequently increased infiltration and good groundwater recharge [16], [39].

Table 4

Classification of groundwater potentiality of thematic layers.
Factors	Attributes	Groundwater potential
Lithology	Younger basalts Recent basalts Older basalts	Moderate to good Poor to very poor Good to very good
Elevation	Very low Low Moderate High Very high	Very good Good Moderate Poor Very poor
Slope	Gentle sloping Moderate sloping Moderate to steep sloping Steep sloping Very steeply sloping	Very good Good Moderate Moderate to poor Poor
Lineament density	Very low Low Moderate High Very high	Poor Moderate Good Good to very good Very good
Drainage density	Very low Low Moderate High Very high	Very good Good Moderate Poor Very poor
Land use	Build up / Developed Forest Herbaceous Planted / Cultivated	Moderate Poor Good Very good

5.2. Validation of the groundwater potential zones map

Validation of groundwater potential zones map using surface geophysics and borehole yield data has been carried out successfully elsewhere [11], [13]. The correlation of the GWPZ map with the seismoelectric (iso-conductivity) map gives agreeable results. The iso-conductivity map (Fig. 8) highlights the very good groundwater potentiality of the southeastern part of the study area and the poor groundwater potentiality of the northwestern part (Fig. 7). It shows that the poor groundwater potential zone has very high apparent resistivity (low conductivity) and further confirms the reason of the dry holes in Bova, Bonakanda and Bokova is due to very low conductivity of the formations. The very low conductivity implies no or little groundwater in this region. Low lineament density, high elevation, high drainage density, young basanites and steep slopes are also characteristic of this zone.

After the correlation between the groundwater potential zones map and the iso-conductivity map, additional validation was done by overlying existing springs, boreholes and static water level data on the GWPZ map (Fig. 9). Understanding the spatial distribution of existing springs and the static groundwater level is very important in understanding the groundwater resource distribution [13]. The reported static water level and spring data show a good correlation with the groundwater potential zones and iso-conductivity map. Springs are an expression of subterranean flow processes, which otherwise are accessible only through boreholes. Springs are common in volcanic terrains, particularly in lowland areas. All the springs in the study area occur in the older basanite formation at lower altitudes within the “very good” groundwater potential zones with high formation conductivity. Almost all the productive boreholes fall within the very good groundwater potential zone at lower elevation and gentle slope angles. From the groundwater potential zones map it is evident that all boreholes with static water level of 10–40 m are in the “very good” zone whereas the dry boreholes occur in the “moderate” zone with high resistivity (very low conductivity). Another notable observation is that as we move from the ‘‘very good’’ to ‘‘poor’’ zones, which also correspond from low elevation to high elevation, the static water level drops deeper.

This study displays the efficiency of employing an integrated approach of remote sensing, geophysical and existing springs and static water level data in prospecting for groundwater in volcanic terrain. Employing geospatial techniques, five groundwater potential zones were determined using six primary groundwater recharge controlling parameters that include lithology, elevation, slope, drainage density, lineament density, and land use. An iso-conductivity map used to verify the GWPZ map highlighted four groundwater potential zones, which are in strong correlation with that of the GWPZ map. More than half of the study area falls in the ‘good’ to ‘very good’ zones, indicating most of the study area has good groundwater potential. These results were validated by existing springs and static water level data overlay. The spatial distribution of the existing springs and the static water level is in close agreement with the GWPZ and iso-conductivity maps. The moderate groundwater potential zone, the very low conductivity zone, the absence of springs and the dry boreholes explain the scarcity in water supply being experienced in Bova, Bonakanda and Bokova. The final map can be used to increase the success ratio of installing groundwater boreholes and as a guide for decision-making regarding sustainable groundwater development and management.

This integrated approach for groundwater prospecting is recommended in both urban and remote areas, in order to obtain effective results. Furthermore, due to the higher cost of geophysical surveys, it is also recommended that remote sensing and geospatial techniques be used as a time- and cost-effective approach to select sites for field geophysical surveys before drilling boreholes.

Funding The authors declare that no funds, grants, scholarship or other supports were received during this research.

Competing Interests The authors have no competing interests to declare that are relevant to the content of this article.

Authors’ contributions All authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by Lionel Junior Ngome Njumbe. Lionel Junior Ngome Njumbe wrote the first draft of the manuscript while Anatole Eugene Djieto Lordon and Christopher M. Agyingi1 revised it critically for important intellectual content. All authors read and approved the final manuscript and agree to be accountable for all aspects of the work.

Data Availability Statement The Landsat 8 and SRTM data are available at https://earthexplorer.usgs.gov/. The shapefiles, borehole and geophysical data generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

W. Ott, “Access to Drinking Water and Stakeholder Action - Drinking Water Governance in Cameroon from a Political-Ecological Perspective Case Study: Upper Mefou Watershed, Cameroon,” Master Thesis, Freie Universität Berlin, 2014.
R. Emile, J. R. Clammer, P. Jayaswal, and P. Sharma, “Addressing water scarcity in developing country contexts: a socio-cultural approach,” Humanit Soc Sci Commun, vol. 9, no. 1, Dec. 2022, doi: 10.1057/s41599-022-01140-5.
A. Gadgil, “Drinking water in developing countries,” Annual Review of Energy and the Environment, vol. 23, no. 1, pp. 253–286, 1998, doi: 10.1146/annurev.energy.23.1.253.
D. N. Chakkaravarthy, “Water Scarcity- Challenging the Future,” International Journal of Agriculture Environment and Biotechnology, vol. 12, no. 3, Sep. 2019, doi: 10.30954/0974-1712.08.2019.2.
T. Gleeson, K. M. Befus, S. Jasechko, E. Luijendijk, and M. B. Cardenas, “The global volume and distribution of modern groundwater,” Nature Geoscience, vol. 9, no. 2. Nature Publishing Group, pp. 161–164, Feb. 01, 2016. doi: 10.1038/ngeo2590.
A. Richts, W. F. Struckmeier, and M. Zaepke, “WHYMAP and the Groundwater Resources Map of the World 1:25,000,000,” Sustaining Groundwater Resources, pp. 159–173, 2011, doi: 10.1007/978-90-481-3426-7_10.
Y. Fan, H. Li, and G. Miguez-Macho, “Global patterns of groundwater table depth,” Science (1979), vol. 339, no. 6122, pp. 940–943, Feb. 2013, doi: 10.1126/science.1229881.
C. Serele, A. Pérez-Hoyos, and F. Kayitakire, “Mapping of groundwater potential zones in the drought-prone areas of south Madagascar using geospatial techniques,” Geoscience Frontiers, vol. 11, no. 4, pp. 1403–1413, Jul. 2020, doi: 10.1016/j.gsf.2019.11.012.
O. Tafila, R. T. Ranganai, D. B. Moalafhi, K. K. Moreri, and J. G. Maphanyane, “Investigating groundwater recharge potential of Notwane catchment in Botswana using geophysical and geospatial tools,” J Hydrol Reg Stud, vol. 40, Apr. 2022, doi: 10.1016/j.ejrh.2022.101011.
A. Gaber, A. K. Mohamed, A. Elgalladi, M. Abdelkareem, A. M. Beshr, and M. Koch, “Mapping the groundwater potentiality of West Qena area, Egypt, using integrated remote sensing and hydro-geophysical techniques,” Remote Sens (Basel), vol. 12, no. 10, May 2020, doi: 10.3390/rs12101559.
K. S. Ishola, A. A. Fatoyinbo, A. I. Hamid-Mosaku, C. J. Okolie, O. E. Daramola, and T. O. Lawal, “Groundwater potential mapping in hard rock terrain using remote sensing, geospatial and aeromagnetic data,” Geosystems and Geoenvironment, vol. 2, no. 1, p. 100107, Feb. 2023, doi: 10.1016/j.geogeo.2022.100107.
R. Thapa, S. Gupta, S. Guin, and H. Kaur, “Assessment of groundwater potential zones using multi-influencing factor (MIF) and GIS: a case study from Birbhum district, West Bengal,” Appl Water Sci, vol. 7, no. 7, pp. 4117–4131, Nov. 2017, doi: 10.1007/s13201-017-0571-z.
R. Amer, M. Sultan, R. Ripperdan, A. Ghulam, and T. Kusky, “An integrated approach for groundwater potential zoning in shallow fracture zone aquifers,” Int J Remote Sens, vol. 34, no. 19, pp. 6539–6561, 2013, doi: 10.1080/01431161.2013.804221.
L. A. Magaia, T. Goto, A. A. Masoud, and K. Koike, “Identifying Groundwater Potential in Crystalline Basement Rocks Using Remote Sensing and Electromagnetic Sounding Techniques in Central Western Mozambique,” Natural Resources Research, pp. 1–14, 2017, doi: 10.1007/s11053-017-9360-5.
F. Z. Echogdali et al., “Application of Analytical Hierarchy Process and Geophysical Method for Groundwater Potential Mapping in the Tata Basin, Morocco,” Water (Switzerland), vol. 14, no. 15, Aug. 2022, doi: 10.3390/w14152393.
R. K. Prasad, N. C. Mondal, P. Banerjee, M. V. Nandakumar, and V. S. Singh, “Deciphering potential groundwater zone in hard rock through the application of GIS,” Environmental Geology, vol. 55, no. 3, pp. 467–475, Aug. 2008, doi: 10.1007/s00254-007-0992-3.
N. S. Magesh, N. Chandrasekar, and J. P. Soundranayagam, “Delineation of groundwater potential zones in Theni district, Tamil Nadu, using remote sensing, GIS and MIF techniques,” Geoscience Frontiers, vol. 3, no. 2, pp. 189–196, 2012, doi: 10.1016/j.gsf.2011.10.007.
A. Shebl, M. I. Abdelaziz, H. Ghazala, S. A. S. Araffa, M. Abdellatif, and Á. Csámer, “Multi-criteria ground water potentiality mapping utilizing remote sensing and geophysical data: A case study within Sinai Peninsula, Egypt,” Egyptian Journal of Remote Sensing and Space Science, vol. 25, no. 3, pp. 765–778, Dec. 2022, doi: 10.1016/j.ejrs.2022.07.002.
M. T. Yossa, A. E. D. Lordon, C. M. Agyingi, J. Agbor-Taku, Y. N. Shandini, and C. E. Bessong, “Remote sensing and geographic information system (GIS)-based high-resolution mapping of potential groundwater recharge zones on the hard rock terrains of the Cameroon volcanic line (CVL),” SN Appl Sci, vol. 5, no. 1, Jan. 2023, doi: 10.1007/s42452-022-05248-w.
M. S. Njome, C. E. Suh, R. S. J. Sparks, S. N. Ayonghe, and J. G. Fitton, “The Mount Cameroon 1959 compound lava flow field: Morphology, petrography and geochemistry,” Swiss J Geosci, vol. 101, no. 1, pp. 85–98, May 2008, doi: 10.1007/s00015-007-1245-x.
J. C. Olivry, “fleuves et rivières du cameroun,” pp. 15–17, 1986, Accessed: Mar. 13, 2023. [Online]. Available: https://books.google.com/books/about/Fleuves_et_rivi%C3%A8res_du_Cameroun.html?hl=fr&id=XqJWAAAAMAAJ
E. L. Molua and C. M. Lambi, “Climate hydrology and water resources in Cameroon,” p. 37, 2006, doi: 10.1007/978-94-009-0207-7_1.
C. E. Suh et al., “The 1999 and 2000 eruptions of Mount Cameroon: Eruption behaviour and petrochemistry of lava,” Bull Volcanol, vol. 65, no. 4, pp. 267–281, 2003, doi: 10.1007/s00445-002-0257-7.
B. Ateba et al., “Eruptive and earthquake activities related to the 2000 eruption of Mount Cameroon volcano (West Africa),” Journal of Volcanology and Geothermal Research, vol. 179, no. 3–4, pp. 206–216, Jan. 2009, doi: 10.1016/j.jvolgeores.2008.11.021.
J. D. Fairhead and C. S. Okereke, “A regional gravity study of the West African rift system in Nigeria and Cameroon and its tectonic interpretation,” Tectonophysics, vol. 143, no. 1–3, pp. 141–159, Nov. 1987, doi: 10.1016/0040-1951(87)90084-9.
G. W. Stuart, J. D. Fairhead, L. Dorbath, and C. Dorbath, “A seismic refraction study of the crustal structure associated with the Adamawa Plateau and Garoua Rift, Cameroon, West Africa,” Geophysical Journal of the Royal Astronomical Society, vol. 81, no. 1, pp. 1–12, 1985, doi: 10.1111/J.1365-246X.1985.TB01346.X.
Y. H. Poudjom Djomani, J. M. Nnange, M. Diament, C. J. Ebinger, and J. D. Fairhead, “Effective elastic thickness and crustal thickness variations in west central Africa inferred from gravity data,” J Geophys Res, vol. 100, no. B11, 1995, doi: 10.1029/95jb01149.
J. G. Fitton, “The Cameroon line, West Africa: A comparison between oceanic and continental alkaline volcanism,” Geol Soc Spec Publ, vol. 30, no. May 2014, pp. 273–291, 1987, doi: 10.1144/GSL.SP.1987.030.01.13.
B. Ateba and N. Ntepe, “Post-eruptive seismic activity of Mount Cameroon (Cameroon) West Africa: a statistical analysis,” 1997.
J.-P. F. Tsafack, P. Wandji, J.-M. Bardintzeff, H. Bellon, H. Guillou, and J. Tsafack, “The Mount Cameroon stratovolcano (Cameroon Volcanic Line, Central Africa): Petrology, geochemistry, isotope and age data,” vol. 47, pp. 65–78, 2009.
“EarthExplorer.” https://earthexplorer.usgs.gov/ (accessed Mar. 15, 2023).
R. P. Gupta, Remote Sensing Geology, Second Edition. New York: Springer-Verlag Berlin Heidelberg, 2003. doi: 10.1007/978-3-662-05283-9.
E. Farahbakhsh et al., “Computer vision-based framework for extracting tectonic lineaments from optical remote sensing data,” Int J Remote Sens, vol. 41, no. 5, pp. 1760–1787, Mar. 2020, doi: 10.1080/01431161.2019.1674462.
“Getting Started with ENVI Tutorial.” https://envi.geoscene.cn/help/Subsystems/envi/Content/Tutorials/Tools/GettingStartedWithENVITutorial.htm (accessed Mar. 19, 2023).
Y. C. Chang, G. S. Song, and S. K. Hsu, “Automatic extraction of ridge and valley axes using the profile recognition and polygon-breaking algorithm,” Comput Geosci, vol. 24, no. 1, pp. 83–93, Jan. 1998, doi: 10.1016/S0098-3004(97)00078-2.
T. L. Saaty, Analytic Hierarchy Process Planning, Priority Setting, Resource Allocation. McGraw-Hill, Inc, 1980.
K. Goepel, “Implementation of an Online software tool for the Analytic Hierarchy Process (AHP-OS),” International Journal of the Analytic Hierarchy Process, vol. 10, no. 3, pp. 469–487, 2018, doi: 10.13033/ijahp.v10i3.590.
R. Peng, J. X. Wei, B. R. Di, P. B. Ding, and Z. C. Liu, “Experimental research on seismoelectric effects in sandstone,” Applied Geophysics, vol. 13, no. 3, pp. 425–436, 2016, doi: 10.1007/s11770-016-0570-0.
P. K. Ghosh, S. Bandyopadhyay, and N. C. Jana, “Mapping of groundwater potential zones in hard rock terrain using geoinformatics: a case of Kumari watershed in western part of West Bengal,” Model Earth Syst Environ, vol. 2, no. 1, Mar. 2016, doi: 10.1007/s40808-015-0044-z.
D. W. O’leary, J. D. Friedman, and H. A. Pohn, “Lineament, linear, lineation: Some proposed new standards for old terms,” Bulletin of the Geological Society of America, vol. 87, no. 10, pp. 1463–1469, 1976, doi: 10.1130/0016-7606(1976)87<1463:LLLSPN>2.0.CO;2.
G. E. Tucker, F. Catani, A. Rinaldo, and R. L. Bras, “Statistical analysis of drainage density from digital terrain data,” Geomorphology, vol. 36, no. 3–4, pp. 187–202, 2001, doi: 10.1016/S0169-555X(00)00056-8.
D. Greenbaum, “Review of remote sensing applications to groundwater exploration in basement and regolith,” Nottingham NG12, Mar. 1985. Accessed: Apr. 01, 2023. [Online]. Available: https://nora.nerc.ac.uk/id/eprint/505150/1/WC_OG_85_1.pdf
G. E. Omolaiye et al., “Integration of remote sensing, GIS and 2D resistivity methods in groundwater development,” Appl Water Sci, vol. 10, no. 6, Jun. 2020, doi: 10.1007/s13201-020-01219-x.

No competing interests reported.

Download PDF

Editorial decision: Major revision
08 Jul, 2023
Reviews received at journal
07 Jul, 2023
Reviews received at journal
12 Jun, 2023
Reviewers agreed at journal
05 Jun, 2023
Reviewers agreed at journal
02 Jun, 2023
Reviewers invited by journal
02 Jun, 2023
Editor assigned by journal
01 Jun, 2023
Submission checks completed at journal
01 Jun, 2023
First submitted to journal
28 May, 2023

You are reading this latest preprint version

Determination of Groundwater Potential Zones on the Eastern Slope of Mount Cameroon using Geospatial Techniques and Seismoelectric Method

Status:

Version 1

Abstract

Figures

Article Highlights

1. Introduction

2. Location and Geology

3. Materials and Methods

3.1. Geospatial techniques

3.2. Analytic Hierarchy Process (AHP)

3.3. Verification and Validation

4. Results and Interpretation

4.1. Factors Influencing Groundwater Recharge

4.1.1. Lithology

4.1.2. Elevation

4.1.3. Slope

4.1.4. Land use

4.1.5. Lineament density

4.1.6. Drainage density

4.2. Groundwater potential zonation

4.3 Model validation

4.3.1. Seismoelectric survey

4.5.2. Existing springs and static water level

5. Discussion

5.1. Groundwater potential zones

5.2. Validation of the groundwater potential zones map

6. Conclusions

Declarations

References

Additional Declarations

Status:

Version 1