Groundwater Vulnerability Assessment in the Metaponto Coastal Plain (Basilicata, Italy)

: This study aims at a groundwater vulnerability assessment of the Metaponto coastal plain, located in the Basilicata region (southern Italy). In the last century, intensive agriculture, zootechnical and industrial activities have signiﬁcantly changed the plain. These changes led to negative impacts on the hydrogeological system intensifying the risk of the aquifer to pollution. The paper presents the assessment of the intrinsic vulnerability of the coastal aquifer carried out by the GIS-based application of the SINTACS method. It considers several aquifer parameters such as water table depth, effective inﬁltration, unsaturated conditions, soil media, aquifer media, hydraulic conductivity and topography. Furthermore, the anthropogenic inﬂuence in the study area was considered by applying the SINTACS-LU method, in which the parameter of land use (LU) was added. The SINTACS and SINTACS-LU vulnerability indexes were provided by summing the product of ratings and weights assigned to each parameter. The analysis of the intrinsic vulnerability map allowed for determining three classes ranging from low to high vulnerability. In both cases, the southeastern part of the coastal plain, closest to the sea, shows the highest vulnerability class, indicating that it is the most vulnerable to contamination due to the hydrogeological intrinsic factors. The wide central part of the study area shows a moderate class of vulnerability and the low class is scattered in small parts in the northern portion of the plain, which represents the areas less contaminable in space and time in the case of potential pollution. In the SINTACS-LU map, some areas classiﬁed as highly vulnerable in the SINTACS method show a minor vulnerability class. These areas are localized in natural and wooded sectors of the Metaponto plain, which are less populated, where human impact on the groundwater is minimal.


Introduction
In the 21st century, water scarcity and its pollution emerged as relevant problems worldwide. The protection of groundwater quantity and quality is therefore one of the main challenges to be undertaken in order to ensure people have availability of freshwater, which is highly threatened by overexploitation and contamination caused by anthropogenic activities and by the evident impacts of the climate-driven changes in water disposability. The authorities and stakeholders have an important role in the integrated planning of water resources, their management and protection, to ensure the best available drinking-water quality [1].
Nowadays, groundwater represents an important water supply resource subject to pollution, especially in agriculture-devoted areas having intensive activity involving the use of fertilizers and pesticides.
In the last decades, many studies of the international scientific community focused their attention on groundwater protection and management based on the groundwater vulnerability assessment. Groundwater vulnerability assessment has been widely recognized for its spatial prediction capability of groundwater to pollution [2][3][4][5][6]. It represents a  [54] and 1:50,000 scale [55]).
The marine-terrace deposits (Middle-Upper Pleistocene), constituted by thin gravel, sand and silt wedges organized in superimposed layers, unconformably overlie the substratum made by marine grey-blue silty clays of the sub-Apennine Clays Formation  [54] and 1:50,000 scale [55]).
The Coastal Deposits defined the sandy beaches, which become sandy or gravely sandy, with pebbly lenses towards the Sinni River [57]. Coastal dunes and marshy areas made up of weakly cemented sands border the beaches inland. Wetlands and coastal dunes The lithostratigraphic setting of the Metaponto Coastal Plain, derived from the crosssection AA' (Figure 2) of the Geological Map of Italy 1:50,000 scale [55] and realized based on several datasets referring to the boreholes present in the study area, highlight different stratigraphic units [55] (Figures 1 and 2). From the bottom to the top, t is composed of the irregular depth of shelf-transition clayey, clayey-silty, and clayey-sandy deposits, locally with intercalations of discontinuous gravel levels, belonging to sub-Apennine Clay Formation (Middle-Pleistocene), reaching about 120 m b.s.l. in correspondence with the paleo river valleys [53], mainly filled by estuarine deposits [49,51]. Fluvial and/or deltaic sandy-gravelly deposits with locally clayey and clayey-silty deposits belonging to marine-terrace deposits (Upper Pleistocene) characterize the second unit.
Water 2022, 14, x FOR PEER REVIEW 6 of 24 conductivity of the plain deposits are 2.28 × 10 −4 and 6.53 × 10 −5 m/s, respectively [59,60]. The hydraulic conductivity can locally decrease to very low values for the silty-clayey levels.
Even if the direct rainfall infiltration in the study area is scarce, both due to the low rate of effective rainfall and the effect of the low-permeability soils [59,60], the groundwater of the coastal aquifer is assured by the discharge from upward aquifers, in which the recharge areas are localized, and by the river leakage. The trend of the piezometric surface, almost striking parallel to the coastline, highlights that the preferential groundwater flow pathways are direct from the marine terraces and the rivers' valley alluvial deposits to the Ionian coast.

SINTACS and SINTACS-LU Methods
The SINTACS method, belonging to the GIS-based overlay and index approaches [3], was selected and applied in the Metaponto coastal plain to assess the intrinsic vulnerability to pollution. The SINTACS method, widely applied for determining groundwater vulnerability, represents the evolution of the DRASTIC method [23] suitable to the hydrogeological Mediterranean environments [45,62].
In this method, the hydrogeological scenarios characterized by the weights attributable to each parameter are more numerous and flexible than in the DRASTIC, which provides strings of weights for only two hydrogeological contexts. The index-based SINTACS method is easy to apply, low-cost, requires limited datasets and low computational time, and provides adequate reliability and accuracy of the results [42].
The acronym "SINTACS" derives from the initials of the Italian words that express the input conditioning factors used by the method to assess the intrinsic vulnerability of the Metaponto coastal plain. These parameters are: (i) S (Soggiacenza) water table depth; (ii) I (Infiltrazione efficace) effective infiltration; (iii) N (Non-saturo) unsaturated zone; (iv) T (Tipologia della copertura) soil media; (v) A (Acquifero) aquifer characteristics; (vi) C The last two lithostratigraphic units constitute the buried and outcropping coastal prisms of Metaponto, both belonging to the Alluvial, transitional, and marine deposits. Sandy, gravelly and silty deposits constitute the buried prism (Upper Pleistocene). They have filled the paleo river valleys during low-stand of sea level, the depth of which ranges between 30 and 40 m below the current sea level. These deposits, laterally discontinuous, lie by erosional surface on the substrate and the sandy-gravelly deposits. The outcropping deposits of the coastal prisms (Upper Pleistocene-Holocene) constitute the upper unit. These medium-to-fine sandy, gravelly and silty deposits of a continental to transitional environment during a lowstand of relative sea-level, filled the incised paleovalleys, and, subsequently, the aggradation took place to form the current Metaponto coastal area.
The peculiar hydrogeological structure of the coastal plain of Metaponto reflects the complex and articulated geological and stratigraphic context. The distinctive boundary conditions of the groundwater flow have been determined by the spatial distribution of the different formations and the geomorphological evolution of the entire area. The grey-blue clays, characterized by very low hydraulic conductivity, constitute the bottom of the entire aquifer system [58,59]. Different mostly unconfined aquifers can be distinguished [60,61]: two aquifer typologies are located in the marine terraced deposits and alluvial deposits of river valleys. The third one is located in the interfluvial coastal plain: the plain deposits create a shallow sandy aquifer, which is widely exploited. The hydraulic conductivity of the marine terraced deposits is generally medium, with maximum values of 10 −3 m/s [59], in correspondence with the sandy and gravelly sediments, and becomes very low for silty-clayey levels. In the coastal plain, alluvial sediments are more distinctly sandy and form transitional deposits consisting of sand, gravel and silt in layers variously distributed in space. On a large scale, the mean and median values of the hydraulic conductivity of the plain deposits are 2.28 × 10 −4 and 6.53 × 10 −5 m/s, respectively [59,60]. The hydraulic conductivity can locally decrease to very low values for the silty-clayey levels.
Even if the direct rainfall infiltration in the study area is scarce, both due to the low rate of effective rainfall and the effect of the low-permeability soils [59,60], the groundwater of the coastal aquifer is assured by the discharge from upward aquifers, in which the recharge areas are localized, and by the river leakage. The trend of the piezometric surface, almost striking parallel to the coastline, highlights that the preferential groundwater flow pathways are direct from the marine terraces and the rivers' valley alluvial deposits to the Ionian coast.

SINTACS and SINTACS-LU Methods
The SINTACS method, belonging to the GIS-based overlay and index approaches [3], was selected and applied in the Metaponto coastal plain to assess the intrinsic vulnerability to pollution. The SINTACS method, widely applied for determining groundwater vulnerability, represents the evolution of the DRASTIC method [23] suitable to the hydrogeological Mediterranean environments [45,62].
In this method, the hydrogeological scenarios characterized by the weights attributable to each parameter are more numerous and flexible than in the DRASTIC, which provides strings of weights for only two hydrogeological contexts. The index-based SINTACS method is easy to apply, low-cost, requires limited datasets and low computational time, and provides adequate reliability and accuracy of the results [42].
The acronym "SINTACS" derives from the initials of the Italian words that express the input conditioning factors used by the method to assess the intrinsic vulnerability of the Metaponto coastal plain. These parameters are: (i) S (Soggiacenza) water table depth; (ii) I (Infiltrazione efficace) effective infiltration; (iii) N (Non-saturo) unsaturated zone; (iv) T (Tipologia della copertura) soil media; (v) A (Acquifero) aquifer characteristics; (vi) C (Conducibilità idraulica) hydraulic conductivity and (vii) S (Superficie topografica) topographic slope.
Based on the structure of the SINTACS method [3,26,44], each hydrogeological parameter is divided into different classes, to which it is assigned ratings and relative weights.
The SINTACS index is provided using the following Equation (1): In this study, one more parameter was considered in the SINTACS method with the aim to improve the accuracy in representing groundwater vulnerability: the land use (LU) was included to take into account the anthropogenic influence in the study area. Land cover plays a key role in relation to the pressure exerted on groundwater quality by industrial, urban and vegetation land uses, posing a serious threat to groundwater resources [19,20]. This parameter was ranked depending on land use and weighted with a 5 value ( Table 1). The modified SINTACS index was computed using the following Equation (2): SI NTACS − LU index = S r S w + I r I w + N r N w + T r T w + A r A w + C r C w + S r S w + LU r LU w (2) where for both equations, S, I, N, T, A, C, S and LU are the above-reported parameters; r is the rating value; w is the weight assigned to each parameter. The data implementation and processing were successfully performed in a GIS environment to assess the vulnerability of groundwater with the SINTACS methodologies. The advantages of GIS are identifiable in the ability to store, organize and analyze uniquely georeferenced spatial geographic data and the efficiency in combining different layers of information. Eight layers (one for each parameter) were prepared and processed to produce the groundwater vulnerability maps and then reclassify the indexes ( Figure 3). is the rating value; w is the weight assigned to each parameter.
The data implementation and processing were successfully performed in a GIS environment to assess the vulnerability of groundwater with the SINTACS methodologies. The advantages of GIS are identifiable in the ability to store, organize and analyze uniquely georeferenced spatial geographic data and the efficiency in combining different layers of information. Eight layers (one for each parameter) were prepared and processed to produce the groundwater vulnerability maps and then reclassify the indexes ( Figure 3).  Each parameter is evaluated depending upon its relative importance to the potential intrinsic vulnerability and indexed by attributing rating ranges from 1 to 10, according to Table 1.
To each parameter is subsequently attributed a multiplying weight ranging from 1 to 5 (Table 1). The SINTACS method proposes six different strings of weights related to the different hydrogeological and anthropogenic conditions, to which the Land Use (LU) parameter was added ( Table 2). The choice of the string of weights "Nitrates" was made considering the mainly agricultural vocation of the study area [44]. To the relevant parameters is assigned a weight of 5, whereas the less important are assigned a weight of 2 (Table 1). Finally, the vulnerability indexes are calculated by Equations (1) and (2). The higher intrinsic vulnerability to pollution is indicated by the higher SINTACS index. The data collection and the surveys were carried out in the overall Metaponto plain. The elaboration and preparation of the SINTACS thematic layers were performed with QGIS software and rasterized with grids of 5 m cell size, using the Kriging approach [63] ( Figure 3). The Kriging technique provides an efficient method for the interpolation of the limited number of observations such as wells, rainfall and temperature data for hydrogeological water balance, hydraulic conductivity determined in punctual sites and for preserving the theoretical spatial correlation [64]. The elaborated data were retrieved as follows.
4.1.1 The water table depth (S) is the distance between the ground surface and the water table and indicates the depth that a contaminant must pass through to reach groundwater [46]. It affects the contaminant in biological and chemical reactions, and the time required to reach groundwater [47]. A low depth to groundwater leads to a higher vulnerability rating. This parameter originated from 42 selected monitored wells' water level measurements located in the study area ( Figure 1). The layer was elaborated based on their localization and the water table depth spatial distribution map was obtained through the Kriging interpolation technique in terms of raster data.
4.1.2 The effective infiltration (I) parameter, is the amount of infiltration that penetrates from the ground surface and reaches the aquifer, and can transport contaminants to the water table through the unsaturated zone. If the effective infiltration increases, the contamination potential of groundwater will rise too. The effective infiltration in the study area was evaluated by the estimation of the hydrogeological budget components applying the inverse hydrogeological water balance approach [65,66]. The effective infiltration results directly related to groundwater contamination. The hydrogeological water balance of the study area, due to the low hydraulic conductivity of the outcropping soils [67] and the unfavorable climatic conditions, shows low values of the effective infiltration rate. 4.1.3 Unsaturated zone (N) characteristics and thicknesses influence groundwater intrinsic vulnerability. The vadose zone included between the bottom of the soil horizon and the water table controls the dilution and the contaminant rate that reach the saturated zone. The ratings assigned to the unsaturated zone classes are given in Table 1. This parameter was elaborated considering lithological characteristics of the geological formations and available stratigraphic data of wells.
4.1.4 Soil media (T) affects the transport mechanisms of the contaminant from the soil surface to the water table. Aquifer potential vulnerability depends upon soil characteristics such as soil organic materials percentage, texture, and permeability [46]. Soil media plays a key role in assessing the vulnerability index of groundwater as it restricts the vertical flow of contamination into the subsurface. In addition, it has a significant impact on the amount of recharge that flows through the soil. Soil data was taken from the pedological map of the Basilicata Region [68].
4.1.5 The aquifer media (A) parameter refers to the characteristics of the saturated zone, such as the porosity rate, particle types and size [46]. The relative layer was prepared using geological data, lithological maps and stratigraphic borehole data available for the study area. The ratings assigned to the aquifer media parameter are reported in Table 1. 4.1.6 Hydraulic conductivity (C) measures the capability of the aquifer to transmit water and controls the rate at which the groundwater will move when subjected to a given hydraulic gradient. This parameter controls the transport and distribution of contaminants from the injection point inside the saturated zone [46]. High permeability allows more infiltration rates and higher concentrations of contaminants to reach the saturated zone. The high permeability of the aquifer system provides a high vulnerability index. The hydraulic conductivity values were obtained from many pumping tests carried out in the study area in the last twenty years [59]. These values were used to elaborate the hydraulic conductivity map.
4.1.7 Topography (S) indicates the slope of the topographic surface and influences the vulnerability assessment in relation to the fact that water and pollutant can run off or stay on the surface for a time long enough to infiltrate. For low topographic slope, the infiltration rate is higher, and, consequently, contaminants have a greater probability of infiltrating the aquifer. With larger slopes, the infiltration capacity decreases and the surface is less vulnerable to groundwater contamination. Thus, by decreasing the slope, the vulnerability of the aquifer increases. The slope map was derived from the digital elevation model data covering the investigated area and downloaded from the Basilicata Region website with a spatial resolution of 5 m and it is expressed as a percentage of the steepness.
4.1.8 Land use (LU) impact on vulnerability was considered by including the land use parameter in the SINTACS approach. The land use data were derived from Corine Land Cover 2018 [69]. In the study area, the groundwater quality is subjected to degradation due to the anthropic impact of the urban, agricultural and industrial pressures. In particular, the Metaponto coastal plain aquifer system is conditioned by intensive agricultural activities and other pollution sources.

Sensitivity Analysis
The selection of ratings applied to the parameters included in SINTACS and SINTACS-LU methods is inevitably linked to subjectivity. Because such a choice significantly impacts on the final vulnerability maps, the sensitivity analysis can be performed to understand the influence of the input parameters on the output, assess the consistency of the results, individuate which variables are more relevant and hence require more thorough information and accuracy [70].
In this research, map removal and single variable analysis were the two types of sensitivity analysis performed with a quick and low-cost GIS-based approach.

Map Removal
The goal of the map removal analysis is to figure out if all the parameters are necessary for the computation of the vulnerability index, or if one or more of the parameters may be deleted without having a significant impact on the outcome. The analysis was performed for both SINTACS and SINTACS-LU methods by deleting one parameter at a time and assessing the effect on the output by calculating the following sensitivity measurement in each cell in which the study area was discretized [71][72][73]: where S is the sensitivity measurement expressed in terms of a variation index, V is the unperturbed vulnerability index obtained using all the parameters, V is the perturbed vulnerability index computed using a lower number of input layers, N and n are the number of parameters considered to calculate V and V , respectively. The average variation index of all cells for each parameter can be used to identify which parameter could have less influence on the outcome if omitted [72]. Any parameter with a lower overall average index value has less influence. A cell with a high or low index for a removed parameter will show a greater or smaller influence on the outcome, accordingly.

Single-Parameter Sensitivity Analysis
The single-parameter sensitivity analysis is a cell-based process, too [70]. The theoretical weight assigned by the SINTACS and SINTACS-LU methods to each parameter was compared to the "effective" or "actual" weight, which was calculated as follow: where P r and P w are the rating and weight of the P parameter, respectively, and V is the vulnerability index as computed in Equations (1) or (2). The goal of this analysis was to determine which parameter had the highest effective weight and therefore provided higher values of the vulnerability indexes when its spatial variation throughout the entire study area was considered [70,73].

Conditioning Indicators
5.1.1 The depth of the water table (S) appears to be one of the parameters that most significantly impacts groundwater vulnerability. Shallow aquifers are commonly more vulnerable than deep aquifers, where infiltrating contaminants require a longer interaction time with unsaturated media, and biogeochemical reactions can reduce their impact. The data on the groundwater level were acquired from 42 wells utilized for agricultural and domestic supply, and monitoring aims (Figure 1). These data were elaborated to determine water table depth (S). The observed data highlight that the water table is relatively shallow and the water depth ranges from 1.5 to 8 m. Kriging method was used for the interpolation of the punctual data. The depth of the water rasterized layer was ranked by assigning four classes of ratings, listed in Table 1, and its spatial distribution is shown in Figure 4a. It is evident that the shallow aquifer will result in a higher potential intrinsic vulnerability, whereas the deep aquifer is relatively less impacted by contamination. In the investigated area, the southeastern sector shows high intrinsic vulnerability due to the higher rating of this parameter (Figure 4a). the degradation phenomena currently underway, mainly in the study area characterized by unfavorable climatic conditions [69,76]. The land use impact on the vulnerability was taken into account by adding the LU parameter to the SINTACS method (Table 1).

SINTACS Vulnerability Index
The SINTACS and SINTACS-LU indexes were determined using Equations (1) and (2), and the resulting intrinsic vulnerability indexes are represented in Figure 5a,b.

Effective infiltration (I)
is the amount of water that infiltrates through the unsaturated zone and reaches the aquifer. The infiltration was estimated using the inverse hydrogeological water balance [24,74]. This methodology consists of the infiltration estimation starting from the analysis of climate data, such as rainfall and temperature, topographic and hydrogeological parameters, such as the altimetry and the hydrogeological complexes characteristics of the study area. The GIS-based method allows to obtain the spatial distribution of the effective infiltration rate in the hydrogeological basin [65,66]. In QGIS software, daily rainfall and air temperatures series were collected and elaborated for the period 2000-2015. In the hydrogeological basin, the mean annual rainfall is about 538 mm. The results of the GIS-based distributed inverse hydrogeological water balance procedure show that the actual evapotranspiration (equal to 465 mm) provides an effective infiltration rate equal to 52 mm/year, which is equivalent to about 1080 L/s. The mean annual value of direct runoff is equal to 21 mm/year [67] (Figure 4b).
5.1.3 The unsaturated zone (N) or vadose zone represents the thickness defined by the soil horizons and the groundwater table. The attenuation capacity of the contaminants is strongly affected by the unsaturated media properties that influence many processes such as biodegradation, diffusion and chemical reactions. The possibility of groundwater contamination can increase with small thicknesses of the vadose zone and high permeability, consequently enhancing the intrinsic vulnerability of the aquifer. The values of the parameter N were obtained from well stratigraphic profiles integrated with geological maps and geophysical surveys [60]. The N parameter ratings range from 2 to 7 ( Table 1). The results highlighted that among the vadose zone classes: Medium-fine alluvial deposits (rating 5), Fine alluvial deposits (rating 4), Coarse alluvial deposits (rating 6), Sand coastal deposits (rating 7) and Clay deposits (rating 2) cover around 42%, 26%, 20%, 11% and 1% of the total study area, respectively (Figure 4c).

Soil media (T)
is the uppermost vadose layer where biological activities are significant. Soil characteristics condition the behavior and the amount of the contaminants that can reach the water table [75]. The presence of organic components, the thickness of soil media, and the grain size can significantly impact the contaminant's attenuation mechanisms. Different textures of soil media (coarse sand, sand, sandy loam, loam, silty clay loam and clay loam) are present throughout the study area (Figure 4d). 5.1.5 Aquifer media (A) describe the components' characteristics of the saturated zone such as the porosity rate, particle types and grain size, and the processes that take place below the piezometric level. Basically, these processes are hydrodynamic dispersion, dilution, absorption and chemical reactions between the aquifer media and the contaminants [61]. The aquifers of the study area are mainly unconfined. The results showed that the medium-fine alluvial complex covers 59% of the study area. Alluvial, sandy and fine sand complexes cover the 29.5%, 11%, and 0.5%, respectively (Figure 4e). The analysis of the thematic layer suggests that the south-eastern portion of the study area (closest to the coast) is under high intrinsic vulnerability. 5.1.6 Hydraulic conductivity (C) is defined as the ability of the fluid to pass through the aquifer media. This physical property expresses the groundwater flow and mobile contaminants rate that can reach the water table [23]. High hydraulic conductivity gives higher potential groundwater contamination and, therefore, greater intrinsic vulnerability. The hydraulic conductivity of the study area ranges from 10 −5 to 10 −3 m/s (Table 1, Figure 4f). 5.1.7 Topography (S) is expressed by the variability of a land surface slope as a percentage. This study used DEM data for extracting the slope information. Low-slope areas can be more vulnerable to pollution as they retain water longer and contaminants are more likely to infiltrate. More than 70% of the study area has a slope class of 0-2%, corresponding to flat areas, which are given the highest rating and, therefore, have a significant impact on intrinsic vulnerability (Figure 4g). 5.1.8 Land use (LU) of the study area is constituted as urban, industrial and vegetation land cover [69]. The quality of groundwater resources in the Metaponto coastal plain is mainly threatened due to population growth, intensive agricultural, industrial and urban activities. The development of agricultural activities has resulted in a general increase in the size of cultivated areas. The vegetation covers areas largely characterized by agricultural crops. The plain is characterized by the presence of forests and pasture located mainly in the southeastern part. Along the coast, the Mediterranean maquis takes on the prevailing vegetational characteristics. The Mediterranean maquis is replaced by the garrigue, a low, soft-leaved scrubland, in the areas with unfavorable climatic and soil conditions and where anthropogenic pressures, such as fire and grazing, are stronger [69,76]. Some areas are characterized by parcels allocated to permanent crops, mainly vineyards, olive trees and sown land. Permanent crops are present on all plain; the coastal area is mainly covered by orchards. Inadequate or intensive farming practices have aggravated the degradation phenomena currently underway, mainly in the study area characterized by unfavorable climatic conditions [69,76]. The land use impact on the vulnerability was taken into account by adding the LU parameter to the SINTACS method (Table 1).

SINTACS Vulnerability Index
The SINTACS and SINTACS-LU indexes were determined using Equations (1) and (2), and the resulting intrinsic vulnerability indexes are represented in Figure 5a According to the procedure [3], the results of the SINTACS parameters elaboration show that the groundwater vulnerability index ranges from 62 to 185; it was classified into low (<105), moderate (106-140) and high (>140) vulnerability classes [3] (Figure 5a). The results even indicated that the low, moderate and high groundwater vulnerability zones cover 7%, 51% and 42% of the study area, respectively. The application highlights that more than half of the study area is at moderate vulnerability in terms of intrinsic susceptibility to pollution (Figure 6a). According to the procedure [3], the results of the SINTACS parameters elaboration show that the groundwater vulnerability index ranges from 62 to 185; it was classified into low (<105), moderate (106-140) and high (>140) vulnerability classes [3] (Figure 5a). The results even indicated that the low, moderate and high groundwater vulnerability zones cover 7%, 51% and 42% of the study area, respectively. The application highlights that more than half of the study area is at moderate vulnerability in terms of intrinsic susceptibility to pollution (Figure 6a).

Map Removal
In the map removal sensitivity analysis, each SINTACS and SINTACS-LU parameter was removed one at a time to see how it affected the output. The effect is expressed in terms of the variation index [72] calculated in each cell of the study area by Equation (3). Tables 3 and 4 give a statistical overview of the variation index resulting from the removal By a combination of the eight thematic layers in the GIS-based SINTACS-LU method, the intrinsic vulnerability index was calculated. The elaboration carried out considering the LU parameter shows that the vulnerability index ranges from 77 to 215 (Figure 5b).
Since the SINTACS and SINTACS-LU indexes do not vary in the same range, they have been normalized to allow the vulnerability zonation. The normalization of the SIN-TACS and SINTACS-LU indexes of the single cells was performed based on the following equation [3,26,44]: where I norm is the normalized index, and I min and I max are the minimum and maximum values of the indexes, respectively. The groundwater vulnerability map resulting in the application and normalization of the SINTACS-LU classifies the majority of the plain (56%) as moderate vulnerability. Seven percent of the study area is classified as low and the 37% of the study area as high vulnerability (Figure 6b).
SINTACS and SINTACS-LU results classified as a very high vulnerability class in any part of the study area.
Low vulnerability areas, situated at a small scale in the upper part of the hydrogeological basin, showed low ratings of the water depth and soil media, despite high ratings of the hydraulic conductivity and effective infiltration. Part of the northwestern and the central sides show diffuse moderate vulnerability. The southeastern side shows high vulnerability due to the effects of the shallow water table, the small thicknesses of the unsaturated zones and the soil media typologies.
The vulnerability classes of both methods pointed out comparable pattern. In the SINTACS-LU approach, a great land use impact was noticed due to the high weight assigned to it.
The use of fertilizers in agricultural activities can be recognized as a source of nitrate contamination in the study area. Hence, the consideration of the land use as a vulnerability conditioning parameter and the "Nitrates" string of weight are significantly representative of the actual land use condition and implies the effective suitability of SINTACS-LU for assessing groundwater vulnerability in the study area. The comparison between SINTACS and SINTACS-LU results, highlights that the SINTACS-LU method takes into account the specific characteristics of the land use that plays an essential role in groundwater vulnerability evaluation.
The southeastern part is classified as highly vulnerable on both maps. The lowest vulnerable portions, located in the north-western and southern zones of the study area, are similarly distributed in both maps. The intensive agricultural, industrial and other anthropogenic activities with the peculiar hydrogeological characteristics of the aquifer system favor the high vulnerability to contamination. In the SINTACS-LU map some areas, classified as highly vulnerable in the SINTACS method, show a minor vulnerability class. These areas are situated in natural and forest sectors of the Metaponto plain, which are less populated, and where human impact on the groundwater is limited.

Map Removal
In the map removal sensitivity analysis, each SINTACS and SINTACS-LU parameter was removed one at a time to see how it affected the output. The effect is expressed in terms of the variation index [72] calculated in each cell of the study area by Equation (3). Tables 3 and 4 give a statistical overview of the variation index resulting from the removal of each parameter. The mean of the variation index in all cells for each variable was used to evaluate the sensitivity of that variable if omitted. For the SINTACS method, the removal of the effective infiltration (I) parameter caused the largest variation index, with a mean of 2.99%. Aquifer media (A) and hydraulic conductivity (C), with mean values of 2.23% and 1.90%, respectively, were the other parameters with higher variation indexes. In descending order, the other parameters are the depth of the water table (S), unsaturated zone (N) and soil media (T). With a variation index of 0.78%, Topography (S) had the lowest means.
In the SINTACS-LU approach, the removal of the land use (LU) parameter caused the highest variation index, with a mean of 2.03%. In descending order, the other parameters are the effective infiltration (I), depth of the water table (S), aquifer media (A), hydraulic conductivity (C), soil media (T), and topography (S). The unsaturated zone (N), with a variation index of 0.27%, had the lowest mean. Tables 5 and 6 show for each variable the statistics of the single-parameter sensitivity analysis: the theoretical weight originally assigned by SINTACS and SINTACS-LU methods, the same weight normalized to 100, the average effective weight computed on the entire study area, the standard deviation and the minimum and maximum values. In both methods, almost all parameters showed an effective weight higher than the theoretical one. The effective infiltration (I) showed effective weight was significantly lower than theoretical one, whereas the unsaturated zone (N) and soil media (T) showed the lowest differences between the effective and theoretical weights.

Single-Parameter Sensitivity Analysis
For SINTACS method, the depth of the water table (S) mostly influenced the final vulnerability index, with an average weight of 27.17%, in agreement with the highest theoretical weight of 19.23%. On the contrary, the effective infiltration (I), which together with the depth of the water table (S) and soil media (T) had the highest theoretical weight, had the lowest effective weight with an average value of 5.17%. In the SINTACS-LU approach, the land use (LU) was the most effective parameter on the final vulnerability index (average weight 23.95%), coherently with the highest theoretical weight of 16.13%, in agreement with the results of the map removal sensitivity analysis. On the other hand, as for the SINTACS method, the effective infiltration (I), which together with the depth of the water table (S), soil media (T), and land use (LU) had the highest theoretical weight and had the lowest effective weight with an average value of 4.41%.

Conclusions
In the present paper, the overlay and index methods SINTACS and SINTACS-LU were applied for assessing the groundwater vulnerability of the Metaponto coastal plain. The land use parameter was added because of its importance in the study area and its impact on groundwater vulnerability. The groundwater vulnerability to pollution of the study area was classified into low, moderate and high classes.
The sector classified as highly vulnerable is closest to the sea and represents the area that most needs opportune management of groundwater resources.
Map removal and single-parameter were the two types of sensitivity analysis performed in this research to understand the influence of the input parameters on the output.
The assessment of the intrinsic vulnerability of the Metaponto coastal plain aquifer will assist in prioritizing groundwater protection measures and help in the selection of the most vulnerable areas for further investigation and monitoring. This is more important, considering that vulnerability assessment is gaining increasing attention because of the key role that groundwater quality plays in the study area.
The knowledge of potential intrinsic vulnerability and the issues related to the quality of groundwater represents a significant task to prevent and control the pollution of the resources. Indeed, the identification of the aquifer zones characterized by high vulnerability represents an important step useful for both groundwater management and communities' awareness concerning the resulting environmental problems.
Furthermore, this study leads to new considerations concerning the link of the aquifer hydrogeological behavior to the anthropogenic activities, related to the release of the pollutants in the subsurface, and climate change. For understanding these issues, future studies should consider the factors influencing the hydrogeological situations in responding to changing conditions, with particular attention to the effects of climate change on groundwater resources that could determine negative impacts in terms of both quantity and quality.
Author Contributions: Investigation, data analysis, methodology application and writing-original draft F.C. and R.M.; writing-review F.C., R.M. and F.S.; supervision F.C. and F.S. All authors have read and agreed to the published version of the manuscript.
Funding: This research received no external funding.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.