Geochemical, Geological and Groundwater Quality Characterization of a Complex Geological Framework: The Case Study of the Coreca Area (Calabria, South Italy)

: Hydrogeochemical characterization and statistical methods were used to investigate the groundwater quality and the origin of constituents (anthropic or natural) in groundwater of the Coreca area (Calabria, South Italy). Coreca is characterized by an articulated geological setting where the three main geological complexes that distinguish the Northern Calabria Peloritan Orogen (CPO) outcrop. This complex asset affects the quality of groundwater mainly exploited for irrigation use. In particular, the presence of ultramafic rocks (e.g., serpentinite and metabasite) promotes the release of harmful elements such as Cr and Ni. In the studied area, two groups of waters were identified: Ca-HCO 3 waters strongly controlled by the interaction with Ca-rich phases (e.g., limestone), and Mg-HCO 3 waters related to the interaction of meteoric water with the metamorphic units. Statistical elaboration allowed to detect, in the Mg-HCO 3 group, a good correlation between Cr and Ni (not observed in Ca waters) and a negative correlation between Cr, Ca and Al, in agreement with direct interaction with ultramafic rocks characterized by low concentrations of CaO and Al 2 O 3 . The concentration of major and trace elements has been compared with the Italian law limit values and the drinking water guidelines provided by the World Health Organization (WHO). Only three samples showed Mn and Ni concentration higher than the Italian law threshold. Furthermore, the assessment of groundwater quality was carried out using salinity and metal indexes. The groundwater quality assessment for irrigation allowed to classify the resource as “excellent to good” and “good to permissible”; nevertheless, a salinity problem and a magnesium hazard were found. Lastly, a metal index (MI) calculation revealed values <1 for almost all samples, pointing to good overall quality. Only a few samples showed a value extremely higher than 1, attributable to prolonged interaction with ultramafic rocks and/or localized anthropogenic pollution. From a general point of view, groundwater showed a generally good quality except for limited areas (and limited to the set of constituents analyzed) and a mild exceedance of the maximum salinity thresholds that must be monitored over time. Through a multidisciplinary approach, it was possible to ascertain the main anomalies attributable to the interaction with the hosting rocks and not (with few exceptions) to anthropic processes.

From a geological point of view, the Mesozoic Apennine Unit outcrops in a faultbounded tectonic window. The Apennine Unit consists of Triassic dolostone and dolomitic limestone (Verbicaro Unit) [42,43] and is overthrusted by the ophiolitic sequence (allochthonous Alpine Liguride) belonging to the Frido (mainly metapelites and slates) and the Gimigliano-Monte Reventino Units (serpentinites, metabasalts, phyllites, and carbonates) [52,53]. The metamorphic units are sealed by the Miocene sedimentary sequences consisting of calcareous sandstones, calcarenites, clays, marls, and Messinian limestone ("Calcare di base"). Pleistocene terraced deposits, consisting of conglomerates and sands, outcrop at the top of the succession.
The hydrogeology of the Coreca area is characterized by two different kinds of aquifers: fissured aquifers in the metamorphic units and partially in the Miocene deposits, and porous multilayer aquifers developed in the Miocene/Quaternary successions. Both types present relatively high volumes of groundwater storage and circulation, with the richest amount located in the fractured-altered superficial portion of the metamorphic units characterized by high permeability. Porous aquifers occur in gravelly-sandy permeable deposits of the quaternary succession, with flowrates variable according to annual rainfall. As reported by [54][55][56][57][58][59][60], Mesozoic limestone (Apennine Unit) represents the main thermal aquifer of the region and one of the main sources of drinking water supplies in Calabria and Southern Italy. However, in the Coreca area, Mesozoic successions do not highlight thermal evidence, with much less developed surface aquifers due to the high fracturing.

Methods
A total of 23 groundwater samples, from 2011 to 2014, were collected and analyzed for major cations, anions, and trace elements. Furthermore, 9 representative rock samples were collected and analyzed.
Water samples were collected from 16 springs and 7 wells (the location for each sample is reported in Figure 1).
Chemical-physical parameters such as pH, Eh, temperature, alkalinity, and specific electrical conductivity were measured in the field by means of portable instruments (HI-9828). Two pH buffers, with nominal pH values of 4.01 and 7.01 at 25 °C, were used for pH calibration. The ZoBell's solution [61] was used to calibrate the mV-meter for Eh measurement. Total alkalinity was determined by acidimetric titration, using HCl 0.05N as titrating agent and methylorange as indicator.
In the laboratory, the concentrations of F, Cl, Br, SO4, NO3, PO4, Na, K, Mg, and Ca were determined by HPLC (DIONEX DX 120). During the same day, it was measured dissolved reactive SiO2 by VIS spectrophotometry upon reaction with ammonium molybdate in acid media (and treatment with oxalic acid) to form a yellow silicomolybdate complex, whose absorbance was read at 410 nm [62]. Trace elements were analyzed by a quadrupole ICP-MS (Perkin Elmer/SCIEX, Elan DRCe) with a collision reaction cell capable of reducing or avoiding the formation of polyatomic spectral interferences. Data quality for major components was estimated by charge balance. Deviation between the sum of cation concentrations and the sum of anion concentrations, both in equivalent units, varied between −5% and +5%. Data quality for minor and trace elements was checked by running the NIST1643e standard reference solution. Deviations from the certified concentrations were found to be lower than 5%. The results of laboratory analyses and field data are shown in Tables 1 and 2. For each sample, the saturation index (SI), with respect to the mineral phases, was performed using PHREEQC Interactive software, version 3.1.1 [63] using the LLNL thermodynamic database.
The mineralogical associations for each main lithotype were determined using an optical microscope and by means of a Bruker D8 Advance XRD Diffractometer.

Mineralogical Characteristics
To characterize the main outcropping lithotypes, samples of metabasalts, serpentinites, and phyllites were collected for mineralogical analyses.
Optical microscopy observations and XRD analyses indicated that: (i) metabasalts, outcropping in the Coreca area, consisted mainly of chlorite, epidote, actinolite, and albite, with small amounts of calcite; (ii) serpentinites were composed of fine-grained crystals of serpentine minerals (antigorite as primary phase, as also reported by [64]), with subordinate magnetite; (iii) phyllites were made up of white mica and albite with significant amounts of chlorite and quartz, and small amounts of calcite.

Physical-Chemical Parameters
Springs (n = 16) and wells (n = 7) collected in the Coreca area ( Figure 1) showed an average temperature of 19.5 ± 2 °C and 20.5 ± 2.4 °C, respectively, which were slightly higher than the yearly mean atmospheric temperature (15.6 ± 5.3 °C). The pH and EC (electrical conductivity) values were comparable between springs and wells and showed average values of 7.47 ± 0.38 and 983 ± 149 (μS/cm) and 7.14 ± 0.17 and 964 ± 198 (μS/cm), respectively. Overall, Eh showed a wide positive range of values both for springs and wells, with only samples S25 and S39 highlighting negative values (Table 1).

Geochemical Characteristics
Water chemistry was investigated by using: (i) triangular plots involving major cations and major anions (Figure 2), both prepared starting from the concentrations in equivalent units; (ii) a correlation graph for SO4 vs. HCO3 + Cl, in which iso-salinity lines were drawn for reference ( Figure 3); (iii) box and whisker plots ( Figure 4); Triangular plots ( Figure 2) allowed to identify two groups of waters. The first group had a Ca-HCO3 composition suggesting a chemism controlled by dissolution of Ca-rich phases and/or other processes such as ionic exchange. Among the Ca-bearing phases, calcite, which occurs both in metabasalts and carbonatic rocks, represents the phase with the highest dissolution rates. Its presence promotes Ca-HCO3 waters [12,65,66]. The second group had a Mg(Ca)-HCO3 composition, probably due to the interaction with ultramafic rocks and/or carbonate-dolomitic successions, where antigorite, actinolite (in serpentinite and metabasalt), and dolomite are the local phases that can promote formation of Mg-HCO3 waters [5,12,67].   Triangular diagrams represent a useful tool to define the main geochemical groups. Unfortunately, no information about salinity was provided. To cope with this, samples were also classified using the TIS diagram (SO4 vs. HCO3 + Cl, Figure 3a), in which comparable salinity values between 12 meq/L to 24 meq/L are evident for the two groups of water. Salinity values and ratios of major constituents obtained in S22, S40, and S35 samples, belonging to Mg-HCO3 waters, suggested a prolonged water-rock interaction with the hosting units (see below and Table 3), increasing their salinity and Mg concentration. Furthermore, the S25 and S39 samples (Ca-deep wells) highlighted a considerable sulphate increase and high Fe (5733 and 7673, respectively) and Mn (112 and 76, respectively) concentration. These samples were representative of a third aquifer (Ca waters, Fe-rich), hosted in the calcareous sandstone directly in contact (tectonic contact) with the metamorphic basement (grey phyllites). This aquifer is isolated from the Ca-rich and Mg-rich systems due to the presence of the normal fault system (N-S), which put in contact (aquiclude) the Tortonian-Messinian filling to the east (the aquifer) with the metamorphic basement to the west. The anomalies were probably due to repeated alternations of reducing and oxidizing conditions that can promote dissolution of sulfides and Fe-Mn oxy-hydroxides.
With the aim to improve the knowledge about processes and evolution undergone by the considered systems, for each sample, the saturation index (SI), with respect to specific mineral phases, was performed using PHREEQC Interactive software, version 3.1.1 [63] using the LLNL thermodynamic database. Geochemical data were elaborated by using of the triangular plot for (HCO3 + CO3)-Mg-SiO2 reported in Figure 3b [68]. Unlike previous triangular diagrams, the triangular plot (HCO3 + CO3)-Mg-SiO2 in Figure 3b was prepared starting from the concentrations in weight units. The diagram allowed to compare the observed compositions with those expected for congruent dissolution of different magnesian minerals, such as serpentine, talc, sepiolite, brucite, and magnesite (and other carbonates); and for incongruent dissolution of Mg saponite and clinochlore (accompanied by precipitation of Al-secondary silicates). The compositions expected for dissolution of these solid phases were represented based on the stoichiometric coefficients of the relevant reactions.
This triangular plot is useful in areas with lithotypes characterized by Mg-bearing phases such as metabasites (chlorites) and serpentinites. With the diagram, an attempt was made to discriminate whether the Mg system belonged to the metabasite or serpentinite aquifer, or at least where the predominant circulation and interaction occurred.
Water-rock interactions with the specific phases are shown in Figure 3b. Mg-HCO3 waters fell between the compositions expected for dissolution of clinochlore and calcite, suggesting a prevailing interaction with rocks holding these two phases (e.g., metabasite). A shift toward phases linked to ultramaphic rocks was only mildly evident. Hydrogeological evidence confirmed the geochemical data that allowed the exclusion of Mg compositions linked to an interaction with the dolomitic successions (see Section 5). Ca-HCO3 fell above the dolomite-diopside line, suggesting an interaction with carbonate phases.
As reported in Table 3, most waters reached oversaturation with calcite, dolomite, clinochlore, tremolite, and albite, and reached values close to the saturation with phases characterizing the main outcropping lithotypes. Samples S22 and S35 showed the highest SI values with respect to clinochlore, albite, and tremolite, directly linked to metabasites, and phases characterizing ultramafic rocks (e.g., antigorite).

Statistical Analysis
Elements concentration and distribution were elaborated by a statistical approach based on box and whisker plots ( Figure 4) and Pearson's correlation coefficients (Tables 4  and 5).
As shown in Figure 4, Ca-HCO3 waters highlighted a relatively high concentration of trace elements (Sr, Mn, Fe and Ba), whereas Mg-HCO3 waters were characterized by relatively higher concentrations of elements such as Cr and Ni. The high concentrations of Fe and Mn, as previously stated, suggested a reducing condition established in a third system, outside the study area, directly connected with the phylladic basement, and where the mobilization of reduced species such as Fe 2+ and Mn 2+ is favored. High Sr concentrations were traced back to the dissolution of carbonates phases. For the Mg-HCO3 group, slight enrichments in Cr and Ni were probably linked to the slight interaction with ultramafic rocks. Moreover, to define the source of elements in waters and assess the expected relationships between elements by virtue of interactions with specific phases, correlations between major and trace elements for each group (Ca and Mg) were evaluated for the entire dataset (Tables 4 and 5) with a p-value of 0.05, including pH, temperature, EC, and Eh as additional parameters. The result of the correlation analysis showed a different constituent association for each group. In Ca-HCO3 waters, EC was strongly associated with Ca, Na, Cl, HCO3, and SO4, indicating high conductivity of groundwater due to the presence of these ions. The presence of evaporitic levels (e.g., gypsum) in the Messinian successions, in addition to CaCO3 as primary or secondary phase, can easily explain the high conductivities. In the Mg-HCO3 group, Mg replaced Ca, highlighting a direct role of Mg-bearing phases during water-rock interaction. Furthermore, the Mg group exhibited a good correlation between Cr and Ni (not observed in Ca waters), in agreement with previous studies carried out in areas with ultramaphic rocks in [11,[69][70][71][72][73][74] and a negative correlation between Cr, Ca, and Al, consistent with a direct interaction with ultramaphic rocks characterized by low concentrations of CaO and Al2O3 [75].
As reported by [11,69,70], for water interacting with ultramafic rocks, it is reasonable to expect a strong correlation between Mg, Cr, Ni, and Fe. Mg-HCO3 waters of the Coreca area showed, as previously highlighted, a strong positive correlation between Cr and Ni and a slight negative correlation between these and Fe and Mg. Poor correlation can be related to secondary processes (precipitation), which tend to vary the relative ratios between each element [76].
Mg vs. Ni, Cr, and Mn diagrams ( Figure 5) display the relationships between constituent directly linked to an interaction with ultramafic rocks. The diagrams show for all samples a general scatter of Ni and Cr more enriched in the Mg groups. In particular, the diagrams highlight a good correlation along the left side of each diagram-low salinity and low Mg concentration -. The positive correlation is lost on the right side -high salinity and high Mg concentration -where the samples are oversaturated in calcite and with the highest salinities. For these groups, calcite precipitation could play an important role in the final composition. In fact, as reported by [76], trace elements are incorporated into the growing carbonate crystal. The authors highlighted that calcite can effectively sequester a variety of toxic cations from solution (e.g., Cd, Pb, Zn, Cu, Mn, Co, Fe, etc.) and trap them into the solid phase [76][77][78].

Groundwater Flow Interpretation
The absence of deep geological and hydrogeological surveys resulted in a speculative reconstruction of the groundwater system of the area. The hypothetic groundwater table was outlined using geological sections (Figures 6 and 7). Well and spring data revealed a multi-aquifer system, suggested by the presence of springs and a deep well (productive at about 30 m) close to each other. The occurrence of a multi-aquifer system was also supported by water chemistry indicating the coexistence of Ca-HCO3 waters coming from a shallow Miocene-Quaternary aquifer, and Mg-HCO3 waters from a deeper metamorphic aquifer (metabasalts and serpentinites).
Schematic cross sections A-A' and B-B' (Figures 6 and 7) show the presence of Ca 2+enriched wells and springs both along the highest portions and the west side of the study area where the sedimentary Miocene succession outcrops. In the southern side, a permeability limit exists between the dolomitic (fractured and without springs) and the metamorphic rocks (with a high number of springs). Moreover, the schematic cross section B-B' (Figure 7) illustrates the multi-aquifer system and the differences between the west and east side of the studied area. A suspended Miocene carbonate aquifer occurs in the west side related to a permeability threshold with the metamorphic basement, whose flowrate is greater to the top of the area. On the east side, the metamorphic aquifer prevails, with the presence of Mg-HCO3 waters. Metamorphic aquifer is well developed in the altered and fractured shallow portions that favor the water circulation. The east side is characterized by a permeability barrier, probably occurring between the metamorphic complex (metabasalts and serpentinites) and the underlying phyllites, and/or at the contact with the Miocene carbonate deposits. As reported in Section 4.3, the east downstream side (southern portion of the studied area) is pouring in Mg-HCO3 waters. In this portion, the Mg aquifer is confined to the east by a normal fault system (N-S) representing a physical barrier to fluid flow (sealing fault zone) between the Tortonian-Messinian filling and the metamorphic basement (underground watershed without racking). Furthermore, the Upper Tortonian silty clays and siltstones, interbedded with the Tortonian sandstones and calcarenites, could play a buffering role not favoring the mixing between the two aquifers characterized by different chemical compositions. Figure 6. Geological-hydrogeological schematic cross-section A-A' through the southern sector of the Coreca Area (modified after [79]). The trace of the section is shown in Figure 1.

Water Quality
Irrigation-water-quality parameters of the analyzed groundwater are shown in Table  6. Lastly, the general metal index (MI) was calculated to assess the water quality compared to metals in solution.

Salinity Hazard
The salinity hazard of irrigation water in terms of sodium adsorption ratio (SAR) is based on the relationship between Na ion and divalent cation [18]. It was determined using the following equation: Water with an SAR value <10 is considered excellent; 10-18 is good; 18-26 is fair; and above 26 is unsuitable for irrigation use [24]. In the study area, all the samples were excellent for irrigation purpose, with maximum values <2.5.
The suitability of groundwater for irrigation depends on the mineralization of water and its effect on plant and soil [80]. If the concentration of sodium is high in irrigation water, Na + tends to be absorbed by clay particles displacing Mg 2+ and Ca 2+ ions, thereby reducing soil permeability [80]. Sodium hazard based on SSP (soluble sodium percentage-SSP%) is useful in characterizing a water, since a high value indicates a soft water, and a low value indicates a hard water. It was calculated using the following equation: SSP = [(Na + + K + ) × 100]/(Ca 2+ + Mg 2+ + Na + + K + ) If SSP is <20, the groundwater is excellent for irrigation; between 20 and 40 is good; between 40 and 60 is permissible; between 60 and 80 is doubtful and may be dangerous; and >80 is unsuitable for irrigation. The calculated values of SSP of the groundwater samples indicated that 100% of the water samples were excellent for irrigation purposes.
The SSP vs. EC plot [24] provides a method for rating irrigation water. In Figure 8a, the Coreca groundwater fell in the "excellent to good" and "good to permissible" fields. However, in the SAR vs. EC diagram (Figure 8b, [20]), the Coreca groundwater was ranked as C2-S1 and C3-S1 (medium to high salinity hazard and low sodium hazard). Excess salinity reduces the osmotic activity of plants, limiting the absorption of water and nutrients from the soil of exchangeable sodium [81].  [20,24]).
Ca 2+ and Mg 2+ ions maintain a state of equilibrium in most groundwater [82]. In equilibrium, Mg 2+ in water affects the soil by making it alkaline, affecting the final crop yield [83]. The measure of the effect of magnesium in irrigated water is expressed as the magnesium adsorption ratio (MAR). [84] developed an index for calculating the magnesium hazard. MAR is calculated using the formula: If the MAR value is <50, then the groundwater may be used for irrigation purposes; on the contrary, if it is >50, it is unsuitable. High magnesium content in groundwater induces an additional alkalizing effect. Usually, Ca 2+ and Mg 2+ are in equilibrium in most waters, [20,85] suggested that MAR values exceeding 50% indicate a magnesium hazard as the soil becomes more alkaline, resulting in decrease in the availability of phosphorous [86,87]. The results of the MAR calculation showed that 100% of the Ca-HCO3 waters were suitable for irrigation purposes, whereas the Mg-HCO3 waters highlighted values of between 50 and 60, indicating a magnesium hazard.
The Kelly ratio (KR) is an important parameter formulated by [19] based on the level of Na against Ca and Mg. It is expressed as: Groundwater with a Kelly ratio <1.0 is considered suitable for irrigation purposes, while a KR value between 1 and 2 is classified as marginal and potentially dangerous for the groundwater, and if the KR is >2, then it is unsuitable. The KR value of the investigated groundwater showed that 100% of the samples were suitable for irrigation.
Lastly, the solubility of salts can be used to determine the potential salinity (PS). It is based on the sum of Cland half of SO4 2concentrations [22], and is expressed as: The Mg-HCO3 and Ca-HCO3 waters displayed PS values below 3, with only 3 cases (2 in Ca waters and 1 in Mg waters) with values between 3 and 4.

Metal Index (MI)
The metal index (MI) was used to assess the influence of overall pollution and illustrate the spatial distribution of heavy-metal concentration and the pollution index in the groundwater of the studied area [88].
The MI for drinking water was defined as: where C is the concentration of each element in solution and i is the ith sample. For elements with concentration < detection limit (d.l.), C = d.l. was assumed. Therefore, MI was defined to evaluate the quality of the water based on the content of 13 metals. The higher the concentration of a metal compared to its respective maximum acceptable concentration (MAC) value, the worse the quality of the water. If the concentration of a specific element is higher than the respective MAC value (MI > 1), the water cannot be used according to law. The presence of several elements with concentrations smaller than but close to the respective MAC values will also decrease the overall quality of water because of an additive effect. Thus, an MI value >1 is a threshold of warning. Almost all samples showed MI values <1, pointing to a good general quality of the resource ( Table 2). Only S18, S25, S39, S20, S36, and S9 showed values of 1.03, 117.6, 155.2, 3.73, 1.06, and 1.2, respectively. S18, S36, and S9 had values close to 1 due to an overall high metal concentration relative to the entire dataset (additive effect). S20, S25, and S39 displayed a very high MI value, owing to the high concentration of a few specific constituents (Mn, Fe, and Ni). S25 and S39 were collected in deep wells where repeated alternations of reducing and oxidizing conditions can promote the release of Fe and Mn in solution. High Ni concentration detected in S20 (spring with MgHCO3 composition) could be attributable both to prolonged interaction with ultramaphic rocks and/or localized anthropogenic pollution.

Conclusions
In the current study, hydrogeochemical characterization and statistical methods were used to investigate the groundwater quality and the origin of constituents (anthropic or natural) in groundwater of the Coreca area (Calabria, South Italy). Rock and water compositions were elaborated following statistical methods combined with hydrogeochemical modeling and conventional plots to investigate groundwater and related geochemical processes.
The mineral assemblage of the outcropping rocks and the multidisciplinary approach allowed the reconstruction of the water-rock interaction processes responsible for groundwater composition. Geochemical data and hydrogeological evidence confirmed the existence of two groups of groundwater: (a) Ca-HCO3, hosted in the shallow Miocene/Quaternary aquifer, and (b) Mg-HCO3 localized in the ultramaphic aquifer (serpentinites and metabasites). Calcite on the one hand and antigorite, tremolite, and clinochlore on the other represented the main phases that favored the formation of Ca and Mg systems.
The water-rock interaction with specific phases was confirmed by calculation of saturation indexes (SI), performed using PHREEQC Interactive software, and statistic elaboration. SI confirmed oversaturation with calcite, clinochlore, tremolite, and albite, and saturation with other phases characterizing the main outcropping lithotypes. Statistical approach allowed us to define, in Ca-HCO3 waters, high concentrations of Sr, Mn, Fe, and Ba, with a strong correlation between EC and Ca, Na, Cl, HCO3, and SO4, indicating high conductivity of groundwater and direct interaction with carbonatic phases. In the Mg-HCO3 group high concentrations in trace elements such as Cr and Ni were highlighted, and Mg replaced Ca, underlining a direct role of Mg-bearing phases during water-rock interaction. Furthermore, the Mg group exhibited a good correlation between Cr and Ni (not observed in Ca waters) and a negative correlation between Cr, Ca, and Al, in agreement with a direct interaction with ultramaphic rocks (e.g., serpentinite) characterized by low concentration in CaO and Al2O3.
Data confirmed a multi-aquifer system, as shown in the schematic cross sections that illustrate the differences between the west and east side of the studied area. A suspended Miocene/Quaternary carbonate aquifer occurs in the west side, related to a permeability threshold with the metamorphic basement, whose flowrate is greater to the top of the area. The east side is characterized by a permeability threshold, probably occurring between the metamorphic complex (metabasalts and serpentinites) and the underlying phyllites, and/or at the contact with the Miocene carbonate deposits. On the east side, the metamorphic aquifer prevails with the presence of Mg-HCO3 waters, and it is well developed in the altered and fractured shallow portions that favor the circulation of water.
Subsequently, major and trace elements were compared with the Italian law limit values and the drinking water guidelines provided by World Health Organization (WHO, 2004). Only the S20, S25, and S39 samples showed Mn and Ni contents higher than the lowest threshold provided by Italian law (50 ppb and 20 ppb, respectively).
The salinity index SSP allowed us to classify the Coreca groundwater as "excellent to good" and "good to permissible"; nevertheless, a salinity problem and a magnesium hazard (east side) were found. This suggests that efforts, including leaching and proper drainage, are needed to control the salinity hazard, especially for those waters, representing a non-negligible part of the Coreca groundwater dataset, having an EC higher than 750 S/cm. Moreover, the presence of a magnesium hazard suggests a need for long-term magnesium monitoring in the area to plan local alkalinity mitigation policies.
Lastly, calculation of the metal index (MI) revealed values <1 for almost all samples, testifying to a good general quality of the resource. Only a few samples had values much higher than 1, attributable to prolonged interaction with ultramaphic rocks and/or localized anthropogenic pollution.
From a general point of view, the data highlighted a good quality of groundwaters in the studied area (limited to the set of constituents analyzed), except for a few localized points and for limited exceedances of the maximum salinity thresholds, which must be monitored over time. Through a multidisciplinary approach, it was possible to ascertain the main anomalies attributable to the interaction with the hosting rocks and not (with few exceptions) to anthropic processes. Data Availability Statement: The manuscript is a data self-contained article, whose results were obtained from the laboratory analysis, and the entire data is presented within the article. However, if any additional information is required, these are available from the corresponding author upon request to the e-mail giovanni.vespasiano@unical.it.