1. Introduction
In many countries, groundwater has a significant contribution to fresh water supply for various anthropogenic activities associated with the industry, irrigation, human consumption and economic development [
1]. During the last few decades, an unprecedented development in the utilization of groundwater resources has been seen, with estimates suggesting that groundwater covers more than 30% of the world’s water demand [
2]. Especially in regions characterized by arid, semi-arid or Mediterranean climates, groundwater often represents the primary potable water source, particularly where scarcity of surface water is substantial.
Notwithstanding, groundwater quality and quantity are at high risk due to the haphazard development of groundwater resources [
3,
4]. In particular, numerous anthropogenic activities, such as population growth, rapid urbanization, deforestation, intensive agriculture, alterations in land use and industrialization, lead to serious groundwater quality degradation and depletion [
5,
6,
7]. Furthermore, climate variability is a climatic factor that has the propensity to stress groundwater and exacerbate the situation by influencing groundwater recharge markedly [
3,
8]. In recent years, in various agricultural areas, the extensive application of chemical fertilizers has caused substantial nitrate pollution of aquifers, which in turn, has crucial effects on public health and the ecosystem [
7,
9,
10,
11]. Hence, the safeguarding of groundwater and pollution prevention are fundamental prerequisites for efficient groundwater management and the socioeconomic advancement of an area [
12,
13].
A cost-effective key tool for the prevention and reduction of groundwater pollution can be the assessment of groundwater vulnerability through the utilization of various techniques and methodologies, as proposed in the literature [
14,
15,
16]. A typical classification of the groundwater vulnerability methods comprises three primary categories, namely: (i) statistical assessment, (ii) process-based simulationand (iii) index and overlay methods [
17]. Among the most commonly used index and overlay methods in porous aquifers are DRASTIC [
18], GOD [
19], AVI [
20], SINTACS [
21] and SI [
22]. Moreover, groundwater vulnerability assessment can be divided into intrinsic, specific and integrated. Regarding the intrinsic, it represents the vulnerability based on the geological and hydrogeological characteristics of an aquifer, whilst the specific describes the vulnerability against specific pollutants. Both types are considered essential and useful for the safeguarding of groundwater resources and the protection of groundwater quality from pollution [
23,
24]. Additionally, the integrated vulnerability assessment, alongside the intrinsic vulnerability, incorporates the presence of danger centers, namely construction material factories, sewage outlets, gas stations, livestock farms, proximity to main roads, rivers and residential areas, as well as landfills, olive mills, animal barns and stockpiles of fertilizers and pesticides [
6,
25,
26].
One of the most popular and extensively used index methods for groundwater vulnerability mapping is DRASTIC, which was originally developed by the US Environmental Protection Agency (EPA) [
6,
27]. This method utilizes seven hydrogeological parameters controlling the flow of the water from the surface to the aquifer, namely depth to water table (D), net recharge (R), aquifer media (A), soil media (S), topography (T), impact of the vadose zone (I) and hydraulic conductivity (C) [
28]. A significant advantage of DRASTIC is its adaptability and appropriateness for regions with varying hydrogeological characteristics and sources of pollution, such as agricultural, urban and industrial [
7,
29,
30,
31]. Despite being the most established method for regional vulnerability assessment, DRASTIC has its limitations, including the subjectivity and uncertainty in the delineation of the rates and the weights of each parameter [
32]. Thus, various techniques have been employed to modify and enhance the reliability and precision of the method and minimize the subjectivity in the evaluation of each parameter’s ratings and weight. These techniques include sensitivity analysis, analytic hierarchy process (AHP), Fuzzy AHP method, binary classifier calibration with a genetic algorithm, multiple linear regression, multi-criteria decision-making approach coupled with a metaheuristic algorithm, as well as additional factors, such as land use (DRASTIC-LU) have been incorporated [
33,
34,
35,
36,
37,
38,
39,
40,
41,
42,
43]. In particular, the impact of land use, such as agriculture, is considered to have a noteworthy impact on groundwater quality and thus warrants inclusion in groundwater vulnerability assessment [
44].
Martínez-Bastida et al. [
45] proposed the Nitrate Vulnerability (NV) index, which is a novel approach based on the standard DRASTIC and the incorporation of land use parameters (as a potential nitrogen source). The new specific vulnerability index is based on a multiplicative model and is concentrated on nitrate pollution by adding an additional parameter denoted as the “potential risk associated with land use” (LU). Results indicated that the NV index achieved greater accuracy in the assessment of specific vulnerabilities compared with DRASTIC-LU [
45,
46].
Another significant and widely employed low-cost vulnerability method for enhancing and adapting DRASTIC to the particularities of Mediterranean regions, such as Italy and Greece, is SINTACS [
21,
47]. The aforementioned method, developed in Italy, employs the same hydrogeological parameters as DRASTIC, albeit it has a more flexible weighting and rating system that is contingent upon regional-specific conditions [
48]. The inherent subjectivity in the parameters’ weight and rating is a major limitation. To overcome this drawback, the SINTACS method can be enhanced by integrating land use parameters into its calculations (SINTACS-LU), reporting greater performance and giving more accurate results over the standard version [
48,
49,
50]. In general, the method can be employed for medium to relatively large-scale porous-media aquifers, whereas in fractured, karstic or carbonate aquifer systems, the results might be debatable.
This study aims to assess and identify the groundwater vulnerability zones to specific pollution in the NeaMoudania aquifer, Chalkidiki, Greece. However, the application of a single groundwater vulnerability method in a given area might, in some cases, lead to imprecise results. Therefore, our study endeavors to conduct a comparative analysis by employing an ensemble of groundwater vulnerability methods, namely Pesticide DRASTIC-LU, Nitrate SINTACS-LU and NV index, in order to evaluate the efficacy and appropriateness of these methods in an agricultural area. The verification of the produced vulnerability maps of the aforementioned methods was appraised through the correlation between the nitrate concentration (NO3−) in groundwater and the vulnerability index. Furthermore, the effective weight of each parameter in Pesticide DRASTIC-LU, Nitrate SINTACS-LU and NV index was evaluated by performing a single parameter sensitivity analysis (SPSA).
2. Study Area
The watershed of Nea Moudania (
Figure 1) is located in the southwestern part of the Chalkidiki peninsula, in the Region of Central Macedonia, Northern Greece, and it is a constituent of a larger area which is known as Kalamaria Plain. This coastal alluvial basin occupies an area of approximately 78 km
2, and it belongs, administratively, to the municipalities of Nea Propontida and Polygyros. Morphologically, the average altitude is relatively low (approximately 200 m above sea level) and features a topography with gentle slopes (mean soil slope approximately 2%). This area forms the primary agricultural land of Chalkidiki [
51]. The climate is typical thermo-Mediterraneanand is characterized by high temperature and aridity during the summer months, while the winters are mild and wet. The mean annual precipitation is approximately 420 mm [
51]. The watershed of the Nea Moudania basin belongs to the Peonia geological subzone, which is part of the Axios geotectonic zone [
52]. The Nea Moudania aquifer is covered mainly by Neogene (upper Pliocene-middle Pleistocene) sediments, which cover the pre-Neogene bedrock formations. The Neogene sediments consist of alternated beds of red silty clay with fine-grained quartz and mica intercalated by sandstones, sands, marl, silts, and conglomerate lances. The pre-Neogene bedrock consists of recrystallized bluish limestone, greenish-brown two-mica gneiss and muscovite gneiss. Moreover, the study area comprises Quaternary formations, such as alluvial deposits (sands, gravels, conglomerates and clay), coastal deposits (beach ridges) and lagoon deposits (sands and clay sand) [
53]. Hydrogeological interest is primarily concentrated on the recent deposits due to their considerable capacity for water storage and their substantial thickness [
54]. On the contrary, the pre-Neogene bedrock formations are considered impermeable with limited fractured zones; therefore, they do not exhibit hydrogeological interest.
The fairly flat topography and the favorable regional climate conditions make the study basin an important agricultural area. Indeed, approximately 80% of the total area is extensively cultivated and irrigated, and there is also a notable presence of touristic and urban development [
55]. The area has a permanent population of approximately 16,000 inhabitants, whilst the maximum population, mostly during the summer season, surpasses 40,000 people [
56]. Water demand for irrigation and domestic purposesis high, especially during the period spanning from May to September. Moreover, the combination of surface water scarcity and generally inadequate yearly rainfall renders groundwater the only viable source of water [
57]. However, this fact has caused a noteworthy decrease in the groundwater level throughout the area. In addition to that, since the 1990s, overexploitation via numerous deep wells has led to quantitative issues. Another major challenge over the years is the qualitative deterioration of the aquifer due to the intensive agricultural activities that result innitrate pollution and seawater intrusion. Undoubtedly, the implementation of a sustainable approach to managing water resources is henceforth absolutely necessary to ensure the quantitative and qualitative protection of the Nea Moudania aquifer.
5. Conclusions
In this study, three different GIS-based modified groundwater vulnerability assessment methods, namely the Pesticide DRASTIC-LU, the Nitrate SINTACS-LU and the Nitrate NV index were implemented in order to evaluate and identify different groundwater vulnerability zones in the porous aquifer of Nea Moudania, Greece. All three methods deploy eight geological and hydrogeological parameters, which predominantly affect the vulnerability of the porous aquifer. Particularly, the areas with the highest vulnerability potential are mainly distributed in the southern and southwestern parts of the study basin. This is attributed to the shallow groundwater depth, the intense anthropogenic activities, such as agriculture and the topography, which is characterized by gentle slopes. Conversely, the northern region of the Nea Moudania aquifer is identified as a low groundwater vulnerability zone, where the topography is slightly steeper, the groundwater table is generally deeper, the vadose zone is thicker and the constituent materials are less permeable, as well as because the agricultural activities are less intensive. Thereafter, a single-parameter sensitivity analysis was performed, revealing the importance of the land use and the topography parameters in vulnerability delineation, emphasizing the necessity for accurate, detailed and representative pertaining to these parameters. On the contrary, the least effective parameter is hydraulic conductivity because of its relatively low value in the study area. As a consequence, the utilization of sensitivity analysis proves to be a handy tool for revising the parameter weights, helping to obtain reasonable results according to the specific hydrogeological conditions of an area. Furthermore, the water quality parameter, namely nitrate NO3− concentration, has been adopted to validate the aquifer vulnerability maps by employing the coefficient of determination (R2). In the Nea Moudania aquifer, the most accurate and suitable method for groundwater vulnerability evaluation was found to be the Nitrate NV index (R2 = 0.69), having scored a marginally superior correlation in comparison to the other two methods. Hereupon, the Pesticide DRASTIC-LU and the Nitrate SINTACS-LU also provided a noteworthy positive coefficient of determination, 0.67 and 0.66, respectively. The aforementioned results may serve as a guideline for the sustainable and effective groundwater management and protection of the Nea Moudania porous aquifer, assisting policymaking authorities in decision-making at a regional scale. Therefore, the design and systematic maintenance of a groundwater quality monitoring network on a regular basis, particularly in highly vulnerable zones, is crucial for the prevention of further groundwater quality deterioration, and thus, it is highly recommended. Furthermore, it is an urgent need to reduce the pumping flow rates of the irrigation wells (regulate the amount of water withdrawal) and enforce strict regulations on agricultural activities, as well as reschedule the general pumping strategy to attenuate groundwater depletion in the studied aquifer. Additionally, the implementation of sustainable agricultural practices, such as precision irrigation and optimal crop allocation and water-saving irrigation methods, is recommended. Moreover, the adoption of advanced technologies, such as Geographic Information Systems (GIS) and remote sensing, can facilitate the systematic monitoring of groundwater level and quality, the identification of potential sources of pollution and the planning of recharge activities. Finally, new data collection, such as groundwater level data, precipitation data, land use changes and hydrochemical data, could enhance the performance and reliability of the groundwater vulnerability assessment methods, ensuring the sustainability of the study area’s aquifer.