Full Terms & Conditions of access and use can be found at https://www.tandfonline.com/action/journalInformation?journalCode=tjag20 NRIAG Journal of Astronomy and Geophysics ISSN: (Print) (Online) Journal homepage: https://www.tandfonline.com/loi/tjag20 Triangulation Approach for Mapping Groundwater Suitability Zones in Coastal Areas Around Lagos, Nigeria. Using Multi-Criteria Decision-Making Technique Akintunde A Oludayo, Ozebo V Chidi & Oyedele F Kayode To cite this article: Akintunde A Oludayo, Ozebo V Chidi & Oyedele F Kayode (2021) Triangulation Approach for Mapping Groundwater Suitability Zones in Coastal Areas Around Lagos, Nigeria. Using Multi-Criteria Decision-Making Technique, NRIAG Journal of Astronomy and Geophysics, 10:1, 423-442, DOI: 10.1080/20909977.2021.1987118 To link to this article: https://doi.org/10.1080/20909977.2021.1987118 © 2021 The Author(s). Published by Informa UK Limited, trading as Taylor & Francis Group. Published online: 26 Nov 2021. Submit your article to this journal View related articles View Crossmark data https://www.tandfonline.com/action/journalInformation?journalCode=tjag20 https://www.tandfonline.com/loi/tjag20 https://www.tandfonline.com/action/showCitFormats?doi=10.1080/20909977.2021.1987118 https://doi.org/10.1080/20909977.2021.1987118 https://www.tandfonline.com/action/authorSubmission?journalCode=tjag20&show=instructions https://www.tandfonline.com/action/authorSubmission?journalCode=tjag20&show=instructions https://www.tandfonline.com/doi/mlt/10.1080/20909977.2021.1987118 https://www.tandfonline.com/doi/mlt/10.1080/20909977.2021.1987118 http://crossmark.crossref.org/dialog/?doi=10.1080/20909977.2021.1987118&domain=pdf&date_stamp=2021-11-26 http://crossmark.crossref.org/dialog/?doi=10.1080/20909977.2021.1987118&domain=pdf&date_stamp=2021-11-26 Triangulation Approach for Mapping Groundwater Suitability Zones in Coastal Areas Around Lagos, Nigeria. Using Multi-Criteria Decision-Making Technique Akintunde A Oludayo a,b, Ozebo V Chidi a and Oyedele F Kayode c aDepartment of Physics, University of Lagos, Lagos, Nigeria; bSchool of Science and Technology, Pan Atlantic University, Lagos, Nigeria; cDepartment of Geoscience, University of Lagos, Lagos, Nigeria ABSTRACT Water quality in Lagos's coastal area is prone to constant seawater intrusion due to its proximity to the Atlantic Ocean. Seawater intrusion in the coastal aquifer is a major challenge to human survival and therefore needs urgent attention due to deficiency in groundwater quality. This study applies integrated techniques to map groundwater suitability zone (GSZ) around the Lagos coastal area for sustainable development. The objectives were achieved by preparing eight different thematic layers that influence groundwater suitability viz; drainage density, geology, geomorphology, lineament density, rainfall, longitudinal conductance, slope, and transverse resistance. Weighting assignment, normalization and, pair-wise comparison of the thematic maps were done using Saaty’s Analytical Hierarchical Process (AHP) model before being integrated through the geographical information system (GIS) spatial analyst model with weight overlay to generate the overall groundwater suitability map. The integrated result reveals that 7% of the research location has very low suitability, 32.2% low, 47.8% moderate, and 13% high groundwater suitability. Thus, the groundwater suitability zone (GSZ) was validated via Electrical conductivity (EC) and total dissolved solids (TDS) and the result reveals a good correlation. The research has proven the efficacy of integrating different surface and subsurface data for effectively mapping the groundwater suitability zone (GSZ). ARTICLE HISTORY Received 1 March 2020 Revised 2 September 2021 Accepted 25 September 2021 KEYWORDS Remote sensing; electrical resistivity; geographic information system; analytical hierarchy process; Groundwater Suitability zone 1. Introduction Population increases in the coastal areas around the world, especially in the Lagos province coupled with increasing urbanisation and industrialisation have war ranted greater demands for potable water supply for socio-economic development. Also, the demand for water has put more pressure on the coastal aquifer sys tems since it is the only source of water used for human consumption and other anthropogenic events in the metropolitan city of Lagos and the surface waters present are somewhat contaminated. Also, improper waste man agement from household sewages disposal and hazar dous industrial chemical wastes have led to degradation in groundwater resources (Bachaer et al. 2018). More so, the coastal regions are prone to seawater intrusion (SWI) because the equilibrium of hydrostatic balance between the groundwater protective capacity and seawater inter face has been broken. The main causes of seawater intru sion into freshwater are high demographic density. This is has a result of over-exploitation and intensive pumping (Bagyaraj et al. 2012). Also, high relief slope, topographi cal depression, precipitation, and high drainage density can increase surface runoff that infiltrates into the coastal aquifer (Demiroğlu and Dowd 2014). However, poor economic situation, inadequate infrastructural facilities from the government like pipe-borne water, and alternative water supplies have made people living around metropolis source for alternative water supply which is solely groundwater. A balance between human demands for water consumption and groundwater sup ply has not been attained for their sustainability. Also, the quantity and quality of groundwater have been changing due to the effect of anthropogenic activities coupled with natural geological influence (Chawla et al. 2010; Akintunde 2018). Groundwater is found in the subsurface of the earth’s crust. It is void occupying water in rock cracks and pore space in sands. It serves as a reservoir for water storage and slowly moves within the geological formation of the subsurface litho-stratification. Groundwater occurrence at any place on the earth is due to consequential interac tion of climate change, geological formation, hydrological and ecological factors (Scibeka et al. 2007). Soil porosity, topography, and permeability of the underlying litholo gies are the dynamic factors controlling groundwater storativity and transmissivity. The study area lithologies are the constituent that forms the stratigraphy of various geomorphic units, porosity media, and permeability thereby causing some noticeable change in the ground water potential (Surajit and Ramakar 2015). Hydro- geological parameters such as topographic slope, CONTACT Akintunde A Oludayo dayodeji@live.com Department of Physics, University of Lagos, Ibeju-Lekki, Km 52 Lekki-Epe Expressway, Lagos, Nigeria NRIAG JOURNAL OF ASTRONOMY AND GEOPHYSICS 2021, VOL. 10, NO. 1, 423–442 https://doi.org/10.1080/20909977.2021.1987118 © 2021 The Author(s). Published by Informa UK Limited, trading as Taylor & Francis Group. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. http://orcid.org/0000-0002-8642-0481 http://orcid.org/0000-0001-5062-0548 http://orcid.org/0000-0002-5711-766X http://www.tandfonline.com https://crossmark.crossref.org/dialog/?doi=10.1080/20909977.2021.1987118&domain=pdf&date_stamp=2021-11-20 geological formation, lineament extraction, drainage net work, rainfall, geomorphological structure, and geophy sical surveys are factors that can be used in groundwater recharge investigation. Recently, geographic information systems (GIS) and remote sensing (RS) are valuable tools for spatial modelling of hydrogeological exploration of a large surface area with more dependability (Gyeltshen et al. 2020). GIS and RS are time-saving tools and eco nomical ways of demarcating groundwater monitoring management. The multi-criteria decision model (MCDM) using the Analytical Hierarchical Process (AHP) technique has become an effective tool adopted by numerous researchers for an accurate description of groundwater suitability studies (Pinto et al. 2015; Fashae et al. 2014). In recent times, researchers have used a different inte grated approach to map and monitor groundwater in some parts of Lagos province (Adeoti et al. 2010; Oyedele and Oladele 2011; Aladejana et al. 2020). In some of these studies, saltwater intrusion has been iden tified in the first aquifer unit, such as Sangotedo, Ajah, Lakowe, Lekki, and Agboyi creek. The degradation of groundwater quality around the Oniru and Eputu- Lakowe environs may be due to the high infiltration of saltwater in the freshwater aquifer zone (Ozebo and Ajiroba 2011; Oyedele and Oladele 2011). (Anthony 2016; Adiat et al. 2012) have successfully combined geo physical methods and geospatial techniques using multi- criteria decision analysis to develop groundwater poten tial and water monitoring models. It is, therefore, necessary to identify and map potential areas around the Lagos coastal area for sustainable development. The present study focuses on applying an integrated approach to identify and map quality groundwater suitability zone (GSZ) around the Lagos lagoon metropolis using geos patial and geophysical techniques. These objectives will be achieved by preparing the thematic maps of contribut ing parameters such as drainage, lineament, geology, geomorphology, rainfall, and slope through the GIS environment. Also, a geophysical investigation was car ried out to corroborate the results from the geospatial data. Then, Groundwater suitability zones (GSZ) map was generated for the study areas. Finally, the GSZ map was validated with the physico-chemical analysis of some collected groundwater samples. 2. Location, hydrology, and geology the study area This study area falls within the Lagos coastal plain which is accessible through waterways within the metropolis. According to Badejo et al. (2014), the Lagos Lagoon is very significant because it not only extends across the southern region of the metropolis but also serves as a linkage to the Atlantic Ocean especially along the western and southern as well as Lekki Lagoon in the eastern part of the lagoon. The study location lies between longitude 3o200Eto4o150E and Latitude 6o23N0to6o440N as shown in Figure 1. The proximity to the Atlantic Ocean and the variation Figure 1. Map of the study area showing borehole points, VES points and elevation (Nigeria Geological Survey Agency 2006 and Agagu 1985). 424 A. A. OLUDAYO ET AL. in both seasonal and semi-diurnals are the main sources of salinity (Medina-Gómez et al. 2014; Mohamed et al. 2014). The weather season in Lagos varies significantly throughout the year, the wetter season lasts 6.8 months, from April to October while the drier season lasts for 5.2 months, from November to March (Adenodi 2018). The study area is charac terised by a multi-aquifer system that consists of three major aquifer horizons differentiated by silty or clay layers. The first aquifer unit is the water table aquifer and is highly prone to pollution due to its closeness to the ground surface (Oyedele and Oladele 2011). The second aquifer unit is a confined aquifer that consists of discontinuous arrangements of sand and clay. These aquifer units are embedded within the continental Ilaro Formation and are harnessed through the borehole. Though some of this aquifer is slightly contaminated, and it is the most exploited (Adeoti et al. 2010; Olorode et al. 2016). The third aquifer consists of fine and medium sand and it is the most reliable aquifer used by industries and water works in Lagos metropolis. The groundwater in this aquifer exists as confined to the semi-confined aquifer (Longe et al. 1987). Dahomey represents one of the sedimentary basins that can be found at the Guinea Gulf continental boundary. The basin is a marginal pull part-system of the West African margin which erupted during the rifting era that separated the African and South American plates and subsequently create the Gulf of Guinea from the late Jurassic era to the early Cretaceous time (Burke et al. 1971; Whiteman 1982). However, Wilson and Willians (1979) suggested that Dahomey is belted to the western part of Ghana escarpment that is essentially an adjunct of the Romanche fracture belt and bounded to the eastern part of the Benin Hinge line. Furthermore, the base ment ridge that differentiates the Niger Delta Basin from the Okitipupa Structure and consequently became the Chain Fracture extension zone of the continent (Billman 1976). The stratigraphic column and the geology map of the study area are present in Table 1 and Figure 1, respectively (GAF et al. 2016 and Agagu, 1985). The major geological features of the study area are the Coastal Plain sand and Coastal Alluvium which belongs to the Quaternary to the Recent Formation with an alternating sequence of clay, clayey sand, pebbly sand, sandy clay, sand, lignite, and silt (Jones and Hockey 1964). While the minor geological units are the Ilaro formation, Ewekoro for mation, Oshosun formation, and Abeokuta complex basin which are the shale, sandstone, limestone, marls gastropod, and siltstone (Omatsola and Adegoke 1981). High explorations of groundwater in the alluvia belt are proportional to the high rate of upward move ment of saline water in this zone (Longe 2011). 3. Materials and methods The method employed in this study is centred on three basic phases; Firstly, is the use geospatial approach to create thematic layers such as lithological layers, drai nage patterns, lineament networks, hydrological cycle, and geomorphological data. Secondly, in-situ data using the geophysical electrical resistivity, digital image processing (DIP) from the remote sensing data were used for spatial classification of groundwater suitability zone (GSZ). In the third phase, a Multi- Criteria Decision Analysis (MCDA) approach based on the analytical hierarchical process (AHP) was used to integrate the multiple data set in an ArcGIS environment. 4. 3.1: Geo-spatial data acquisition Drainage Network: The Shuttle Radar Topographic Mission Digital Elevation Model (STRM DEM) 30 m resolution (1 arc-second) data was acquired from USGS Earth Explorer. The topographic data was geo- referenced to UTM coordinate Systems in Minna Datum 32 N using a geographic information system (GIS). The Arc-Hydro tool of ArcMap employed involves a start from the sink fills up to the stream to feature function for extracting stream channels from the topographic map in Figure 4. Researchers like (Prasad et al. 2008; Singh et al. 2013; Olorunfemi and Table 1. Eastern Dahomey basin stratigraphic units. Jones and Hockey (1964) Omotsola and Adegoke (!881) Agagu (1985.) GAF et al. 2016 Era Age Formation Age Formation Age Formation Lithology Quaternary Recent Alluvium Alluvium Silt Clay Sand Tertiary Pleistocene - Oligocene Eocene Palaeocene Coastal Plain sands Ilaro Ewekoro Pleistocene - Oligocene Eocene Palaeocene Coastal Plain sand. Ilaro Oshosun Akinbo Pleistocene - Oligocene Eocene Palaeocene Coastal Plain sand Ewekoro Oshosun Akinbo Ewekoro Sandstone Shale Limestone Cretaceous Late Senonian Abeokuta Maastrichtian – Neocomian Araromi Afowo Ise Maastrichtian – Neocomian Araromi Afowo Ise Conglomerates, Sandstones, Shales and Minor Carbonate Precambrian Crystalline Basement Rocks. NRIAG JOURNAL OF ASTRONOMY AND GEOPHYSICS 425 Oni 2019) have successfully used the fill tools for different spatial analysis modules. The drainage net work was used to indicate hydrological conditions. The density of the drainage network is restrained by its lithology and the precipitation infiltration rate of the subsurface runoff (Edet et al. 1998). The Drainage density is the total length of the drainage line via the sub-basin zone and was expressed in terms of km/km2. For the purpose of the research, a high drainage den sity zone is allocated with higher weight while a low drainage density zone is assigned with lower weight as presented in Table 2 (Shaban et al. 2006). The drainage network is distinctively classified with a natural break (Jenks) into five categories as presented in Table 2. Lineament extraction: Lineament analysis shows the relationship between the surface water penetration and weathering zones that control the water mobility and infiltration rate (Recep 2019). The Landsat 8 (OLI_TIRS) enhanced thematic mapper (ETM+) sen sor of path 190 and raw 055 and 056 acquired in May 2019 from USGS Earth Explorer was employed to extract the lineament pattern. The Landsat data reflectivity was enhancement via radiometric correc tion and the principal component analysis (PCA) using forward PC rotation in the ENVI 5.1 software. The automated lineament extraction was completed with the PCI Geomatica 2018 software using the LINE extraction tool from the algorithm librarian. The algo rithm of the LINE extraction modular consists of edge detection, thresholding, and curve extraction, which were the three stages of extracting the Lineaments pattern. Rockwork 17 software was used to identify the bearing of the lineament lines using the rose dia gram. The rose diagram tools are used to analyse lineament orientation and the directional frequency of the extracted lineament over the area. Layers of the lineament features were enhanced and improved using the panchromatic band of the ETM+ to increase details on the GIS model, especially layer grouping and classification. The lineament density was calcu lated by dividing the total lineament length by the cross-sectional area of the basin and its expressed in km/km2. The lineament density was classified into five categories as presented in Table 2 with natural break (Jenks) intervals. Geological structures: The geological map was geo-referenced, digitised, and classified from the 1: 50,000 toposheet 279 SW of the Nigeria Geological Survey Agency (NGSA) using ArcGIS Tools. More so, the geology of the study area was reclassified (Table 2) by their influence on saltwater intrusion. The Abeokuta complex basin is given the magnitude of 5, the Ewekoro formation is given the magnitude of 4, the Ilaro and Oshosun formation were classed as 3 while the coastal plain and the coastal Alluvia were classified as 2 and 1 respectively (Figure 6). Hydro-geomorphological map: the geomorphology features are the topographic features that illustrate the various landform and landscapes of the study environ ment. The geomorphology map was extracted from the world geological map (One-Geology database, which is affiliated with the British Geological Survey (BGS)). The extracted data was geo-referenced using the ArcGIS referencing tool, then ratified into the tiff format for layer integration before being digitised for classification. Also, the geomorphology map was enhanced by differ ent geomorphology features presented by (Allen 1965) and updated from the field survey. Weight was assigned to each geomorphological feature regarding their influ ence on saltwater intrusion (Offodile 2014). Features with higher influence were assigned the lowest value and vice versa (Table 2). Slope and elevation: To determine the groundwater controlling factors, the slope map plays a significant role (Kadam et al. 2012). The research zone slope map was Table 2. Categorisation of factors influencing of groundwater suitability zones. Factors Assigned Weight Normalised Weight (w) Domain Effect Rank Geology (Ge) 7 0.175 Abeokuta Basin 5 Ewekoro formation 4 Oshosun formation 3 Ilaro formation 3 Coastal plain sand 2 Coastal Alluvia 1 Geomorphology (Gm) 2 0.050 Coastal Alluvium 4 Beach ridge 2 Coastal creek 2 Mangrove Swamp 3 Floodplain 4 Lagoon 1 Pedi plain 5 Lineament Density (Ld.) 6 0.150 0–0.319 5 0.320–0.562 3 0.56–1.76 1 Drainage density (Dd) 5 0.125 4 5 3 4 2 3 1 2 0 1 Slope 4 0.100 0.001–0.262 5 0.262–1.442 4 1.442–2.360 3 2.360–3.541 2 3.541–33.44 1 Longitudinal conductance (Lc.) 3 0.075 0.03208–0.4157 6 0.4158–0.7994 5 0.7995–1.617 4 1.618–2.457 3 2.458–5.762 2 5.763–14.010 1 Transverse Resistance (Tr.) 8 0.200 7.006–25.00 1 25.01–100.00 2 100.01–250.00 3 250.01–500.00 4 500.01–1000.00 5 1000.01–2000.00 6 Rainfall (Rf) 5 0.100 178.888 1 135.039 2 219.779 3 Total Sum 40 1 37 426 A. A. OLUDAYO ET AL. generated from the STRM DEM topographic map 30 m resolution (1 arc-second) data downloaded from USGS Earth Explorer and was processed with the ArcGis spatial tools. The slope was classified based on the SOTER model in the ArcGis environment and then reclassified into a close level, very gentle, somewhat gentle, somewhat steep, and highly steep as presented in Figure 6 (European Commission 1995). The lowest classes slope is ranked as five and one is given the highest step of the slope as presented in Table 2. Rainfall: Rainfall is a vital factor that affects the hydrological cycle and represents a direct recharge source in groundwater. Rainfall data can be used to determine the quantity of water present in ground water. Rainfall data were collected from the Nigerian Meteorological Agency (NIMET) stations in Lagos state. The annual rainfall data used covers from 2003 to 2017 and the average annual rainfall data for 15 years is plotted in Figure 2. Rainfall distribution is affected by the slope gradient because it increases the infiltration rate of water runoff and groundwater rechargeability (Gyeltshen et al. 2020). The rainfall data were processed in the ArcGis platform using the Thiessen polygons tools to get the rainfall spatial dis tribution map. Researchers have used Thiessen poly gon to analyse climate data when there are fewer weather stations in the state. The Thiessen polygon contains a single influential point input feature so that any location within the study area is closer to its associated point than to any other input point (Heywood et al. 1998). 4.1. Geophysical field survey A geo-electrical resistivity method was employed to determine the physical parameters of the sub-surface lithology based on its apparent resistivity. A total of one hundred and ninety-seven (197) Vertical Electrical Sounding (VES) were conducted at Twenty- five locations around the study area using the Schlumberger arrangement with maximum current electrode spacing (AB/2) of 200 m at each point (Figure 1). A period of four months (May to August 2019) was used in the acquisition of the geo physical data. The equipment used in acquiring data was Pasi Terameter (16-GL). The accessories were electrodes, hammers, tape measures, and Global Positioning System (GPS) for geo-referencing the study points. The measured apparent resistivity values were interpreted quantitatively and qualitatively using the partial curve matching technique and computer iterations (Win-Resist software) (Vander Velpen 2004). Out of the 197, VES stations sounded, 100 VES points were presented in Table 3 showing the longitude, latitude, layer resistivities, and layer thick ness of each station. Four to six subsurface geo- electrical layers were delineated from the interpreted VES. The distribution of their curve types is presented in Figure 3. Due to the ambiguity in the interpretation of resistivity, and the similarities in geophysical sub surface lithology, the Dar Zarrouk parameter was employed. Geoelectrical parameter: The importance of using the Dar Zarrouk parameter is to solve the problem of overlapping character and reduce the uncertainty in the interpretation of the VES result (Henriet 1976; Batayneh 2013). The Dar Zarrouk parameters were used to determine the resistivity regime of the water- bearing formations or water bodies while dealing with complicated situations of resistivity ranges intermix ing for saline, brackish, and freshwater aquifers. Dar Zarrouk parameters were used to identify the long itudinal conductance and transverse resistance from the VES data. In demarcating groundwater suitability, Figure 2. Average monthly rainfall data for the past 15 years (2005–2019) at the three stations. NRIAG JOURNAL OF ASTRONOMY AND GEOPHYSICS 427 Table 3. 100 vertical electric sounding results and the estimated Dar Zarrouk parameters. NAME? LONG LAT R1 R2 R3 R4 R5 R6 H1 H2 H3 H4 H5 TR LC Owode-Ibese 3.473067 6.552704 178.9 636.6 1586.7 4404.3 588.4 0.7 6.8 22.8 75 74,190.67 0.045993 3.475682 6.555943 215.4 510.2 1218.9 3155.6 662.1 0.6 3.7 10.4 23.2 17,580.69 0.025922 3.480382 6.55546 309.4 506.1 1045.9 6132.4 394.6 0.7 4.6 15.1 83 105,465.4 0.039324 Ofin Oreta 3.522196 6.537954 343.2 609.9 1872.8 6081.6 467.8 0.7 5.1 16 57.9 77,088.03 0.028466 3.51591 6.597493 15.3 2.5 0.7 5 17.2 0.6 2.6 1.1 19.7 22.99 6.590644 3.508483 6.553166 4.8 1.2 0.5 2.3 12.2 0.6 2.8 3.7 35.4 17.902 25.24964 3.500127 6.541559 1.3 1.7 2.7 8.8 0.6 6.2 16.4 11.12 10.18267 Bayiku 3.554542 6.546756 186.3 711.6 2443.7 718.5 0.7 5.7 56.6 28,499.99 0.034929 3.566826 6.567844 337.2 1291.9 3802.5 1290.1 528.1 0.8 4 22.2 100.1 43,798.37 0.088898 3.53923 6.553384 230.2 2161.7 3916.4 669.8 3334.1 0.7 4.4 17.5 173.1 38,830.4 0.26798 3.557059 6.565799 367.6 1799.8 4490.6 1246.9 516.1 0.7 3.9 11.6 79.7 31,749.09 0.070573 Ijede 3.607656 6.58398 417.6 718.4 1265.4 795.3 551.5 0.6 5.7 20.4 63.2 16,084.51 0.104959 3.691796 6.658793 371.4 415.8 1257.1 900.3 318.2 0.7 5.4 18.7 61.3 16,240.29 0.097836 3.588603 6.572281 404.3 587.9 1189.3 872.1 257.4 0.6 4.7 16.3 57.9 14,577.18 0.089576 Lagasa (orubga) 3.774026 6.623492 380.1 625 412.2 2017.3 744.2 0.7 3.7 10.4 23.3 10,773.71 0.044542 3.773612 6.633256 312.2 612.5 321.7 2124.5 644.2 0.6 3.1 11.1 19.8 9544.408 0.050807 3.771944 6.645947 298.1 557.9 204.6 1981.8 664.2 0.5 4.1 11.3 20.5 9075.064 0.0746 Ikosi-Agbowa 3.676676 6.604238 312 308.2 1012.6 92.1 277.7 0.8 2.5 4.8 29.2 1713.98 0.332463 3.707469 6.627901 167.1 901 147.9 1679.4 1.1 5.9 52.9 3330.905 0.370805 3.71195 6.639205 559.5 155.8 61.4 468.5 0.9 3 19.1 535.9225 0.331939 3.706407 6.613382 329.9 155.1 45.3 353.3 1.2 2.3 10 301.4025 0.239217 Ejirin-Ikosi 3.875877 6.633158 412.4 799.2 1161.7 1308.2 286.7 0.6 2.9 10.6 82.6 24,587.29 0.077348 3.899557 6.613982 362.3 507.6 1076.5 1873.8 486.1 6059.5 0.6 2.2 6.6 12.7 41.9 8767.308 0.105095 3.884937 6.618626 307.5 868.6 1257.9 1060.6 834.1 2881.9 0.7 2.8 8.5 30 60.7 15,964.56 0.113316 3.72148 6.651098 380.1 625 2317.4 3638 8118 744.2 0.5 4 10.4 23.2 44.6 78,875.9 0.024074 Orinsha 3.932677 6.609692 473.8 747.2 787.6 6004.3 468.8 0.6 3.1 8.9 34.3 43,111.55 0.022428 3.955684 6.609008 471.6 842.4 2026.7 7884.6 881.6 0.6 8 15 56.3 96,265.13 0.025311 3.924728 6.61318 408.2 361.2 1626.3 6390.5 535.5 0.7 4 11 60.9 81,760.26 0.029083 Epe 4.002733 6.589904 589.5 70.2 11.8 63.2 308.9 0.6 2.2 13.5 43.4 682.064 1.863134 3.964254 6.664254 329.3 86 19.9 89.6 357 0.7 2.9 19.3 55 1158.396 1.619535 3.968749 6.62913 204.7 91.9 20.1 7.5 188.3 0.8 2.8 9 18.5 148.146 2.948804 3.902899 6.617087 220.1 232.9 32.7 9.6 121.4 217 0.6 1.3 6.8 9.5 56.6 1269.938 1.67207 Ibeju-Lekki 3.862969 6.506399 43.3 13 53.3 15.5 178.3 0.5 2.2 3.9 22.9 122.614 1.731368 3.922679 6.467515 14.6 17.3 54.4 15.5 102.6 0.6 1.9 4.9 34.7 169.208 2.479706 3.716066 6.444968 18.4 63.2 142.5 7.3 69.2 0.8 2.5 5 8.1 188.87 1.227712 3.733927 6.493425 7.9 30.6 42.1 21.1 108.4 0.8 3.1 8 21.5 178.326 1.411554 Eputu-Lakowe 3.691344 6.481131 3419.2 388.9 140.6 53.7 428.4 0.6 2.9 17 46 1607.946 0.985154 3.612471 6.510547 5703.9 651.6 92.6 557.7 113.5 609.1 0.7 1.9 1.8 8.4 62.6 2864.538 0.589081 3.581578 6.491106 3920.6 521.6 59.2 720.8 35.4 587 0.6 1.6 2.1 8.6 25.1 1733.11 0.759664 3.565397 6.471926 2034.8 233.4 37.6 11.9 97.4 0.7 2.8 9.9 5.8 503.828 0.763033 Lambasa 3.634962 6.44286 3429.1 640.4 123.9 1359.1 58.2 0.7 2.4 19.1 56.1 16,509.87 0.199386 3.624312 6.483479 3282 376.5 69.3 476.8 87.2 0.7 2.9 18.4 63.4 6978.698 0.406398 3.670724 6.472587 2860.1 465.5 147.3 188.8 29.1 1 3 25.5 110.9 5790.134 0.767304 3.573707 6.477972 2135.2 814.2 136.4 63.2 59.5 53.1 0.7 2.6 12.1 29.3 1422.752 0.555838 Sangotedo 3.634962 6.44286 2426.2 447.9 176.1 235.5 67.8 0.7 2.8 7.3 88.2 5001.818 0.422516 3.624312 6.483479 3000.2 332.8 58.3 223.7 49.9 0.7 3 15.4 75.6 4181.616 0.611351 3.670724 6.472587 4211.7 630 128.4 195.9 33.1 0.6 3.2 11.1 62.1 3626.73 0.408669 3.640231 6.462005 3214.2 578.8 95.9 254.1 33.6 0.6 3 45.1 96.4 6497.05 0.85503 Ajah 3.573707 6.477972 1424.7 314.4 102.6 34.4 327 0.8 4.2 21.3 67.1 1390.772 2.172104 3.557353 6.435873 999.7 271.2 79.2 24.8 117.9 0.7 3.1 14.9 69.4 888.342 2.998649 Eti-Osa 3.482543 6.432149 1866 3822 87.6 58.1 2371.1 0.7 3.5 6.2 61.9 3764.542 1.137472 3.57924 6.438695 1636.3 550.5 214.1 126.7 929.1 0.7 3.6 11.3 49.4 2361.104 0.449644 3.599389 6.443611 2371.5 1130.3 175.9 65.4 1378.4 0.9 3.3 12.5 139.5 3437.278 2.20739 3.639255 6.428793 1008.8 345.1 135.4 35.1 800.8 1.1 5.1 14.2 187.1 2271.916 5.451228 3.565304 6.430002 1966.5 405.1 137.3 67.8 1439.1 0.7 4.2 21.6 81.1 2308.446 1.364209 Ikoyi Airtel 3.38397 6.461299 83.7 364.7 24.3 4.4 457.4 0.6 3.7 15.9 75.1 423.284 17.73982 3.460691 6.4533 74.9 312.9 23.1 5.7 421.3 0.6 4.1 19.7 63.7 429.198 12.04937 3.430051 6.444843 87.3 402.4 32.1 6.1 541.8 0.6 2.9 17.5 70.7 442.472 12.14941 3.41324 6.456833 81.9 360.7 26.5 5.4 473.5 0.6 3.5 17.8 69.8 432.042 13.61465 Iwopin 4.163984 6.489433 424.4 50.1 8.4 37 151.4 0.7 2.8 29.7 83.6 756.008 5.852711 4.199152 6.520966 210.9 72.6 558.1 79.3 458.7 0.8 2.7 12.6 34.7 2029.702 0.501139 4.222945 6.53625 213.6 32.6 6.4 52.7 127.6 0.8 3.7 18.9 59.9 713.838 4.20699 4.21158 6.566172 318 39 7.2 2 27.5 0.7 3.6 14.3 18 100.392 11.08062 Oworoshoki 3.403984 6.548789 25 28 10.7 247.2 0.4 5.3 21 95.775 2.167903 3.403101 6.559842 52.3 82.7 7.1 315.3 751.1 0.6 4.1 19.7 23.9 1609.198 2.911498 3.399997 6.569001 43 15.8 63.7 6.7 106.1 0.5 2.6 2.4 19.4 69.088 3.109384 3.404661 6.548052 53.1 30.7 117.3 21.8 214.5 0.8 2.4 4.6 25.1 240.584 1.283834 Tin can Island 3.357869 6.432569 1424 3907.6 903.2 115.2 40.2 369.3 1.2 0.3 8.8 31.4 73.5 3480.244 2.11159 3.346282 6.439306 247.5 37.6 10.7 191 289.2 1290.7 0.6 2.2 1.4 8.4 51.1 3325.744 0.412449 3.354608 6.449967 233.3 36.8 7.6 3.3 51.5 349.4 0.7 2.7 9.6 5.2 46.1 545.388 3.810431 3.337184 6.432654 285.1 47.2 8.8 1.9 22.8 96.1 0.8 3.1 10.9 23.3 113.7 621.39 18.55712 makoko-Okobaba 3.381124 6.499842 15.9 44.8 3.8 1.1 12.4 0.5 1.5 7.1 15.1 23.748 15.66062 3.367749 6.478124 48.9 112.8 15.2 2.1 197.2 0.7 2.1 9.3 18.4 90.222 9.406679 3.376114 6.493259 40.2 10.8 5 108.1 0.7 2.2 13.4 29.725 2.901117 3.372416 6.492654 87.1 13.1 7.8 211.1 510.3 0.8 2.7 5.6 14.9 658.824 1.003823 Ilaje Bariga 3.383287 6.535416 25.2 3 0.8 18.8 42.9 0.7 1.9 6.9 14.5 60.292 10.05739 3.393219 6.543001 37.9 7.6 78.2 123.1 412.2 0.8 1.9 9.1 23.7 734.77 0.580003 (Continued) 428 A. A. OLUDAYO ET AL. the transverse resistivity (Tr) and the Longitudinal conductance (Lc) were used in this study. Tr was used to identify target areas with good groundwater and is directly associated with groundwater transmis sivity. Tr was obtained from the first order geoelectric parameters, which are the layer’s apparent resistivity and the thickness. Equation 1 was used in calculating the transverse resistance (Tr) values with the layer’s resistivity and h the layer thickness for n-layers (Zohdy et al. 1974). Tr ¼ Xn i� 1 h1 � ρ1 þ h2 � ρ2 þ . . .þ hn � ρn Ωm2� � (1) Also, the longitudinal conductance (Lc) value was obtained from the acquired VES results (ρnadhÞ: The Lc was estimated by using equation 2. Lc ¼ Xn i� 1 h1 ρ1 þ h2 ρ2 þ . . .þ hn ρn mhosð Þ (2) The result from both Tr and Lc was interpolated with the spatial tools in the ArcGIS Platform using the inverse distance weight (IDW) interpolation scheme and their classification is presented in Table 2. Moreover, the map was classified according to its influence on groundwater quality. Areas with high Tr were given a high class and vice versa as presented in Table 3. (Continued). NAME? LONG LAT R1 R2 R3 R4 R5 R6 H1 H2 H3 H4 H5 TR LC 3.40099 6.513157 10.2 1.6 0.5 5.4 24.1 0.7 2.3 5.9 8.6 12.042 14.89872 3.397736 6.536293 48.7 10.1 2.1 178.2 552.7 0.9 2.4 12.7 18.2 667.596 6.405856 Orile iganmu 3.36368 6.47843 40.6 79.2 11.3 87.3 293.1 1 2.3 38.2 72 1388.004 4.258944 3.347323 6.480751 71.1 140.2 11.9 146.9 737.4 0.9 2.1 21.8 114.1 3475.824 2.636288 3.343693 6.48483 41.3 97.2 10.8 239 688.3 1 2.3 41 94.2 4644.292 4.238314 3.350823 6.478089 78.9 136.4 8.9 335 7.7 0.8 1.2 7.7 127.5 8601.566 1.264703 Erun Kan 3.412928 6.60654 75.3 18 4.5 34.1 215.7 0.7 2.9 2.5 30.7 232.606 1.626256 3.430786 6.61202 61.6 9.7 32.5 105.1 0.8 6.1 27.5 200.44 1.488007 3.425746 6.61467 54.8 37.8 26.1 97.7 350.8 0.7 4.3 21.2 42.9 989.11 1.37789 3.422636 6.617255 90 25.8 9.5 44.5 156.6 0.7 2.6 9.6 37.2 375.336 1.955034 Akoka Unilag 3.388017 6.512451 30.3 16.2 51.6 99.4 24.9 0.7 2.6 6.3 27.6 626.37 0.583355 3.389232 6.513717 14.6 15.7 57.1 10.1 286.4 0.6 3 14.4 48.3 273.186 5.266546 3.396077 6.513099 92.9 26.9 189.4 27.6 8.5 0.7 4.3 8.9 77.5 801.072 3.022348 3.38796 6.518927 308.8 15.4 98 115 45 0.7 1.5 6 53.6 1398.252 0.626981 Ita-Oko Island 4.192854 6.432275 60.7 98.9 761.3 1716.9 3006.4 555.4 0.6 3.4 13 53.1 105.8 69,919.02 0.127458 4.207758 6.433615 51.9 236.7 1554.6 4309.6 3386.7 2298.9 0.6 2.9 6.7 26.2 166.8 114,824.4 0.083453 4.159693 6.427474 170.4 96.4 1253.1 1618.9 6594.8 0.7 2.7 7.8 19.3 8279.702 0.050263 4.15241 6.41178 47.5 436.8 201.1 1891.4 6674.5 0.6 1.9 12.5 52.2 20,420.65 0.106738 Ikosi-ketu 3.38058 6.61595 190.1 724 65.4 343.6 857.8 0.6 1.1 3.3 35.9 2692.304 0.159616 3.385437 6.625939 111.5 934.9 130.1 3396.8 3797.8 0.4 1.5 7.2 48.2 33,221.89 0.074724 3.401268 6.639044 565.4 36.1 403.6 885.6 4439.3 0.5 3 16.7 65.1 12,956.74 0.198874 3.388452 6.630058 685.9 415.8 468.7 871.3 765.9 0.7 3.2 13.9 49.9 10,360.7 0.095644 3.384336 6.639651 612.8 341.6 487.3 944.8 697.3 0.6 3.1 15.8 50.9 11,443.26 0.096351 Figure 3. Frequency distribution for qualitative analyses of the VES curve types. NRIAG JOURNAL OF ASTRONOMY AND GEOPHYSICS 429 (Table 2). Finally, the calculated parameter of the longitudinal conductance (Lc) and transverse resis tance (Tr) is presented in Table 3. 4.2. Analytical Hierarchy Process (AHP) The AHP technique was used as a measurement for both subjective and objective evaluation of decision-makers, coupled with a provisional technique that checks and validates the consistency of the evaluated measurements. For this study, the eight contributing thematic para metric layers that influence groundwater suitability were ranked, then assigned weights by their importance. The AHP selects the likelihood of uncertainty in assign ments via the principal Eigenvalue and the consistency index (Saaty 2004). A parameter weighted high was associated with the layer of high impact and the para meter weighted low was associated with a small ground water impact. The weightage of each parameter and the pair-wise comparison was assigned based on standard Saaty’s scale 1–9 (Table 4) of relative importance value, with each factor, scored in accordance to its relative influence (Saaty 1980). The Saaty’s scale of relative importance presented in Table 4 show that each pair with a maximum difference is having one factor that is nine (9) times important than the other. The pair-wise comparison matrix for the selected groundwater para meter was developed by estimating the bulk value in every column in Table 5. The Sum of the columns was divided by every element in the matrix in other to nor malise the matrix in Table 5. Also, in furtherance of validating the judgement consistency, the principal Eigenvalue was derived by multiplying the right matrix of the judgement by the normalised Eigenvector. The result of the pair-wise comparisons matrix (Weights, λmax) is summarised in Table 6. Then, the consistency index was evaluated using equation 4. λmax ¼ X1 n λ n ¼ 7:333þ 9:003þ 7:948þ 7:820 þ9:754þ 10:377þ 9:826þ 8:529 8 ¼ 8:8239 (3) Where n is the number of factors (i.e. 8) and λ is the average eigenvalue of the consistency vector. CI ¼ ðλmax � nÞ n � 1 ¼ 8:8239 � 8 8 � 1 ¼ 0:1177 (4) Also, the Consistency Ratio (CR) was used to measure the consistency of the pairwise matrix by using equation 5. Figure 4. Drainage map of area around the Lagos Lagoon. Table 4. Saaty’s 1–9 scale of relative importance (Saaty 1980). Scale Interpretation 1 Requirements i and j are equal value 3 Requirement i has a slightly higher value than J 5 Requirement i has a strongly higher value than j 7 Requirement i has a very strong higher value than j 9 Requirement i has an absolutely higher value than j 2,4,6,8 These are intermediate scales between two adjacent judgements Reciprocals If requirement 1 has a lower value than j Table 5. Pair wise comparison matrix of 7 criteria for the AHP process. i/j Tr Ge Ld Rf Dd Sl Lc Gm Tr 1 2 3 4 4 3 3 5 Ge 1/2 1 1 3 3 1/3 1 1 Ld 1/3 1 1 2 1 2 1 2 Rf 1/4 1/3 1/2 1 1 1/2 1/5 1 Dd 1/4 1/3 1 1 1 1/4 1/3 3 Sl 1/3 3 1/2 5 4 1 1/2 2 Lc 1/3 1 1 5 3 2 1 5 Gm 1/5 1 1/2 1 1/3 1/2 1/5 1 Sum 3.200 9.667 8.500 19.000 17.333 9.583 7.233 20.00 430 A. A. OLUDAYO ET AL. CR ¼ CI RI ¼ 0:1177 1:41 ¼ 0:0835 (5) Where RI is the Ratio index calculated by Saaty (2004), using Table 7. For the comparison value to be consid ered acceptable, the consistency ratio must be smaller than or equal to 0.1. This implies that the judgement is consistent, and the derived weights can be used. For this study, n is 8, and CR is 0.0835. 4.3. Delineation of groundwater suitability zones The groundwater suitability zones (GSZ) were deli neated by the weighted linear combination (WLC) tech nique. (WLC) is a deterministic measure used to integrate non-parametric physical criteria in the ArcGIS module. The WLC was to calculate the ground water suitability index by assigning weights of relative importance and then combine the standardised weights to design a suitability map. In order to generate a map that will portray a suitable area for potable water, all the contributing eight thematic parametric layers were inte grated into the ArcGIS module and the overall weight of every map was computed using Equation 6: GPZ ¼ GmW � Gmrð Þ þ Lcw � Lcrð Þ þ Slw � Slrð Þ þ Ddw � Ddrð Þ þ Rfw � Rfrð Þ ¼ Ldw � Ldrð Þ þ Gew � Gerð Þ þ Trw � Trrð Þ (6) Where w is the normalised weight and r is the influen cing factor ratings. The weighted average was reclassi fied into five values to create a raster output file. A model builder was used to model all the influencing parameters. Each influencing parameter was repre sented with a band in the ArcGIS environment by using the normalised weights of each thematic layers. The vector thematic layer was transformed to raster thematic layers which later was reclassified for the optimisation of the integrated thematic map, then the groundwater suitability zone was estimated (Table 2). Moreso, the Groundwater Suitability Index Map was produce using the weight overlay model as presented in Figure 13. The groundwater suitability index of the study area was classified into Very Low, Low, Moderate, and High suitability using the natural breaks classification method. 4.4. Validation of groundwater suitability zones Finally, the GSZ model was validated with physico chemical data. Seventy-two (72) groundwater sam ples were collect sporadically within the study area with depths ranging from 28 to 87 m (Figure 1) to test for their electrical conductivity (EC) and total dis solved solids (TDS) concentrations. The water sam ples were taken to the chemistry laboratory for EC and TDS analysis using the American Public Health Association (APHA) (2005) standard. Salinity in groundwater is associated with the amount of dis solved salt in water. When salt dissolves in water, their ions dissociate and increases the concentration of TDS and EC in groundwater (EPA. 2012). The statutory standard measures from World Health Organization (World Health Organization (WHO) 2018) and Nigerian standard for drinking water qual ity (NSDWQ 2015) for guidance on drinking water (Table 8) were used to corroborate the result of the EC and TDS. The descriptive statistical result of the EC and TDS was presented in Table 8 and the elec trical conductivity (EC) and total dissolved solids (TDS) values were interpolated with inverse distance weight to produce the EC and TDS maps. Then, the EC and TDS maps were compared and correlated with the groundwater suitability map to validate the GSZ map, and to verify the effectiveness of the AHP method in demarcating GSZs in the study area. A very low electrical conductivity in groundwater is Table 6. Determining the relative criterion weights. i/j Tr Ge Ld Rf Dd Sl Lc Gm Normalised value Eigen value Tr 0.312 0.207 0.353 0.211 0.231 0.313 0.415 0.250 0.286 7.333 Ge 0.516 0.103 0.118 0.158 0.173 0.034 0.138 0.050 0.116 9.003 Ld 0.104 0.103 0.117 0.105 0.058 0.209 0.138 0.100 0.117 7.948 Rf 0.078 0.035 0.059 0.052 0.058 0.052 0.027 0.050 0.051 7.820 Dd 0.078 0.035 0.117 0.052 0.058 0.026 0.046 0.150 0.070 9.754 Sl 0.104 0.310 0.059 0.105 0.231 0.104 0.069 0.100 0.135 10.377 Lc 0.104 0.103 0.118 0.263 0.173 0.209 0.138 0.250 0.170 9.826 Gm 0.063 0.103 0.059 0.052 0.019 0.052 0.028 0.050 0.530 8.529 Table 7. Saaty’s ratio index for different values of n. n 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 RI 0.00 0.00 0.58 0.90 1.12 1.24 1.32 1.41 1.45 1.49 1.51 1.48 1.56 1.57 1.59 Table 8. Descriptive statistical result of the physico-chemical analyses. Parameter Minimum Maximum Mean WHO (2018) SON (2006) EC (μS=cmÞ 126 2380 925.79 1000 1000 TDS ðmg=lÞ 113 2264 481.03 500 500 NRIAG JOURNAL OF ASTRONOMY AND GEOPHYSICS 431 caused by the absence of impurities in groundwater. For groundwater to properly conduct electricity, there must be ions contained within it (EPA. 2012). Also, natural impacts like evaporation and rain, and human impacts which include agricultural runoff, road salt, and septic leachate) are factors that dictate the fluctuation of EC in groundwater (EPA. 2014). 5. Results and discussion Drainage Network The estimated drainage density in the study area varies from 1 to 318,900 km/km2 and the drainage patterns are dendrites (i.e. a tree root-like). High drainage density will make the surface runoff high and lower the infiltration rate, while low drainage density zones indicate a high rate of infiltration and probably high recharge zones (Vasanthavigar et al. 2011; Shailesh et al. 2019). The contributing influence of drainage density is low, with an overall percentage of 12.5% as presented in Table 2. The drainage density map presented in Figure 4, shows that lower density values are poor zones for groundwater accumulation while higher density value zone indicates good subsur face permeable materials and therefore enhances the rechargeability of groundwater Lineament extraction: The analysed lineament result is presented in the lineament density map with some lineament lines in Figure 5. Concentrated linea ment density zone indicates high permeable and porous zones, which enhances groundwater accumulation (Fashae et al. 2014). The lineament networks are criss crossed by minor to major lineaments lines varying in length. The lineament rose diagram has a bearing of 284:3o with a population density of 2398 that covers a total polygon area of 68.69 Km2. The rose diagram reveals three noticeable directions (NW-SE, N-S, and NE-SW) but in different ranges. The dominant orienta tion of the lineament direction tends towards the NW- SE direction. The lineament density map is classified into five categories: very low to low (lineament density less than 0.104 km/Km2), moderate (0.320–0.562 km/ km2), and high to very high (0.56 − 1.76 km/km2) as shown in Table 2. Zones with a higher value of linea ment density are weighted with high value and vice versa. The orientations of the lineament lines are indi cative of the preferential flow path (Khalil et al. 2014). Geological structures: The ability of the subsurface media to transports and accumulates groundwater is built upon its geological structures, mineral content, and lithological texture (Salako et al. 2019). The major geological features of the study area are the Coastal Plain sand and Coastal Alluvia which belongs to the Quaternary to the Recent Formation that has an alter nating sequence of sand, clay, lignite, and silt/clay (Jones and Hockey 1964). While the minor geological units are the Ilaro formation, Ewekoro formation, Oshosun formation, and Abeokuta complex basin which is the shale, sandstone, limestone, marls gastro pod, and siltstone. High explorations of groundwater in the alluvia belt are proportional to the high rate of upward movement of saline water in this zone (Longe 2011). The geological map of the study area is presented in Figure 6. The contributing influence of the geological structure has an overall percentage of 17.5% as pre sented in Table 2. More so, the geology of the study area was reclassified (Table 2) by their influence on saltwater intrusion. Figure 5. Lineament map and lineament density map of areas around the Lagos Lagoon. 432 A. A. OLUDAYO ET AL. Hydro-geomorphological map: Geomorphology reflects the various structural behaviour, nomencla ture processes of landforms and landscapes. Geomorphological features control the distribution of a surface water body and groundwater situation. The geomorphology of the research zone is cate gorised into coastal alluvium, beach ridge, mangrove swamp flats, pediplain, floodplain, a barrier beach complex, coastal creeks, and lagoons as shown in Figure 7 (Offodile 2014). The contributing influence of the geomorphology is given an overall percentage of 5 and the geomorphology features are classified and weighted according to their increasing significant influence on groundwater suitability as presented in Table 2. The southern zone of the study area is domi nated by Barrier beach, beach ridge Mangrove swamp, and Alluvium plain which are lowland areas with low drainages and high infiltration rate while the Flood plain and pediplain are the higher land region with high drainages and low infiltration rate. Figure 6. Geology map of area around the Lagos Lagoon. Figure 7. Geomorphology map of areas around the Lagos Lagoon. NRIAG JOURNAL OF ASTRONOMY AND GEOPHYSICS 433 Slope and elevation: The slopes are terrain para meters that give a clue on the geologic nature and geodynamic process (Arulbalaji et al. 2019). It also expresses the steepness on which surface runoff and infiltration depend. The larger the slope, the smal ler the recharge. This is because the flow of water is rapid down the steep slope during rainfall, thereby giving rise to low infiltration. The geospatial dis tribution of the slope steepness is presented in Figure 8. Areas with a steep slope (2:5o � 34:0o) are characterised by high runoff and low rainwater infiltration rate while, areas with a gentle slope (0:01o � 2:5o) allow a high infiltration rate due to slow runoff and promote appreciable groundwater recharge. This implies that the slope steepness is inversely proportional to the infiltration of surface runoff as presented in Table 2. Rainfall: Rainfall is one of the most influencing hydrological parameters that contribute to the groundwater recharge zone The average rainfall in the study location ranges from 135.039 to 219.778 mm. As presented in Figure 9, the amount of rainfall decreases towards the shore while the slope gradient directly influences the infiltration rate of water runoff. The rainfall data were reclassified into three classes as presented in Table 2. Figure 8. Slope map of areas around the Lagos lagoon. Figure 9. Mean annual rainfall map of areas around the Lagos lagoon. 434 A. A. OLUDAYO ET AL. Geoelectric Resistivity Result: The iterated VES result of some selected areas around the Lagos Lagoon is presented in Figure 10. The RMS errors were less than 4.6. The maximum and the minimum values at the study location are between 7469 Ωm northern zones and 1.1 Ωm at the southern zone respectively with a mean resistivity value of 469.79 m. The representative curve types (Figure 10) of the geoelectric layer signifies the alternating sequence of conductive to resistive layers of the subsurface and reflecting the unconsoli dated nature of their lithologic composite and its decomposition characteristics. Seventeen curve types are obtained after plotting 197 partial curves. Figure 3 shows the frequency of the occurred curve types with QHA, QQH, and QKH having the highest frequency 17.65%, 11.76%, and 10.59% distribution respectively for the five-layer types and KQ and AK and KH are the four-layer curve types and the six-layer curve type is QHKQ. For the study area, four to six geoelectric layers are delineated. The first layer is the topsoil having a resistivity value that varies from 0.8 to 8118 Ωm, with thicknesses between 0.4 and 1.6 m. The second layer is characterised as a shallow unconfined aquifer (Asseez 1972; Longe 2011; Oyedele 2019) due to litho- sequence of sand, sandy silt, and very low imperme able clay. The second layer has a resistivity value ran ging from 1.2 to 3907.6 Ωm, with thicknesses between 0.3 and 8 m. Lower resistivity (<30 Ωm) is observed around Bariga, Agboyi, Irawo, Ibeju Lekki, Oworoshoki, and Makoko regions while high resistiv ity value (>500 Ωm) is observed at Bayiku, Agbowa, Figure 10. Representative inverted resistivity curve types of selected area around the Lagos Lagoon. NRIAG JOURNAL OF ASTRONOMY AND GEOPHYSICS 435 Ejirin, Tin-can, and Ajaganbe environments. The higher resistivity values demarcated free polluted aquifer zones and the lower resistivity value shows the aquifer contamination description (Oyedele 2019). Although, some of the low resistivities areas towards the northern zone of the study location are characterised as clay and peat (Ejirin ikosi, Epe, and Orinsha). The third geo-electric layer has resistivity values in the range of 0.8 to 4490.6 Ωm with a mean thickness of 20.45 m indicating the mixture of sand and lateritic clayey sand. The fourth geo-electric sequence is characterised as the confined coastal plain aquifer, separated by impermeable sandy clay to silty clay and clay (Oyedele 2019). The resistivities value varies from 1.1 to 7884.6 Ωm with a mean thick ness of 51.2 m at an average depth of 72.8 m. (Longe et al. 1987) noted that the aquifer layer at this depth is a multi-layer system with three distinctive sub-zones. The fifth and sixth geoelectric layers were delineated with resistivity values ranging from 53.1 to 8118 Ωm with a signature of the basement complex aquifer zone (Coode et al. (1996), Nton and Adeyemi 2014), although, the depth at this point could not be deter mine due to the termination of the current. Therefore, their depths and thicknesses at that point could not be determined. Most of the confined aquifer is located at the second to third strata of the subsurface and is dominant towards the north-eastern zone of the research location. The aquifer in this region shows high protective capacity as the thickness of clay in this region is high compared with the southwestern zone, which has low clay thickness. Geoelectrical parameter: Environmental degrada tion such as surface water depletion and saltwater intrusion could be due to excessive abstraction of groundwater within the belt of the Lagos Lagoon and its proxy to the Atlantic Ocean. Furthermore, geo electric parameters such as the longitudinal conduc tance, transverse resistance, thickness, and depth of the geological strata are essential in evaluating groundwater suitability. The result from Dar Zarrouk parameters was interpolated in the ArcGIS platform with the help of an inverse distance weight (IDW) module to generate the distributed spatial map. Also, the result of both Tr and Lc was presented in Table 3. Transverse resistance: In this study, the transverse resistance values vary from 11.12 Ωm2 to 137,789 Ωm2 with a mean value of 13,419.78 Ωm2 (Table 3). However, Zohdy et al. (1974), Suggest that an increase in Transverse resistance may be due to an increase in the resistive thickness of the subsurface. High trans verse resistance (Tr) indicates regions with high aqui fer transmissivity and porous media which can enhance the prospect contaminant flow into the groundwater (Tsepav et al. 2015; Senthikumar et al. 2019). Thus, the transverse resistance spatial distribu tion map of the research area is presented in Figure 11. A high value of Tr is observed (>1000 Ωm2) at the northern and northwest zone around the lagoon shore, as well as a low value of Tr less than 500 Ωm2 is noticed at the eastern and southern zone. Longitudinal conductance: The longitudinal con ductance (Lc) values span between 0.025 and 25.54 mhos with an average value of 2.85 mhos (Table 3). The contrasting longitudinal conductance along the coastal belt is indicative of Atlantic marine sediments (Evans et al. 2017). Relatively, high longitudinal con ductance values (above 10 mhos) are observed at Ilaje Bariga, Agboyi Ketu Owode-Ibese, and Erun-kan which scrutinise the presence of Brackish to saline Figure 11. Transverse resistance map of areas around the Lagos Lagoon. 436 A. A. OLUDAYO ET AL. water and clay to slit dominated regions (Ayolabi et al. 2013). The average value of Lc was observed around areas close to the lagoon. This could indicate the brackish to saline groundwater nature as presented in Figure 12 (Laouini et al. 2017). The low conduc tance values, noticed around the central north and south zones, could indicate the influence of the Ewekoro and Ilaro formation (Yusuf et al. 2018; Akinlalu and Afolabi 2018). The high value of Tr and low Lc around the central northern part of the research location reveals the existence of a coastal plain sand formation. Groundwater Suitability Zone (GSZ): The GSZ, is a dimensionless model that can identify the suitability of the groundwater zone within the study area map by the assignment and normalising weight to each influ encing thematic layer (Figure 13). The groundwater suitability index map is classified into Very Low, Low, Moderate, and High suitability. The groundwater Suitability map shows that an excellent GSZ is con centrated around the central-western part (Ikosi– Agbowa, Lagasa, Ejirin–Ikosi Ajaganabe, and Ijede) of the study area due to the geological influence of Ilaro, Ewekoro, and the Oshosun formation with the high transverse resistance (Coode et al. 1996). In this research, 7% of research locations have very low suit ability, 32.2% were low, 47.8% e moderate and 13% have high groundwater suitability. The transverse resistance and the geology feature play a very signifi cant part in the suitability of groundwater in this Figure 12. Longitudinal conductance map of area around the Lagos Lagoon. Figure 13. Groundwater suitability index map of area around the Lagos Lagoon. NRIAG JOURNAL OF ASTRONOMY AND GEOPHYSICS 437 research with an influencing capacity of 37.5%. The drainage density, lineament density, rainfall, and slope also help in the rechargeability of the aquifer, aquifer mobility rate, and the steepness in which surface run off and the infiltration capability of the groundwater resources. The geomorphology and longitudinal con ductance have a reasonable effect of 15% on the phy sical features of geological structures and the aquifer protective capacity of saltwater infiltration (Jaroslav and Balthazar 2011). Groundwater Suitability Zones Validation: A validation assessment that demonstrates that the model is a reasonable representation of the actual aquifer system of the study area. In doing this, the groundwater suitability zone was validated with results obtained from the physicochemical analysis of groundwater. The descriptive statistical result of the EC and TDS is presented in Table 8. The EC and TDS results show different variations across the study area with a minimum value of (126 (µS/cm) and 113 (mg/ l)) and maximum value of (2380 (µS/cm) and 2264 (mg/l)) respectively. The permissible limit of EC and TDS stipulated by World Health Organization (WHO) (2018) and NSDWQ (2015) are 1000 (µS/ cm) and 500 (mg/l) respectively. The distribution maps of EC and TDS are presented in Figures 14 and Figure 14. Electrical conductivity map of areas around the Lagos Lagoon. Figure 15. Total dissolve solid map of areas around the Lagos Lagoon. 438 A. A. OLUDAYO ET AL. 15 respectively. High EC and TDS are observed towards the southern zone and, the western north zone of the study location. The high EC and TDS might be due to the influence of the Ogun river at the western north and the proximity to the Atlantic Ocean in the south zone of the study location respec tively. It was observed that low EC and TDS are situated around the central north-eastern zone (Ikosi–Agbowa, Lagasa, Ejirin–Ikosi Ajaganabe, Orinsha, and Epe), probably due to high altitude and presumably lower water table with lesser groundwater abstraction (UNICEF 2008; Oyem et al. 2014). Similarly, the low concentration EC and TDS within this location ascertained that the GSZ model imple mented is reasonable and the observation from the geological and geomorphological analysis have been correctly verified. The existence of the Ilaro and Oshosun formation has a high protective capacity that produces a shield that protects the aquifer from saltwater intrusion (Ayolabi et al. 2013). The EC and TDS measure of groundwater around this zone (southern to the eastern) corroborate to an agreeable degree to the GSZ model. Though there is a trivial significant difference between the GSZ model and the EC /TDS, yet both maps agree to a large extent. 6. Conclusion and recommendation Both the AHP and GIS approaches have been used to integrate different geophysical, hydrological, and remote sensing maps to produce groundwater suit ability zone around the Lagos coastal area. The inte grated groundwater suitability zone (GSZ) map was categorised by assigning normalised weights to differ ent thematic maps. The AHP was used for the pairwise comparison matrix for the various thematic layers like geology, geomorphology, lineament density, drainage density, slope, rainfall, transverse resistance, and long itudinal conductance. The integrated groundwater suitability zone map of the research locations reveals that 7% of research locations have very low suitability, 32.2% were low, 47.8% e moderate and 13% have high groundwater suitability. Subsequently, EC and TDS maps generated from the 72-borehole samples were used to validate the groundwater suitability model and to check the quality of the water. The GSZ model corroborated with EC/TDS maps. The GSZ model shows that the northern zone of the study locations is a suitable area for groundwater exploration. The GSZ map will assist in the decision-making of the state to plan a sustainable groundwater monitoring scheme. Finally, the study has been used to establish the efficacy of integrating approach between, MCDA (AHP), GIS-based module, remote sensing, rainfall precipitation, and geophysical method in effectively mapping groundwater suitability zone for human consumption. Disclosure statement No potential conflict of interest was reported by the author(s). ORCID Akintunde A Oludayo http://orcid.org/0000-0002-8642- 0481 Ozebo V Chidi http://orcid.org/0000-0001-5062-0548 Oyedele F Kayode http://orcid.org/0000-0002-5711- 766X References Adenodi RA. 2018. A centurial analysis of rainfall variability in Nigeria. Nigerian Journal of Technology (NIJOTECH). 37(2):543–547. doi:10.4314/njt.v37i2.34. Adeoti L, Alile OM, Ughegbulam O. 2010. Geophysical investigation of saline water intrusion into freshwater aquifers: a case study of Oniru, Lagos State. Scientific Research and Essays. 5(3):248–259. (1992-2248). Adiat KAN, Nawawi MNM, Abdullah K. 2012. Assessing the accuracy of GIS-based elementary multi criteria deci sion analysis as a spatial prediction tool—a case of pre dicting potential zones of sustainable groundwater resources. J Hydrol. 440-441:75–89. doi:10.1016/j. jhydrol.2012.03.028 Agagu, OA. 1985. A geological guide to bituminous sedi ments in Southwestern Nigeria. Department of geology, University of Ibadan Akinlalu AA, Afolabi DO. 2018. Borehole depth determination to freshwater and well design using geophysical logs in coastal regions of Lagos, southwestern Nigeria. Applied Water Science. 8(6):152–169. doi:10.1007/s13201-018- 0798-3. Akintunde OA 2018. Hydro-geophysical analyses and estimation of aquifer hydraulic conductivity, overbur den vulnerability and transmissivity at Ejinrin-Ijebu, Lagos. Presented at University of Lagos 13th Annual Research Conference and Fair. (163–171) Lagos, Nigeria. Aladejana JA, Robert M, Kalin RM, Sentenac P, Hassan I. 2020. Hydrostratigraphic characterisation of shallow coastal aquifers of Eastern Dahomey Basin, S/W Nigeria, using integrated hydrogeophysical approach; implication for saltwater intrusion. Geosciences. 10 (2):65. doi:10.3390/geosciences10020065. Allen JRL. 1965. A review of the origin and characteristics of recent alluvial sediments. J Int Assoc Sedimentology. 5 (2):89–191. doi:10.1111/j.1365-3091.1965.tb01561.x. NRIAG JOURNAL OF ASTRONOMY AND GEOPHYSICS 439 https://doi.org/10.4314/njt.v37i2.34 https://doi.org/10.1016/j.jhydrol.2012.03.028 https://doi.org/10.1016/j.jhydrol.2012.03.028 https://doi.org/10.1007/s13201-018-0798-3 https://doi.org/10.1007/s13201-018-0798-3 https://doi.org/10.3390/geosciences10020065 https://doi.org/10.1111/j.1365-3091.1965.tb01561.x American Public Health Association (APHA). 2005. Standard methods for the examination of water and wastewater. 21st. Washington (DC., USA): American Public Health Association/American Water Works Association. Anthony MK. 2016. Combining geophysical techniques and multi-criteria GIS-based application modelling approach for groundwater potential assessment in southwestern Nigeria. Environ Earth Sci. 75(16):1181–1201. doi:10.1007/s12665-016-5897-6. Arulbalaji P, Padmalal D, Sreelash K. 2019. GIS and AHP techniques based delineation of groundwater potential zones: a case study from Southern Western Ghats, India. Sci Rep. 9(1):2082. doi:10.1038/s41598-019-38567-x. Asseez LO. 1972. Hydrogeology of southwestern Nigeria. The Nigerian Engineers. 7(1):22–44. Ayolabi EA, Folorunso AF, Odukoya AM, Adeniran AE. 2013. Mapping saline water intrusion into the coastal aqui fer with geophysical and geochemical techniques: the University of Lagos campus case (Nigeria). Springerplus. 2(1):433–449. doi:10.1186/2193-1801-2-433. Bachaer A, Ikram J, Samir S, Salem B. 2018. The seawater intrusion assessment in coastal aquifers using GALDIT method and groundwater quality index: the Djeffara of medenine coastal aquifer (Southeastern Tunisia). Arabian Journal of Geosciences. 11(20):609. doi:10.1007/s12517- 018-3966-8. Badejo OT, Olaleye JB, Lademomi AS. 2014. Tidal charac teristics and sounding datum variation in Lagos State. International Journal of Innovative Research and Studies. 3: 1–23. 2319-9725. Bagyaraj M, Ramkumar T, Venkatramanan S, Gurugnanam B. 2012. Application of remote sensing and GIS analysis for identifying groundwater potential zone in parts of Kodaikanal Taluk, South India. Frontiers of Earth Science. 7(1):65–75. doi:10.1007/s11707-012- 0347-6. Batayneh AT. 2013. The estimation and significance of Dar-Zarrouk parameters in the exploration of quality affecting the Gulf of Aqaba coastal aquifer systems. Journal of Coastal Conservation. 17(3):623–635. doi:10.1007/s11852-013-0261-4. Billman HG 1976. Offshore stratigraphy and palentology of the Dahomey Embayment. Proceeding of the 7th Afican Micropaleontology Colloquium, Ile-Ife. Burke K, Dessauvagie TE, Whiteman A. 1971. Opening of the Gulf Of Guinea and geological history of the Benue depression and Niger-Delta. Nature, Phys Sci. 233 (38):51–55. doi:10.1038/physci233051a0. Chawla JK, Khepar SD, Sondhi SK, Yadav AK. 2010. Assessment if long-term groundwater behaviour in Pujab, India. Water Int. 35(1):63–77. doi:10.1080/ 02508060903513502. Coode B, Oteri AU, Rofe KL (1996) Hydrogeological inves tigation of Lagos State. Final report, vol. I, II. Demiroğlu M, Dowd J. 2014. The utility of vulnerability maps and GIS in groundwater management: a case study. Turkish Journal of Earth Science. 23:80–90. doi:10.3906/yer-1205-6 Edet AE, Okereke CS, Teme SC, Esu EO. 1998. Application of remote-sensing data to groundwater exploration: a case study of the cross-river state, southeastern Nigeria. Hydrogeol J. 6(3):394–404. doi:10.1007/s100400050162. EPA. 2012. Conductivity. In Water: Monitoring and Assessment. Retrieved from http://water.epa.gov/type/ rsl/monitoring/vms59.cfm EPA. 2014. Sediments. In Water: Pollution Prevention & Control. Retrieved from http://water.epa.gov/polwaste/ sediments/ European Commission (1995) Soil terrain database. Land management and natural hazards unit. IES and JRC, European Commission, Brussels. Available at http:// eusoils.jrc.ec.europa.eu/projects/SOTER/Soter_Model. html. Accessed Dec 2016. Evans UF, Abdulsalam NN, Mallam A. 2017. Natural vul nerability estimate of groundwater resources in the coastal area of Ibaka community, using dar zarrouk geoe lectrical parameters. J Geol Geophys. 6:295–303. doi:10.4172/2381-8719.1000295 Fashae OA, Tijani MN, Talabi AO, Adedeji OI. 2014. Delineation of groundwater potential zones in the crystal line basement terrain of SW-Nigeria: an integrated GIS and remote sensing approach. Applied Water Science. 4 (1):19–38. doi:10.1007/s13201-013-0127-9. GAF D, Kaki C, Adeoye JA. 2016. Benin and Western Nigeria Offshore basins: a stratigraphic nomenclature comparison. International Journal of Geosciences. 7 (2):177–188. doi:10.4236/ijg.2016.72014. Gyeltshen S, Tran TV, Gunda GKT, Kannaujiya S, Chatterjee RS, Champatiray PK. 2020. Groundwater potential zones using a combination of geospatial tech nology and geophysical approach: case study in Dehradun, India. Hydrological Science Journal. 65 (2):169–182. doi:10.1080/02626667.2019.1688334. Henriet JP. 1976. Direct application of Dar-Zarrouk para meters in ground water surveys. Geophysical Prospecting. 24(2):344–353. doi:10.1111/j.1365-2478.1976.tb00931.x. Heywood I, Cornelius S, Carver S. 1998. An introduction to geographical information systems. New Jersey: Prentice Hall. Jaroslav V, Balthazar TV. 2011. Groundwater of emergency situations: a methodological guide. United Nations Educational, Scientific and Cultural Organization. 3:173–185. Jones HA, Hockey RD 1964. The geology of Part of South-Western Nigeria. Geological survey of Nigeria Bulletin No. 31, Geological Survey of Nigeria, Kaduna. Kadam AK, Kale SS, Pande NN, Pawar NJ, Sankhua RN. 2012. Identifying potential rainwater harvesting sites of a Semi-arid, Basaltic region of Western India, Using SCS-CN Method. Water Resour Manage. 26 (9):2537–2554. doi:10.1007/s11269-012-0031-3. Khalil M, Ahmed K, Elnahry A, Hasan A. 2014. Integrated geophysical, remote sensing and GIS studies for ground water assessment, Abu Zenima Area, West Sinai, Egypt. International Journal of Geosciences. 5(9):882–907. doi:10.4236/ijg.2014.59078. Laouini G, Etuk SE, Agbasi OE. 2017. Delineation of aqui fers using Dar Zarrouk parameters in parts of Akwa Ibom, Niger Delta, Nigeria. Journal Hydrogeol Hydrol Eng. 6(1):18–42. doi:10.4172/2325-9647.1000151. Longe EO. 2011. Groundwater resources potential in the coastal plain sands aquifers, Lagos, Nigeria Research. Journal of Environmental Earth Science. 3(1):1–7. 440 A. A. OLUDAYO ET AL. https://doi.org/10.1007/s12665-016-5897-6 https://doi.org/10.1038/s41598-019-38567-x https://doi.org/10.1186/2193-1801-2-433 https://doi.org/10.1007/s12517-018-3966-8 https://doi.org/10.1007/s12517-018-3966-8 https://doi.org/10.1007/s11707-012-0347-6 https://doi.org/10.1007/s11707-012-0347-6 https://doi.org/10.1007/s11852-013-0261-4 https://doi.org/10.1038/physci233051a0 https://doi.org/10.1080/02508060903513502 https://doi.org/10.1080/02508060903513502 https://doi.org/10.3906/yer-1205-6 https://doi.org/10.1007/s100400050162 http://eusoils.jrc.ec.europa.eu/projects/SOTER/Soter_Model.html http://eusoils.jrc.ec.europa.eu/projects/SOTER/Soter_Model.html http://eusoils.jrc.ec.europa.eu/projects/SOTER/Soter_Model.html https://doi.org/10.4172/2381-8719.1000295 https://doi.org/10.1007/s13201-013-0127-9 https://doi.org/10.4236/ijg.2016.72014 https://doi.org/10.1080/02626667.2019.1688334 https://doi.org/10.1111/j.1365-2478.1976.tb00931.x https://doi.org/10.1007/s11269-012-0031-3 https://doi.org/10.4236/ijg.2014.59078 https://doi.org/10.4172/2325-9647.1000151 Longe EO, Malomo S, Olorunniwo MA. 1987. Hydrogeology of Lagos Metropolis. Afr J Earth Sci. 6(2):163–174. Medina-Gómez I, Kjerfve B, Mariño I, Herrera-Silveira J. 2014. Sources of salinity variation in a coastal Lagoon in a Karst landscape. Estuaries and Coasts. 37(6):1329–1342. 10.1007/s12237-014-9774-9. Retrieved from http://www. jstor.org/stable/44851210 Mohamed HK, Khalid SA, Alaa Eldin HE, Alaa NH. 2014. Integrated geophysical, remote sensing and GIS studies for groundwater assessment, Abu Zenima Area, West Sinai, Egypt. International Journal of Geosciences. 5 (9):882–907. doi:10.4236/ijg.2014.59078. Nigeria Geological Survey Agency (2006) Published by the authority of the federal Republic of Nigeria. NSDWQ. 2015. Nigerian standard for drinking water quality. Niger Ind Stand. 554:13–14. Nton ME, Adeyemi MO. 2014. Petrography, compositional characteristics and stable isotope geochemistry of the Ewekoro formation from Ibese Corehole, eastern Dahomey basin, southwestern Nigeria. Global Journal of Geological Sciences. 13(1):35–52. doi:10.4314/gjgs.v13i1.5. Offodile ME (2014) Hydrology: groundwater study and development in Nigeria. ISBN:978-30956-4-1. Olorode DO, Adedayo AS, Akintunde AO. 2016. Application of geotechnical and geophysical methods to investigate tilt buildings at Lagos State, Nigeria. IOSR Journal of Applied Geology and Geophysics (IOSR-JAGG). 4(5):21–28. doi:10.9790/0990-0405022128. Olorunfemi MO, Oni AG. 2019. Integrated geophysical methods and techniques for siting productive bore holes in basement complex terrain of Southwestern Nigeria. Ife Journal of Science. 21(1):013. doi:10.4314/ijs.v21i1.2. Omatsola ME, Adegoke OS. 1981. Tectonic evolution and cretaceous stratigraphy of the Dahomey Basin. Journal of Mining and Geology. 18:130–137. Oyedele AA. 2019. Use of remote sensing and GIS tech niques for groundwater exploration in the basement complex terrain of Ado-Ekiti, SW Nigeria. Applied Water Science. 9:51. doi:10.1007/s13201-019-0917-9. Oyedele KF, Oladele S. 2011. Geoelectrical assessment of groundwater potential in the coastal aquifer of Lagos, Nigeria. Advances in the Research of Aquatic Environment. 2:29–36. doi:10.1007/978-3-642-24076- 8_4 Oyem HH, Oyem IM, Ezeweali D. 2014. Temperature, pH, electrical conductivity, total dissolved solids and chemical oxygen demand of groundwater in Boji-Boji Agbor/Owa area and immediate suburbs. Research Journal of Environmental Sciences. 8(8):444–450. doi:10.3923/ rjes.2014.444.450. Ozebo VC, Ajiroba SO. 2011. Groundwater assessment in apapa coast-line area of Lagos using electrical resis tivity method. Journal of Emerging Trends in Engineering and Applied Sciences (JETEAS). 2 (4):673–679. (2141-7016). Pinto D, Shrestha S, Babel MS, Ninsawat S. 2015. Delineation of groundwater potential zones in the Comoro watershed, Timor Leste using GIS, remote sensing and analytic hierarchy process (AHP) technique. Applied Water Science. 7(1):503–519. doi:10.1007/s13201-015-0270-6. Prasad RK, Mondal NC, Banerjee P, Nandakumar MV, Singh VS. 2008. Deciphering potential groundwater zone in hard rock through the application of GIS. Environ Geol. 55(3):467–475. doi:10.1007/s00254-007-0992-3. Recep Ç. 2019. Evaluation of groundwater potential by GIS-based multicriteria decision making as a spatial pre diction tool: case study in the tigris river Batman-Hasankeyf Sub-Basin, Turkey. Water. 11 (12):2630. doi:10.3390/w11122630. Saaty TL 1980. The analytic hierarchy process. McGraw- Hill, New York. reprinted by RWS Publications, 4922 Ellsworth Avenue, Pittsburgh, PA, 15213, 2000a. Saaty TL. 2004. Decision making the analytic hierarchy and network processes (AHP/ANP). Journal of Systems Science and Systems Engineering. 13(1):1–35. doi:10.1007/s11518-006-0151-5. Salako AO, Osotuyi AG, Adepelumi AA. 2019. Seepage investigations of heterogeneous soils beneath some build ings using geophysical approaches: example from south western Nigeria. Geo-Engineering. 10(11):17–29. doi:10.1186/s40703-019-0107-5. Scibeka J, Allena DM, Cannonb AJ, Whitfield PH. 2007. Groundwater-surface water interaction under scenarios of climate change using a high-resolution transient groundwater model. Journal of Hydrology. 333(2– 4):165–181. doi:10.1016/j.jhydrol.2006.08.005. Senthikumar S, Vinodh K, Johnson Babu G, Gowtham B, Arulprakasam V. 2019. Integrated seawater intrusion study of coastal region of Thiruvallur district, Tamil Nadu, South India. Applied Water Science. 9 (5):124–144. doi:10.1007/s13201-019-1005-x. Shaban A, Khawlie M, Abdallah C. 2006. Use of remote sensing and GIS to determine recharge potential zones: the case of Occidental Lebanon. Hydrogeol J. 14 (4):433–443. doi:10.1007/s10040-005-0437-6. Shailesh KS, Malte Z, Ude S, Griffiths GA. 2019. Potential groundwater recharge zones within New Zealand. Geoscience Frontiers. 10(3):1065–1072. doi:10.1016/j. gsf.2018.05.018. Singh P, Thakur JK, Kumar S. 2013. Delineating ground water potential zones in a hard-rock terrain using geos patial tools. Hydrol Sci J. 58(1):213–223. doi:10.1080/ 02626667.2012.745644. SON 2006. Nigerian Industrial Standard for Potable Water developed by Standards Organisation of Nigeria and National Guidelines and Standards for Water Quality in Nigeria. ICS 13.060.20 Surajit M, Ramakar J 2015. Identification of groundwater potential zones using remote sensing and gis in a mine area of Odisha, Odisha, India. National Conference on recent approach to water resource management. https:// DOI:10.13140/RG.2.1.3374.6644 . Tsepav MT, Ibrahim SI, Bayegun FA. 2015. Geoelectrical char acterization of aquifer precincts in parts of lapai, north central nigeria. journal of applied science and environmental. Management. 19(2):295–301. doi:10.4314/ jasem.v19i2.17. UNICEF. 2008. UNICEF handbook on water quality. New York (USA): United Nations Children’s Fund (UNICEF). Vander Velpen BPA (2004). Win RESIST Version 1.0. M. Sc Research Project. ITC, Deft (Netherlands). NRIAG JOURNAL OF ASTRONOMY AND GEOPHYSICS 441 https://doi.org/10.1007/s12237-014-9774-9 https://doi.org/10.1007/s12237-014-9774-9 http://fromhttp://www.jstor.org/stable/44851210 http://fromhttp://www.jstor.org/stable/44851210 https://doi.org/10.4236/ijg.2014.59078 https://doi.org/10.4314/gjgs.v13i1.5 https://doi.org/10.9790/0990-0405022128 https://doi.org/10.4314/ijs.v21i1.2 https://doi.org/10.1007/s13201-019-0917-9 https://doi.org/10.1007/978-3-642-24076-8_4 https://doi.org/10.1007/978-3-642-24076-8_4 https://doi.org/10.3923/rjes.2014.444.450 https://doi.org/10.3923/rjes.2014.444.450 https://doi.org/10.1007/s13201-015-0270-6 https://doi.org/10.1007/s00254-007-0992-3 https://doi.org/10.3390/w11122630 https://doi.org/10.1007/s11518-006-0151-5 https://doi.org/10.1186/s40703-019-0107-5 https://doi.org/10.1016/j.jhydrol.2006.08.005 https://doi.org/10.1007/s13201-019-1005-x https://doi.org/10.1007/s10040-005-0437-6 https://doi.org/10.1016/j.gsf.2018.05.018 https://doi.org/10.1016/j.gsf.2018.05.018 https://doi.org/10.1080/02626667.2012.745644 https://doi.org/10.1080/02626667.2012.745644 https://DOI:10.13140/RG.2.1.3374.6644 https://DOI:10.13140/RG.2.1.3374.6644 https://doi.org/10.4314/jasem.v19i2.17 https://doi.org/10.4314/jasem.v19i2.17 Vasanthavigar M, Srinivasamoorthy K, Vijayaragavan K, Gopinath S, Sarma S. 2011. Groundwater potential zon ing in Thirumanimuttar sub-basin Tamilnadu, India (A GIS and remote sensing approach. Geo-Spatial Inf Sci. 14 (1):17–26. doi:10.1007/s11806-011-0422-2. Whiteman A (1982) Nigeria: its petroleum geology, resources and potential. Graham and Trotman, 394. Wilson RCC, Willians CA. 1979. Oceanic transform struc tures and the developments of Atlantic continental mar gin sedimentary basin a review. Journal of Geological Society of London. 136(3):311–320. doi:10.1144/ gsjgs.136.3.0311. World Health Organization (WHO). 2018. Guidelines for drinking water quality. First Addendum 3rd1 491–493 Geneva. Yusuf MA, Abiye TA, Butler MJ, Ibrahim KO. 2018. Origin and residence time of shallow groundwater resources in Lagos coastal basin, south-west Nigeria: an isotopic approach. Heliyon. 4(11):e00932. doi:10.1016/j.heli yon.2018.e00932. Zohdy AAR, Eaton GP, Mabey DR (1974) Application of surface geophysics to groundwater investigations, techni ques water resources investigations of the US geological survey. Washington, 195–205. 442 A. A. OLUDAYO ET AL. https://doi.org/10.1007/s11806-011-0422-2 https://doi.org/10.1144/gsjgs.136.3.0311 https://doi.org/10.1144/gsjgs.136.3.0311 https://doi.org/10.1016/j.heliyon.2018.e00932 https://doi.org/10.1016/j.heliyon.2018.e00932 Abstract 1. Introduction 2. Location, hydrology, and geology the study area 3. Materials and methods 4. 3.1: Geo-spatial data acquisition 4.1. Geophysical field survey 4.2. Analytical Hierarchy Process (AHP) 4.3. Delineation of groundwater suitability zones 4.4. Validation of groundwater suitability zones 5. Results and discussion 6. Conclusion and recommendation Disclosure statement ORCID References