<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">IJG</journal-id><journal-title-group><journal-title>International Journal of Geosciences</journal-title></journal-title-group><issn pub-type="epub">2156-8359</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ijg.2023.1410053</article-id><article-id pub-id-type="publisher-id">IJG-128725</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Earth&amp;Environmental Sciences</subject></subj-group></article-categories><title-group><article-title>
 
 
  Assessment of Groundwater Potential and Prediction of the Potential Trend up to 2042 Using GIS-Based Model and Remote Sensing Techniques for Kiambu County
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Mark</surname><given-names>Boitt</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Patricia</surname><given-names>Khayasi</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Catherine</surname><given-names>Wambua</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Department of Research and Development, GIS and Remote Sensing Section-Mapinfotek, Nairobi, Kenya</addr-line></aff><aff id="aff1"><addr-line>Institute of Geomatics, GIS and Remote Sensing, Dedan Kimathi University of Technology, Nyeri, Kenya</addr-line></aff><pub-date pub-type="epub"><day>12</day><month>10</month><year>2023</year></pub-date><volume>14</volume><issue>10</issue><fpage>1036</fpage><lpage>1063</lpage><history><date date-type="received"><day>8,</day>	<month>June</month>	<year>2023</year></date><date date-type="rev-recd"><day>28,</day>	<month>October</month>	<year>2023</year>	</date><date date-type="accepted"><day>31,</day>	<month>October</month>	<year>2023</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Groundwater is one of the important necessary renewable resources of the world. It forms part of the natural water cycle that is present in the underground strata with the principal sources being precipitation and streamflow. Traditionally, information on the potential occurrence of groundwater was 
  obtained using techniques such as drilling, geophysical, geological, hyd
  ro-
  geological and geo-electrical which are time-consuming, costly and lacked full coverage. This study shows that remote sensing and GIS techniques can be utilized to map groundwater potential using a GIS-based model, the Modified DRASTIC Model, which incorporates factors that influence groundwater occurrence. These factors are the surface attributes that infer groundwater potentials and they include geology, soil texture, land use, lithology, landforms, slope steepness, lineaments and drainage systems. A prediction of the groundwater prediction was done by utilizing the MOLUSCE tool, a plugin in Qgis that utilizes ANN, multicriteria evaluation, weights of evidence and LRs algorithms in predicting land changes. The kappa value for prediction was 0.83. The results showed areas in the Southwest region had low to very low potential and the central region had high to very high potential for all the years and there were little changes between the years. The prediction showed that by 2042, the eastern region of Kiambu County will have a decline in groundwater potential.
 
</p></abstract><kwd-group><kwd>Groundwater</kwd><kwd> DRASTIC</kwd><kwd> MOLUSCE</kwd><kwd> Remote Sensing</kwd><kwd> ANN</kwd><kwd> LR</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Groundwater is one of the important necessary renewable resources of the world. It forms part of the natural water cycle that is present in the underground strata with the principal sources being precipitation and streamflow. Groundwater becomes usable when the water-bearing formations are permeable enough to allow water to infiltrate through them, to yield an adequate amount of water for use through boreholes, hand dug-well and springs, which can be replenished from recharge sources to permit continued exploitation [<xref ref-type="bibr" rid="scirp.128725-ref1">1</xref>] .</p><p>Groundwater is a valuable resource that serves as a significant source of water for communities, agricultural and industrial purposes. It plays a critical role in ecosystem sustainability and human resiliency in the face of catastrophic and unpredictable global climate change [<xref ref-type="bibr" rid="scirp.128725-ref2">2</xref>] . Compared to other sources of water, groundwater is less vulnerable to climate fluctuations in an undisturbed aquifer system and can therefore, act as a critical buffer against drought and variations in rainfall [<xref ref-type="bibr" rid="scirp.128725-ref3">3</xref>] . This brings out the need to identify groundwater potential zones for groundwater development and effective water resource management. The potential of groundwater in a region depends on different facts and varies from place to place. Since groundwater cannot be seen directly from the earth’s surface, its occurrence, distribution and movement depend upon the geological and hydro-geomorphological features of the area [<xref ref-type="bibr" rid="scirp.128725-ref4">4</xref>] .</p><p>A variety of techniques have been put forward to give information on the potential occurrence of groundwater for easy exploration of the resource. Some of the techniques comprise: drilling, geophysical, geological, hydro-geological and geo-electrical but are deemed expensive and time-consuming. For that reason, the exploration calls for the putting into practice of actual approaches of precision and that saves both time and money [<xref ref-type="bibr" rid="scirp.128725-ref5">5</xref>] . Such approaches include Remote Sensing (RS) and Geographical Information System (GIS) technologies that have been applied extensively in ground studies. RS and GIS techniques provide access to large coverage including inaccessible areas. GIS offers spatial data management and analysis tools that assist in organizing, storing, editing, analyzing and displaying positional and attribute information about geographical data [<xref ref-type="bibr" rid="scirp.128725-ref3">3</xref>] . In the recent past, several researchers have utilized these techniques in the identification of groundwater potential zones around the world. [<xref ref-type="bibr" rid="scirp.128725-ref6">6</xref>] explored groundwater potentials by integrating GIS and remote sensing technologies. [<xref ref-type="bibr" rid="scirp.128725-ref2">2</xref>] stacked a few algorithms to create a hybrid model and used the model to perform groundwater mapping. On the other hand, [<xref ref-type="bibr" rid="scirp.128725-ref5">5</xref>] integrated remote sensing and vertical electrical sounding in mapping groundwater potential in Asals areas. [<xref ref-type="bibr" rid="scirp.128725-ref7">7</xref>] and [<xref ref-type="bibr" rid="scirp.128725-ref8">8</xref>] were keen to employ a resistivity survey technique in combination with GIS and remote sensing in prospecting groundwater potential zones.</p><p>In this study, an endeavour was made towards delineating and classifying underground water potential areas by using an overlay GIS model, the modified DRASTIC model which is a substitute of the initial D—depth to groundwater, R—recharge rate, A—aquifer, S—soil, T—topography, I—vadose zone’s impact, and C—aquifer’s hydraulic conductivity (DRASTIC) model that was used to map groundwater vulnerability zones [<xref ref-type="bibr" rid="scirp.128725-ref9">9</xref>] . This model is one of the better-known and extensively used in various studies to assess the intrinsic vulnerability of groundwater by considering known properties of the aquifers [<xref ref-type="bibr" rid="scirp.128725-ref10">10</xref>] . For this research, the factors used in the initial model are replaced with the factors that influence groundwater occurrence. These factors are the surface attributes that infer to groundwater potentials and they include: geology, soil texture, land use, lithology, landforms, slope steepness, lineaments and drainage systems. In addition, the climate conditions especially rainfalls also control movement and storage of groundwater.</p><p>Kiambu County is one of the 47 counties in Kenya. The county experiences high population projections due to the influx of people from other counties in search of better livelihood. Kiambu is situated next to the capital city, Nairobi and it is perceived that people who work in the city prefer to live in Kiambu County and surrounding areas where there are less congestion and well-developed infrastructure [<xref ref-type="bibr" rid="scirp.128725-ref11">11</xref>] . Furthermore, there has been a general shift in land use from agriculture to residential and commercial buildings in the area due to the increase in population. For that reason, there is need to diversify water sources to meet the high demand for water in Kiambu County. Some of these sources include groundwater resource, which is highly ignored in Kiambu and Kenya in general [<xref ref-type="bibr" rid="scirp.128725-ref3">3</xref>] . Therefore, this resource can be developed to supplement the water from other sources which is not enough for the growing populations in Kiambu county. Additionally, mapping groundwater potential zones will provide information on the water potential of Kiambu County and prevent the effect of sudden drying of boreholes or minimal yields from the boreholes during the dry seasons.</p><p>This paper, therefore, aims to provide information about groundwater potential zones and predict the future of groundwater potential in Kiambu County using the modified DRASTIC model together with Remote sensing techniques for further groundwater exploration, proper planning, sustainable utilization and management of groundwater resources.</p></sec><sec id="s2"><title>2. Description of Study Area</title><p>Kiambu County lies between latitudes 10˚20' and 00˚25'S and longitudes 36˚31' and 37˚15'E. As per the 2019 census report, the county’s population is 2,417,735 with a population density of 952.4/km<sup>2</sup>. The county covers an area of approximately 2538.6 km<sup>2</sup> (<xref ref-type="fig" rid="fig1">Figure 1</xref>).</p></sec><sec id="s3"><title>3. Materials and Method</title><sec id="s3_1"><title>3.1. Data</title><p>To identify the groundwater potential zones in Kiambu County, different data were used to prepare various thematic layers of the study area. The used data incorporated satellite imagery in conjunction with auxiliary data such as geology, soil, lithology, geomorphology and rainfall as shown in <xref ref-type="table" rid="table1">Table 1</xref> below.</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Data sources</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Data</th><th align="center" valign="middle" >Source</th><th align="center" valign="middle" >Spatial Resolution (m)</th><th align="center" valign="middle" >Purpose</th></tr></thead><tr><td align="center" valign="middle" >Landsat 4 (TM)/Landsat 7 (ETM+)/ Landsat 8 (OLI) (1992, 2002, 2012, 2022)</td><td align="center" valign="middle" >Earth Explorer</td><td align="center" valign="middle" >30 m</td><td align="center" valign="middle" >To extract lineament To perform LULC classifi-cation</td></tr><tr><td align="center" valign="middle" >SRTM DEM</td><td align="center" valign="middle" >Earth Explorer</td><td align="center" valign="middle" >30 m</td><td align="center" valign="middle" >Obtaining slope and drainage density</td></tr><tr><td align="center" valign="middle" >Soil</td><td align="center" valign="middle" >ISRIC</td><td align="center" valign="middle" ></td><td align="center" valign="middle" >Obtaining soil types and texture</td></tr><tr><td align="center" valign="middle" >Lithology &amp; Geomorphology</td><td align="center" valign="middle" >RCMRD</td><td align="center" valign="middle" ></td><td align="center" valign="middle" >To produce lithology and geomorphology maps</td></tr><tr><td align="center" valign="middle" >Geology</td><td align="center" valign="middle" >USGS</td><td align="center" valign="middle" ></td><td align="center" valign="middle" >To produce a geology map</td></tr><tr><td align="center" valign="middle" >Rainfall-CHIRPS (1992, 2002, 2012 &amp; 2022)</td><td align="center" valign="middle" >USGS FEWS NET DATA Portal</td><td align="center" valign="middle" >5 Km</td><td align="center" valign="middle" >To give rainfall estimates</td></tr></tbody></table></table-wrap></sec><sec id="s3_2"><title>3.2. Methodology</title><p>The potentiality of groundwater in Kiambu County was evaluated using the modified DRASTIC approach. The approach considers the intrinsic characteristics of an aquifer, in hand with anthropogenic activities on the earth surface [<xref ref-type="bibr" rid="scirp.128725-ref10">10</xref>] . The factors that were employed in the model included: soil, land use land cover, lineaments density, drainage density, geology, geomorphology, lithology, slope and rainfall. The processing of these factors is given in detail in the subsections below.</p><sec id="s3_2_1"><title>3.2.1. Land Use Land Cover Classification</title><p>The three Landsat Imageries, Landsat 4 Thematic Mapper (TM)/Land-sat 7 Enhanced Thematic Mapper Plus (ETM+)/ Landsat 8 Operational Land Imager (OLI) acquired for the years 1992, 2002, 2012 &amp; 2022 were pre-processed and later classified to obtain the land use land cover of the study area. The images were classified into six classes namely: Built-up, Cropland, Grassland, Bareland, Water and Forest land.</p></sec><sec id="s3_2_2"><title>3.2.2. Lineaments Extraction</title><p>The Landsat images were used to extract lineaments. Principal Component Analysis (PCA) an image enhancement technique was performed where principal component statistics were computed in ENVI software. PCA was performed in order to minimize the dimensions of the data and to collect the most information in a band. The PCA image is used as an input image in PCI Geomatica for automatic lineament generation. The algorithm for lineament extraction in PCI Geomatica software consists of edge detection, thresholding and linear extraction steps which were carried out.</p></sec><sec id="s3_2_3"><title>3.2.3. Soil</title><p>The KENSOTER database which contained the soil components was processed. This was done in ArcGIS environment. The downloaded soil texture data was processed to obtain the soil texture map.</p></sec><sec id="s3_2_4"><title>3.2.4. Geology</title><p>A geology map was obtained from the downloaded Africa’s geology data. The data was clipped to the area of interest and projected to be able to identify the various rocks in the study area. It was converted to raster for the model input.</p></sec><sec id="s3_2_5"><title>3.2.5. Geomorphology &amp; Lithology</title><p>Data downloaded on geomorphology and lithology was processed to identify the various geomorphic and lithological units in Kiambu County.</p></sec><sec id="s3_2_6"><title>3.2.6. Rainfall</title><p>Climate Hazards Group InfraRed Precipitation with Station Data (CHIRPS) is a quasi-global rainfall dataset that is obtained by combining data from real-time observing meteorological stations with infra-red data to estimate precipitation [<xref ref-type="bibr" rid="scirp.128725-ref12">12</xref>] . To obtain the annual rainfall estimates for each of the study years, the data was processed. The processes consisted of projections, resizing and spatial analysis. Spatial analysis entailed the computation of cell statistics and Inverse Distance Weighting (IDW) interpolation.</p></sec><sec id="s3_2_7"><title>3.2.7. Slope and Drainage</title><p>The Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM) was processed in the ArcGIS Spatial Analyst environment. The slope tool was applied in acquiring the slope map of the area. The hydrology tools were used to obtain the watershed of the area. The drainage properties of the watershed were analyzed to determine the drainage density.</p><p>These contributing factors were weighted according to the weights in <xref ref-type="table" rid="table2">Table 2</xref> to obtain the final groundwater potential map</p><p>The above discussed processes are captured <xref ref-type="fig" rid="fig2">Figure 2</xref> below. To project the groundwater potential to 2042, Artificial Neural Network (ANN) technique inside the Modules for Land Use Change Evaluation (MOLUSCE) plugin was used to model transition potentials and simulate the future. ANN approach is more effective than Linear Regression (LR) [<xref ref-type="bibr" rid="scirp.128725-ref13">13</xref>] . MOLUSE, a plugin that is employed in Quantum GIS (Qgis) environment, is designed to analyze, model and simulate land use changes. The plugin incorporates well-known algorithms, which can be used in land change analysis, urban analysis as well as forestry applications and projects.</p><p>Based on the groundwater potential data for 2002, 2022, explanatory variables and transition matrices, we projected the groundwater potential for 2042. To validate the model and prediction accuracy, the plugin offered a kappa validation technique and comparison of actual and projected images, and a kappa value of 0.83 was obtained. In the ANN learning process, 100 iterations were chosen for the projection. <xref ref-type="fig" rid="fig3">Figure 3</xref> outlines the steps used in implementing the prediction in the MOLUSCE model.</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Weights applied in the DRASTIC based overlay scheme</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Factors</th><th align="center" valign="middle" >Weights (%)</th></tr></thead><tr><td align="center" valign="middle" >LULC</td><td align="center" valign="middle" >12</td></tr><tr><td align="center" valign="middle" >Slope</td><td align="center" valign="middle" >12</td></tr><tr><td align="center" valign="middle" >Lineament density</td><td align="center" valign="middle" >15</td></tr><tr><td align="center" valign="middle" >Drainage density</td><td align="center" valign="middle" >15</td></tr><tr><td align="center" valign="middle" >Geology</td><td align="center" valign="middle" >10</td></tr><tr><td align="center" valign="middle" >Rainfall</td><td align="center" valign="middle" >8</td></tr><tr><td align="center" valign="middle" >Geomorphology</td><td align="center" valign="middle" >10</td></tr><tr><td align="center" valign="middle" >Lithology</td><td align="center" valign="middle" >10</td></tr><tr><td align="center" valign="middle" >Soil types</td><td align="center" valign="middle" >4</td></tr><tr><td align="center" valign="middle" >Soil texture</td><td align="center" valign="middle" >4</td></tr><tr><td align="center" valign="middle" >Total</td><td align="center" valign="middle" >100</td></tr></tbody></table></table-wrap></sec></sec></sec><sec id="s4"><title>4. Result and Discussion</title><sec id="s4_1"><title>4.1. Land Use Land Cover</title><p>Land use land cover distribution within the area was obtained for 1992, 2002, 2012 and 2022 as shown in <xref ref-type="fig" rid="fig4">Figure 4</xref>. From the analysis, it was realized that built -up areas increased gradually from 1992 to 2022. In 1992, built up areas covered an area of approx. 164.96 km<sup>2</sup> to 347.28 km<sup>2</sup> in 2022. This effect was attributed to the increase of population in the area that led to the development of towns and growth of settlements. Contrary to built-up areas, forest areas reduced from 688.48 km<sup>2</sup> in 1992 to 370.94 km<sup>2</sup> in 2022. The reduction was assigned to deforestation to create more land for settlement and agricultural use. Croplands increased immensely from 799.65 km<sup>2</sup> in 1992 to 990.11 km<sup>2</sup> in 2022 owing to the high demand of food to feed the growing populations in the area. Water followed the same trajectory. It increased from 11.49 km<sup>2</sup> in 1992 to 65.34 km<sup>2</sup> in 2022 as a result of increased rainfall and run off water from upstream. Grasslands increased slightly from 469.71 km<sup>2</sup> in 1992 to 488.70 km<sup>2</sup> in 2022. The increase was as a consequence of bareland being converted to grassland and other land. This resulted in reduction of bareland from 410.40 km<sup>2</sup> in 1992 to 282.32 km<sup>2</sup> in 2022. The effect of land use land cover is manifested either by reduced runoff or by trapped water. Water droplets trapped go down to recharge groundwater hence vegetal cover reduces evaporation and runoff, increasing infiltration and chances of groundwater recharge as compared to barren soil. Furthermore, the root system makes the soil pervious. In settlements and built-up areas, infiltration is low due to roads, pavements and buildings covering the soil surface hence low groundwater potentials. Due to the presence of built-up structures in these areas, they are largely avoided in locating suitable sites for groundwater potential.</p></sec><sec id="s4_2"><title>4.2. Lineaments Density</title><p>Lineaments are structurally controlled linear and curvilinear features identified from the satellite imagery by their relatively linear alignments. They express the surface topography showing the zones of faulting and fracturing that increase the porosity and permeability. Lineaments were obtained for the years 1992, 2002, 2012 and 2022, and from their analysis, it was revealed that there was variability in lineament structures for each of these years. [<xref ref-type="bibr" rid="scirp.128725-ref14">14</xref>] found out that the number and orientation of the lineaments can change significantly and gradually return to its initial state after sometime. Lineaments densities were obtained for each of the years as shown in <xref ref-type="fig" rid="fig5">Figure 5</xref> below. Generally, the lineaments density was very high on the western region which was an indication of many faults and sharp change in linear alignment. Water drains through the faults to the permeable rock hence they are suitable sites of groundwater potential. The eastern part of the county comprised of a mixture of high, moderate to low lineaments density in 1992, 2002 and 2022. In 2012, the central area consisted of low lineaments density thus hinders the percolation of water into the earth surface. The information obtained from the lineaments show the movement and storage of groundwater. Potential sites for productive water are usually located around these features because, they are responsible for infiltration of surface runoff into subsurface and also for movement and storage of groundwater. Therefore, areas with high lineament density are good for ground water potential zones.</p></sec><sec id="s4_3"><title>4.3. Rainfall</title><p>CHIRPS data was processed to give the annual rainfall values for 1992, 2002, 2012 and 2022 in Kiambu county as shown in <xref ref-type="fig" rid="fig6">Figure 6</xref>. The annual rainfall ranged from 596 mm to 1488 mm in 1992, 775 mm to 1899 mm in 2002, 860 mm to 2266 mm in 2012 and 428 mm to 1396 in 2022. The highland region of</p><p>the county receives the highest rainfall through the study period. The areas towards the south, south east and south west, received the lowest annual rainfall throughout the study period. The central regions receive moderate to low rainfall as per the annual average scale. The regions towards the east and west of the county are characterized by moderate rainfall. Rainfalls are the primary sources of groundwater and dominantly recharge the groundwater. Rainfall determines the amount of water that would be available to percolate into the groundwater. High rainfall is favourable for high groundwater potentials hence it is assigned higher priority during weightings.</p></sec><sec id="s4_4"><title>4.4. Slope</title><p>The slope map of the area is shown in <xref ref-type="fig" rid="fig7">Figure 7</xref>. The area is relatively flat with the majority of the region having slope of less than 5%. The slope on the North East part of the area ranges between 5% - 25% while the slightly steep areas of &gt;25%, covers the least area. Slope influences groundwater infiltration and recharge, in sense that, areas with steep slopes cause more runoff, less infiltration and have low groundwater prospects compared to areas with gentle slope. Gentle slopes cause less runoff, high infiltration rate and have good ground water prospects hence slope is a proxy for groundwater potential analyses.</p></sec><sec id="s4_5"><title>4.5. Drainage</title><p>Stream classification is very important because it gives an understanding of a stream ecosystem [<xref ref-type="bibr" rid="scirp.128725-ref15">15</xref>] . Ordering the streams in the area showed that the region is largely endowed with the third order streams draining into the region as shown in <xref ref-type="fig" rid="fig8">Figure 8</xref> below. The fifth to sixth order streams were few same as the first and second order streams. The capacity of the fourth and third order streams were higher hence suitable for storage capacity in the area.</p><p>The stream network below was used to generate the drainage density as shown in <xref ref-type="fig" rid="fig9">Figure 9</xref>. Low drainage density areas were observed on the exterior periphery of the area since these areas are characterized with first order stream. High drainage density areas were seen on the North Eastern part towards the central part of the area. This was so because these areas have a high concentration of third and fourth order streams. When the drainage density is high is an indication of high runoff and consequently low infiltration rate. Low drainage density implies low runoff and high infiltration hence the areas have high groundwater potential.</p></sec><sec id="s4_6"><title>4.6. Soil</title><p>Infiltration of water is highly dependent on the type of soil and soil texture. The distribution of soil classes within the county were obtained as shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>0 below. Nitisols had the highest coverage while the Fluvisols, Gleysols and Phaeozems had the lowest coverage. Andosols have excellent internal drainage because of their high porosity. Planosols are of alluvial horizon having loamy or coarser textures. Ferralsols are deeply weathered and mostly clayey whereas Regosols are deep, well drained, medium textured and sensitive to erosion. They have low water holding capacity and highly permeable.</p><p>The soil texture was represented as given in <xref ref-type="fig" rid="fig1">Figure 1</xref>1. Very clayey soil texture had the highest coverage, followed by clayey and loamy while sandy soil texture had the least coverage. Soils with sandy texture have large particle constituents, which makes it to have high transmissivity and high infiltration values.</p><p>On that account, areas with this kind of texture have high groundwater recharge and potential as opposed to areas with clayey texture soils. Very clayey and clayey texture have soils with small particle constituents, resulting in low infiltration rates. Soils with loamy texture have a medium weight because they neither have a low infiltration</p></sec><sec id="s4_7"><title>4.7. Geomorphology</title><p>Geomorphology plays an essential role in the groundwater conditions of an area. Geomorphological features of a given area controls not only the occurrence but the surficial distribution of a surface water as well as the groundwater conditions. Geomorphological units of Kiambu county are captured in <xref ref-type="fig" rid="fig1">Figure 1</xref>2 and they include: Escarpments, footslope, hills and mountain footridges, mountains, plain and plateau. Majority of the area is covered by the plain and hills and mountains.</p></sec><sec id="s4_8"><title>4.8. Geology</title><p>Geology also plays an important role in the occurrence of groundwater in a given area. Geologically, the study area is underlined by formations of quaternary, tertiary and Precambrian rocks as shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>3. Tertiary rocks are the dominantly rock in the area while Precambrian rocks had the least coverage in the area.</p></sec><sec id="s4_9"><title>4.9. Lithology</title><p>Lithology is the rock composition and texture. The serial arrangement of different rocks and their interaction determines the total infiltration capacity of an area. Various landform and drainage characteristics that have a direct control on the occurrence and flow of groundwater were mapped and delineated. <xref ref-type="fig" rid="fig1">Figure 1</xref>4 shows the lithological units in the area and they comprise of volcanic, colluvium, extrusive volcanic and metaigneous. Volcanic rocks covered the largest part of Kiambu contrary to metaigneous rocks that covered the least area. The volcanic aquifer stored in the fractured and weathered parts of volcanic rocks are sources of groundwater [<xref ref-type="bibr" rid="scirp.128725-ref16">16</xref>] . Porous and permeability of lithology units translate to the storage and transmitting capacity which supports the groundwater occurrence and distribution in an area.</p></sec><sec id="s4_10"><title>4.10. Groundwater Potential Zonation</title><p>From <xref ref-type="fig" rid="fig1">Figure 1</xref>5, we observed that the suitable sites kept on changing form one year to another. In 1992, the central towards the eastern region had high potential of groundwater. The south west region had potential ranging from medium to very low, and a few areas had high potential. In 2002, the northern region began to have high potential as a result of increased rainfall and conversion of the forest land to croplands. On the other hand, potential on the eastern region reduced and this was attributed to the emergence of settlements and low rainfall. The potential on the south west region continued to reduce as people began to build settlements in those areas. In 2012, the northern region maintained high potential since the rainfalls continued to increase on those sides. The south west continued to have low potentials as settlements continued to grow in the areas. In 2022, there was a tremendous decrease in the groundwater potential for most of the regions. The potential on the northern region began to reduce to medium potential and the same applied to the eastern region. However, the central region had high potentials and the south west maintained the low potentials.</p></sec><sec id="s4_11"><title>4.11. Prediction of Groundwater Potential</title><p>This section is an important part of the study where groundwater potential of 2042 was predicted with the help of artificial neural model in QGIS. The predicted groundwater potential reveals that majority of the north side will have medium to very high potential unlike the west and some parts of the south side which will have low to very low potentials. The potential on the eastern region will reduce greatly as compared to the potential in 2022. Very high potential areas will be in regions around Ndumberi, Kijabe, Kagwe, Gathungu, Komothai and Gatundu.The ground water potential in Thika, Wangige, Kamae, Gatukuyu, Kimunyu and Ruiru will be medium. Juja, Ngenya and Chomo will have high to medium potential whereas areas around Limuru, Nachu, Kilima mbogo, Karuri, Uplands and Ngecha will have potential ranging from low to very low. The discussed above information is captured in <xref ref-type="fig" rid="fig1">Figure 1</xref>6 which shows how the groundwater potential zones will be in 2042.</p></sec><sec id="s4_12"><title>4.12. Change Analysis</title><p>To obtain the change for each groundwater potential class, the areas for each year were obtained as shown in <xref ref-type="table" rid="table3">Table 3</xref>, from which the difference between the years were computed, which represents magnitude of change between corresponding years.</p><p>Changes in groundwater potential are inextricably tied to geography, climate change, physical and socioeconomic factors. During the study period, we observed an uneven shift in land use due to rapid urban expansion and deforestation. Urban areas are given low weights in determining groundwater potentials hence the change in ground water potential between the years. Furthermore, the effect of climate change has been adverse resulting to reduced rainfalls i.e., low rainfall declines water level thus the change in groundwater potential. The changes over the years were derived as indicated in <xref ref-type="table" rid="table4">Table 4</xref> below.</p><p>The changes were little among the groundwater potential more so, for the low and very low ground water potential class. However, the changes are expected to become greater in the future if the factors affecting the occurrence of ground water continue with the same trend.</p></sec><sec id="s4_13"><title>4.13. Trend Analysis</title><p>Trend analysis of the ground water potential zones was done where, we saw the very low potential zones increase at a very low rate. The low potential zones kept on varying but it was expected to reduce by 2042. The same applied to medium potential zones which are expected to increase immensely by 2042. The high potential zones will reduce significantly while the very high potential zones increased consistently. This information is illustrated in <xref ref-type="fig" rid="fig1">Figure 1</xref>7 below.</p><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Areas for each ground water potential class in sq&#183;km</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  colspan="7"  >Areas of Groundwater potential zones (sq&#183;km)</th></tr></thead><tr><td align="center" valign="middle" >Years</td><td align="center" valign="middle" >Very low</td><td align="center" valign="middle" >Low</td><td align="center" valign="middle" >Medium</td><td align="center" valign="middle"  colspan="2"  >High</td><td align="center" valign="middle" >Very High</td></tr><tr><td align="center" valign="middle" >1992</td><td align="center" valign="middle" >511.519</td><td align="center" valign="middle" >502.237</td><td align="center" valign="middle" >504.837</td><td align="center" valign="middle" >515.558</td><td align="center" valign="middle"  colspan="2"  >485.647</td></tr><tr><td align="center" valign="middle" >2002</td><td align="center" valign="middle" >504.208</td><td align="center" valign="middle" >505.508</td><td align="center" valign="middle" >532.652</td><td align="center" valign="middle" >487.615</td><td align="center" valign="middle"  colspan="2"  >489.815</td></tr><tr><td align="center" valign="middle" >2012</td><td align="center" valign="middle" >505.762</td><td align="center" valign="middle" >503.552</td><td align="center" valign="middle" >514.447</td><td align="center" valign="middle" >505.781</td><td align="center" valign="middle"  colspan="2"  >490.256</td></tr><tr><td align="center" valign="middle" >2022</td><td align="center" valign="middle" >506.804</td><td align="center" valign="middle" >508.767</td><td align="center" valign="middle" >498.699</td><td align="center" valign="middle" >509.473</td><td align="center" valign="middle"  colspan="2"  >496.055</td></tr><tr><td align="center" valign="middle" >2042</td><td align="center" valign="middle" >507.910</td><td align="center" valign="middle" >427.911</td><td align="center" valign="middle" >731.893</td><td align="center" valign="middle" >318.925</td><td align="center" valign="middle"  colspan="2"  >533.159</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr></tbody></table></table-wrap><table-wrap id="table4" ><label><xref ref-type="table" rid="table4">Table 4</xref></label><caption><title> Change detection in sq&#183;km</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  colspan="8"  >Change Areas of Groundwater potential zones (sq&#183;km)</th></tr></thead><tr><td align="center" valign="middle" >Years</td><td align="center" valign="middle"  colspan="2"  >Very low</td><td align="center" valign="middle" >Low</td><td align="center" valign="middle" >Medium</td><td align="center" valign="middle"  colspan="2"  >High</td><td align="center" valign="middle" >Very High</td></tr><tr><td align="center" valign="middle"  colspan="2"  >2002-1992</td><td align="center" valign="middle" >−7.311</td><td align="center" valign="middle" >3.271</td><td align="center" valign="middle"  colspan="2"  >27.815</td><td align="center" valign="middle" >−27.943</td><td align="center" valign="middle" >4.168</td></tr><tr><td align="center" valign="middle"  colspan="2"  >2012-2002</td><td align="center" valign="middle" >1.554</td><td align="center" valign="middle" >−1.956</td><td align="center" valign="middle"  colspan="2"  >−15.059</td><td align="center" valign="middle" >18.166</td><td align="center" valign="middle" >0.441</td></tr><tr><td align="center" valign="middle"  colspan="2"  >2022-2012</td><td align="center" valign="middle" >1.043</td><td align="center" valign="middle" >5.215</td><td align="center" valign="middle"  colspan="2"  >−15.748</td><td align="center" valign="middle" >3.692</td><td align="center" valign="middle" >5.799</td></tr><tr><td align="center" valign="middle"  colspan="2"  >2042-2022</td><td align="center" valign="middle" >1.106</td><td align="center" valign="middle" >−80.857</td><td align="center" valign="middle"  colspan="2"  >233.194</td><td align="center" valign="middle" >−190.548</td><td align="center" valign="middle" >37.104</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr></tbody></table></table-wrap></sec><sec id="s4_14"><title>4.14. Water Quality</title><p>Not all the types of water are safe for human being. With every kind of development happening in this age, it is difficult to get naturally purified water because of industrialization that pollutes groundwater. Groundwater pollution often results from improper disposal of wastes on land. The sources of these pollutions include: industrial and household garbage, landfills, leaking underground oil storage tanks and pipelines, sewage sludge and septic systems for liquid waste [<xref ref-type="bibr" rid="scirp.128725-ref17">17</xref>] . To check the quality of water, a few parameters are investigated to determine their concentration in the water [<xref ref-type="bibr" rid="scirp.128725-ref18">18</xref>] . These parameters include manganese, fluoride, cadmium, iron, turbidity, total coliforms and Escherichia coli (E. coli) bacteria [<xref ref-type="bibr" rid="scirp.128725-ref19">19</xref>] . The contamination majorly affects the shallow well/boreholes, making the water unsuitable for use in drinking and food processing due to the presence of E. coli bacteria. [<xref ref-type="bibr" rid="scirp.128725-ref20">20</xref>] observed that shallows wells in Ruiru had faecal coliform contamination which is harmful to human health.</p></sec><sec id="s4_15"><title>4.15. Validation</title><p>Field survey was carried out to confirm the image processing results. A questionnaire was generated to help conduct the verification. The major factors that were considered in the questionnaire were: water quality, reliability and sustainability of the underground water. Data was collected from 32 boreholes/wells and the following observations were made:</p><p>1) Out of the samples collected, water was sustainable in 19 boreholes/wells since it was available throughout even during the dry seasons.</p><p>2) Water provided by 24 boreholes/wells was sufficient for the residents in the area in which the borehole/well serves.</p><p>3) The quality of water was observed to be good for 25 boreholes/wells and the water in the boreholes/wells is regularly tested to ensure it is fit for human use.</p><p>The areas that were classified to have low to very low groundwater potential on the maps, were observed to have boreholes/wells that run dry occasionally especially in times of low to no rainfalls.</p></sec></sec><sec id="s5"><title>5. Conclusion</title><p>The aim of the paper was to provide information about groundwater potential zones and predict the future of groundwater potential in Kiambu county using the modified DRASTIC model together with Remote sensing techniques for further groundwater exploration, proper planning, sustainable utilization and management of groundwater resources. The study affirms the potential of groundwater in Kiambu county mostly for the central and some parts of the eastern region. There was a consistent trend in South west region where the potential continued to decline. Additionally, we observed little changes in the groundwater potential between the years. It was noted that most the changes arose from geography, climate change, physical and socioeconomic factors. The prediction showed that by 2042, the eastern region of Kiambu county will have a decline in ground water potential and areas that initially had high potential will shift to medium.</p></sec><sec id="s6"><title>Acknowledgements</title><p>We would like to acknowledge the support received from Dedan Kimathi University of Technology and Mapinfotek Geomatiks Ltd to work on the research and development.</p></sec><sec id="s7"><title>Conflicts of Interest</title><p>The authors declare no conflicts of interest regarding the publication of this paper.</p></sec><sec id="s8"><title>Cite this paper</title><p>Boitt, M., Khayasi, P. and Wambua, C. (2023) Assessment of Groundwater Potential and Prediction of the Potential Trend up to 2042 Using GIS-Based Model and Remote Sensing Techniques for Kiambu County. International Journal of Geosciences, 14, 1036-1063. https://doi.org/10.4236/ijg.2023.1410053</p></sec></body><back><ref-list><title>References</title><ref id="scirp.128725-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Kuria, D.N. (2012) Mapping Groundwater Potential in Kitui District, Kenya Using Geospatial Technologies. The International Journal of Water Resources and Environmental Engineering, 4, 15-22. https://doi.org/10.5897/IJWREE11.119</mixed-citation></ref><ref id="scirp.128725-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Sarkar, S.K., Alshehri, F., Shahfahad, Rahman, A., Pradhan, B. and Mohamed, A. (2022) A National-Level Study on Groundwater Potentiality Mapping Using a Hybrid Machine Learning Models under the Scenario of Climate Change. https://doi.org/10.21203/rs.3.rs-1818227/v1</mixed-citation></ref><ref id="scirp.128725-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Chepchumba, M.C. (2010) Exploration of Groundwater Potential Using Gis and Remote Sensing in Embu County, Kenya. Doctoral Dissertation, JKUAT-IEET.</mixed-citation></ref><ref id="scirp.128725-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Rajaveni, S.P., Brindha, K. and Elango, L. (2017) Geological and Geomorphological Controls on Groundwater Occurrence in a Hard Rock Region. Applied Water Science, 7, 1377-1389. https://doi.org/10.1007/s13201-015-0327-6</mixed-citation></ref><ref id="scirp.128725-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Nyaberi, D.M., Basweti, E., Barongo, J.O., Ogendi, G.M. and Kariuki, P.C. (2019) Mapping of Groundwater through the Integration of Remote Sensing and Vertical Electrical Sounding in ASALs: A Case Study of Turkana South Sub-County, Kenya. Journal of Geoscience and Environment Protection, 7, 229-243. https://doi.org/10.4236/gep.2019.711017</mixed-citation></ref><ref id="scirp.128725-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Mwega, W.B., Bancy, M.M., Mulwa, J.K. and Kituu, G.M. (2013) Identification of Groundwater Potential Zones Using Remote Sensing and GIS in Lake Chala Watershed, Kenya. Proceedings of 2013 Mechanical Engineering Conference on Sustainable Research and Innovation, Kitui, 24-26 April 2013, 42-46.</mixed-citation></ref><ref id="scirp.128725-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Jhariya, D.C., Khan, R., Mondal, K.C., Kumar, T., Indhulekha, K. and Singh, V.K. (2021) Assessment of Groundwater Potential Zone Using GIS-Based Multi-Influencing Factor (MIF), Multi-Criteria Decision Analysis (MCDA) and Electrical Resistivity Survey Techniques in Raipur City, Chhattisgarh, India. Water Infrastructure, Ecosystems and Society, 70, 375-400. https://doi.org/10.2166/aqua.2021.129</mixed-citation></ref><ref id="scirp.128725-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Indhulekha, K., Chandra Mondal, K. and Jhariya, D.C. (2019) Groundwater Prospect Mapping Using Remote Sensing, GIS and Resistivity Survey Techniques in Chhokra Nala Raipur District, Chhattisgarh, India. Journal of Water Supply: Research and Technology-Aqua, 68, 595-606. https://doi.org/10.2166/aqua.2019.159</mixed-citation></ref><ref id="scirp.128725-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Musa, K.A., Akhir, J.M. and Abdullah, I. (2002) The Effect of Major Faults and Folds in Hard Rock Groundwater Potential Mapping: An Example from Langat Basin, Selangor.</mixed-citation></ref><ref id="scirp.128725-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Garewal, S.K., Vasudeo, A.D., Landge, V.S. and Ghare, A.D. (2019) Groundwater Vulnerability Mapping Using Modified DRASTIC ANP. Journal of the Croatian Association of Civil Engineers, 71, 283-296. https://doi.org/10.14256/JCE.1951.2016</mixed-citation></ref><ref id="scirp.128725-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Mugo, G.M. and Odera, P.A. (2019) Site Selection for Rainwater Harvesting Structures in Kiambu County-Kenya. The Egyptian Journal of Remote Sensing and Space Science, 22, 155-164. https://doi.org/10.1016/j.ejrs.2018.05.003</mixed-citation></ref><ref id="scirp.128725-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Funk, C., et al. (2015) The Climate Hazards Infrared Precipitation with Stations—A New Environmental Record for Monitoring Extremes. Scientific Data, 2, Article No. 150066. https://doi.org/10.1038/sdata.2015.66</mixed-citation></ref><ref id="scirp.128725-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Muhammad, R., Zhang, W., Abbas, Z., Guo, F. and Gwiazdzinski, L. (2022) Spatiotemporal Change Analysis and Prediction of Future Land Use and Land Cover Changes Using QGIS MOLUSCE Plugin and Remote Sensing Big Data: A Case Study of Linyi, China. Land, 11, Article 419. https://doi.org/10.3390/land11030419</mixed-citation></ref><ref id="scirp.128725-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Sichugova, L. and Fazilova, D. (2021) The Lineaments as One of the Precursors of Earth-Quakes: A Case Study of Tashkent Geodynamical Polygon in Uzbekistan. Geodesy and Geodynamics, 12, 399-404. https://doi.org/10.1016/j.geog.2021.08.002</mixed-citation></ref><ref id="scirp.128725-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">McManamay, R.A. and DeRolph, C.R. (2019) A Stream Classification System for the Conterminous United States. Scientific Data, 6, Article No. 190017. https://doi.org/10.1038/sdata.2019.17</mixed-citation></ref><ref id="scirp.128725-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Fenta, M.C., Anteneh, Z.L., Szanyi, J. and Walker, D. (2020) Hydrogeological Framework of the Volcanic Aquifers and Groundwater Quality in Dangila Town and the Surrounding Area, Northwest Ethiopia. Groundwater for Sustainable Development, 11, Article ID: 100408. https://doi.org/10.1016/j.gsd.2020.100408</mixed-citation></ref><ref id="scirp.128725-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Kirori, P., Matiru, V. and Mutai, J. (2022) Factors Associated with Bacterial Contamination of Shallow Well Water Sources. Case Study of Juja Hostels Kiambu County. Journal of Agriculture, Science and Technology, 21, 35-43. https://doi.org/10.4314/jagst.v21i4.4</mixed-citation></ref><ref id="scirp.128725-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Makokha, J.K. (2017) Analysis of Groundwater Quality and Identification of Abstraction Points in Kahawa Wendani, Kiambu County. Master’s Thesis, University of Nairobi, Nairobi.</mixed-citation></ref><ref id="scirp.128725-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Kilonzo, W., Home, P., Sang, J. and Kakoi, B. (2019) The Storage and Water Quality Characteristics of Rungiri Quarry Reservoir in Kiambu, Kenya, as a Potential Source of Urban Water. Hydrology, 6, Article 93. https://doi.org/10.3390/hydrology6040093</mixed-citation></ref><ref id="scirp.128725-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Otieno, R.O. (2016) Seasonal Assessment of Groundwater Quality in Ruiru Kiambu County, Kenya.</mixed-citation></ref></ref-list></back></article>