Sustainable Groundwater Potential Zoning with Integrating GIS, Remote Sensing, and AHP Model: A Case from North-Central Bangladesh

: Groundwater is one of the most valuable natural resources, and the most dependable source of fresh water. For sustainable groundwater management, the present study aimed to model groundwater potential zones in the north–central region of Bangladesh using GIS, remote sensing, and the analytical hierarchy process. The present study included eight thematic layers: lineament density, geomorphology, soil types, slope, land use/land cover, drainage density, elevation, and rainfall features to delineate a groundwater potential zone of the area. Integration of the eight thematic layers was performed through weighted overlay analysis, which assisted in delineating groundwater potential zones. This simple and systematic method successfully provides a satisfactory result concerning the delineation of groundwater potential zones. The study resulted in a groundwater potential zone map, which identiﬁes about 11.51% of the study area as being under a very high groundwater potential zone, covering an area of 504.09 km 2 . The AHP analysis shows that the physiographical parameters, such as lineament density, slope, and drainage density, and meteorological factors such as annual rainfall, have greater inﬂuence over groundwater potentiality. The result obtained from the weighted overlay analysis was veriﬁed with actual well yield and groundwater depth data, which show a signiﬁcant positive correlation. The outcome of the study will help in taking effective measures to ensure sustainable use and extraction of groundwater in this region. Contributions: Conceptualization, M.A.S. M.A.I.; methodology, M.A.S., M.A.I., U.P. M.N.-E.-A.; software, M.A.I., M.N.-E.-A., S.K.S. U.P.; validation, M.A.S., A.R.M.T.I. M.F.U.; formal analysis, M.A.I., M.N.-E.-A., S.K.S. U.P.; investigation, M.A.I. and U.P.; resources, S.I.I., A.R.M.T.I. M.F.U.; data curation, M.A.S., M.A.I., A.R.M.T.I. S.K.S.; writing—original preparation, U.P., M.A.I., M.N.-E.-A. and M.A.S.; writing—review and editing, A.R.M.T.I., M.F.U., A.E.R. and S.K.S.; visualization, M.A.S., A.R.M.T.I., A.E.R. and S.K.S.; supervision, M.A.S.; project administration, M.A.S.; funding acquisition, M.A.S., A.E.R., M.A.I., M.N.-E.-A., U.P.


Introduction
Nature's greatest gift, water, is essential to human civilization, primarily for drinking, irrigation, industry, and other purposes. Almost three-quarters of the Earth's surface is covered with water, but freshwater resources are limited [1]. Groundwater plays an important role in meeting freshwater demands, particularly for drinking water [2]. As a primary source of clean, potable water in all climatic regions, groundwater is recognized as one of the most valuable natural resources [3]. Therefore, groundwater supplies irrigation The study area of Mymensingh is in the north-central zone of Bangladesh. Mymensingh district is located at latitudes of 24 • 15 to 25 • 12 N and longitudes of 90 • 04 to 90 • 49 E. The district has an area of approximately 4396.53 km 2 , with several small valleys between high forests, and accounts for approximately 2.96% of Bangladesh's total area [37], as shown in Figure 1. The city of Mymensingh is located on the bank of the Old Brahmaputra River. In its north-western corner, near Karaibari, the Brahmaputra River enters Mymensingh, then flows towards the south-east and south to the Meghna, near Bhairab Bazar [38]. At a width of 94 m, the Jamuna forms the western edge of Mymensingh, and during the rainy season, this river widens to 5 to 6 m in many areas, while the Meghna flows just over the south-eastern part of the district for a limited distance [39]. Mymensingh has a moderate climate with a humid, oppressive, and mostly cloudy wet season, and a warm, mostly clear dry season [40]. It has a subtropical monsoon climate, with rainfall from May to August and the average annual rainfall is 2541 mm [41]. The region is dominated by recent alluvial flood plain deposits, with underlying Pleistocene Madhupur clay deposits [39]. Mymensingh has a total of 635 perennial bodies of water, covering an area of 49.43 km 2 [42]. Here, the shallow aquifers are characterized by a near-neutral pH, moderate COD, and high arsenic (As) and iron (Fe) levels [4]. The district of Mymensingh is heavily reliant on groundwater. Agricultural irrigation, and other domestic and industrial activities, in this area meet their water needs using groundwater harvesting. Therefore, the decline in the groundwater level has become a vital problem in this area [43]. and suitable sites for tube well installation. The study's findings will also aid in the implementation of efficient strategies to ensure the sustainable use and extraction of groundwater in this region.

Description of the Study Area
The study area of Mymensingh is in the north-central zone of Bangladesh. Mymensingh district is located at latitudes of 24°15′ to 25°12′ N and longitudes of 90°04′ to 90°49′ E. The district has an area of approximately 4396.53 km 2 , with several small valleys between high forests, and accounts for approximately 2.96% of Bangladesh's total area [37], as shown in Figure 1. The city of Mymensingh is located on the bank of the Old Brahmaputra River. In its north-western corner, near Karaibari, the Brahmaputra River enters Mymensingh, then flows towards the south-east and south to the Meghna, near Bhairab Bazar [38]. At a width of 94 m, the Jamuna forms the western edge of Mymensingh, and during the rainy season, this river widens to 5 to 6 m in many areas, while the Meghna flows just over the south-eastern part of the district for a limited distance [39]. Mymensingh has a moderate climate with a humid, oppressive, and mostly cloudy wet season, and a warm, mostly clear dry season [40]. It has a subtropical monsoon climate, with rainfall from May to August and the average annual rainfall is 2541 mm [41]. The region is dominated by recent alluvial flood plain deposits, with underlying Pleistocene Madhupur clay deposits [39]. Mymensingh has a total of 635 perennial bodies of water, covering an area of 49.43 km 2 [42]. Here, the shallow aquifers are characterized by a nearneutral pH, moderate COD, and high arsenic (As) and iron (Fe) levels [4]. The district of Mymensingh is heavily reliant on groundwater. Agricultural irrigation, and other domestic and industrial activities, in this area meet their water needs using groundwater harvesting. Therefore, the decline in the groundwater level has become a vital problem in this area [43].

Description of Data Processing
Various analogue maps were obtained from various organizations during the early stages of GIS spatial database development (Table 1). Rainfall data was collected from ten Sustainability 2022, 14, 5640 4 of 24 (10) Bangladesh Water Development Board (BWDB) sampling stations for the year of 2020, and interpolated into the GIS system. The soil map was prepared from the "digital soil map of the world" published by the Food and Agriculture Organization (FAO). Land use, land cover and lineament density data were generated from downloaded Landsat-8 imagery from the USGS Earth Resource Observation System Data Center [https: //earthexplorer.usgs.gov (accessed on 7 February 2021)]. The map of land use and land cover was generated using supervised classification based on non-parametric rules, and then a theoretical confusion matrix (error matrix) analysis [44] was performed to confirm accuracy using Google Earth and Google Maps, as shown in Table 2.   (5) Waterbodies  0  0  0  13  1  14  1  Total  14  11  11  30  1  67  53 Final accuracy of LULC classification was evaluated by the following Equation [44]: Accuracy (%) = number of correct sampled/total number of samples = (53/67) × 100 = 79.10% Multi-software processing of Landsat-8 data was used to create a lineament density map. The lineament of the research area was initially derived from Landsat-8 data (Band-04) in PCI Geomatics (demo version). Then, using ArcGIS' spatial analysis tool, a map of lineament density (km/km 2 ) was created. All images were downloaded at level 1 geotiff format, and have a 30 × 30 m spatial resolution, without the presence of cloud cover (4.17%). A geomorphological map of the study area was prepared from the "Quadrangle Map" published by the Geological Survey of Bangladesh (GSB), by digitizing and geo-referencing.
The drainage density map, elevation map, and slope map were prepared from the ASTER Global Digital Elevation Model V3 (DEM), with 30 × 30 m spatial resolution, available on the Earth data site of NASA. The model was validated using well yield data from the Bangladesh Water Development Board's (BWDB) eight (08) sampling stations. Various analyses and interpretations (digitization, conversion, interpolation, classification, reclassification, enhancement, filtering, and other GIS processing) of satellite data were performed to produce thematic maps in a GIS environment. Lineament density, drainage density, slope, soil, rainfall, elevation, geomorphology, and land use/land cover were all the layers derived from this 3-stage method of groundwater potential zone mapping.
For groundwater potential zone mapping, at first, all data layers were transformed to raster data, with the same 30 × 30 m spatial resolution. A pair-wise comparison was carried out among all the parameters, to determine the relative importance for groundwater potentiality of each parameter [45]. Each map was reclassified based upon those weights generated. Finally, a groundwater potential map was created by combining all the thematic layers using Arc GIS 10.8's weightage overlay module. A detailed examination of procedures used in this study is illustrated in Figure 2.
the Bangladesh Water Development Board's (BWDB) eight (08) sampling stations. Various analyses and interpretations (digitization, conversion, interpolation, classification, reclassification, enhancement, filtering, and other GIS processing) of satellite data were performed to produce thematic maps in a GIS environment. Lineament density, drainage density, slope, soil, rainfall, elevation, geomorphology, and land use/land cover were all the layers derived from this 3-stage method of groundwater potential zone mapping.
For groundwater potential zone mapping, at first, all data layers were transformed to raster data, with the same 30 × 30 m spatial resolution. A pair-wise comparison was carried out among all the parameters, to determine the relative importance for groundwater potentiality of each parameter [45]. Each map was reclassified based upon those weights generated. Finally, a groundwater potential map was created by combining all the thematic layers using Arc GIS 10.8's weightage overlay module. A detailed examination of procedures used in this study is illustrated in Figure 2.

Assigning Weight Using Multi-Criteria Decision Making Approach
Using AHP, a multi-criteria decision-making approach, it is possible to establish the relative importance of various thematic layers and their associated attributes [35,46]. This AHP approach employs the hierarchy to look at many attributes individually, in order to identify possible groundwater flow and storage regulators from competing parameters [36]. The importance of each parameter was graded on a scale of 1 to 9, using a pairwise comparison matrix. The weighting of all the influencing parameters was based on their reaction to the occurrence of groundwater and expert opinion [28]. A parameter with a high weight indicates a large influence on groundwater potential, whereas a parameter with a low weight represents a minor impact [47]. Each parameter's weightage was allocated using Saaty's scale of relative importance, with values of 1 to 9 [48][49][50][51]. Then, the weights were allocated based on previous studies, together with field experience. According to Saaty's relative importance scale, a value of 1 means equal importance, 2 weak, 3 moderate, 4 moderate plus, 5 strong, 6 strong plus, 7 very strong, 8 very strong plus, and 9 indicates extreme importance [50,51]. Relative weights are ascribed to each thematic layer, based on their categorization, as shown in Table 3.

Weight Normalization
Each row's weights were averaged to find the appropriate ranking, resulting in the normalized weight of each parameter, as shown in Table 4. Based on the study's findings, observed lineament density has a higher value than any other metric. Substantial lineament density indicates the prospect for high groundwater recharge, while low lineament density implies low groundwater recharge, and a low groundwater potential zone. There are more elements that affect GWPZ, listed in increasing order: lineament density > slope > rainfall > drainage density > land use/land cover > geomorphology > soil > elevation.

Principal Eigen Vector
Using the algorithm for the consistency index, the maximum normalized main Eigen vector values (max) were calculated for each parameter, given in Table 5. As examples, lineament density = 34.53 was multiplied by the total value determined within a pairwise comparison matrix (e.g., lineament density = 2.50), shown in Table 4, to obtain the weight of the first criteria. The remaining eight parameters were treated in the same way. Table 5 shows the major eigen vector values (max = 8.527) for determining the consistency index (CI), by adding these values together.

Consistency Index (CI) and Consistency Ratio (CR) Calculation
The consistency ratio was calculated to evaluate the appropriateness of the weights applied to each parameter. For the weight to be considered constant, the consistency ratio must be less than 0.10. Consistency ratio (CR) was calculated using the following Equation (1) [50]: where CI is the consistency index, and RI is the random index. For the present research, the random index value is 1.41, since eight parameters were employed, as per Table 6 [50]. The consistency index was calculated using the following Equation (2): where λmax is the principal eigenvalue of the comparison matrix, and n is the total parameters used in this study. The consistency ratio was calculated to evaluate the appropriateness of the weights applied to each parameter. For the weight to be considered constant, the consistency ratio must be less than 0.10. From the calculation, the consistency index (CI) value is 0.07532. The consistency ratio (CR) result is 0.0538, which is less than 0.1, which means that the given weight for each parameter is valid for further analysis.

Assigning Ratings for the Weighted Overlay Analysis
Appropriate ratings for the overlay approach were assigned based on the influence of each thematic layer, and their classes on groundwater potentiality. The ratings were developed using expert opinion, field study, and previous groundwater potential zone mapping studies. Ratings were assigned on a scale of 1 to 5, with 1 indicating very little influence, and 5 indicating very significant influence.
If lineament density is high in a region, groundwater potential is also high. As a result, regions with lineament density of 0.96 to 1.78 km/km 2 are assigned a 5 (very high) rating, as these areas have long, major faults. Another rating, 4 (high), is given to the lineament density showing local faults that are linked together (frequent faults). A Sustainability 2022, 14, 5640 8 of 24 rating of 3, 2, and 1 is given for 0.49-0.96 km/km 2 , 0.28-0.49 km/km 2 , 0 to 0.28 km/km 2 lineament density, respectively. In the research region, slopes are classified as very low (0-10), low (1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(14)(15)(16)(17)(18)(19)(20), moderate (2-30), high , and very high (>50) slopes. Flat slopes receive the highest rating (very low) because they hold water for a long period, allowing groundwater to easily refile through infiltration. The steep (very high) slope is assigned the lowest ranking, because it causes a lot of runoff and minimal groundwater recharge. The availability of water for infiltration is increased when the yearly rainfall is high. So, locations with 2443-2634 mm/year rainfall receive a high grade of 5. As a result, areas with annual rainfall between 2316 and 2442 mm, 2208 to 2315 mm, 2106 to 2207 mm, and 1964 to 2106 mm are assigned a rating of 4, 3, 2, and 1, respectively. As a high drainage density area causes more surface runoff, the region with the highest drainage density (1.98 to 3.18 km/km 2 ) receives the lowest rating of 1. As a result, high, moderate, low, and very low drainage density receive ratings of 2, 3, 4, and 5, respectively, as shown in Table 6. The ratings for land use and land cover are based on the quantity of water available for recharging, depending on the landform. Agricultural land and bodies of water are assigned the highest rating of 5, paved settlement areas receive the lowest grade of 1, and vegetation area has a rating of 3. For its low permeability, fallow land receives a rating of 2. There are several geomorphologies in this research area. Madhupur clay residuum is given a rating of 2, as these areas are less suitable for groundwater recharge. Young gravelly sand, lakes, and the Chandina alluvium floodplain all obtain a rating 5, since these formations are well-known for effective groundwater recharging. Alluvial silt, and alluvial silt and clay are rated 4 and 3, respectively. Marsh clay and peat, and the Dihing and Dupi Tila formation undivided are assigned a rating of 1, due to relatively lower permeability. Depending on the characteristics, the soil groups of Dystric Nitosols (Nd) are assigned the maximum rating of 5, while Eutric Gleysols (Ge), with greater silt and clay, receive the lowest rating of 4. The elevation value ranges from 3 to 12 m, receiving the maximum rating of 5, indicating that it is more likely to recharge groundwater. Based on their appropriateness for groundwater potential, elevations ranging from 12 to 15, 15 to 20, 20 to 27, and 27 to 75 m are assigned ratings of 4, 3, 2, and 1, respectively. Table 6 represents the weights and the ratings (rank) of the eight criteria with their area coverage (km 2 ). Table 7 depicts the list of parameters, together with their assigned weights, ratings, and area coverage.  Weighted overlay analysis [52] is an approach for generating an extensive analysis, by assigning a common scale of values to input components based on AHP analysis. Each input layer was reclassified to a common ratio scale of five specific classes. Using the weighted overlay tool of ArcGIS 10.8 software, the reclassified layers of lineament density, drainage density, slope, soil, rainfall, elevation, geomorphology, and land use/land cover, and their corresponding percentage influence on recharge, were integrated to produce a map of the spatial distribution of groundwater potentials within the Mymensingh district. By multiplying the cell values of each factor class by the factor weight, and summing the resulting cell values together, the weighted overlay analysis tool reclassified values in the input raster layers into a common evaluation scale of 1, 2, 3, 4, and 5; very low, low, moderate, high, and very high, respectively. The following Equation (1) [28,29,31,47,52] was used to measure ground water potential zones of the study area: where Rf = rainfall; Ge = geomorphology; Sl = slope; DD = drainage density; El = elevation; S = soil; LULC = land use/land cover; and LD = lineament density. The symbol "w" expresses the weightage of a thematic layer, and the symbol "r" expresses the rating of subclasses (the rank) in each layer.

Validation of Results
One of the most important criteria in scientific study is validation. The defined GWPZs were validated using the well yield of pump test data obtained from the Bangladesh Water Development Board (BWDB). Cross-verification analysis [53] was performed between the groundwater potential zone values and well yield data. A total of 8 well yield data (liter per second) were found and used for the investigation.
Again, we used groundwater depth data acquired from the monitoring wells of the Bangladesh Water Development Board (BWDB), and the Water Resources Planning Organization (WARPO), for the years 2019 and 2020, to validate our findings. A total of 68 monitoring well data points were obtained and interpolated in GIS, using the kriging interpolation method. Then, an accuracy assessment was performed using confusion matrix (error matrix) analysis between the resulting groundwater potential zones map and the observation wells data. The observation wells data were used as a point of reference for determining the accuracy of classification. The total accuracy was calculated using the following equation [54,55]: Accuracy (%) = (number of correct sampled/total number of samples) × 100 Then, kappa (K) analysis [56] was used, which is a multivariate way to evaluate accuracy that yields a Khat statistic. It was determined by applying the following equation [56,57]: (5) where TS stands for total samples, and TCS stands for total correct sampled.

Description of Influential Factors
Thematic maps were prepared at a scale of 1:50,000, with a spatial resolution of 30 × 30 m from satellite imagery, geomorphological mapping, soil mapping, and other hydrogeological field data. Thematic maps for each parameter are prepared as follows: 3.1.1. Lineament Density Lineament density in the study area ranges from 0 km/km 2 to 1.78 km/km 2 . The region is divided into five classifications to make it easier to provide a rating. Potentiality of groundwater depends largely on lineament density because it helps in groundwater recharge [57]. When the lineament density in an area is high, the potentiality of groundwater is also high. So, the rating of 5 (very high) is given to areas with a lineament density of 0.96 to 1.78 km/km 2 . This is because these areas have long, major faults. A little more than 6.05% of the area is covered by this class. Another rating, 4 (high), is given to the lineament density showing local faults that are linked together. It covers 17.03% of the study area, and has a value of 0.70 to 0.96 km/km 2 (frequent faults). In the middle of the scale, the lineament density ranges from 0.49 km/km 2 to 0.70 km/km 2 . Almost 25.68% of the research region is in this classification, and so this class receives a rating of 3, as it shows local faults with a lot of fractures. A total of 30.27% of the area has lineament density of 0.28-0.49 km/km 2 or less, which is assigned a rating of 2, because it shows fractures with short lines. Lineament density in these areas from 0 to 0.28 km/km 2 , covering about 21% of the total study area, receives a rating of 1. Lineament density is assigned 34.53% of the weight. Figure 3A shows lineament density in the study area. lineament density ranges from 0.49 km/km to 0.70 km/km . Almost 25.68% of the research region is in this classification, and so this class receives a rating of 3, as it shows local faults with a lot of fractures. A total of 30.27% of the area has lineament density of 0.28-0.49 km/km 2 or less, which is assigned a rating of 2, because it shows fractures with short lines. Lineament density in these areas from 0 to 0.28 km/km 2 , covering about 21% of the total study area, receives a rating of 1. Lineament density is assigned 34.53% of the weight. Figure 3A shows lineament density in the study area.

Slope
Slope categorization in the study area is categorized as very low (0-1 • ), low (1-2 • ), moderate (2-3 • ) slopes, as well as high (3-5 • ) and very high (>5 • ) slopes. Figure 3B shows the reclassified slope map of the study area. The highest rating is assigned to flat slopes (very low), because they hold water for a long time, meaning that groundwater is easily refiled by infiltration [58]. The lowest rank is given to steep (very high) slopes because they cause a lot of runoff and little groundwater recharge. A total of 39.38% of the study area is flat land with very little slope, 15.10% is medium slope, and 6.81% is very steep slope.

Rainfall
Rainfall is the most important source of water in the water cycle, as well as one of the most influential factors for groundwater availability. Rainfall intensity and duration have a big impact on how much water is infiltrated, and how much water runs off [59]. The hydrology of the study area mainly depends on monsoon rainfall [12]. Table 8 depicts the study area's yearly rainfall for the year 2020. The amount of rain that falls inside the north and north-east highland parts of the study area directly affect the rate at which groundwater is infiltrated in the central downstream area. The rainfall map of the study area is shown in Figure 3C.

Drainage Density
Drainage density is connected to permeability, which affects how much water runs off, and how much water soaks into the ground [60]. The drainage density is classified into five categories: very high, high, moderate, low, and very low drainage density zone, as illustrated in Figure 3D. The drainage density in the study area ranges from 0 to 3.18 km/km 2 . Approximately 9.56% of the overall land has a relatively low drainage density. Only 2.53% of the study area has a very high drainage density, ranging from 1.98 to 3.18 km/km 2 . Compared to a low drainage density region, a high drainage density area increases surface runoff. The maximum drainage density number reflects the greatest possibility of runoff, which eventually leads to less percolation. As a result, the lowest rating of 1 is given to the region with the highest drainage density, and a rating of 5 is given to the region with the lowest drainage density, 0 to 0.83 km/km 2 .

Land Use/Land Cover
The land use type of the study region determines information on soil moisture, permeability, amount of runoff, and percolation [61]. Mymensingh is a semi-urban region comprised of 43.38% of agricultural land, 3.36% of fallow land, 13.37% of water body mask, 7.64% of settlements, and 32.24% of vegetation, according to the LULC map, shown in Figure 3E. According to the results of the theoretical confusion matrix (error matrix) analysis, the LULC map has an accuracy of 79.10%, suggesting that the classified image is appropriate for further investigation. The ratings are assigned based on the amount of water available for recharging, depending on the landform. Agricultural land and waterbodies receive the maximum rating of 5, while paved settlement areas receive the lowest rating of 1, and vegetation area has a rating of 3. Fallow land receives a rating 2 due to its limited permeability. This layer is assigned a weight of 6.89%.

Geomorphology
The features of various water-bearing geomorphological formations influence the occurrence and flow of groundwater [62]. In this area, there are nine major geomorphologies: young gravelly sand, marsh clay and peat, Madhupur clay residue, lakes, Dihing and Dupi Tila formation undivided, Chandina alluvium, alluvial silt, alluvial silt and clay. Here, Madhupur clay residuum underlays 16.02% of the entire area, which is assigned a rating of 2, since these locations do not favor groundwater recharging. Young gravelly sand (7.4%), lakes (0.39%), and the Chandina alluvium floodplain (27.31%) each receive a rating of 5, as these formations are well-recognized for efficient groundwater recharge, and the Chandina alluvium floodplain occupies the majority of the research region. Alluvial silt (20.03%) and alluvial silt and clay (12.85%) have ratings of 4 and 3, respectively. Marsh clay and peat (14.68%) and Dihing and Dupi Tila formation undivided (0.94%) receive a rating of 1, due to extremely low permeability. This thematic layer is assigned a weight of 4.49%. Figure 3F shows the geomorphology map of the study area.

Soil Texture
The rate of infiltration is determined by the uppermost soil's permeability and water holding capacity [63]. Figure 3G shows the different types of hydrogeological soils classes in the study area. The largest part (73.19%) of the study area, is made up of Eutric Gleysols (Ge), with silty characteristics (containing mainly silt and clay). Another 26.81% of the region consists of Dystric Nitosols (Nd), with sandy characteristics. Depending on the characteristics, the highest rating of 5 is assigned to the soil groups of Dystric Nitosols (Nd), while the lower rating of 4 is assigned to Eutric Gleysols (Ge) with higher silt and clay. The overall weight for this layer is 2.87%.

Elevation
The elevation of the land is an essential aspect of groundwater recharge. Plains with lower altitudes tend to retain water longer, resulting in higher groundwater recharge [64]. Inside the northern and southern parts of the area, there are many big hills. Figure 3H shows that the elevation inside the study ranges from 3 to 75 m. The elevation is higher in the northern and southern parts of the study area. A total of 50.80% of the study area has elevation values of 3 to 12 m, and receives the highest rating, 5, as it is more likely to allow groundwater recharge. As follows, 12 to 15, 15 to 20, 20 to 27, and 27 to 75 m elevation, covering 30.18%, 15.66%, 3.04%, and 0.33% of the study area, respectively, are assigned the ratings of 4, 3, 2, and 1, based on their suitability for groundwater potential. The overall weight for this layer is 2.03%.

Groundwater Potential Zone
The eight thematic layers were generated in a GIS context, and the applicable rankings and weights, decided by the AHP technique, were allocated to them for the overlay analysis by combining all thematic maps. The findings are classified into five groups, based on the total values generated from Equation (3), ranging from very low potential zones to very high potential zones (Figure 4). The values of groundwater potential zones range from 4.49 to 338.50. Natural breaks classification is used to classify this groundwater potential zones (Table 9).

Validation of Groundwater Potential Zone
Cross-verification [53,64] was performed between the groundwater potential zone map and well yield data from the BWDB. A total of eight well yield data were found and used for the investigation, and their water production (liter per second) was extrapolated from the measured data. Well yield production (liter per second) varies from 25.21 to 37.56. This groundwater potential zone is in close (87.5%) accord with the results of the point source data provided in Table 10. High-yielding groundwater is found in several locations of the country. Rift faults in the region are likely responsible for this, since they may have created a range of degrees of displacement within rock formations that are now in touch with high-permeability rock types [62]. It is clear from the table that the findings achieved using this technique are in excellent agreement with real-time data. Therefore, the findings produced are more in line with the real field data.  Due to the flat surface of the research area, there is a higher likelihood of groundwater recharge. The northern, north-eastern, and central parts of the research area, all of which are inundated by water, have a greater possibility of having a zone of very high to high potential, covering 11.51% (504.09 km 2 ) and 22.74% of the study area, respectively (995.58 km 2 ). Areas where lineament density is higher, but the drainage density is lower, result in increased groundwater potentiality. Furthermore, in areas with low rainfall intensity and high drainage density, all of the water flows away, with little time to recharge the groundwater. In these locations, 5.375% (235.17 km 2 ) and 26.88% (1176.84 km 2 ) are classified as having very low to somewhat low potential zones. The eastern half of the research area has a low to moderate potential zone, with moderate rainfall, lineament, slope, and drainage density. A moderate potential zone, with an extent of up to 33.50%, covers the majority of the research area (1466.46 km 2 ).

Validation of Groundwater Potential Zone
Cross-verification [53,64] was performed between the groundwater potential zone map and well yield data from the BWDB. A total of eight well yield data were found and used for the investigation, and their water production (liter per second) was extrapolated from the measured data. Well yield production (liter per second) varies from 25.21 to 37.56. This groundwater potential zone is in close (87.5%) accord with the results of the point source data provided in Table 10. High-yielding groundwater is found in several locations of the country. Rift faults in the region are likely responsible for this, since they may have created a range of degrees of displacement within rock formations that are now in touch with high-permeability rock types [62]. It is clear from the table that the findings achieved using this technique are in excellent agreement with real-time data. Therefore, the findings produced are more in line with the real field data. Again, we used groundwater depth data acquired from the Bangladesh Water Development Board (BWDB), and the Water Resources Planning Organization (WARPO), for the years 2019 and 2020, to validate study's findings, as shown in Figure 5. A total of 68 monitoring well data were obtained and interpolated in GIS, using the kriging interpolation method. The average depth of the monitoring well ranges from 3.95 m to 10.37 m. Again, we used groundwater depth data acquired from the Bangladesh Water velopment Board (BWDB), and the Water Resources Planning Organization (WAR for the years 2019 and 2020, to validate study's findings, as shown in Figure 5. A tot 68 monitoring well data were obtained and interpolated in GIS, using the kriging inte lation method. The average depth of the monitoring well ranges from 3.95 m to 10.37 An accuracy assessment was performed between the resulting groundwater po tial zones map and the observed well data. The observation well data were used as a p of reference for determining accuracy of classification. A confusion matrix (error matrix) analysis was used to measure accuracy, as sh in Table 11. The total accuracy was calculated using the Equation (4): Accuracy (%) = ((59/68) × 100 = 86.76% After that, kappa coefficient (T) analysis was used, which is a multivariate wa evaluate accuracy that yields a Khat statistic. It was determined by applying Equation  An accuracy assessment was performed between the resulting groundwater potential zones map and the observed well data. The observation well data were used as a point of reference for determining accuracy of classification.
A confusion matrix (error matrix) analysis was used to measure accuracy, as shown in Table 11. The total accuracy was calculated using the Equation (4): Accuracy (%) = ((59/68) × 100 = 86.76% Table 11. Theoretical confusion matrix for total accuracy and Kappa coefficient analysis.

GWPZs
(1) After that, kappa coefficient (T) analysis was used, which is a multivariate way to evaluate accuracy that yields a Khat statistic. It was determined by applying Equation (5) × 100 The total accuracy and kappa coefficients are 86.76% and 83.28%, respectively, with a considerable level of agreement [65], indicating a good connection between groundwater potential zones and actual monitoring well information.

Discussion
Groundwater is the most abundant source of potable freshwater, as surface water sources are contaminated, or being contaminated, with anthropogenic contaminants [66][67][68][69]. Groundwater depletion is one of the world's most pressing worries, with the effects felt most acutely in underdeveloped countries [70]. Water consumption has risen dramatically in recent years, as a result of population increase and socioeconomic development [61]. Groundwater modeling is now regarded as a very efficient technique for managing groundwater resources. Groundwater recharge has considerably diminished, almost everywhere in the world, as a result of numerous anthropogenic activities and imbalanced development [71][72][73]. Knowing where groundwater potential zones are can aid in the more precise study of groundwater resources for sustainable use, as well as in the establishment of artificial groundwater recharge systems [29,74,75]. Groundwater is the principal source of fresh water in Bangladesh, and its agricultural sector is largely dependent on groundwater [76]. Within study region, a divisional city of Bangladesh, great strain on groundwater is experienced, due to the fast population growth, urbanization, and industrialization [77]. In the north-central region of Bangladesh, the majority of the tube wells were installed without regard to the hydrogeological features of the aquifers underneath them. This ignorance is currently hampering the performance of tube wells, because the level of groundwater has dropped below the suction limit [43]. This study analyzed hydrologic and geographic attributes of the north-central region of Bangladesh, covering eight major factors influencing groundwater potential: geomorphology, land use/land cover, drainage density, lineament density, soil type, slope, elevation, and annual rainfall. The hydrogeological information of the study area will help future researchers and policymakers in groundwater resource management. Similar studies for assessing groundwater potential zones were conducted worldwide [23,24,26,[28][29][30][31][32][33][34]74,[78][79][80], as it is a quick, cost-effective, and reliable way of predicting groundwater potentiality. In Bangladesh, similar studies covering the Atrai-Sib River basin [81], Dhaka city [82], agricultural area in Dinajpur district [83], the hilly terrain of Bangladesh [84], and drought-prone areas of Bangladesh [85] were performed, and no study was conducted in the Mymensingh district for groundwater potentiality. During the last few years, this district experienced intensive industrial activities. Many industries were built in this district, which created excessive pressure on ground water consumption. This study is useful in groundwater management, decision making, and as a guide for future detailed investigations in the study area. Based on the influence of groundwater potentiality, all the thematic layers were quantitatively placed together using GIS and multicriteria decision making approaches (AHP), and categorized the study area into very high, high, moderate, low, and very low groundwater potential zones. The modeled map was validated by groundwater depth data using a confusion matrix, and kappa coefficient analysis. This study has certain limitations, such as model validation using cross-verification relying on groundwater yield data from BWDB's eight monitoring stations, as there were only eight monitoring stations found from a reliable data source (BWDB) in the study region. To meet this gap, we also employed groundwater depth data from 68 monitoring wells in the study area for validation, using a confusion matrix and kappa coefficient. The total accuracy and kappa coefficients are 86.76% and 83.28%, respectively, with a considerable level of agreement [68,69], indicating a good connection between groundwater potential zones and actual monitoring well information. As per the study's findings, areas with flat surface, higher lineament density, lower drainage density, and higher annual rainfall show higher groundwater potentiality. The findings of the study can be used as a working document to address groundwater concerns in the study area. This document is useful for the industrial and agricultural sectors in selecting work sites and suitable sites for tube well installation. Based on the findings of this study, and earlier research, it is possible to conclude that integrated GIS and remote sensing techniques are immensely helpful, time-efficient, and cost-effective tools for modeling groundwater potential zones. Integrated remote sensing and GIS approaches are highly recommended for rapid, cost-effective, and reliable assessment in groundwater potential zone research.

Conclusions
At present, groundwater is an essential and valuable resource. The geographic information system (GIS) and remote sensing are remarkable and practical tools, as they are less time consuming and more cost-effective ways to discover groundwater potential zones in any region. In the present study, the analytical hierarchy process (AHP), along with GIS and remote sensing, was used to delineate the groundwater potential zones of the Mymensingh district (the north-central region of Bangladesh), based on the influential factors for groundwater potential zones. This simple and systematic method successfully provides a satisfactory result concerning the delineation of groundwater potential areas. The study results in a groundwater potential zone map for Mymensingh district that identifies that there is a higher potential for about 11.51% of the study area to be a high groundwater potential zone, which is an area of about 504.09 km 2 . The AHP analysis reveals that physiographical parameters, such as lineament density, slope, drainage density, geomorphology, and meteorological factors such as annual rainfall, have greater influence over groundwater potentiality. The groundwater level is depleting alarmingly in the study area, and making a lot of deep and shallow tube wells useless. To manage these valuable resources, we need to reduce our dependency on groundwater, by searching for alternatives. Proper management of groundwater recharge zones is also a very reliable way to enhance these resources. Agricultural activities in the study area directly depend on groundwater. Searching for alternative water resources is keenly needed for sustainable management of groundwater. Dredging the Old Brahmaputra River, and its adjacent canals, could be an important alternative in this regard. This study will help decision making for the sustainable management of groundwater resources, site determination for groundwater investigation, and exploitation. This study will also help develop natural or artificial recharge structures, for the prospect of groundwater use in the Mymensingh district to fulfil their water crisis during the dry season.