Delineating springs safeguard zones with DISCO-URBAN index: a valuable tool for groundwater vulnerability mapping in local-scale urban areas

Assessing groundwater vulnerability to contamination is vital worldwide, particularly in sustainable water resources management. That is mainly a concern in fractured media in urban areas due to a large diversity of contaminant sources and the complexity of recharge pathways. Thus, groundwater vulnerability assessment is essential to delineate groundwater protection zones around springs or wellheads. Furthermore, it considers the groundwater system’s heterogeneity and the surrounding hydrogeological conditions, as well as provides suitable solutions to protect the resource and mitigate potential hazards. DISCO-URBAN index focused on urban areas was applied to evaluate the intrinsic vulnerability in fractured media in the surroundings of Penafiel city (NW Portugal). The analytic hierarchy process (AHP) was used to determine the weight of each evaluation factor. Furthermore, multi-criteria indexes were applied: GOD‐S, DRASTIC‐Fm, SI DISCO and DISCO-URBAN. Low to moderate vulnerability classes dominate in the combined approach of the vulnerability indexes. However, very high vulnerability classes occur in DISCO and DISCO-URBAN, corroborated by the water-enriched nitrates. Therefore, the DISCO-URBAN method highlights a better delineation of groundwater safeguard zones. In fact, the DISCO-URBAN index is reliable in urban areas to be integrated as a tool to develop local site hydrogeological investigations related to springs safeguard zones.

Location of Santa Marta springs and underground water mine, Penafiel urban area (NW Portugal): a a general background and hydrogeological inventory; b local geological setting; c land use framework; d climatology features Table 1 Synthesis of the hydrogeological features of the groundwater systems in Santa Marta springs surroundings (regional and local data updated from [15] na not applicable; X outcrop Hydrogeological   Figure 2 presents the Santa Marta water mine setting, which can be crucial for water supply to the main green spaces and public fountains in the urban area of Penafiel (Fig. 2a). The area's groundwater flow is related to the Santa Marta springs [15,49]. The Santa Marta springs are located at about 380 m, close to the water deposit of the Municipality of Penafiel, namely in a hillslope "Monte do Fogo" site. The underground galleries mainly trend an E-W orientation, then inflect to NNE-SSW. Historical and current data related to the springer flow was collected to assess the potential of groundwater sustainability (Fig. 2a, b). The measured flow in the summer at the Santa Marta mine is close to the recorded historical values, i.e., about 0.40 L/s (details in [15]).
The regional geomorphology and geotectonic framework is part of the Central-Iberian Zone (ZCI) of the Iberian Massif [56]. Locally, the study area is located in dominant tectonic lineaments towards NE-SW, mainly in the flattened area (Santa Marta) (Fig. 2c). That will be responsible for the regional arrangement of the relief in this region, where the entire urban area of Penafiel develops.
Moreover, the drainage network refers to this tectonic control, which gives the area morphostructural features. The mine's entrance is at the height of 352 m (Fig. 2c) and should be found, entirely excavated in granitic rocks, namely in the porphyritic granite, coarse-grained, slightly to moderate weathered (W 2 -W 3 ) (Fig. 2e). The principal fracture network highlights the following directions: NE-SW, NNW-SSE, and ENE-WSW. In addition, at the local scale, outcrops quartz veins and doleritic rocks, with N50°E orientation (at the surface, dolerites are mainly found highly weathered [W [4][5], with orangey colour).

Materials and methods
To achieve a comprehensive analysis and conceptualisation of the Santa Marta springs (Penafiel urban area, NW Portugal), a multidisciplinary approach using field and laboratory techniques were performed to access urban groundwater vulnerability to contamination and delineate safeguard zones in urban areas (Fig. 3). Using a GIS-based mapping, a collection of several urban data was successfully organised, namely geology, hydrogeology, hydrogeomorphology, land use/cover, and urban hydraulics and sanitation. Also, GIS overlay techniques were used to integrate and analyse the urban geodatabases and produce all thematic maps [50,52] and geovisualisation techniques [57]. As a result, the mapmaking procedures and design and the conceptualisation of maps have developed immensely. Furthermore, new trends exploring geovisualisation analysis integrate approaches from different disciplines, including scientific visualisation, image analysis, information visualisation, exploratory data analysis and GIS science (e.g., [50,57]). According to Kraak [58], geovisualisation combines the strength of the computer (automated analysis techniques and geo-computation) and the human (interactive visualisations for adequate understanding, reasoning and decision making).
Portugal land use/cover map (COS 2018) was used to define the land use and cover of the urban study area. Furthermore, to improve the high-resolution mapping and maintain the base structure of the COS2018, it was cross-checked and updated using Maxar satellite imagery (0.5 m resolution, available through the ArcGIS base map service) and fieldwork, following the guidelines suggested by [59]. Additionally, the inventory of potential contamination activities was mapped from extensive PhD research surveys during 2014-2021 [60,61]. The geology map background was revised from the Portuguese Geological Survey and was updated during the fieldwork. The hydrogeotechnical surveys have followed the recommendations of the [62].
Moreover, a multi-scale urban groundwater vulnerability to contamination assessment in the Santa Marta water springs was delineated based on comparison and integration of three indexes to reach more accurate results, namely GOD-S [27,41] DRASTIC-Fm [42], SI [46]. The previous indexes are significant for regional to local purposes. The data was processed, and the indexes were calculated using raster format and overlay tools in the GIS environment. The information that served as a basis for the vulnerability assessment corresponds to the hydrogeological background. The overlay analysis was made by converting each layer into a raster of 4 × 4 m (16 m 2 ). The international colour code for DRASTIC was applied to the vulnerability representation.
DISCO index [30,43,63] and DISCO-URBAN index [37] were calculated to delineate groundwater protection zones mapping around the Santa Marta springs and underground water mine site. The DISCO method has successfully proven its ability to assess groundwater's vulnerability to contamination, especially in highly heterogeneous fractured aquifers, by its dependence on DIScontinuities and protective COver parameters (e.g., [30,43]). Conversely, the DISCO-URBAN index is a new approach developed for urban areas with an excellent asset for local-scale investigations [37]. Table 2 shows the strengths and weaknesses of the groundwater vulnerability indexes in fissured rock media focused on urban areas.

Step 1. Assessment of the discontinuities (D) and protective cover (P) parameters
This step aims to assess and map the parameter's discontinuities and a protective cover over the entire catchment region. The discontinuities parameter (D) is based on the tectonic lineaments analysis. It considers the groundwater flow velocity within the fractured aquifer between an infiltration point in the water catchment region and the springs under
The maps use diverse legends and explanatory notes to promote a uniform mapping methodology (e.g., [38]) Local-scale: Large-scale maps to high-resolution vulnerability maps for delineating safeguard zones (local framework: 1:1000-1:10,000; detailed surveys: ranging 1:100-1:250). The maps use diverse legends and explanatory notes to promote a uniform mapping methodology (e.g., [38]) Water sources Focused on a large number of water occurrences Focused on delineation safeguard zones for singular water occurrence (e.g., spring) Watershed features Do not consider the watershed limit

Considers the watershed limits
Urban features Do not take into account the land use and impervious areas Take in account the land use/cover Do not take into account the land use and impervious areas Take into account the land use/ cover and impervious areas Do not take into account the hydraulic and sanitary features Hydraulic and sanitary features: water supply, stormwater, sewer consideration. The "D" rating values are in the range 0-3, with increasing values corresponding to higher residence time and attenuation processes. The protective cover parameter (P) is based on soil analysis. It considers the protective effect of water flow directly through the soil and geological formations overlying the fractured aquifer. The rating values of "P" are in the range 0-3, with increasing values corresponding to each of the higher protective cover thickness and lower permeability of the deposits.

Step 2. Determination of the intermediate protection factor (F int )
The protection factor (F int ) map is estimated by combining all the parameter maps, and it ranges from 1 to 9 related to the vulnerability, whereas an extremely low protection factor corresponds to a very high vulnerability. To calculate the intermediate protection factor (F int ), the following formula was applied: F int = 2 × D + P.

Step 3. Protection zone delineation (F)
According to the DISCO methodology [43,63] and considering the topography's slope and the runoff area, achieving the final protection factor "F" is possible. Considering the watershed limits is crucial since the runoff may transport some contaminant load for the spring area. This contaminant load will be most damaging if it reaches areas of high vulnerability. Therefore, the final protection factor (F) map is converted into safeguard zones, using the inverse relationship between the value of the final protection factor "F" and safeguard zones "S".
Meerkhan et al. [37] proposed that the DISCO-URBAN Index considers four main parameters. Consequently, this Index contemplates the same two parameters of the DISCO methodology: DIScontinuities "D" and protective Cover "P", which are used in the original DISCO index and they take the same classes, and also two new parameters: Land Use/Cover parameter "LUC" and Urban hydraulics and sanitary features "HS" including water supply, stormwater and sewer networks. The general methodology and tables with weights assigned to each parameter could be found in the work of Meerkhan et al. [37]. First, GIS-based mapping was used to achieve the urban hydraulic and sanitary features by applying a line density tool considering different weights for each network: 1 × water supply; 2 × stormwater, and 4 × sewer network. LUC ranged from 0 (high to very high) to 3 (very low to low), and HS is also in the range of 0 (very high-lower protection) to 4 (very low to null-higher protection). Next, the relative weight/score for each parameter used in the DISCO-URBAN was calculated using the analytical hierarchy process (AHP). The AHP multi-criteria technique compares all parameters, assigning a weighting to each parameter according to its importance (e.g., [64][65][66][67]). Next, a proper Consistency Ratio (CR) is calculated by dividing the Consistency Index for the set of judgments by the Index for the corresponding random matrix [65]. Saaty [64] suggests that if that ratio exceeds 0.1, the set of judgments may be too inconsistent to be reliable. In practice, when CR < 0.1, the consistency of the judgment matrix is acceptable; when CR ≥ 0.1, the judgment matrix needs to be tested again (details in [65][66][67][68]).
Therefore, the understanding collected in the hydrogeological fieldwork was essential for the hierarchy of parameters. Then, the following formula was applied to calculate the urban intermediate protection factor (F int-Urb ): F int-Urb = 2 × D + P + 2 × LUC + HS. Finally, the safeguard zones were computed, applying the inverse relation between the urban protection factor (F) and the safeguard zones (S). Once again, GIS technology was used to calculate these indexes and vector operations were performed using the intercept function to overlay the vector data (details of the methodology in [37]. Figure 3 presents a general methodological flowchart with the main steps of delineating groundwater vulnerability and safeguard zones mapping in urban areas and applying the DISCO and the DISCO-URBAN indexes to be effective on a local scale (urban areas). Finally, the integrative approach improves the urban hydrogeological conceptual site model in the Santa Marta springs, improving water resources management and developing environmental sustainability. legislation defining strategies for determining protection zones around urban groundwater recharge areas (e.g., Decree-Law 382/99; Decree-Law 58/2005; Ministerial Order P.702/2009). Therefore, several indexes were applied using GIS overlay techniques. By integrating and comparing these indexes' results, it will be possible to analyse the vulnerability accurately to groundwater contamination in the Santa Marta water mine (Penafiel urban area).

Intrinsic vulnerability indexes: GOD-S, DRASTIC-Fm and SI
The Santa Marta water mine area's intrinsic vulnerability assessment was evaluated based on GOD-S, SI, and DRASTIC-Fm indexes. Table 3 presents the synthesis of the parameters' description and the classification adopted for each method. In addition, in the supplementary materials, a synthesis of the parameters' description is presented, along with the classification adopted for the three methods (Tables S1-S3).
The GIS-based urban maps were created to assess the groundwater vulnerability to contamination in the study area. Figure 4 presents the urban groundwater vulnerability maps based on the study area's GOD-S, SI, and DRASTIC-Fm indexes. In addition, an inventory of surface potential contamination activities was surveyed over 12 potential sources (details in Table S4 in supplementary material). The potential contamination sources are located in the vicinities of the Santa Marta water mine.
The groundwater vulnerability assessment based on the GOD-S Index illustrates that most of the Santa Marta water mine area fits in a low vulnerability category. Besides the moderate vulnerability category, which is related to the alluvial deposits (Fig. 4a). The DRASTIC-Fm index, which is the updated version of the DRASTIC method, taking into account the fractured media parameter, shows that the most significant space of the study area fits in a low vulnerability category. The small granitic flattened area fits low to moderate vulnerability. There is also a moderate vulnerability category related to the alluvial deposits in the DRASTIC-Fm index (Fig. 4b). Finally, the SI index indicates that most of the study area fits in the low vulnerability category, where porphyritic granite and forested areas are combined. The rest of the area fits in the low to moderate vulnerability category. In addition to the existence of a moderate vulnerability category related to alluvial deposits (Fig. 4c).
In comparing the three methods, it was recognised that the GOD-S, DRASTIC-Fm, and SI indexes indicated that Santa Marta mine water's surrounding area fits in a predominantly low to moderate vulnerability category. However, two areas should be highlighted with moderate to high vulnerability categories, corresponding to sedimentary deposits (SW and SSE of Crasto de Cima site).

The application of DISCO and DISCO-URBAN indexes
DISCO method allows assessing vulnerable springs' intrinsic vulnerability in highly heterogeneous fractured media, especially on a local scale [43]. This study highlights the application and the importance of the DISCO methodology in

Research
Discover Water (2022) 2:2 | https://doi.org/10.1007/s43832-022-00011-z 1 3 spring protection areas. In addition, it aims to revalidate the development of the DISCO-URBAN for local-scale investigations purposes in urban areas [37]. Figure 5 presents the thematic maps of the parameters used to apply the DISCO index. Figure 5a presents the discontinuities parameter in the study area. D0 is related to the highly permeable discontinuities with preferential connection to the spring, which fits the primary tectonic lineaments in the NE-SW orientation. D1 is related to the discontinuities with a relatively connected to the spring. D2 is related to the discontinuities with a relatively slow connection to the spring. Finally, D3 is related to the low permeability zone. According to the distribution of the classes ( Table 4), most of the area (87%) is characterised by low permeability, which means that the discontinuities have a slow connection to the springs. As such, purification processes are significant and occur in most studied areas. It should also be noted that in the proximity of the spring, these processes are limited due to the preferential connection between the discontinuities and the spring (D0). According to the protective cover parameter, the study area is divided into three classes (Fig. 5b). P1 was related to moderate permeability soil, which fits the alluvial deposits unit. P2 is related to moderate to low permeability soils, consists of porphyritic granite and represents most of the study area (97%). Finally, the P3 class was related to low permeability soil, where the crystalline bedrock is interrupted by dolerite dyke. Figure 5c presents the study area's intermediate protection factor (F int ). Most of the study area fits the high intermediate protection factor (8)(9)(10), which can be justified by the slow connection between the discontinuities (D3) and the spring and the occurrence of soils with moderate to low permeability constituted by porphyritic granite (P2). The low intermediate protection factor (2)(3)(4) is related to the main tectonic lineaments (D0 e D1). The moderate intermediate protection factor (5-7) is related to the alluvial deposits and discontinuities with a low connection with the spring. The class with the lowest representation is the very-low intermediate protection factor, occurring due to the overlap between alluvium and discontinuities with preferential connection to the spring. Figure 5d presents the groundwater protection zones or safeguard zones that classify the study area into four classes according to the vulnerability (Table 4). S 1 zone is related to the very high vulnerability, where the main tectonic lineaments coexisted with porphyritic granite characterised for a moderate to low permeability. S 2 zone represents the high vulnerability class, and it is associated with the occurrence of secondary tectonic lineaments. S 3 zone is related to the moderate vulnerability, the most representative class in the study area, where discontinuities have a relatively slow connection to the spring (D2/D3) coexisting with the porphyritic granite a moderate to low permeability. S 4 zone related to the low to very low vulnerability, which fits the alluvial deposits. Although the sedimentary deposits unit has moderate permeability, the discontinuities have no connection with the Santa Marta springs.

The application of the DISCO-URBAN index
Considering the original DISCO methodology and the used parameters, the newly DISCO-URBAN method is helpful with more accurate results in the urban contexts, especially in local-scale investigations. Furthermore, DISCO-URBAN plays an essential role in managing and protecting groundwater resources in urban environments. So, it is essential to calculate the DISCO-URBAN index to determine the parameters affecting urban groundwater vulnerability. Figure 6 and Table 4 present the application of DISCO-URBAN in the Santa Marta water mine. Figure 6a, b correspond to the parameters presented and used to calculate the DISCO index. Figure 6c presents the Land Use/Cover parameter, which ranges between 0 and 3 classes. Class LUC0 is related to the urban continuous, with an impervious surface grade high to very high. Class LUC1 is related to urban discontinuous, in which moderate to high impervious surface is observed. The rest of the study area is characterised by a low to a very low degree of imperviousness, namely those with agricultural (LUC2) and forestry (LUC3) areas. Figure 6d shows the study area's urban hydraulics and sanitary features parameter (HS). It ranged from 0 (very high) to 3 (low) based on the distribution of the area's water supply, stormwater, and sewer networks. According to the HS parameter, most of the area has a moderate (HS2) to low (HS3) density (84.7%). Figure 6e, f present the urban intermediate protection factor (F int-Urb ) and groundwater protection zones or safeguard zones in the Santa Marta water mine. Analysing the F int-Urb (Fig. 6e) distribution, it can be concluded that the D and LUC parameters stand out from the others. Most of the area (93.5%) fits in the moderate to high classes. In fact, it can be observed in these areas that the coexistence between the highly permeable discontinuities with preferential connection to the spring (D0 and D1) and areas with an impervious surface grade vary from very high to low.
The groundwater protection zones or safeguard zones in the Santa Marta springs (Fig. 6e) split the area into 4 classes. The S 1 class is related to the very high vulnerability and is characterised by the highly permeable discontinuities with preferential connection to the Santa Marta springs. The area in the secondary tectonic lineaments is related to the high vulnerability class (S 2 ). The remaining area, which fits in the agricultural and forest area without any hydraulic connection to the Santa Marta spring, is related to the moderate and low to very-low vulnerability classes (S 3 and S 4 ). processes; D2-discontinuities with a relatively slow connection to the spring (residence time of approximately ten days)/significant purification processes; D3-low permeability zone or discontinuities with a slow connection to the spring (residence time of several tens of days)/efficient purification processes P: P1-Moderate permeability soil (silt, loam), thickness > 1 m; P2-Moderate to low permeability soil (silt, loam, clay), thickness 0.5-1 m; P3-Low permeability soil (loam, clay), thickness > 1 m LUC: 0-Urban continuous; 1-Urban discontinuous, commerce areas and roads; 2-Agricultural areas; 3-Forest areas HS: 0-High; 1-Moderate; 2-Low, 3-Very low Comparing DISCO and DISCO-URBAN indexes shows that the DISCO-URBAN index results are more accurate and give more details regarding the groundwater vulnerability and protection zones. The application of the DISCO and DISCO-URBAN indexes allowed us to conclude that the DISCO-URBAN index reflects more appropriately the hydrogeological conditions observed in urban areas. Table 4 highlights the most noticeable difference between DISCO and DISCO-URBAN indexes are in the S3 and S4 safeguard classes. Also, some S1 areas (very high vulnerability) were more accurately classified in DISCO-URBAN since they do not contribute to spring flow. The DISCO-URBAN methodology also allows distinguishing the S3 and S4 classes, incorporating the LUC and HS parameters.
In summary, it can be concluded that the purification process is limited or insignificant near the Santa Marta spring. Near Santa Marta spring, the area is characterised by: highly permeable discontinuities with preferential connection to the Santa Marta spring (D0), soils with moderate to low permeability, mainly constituted by porphyritic granites (P2); the land use is occupied for urban discontinuous with moderate to a high degree of imperviousness and the urban hydraulic and sanitary features have a moderate density (HS2). Given these conditions, the protection factor is low to very low, which indicates that the Santa Marta springs area has a very high (S 1 ) vulnerability.
The EU legislation defines a water enriched nitrate threshold of 50 mgL −1 to define Nitrate Vulnerable Zones [69,70]. Moreover, the regulation defines action programmes to reduce nitrate loss from agricultural uses, human health protection and water resources contamination [70]. Several works pointed out the correlation of the vulnerability indexes with contaminant parameters or land use (e.g., [71][72][73]). In the present study, the preliminary analytical NO3-data pointed out a value range of 65-95 mgL −1 [15]. In addition, there are recorded nitrate values (lower 40 mgL −1 ) in other springs in the surroundings of the study area included in the same hydrogeological unit (porphyritic granite coarse-grained, twomica). Comparing the nitrates values with the regional vulnerability mapping (GOD-S, DRASTIC-Fm, SI), the higher values match the low to moderate classes (Fig. 4). Conversely, in the local vulnerability mapping (DISCO and DISCO-URBAN), these groundwater enriched nitrates values tie in the very high vulnerability classes (Figs. 5 and 6).
Consequently, The DISCO method allows, at a local scale, a consistent evaluation to delineate groundwater vulnerability protection zones of springs and/or wellhead zones. Thus, DISCO or DISCO-URBAN methodologies' main advantage is related to the local-scale of application, i.e., it is not the best methodology for regional approaches (details in [30,37,43]). That demonstrates the applicability of the DISCO-URBAN approach underlining the local-scale to high-resolution vulnerability mapping to delineate safeguard zones (e.g., springs, wellheads) and urban groundwater management activities.
Comparing the several vulnerability indexes shows the benefit of detailed analysis, offering a better insight into intrinsic and specific vulnerability's existence and spatial distribution, as well as contributing to the site hydrogeological conceptual model [47,74]. Saida [75] pointed out several different models to assess groundwater vulnerability, i.e., there is no universal model for groundwater vulnerability assessment. The selection of models with which the groundwater vulnerability will be assessed depends mainly on the hydrogeological settings of the natural system and data availability [75].

Delineating safeguard protection zones of Santa Marta springs
The Santa Marta water mine may constitute one of the primary sources to supply the central gardens and fountains in the urban area of Penafiel [15]. Mapping groundwater vulnerability and protection zones improved the conceptual hydrogeological site model (Fig. 7). The Santa Marta mine's recharge area is related to the flattened areas of Santa Marta, which are mainly occupied by forests. The precipitation values range between 1300 and 1350 mm/year, with infiltration rates of around 6%. According to the DISCO-URBAN index, the Santa Marta springs is characterised by a very high vulnerability. European and Portuguese legislation defines strategies to determine the protection zones around urban springs (Decree-Law 382/99; Decree-Law 58/2005; Ministerial Order P.702/2009). Therefore, DISCO-URBAN has been proved to be a valuable tool for refining the delimitation of these areas (e.g., [37]). In the Santa Marta springs area, two types of aquifer systems were identified that are directly connected to the springs (details in [15]): i. A superficial unconfined unit (5-10 m), corresponding to the alluvial deposits unit, a shallow water table (< 5 m), a moderate transmissivity (10-20 m 2 /d), a moderate permeability (< 4 m/d), and a Na-Cl to Na-NO 3 hydrochemical facies. Thus, the vulnerability to contamination is moderate. ii. An unconfined to confined unit, corresponding to the highly weathered to decomposed two micas granite, coarsegrained, with a thickness ranging from 0 to 20 m, a shallow water table (5-10 m), a low transmissivity (1.5 m 2 /d), a low permeability (< 1 m/d), and a Na-Cl to Na-Cl-SO 4 hydrochemical facies. Wherein the vulnerability to contamination is low to moderate. iii. An unconfined to confined unit, corresponding to the fresh to moderately weathered two micas granite, coarsegrained, with a thickness ranging from 0 to 10 m, a shallow water table (5-10 m), a very low transmissivity (0.5-1.5 m 2 /d), a very low to low permeability (< 0.4 m/d), and a Na-Cl to Na-Cl-SO 4 hydrochemical facies. In which the vulnerability to contamination is low to moderate.

Conclusions
Evaluating groundwater vulnerability to contamination provides visual analytics and quantitative tools to help planners and decision-makers face the increasing pressure of development and residential dependency on groundwater systems. This study performs an integrative assessment of groundwater vulnerability to contamination in a local-scale urban area, namely in the surroundings of Penafiel city (NW Portugal). A multidisciplinary integrative methodology suitable for local-scale urban areas was used by combining local data from geology, hydrogeology, hydrogeomorphology, hydrogeotechnics, land cover/use and urban hydraulics. GIS-based mapping technology allows the development of an accurate assessment of groundwater vulnerability to contamination. Furthermore, several intrinsic vulnerability indexes (e.g., GOD-S, DRASTIC-Fm, SI, DISCO and DISCO-URBAN) were applied within a combined approach to define the vulnerability safeguard zones. Understanding the source protection areas is crucial for the urban management system by supporting sustainable groundwater systems. The combined assessment of the intrinsic vulnerability maps allows us to conclude that the study area has a low to moderate vulnerability. In addition, the DISCO and DISCO-URBAN index application enables us to delimitate the groundwater protection areas around the Santa Marta springs. The comparison between the two indexes shows that DISCO-URBAN is more accurate for urban areas since it uses some urban parameters that significantly impact groundwater quantity and quality. The study also underlines the DISCO-URBAN index as a consistent method for local-scale delineating groundwater protection areas for groundwater systems management in urban areas.