Assessing Groundwater Vulnerability: DRASTIC and DRASTIC-Like Methods: A Review

: Groundwater vulnerability studies are sources of essential information for the management of water resources, aiming at the water quality preservation. Di ﬀ erent methodologies for estimating the groundwater vulnerability, in general, or of the karst aquifer, in particular, are known. Among them, DRASTIC is one of the most popular due to its performance and easy-to-use applicability. In this article, we review DRASTIC and some DRASTIC-like methods introduced by di ﬀ erent scientists, emphasizing their applications, advantages, and drawbacks.


Introduction
In recent decades, water scarcity and its pollution became a major issue all over the world.
Preserving the groundwater quality is very important for assuring the drinking water resources, given that billions of people all over the world do not have access to water or suffer from water scarcity [1].
Since 1968, when Margat [2] introduced the concept of groundwater vulnerability, many definitions were proposed for this concept. For example, Hirata and Bertolo [3] defined the groundwater vulnerability as "the property of a groundwater system that depends on the sensitivity of the material in permitting the degradation of the saturated zone by pollutant substances originating from human activities", while the National Research Council [4] defined this term as "the relative ease with which a contaminant (in this case a pesticide) applied on or near the land surface can migrate to the aquifer of interest under a given set of agronomic management practices, pesticide characteristics, and hydrogeological sensitivity conditions".
The intrinsic vulnerability describes the water vulnerability to different pollutants (independent of their nature) resulted from human activities and is related to the hydrological, geological, and hydrogeological aquifer's characteristics. Given that the aquifers have different reactions to the same contaminant due to their physicochemical characteristics, the specific vulnerability shows the groundwater vulnerability to a pollutant (or a group of pollutants), determined by the pollutant's properties, taking into account the time of impact, its intensity, and the interaction between the intrinsic vulnerability components and the contaminant [5,6].
Adams and Foster [7] emphasized that the aquifer vulnerability depends on the properties of the layers situated above the saturated zone to attenuate the pollutants' effect, by retention or neutralization by chemical reactions.
Gogu and Dassargues [6] divided the approaches of assessing the groundwater vulnerability in three groups, as a function of the groundwater protection. The first group takes into account only the soil and unsaturated zone, the second one takes into consideration the groundwater flow and the contaminant transfer to some extent [8], whereas the third focuses on the soil, the unsaturated medium, and the aquifer. Different approaches are used for estimating groundwater vulnerability. They can be grouped into three categories. The first group is formed by the index-based methods, which take into consideration only the characteristics of soil and unsaturated zone. They are divided into Hydrogeological Complex and Settings methods (HCS) [9]; Matrix Systems [10], approaches based on the combination of two parameters, and Rating Systems [11][12][13]. They work by building water vulnerability groups using different ratings associated with the physical characteristics of the study media. The second group contains the statistical approaches that assess the groundwater vulnerability through statistical analysis or regression models [14][15][16]. The third one contains the methods based on simulation, which uses simulation techniques for forecasting the processes related to contaminant transport [17][18][19][20]. The index-based techniques have the advantage that they do not depend on data availability or similarities [21].
The procedures that belong to the first and second categories are used for studying the intrinsic vulnerability of large areas [22].
Introduced in 1985, DRASTIC is among the most popular approaches used in groundwater vulnerability estimation due to its capability and easy-to-use. In the following, we shall focus on reviewing this method, and some of the DRASTIC-like procedures that aim to improve the performance of the groundwater vulnerability estimation, emphasizing the differences between them. We shall not focus on the methods assessing the groundwater vulnerability for the karst aquifer because of the extensive literature for the general case and the lack of space. Some classifications of the methods that will be presented in next sections are: 1.
Based on the extent of their use: a.
With general applicability-DRASTIC, GOD b.
We shall indicate the references to the articles treating these methods in the next sections, together with a description of approaches.
The methods (and corresponding parameters) for groundwater vulnerability assessment discussed in this article are summarized in Table 1.

DRASTIC
DRASTIC is a model that considers the main hydrological and geological factors with a potential impact on aquifer pollution. Its acronym stands for D-depth to groundwater, R-recharge rate, A-aquifer, S-soil, T-topography, I-vadose zone's impact, and C-aquifer's hydraulic conductivity [38]. The depth to water table (D) [m] is the thickness of the layer crossed by the pollutant before reaching the aquifer. The aquifer vulnerability is inverse proportional to the depth to the water table.
The land use parameter characterize the human activity as effect on the runoff coefficient, not as the contaminants' nature. ## Refers to the impact of the human activity as impact of the built environment or the nature of pollutant. * replaced by the vadose zone lithology. ** replaced by the contaminant impact. *** replaced by the ratio of the clay layers' thickness to the vadose zone thickness.
Water 2020, 12, 1356 4 of 22 The net recharge (R) [mm/year] represents the volume of infiltrated water that reaches the aquifer. The contamination possibility increases if the net recharge increases. Three types of recharges can be distinguished: direct, indirect, and localized [39,40]. The aquifer media (A) consists of different types of rocks serving as an aquifer.
The upper part of the vadose zone, with intense biological activity, is defined to be the soil media (S). The topography (T) (%) is defined by the terrain slope, together with its variation. A low slope will determine a small surface flow and a high pollution risk.
The vadose zone's impact (I)-The unsaturated or discontinuously saturated layer situated above the water table is called vadose. The pollutant's transfer is influenced by the vadose zone's lithology.
The aquifer hydraulic conductivity (C) is the aquifer materials' capacity to leave the water to pass through it. The aquifer vulnerability is low for reduced hydraulic conductivities.
The hypotheses of the DRASTIC models are: • The pollutants are produced at the surface of the Earth; • The pollutants are transported into the soil by precipitation; • The pollutants' travel velocity is that of the water; • The affected area must be big enough.
Firstly, a rate from 1 to 10 is assigned to each parameter, 1 being the least important [38]. Then, the DRASTIC index score is built, using the weights fixed for each parameter. The formula for DRASTIC index is: where R is the rate and w is the parameter weight. The weights have been set up by EPA (the United States Environmental Protection Agency) based on the experts' knowledge after studying different regions. In the original DRASTIC algorithm the weights range from 1 to 5 (1 being the least important), the smallest possible index score is 23 and the highest, 230. Tables 2 and 3 contain the weights and ratings of the components, firstly provided in [38]. Lower groundwater vulnerability is described by a lower index score.  Different authors [41,42] pointed out the DRASTIC results' accuracy, the small amount of input data, its application's low cost [38,43], reduced computational time, and simple computational procedure [44]. DRASTIC proved to be useful for evaluating the aquifer vulnerability in priority monitoring areas and as a valuable indicator where detailed hydrogeological evaluation is necessary. Other authors emphasized the limited validation procedure of the DRASTIC methodology [45,46] and a low correlation between the experimental data and the model's output [47,48]. Wang et al. [49] remarked on the necessity of procedure adaption for urban areas, while the parameters' weight choice in the DRASTIC index was criticized by Merchant [50]. Therefore, several approaches were proposed for improving the groundwater vulnerability estimation accuracy, each of them involving a different number of parameters. In the following, we shall present some of these methods and the rationale for their use.
Although the DRASTIC model was intended to be used in mapping applications, it was not expressly designed for use in a GIS, its initial applications employing a manual map overlay and computation procedure [50]. The main importance of vulnerability maps is that their analysis can provide effective information for making informed decisions for water management [51].
The main GIS advantage is its efficiency of combining data layers and changing the parameters used for the vulnerability classification [49]. For producing a groundwater pollution risk map (Figure 1), it is necessary to prepare the seven individual maps (one for each component in the model). Therefore, all data should be available, accurate enough [50], and introduced in a GIS database.
The main GIS advantage is its efficiency of combining data layers and changing the parameters used for the vulnerability classification [49]. For producing a groundwater pollution risk map (Figure 1), it is necessary to prepare the seven individual maps (one for each component in the model). Therefore, all data should be available, accurate enough [50], and introduced in a GIS database. The D parameter layer is generated using topographic maps. Then, the IDW method is applied for interpolating the water level data and obtaining the depth to the water layer. The general water balance equations are used for generating the recharge layer, R. The lithology, and the type of aquifer media are considered for the estimation of the A factor. The soil media, S, is determined by using the soil textural classification chart. The topographic slope, T, is determined by using a digital elevation model and the data extraction from the topography layer. Then the slope layer (%) is generated, and the range is reclassified taking into account the DRASTIC ranges. The impact of the V factor is determined using the depth to the water layer and the well logs report. The hydraulic conductivity data is retrieved by experimental measurements. Finally, the pollution risk layer is produced using the seven layers previously built by GIS, and all the DRASTIC thematic layers are combined [62]. The D parameter layer is generated using topographic maps. Then, the IDW method is applied for interpolating the water level data and obtaining the depth to the water layer. The general water balance equations are used for generating the recharge layer, R. The lithology, and the type of aquifer media are considered for the estimation of the A factor. The soil media, S, is determined by using the soil textural classification chart. The topographic slope, T, is determined by using a digital elevation model and the data extraction from the topography layer. Then the slope layer (%) is generated, and the range is reclassified taking into account the DRASTIC ranges. The impact of the V factor is determined using the depth to the water layer and the well logs report. The hydraulic conductivity data is retrieved by experimental measurements. Finally, the pollution risk layer is produced using the seven layers previously built by GIS, and all the DRASTIC thematic layers are combined [62].

Modified DRASTIC (DRASTICM)
Scientific studies pointed out that geologic structures have a significant impact on highly fractured environments' vulnerability. Therefore, in a study performed for a region form Nicaragua, Mendoza and Barmen [63] modified the DRASTIC index by including the influence of the length, connectivity, and lineament density. They introduced the lineament influence, denoted by M, in the new model, called Modified DRASTIC, whose index, MDI, is defined by where R is the rating, M is the lineament factor. A rate between 0 and 3 was assigned to the influence of the lineament. Data collected from the field and photographic interpretation were normalized and combined in a map to assess the lineament influence. This map and the other seven (from DRASTIC) contributed to building the Modified DRASTIC map.
The results show that D and T are the factors with a significant influence on vulnerability prediction. Compared with DRASTIC, the modified DRASTIC gave a better estimation of the contamination risk in zones with high fractured structures.

DRIST and Modified DRASTIC
Introduced for investigating the underground water vulnerability in Grombalia, the DRIST model was adapted to the hydrogeological system properties from this region. DRIST considers only parameters related to the unsaturated aquifer zone, while DRASTIC works with the aquifer saturated zone characteristics [40]. The calculation of the DRIST vulnerability index is similar to that for DRASTIC (but ignoring A and C parameters).
In the same article, Chenini et al. [40] proposed a Modified DRASTIC method. The difference between these approaches resides in the estimation of the factors A and I. In the new one the lithology is substituted by the permeability, as suggested in [64]. The other maps are created by the same procedure as in DRASTIC.
The permeability map of the vertical vadose zone is realized based on the vertical permeability formula: where K 1 is the vertical average permeability (m/s), H-the unsaturated zone total thickness (m), h i -the thickness of the i th layer (m), k i -the permeability of the i th layer (m/s), and p-the number of layers [64]. The saturated zone's permeability map is determined using the formula of the horizontal permeability [65]: where K 2 is the average horizontal permeability (m/s), h i , k i and p have the same significance as in h i at the denominator of formula (10) is the saturated zone total thickness (m).
Comparing the two vulnerability maps, Chenini et al. [40] showed that there are differences between them. The area with medium vulnerability is more significant in the Modified DRASTIC due to the minimization of the saturated zone effect, as an effect of the permeability replacement by lithology in the process of parameters' estimation. DRIST map reflects the effect of removing the factors related to the saturated zone.
Sakala et al. [66] used the same model and a neural network approach to generate a groundwater vulnerability model. The network used as input the DRIST parameters, and as the training dataset, the sulfate and Total Dissolved Solids (TDS) concentrations retrieved from five groundwater samples. The groundwater vulnerability model was finally obtained by applying a fuzzy operator for combining the training and classification results. The model's results are well correlated with the available data and the output of the DRIST model.

DRAV
DRAV is a model designed by modifying DRASTIC for taking into account the groundwater characteristics from the arid zones [67]. Since generally, in those areas, there is no horizontal runoff, the DRASTIC T term was removed, and S was replaced by V (vadose zone's lithology). The factors D, R, and A were kept in the new model.
The DRAV index is a linear combination of the factors D, R, A, and V with the normalized weights 0.20, 0.15, 0.31, and 0.34, respectively.
DRAV was used to analyze the pore groundwater in the northwestern part of China, but no comparison with other methods is provided. Therefore more studies are necessary to validate this approach.

DRAMIC
Many scientists emphasized the limitations of DRASTIC's application for urban areas [49,68], as follows. (1) The terrain where the cities are situated is mostly flat, so the T factor in the DRASTIC model is not relevant. (2) The values of the soil media can be hardly obtained because the ground surface is mostly covered by concrete. (3) The hydraulic conductivity is not relevant. Therefore, they built the DRAMIC index, by replacing in DRASTIC the S factor by the thickness of the aquifer (M), and the C factor by the contaminant impact (denoted by C as well). It must be noticed that DRAMIC does not consider the pollutants' properties, but its stability and infiltration capacity into the aquifer. The parameters (and ratings) in DRAMIC are [49]: • Aquifer thickness (m): 0-6 (9), 6-15 (7) Stability, infiltration easiness (9) Stability, infiltration relative easiness (7) Stability, infiltration uneasiness, and Relative stability, infiltration easiness (5) Relative stability, infiltration relative easiness (4) Relative stability, infiltration uneasiness, and Instability, infiltration easiness (3) Instability, infiltration relative easiness (2) Instability, infiltration uneasiness (1) The DRAMIC Index is computed by the relation where R is the rating. The main factors considered in DRAMIC are the stability of the pollutant and the easiness of the pollutant infiltration. The results of this model applied in a study from China (Wuhan region) were compared with the field data, showing a good correlation. Despite promising results, other studies are needed to validate this method for other urban areas.

DRASTICA
DRASTICA is a modified DRASTIC model, which includes the anthropogenic influence in urbanized environments [68,69]. A new factor (A-anthropic factor) was introduced, with the weight equal to 5. The index is computed as in DRASTIC, adding the new term, A R A w , where A R is the rating and A w the weight. The anthropic factors and the rating assigned are the following [68]: Effluents/sewage/industrial waste (untreated), Oil spillage/gas flaring and E-wastes -9, Open dumpsites (non-sanitary landfill) and Emissions from automobiles/generators -8, Cementary/soakaway/pit latrine (unlined) and Fertilizer/agrochemicals-7, Domestic waste (organic/degradable) -6, Effluents/sewage/industrial waste (treated) and Sanitary landfill -5, Cementary/soakaway/pit-latrine (lined) and Bush burning -4. The rating and weighing of the other parameters were kept as in DRASTIC.
Four vulnerability categories were built (low, moderate, high, and very high), corresponding to values of vulnerability indexes in the intervals 140-159, 160-179, 180-199, 200-215. In a study of the water pollution impact in Lucknow, India, DRASTICA better performs by comparison to DRASTIC, when the models were validated using field data. The sensitivity study emphasized that the less sensitive factors were A (aquifer), followed by S and T. The parameters with the highest impact are D, followed by A (anthropogenic factor) and C.
Another research concerning the groundwater vulnerability in the Niger Delta [69] concluded that the anthropic activity (incorporated in the A factor) had a consistent impact on the groundwater contamination.

DRASTIC-LU
Studies concerning the groundwater vulnerability showed an increasing impact of land use on water contamination [51,70,71]. Alam et al. [51] indicated that industrial and sewage pollution, pesticides, and fertilizers alter groundwater quality. They proposed a new index, DRASTIC-LU, adding "the land use pattern" (LU) parameter. The land use categories considered (and the rating) are respectively: urban and industrial (10), rural and industrial (9), rural and agriculture (8), with a weight of 5.
The DRASTIC-LU index is computed by: where the land use rating and weight are L R and L w , respectively. The other acronyms have the same significance as in the DRASTIC index. The parameter of the vadose zone impact (I R ) is computed by [70]: where T is the vadose zone total thickness, T i is the ith layer thickness and I r i is the ith layer rating. Since this approach considers many layers of the vadose zone, it is expected to provide more accurate results.
The values of the DRASTIC-LU index are situated in the interval [158,190], divided into subintervals as follows: less than 160 (corresponding to low vulnerability zone), 160-170 (medium vulnerability zone), 170-180 (corresponding to high vulnerability zone), greater than 180 (very high vulnerability).
Some research [51,70,72,73] investigated the groundwater vulnerability in different regions of India. In a study related to a zone of Central India, Alam [50] showed that the most significant parameters in the model DRASTIC-LU model are D, I, C, and LU.
In a vulnerability analysis in the Basin of Damodar River, Kumar and Khrisma [72] compared the performance of DRASTIC and DRASTIC-LU, emphasizing the significant impact of the LU component. The sequence of impact intensities I > D > C > LU > S > T > R > A resulted after investigating the map sensitivity. At the models' validation stage, a better correlation between the field data and the estimated ones resulted in the DRASTIC-LU model (0.893 against 0.781 for DRASTIC). Therefore one can conclude that the essential factors that should be taken into account for assessing the vulnerability in the study zone are A, T, I, and LU.
In the sensitivity analysis by map removal in a DRASTIC-LU approach for Karun Basin, Sinha et al. [73] found a different sequence of impact intensities by comparison with [72] (LU > S > T > D > I > A > R). Therefore, the LU and S factors have the main effect on the DRASTIC-LU index. This result is concordant with the field reality (the aquifer' shallow waters). Sensitivity analysis revealed that depth of water table, land use, and topography produce large variations of vulnerability index by comparison with other parameters.

DRASIC-LU
DRASIC-LU is a version of DRASTIC, initially used for assessing the groundwater pollution risk in some sub-regions of India (Ganga Plain) [71]. Due to the topographic small variation, the parameter T was removed from the DRASTIC index and was replaced by the parameter L (land use), which reflects the land use impact on the water quality [74]. The land use categories are the same as in the DRASTIC-LU, and the vadose zone impact parameter is computed by the relation (7). Qinghai et al. [75] introduced the hydraulic conductivity values in concordance with the experimental data. They are respectively:

•
The ratings for D (depth to the water table) are 2, 3 and 5, while the weighting factor is 5; • The rating for net recharge (R) is 9, and the weight scale is 4; • The rating for aquifer media (A) is 8, and the weight scale is 3; • The ratings for soil media (S) are 5 and 6, and the weight scale is 2; • The ratings for vadose zone impact (I) are 1 and 2, while the weight scale is 5; • The ratings for hydraulic conductivity (C) are 4, 8 and 10, while the weight scale is 3; • The ratings for land use (L) are 8, 9, and 10, and the weight scale is 5.
The new index is defined by: The terms have the same significance as in Equations (1) and (6) Studying an aquifer from the Ganga Plain, Umar et al. [74] concluded that D, C, I, and LU are the main factors to be considered for vulnerability mapping.

SI Index
Ribeiro [76] introduced the SI method for the estimation of the groundwater vulnerability to pollutants generated in areas at medium and large in Portugal. SI is obtained by removing S, I, and C from DRASTIC an including the land use parameter (LU) that incorporates the agricultural activities' impact (especially nitrates) on the water quality [77]. Therefore this method assesses the specific vulnerability of groundwater.
The SI index is computed by: where the parameters' weights are [77]: The essential land use activities classes and the corresponding rating values (displayed inside the brackets) [77] varies between 0 (for semi-natural zones and forest) and 100 [59,[76][77][78][79][80][81][82][83][84]. Hamza and Added [82] show that DRASTIC does not consider the contaminant's nature and gives great weight to the hydrogeological factors. The case study supports the idea that the SI method was designed for taking into account the nitrates properties and the relations between them and the intrinsic vulnerability. LU factor integrates the land use types, allowing the integration of different particular characteristics.
The results of Stigter et al. [77] and Hamza et al. [83] show that permeable aquifer and high recharge are responsible for the pollution vulnerability increase. For chloride or nitrate contaminants in specific conditions, the dilution potential may have a significant role in the determination of contamination degree [81]. Validating the vulnerability maps using the measured nitrites concentration, Stigter et al. [77] emphasized the groundwater vulnerability underestimation when DRASTIC was used instead of SI. Another comparative analysis of these methods validated in the field showed a better concordance when using the SI approach [83].

DRARCH
This model was introduced for studying the water vulnerability at arsenic in the Taiyuan basin and is based on simulation of the solute transport. The procedure can be summarized as follows [75]: The ratio of the clay layers' thickness to the vadose zone thickness (R), introduced for emphasizing that the clay has a specific surface area and an adsorption capacity greater than other sediments; • The coefficient of pollutant's adsorption by the sediment in the vadose zone (C); • Aquifer hydraulic conductivity (H).

of 22
where V is the DRARCH score, R-the rating value, w-the parameter weight. The vulnerability index values are between 31 and 310 and five vulnerability classes were adopted: very low , low , moderate (143-198), high (199-254), and very high (255-310).
Other approaches of the aquifer vulnerability to arsenic used a GIS-based DRASTIC [85], with the vulnerability classifications and the indices values given in [86].

SINTACS
SINTACS was proposed and developed by Civita [26] and Civita and De Maio [27,28] for improving and adapting the DRASTIC model to the particularities of Italy. The letters in SINTACS are the first letters of the Italian words that define the models' factors. They are the depth to the water table (Soggicenza), effective infiltration (Infiltrazione), attenuation capacity of the unsaturated zone (Nonsaturo), type of the soil media (Tipologia della copertura), characteristics of the saturated zone (Acquifero), hydraulic conductivity (Conducibilità), and topographic slope (Superficie topografica) [25,29].
Civita [29] remarked that for using one or another method for assessing the groundwater vulnerability, one should consider the density of the observation points, the data availability, its completeness, and reliability, the homogeneity of the study region. In a critical review of some methods he presents the reasons for searching a better approach for the evaluation of groundwater vulnerability:

•
The soil action is isolated from the action of the embedding system.

•
The climatic factors and their influence on the water system is not considered • Most methods have only a local application • The use of vulnerability maps for the prevention of the groundwater quality deterioration should be supported by a deep insight into the mechanism of the contaminant production and its risk level [21].
Based on the use of the same parameters, the SINTACS structure has a higher complexity than the DRASTIC one.
For a complete and reliable database, the SINTACS procedure is the following [25,29].
• Select the factors used in the study • Divide the factors into types or subintervals containing the factors' values • Assign a rating, P, between 1 and 10, to each subinterval, in concordance with its importance in the last step of the algorithm ( Figure 2) Water 2020, 12, x FOR PEER REVIEW 12 of 22 strings of weights (Table 4) to describe the environmental conditions [6]. The most difficult task remains the range selection and the assignation of weight and ratings.   [29]).

Parameter
P i being the rating value and w i is the corresponding weight.
One of the SINTACS advantages is the possibility of simultaneous use in different zones since each situation has assigned a specific weighting rate. Notice the differences between the DRASTIC and SINTACS procedures of weighting and rating, the last one operating in parallel with different strings of weights (Table 4) to describe the environmental conditions [6]. The most difficult task remains the range selection and the assignation of weight and ratings.
For extending the applicability of SINTACS to the entire Italian territory, a new approach was introduced, by combining the SINTACS Release 5 [28] with the GNDCI_CNR Basic Method [29]. Since its release, SINTACS became one of the most used methods for the assessment of groundwater vulnerability in countries as Algeria, Italy, Jordan, Morocco, Thailand [87][88][89][90][91][92][93].
Corniello et al. [89] remarked that lithological and morphological settings play an important role in the process of generating SINTACS vulnerability maps. In a comparative study of three methods [90] on sites situated in a Mediterranean region, it is shown that the climatic conditions have a significant influence on the methods' performance, DRASTIC providing better results than SINTACS and AVI. A comparative study of the vulnerability maps produced by DRASTIC and SINTACS for an aquifer situated in Algeria [91] shows that the results are statistically concordant. Luoma et al. [92] emphasize in their research on a coastal aquifer that the SINTACS vulnerability maps are concordant with the field reality.
From the comparative analysis of the results provided by performing DRASTIC, SINTACS, and GOD methods on a database from Central Romania [93], one can remark on the similarity of the maps generated by the first two methods, there are few differences in the extent of the class of low vulnerability. In zones with small vulnerability variations, GOD performed worst. Therefore this method should be used only for regions with big vulnerability variations.
Aiming at detecting the capabilities of five groundwater vulnerability approaches, Civita and De Regibus [11] developed their research in three zones (mountains, hills, and flat). SINTACS and DRASTIC could adapt to the various situations, by comparison to the other competitors (GOD being among them) due to their flexibility.
Secunda et al. [94] and Noori et al. [95] used a SINTACS-LU approach in their research. The new factor, LU, was introduced by analogy to DRASTIC-LU, for considering the land use effect on the groundwater vulnerability. The new approaches better performed than SINTACS in case studies from Israel and Iran. Both SINTACS-LU and DRASTIC-LU vulnerability maps delimited the zones highly affected by human activity [94]. The sensitivity analysis for SINTACS-LU [95] showed that the parameter with the highest impact was the vadose zone, followed by the land use. The analysis of the correlation between the vulnerability index and the nitrate values (recorded on-field) was the highest for SINTACS-LU (0.75), by comparison to those of DRASTIC-LU (0.68) and SI (0.64).

Groundwater Vulnerability Assessment to Specific Pollutants
Even if all the described methods could be applied to assess the groundwater vulnerability to contamination, new approaches have been proposed to account for the specific properties of some pollutants. These include Pesticide DRASTIC, Pesticide DRASTIC-LU [38,96], Modified Pesticide DRASTIC [80,81,96], Modified DRASTIC for nitrate [20,46,[96][97][98][99]. The factors weights in Pesticide DRASTIC and Pesticide DRASTIC LU differs from those of DRASTIC, the rating being preserved. They are presented in Table 5. The ratings of the land use in Pesticide DRASTIC-LU are 1, 5, 7, or 8. DRASTIC, Pesticide DRASTIC, and Pesticide DRASTIC-LU were used for a study in a part of the Gangetic Plain with intense agricultural activities. Statistical analysis of the average values revealed that the most significant contribution to calculation the vulnerability indexes were I, T (in DRASTIC and Pesticide DRASTIC), and D, followed by R (in Pesticide DRASTIC LU). The sensitivity analysis found that A and R factors had the highest impact on all the models. The less significant parameters were S, and T-in DRASTIC, I and C-in Pesticide DRASTIC. Pesticide DRASTIC was the best model point of view of the correlation between the field data and prediction [96].
DRASTIC, Pesticide DRASTIC, and SI were applied in a case study of an aquifer from Tunisia. SI and Pesticide DRASTIC better detected the pollution risk. The concordance between the categories of vulnerability determined by these approaches was 64%. The authors [80] recommend the use of these two approaches for different purposes; the first one for monitoring, whereas the second one as a part of a multicriteria decision tool for allocating different zones to specific anthropic activities.
The performance of the same three models, together with Modified DRASTIC were compared on an aquifer in India (Telangana) [81]. The D factor has a considerable impact, followed by soil, the smaller one being that of R. The vulnerability classes are almost the same in SI, modified DRASTIC, and modified Pesticide DRASTIC because of the effect of LU inclusion. The Modified Pesticide DRASTIC map contains a higher area with high vulnerability, compared to Pesticide DRASTIC. The scientists remarked [81] the DRASTIC vulnerability underestimation and SI overestimation. All the models (but SI) are in concordance by at least 60%. It seems that the modified Pesticide DRASTIC provided the best predictions.
The nitrate is not a natural compound of soil, being the result of human activities, like the fertilizers used in agriculture or defecation [100]. While some authors used the Pesticide DRASTIC or Pesticide DRASTIC-LU to study the soil contamination with nitrates [46], other scientists [97][98][99] developed new approaches for improving the weights assigning in DRASTIC. Antonakos and Lambrakis [98] proposed DRASTIC-based hybrid methods, Huan et al. [99] adjusted the DRASTIC rating and weighting system. They validated their models on study cases from Greece and China, respectively. Kazakis and Voudouris [97] replaced the A, S, and I factors of DRASTIC by the thickness of the aquifer, losses of nitrogen from the soil, and hydraulic resistance. The second factor was estimated by the GLEAMS model [101]. The parameters' range and weights were also modified. The two new methods, named DRASTIC-PA and DRASTIC-PAN, were compared to DRASTIC and LOSN-PN. Their performance were validated by the sensitivity analysis.
As already mentioned, one of the main criticism of DRASTIC was that it does not take into account the study particular characteristics of each study region, and does not adapt the ratings and weights [46]. To surpass this inconvenience, other techniques have been proposed, as follows: • Approaches that use genetic algorithms and neural networks [45,[115][116][117].

•
Correspondence Analysis [118], aiming at minimizing the redundancy between factors using a multivariate statistical method • Calibration techniques, proposed in [98,119], used in [99,100]. • For a deep insight into these approaches, the readers can access the cited articles.

Models' Validation
Different authors used many models to validate the results of the vulnerability maps [25,90,120]. Kumar et al. [21] emphasize that this comparison is not advisable because various approaches use different parameters, so the vulnerability maps might not be similar. The benefit of such procedures is to offer an insight into the existence and spatial distribution of groundwater pollution. Therefore, other techniques should be used, as the validation of the vulnerability maps on contaminants data sets collected on-site from wells distributed in the study region. This is usually done using the concentrations of nitrates in the collected samples. A method whose results are most contrasting could be considered most sensitive, so it can be used [6].
Napolitano-Fabri [121] proposed the single-parameter sensitivity analysis (SPSA), which is the most frequently used technique for evaluating the significance of the parameters in the vulnerability models [46,57,58,64,99,[122][123][124][125][126][127]. SPSA provides information on the rating and weighting assigned to each parameter, enabling its evaluation by the researcher. SPSA compares the theoretical and effective weights assigned to each parameter in a model.
Lodwick et al. [128] introduced the map removal sensitivity analysis (MRSA) for assessing the uncertainty degree of the models' output. MRSA consists of removing one map from the analysis of vulnerability and computing a variation index.

Conclusions
In this article, we reviewed DRASTIC and the main DRASTIC-like approaches proposed by scientists for improving the initial algorithm. The methods for assessment of the groundwater vulnerability in karstic regions were not discussed here. DRASTIC is among the most used tools for groundwater vulnerability evaluation. Generally, it uses readily available geodata with no experimental data. DRASTIC employs numerous parameters, and its outputs are only sometimes compared with field-collected data. Therefore DRASTIC-based forecast should be rigorously checked before making management decisions.
The groundwater vulnerability maps are important tools for assessing the groundwater vulnerability and planning future land use. No method developed for creating vulnerability maps is the most reliable, each of them depending on the aquifer characteristics, the land use, the data availability, the parameters involved in the model, the weightings, and rating assigned to each parameter.
Some aspects should be addressed in the next studies: • Development of analytical methods for choosing and validating the ratings and weightings attached to each parameter in the models • Integrating the models of water flow and pollutants' transport in different soils types in the methodology of choosing the weighting values of different parameters • Detecting the relationships between the parameters used in the models by statistical methods and removing the effect of this correlation by adjustment of the ratings and weightings attached to the corresponding parameters • Development of unified models that should include the soil and geological characteristics • Development of hybrid models to reduce the influence of subjectivity in the parameters' settings and use the statistical methods for the results' validation. • Improvement of the databases containing hydro-chemical elements and their integration into GIS software • Improvement of GIS software by integrating analytical methods with groundwater vulnerability methods • Development of spatio-temporal methods for the groundwater vulnerability assessment.