Comparative analysis of groundwater potentiality zone using fuzzy AHP, frequency ratio and Bayesian weights of evidence methods

Groundwater resource management in the urban area is one of the important aspects because of growing population demand and having inadequate water supply. So, proper information is needed to manage the future urban planning for such kinds of areas. This study emphasizes groundwater potentiality zone (GPZ) assessment in the Asansol urban agglomeration (AUA) region, West Bengal, India. For this purpose, we have incorporated eight conditioning factors namely LULC, lithology, slope, elevation, rainfall, drainage density, lineaments density, and soil map using conventional and remote sensing data in GIS software. All these conditioning factors have been reclassified in ArcGIS and processed by the fuzzy analytical hierarchy process (FAHP), frequency ratio (FR), and Bayesian weights of evidence (BWOE) statistical methods. Then, the groundwater potential index has been formulated, and finally, GPZ maps are generated based on of selected three models. The result shows that very high area of GPZ, e.g. 9.13%, 11.62% and 7.43% are under BWOE, FR and FAHP models, respectively. The receiver operating characteristic validation curves show that FR method (AUC = 96.4%) is well obtained for GPZ in comparison with both BWOE (AUC = 83.8%) and FAHP methods (AUC = 82.9%). Therefore, this statistical method is highly recommended for the study of groundwater potential assessment and this outcome is very suitable for the groundwater resources management in future land use planning. Precautionary works in low potential areas should be given priority for long-term planning. Thus, this study can be considered as a good document for decision support in water exploitation planning and sustainable groundwater management in AUA region.


Introduction
Groundwater is one of the precious entities in nature that occupying the pore and fracture space under the geological stratum (Naghibi et al. 2015). It plays a significant role in terms of environmental, economic and human wellbeing throughout the world (IPCC 2001;Mallick et al. 2019). About 30% of the world freshwater is reserved as groundwater, whereas surface water contributes only 0.3% of freshwater (Senanayake et al. 2016). Since the last century, the demand for freshwater is immensely increasing due to the agricultural revolution, industrialization, and urbanisation (Manap et al. 2014). According to the census of India, in 1991, about 1/10 people lived in urban areas but now (2011) it has become 1/3 (Chandramouli and General 2011). The increasing trends of urbanisation create an impact on the quantity and quality of available groundwater in various ways such as low recharge due to soil sealing, water shortage due to increasing demand, point and a non-point source of water pollution (Rogers 1994;Strohschon et al. 2013;Jha, Singh and Vasta 2008;Wakode et al. 2018;). Consequently, the proper supply of clean and fresh water has become a challenging issue. Near future urban areas have to be self-sustain in water through proper management of water supply and consumption. The most exigent part of groundwater resources management and protection are to map and delineate the groundwater recharge potential (Pathak 2017; Chen et al. 2018). The groundwater potential aquifer has supplied good quality water to the urban areas throughout the year (Singh et al. 2013). So, GPZ delineation is a prerequisite to identify the site suitability for the future urban residential zone . The availability of groundwater in any region is influenced by the interaction among hydrogeological, climatic and biological factors (Das 2019). Hence, geology, slope, topographical and lithological variation, elevation, soil type, drainage, and lineament, rainfall, distance from the river, land use and land cover, etc. identify the groundwater store and movement under the geological stratum (Ozdemir 2011;Rahmati et al. 2014;Tiwari et al. 2019;Oh et al. 2011;Mallick et al. 2014).
The invention of integrated remote sensing (RS) and GIS techniques is now a paradigm shift in groundwater exploration research through accurate access and monitoring. The satellite-based RS techniques have provided freely, accurate and neutral data about the earth surface, which helps to create the attribute data layer for delineating GPZ. GIS-based multi-criteria analysis, geo-information techniques, and integrated multi-criteria decision analysis (MCDA) are very useful methods for decision making for problems involving too many influencing factors in identifying groundwater recharge zone (Machiwal et al. 2011;Kumar et al. 2014;Fenta et al. 2015;Singh et al. 2017). Besides the GIS-based weighted overlay analysis using thematic maps, several statistical approaches like analytic hierarchy process (AHP) (Shekhar and Pandey 2015; Rahmati et al. 2015), frequency ratio (FR) (Oh et al. 2011;Elmahdy and Mohamed 2015;Zeinivand and Nejad 2018), fuzzy analytic hierarchy process (FAHP) (Aryafar et al. 2013;Sener et al. 2018;, certainty factor (Razandi et al. 2015) and Bayesian weights of evidence (BWOE) (Al-Abadi 2015; Chen et al. 2018;Kordestani et al. 2019) are also used to delineate the GPZ for better representation. These methods can be validated through ROC curves .
The Asansol urban agglomeration (AUA) becomes an economic and industrial hub and the second largest urban agglomeration of West Bengal (Maity et al. 2020). Consequently, the population pressure and urban built-up area are increasing rapidly in these cities (Shikary and Rudra 2020;Maity et al. 2021). As this is a mining and industrial region, water supply and management will be a challenging issue in future. Therefore, the aim of this research is to delineate the groundwater potentiality zone for sustainable urban planning. In this regard, this study ensemble fuzzy AHP, FR, and BWOE models to get better results than single models, while delineating groundwater potentiality zones based on eight thematic layers, e.g. DEM, slope, drainage density, lineament density, soil, lithology, rainfall, and LULC, that could have served water resource planners, decision-makers, and urban planners for proper urban planning and sustainable use of water resources.

Study area
Asansol region is one of the significant urban bodies in West Bengal in terms of economic activities and urbanisation. Consequently, the landscape pattern has been changed within the last few decades. Recently, this study area took the attention of the urban geographer. Asansol urban agglomeration, which lies in 23°34'30" N to 23°48'00'' N and 86°8'00" E to 87°8'00'' E ( Fig. 1), is consisting of 3 municipalities and 1 Municipal corporation covering an area of about 461.20 km 2 , and 11.5 lakhs population. Physio-graphically, this area is located in the lower part of Chotanagpur Plateau on the left bank of Damodar River with an average elevation of 111 m from mean sea level and it consists of the meta-sedimentary rocks of Precambrian age, Gondwana sedimentary rocks, Rajmahal basalts, and upper tertiary sediments. This area has coarse gritty soil blended with rock fragments that are formed from the weathering of pegmatite, quartz veins, and conglomeratic sandstones, whereas sandy, Red, and Yellow Ultisols soil with low in nitrogen, calcium, phosphate, and mineral resource such as coal, iron ores, calcium carbonate, abrasives, silica bricks, glass sands, building materials, moulding sands, manganese, bauxite, laterite, etc. are found in this region. A very mixed types of natural vegetation are found in this region like Sal (Shorea robusta), Palas (Butea monosperma), Mohua (Madhuca longifolia), Kendua (Diospyros melanoxylon), Shirisha (Albizia lebbek), Bans (Bambusa arundinacea), Arka (Calotropis gigantea), Arjun (Terminalia arjuna) and Ashan (Tilia tomentosa). The climate of this region is tropical savanna types. The April to June months are the warmest month and December to February months are the coolest month with a mean annual temperature is 26.2 °C and mean annual rainfall is 1430 mm.

Database
To examine the groundwater potentiality zone (GPZ), various types of data are used and acquired from different sources. Rainfall data have been collected from the Indian Meteorological substation in Asansol. The cloud-free satellite images from multi-spectral Landsat -8 (OLI/TIRS) and Digital Elevation Model (DEM) have been acquired from the United States Geological Survey (USGS) (http:// earth explo re. usgs. gov). All the detailed specifications of the satellite images are shown in Table 1. The soil samples collection and the location of existing wells are identified from different sites of the study area during field surveys using GPS. The laboratory experiment of these soil samples has been done for grain size distribution and permeability. The data processing of the whole study was done by ArcGIS, ENVI, MS-office, and SPSS software. All the methods used in this study have been indicated in Fig. 2.

Data processing for groundwater potential zones mapping using GIS
In order to evaluate the groundwater potentiality zone of Asansol Urban Area (AUA), eight thematic maps, namely LULC, lithology, slope, elevation, rainfall, drainage density, lineaments density, and soil characteristics, have been generated using conventional and remote sensing data in GIS software.

Land use/land cover (LU/LC)
Globally, anthropogenic activities (directly or indirectly) are responsible for damage to the environment such as groundwater depletion, deforestation, soil erosion, and loss of soil quality. Land use/land cover change is one of the crucial factors for the change of groundwater. In this study, supervised land use/land cover classification has been done by using the Gaussian Maximum Likelihood Classifier Algorithm (GMLCA) in ArcGIS (v.10.1) software. The classification result has been validated through Kappa coefficient. As a result, the six different LULC classes ( Fig. 4.h) have been found and the proportionate share of each class is about agricultural land (36.85%), built-up area (21.20%), dry/barren land (28.20%), mining and industrial area (4.91%), vegetation cover (7.20%), and waterbodies (1.65%). Consequently, with the increasing trend of urban areas, groundwater availability becomes a question to the Asansol urban area.

DEM
The groundwater prospects are likely to be affected by topographic elevation and regulated by various hydrogeological and geomorphological processes. The topographical elevation in this study has been developed by the digital elevation model (DEM) (Fig. 3a) and extracted from ASTER DEM data. The higher elevation is found in the north-western part and the lower elevation is found in the south and southeastern part of AUA along the Damodar Riverbank. For

Slope
The availability and flow of groundwater are highly controlled by the slope (Yeh et al. 2016). The slope was derived from ASTER DEM data of AUA (Fig. 3b). In this study, the slope has been reclassified into five classes according to the value of slope (Boughariou et al. 2021) and represented by percentage value. The slope area of 0-6% is covering almost entire area, which refers to the gentle slope in the favour of highest rate of infiltration and 6-12%, 12-18%, 18-24% and > 24% covering the remaining part of the study area. The infiltration rate of surface water into subsurface water is high in gentle slope area which allows more time to percolate and vice versa. (Das et al. 2019).

Lineaments Density
The lineaments are the area of weakness surface with some linear to curvilinear features in the geological structure, such as fracture, fault, and joint. The groundwater intensity is highly influenced by the lineament density (Al-Ruzouq et al. 2019;Yeh et al. 2016). The geological map was used to extract the lineament density in this study. The high weight assigned for high lineament density represents the high recharge zone and vice versa. In this study, five lineament density zones (Fig. 3d) were identified which are 0-0.13 (very low), 0.13-0.35 (low), 0.35-0.58 (moderate), 0.58-0.85 (high), 0.85-1.43 km/km 2 (very high).

Drainage density
The drainage density is the sum of all stream lengths in a drainage basin divided by the total area of the basin. The structural analysis of a drainage network helps to identify the groundwater recharge zone, and the quality of drainage system relies on a significant index on percolation rate (Yeh et al. 2016). The drainage density is significantly correlated with the groundwater recharge, higher the drainage density refers to the high level of groundwater recharge (Yeh et al. 2016). The drainage density map (Fig. 3c) is extracted from the ASTER DEM in ArcGIS. In this study, higher the weight attributed to high drainage density and vice versa.

Rainfall
Rainfall is one of the essential weather variables in respect to the delineation of groundwater potentiality zone. Therefore, it is necessary to understand the spatio-temporal phenomena of rains for the study of water resources. The rainfall data were acquired from the meteorological substation in Asansol Municipality. The average annual rainfall of AUA is around 1200-1400 mm and 80-85% of the total rainfall happens during July-September due to south-west Monsoon. The eastern and western part of AUA receives 1200-1400 mm/year rainfall ( Fig. 3g) and central part of AUA receives > 1400 mm/year rainfall. This little difference of rainfall around the whole study area has no major influence on the variability of the groundwater potentiality zone. However, rainfall is still an important controlling factor, because rainfall is input on surface water storage and created a balance among all components of hydrological cycle.

Soil
The availability of groundwater is highly influenced by the infiltration capacity of the topsoil of any region. Soil porosity and permeability are directly controlled by the particle size of the soil. Therefore, the characteristics of soil is an essential controlling factor for delineating the GPZ. The soil map (Fig. 3e) of AUA is collected from the National Bureau of Soil Survey and Land Use Planning. The different types of soil are available in this study area, e.g., loamy soil covers 0.81% of the total area, sandy loam covers 27.16%, sandy loam-gravelly sandy loam 0.16%, and sandy clay loamy covers 71.87% of the total area. The porosity and permeability are higher in sandy and gravelly sandy loam soil and much lower in case of clay soil. The characteristics of the soil found in this study area are more or less similar. However, the weight assigned for GPZ is higher for courser grain and less weight is assigned for finer grain soil.

Lithology
Lithological rock types exposed to the surface has a significant impact on groundwater recharge. In this study, the lithological map (Fig. 3f) has been collected from Geological Survey of India (GSI). This map illustrates five classes of lithological units, e.g. flood plain of Damodar River, toe slope, alluvium plain, gently sloping land, moderately sloping land. The flood plain of Damodar River and alluvium plain has the highest water holding capacity as compared to other lithological units.

Fuzzy analytic hierarchy process (FAHP)
The analytical hierarchy process (AHP) is very useful for the multi-parametric evaluation (Saaty 1980) ( Table 2). FAHP is the analytical method which is very useful to encounter the unstructured issues regarding geomorphological, socioeconomic, and environment (Özdağoğlu and Özdağoğlu 2007). Moreover, fuzzy AHP is one of the hierarchical structures which envisages pairwise comparison of thematic layers with triangle fuzzy number (TFN) and it provides more accurate results compared to the decision-making traditional method AHP (Boughariou et al. 2021).

Frequency ratio (FR)
Frequency ratio (FR) model is one of the bivariate statistical method that helps to measure the probability between dependent variables and independent variables of the groundwater potential zones (Balamurugan et al. 2017;Boughariou et al. 2021). In this study, groundwater conditioning factors are considered as independent variables and bore well data are considered as dependent variable. The FR can be calculated as follows: where W denotes the number of pixels of the bore well locations for each conditioning factor; BW denotes the number of total bore well pixels in study area; P represents the number of pixels of each class of the conditioning factor; TP denotes the number of total pixels in the study area. In FR calculation, the FR model value is obtained from the thematic layers of the conditioning factors that determine the groundwater potentiality (Balamurugan et al. 2017).

Bayesian weights of evidence (BWOE)
Bayesian weights of evidence (BWOE) model is a widely accepted method to calculate the groundwater potentiality (Lee et al. 2012;Kordestani et al. 2019;Boughariou et al. 2021). Based on relative weightage of the conditioning factors, BWOE model is calculated in a study area. BWOE can be expressed as follows: where W + and W − denote the positive and negative weights of the conditioning factors; P is the probability, B stands for conditioning factor, B ′ stands for absence of conditioning factors, A stands for bore well data and A ′ stands for absence of bore well data.
Then, the standard deviation (SD) of the contrast value can be calculated as follows: where S 2 denotes the influencing weight of the conditioning factors (Zeinivand and Nejad, 2018).
The standardised contrast is considered as measure of confidence can be calculated as follows: After assigning all the conditioning parameters with the bore well data, the value of W + , W − , SD c , and has been calculated.

Delineation of groundwater using FAHP
All the thematic layers are integrated with ArcGIS to generate groundwater potential index (GPI), to evaluate the groundwater potentiality zone (GPZ). The GPI has been computed by using the weighted linear combination method (Malczewski 1999;Mallick and Rudra 2021a, b) as follows: where wi is the normalized weight of the ith thematic layer, xj is the rank value of each class with respect to the j layer, m is the total number of themes, and n is the total number of classes in a theme. The GPZ for each factor is calculated using FAHP Eq. 12: where w and r represent the normalized weight index and rating of FAHP of the individual classes, respectively. E is the DEM, SL denotes slope, S represents the distribution of soil, DD represents drainage density, LULC is denoted land use land cover of the study area, P indicates volume of rainfall, LD is the lineament density, and L represents lithological condition of the study area.

Delineation of groundwater using FR
Contrasting the FAHP, in FR, the weightage of each class is not put based on the properties of the conditioning factor but given in the form of spatial occurrence of the wells in each class. Similarly, the FR is calculated for all the conditioning factors. Finally, the GPI has been computed using Eq. 13: where FR is the final weight of the conditioning factors. Then, GPZ is calculated based on frequency ratio-based GPI value of each conditioning factors as follows (Eq. 14): where FR represents the weights of frequency ratio of the individual classes. E is the DEM, SL denotes slope, S represents the distribution of soil, DD represents drainage density, LULC denotes land use land cover of the study area, P indicates volume of rainfall, LD is the lineament density and L represents lithological condition of the study area.

Delineation of groundwater using BWOE
BWOE is similar like FR; it is calculated and mapped according to the τ values using Eq. 15: where τ is the final weight that comes through the calculation of W+, W−, C, and SD c . Then, the GPZ is calculated according to this GPI value as follows (Eq. 16): where E is the DEM, SL denotes slope, S represents the distribution of soil, DD represents drainage density, LULC denotes land use land cover of the study area, P indicates volume of rainfall, LD is the lineament density and L represents the lithological condition of the study area.

FAHP
The thematic layers of conditioning factors were classified and rank was attributed to each class using Saaty's scale (Table 3) to establish the Fuzzy AHP-based groundwater potentiality map. The rank for DEM classes is decreasing following the numerical values of the thematic layer. For the lowest altitude (43-93 m.), the rank is 5; for the intermediate altitude, classes are given the rank 4, 3, and 2, respectively, while the highest altitude (137-209 m.) class has the lowest rank 1. The rating of slope classes also has decreasing rank while the value of slope is increasing, and the highest rank (5) is attributed to the lowest slope value (0-6), whereas the lowest rank (1) is given to the highest slope value (> 24). The presence of lineaments influences the infiltration rate. Therefore, highest rank (5) is given to highest class (0.85-1.43) and rank 1 is provided for the lowest class (0-0.13). The characteristics of soil also regulate the infiltration rate. Here, highest (5) and lowest rank (2) are assigned to loamy and sandy loam-gravelly sandy loam soil. The rating for lithological unit is provided based on the conductivity. For lithological class, highest rank (5) is given to Flood plain of Damodar River and lowest rank (2) is assigned for moderately sloping land. The amount of rainfall has an important role in groundwater recharge. In this study, two rainfall zones have been found and highest rank is given to 1400 mm/year and the lowest rank for 1200 mm/ year. The rate for the LULC class has been assigned based on the infiltration capacity of each class. The highest rank (5) is given for water bodies, rank 4 for vegetation cover and agricultural field, rank 3 for dry/barren land, and rank 2 is assigned for built-up and industrial areas. The FAHP is used to compare the eight thematic layers and normalised weights have been calculated (Table 3). As a result, the highest weight is attributed to lineament density (0.18), and the lowest weight is 0.08 for DEM parameter. The GPZ map has been generated by overlaying a thematic map based on the FAHP weights. The GPZ map (Fig. 4) obtained using FAHP model is specified four groundwater potential zones, like very high zone, covering 7.43% to the total area of AUA and found along Damodar Riverbank, a high zone covering 34.60% area, moderate zone (37.08% area) and low potential zone covering (20.90% area) ( Table 5).

FR
The groundwater potentiality zone estimation using FR model is very significant. This model executed the GPZ through correlation between conditioning factors and location of bore wells. Moreover, higher correlation value indicates the greater groundwater potential and vice versa (Manap et al. 2014). In this study, eight conditioning factors (Fig. 3) and twelve bore wells (WRIS) have been used to create GPZ map. The correlation results (Table 4) of FR model depict that lower slope per cent (0-6) is indicated by the higher value of FR (2.253), while the absent of wells in higher the slope per cent (> 24) area is indicated lower FR (0.00) reflecting the low groundwater potential. Meanwhile, higher value of lineaments density (0.85-1.43) shows high FR (9.328) value which indicates the greater chance to groundwater potential while low FR value (0.522) was found for the lineaments class 0.35-0.58. For lithological conditions, the result shows that the alluvium plain has the greater ability to groundwater potential with FR value (1.951). The soil characteristics are also significant for delineating GPZ. According to FR model, the sandy loam soil (FR = 4.149) has more groundwater potential than other soil found in this region. The amount of annual rainfall is also very significant for groundwater recharge. Here, two different rainfall zones have been found, i.e. > 1400 mm/year with FR value 1.457 and 1200-1400 mm/year with FR value 0.389. Finally, the impact of urbanisation on groundwater potentiality is highly influenced by LULC classes. In this study, water body contributed the high potential ability with 9.923 while mining/industrial area and vegetation cover area found insignificant FR values because FR is analysed by conditioning factor and bore wells data. Due to lack of bore wells data of these classes, the result indicated the low FR value while vegetation cover always greatly influenced the infiltration rate. The GPZ map (Fig. 5) obtained using FR model is specified four groundwater potential zones like very high zone, covering 11.62% to the total area of AUA and found along Damodar river bank and south-eastern part; high zone (32.45% area); moderate zone (33.84% area) and low potential zone (22.08% area) which found in the built-up and industrial areas (Table 5).

BWOE
The groundwater potentiality zone estimation using BWOE model (Fig. 6) is also useful for this study. According to BWOE (τ) values (Table 4), it is noted that for the lineament, the highest τ value is given which indicates a high groundwater potentiality. Then, rainfall, drainage density, lithology thematic maps are given next level of τ value. Similarly, the LULC (specifically vegetation and agricultural land) has been given a positive standard value. For DEM, highest weight of τ is given to the lowest altitude in the Damodar riverbank area. The lowest weight of DD has also a higher value of τ, while the highest value of this parameter is detected in the highest class of lineament. For BWOE model, the quintile method was also adopted to classify the groundwater potential map into three classes: low, medium, and high. According to the final GPZ map, the low GPZ is found mostly on the Damodar riverbank area. The delineation of the BWOE-based GPZ map shows a similar surface cover for the high and moderate groundwater potential. The low GPZ class covered slightly higher area (38.33%) than the high and moderate potential zones ( Table 5).

Validation of GPZ
The FAHP, FR, and BWOE models have been validated through receiver operating characteristics (ROC) curve. The area under the ROC curve (AUC) defines the prediction or classification accuracy (Boughariou et al. 2021). In this study, we have run and validated three models based on the classification result of the GPZ. AUC value is ranged between 0 and 1. If the value is given below 0.5, then it indicates that the result of the model is not suitable, and it needs to classify once again and close to 1 implies that the result is well delineated ). The ROC curves of the GPZ maps were generated for the validation of FAHP, FR, and BWOE methods (Fig. 7). The result demonstrating that the outcome of FR method (AUC = 96.4%) is well obtained for GPZ in comparison to the both BWOE (AUC = 83.8%) and FAHP method (AUC = 82.9%). However, all the obtained results were validated and well delineated . But, it can be stated that the FR is better representative for this study area to denote that the spatial distribution of the GPZ with respect to the fuzzy AHP and BWOE methods. Therefore, this  validation method is highly recommended for the study of the groundwater potential assessment.

Discussion
In this study, three models, e.g. FAHP, FR, and BWOE have been ensemble for groundwater potential zone assessment for getting better result (Nair et al. 2019;Werner et al. 2019). Eight conditioning factors of thematic maps were used as indispensable input parameters to calculate the GPZ. The GPZ value of the three models was calculated based on groundwater potential index values. The FAHP method was developed based on the literature and experts' knowledge which influenced the weightage of the conditioning factors on the study area . Moreover, the FR and BWOE have imposed the distribution of wells that helps to estimate the weights of each class of the conditioning factors. The output maps were classified into four zones, i.e. low potential, moderate potential, high potential, and very high potential. The results of FAHP model demonstrate that very high and low groundwater potentiality zone covers 7.43% and 20.90% to the total area, respectively. However, the maps using FR and BWOE methods revealed better representation where 11.62% and 9.13% of high groundwater potentiality including 38.88% and 22.08% of low groundwater potentiality, respectively. That means this area has high water stress than availability of potential groundwater. For these three adopted methods, the high groundwater potentiality area was located on the riverine area on the southern part of the AUA and low potentiality zones were found in the urbanised areas of the AUA. Moreover, similar method was found for groundwater potentiality validation that was area under the ROC curve (Rahmati et al. 2016).
However, it is well known that the weights of conditioning factor in groundwater potential mapping are significantly influenced by the LULC characteristics and the lithological characterises (Naghibi et al., 2015). The findings of this study support the results of Lee et al (2012) that high lineament density locations are advantageous for groundwater potential. According to Kumar et al. (2014) and Fenta et al. (2015), since lithological character and lineament density facilitate percolation and hence improve groundwater recharge. The groundwater potential increases from 1200 to 1400 mm in case of rainfall and then steadily decreases from 1400 mm onward. This could be attributed to a decrease in impermeable surface infiltration in the upland area, as well as geomorphological and topographical factors that lead to excessive runoff (Pradhan, 2009). Nevertheless, analysis of GPZ result is also helpful to identify the future urban land use development area and to maintain the urban resilience (Mallick 2021).

Conclusion and policy recommendation
Asansol urban area is a rapidly rising urban agglomeration in West Bengal, India. The WWF-India reported that the 30 cities in India will face a grave water risk by 2050 due to sharp increase in population. The Dhanbad city in Bihar is in this list, which is located only 63.5 km away from Asansol region. In future, AUA might be facing the same crisis because both cities have the same physiographic characteristics. Therefore, proper urban built-up area planning is prerequisite for sustainable urban development. This study emphasised the identification of groundwater potentiality zone for predicting built-up area development sites around the AUA. The results of this study will help the AUA and Asansol municipality for sustainable future water resource management, including regional land use planning, future wells construction, and groundwater protection. Although the shallow aquifer of AUA was overexploited, particularly near the Damodar River area, this outcome helps them to take further necessary action in the future days. Thus, this study can be considered as a good document for decision support in water exploitation planning and water management in AUA region. However, a few steps can be taken regarding water supply and increase in groundwater recharge like supply of water to the urban dwellers from Damodar River and Maithon water reservoir after purification, implementation of rainwater harvesting projects in the urban households, etc. Additionally, there need some more inclusive innovative methods that can be demonstrated regarding better management of water resources from the over-exploitation and changing climate scenarios.
Funding The present research has been funded by the University Grant Commission (UGC), Govt. of India, and it is received by authors Maity, B., Mallick, S.K., and Das, P.