Integrated Geospatial and Geophysical Approaches for Mapping Groundwater Potential in the Semi-Arid Bukombe District, Tanzania

: The rapid growth of civil societies coupled with population inﬂux due to the artisanal mining industry in the Bukombe district (BD) has triggered a high demand for water resources. The daily consumption of water resources in the district surpasses the supply from available surface water sources. Thus, the situation has raised the demand for groundwater resources as an alternative. Despite the importance of groundwater resources, no current studies have spatially assessed groundwater potential to locate optimal points for borehole development. This study intended to investigate and map the groundwater potential areas (GWPAs) in the semi-arid BD using remote sensing (RS), the geographic information system (GIS), and the analytic hierarchy process (AHP) to help local communities access clean and safe water. Rainfall, geology, slope, drainage density, land use/land cover and lineament density were prepared to delineate the map of GWPAs. The map was categorized into poor (0.21%), moderate good (51.39%), good (45.70%) and very good (2.70%). Finally, the GWPA map was validated using Vertical Electrical Sounding (VES), 2-D sections and a drilled borehole. The validation results conﬁrmed that the applied approach provides signiﬁcant results that can help in planning the sustainable utilization of groundwater resources.


Introduction
Globally, water is recognized as a basic irreplaceable need of human beings. Its accessibility, however, is a major challenge, especially in developing countries [1], because of the dwindling surface water sources due to climate change and anthropogenic pollution [2][3][4][5]. In recent years, climate change, rapid population growth and surface water losses within the conveying pipes during transportation from sources to storage reservoirs have made groundwater the most demanded resource [6,7]. Moreover, some communities perceive the quality of groundwater to be better and more reliable than other sources of water for domestic uses [6]. Thus, it is important to locate dependable sources of groundwater supply in all climatic regions across the world [8].
Demographic data show that about 11% of the world's population, particularly in Sub-Saharan Africa and Oceania, do not have access to clean and safe water [9]. The majority have aligned themselves to an alternative and cost-effective source, which is groundwater [10]. The lack of adequate hydrogeological and updated geospatial data in place to identify groundwater potential areas has led to the drilling of boreholes and the digging of wells in areas where groundwater potentiality is suboptimal [11][12][13].
Economic water scarcity is a common problem in developing counties such as Tanzania [3,4]. About 25% of Tanzania's population depends on groundwater sources for domestic and Earth 2023, 4 242 other uses [14]. In spite the significance of this concern, groundwater resources are not thoroughly explored in terms of aquifer extent, aquifer characteristics and water quality to sufficiently cater to users [15,16], and the most affected areas where groundwater is poorly explored is the western part of Tanzania. As a result, individuals, communities and institutions have decided to pass on important indigenous knowledge on exploiting and managing shallow aquifers [17]. Due to current climate variability and change, high levels of resource exploitation, increased urbanization and population growth, the determination of groundwater potential areas in developing countries, particularly in semi-arid areas, is inevitable [18]. However, this is challenging as it needs reliable information for locating the boreholes within a groundwater potential area [18][19][20].
In semi-arid region such as the Bukombe district in western Tanzania, groundwater is the major source of water for domestic uses [21]. Due to water scarcity, groundwater projects initiated by the government and NGOs to narrow the water supply gap experienced failure a few years after inauguration. This was largely attributed to the paucity of technical knowledge in the field of satellite data use that hampered the assessment of tangible groundwater potential areas [3]. Research conducted by SMEC (2015) identified Bukombe as one of the areas in the Malagarasi catchment with severe groundwater risk. Furthermore, their study revealed that Bukombe is the area where groundwater recharge takes place and is proposed to be legally protected to safeguard not only the welfare of the people but also the Malagarasi-Muyovosi Ramsar Site [15]. However, the spatial distribution of the groundwater potentiality, aquifer characteristics and groundwater quality are vaguely known.
Therefore, this research intended to assess groundwater potential areas in the Bukombe semi-arid area using geospatial and geophysical approaches which digitally demarcated and categorized the area for further groundwater development. Our findings agree with results from previous researchers which demonstrated that RS and GIS provide potentially powerful tools to study groundwater resources [22][23][24].

Location and Accessibility
Bukombe is one of the districts in the Geita region, northwest Tanzania, covering an area of about 8000 sq. km between latitudes 3.11-4.47 • S and longitudes 31.15-32.167 • E ( Figure 1). The district is bordered to the south by the Kaliua district, to the west by the Kibondo district and to the north by Biharamulo, Chato and Geita districts; in the eastern part, it is bordered by Mbogwe and Kahama districts. A large area in the southern part falls under the Kigosi National Park (Figure 1). The study area is accessible through the tarmac road of the Isaka-Rusahunga highway in almost all town centers. Most of the villages are accessible through rough roads in both wet and dry seasons.
In 2012, the population of the Bukombe district was 224,542 [25]. It is projected to be 333,934 in 2025 using the population projection model as per Equation (1).
where P f , P p , gr and t are the future population, previous population, growth rate and time in years, respectively (https://www.youtube.com/watch?v=ygMCKsb4aAU, accessed on 25 January 2023). The land use in the Bukombe district includes the forest where there is the Kigosi National Park and the community and individual forest reserves [26], along with agricultural land, especially paddy pans and corn cultivated areas. Settlement areas ranging from towns to village housing are also present. The curve number for town centers in general is higher than for villages as a result of accelerating surface runoff. The mining industry also accounts for a small portion of land use in the district. Water demand in the Bukombe district is not explicitly displayed but rather could be inferred from the water demand, estimated on the Malagarasi catchment of which the Bukombe district is a part. SMEC (2015) estimated the water demand in the Malagarasi catchment as 1299 MCM/year in 2015, 1437 MCM/year in 2025 and 1704 MCM/year in 2035 [15]. With a simple areal-and population-based ration, Bukombe's water demand can be downscaled to 208 MCM/year in 2025 that will be used for domestic purposes, irrigation, tourism, wildlife, mining, industry, livestock keeping and recreation. According to BDAR 2021, the district is planning to have 5 recreation sites and 2 sports grounds; develop 1850 ha for irrigation in Bugelenga, Mjimwema and Nyampangwe, where 250 ha and 100 ha in Bugelenga and Mjimwema, respectively, are locally being irrigated; and raise the water supply coverage that will cater to over 130 lodges, industries, mining activities, livestock, etc. [26].

Climate
The climate of the study area is that of a tropical savannah with distinct rainy and dry seasons [16]. The district experiences a unimodal rain pattern with total annual rainfall ranging between 600 and 1000 mm [15]. Rainfall occurs between October and May with a short break in January or February [27]. The rest of the months are the dry season [15]. The temperature ranges between 20 and 35 • C, with annual evapotranspiration of 1600 mm per year.

Hydrogeological Settings
Based on Tanzania's water basins, the Bukombe administrative district falls under the Lake Tanganyika Basin, the second largest basin after the Rufiji Basin, with an area of 160,426 sq. kilometers (Figure 2a,b). BD lies at the upper Malagarasi catchment consisting of three subcatchments made of the Mlera, Nikonga and Kigosi rivers [16]. It is dominated by NE-SW, N-S and NW-SE lineaments that might have led to the formation of the Mlera and Nikonga rivers in the middle, the Moyowosi river in the west and the Kigosi stream in the east (Figure 1). The Moyowosi and Nikonga rivers are perennial, whereas the Mlera and Kigosi rivers are ephemeral in their upstream but become active as they approach the Earth 2023, 4 244 confluence with the Moyowosi river. Groundwater potentiality is mainly dependent on geological structures and alluvial basins ( Figure 3). The available boreholes in the study area are located in an alluvial basin ( Figure 3).
Based on Tanzania's water basins, the Bukombe administrative district falls under the Lake Tanganyika Basin, the second largest basin after the Rufiji Basin, with an area of 160,426 sq. kilometers (Figure 2a,b). BD lies at the upper Malagarasi catchment consisting of three subcatchments made of the Mlera, Nikonga and Kigosi rivers [16]. It is dominated by NE-SW, N-S and NW-SE lineaments that might have led to the formation of the Mlera and Nikonga rivers in the middle, the Moyowosi river in the west and the Kigosi stream in the east (Figure 1). The Moyowosi and Nikonga rivers are perennial, whereas the Mlera and Kigosi rivers are ephemeral in their upstream but become active as they approach the confluence with the Moyowosi river. Groundwater potentiality is mainly dependent on geological structures and alluvial basins ( Figure 3). The available boreholes in the study area are located in an alluvial basin (Figure 3).

Geological Setting
The Bukombe district is part of the Lake Victoria Gold Field (LVGF). LVGF is the area comprising large gold mines and has been documented as an important site for gold mining and exploration in East Africa [28,29]. Rocks within the gold fields have been subdivided into the Nyanzian System and the Kavirondian System. The Nyanzian System comprises an assemblage of mafic metavolcanics, chemical sediments (oxide and silicate facies iron formation) and lesser epiclastic sediments [30]. Minor ultramafic and felsic volcanics plus volcano clastics form basal and capping members, respectively, of the largely tholeitic volcanic pile. These Nyanzian rocks are unconformably overlain by the Kavirondian System, which consists of epiclastic sediments and minor tuffs. Both systems are intruded by large granitic plutons. The Banded iron formation and Banded ferruginous Earth 2023, 4, FOR PEER REVIEW 5 chert are narrow zones that overlie the basic volcanic piles [30]. Visual observations by the researcher confirmed the weathered sedimentary rock locally known as the laterites in both Katente and Bulangwa wards ( Figure 3).

Types of Data and Data Collection Procedures
Geospatial data collection and processing included the satellite image and digital elevation model, the total annual rainfall (TAR) data and the geological map. Fieldwork included a geophysical survey for subsurface lithological analysis.  The Bukombe district is part of the Lake Victoria Gold Field (LVGF). LVGF is the area comprising large gold mines and has been documented as an important site for gold mining and exploration in East Africa [28,29]. Rocks within the gold fields have been subdivided into the Nyanzian System and the Kavirondian System. The Nyanzian System comprises an assemblage of mafic metavolcanics, chemical sediments (oxide and silicate facies iron formation) and lesser epiclastic sediments [30]. Minor ultramafic and felsic volcanics plus volcano clastics form basal and capping members, respectively, of the largely tholeitic volcanic pile. These Nyanzian rocks are unconformably overlain by the Kavirondian System, which consists of epiclastic sediments and minor tuffs. Both systems are intruded by large granitic plutons. The Banded iron formation and Banded ferruginous chert are narrow zones that overlie the basic volcanic piles [30]. Visual observations by the researcher confirmed the weathered sedimentary rock locally known as the laterites in both Katente and Bulangwa wards ( Figure 3).

Types of Data and Data Collection Procedures
Geospatial data collection and processing included the satellite image and digital elevation model, the total annual rainfall (TAR) data and the geological map. Fieldwork included a geophysical survey for subsurface lithological analysis.
The Digital Elevation Model (DEM) was downloaded from the Earth Explorer website (Shuttle Radar Topography Mission (SRTM-DEM)) of 30 m resolution (https://earthexplorer. usgs.gov, (accessed on 25 January 2023)). The satellite images (20% cloud cover) were downloaded from Landsat8, particularly using the OLI sensor data. Three tiles were needed, 171-063, 172-062 and the 171-062, where the first and last numbers represent the path and row of the satellite.
Conventional data were individually handled in accordance with their respective standard procedures. Rainfall data, for example, from 21 rainfall stations were collected from Lake Tanganyika Basin Water Board, Kahama Sub-Office and remotely from Lake Victoria Basin Water Board. The geology maps from the Tanzania Geological Survey were collected in picture mode quarter degree sheets 44, 45, 60, 61, 76 and 77.
Geophysical data were collected in two phases; the first phase was 2-D and the second phase was 1-D. In the 2-D geophysical survey, a 32-channel McOhm Profile-4 resistivity system utilizing a pole dipole electrode array type was used to acquire profiles each 280 m long with an electrode interval of 10 m. Pole dipole array enables one to obtain the quantitative depth, lithological information and resistivity information below the measurement station [31]. Fourteen profiles were measured at the locations of Ilyamchele, Bulangwa, Katente, Shikalibuga, Roman Catholic Mission and Runzewe.
The vertical electrical sounding (1-D) resistivity measurements were carried out in the anomalous points interpreted from the 2-D pseudo-sections. It utilized ABEM Terrameter LS with Schlumberger configuration at the minimum and maximum current electrode spacing (AB) of about 3 and 600 m, respectively, which is suitable for a subsurface investigation depth of about 300 m [13,32]. A total of 80 VES stations were included in the study area and out of the 2-D profiles aiming at widening the scope of comparing the two approaches.

Geospatial Data Analysis
The preparation of thematic maps and later the groundwater potential area map was entirely dependent on GIS software. A total of six thematic layers were generated from DEM, satellite images, climatological data and geology as described below. The DEM from the Shuttle Radar Topography Mission generated two thematic maps, one of the drainage density and one of the slope (Figure 4).
Drainage density was generated by following sequential algorithms: fill, flow direction, flow accumulation, threshold dataset, stream order and stream to feature, then drainage density. The slope is a measure of elevation change and termed as a topographic parameter [33]. The slope was produced in GIS through the conversion of DEM into a slope dataset and reclassification into five classes according to Van Zuidam, who divided the slope as follows: 0-2% flat, 2-15% moderate, 15-25% moderate steep, 25-40% steep and greater than 40% very steep [33]. The rainfall map of the study area was prepared by interpolating the TAR from 21 stations (Appendix C). The data were processed using GIS and interpolated using Kriging interpolation to produce the rainfall map ( Figure 5b). measurement station [31]. Fourteen profiles were measured at the locations of Ilyamchele, Bulangwa, Katente, Shikalibuga, Roman Catholic Mission and Runzewe.
The vertical electrical sounding (1 -D) resistivity measurements were carried out in the anomalous points interpreted from the 2 -D pseudo-sections. It utilized ABEM Terrameter LS with Schlumberger configuration at the minimum and maximum current electrode spacing (AB) of about 3 and 600 m, respectively, which is suitable for a subsurface investigation depth of about 300 m [13,32]. A total of 80 VES stations were included in the study area and out of the 2 -D profiles aiming at widening the scope of comparing the two approaches.

Geospatial Data Analysis
The preparation of thematic maps and later the groundwater potential area map was entirely dependent on GIS software. A total of six thematic layers were generated from DEM, satellite images, climatological data and geology as described below. The DEM from the Shuttle Radar Topography Mission generated two thematic maps, one of the drainage density and one of the slope (Figure 4). Drainage density was generated by following sequential algorithms: fill, flow direction, flow accumulation, threshold dataset, stream order and stream to feature, then drain- Earth 2023, 4, FOR PEER REVIEW 7 age density. The slope is a measure of elevation change and termed as a topographic parameter [33]. The slope was produced in GIS through the conversion of DEM into a slope dataset and reclassification into five classes according to Van Zuidam, who divided the slope as follows: 0-2% flat, 2-15% moderate, 15-25% moderate steep, 25-40% steep and greater than 40% very steep [33]. The rainfall map of the study area was prepared by interpolating the TAR from 21 stations (Appendix C). The data were processed using GIS and interpolated using Kriging interpolation to produce the rainfall map ( Figure 5b). Three satellite images were used to classify the study area into different types of land use and land cover (LULC). The procedures and steps undertaken were image enhancement through dark object subtraction DOS, conversion into false color composite FCC,  Three satellite images were used to classify the study area into different types of land use and land cover (LULC). The procedures and steps undertaken were image enhancement through dark object subtraction DOS, conversion into false color composite FCC, mosaicking and clipping into study area jurisdiction. Training samples' definition was the next step after clipping, during which macro and micro classes were developed [34]. Finally, the five LULCs (water body, forest, built up, agricultural land and the bare land) were generated through interactive supervised classification.
The lineament density map was generated from mosaicked satellite image in QGIS. A single red band was selected, from which hill shades of azimuth angles 45 • , 135 • , 225 • and 315 • were developed. A polyline layer was also prepared, onto which lineaments were manually digitized from the hill shade layers at a 1:100,000 scale [35]. Cultural and/or manmade linear features such as roads, canals and rail lines were eliminated by projecting the lineament layer onto Google Earth to ensure that all the digitized lineaments are purely geological features [35]. The final lineament density map was created.
The geological maps from the Tanzania Geological Survey (QDS 44, 45, 60, 61, 76 and 77), the Lake Tanganyika Basin geological map from SMEC report and the LVGF map were used to extract geological formations in the study area [15,30]. A total of six geological formations were identified ( Figure 3): basalt, recent sediments (limestone, calcrete, gravel, sand, silt and clay), alluvial deposit, granites, metasediments of Archean age (phyllite, schist, gneiss, amphibolite and marble) and the Nyanzian Kavirondian in the form of banded iron formation.
Thematic layers were then input into AHP excel software for pairwise comparison weight calculation. At this stage, each thematic layer was compared to other layers one at a time to choose which of the two is more important than the other with regard to groundwater potentiality. The choices made assigned ranks to themes and their classes from 1 (equally important) to 9 (extremely important) depending on how important a factor is over the others, as per Saaty's (1980) classification [10,36] as Table 1 elaborates. Of equal importance, the consistency ratio was monitored at all times, making sure it was an acceptable value of at most ten percent (see Appendix A) [12]. Table 1. Preference scale for two parameters in AHP (Saaty 1980).

Scale
Intensity of Importance Explanation 1 Equally Two activities contribute equally to the objective 3 Moderately Experience and judgment slightly to moderately favor one activity over another 5 Strongly Experience and judgment strongly or essentially favor one activity over another 7 Very strongly An activity is strongly favored over another and its dominance is shown in practice 9 Extremely The evidence of favoring one activity over another is of the highest degree possible of an affirmation 2,4,6,8 Intermediate values between two judgments Used to represent compromises between the preferences in weights 1, 3, 5, 7 and 9 Ranks for aggregates of each theme were set based on the literature, field experience and groundwater science expertise [37]. The thematic layers in raster format were overlaid using the weighted overlay command as per the following equation to generate groundwater potential area (GWPA) [12,21,38]: where w tm is the weight of the thematic factor and r tm is the scale value or rank each aggregate is given by the researcher. Detailed clarification on weight values and ranks for this study are provided in Appendices A and E. The generated GWPA was used to select potential areas for conducting geophysical measurements. The selection of the traverses during the resistivity survey was based Earth 2023, 4 248 on geomorphological settings, indigenous historical narration and the existence of local water sources, aiming at capturing key areas that seemed to be prominent for groundwater development. The 2-D survey was conducted first to map the subsurface geology using Marcohm Profile 4 machine followed by the 1-D Terrameter Ls on promising points for confirming and locating drilling sites (Figure 4).

Geophysical Data Analysis
Modeling and inversion of the 1-D resistivity data were performed using IP2WIN/DCINV software to generate resistivity models with charts and tables. The 2-D resistivity data were analyzed using RES2-DINV software to produce sections' interpretation of geological layers in relation to groundwater occurrence. The final output was the identification of the number of aquifers, depth to the aquifer, thickness of the aquifer and resistivity of the aquifer [39].

Results
This section portrays two main types of results: the geospatial results and the geophysical results. The geospatial part includes thematic layers which were integrated to produce the GWPA map. The geophysical part accounts for both VES data and 2-D data in terms of resistivity charts and resistivity sections, respectively. The latter part was carried out for the sake of confirming the validity of the former part's results. Lastly, the general conclusion was made that delineated groundwater potential areas can be used to undertake further water projects.

Slope
Geomorphology in any particular area governs the movement of water. Steep slopes promote runoff, while topographical depressions favor infiltration [21]. The slope, like other themes, plays a significant role in groundwater recharge leading to groundwater potentiality. With reference to Van Zuidam's slope classification, BD has only two classes: flat and moderate. However, the slopes in the study area were categorized into five classes: (i) very low slope (<1%), (ii) low slope (1 to 3%), (iii) medium slope (3 to 5%), (iv) high slope (5 to 10%) and (iv) very high slope (>10%). With reference to Figure 5a, a large area falls under low and very low classes, accounting for more than 60% of the entire area. This means that a large amount of water from rain is given ample time to infiltrate, hence recharging the aquifer and increasing the groundwater potentiality, if and only if other factors remain constant.

Rainfall
Rainfall plays major role in recharging groundwater. Therefore, areas with higher rainfall are generally considered as potential zones for groundwater exploitation. The rainfall map of the study area (Figure 5b) was divided into five classes: (1) low rainfall zone with rainfall less than 1043 mm, (2) moderate rainfall zone with rainfall between 1043 and 1068 mm, (3) good rainfall zone with rainfall between 1068 and 1091 mm, (4) very good rainfall zone with rainfall between 1091 and 1123 mm and (5) excellent rainfall zone with rainfall more than 1123 mm. The northeast part receives the least amount of rainfall; the precipitation sequentially increases through to the far southwest end, which has an excellent amount of rainfall (Figure 5b). This is attributed to the Kigosi National Park, which accounts for a large portion of the district. The range of rainfall is not that substantial, and areal percent for all classes are 23.95%, 24.39%, 22.10%, 16.57% and 12.99%, respectively.

Land Use/Land Cover
Land use/land cover of an area is largely governed by the groundwater resources and, at the same time, plays an important role in controlling these resources. It influences Earth 2023, 4 249 many hydrogeological processes in the water cycle such as evapotranspiration, infiltration, surface runoff, etc. [40,41]. In the forest and agricultural land, runoff is generally less prevalent and infiltration is more common, whereas, in settlement areas and other built up parts such as tarmac roads, the rate of infiltration usually decreases due to the high curve number possessed by such areas.
Therefore, agricultural land and vegetation cover area were ranked higher than builtup area [38]. In the Bukombe district in particular (Figure 6a), a major portion of land use is forest, covering 5984.133 sq. km (74.3%), along with bare land covering 1456.541 sq. km (18.1%), agricultural land covering 387.404 sq. km (4.8%) and water and settlements each covering 1.4% with areas of 116.506 sq. km and 109.45 sq. km, respectively. The correctness of these results was checked by creating the classification confusion matrix, as displayed in Table 2. Recent research pointed out that the minimum acceptable accuracy assessment is 70% [42]. The accuracy assessment of this research is more than sufficient, as given in Table 2.  The quantification of drainage density and its types leads to runoff, infiltration, relief and permeability information acquisition [38]. It also reflects the proximity of spacing of  The quantification of drainage density and its types leads to runoff, infiltration, relief and permeability information acquisition [38]. It also reflects the proximity of spacing of streams as well as surface characteristics [43]. Further, drainage patterns give information Earth 2023, 4 250 related to surface materials and subsurface formation. For instance, dendritic drainage indicates homogenous rocks, whereas trellis, rectangular and parallel drainage patterns are suggestive of structural and lithological controls [38]. The Bukombe district in particular is characterized by a dendritic drainage pattern. Geospatial analysis categorized five classes of drainage density: very low, low, medium, high and very high, with areal percentages of 58.2%, 15.6%, 17.9%, 6.7% and 1.6%, respectively (Figure 6b). With these data, Bukombe can be presumed to be gently sloping to flat land that promotes more infiltration. It is of a lower drainage density, which encourages more groundwater potential than high-drainagedensity regions [22,44]. Low drainage density yields a coarse drainage texture and high drainage density leads to a fine drainage texture [44]. So, drainage density has an inverse effect on permeability. Therefore, it is an important feature in evaluating groundwater potential zones.

Lineament Density
The linear faults, accompanied by the fissure, provide space for the occurrence of groundwater [45]. In the stratum with the same lithology, the intersection of the faults leads to the development of the fissure, which tends to be the groundwater enrichment zone with the connectivity enhancement [6]. The linear structures extracted from Landsat 8 OLI images with the aid of the GeoTrace plugin in QGIS led to the formation of five lineament density classes in the study area (Figure 7a). The very low class has an area of 3358.585 sq. km (41.70%), the low lineament has an area of 1184.481 sq. km (14.71%), the medium lineament class has an area of 1401.553 sq. km (17.40%), the high lineament class has an area of 1605.99 sq. km (19.94%) and the very high lineament class has an area of 503.439 sq. km (6.25%). With these lineament data, structural controlled formations in the Bukombe district are limited.
Earth 2023, 4, FOR PEER REVIEW 11 has an inverse effect on permeability. Therefore, it is an important feature in evaluating groundwater potential zones.

Lineament Density
The linear faults, accompanied by the fissure, provide space for the occurrence of groundwater [45]. In the stratum with the same lithology, the intersection of the faults leads to the development of the fissure, which tends to be the groundwater enrichment zone with the connectivity enhancement [6]. The linear structures extracted from Landsat 8 OLI images with the aid of the GeoTrace plugin in QGIS led to the formation of five lineament density classes in the study area (Figure 7a). The very low class has an area of 3358.585 sq. km (41.70%), the low lineament has an area of 1184.481 sq. km (14.71%), the medium lineament class has an area of 1401.553 sq. km (17.40%), the high lineament class has an area of 1605.99 sq. km (19.94%) and the very high lineament class has an area of 503.439 sq. km (6.25%). With these lineament data, structural controlled formations in the Bukombe district are limited.  Geological formations are vessel-like features that hold groundwater in their pore spaces. The knowledge of how these earth materials formed and the evolution they have undergone opens understanding of the distribution of geologic parameters such as permeability and porosity [46]. Higher porosity contributes to higher groundwater storage and higher permeability contributes to higher groundwater yields [47].
The study area comprises a variety of geological formations (Figure 3), including metasediments of Archean age with areal coverage of 7.6%, and the granites and syenites family (41.72%). They are not intensively weathered, which hinders the infiltration of meteoric water, leading to low groundwater recharge. Recent terrestrial sediment accounts for 45.61% of the total area. Such areas have good hydraulic conductivity due to the presence of loose and coarse-grained material. The Nyanzian accounts 3.16% of the total area. The volcanic lava in the northwest corner of the study area has an areal extent of 1.92%. Outcrops of banded iron formation currently serve most of the shallow dug wells and springs for various uses. Groundwater occurrence is higher in metasediments ( Figure 7b) and the least in granites due to their aforementioned properties. The degree of metamorphism and weathering, the presence of geologic structures such as folds, nonconformability and faults escalate significant storage of water.

Groundwater Potential Map
The groundwater potential map area is presented in Figure 7b. Groundwater potentiality was demarcated, ending up with four areas (Table 3): poor, moderate, good and very good. The area with very good potential covers 215.56 sq. km, equivalent to 2.70% of the total district area. The spatial location of this area is mostly in the Kigosi National Park, with a few portions in the western part of Runzewe and Uyovu wards, as per Figure 7b. The good category lies mostly in the forest zone and partly in Runzewe west and Runzewe east wards. The potentiality demarcation by the software concurs with the true situation found on the ground, that the aforementioned areas have better groundwater potentiality than the eastern part of the district. The area with good potential covers 3649.70 sq. km, equivalent to 45.70%. The moderate groundwater potential zone is the largest area of all, occupying most of the residential area in the east, north and west. It covers an area of 4104.07 sq. km, corresponding to 51.39% of the total area. The poor groundwater potential zone is the smallest category, with an area of 16.74 sq. km, equal to 0.21% of the total area, and it appears as a spot on the map (Figure 7b). This area lies in the granitic zone, where little water is able to percolate into ground aquifers. The central to western part corresponds to alluvial plains and lacustrine sediments, which coincide with the low-slope and highlineament areas.

Vertical Electrical Sounding Results
The inversion of the 80 VES measurements was calibrated using control VES data that were taken from existing boreholes. The results and the inferred interpretative curves (Table 4 and Appendix D) reveal that the subsurface structure of the study area comprises a complex sequence of lithological layers. The number and the sequence of different layers vary from one site to another (Figure 8), reflecting the lithological complexity and heterogeneity. The succession of lithological layers and the vertical distribution of resistivity values, as shown in interpretative curves, allow us to distinguish two patterns of subsurface structures, namely the weathering profile of Precambrian basement rocks and the presence of alluvial deposits. The sequence of interpreted lithological layers illustrated that the subsurface is composed of three to seven lithological layers with a generally high-low-highest depth trend of resistivity values (Appendix D).   The VES measurements undertaken in the study area were dominated by several curve types, of which 40% belong to KH type, 19% fall under HA type, 15% were categorized as H type and 6% are of QH type. AA, KHKH, QHKH, HAA, AKH, QHKA, HKHA, KAA, KHA, QKH, QH and HAKH share the remaining percentage. Experience through practical fieldwork conducted by the researcher in the study area revealed that H and/or other combined types with H as one of the aggregates and with resistivity less than or equal to 100 Ωm; such VES sections have more than 70% of striking groundwater.
The uppermost layer of almost all VES measurements conducted has resistivity above 100 Ωm, suggesting the presence of the duricrust, especially laterite with thickness ranging from ~0.4 m to ~1.5 m. The second layer is relatively thicker than the first layer (~8.4 m thick), with resistivity of about 20 Ωm corresponding to loam soil intercalated with fine sand. It is within this formation where shallow wells are taped. The third layer (~24 m, thick) with the resistivity of 26 Ωm reveals more moisturized clays. The fourth layer (~23 m) is characterized by resistivity of about 85 Ωm, suggesting the presence of fresh groundwater formation hosted in gravels, sand and alluvial sediment formation. The fifth through seven layers are characterized by resistivity above 1085 Ωm, suggesting massive felsic rock (Table 4). The VES measurements undertaken in the study area were dominated by several curve types, of which 40% belong to KH type, 19% fall under HA type, 15% were categorized as H type and 6% are of QH type. AA, KHKH, QHKH, HAA, AKH, QHKA, HKHA, KAA, KHA, QKH, QH and HAKH share the remaining percentage. Experience through practical fieldwork conducted by the researcher in the study area revealed that H and/or other combined types with H as one of the aggregates and with resistivity less than or equal to 100 Ωm; such VES sections have more than 70% of striking groundwater.
The uppermost layer of almost all VES measurements conducted has resistivity above 100 Ωm, suggesting the presence of the duricrust, especially laterite with thickness ranging from~0.4 m to~1.5 m. The second layer is relatively thicker than the first layer (~8.4 m thick), with resistivity of about 20 Ωm corresponding to loam soil intercalated with fine sand. It is within this formation where shallow wells are taped. The third layer (~24 m, thick) with the resistivity of 26 Ωm reveals more moisturized clays. The fourth layer (~23 m) is characterized by resistivity of about 85 Ωm, suggesting the presence of fresh groundwater formation hosted in gravels, sand and alluvial sediment formation. The fifth through seven layers are characterized by resistivity above 1085 Ωm, suggesting massive felsic rock (Table 4).

Two Dimensional (2-D) Geophysical Survey Results
The two-dimensional results tested on several points of the study area are pictorially presented in Figure 9. Examinations showed that the sections possess at least three layers, starting with the upper thin layer with resistivity between 250 and 300 Ωm. Its thickness ranges between 1 and 10 m with discrete to continuous extension in different locations. It represents the lateritic layer of the weathered BIF, as experienced during the field visit. The second layer is covered by the aquifer layer in most of the VES measurements undertaken, having low resistivity ranging from 20.0 to 100 Ωm with the greatest thickness in the western part of the study area ( Figure 9C) and the least ( Figure 9A) to almost no thickness in the eastern part ( Figure 9B). It is marked by the blue water drop symbol in Figure 9. Currently, operating boreholes tap their water in this layer; for example, the Bulangwa and Katente boreholes are owned by the Ushirombo town water supply, Isemabuna is used for the Butinzya water scheme, and Kazibila and Nuru boreholes are for the Runzewe town water supply scheme. The third and fourth layers cover the slightly hard to massive granite with resistivity greater than 250 Ωm, marked by a gray water drop symbol. Again, the eastern part comprises profiles ( Figure 9A,B,D,E) with thicker massive rock than the western part ( Figure 9C), concurring with the geological map, which shows the existence of the granite formation in those particular areas. Structural controlling groundwater occurrences were marked by thick downward-facing arrows ( Figure 9A-F), suggesting points prone to groundwater development.

Discussion
Demarcation of groundwater potential areas of the Bukombe semi-arid complex basement aquifer has been undertaken using both geospatial and geophysical methods with the aid of remote sensing and GIS. An overall perspective is that rainfall, geology, lineaments and land use/land cover pattern play important roles in influencing groundwater potentiality (Appendix A). Past studies [22,24,45] concluded more or less similar findings, with slight differences. Allafta pointed out lithology, geomorphology, land use/land cover, rainfall and lineament density as the main controlling factors for groundwater recharge and/or groundwater potentiality; Murasingh suggested that slope, lineament density and drainage density are the major factors influencing groundwater potential, with land use/land cover having medium influence and soil and geomorphology having the least; Rahmat et al. posited that rainfall, lithology and lineament density be the most influential factors for groundwater potentiality. Such conclusions are in line with the results of this study in most cases. In this study, however, it was found that rainfall, geology, lineament and land use/land cover, in order, are the most important factors. Deviations from previous findings are due to a number of reasons, e.g., the number and type of thematic layers used, the difference in aridity index of the study area and the type of geology/lithology of the study area, to mention few.
Locations and areal size of the very good groundwater potential area in the study are relatively bigger than those of the areas identified as poor by 2.5%. The good groundwater potential areas are principally characterized by high rainfall, gentle to flat slopes, a good distribution of lineaments leading to high lineament density, forested land cover and geological formations with good hydraulic conductivity due to the presence of loose and coarsegrained material which allows infiltration and recharge of the aquifer [7,24,38,41,[43][44][45]47]. Meanwhile, the poor groundwater potential area is characterized by high drainage density, low rainfall, few to no lineaments and steep slopes, as reported by other researchers [45].
Geology and rainfall as portrayed by AHP decision support played important roles in divulging the groundwater potential areas. Lineament, in particular, was the third most important theme in this study, concurring with note put forth by recent researchers on providing important information on subsurface fractures which control the movement and storage of groundwater [48]. Land use/land cover is the second least important, and drainage density and slope follow with about 4% weight each (Appendix A). Bukombe district aquifers with poor groundwater potential areas coincide with massive granitic formation, low rainfall, few lineaments and high-drainage-density areas, which accelerate high runoff and low infiltration rates, while the reverse leads to very good groundwater potential areas, as also reported by recent researchers [24].
In linking the geospatial findings to the geophysical survey carried out, it can clearly be observed that groundwater potential is based on structures and alluvial deposits, as prescribed earlier in this paper (Figures 3 and 8). The two-dimensional geophysical results portrayed in Figure 9 and the one-dimensional inversion charts illustrated in Appendix D agree with this statement. Structural controls related to groundwater occurrence are marked by thick downward-facing arrows ( Figure 9A-F) suggesting points prone to groundwater development. Arguably, the 1-D results given in Appendix D display evidence of alluvial deposit resistivity of fresh water ranging from 10 to 100 Ωm. Other VES stations observed to have low resistivity of clay or saline nature (>10 Ωm) in some layers, especially the second, third and fourth ones, with the third layer leading.

Results Validation
Validation is a crucial part of any modeled data in order to render the simulated results realistic. Past studies of a similar kind used different approaches to validate their final results. Among those which utilized borehole yields, the ones with low yields were mostly found in the poor groundwater potential zones, whereas those with high yields were located in the very good groundwater potential zones [18,40,43,45]. Others made use of the number of dug wells or boreholes classified into perennial and non-perennial; perennial dug wells were located in very good and good groundwater potential zones, while non-perennial dug wells were found in the poor and very poor groundwater potential zones [12,21,38]. Moreover, the 1-D geophysical survey results were also used to crosscheck the delineated groundwater potential zones [22]. This approach is supported by the results of this study but with a bantam modification. In this study, however, both 1-D and 2-D results were used to analyze the GWPAs with the aid of a drilled borehole. Thus, good GWPAs have VES curves, with a slightly wider H-curve type signifying thick aquifer formation, with medium to relatively deeper depth; the Runzewe Mcohm profile 01 ( Figure 9) and VES BKE63VES04 in Figure 8 and Appendix D are some representatives. On the contrary, the moderate to poor GWPAs are matched with VES curves characterized by shallow aquifer formation, deep and narrow H-curve type and massive formations underlying the water bearing formation, as depicted by Shikalibuga and Bulangwa Mcohm profiles in Figure 7 and VESs BKE45VES01 and BKE59VES04 in Appendix D. Additionally, the borehole drilled in Nyampangwe on 24 October 2022 helped validate the delineated GWPAs. Geophysical survey at that particular survey station indicated the presence of a thick aquifer at a shallow depth (~15 m), characterized by the resistivity of~30 Ωm. Lithological logs captured during drilling period depicted the aquifer formation as partly weathered with conglomeratic deposits at depths of 20 m to 50 m, as illustrated in Figure 10. The borehole depth was 130 m and it is currently used by villagers for domestic and mining purposes.

Conclusions
The present study demonstrated that remote sensing, GIS and AHP approaches are feasible tools for demarcating GWPAs in the Bukombe district, thereby providing a solid preliminary assessment for the groundwater resources in this basin, saving much money and time. Remote sensing data and conventional data were applied to compose the thematic layers that were then allocated appropriate weightage through the AHP technique. Based on the GWPA map generated, the study area was categorized into four different areas, viz., poor groundwater potential area (16.74 sq. km), moderate groundwater potential area (4104.07 sq. km), good groundwater potential area (3649.92 sq. km) and very good groundwater potential area (215.56 sq. km). Poor GWPAs that cover the eastern parts of Generally, the study area is envisaged with groundwater potential in most parts, varying in depth and yield. Field experience confirms that VES curves with H-type, KHtype and HKH-type having resistivity values ranging from 10 to 100 Ωm have more than 70% of striking groundwater. This finding is supported by other research [22]. Therefore, the delineated GWPAs (Figure 7b) are essential for locating other borehole sites within the study area.

Conclusions
The present study demonstrated that remote sensing, GIS and AHP approaches are feasible tools for demarcating GWPAs in the Bukombe district, thereby providing a solid preliminary assessment for the groundwater resources in this basin, saving much money and time. Remote sensing data and conventional data were applied to compose the thematic layers that were then allocated appropriate weightage through the AHP technique. Based on the GWPA map generated, the study area was categorized into four different areas, viz., poor groundwater potential area (16.74 sq. km), moderate groundwater potential area (4104.07 sq. km), good groundwater potential area (3649.92 sq. km) and very good groundwater potential area (215.56 sq. km). Poor GWPAs that cover the eastern parts of the study area can be ascribed to the accumulated influence of poor hydrogeological-environmental parameters, including low rainfall compared to the western and southern parts and the existence of hard rock regions, particularly the massive granitic rock and low lineament densities. The very good GWPA together with the good GWPA demarcated in the south, southwest and western parts of the study area are contributed to by the favorable geology (unconsolidated sediments), the low slope and drainage density, high lineament density and rainfall, and the availability of dense forest (Kigosi National Park) that favors more infiltration. The generated map is suitable for use in groundwater resources management in the Bukombe district and recommends further sites for exploration through drilling, as shown in Table 5 and Figure 11.  Acknowledgments This research was carried out as part of the first author's M. Sc. thesis at Sokoine University of Agriculture, Morogoro, Tanzania. The author would like to give his sincere thanks to the employer who gave permission to pursue the master program. He is also indebted to esteemed supervisors who provided suggestions and guidelines to shape the research. Lastly, the author expresses his profound sense of gratitude to the anonymous reviewers and editors for their helpful comments on the past drafts of the manuscript.

Conflicts of Interest:
The authors declare that there is no conflict of interest.