Next Article in Journal
Water, Sanitation, and Hygiene Vulnerability among Rural Areas and Small Towns in South Africa: Exploring the Role of Climate Change, Marginalization, and Inequality
Next Article in Special Issue
Groundwater Vulnerability Analysis of Tirnavos Basin, Central Greece: An Application of RIVA Method
Previous Article in Journal
Discharge Estimation with the Use of Unmanned Aerial Vehicles (UAVs) and Hydraulic Methods in Shallow Rivers
Previous Article in Special Issue
Aquifer Parameters Estimation from Natural Groundwater Level Fluctuations at the Mexican Wine-Producing Region Guadalupe Valley, BC
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrogeochemical Processes and Natural Background Levels of Chromium in an Ultramafic Environment. The Case Study of Vermio Mountain, Western Macedonia, Greece

1
School of Mining and Metallurgical Engineering, Division of Geo-Sciences, National Technical University of Athens, 9 Heroon Polytechniou St., 15773 Zografou, Greece
2
Researcher, Orologa 8, 11521 Athens, Greece
*
Author to whom correspondence should be addressed.
Water 2021, 13(20), 2809; https://doi.org/10.3390/w13202809
Submission received: 2 August 2021 / Revised: 27 September 2021 / Accepted: 2 October 2021 / Published: 9 October 2021

Abstract

:
The hydrogeochemical processes and natural background levels (NBLs) of chromium in the ultramafic environment of Vermio Mountain, Western Macedonia, Greece, were studied. Seventy groundwater samples were collected from 15 natural springs between 2014–2020, and an extensive set of physical and chemical parameters were determined. The ultramafic-dominated environment of western Vermio Mt. favors elevated groundwater concentrations of dissolved magnesium (Mg2+), silicon (Si), nickel (Ni), and Cr in natural spring waters. Chromium was the principal environmental parameter that exhibited a wide range of concentrations, from 0.5 to 131.5 μg/L, systematically exceeding the permissible limit of 50 μg/L for drinking water. Statistical evaluation of hydrogeological, hydrochemical, and hydrological data highlighted the water-ultramafic rock process as the predominant contributor of Cr in groundwater. The NBL assessment for Cr and Cr(VI) was successfully applied to the typical ultramafic-dominated spring “Potistis” that satisfied all the methodology criteria. The NBLs of Cr and Cr(VI) were defined at 130 μg/L and 100 μg/L, respectively, revealing that a natural ultramafic-dominated environment exhibits the geochemical potential to contribute very high concentrations of geogenic Cr to groundwater. The holistic methodology, proposed herein, could be implemented in any catchment scale to assess geogenic and anthropogenic Cr-sources that degrade groundwater quality.

1. Introduction

Natural background levels (NBLs) are defined as “the concentration of a substance or the value of an indicator in a body of groundwater corresponding to no, or only very minor, anthropogenic alterations to undisturbed conditions” according to the Groundwater Daughter Directive (GDD) (Directive 2006/118/EC) [1]. Broadly, the term of NBLs is synonymous with the terms of environmental geochemistry “natural/geochemical background values” or “geochemical baseline” used in the past [2]. The NBLs are a set of several varying hydrogeological (i.e., the residence time of groundwater in the saturated zone, recharge by precipitation, hydraulic connection with other aquifer systems) [3,4,5], and hydrogeochemical (i.e., water–rock interaction, pH/redox conditions, chemical, and biological processes in the unsaturated zone) [5,6,7] factors. The determination of NBLs requires in-depth knowledge of geological/hydrogeochemical processes [8] and the distinguishment of natural and anthropogenic factors that affect the groundwater systems [9]. The need to separate NBLs from the anthropogenic impacts (e.g., urbanization, industrialization, agricultural activity) is frequently satisfied through statistical and pre-selection (PS) methods [10]. Such methods were applied within the EU-Specific Targeted Research Project BRIDGE (Background cRiteria for the iDentification of Groundwater thrEsholds), the objective of which was to develop a comprehensive methodology to evaluate threshold values (TVs) and NBLs of various qualitative parameters in the groundwater resources [11]. The first stage of this approach includes the PS method, which assumes that the groundwater samples represent pristine groundwater not affected by anthropogenic pressures [9]. It constitutes the most frequent method to exclude samples influenced by anthropogenic activities based on specific criteria such as concentrations of Cl, Na+, NO3, NH4+, and DO [6,10,12]. The PS method has been successfully applied to establish NBLs for different physical and chemical parameters, including EC, Cl, SO42−, F, As, Cr, Cr(VI), Mn, Ni, Fe, and V in many European water bodies [4,5,7,13,14,15,16,17]. The next stage contains statistical tools such as box plots and normality tests for assessing the NBLs of the target chemical parameter. An approach that incorporates both PS and statistical methods has been performed by many researchers [7,8,9,10,14,15], providing a comprehensive methodology to boost the validity of the assessment, mainly when the geochemical and geological features are adequately considered [16]. Thus, the challenging assessment of NBLs in an environment in which the prevailing geochemical conditions favor the occurrence and mobilization of naturally occurring chemical elements could provide essential information regarding the controversial geogenic and anthropogenic inputs in a complex environmental setting.
The water–ultramafic rock interaction is of great scientific interest due to the high content of the latter in Cr (1000–3000 mg/kg), and other potentially toxic elements (PTEs) such as As, Co, Fe, Mn, and Ni compared to the Earth’s crust composition [8,9,10,11,12,13,14,15,16,17,18,19,20,21] and to other rock types [22]; it constitutes the principal geogenic source of Cr in the environment [18]. Chromium is mainly hosted in spinels (e.g., chromite and magnetite) and silicates (e.g., pyroxene, serpentine, chlorite, olivine, talc). Serpentine group minerals can be highly enriched in Cr because it substitutes for magnesium (Mg) and/or iron (Fe) [18]. In the crystal lattice of most minerals, Cr occurs in the trivalent valence state [Cr(III)]. However, the geochemically immobile Cr(III) is oxidized into the mobile and toxic for the living organisms hexavalent chromium [Cr(VI)] in the presence of natural manganese oxides (MnO2), specifically pyrolusite (b-MnO2), in the typical range of groundwater pH (6.5–8.5) and under oxidizing redox potential (Eh) conditions [23,24,25,26,27,28]. Although an increasing number of studies focus on the occurrence and fate of Cr in the environment [29,30,31,32,33,34,35,36,37,38] only a few have systematically examined the geochemical fingerprint of water–ultramafic rock interaction in natural springs [29,33,35,39]. Typical worldwide examples of ultramafic springs with elevated groundwater concentrations of Cr(VI) have been recorded in the Province of La Spezia, Italy (up to 73 μg/L) [29], the Pollino massif, Italy (up to 30 μg/L) [40], the Gerania springs, Greece (up to 17.2 μg/L) [33], the Euboea Island, Greece (up to 37 μg/L) [41], and the Lesvos and Rhodes Islands, Greece (10–15 μg/L) [42]. The water–rock interaction constitutes a crucial and controlling factor concerning groundwater evolution. The geochemical reactions between the recharging water and the minerals of the host rocks affect the groundwater quality [43,44]. Hydrogeological and hydrogeochemical conditions such as pH, Eh, dissolved oxygen (DO), and groundwater flow path play a significant role in elevated concentrations of PTEs, including Cr, in the aquifer systems. Geochemical reactions such as ion exchange, weathering, precipitation/dissolution, and sorption process control the groundwater’s composition considerably. During chemical weathering, some major ions, PTEs, and other trace elements become mobile and release from the parent rocks to the groundwater along the flow path. In addition, the mobility and solubility of these elements are controlled by water–rock contact time, Eh–pH conditions, and chemical reactions with organic matter [36]. Ionic ratios, saturation indices (SIs), and geochemical bivariate plots are usually evaluated to determine the intensity of water–rock interaction and chemical reactions [43]. Hence, the primary target of studying the mechanism of water–rock interaction is to elucidate the indissoluble association between the geological environment and the qualitative characteristics of groundwater.
Ιn this work, we study the geochemical fingerprint of the water–ultramafic rock interaction process in the western Vermio Mt., Western Macedonia, Greece, and determine the NBLs of Cr in groundwater from natural springs. At the catchment scale of the Sarigkiol Basin, elevated groundwater concentrations of Cr (up to ~140 μg/L) have been recorded in irrigation wells in the lowland [45]. Based on geospatial and multivariate statistical analyses of data from selected natural springs, irrigation wells, and surface waters the increased concentrations of Cr were attributed mainly to geogenic origin with the synergistic contribution of anthropogenic factors [45]. Challenged by the leaching potential of Cr of the ultramafic rocks in the area, we focus, herein, exclusively on the natural springs located in the ultramafic environment of western Vermio Mt, assessing hydrogeochemical data of a 7-year monitoring period (2014–2020). The springs are ideal for setting the NBLs at the catchment scale of the Sarigkiol Basin, because: (a) they record a strong ultramafic footprint, (b) they are located at a high altitude (>1300 m), (c) they exhibit unique worldwide high to very high concentrations of Cr (up to ~130 μg/L) [45], and (d) they are not affected by anthropogenic activities. Defining the NBLs of Cr in western Vermio Mt., will facilitate the identification of Cr origin in groundwater in the Sarigkiol Basin. This is the first systematic study of the natural springs of western Vermio Mt. and provides important hydrogeochemical data for the geogenic footprint of a natural ultramafic environment on the groundwater quality.

2. Materials and Methods

This section contains basic information about the: (i) study area, (ii) geological and hydrogeological setting, and (iii) sampling, chemical analysis procedures, and data processing.

2.1. Case Study

The present study is focused on the western Vermio Mt., located in the eastern part of the Sarigkiol Basin, Western Macedonia, Greece. The altitude of Vermio Mt is 2025 m and the average altitude of the basin is 650 m. The study area lies between the latitudes 40°25′00″ and 40°28′00″ E and the longitudes 21°56′00″ and 22°59′00″ N (Figure 1). In this area, any extensive anthropogenic activities lack except for local livestock farming and sporadically logging.

2.2. Geological and Hydrogeological Setting

The western Vermio Mt. is composed of (Figure 1 and Figure 2) [46,47]: (a) alluvial deposits, (b) clastic conglomerates, talus cones, and breccias, (c) upper Cretaceous flysch, (d) a complex of schists and cherts formations, (e) ultramafic rocks (serpentinites and peridotites), (f) Triassic–Jurassic limestones, (g) Cretaceous limestones.
The aquifer systems in Vermio Mt. are:
(a)
The deep karstic aquifer of the Triassic‒Jurassic limestones, which form the mountainous boundaries and the basement of the Sarigkiol Basin,
(b)
Perched aquifer systems that are developed in the highly fractured serpentinites of Vermio Mt. due to secondary porosity,
(c)
Small in size and capacity, karstic aquifers developed in the scattered Cretaceous limestones. There are many aquifers in which the water table varies from +700 up to +900 m. They are hydraulically connected and recharge the groundwater of the screes and talus cones in the ridges of the basin. The general flow direction of the groundwater is from the mountainous area to the center of the basin, i.e., NE–SW.
Different types of natural springs flow out in Vermio Mt. and specifically (Figure 1):
Contact springs formed where permeable formations (limestones, breccia, conglomerates) overlay formations of low permeability or impermeable (altered ultramafic rocks/serpentinites). Contact springs studied here were: the springs S19, S10, S15 in the Agio Pnevma area, the “Potistis”‒W13 spring, the spring S1 in the Agios Dimitrios area, the spring “Mouratidis”‒S2, the springs S13, S14 in the Agios Panteleimonas area, and the springs S5 and S6 in the Vazelona area.
Fault springs formed where impermeable rocks such as ultramafic rocks are in contact with an unconfined aquifer due to faulting; the springs “Elafakia”‒W14 and W21 in the Agios Panteleimonas area belong to this type.
The spring “Potistis”‒W13, presents great interest, because of the very high concentrations of Cr it exhibits [45]. The spring “Potistis”‒W13 flows out at an elevation of 1300 m in an ultramafic environment characterized by the absence of any anthropogenic activities. It constitutes a contact-type spring in the contact of conglomerates with ultramafic clastic material and limestones and ultramafic rocks. The aquifer, which discharges via the spring, flows through a weathered zone in serpentines. The natural recharge comes mainly from the seasonal precipitations via the permeable upper unsaturated zone (conglomerates, clastic material of ultramafic rocks, and limestones). The recharge water is mainly enriched with released PTEs (mainly Cr) from the ultramafic rocks, as the rainfall infiltrates through the weathered fractured ultramafic rocks. An additional lateral recharge takes place due to secondary porosity in the fractured ultramafic rocks. The high permeability of the unsaturated zone due to the presence of conglomerates in this area facilitates the direct recharge of the aquifer in a short time. The range of discharge was calculated from 205 L/h up to 1200 L/h, with an average value of 482 L/h. In Figure 3, the simplified hydrogeological section describes the natural recharge and the operation mechanism of the spring “Potistis”‒W13.
In western Vermio Mt. ultramafic rocks, mainly serpentinites, carbonates, schists and cherts occur [48]. The main mineral phases of the ultramafic rocks, depending on the degree of serpentinization, are serpentine [(Mg, Mn, Fe, Co, Ni)3−xSiO2O5(OH)4], olivine [(Mg,Fe2+)2(SiO4)]), pyroxene [(Mg,Fe2+)(Si,Al)2O6], talc [Mg3Si4O10(OH)2], chlorite [(Mg,Fe2+)5Al(Si3Al)O10(OH)8], tremolite [Ca2(Mg,Fe2+)5Si8(OH)2O22(OH)2], magnetite (Fe2+Fe23+O4) and Cr-rich magnetite [Fe2+(Fe3+,Cr)2O4] and chromite (FeCr2O4).

2.3. Sampling, Chemical Analyses, and Data Treatment

Springs to be systematically studied herein were selected based on their hydrogeochemical characteristics [45]. A total of 70 representative groundwater samples were collected from 15 natural springs during wet and dry periods from March 2014 to September 2020, following the groundwater sampling guidelines [49]. The 15 sampling sites (Figure 1) were classified into seven groups according to their location, lithology, and type: (i) twenty-three (23) samples were collected from the spring “Potistis”‒W13, (ii) seven (7) samples from the Agio Pnevma area (S18, S16, S19, S10, S17, S15), (iii) twenty-three (23) samples from the spring “Elafakia”‒W14, (iv) three (3) samples from the Agios Panteleimonas area (S13, S14, W21), (v) ten (10) samples from the spring “Mouratidis”‒S2 (vi) two (2) samples (S5 and S6) from the Vazelona area, and (vii) two (2) samples from the spring S1 in the Agios Dimitrios area. Considering that the number of the water samples differs between the seven groups, each group has been treated and evaluated separately (the statistical and geochemical analysis), so the analyses are classified as reliable.
The analytical methods for the determination of physical [i.e., temperature (T), pH, oxidation-reduction potential (ORP), DO, and electrical conductivity (EC)] and chemical parameters (i.e., major ions, PTEs, and other trace elements) are provided in detail in Papazotos et al. [33]. The calculations of Eh and the total dissolved solids (TDS) values were carried out by converting ORP measurements (i.e., adding 200 mV) and the summation of major ions in each collected water sample, respectively.
AquaChem 5.0 software was used to elaborate chemical analyses, develop a Piper diagram, and calculate alkalinity. The statistical analyses of the chemical data were performed with SPSS 22.0 software.

2.4. Spearman’s Rank Correlation Coefficient

Spearman’s rank correlation coefficient (ρ, also signified by rs) measures the strength and direction of association between two ranked variables, evaluating the degree of linear association or correlation between these independent variables. It presents many similarities to Pearson’s coefficient except that it operates on the ranks of the data rather than the raw data [50].
The Spearman’s rank correlation coefficient is calculated according to the following Equation (1) [51]:
r s = 6 i = 1 n d i 2 n ( n 2 1 )
where di = difference in paired ranks, n = number of cases, xi and yi = data pair.
The formula to use when there are tied ranks is Equation (2):
ρ = i ( x i x ) ( y i y ) i ( x i x ) 2 i ( y i y ) 2
The Spearman’s rank correlation coefficient, rs, can get values from −1 to +1. The equation for the calculation is developed so that it gives rs = +1 when the data pairs have a perfect positive correlation (di = 0) and rs = −1 for the perfect negative correlation, whereas rs = 0 indicates no association between ranks. The closer rs is to zero the weaker the association between the ranks is.
The values of the correlation coefficient are classified as very strong (0.80–1), strong (0.60–0.79), moderate (0.40–0.59), weak (0.20–0.39), and very weak (0.00–0.19) [52]. The correlation coefficient is highly statistically significant, marginally statistically significant when the p-value is p < 0.01, p < 0.05.

2.5. Shapiro‒Wilks Test

Shapiro–Wilks is a test of normality in frequentist statistics. The null hypothesis of this test is that the dataset is normally distributed. Thus, if the p-value is less than the chosen alpha level (0.05 in this case), then the null hypothesis is rejected and the data tested are not normally distributed. If the p-value is greater than the selected alpha level, then the null hypothesis cannot be rejected (Equation (3)) [53].
W = i = 1 n ( a i x ( i ) ) 2 i = 1 n ( x i x ) 2
where x(i) is the i-th largest order statistic, x- is the sample mean, and n is the number of observations.

2.6. Quantile–Quantile Plot

The quantile‒quantile (q–q) plot is a graphical tool for defining if two datasets come from populations with a common distribution [54], basically tests the conformity between the empirical distribution and the given theoretical one. On a Q–Q plot normally distributed data, the points in a Q–Q plot will fit on a straight diagonal line.

2.7. Geochemical modeling

The geochemical software PHREEQC version 3.1.2 [55] coupling with the MINTEQ database was used to calculate the saturation indices (SIs) of natural spring samples. Mineral SIs employed to define mineral dissolution and precipitation processes in the natural springs of western Vermio Mt. Saturation index is calculated by the Equation (4):
SI = Log IAP K sp = LogIAP LogK sp
where IAP = ion activity and Ksp = solubility product constant.
A positive SI indicates that the mineral is oversaturated or supersaturated with respect to the solution [56]; thus, the mineral could precipitate. Conversely, a negative SI indicates that the solution is undersaturated with respect to the selected mineral, suggesting that the mineral is dissolved in groundwater to reach equilibrium.

2.8. Calculation of NBLs of Cr

The assessment of NBLs for the target parameter was implemented based on the BRIDGE methodology [11]. The applied modified multi-method was separated into three steps: (a) the hydrogeochemical (bivariate plots, Piper, SI), (b) the PS method, and (c) the statistical analysis for estimating the NBLs (box plots for outliers, Q–Q (quantile–quantile) plots, and normality tests).
The applied methodology for the assessment of NBLs of Cr is described in detail in Figure 4. The pre-selection (PS) method, which is widely applied worldwide, was employed to select the suitable spring water samples for the NBLs assessment [10,16,57,58,59,60]. The PS method constitutes the methodology geochemical approach to validate the dataset according to similar geochemical characteristics and recognize the water samples that are affected by anthropogenic activities. In the first stage, the hydrochemical facies were selected based on the DO concentrations and Eh (ORP) [61]. The first dataset group included the water samples with ORP > 100 mV and DO > 3 mg/L. All water samples from the natural springs satisfied this criterion. The next criterion included consideration of redox conditions; if the prevailing conditions were oxidizing, then the concentrations of NO3 < 10 mg/L would be considered and if the conditions were reducing, then the NH4+ < 0.5 mg/L would be considered [9,61,62] to exclude the samples affected by anthropogenic activities [62,63]. Based on this criterion, the samples with NO3 > 10 mg/L were considered to be affected by anthropogenic activities and thus, were excluded from the new dataset. The next criterion required eight measurements per year for two years or two measurements per year for at least four years to exist for each spring [16].

2.9. Threshold Values (TVs) Derivation

The assessment of TVs was based on three scenarios [11,58] (Figure 5). The reference value was set equal to the water drinking acceptable limit (i.e., World Health Organization (WHO) guideline value).

2.10. Meteoric Genesis Index (MGI)

The meteoric genesis index (MGI) was also employed to classify the groundwater sources based on the depth of the meteoric water [64]. This index was calculated using the following Equation (5):
r 2 = N a + + K + C l S O 4 2 ,   all   concentrations   are   expressed   in   meq / L
when r2 < 1, the groundwater source is of deep meteoric water percolation type whereas when r2 > 1 the groundwater is of shallow meteoric water percolation type [65].

2.11. Meteorological Data

Daily and monthly rainfall data for the period 2014–2019 were evaluated from the meteorological station in the Ermakia village (40°30′325″ N, 21°51′233″ E) which is the most representative and the nearest one, located on Vermio Mt., at an elevation of 1100 m. The annual precipitation for the period 2014–2019 was estimated at 985 mm.

3. Results

3.1. Chemometric Analysis

All chemical analyses were grouped into the above-mentioned seven categories according to their location, lithology, and type to evaluate the results. The descriptive statistics (max, min, median) of the physical and chemical parameters, for the studied spring waters from March 2014 until October 2020, are summarized in Table 1 and Table 2. Some concentrations above the detection limit (DL) and below the quantification limit (QL) were not excluded because their very low concentrations do not affect the data processing of this work. All concentrations determined by ICP–MS presented a QL equal to the DL. The parameters NO2, NH4+, PO43−, Fe, Cd, Co were measured below the detection limit (BDL) in the majority of the samples. Silver, Au, Be, Bi, Cs, Pt, Re, Ga, Ge, Hf, Hg, In, Mo, Nb, Ta, Ti, Th, Tl, W, Zr were detected BDL in all water samples. Hence, these elements were excluded from the statistical analyses. In Table 3, the dataset for Cr-Cr(VI) and the geographical coordinates for the sampling points are given.
In the Agios Panteleimonas area, the natural spring “Elafakia”‒W14 and in the Agio Pnevma area, the natural spring “Potistis”‒W13 were examined separately from the other springs in these areas because of their elevated groundwater concentrations of Cr.
Box plots depicting the variation of eight groundwater quality parameters (pH, Eh, DO, EC, Ca2+, Mg2+, HCO3 and Si) are given in Figure 6 and Figure 7. The alkalinity of the water samples was calculated in a range of 1.95 × 10−3 meq/L (S18) up to 6.71 × 10−3 meq/L “Potistis”‒W13, with an average value of 4.60 × 10−3 meq/L.
The abundance of major ions varied significantly among the natural springs (Table 4). Concentrations of Cr and Cr(VI) exhibited a wide range of values (Table 1, Table 2 and Table 3, Figure 8, Figure 9, Figure 10 and Figure 11). Most physical and chemical parameters (e.g., EC, NO3, SO42−, As, B, Ba, Cd, Cu, Ni, Pb, Sb, Se, Zn, etc.) in the natural springs of Vermio Mt. exhibited lower concentrations than the desirable and permissible limits for drinking water. All samples were compared with the guideline of WHO for drinking water [66]. Exceedances were recorded in the natural springs “Potistis”‒W13, “Elafakia”‒W14, S1‒Agios Dimitrios area, and S18- Agio Pnevma area. Specifically, in the spring S18‒Agio Pnevma area, concentrations of K+ exceeded the guideline value of 12 mg/L for drinking water (WHO, [66]). Two exceedances of As, above the permissible limits of 10 μg/L, were recorded in the natural spring S1 of the Agios Dimitrios area. In the spring “Elafakia”‒W14 seven water samples exceeded the WHO guideline value of 50 μg/L for Cr concentration for drinking water [66], while in the spring “Potistis”‒W13 23 water samples exceeded this value. Finally, in total, three samples (one from the spring “Elafakia”‒W14, and two from the spring “Potistis”‒W13) exceeded the limit of 20 μg/L for the concentration of Ni for drinking water [66]. The systematic exceedances of Cr in the natural springs of Vermio Mt. were the principal reason for selecting Cr as the target parameter for NBLs assessment in this area.

3.2. Correlation Analysis of Water Samples

Spearman’s rank correlation coefficient was applied to selected parameters (Table 5); the rest of the measured parameters were excluded because most of their values were BDL. The most remarkable features of the Spearman’s rank correlation coefficients are the following: Cr presented a statistically significant (p < 0.01) very strong positive correlation coefficient with Cr(VI) (rs = 0.98), strong positive correlation coefficients with Mg2+ (rs = 0.76), Si (rs = 0.75), EC (rs = 0.71) and Ni (0.61), and moderate positive correlation coefficients with HCO3 (rs = 0.55) and alkalinity (rs = 0.55). Hexavalent chromium exhibited a statistically significant (p < 0.01) very strong positive correlation coefficient with Mg2+ (rs = 0.8), a strong positive correlation coefficient with Si (rs = 0.76) and Ni (rs = 0.67), and moderate positive correlation coefficients with HCO3 (rs = 0.59) and alkalinity (rs = 0.59). Magnesium exhibited a statistically significant (p < 0.01) strong positive correlation coefficients with EC (rs = 0.68), HCO3 (rs = 0.68), alkalinity (rs = 0.68) and Ni (rs = 0.60), while very strong correlation coefficient with Si (rs = 0.91). Bicarbonates presented a statistically significant (p < 0.01) strong positive correlation coefficient with EC (rs = 0.62), and moderate positive correlation coefficients with Ni (rs = 0.51), Si (rs = 0.48). Sulfates had a statistically significant (p < 0.01) very strong positive correlation coefficient with U (rs = 0.82), strong positive correlation coefficients with Sr (rs = 0.74), Br (rs = 0.70), Ba (rs = 0.63), and Ca2+ (rs = 0.67), and a moderate positive correlation coefficient with Na+ (rs = 0.52). Arsenic exhibited statistically significant (p < 0.01) strong positive correlation coefficients with U (rs = 0.72), Ba (rs = 0.67) and V (rs = 0.70), and a moderate positive correlation coefficient with K+ (rs = 0.49).

4. Discussion

4.1. Hydrogeochemical Characterization of the Natural Springs of Western Vermio Mt. the Ultramafic Fingerprint

The dominant hydrochemical types of the studied natural springs of western Vermio Mt. were Ca-Mg-HCO3 (40% of the water samples), Mg-Ca-HCO3 (33% of the water samples), and Ca-HCO3 (21%) (Figure 12). Other transitional water types in the study area comprised Ca-K-HCO3-Cl (3%) and Ca-HCO3-SO4 (3%). Based on the type and the geological environment of the springs, the Ca–HCO3 waters are considered to originate through the interaction of meteoric water with rocks containing Ca-bearing minerals, whereas water types enriched in Mg, were derived from the dissolution of ultramafic rocks [29]. The mixed Ca-Mg-HCO3 type indicated fresh recharge waters mainly related to carbonate rocks and less to ultramafic rocks. The Mg-Ca-HCO3 water type represents recharge waters related to Mg-rich rocks, suggesting the strong interaction with the ultramafic rocks of the area [67]. The springs “Potistis”‒W13, S1, and W21 that belong to this type are associated with fissured aquifers in strongly serpentinised ultramafic rocks, or they are in hydraulic connection with them.
The water–ultramafic rock interaction typically produces Mg-HCO3 water type [33,40,68] and slightly alkaline to strongly alkaline pH conditions [69] because of the absorption of dissolved CO2 from atmospheric water in the serpentine and olivine according to the Equations (6) and (7) [70]:
Mg 3 Si 2 O 5 ( OH ) 4 + 6 CO 2 + 5 H 2 O 3 Mg 2 + + 6 HCO 3 + 2 H 4 SiO 4
Mg 2 SiO 4 + 4 H 2 O + 4 CO 2 2 Mg 2 + + 4 HCO 3 + H 4 SiO 4
The pH values that characterised the studied springs varied from 7.3 up to 8.5, indicating near-neutral to slightly alkaline conditions that are typical of groundwater interacting with ultramafic and carbonate rocks [40,71]. Redox potential conditions were oxidizing up to strong oxidizing, as indicated by Eh, ranging from 300 mV up to 410 mV. The pH and Eh conditions in the studied springs favoured the release and solubility of the Cr oxyanion in groundwater since the solubility of oxyanions such as HCrO4, CrO42−, Cr2O72−, H2AsO4, and HAsO42− is enhanced with increasing pH [72].
Various tools are usually employed to define and evaluate the water–rock interaction processes in an area [73]. Bivariate plots of major ions and ionic ratios were used to study the hydrogeochemical evolution processes in the studied springs (Figure 13). In the bivariate plot of Ca vs. Mg, the water samples were grouped into three classes based on their Ca/Mg ratio (Figure 13a). In the first class belong the water samples with a Ca/Mg ratio below the 1:3 line. This class included two seasonal water samples from the spring “Potistis”‒W13 (W13_10_19 and W13_4_19) with a Mg-Ca-HCO3 water type, revealing that the flow path was mainly through serpentinites. The second class contained the water samples of the natural spring S1‒Agios Dimitrios area and “Potistis”‒W13. They are all of Mg-Ca-HCO3 type and characterized by Ca/Mg ratios plotted below the 1:1 and above the 1:3 lines; this suggests a mixture of Mg-HCO3 and Ca-HCO3, indicating that these waters were derived from interaction with serpentinites and Ca-rich rocks. The third class comprised the rest of natural springs, characterized by mixed water types and a Ca/Mg ratio above the 1:1 line, indicating a limited influence of serpentinites. The bivariate plot of Ca + Mg vs. HCO3 (Figure 13b) suggests an excess of (Ca + Mg) over HCO3 reflecting an additional non-carbonate source of Ca2+ and Mg2+ ions, such as the dissolution of silicate minerals [38,57]. Iron-Mg-silicates of ultramafic rocks, such as olivine, pyroxene, and amphibole are transformed to serpentine group minerals during the serpentinisation process. Dissolution reactions favour the Mg2+ and HCO3 release of the Mg-rich minerals (Equations (6) and (7)) [68].
During water–rock interaction, various chemical processes (e.g., fluctuation of ionic concentrations, mobilization of the dissolved components, and change in pH) are fingerprinted on the groundwater quality [74,75]. Gibbs diagrams are generally used to identify the hydrogeochemical evolution of groundwater, which involves precipitation, water–rock interaction, and evaporation–crystallization processes, based on TDS vs. Na+/(Na+ + Ca2+), and TDS vs. Cl/(Cl+ HCO3) scatter diagrams [76]. Herein, Gibbs diagrams were employed to assess hydrogeochemical processes that affect the water chemistry in the natural springs of western Vermio Mt. Figure 14 illustrates that all samples from natural springs fall into the water–rock interaction field, suggesting weathering of carbonate and silicate minerals. Although the use of Gibbs plots for groundwater has been disputed [77], the case study discussed herein exhibits none of the characteristics that could result in misuse of these plots (e.g., high SO42− concentrations, salinity sources, evolutionary flow paths, etc.). The implications of the Gibbs diagrams are in accordance with the calculated MGI index, according to which the waters from the natural springs are characterized as deep percolation types.

4.2. Hydrogeochemistry of Cr in Natural Ultramafic Springs

To further study the hydrogeochemistry of Cr in the studied springs, the average concentration of Cr was plotted vs. the water type of each spring (Figure 15). As shown, each water type is characterised by a wide range of concentrations of Cr, attributed to the different operation mechanisms of the spring and the weathering degree of the host geological formations. The mixed Mg-Ca-HCO3 water type ranges from very high concentrations of Cr, in the spring “Potistis”‒W13 (>100 μg/L), to much lower values (<20 μg/L) in the springs S1‒Agios Dimitrios area and W21‒Agios Panteleimonas area. The mixed Ca-Mg-HCO3 water type is related to a range of concentrations of Cr from 17 to 48 μg/L. On the other hand, all springs that are characterised by a Ca-HCO3 water type exhibit very low Cr concentrations (<5 μg/L) since mostly the carbonate rocks influence their hydrochemistry. In all springs, the dominant anion is HCO3, the principal source of which is the dissolution of carbonate and silicate minerals [33].
An interesting feature of the spring “Potistis”‒W13, derived from the evaluation of hydrogeochemical, hydrological, and meteorological data, is the decrease in concentrations of Cr in a very short time after rainfall; this is further supported by the strong linear regression of Cr vs. discharge (coefficient of determination R2 = 0.85) (Figure 16 and Figure 17). Low discharge results in increased water–ultramafic rock contact time and thus, in elevated concentrations of Cr.
In Figure 18a, the Mg/Si ratio vs. Cr in the springs is presented. The diagram is divided into three sub-groups according to the concentrations of Cr. The water samples with low concentrations of Cr (<30 μg/L) constituted 45.5% of the total samples, 73.3% of which exhibited a Mg/Si ratio lower than 2. Concentrations of Cr from 30 μg/L up to 50 μg/L, comprised 25% of the total water samples, 77.78% of which exhibited a Mg/Si ratio lower than 2.3; only the samples W13_6_18 and W13_7_18 which correspond to the lowest concentrations of Cr recorded in the spring “Potistis”‒W13 exhibited a Mg/Si ratio higher than 2.3. Of the total water samples, 42% exceeded the permissible limit of 50 μg/L for drinking water [66], with most of them corresponding to samples from the spring “Potistis”‒W13. Most samples presented a Mg/Si ratio higher than 2.3. Respective Mg/Si ratios have been reported for groundwater in other natural ultramafic environments [29,33].
The strong fingerprint of the water–rock interaction on the spring water chemistry and the geogenic origin of Cr in groundwater are indicated by the statistically significant very strong positive correlation coefficient of Cr with Si, the strong positive correlation coefficients of Cr with Mg2+, EC, and Ni, and the moderate positive correlation coefficients of Cr with HCO3 and alkalinity. Magnesium and alkalinity are two parameters usually increased with increasing degree of weathering; the latter has been reported to relate to elevated concentrations of Cr in groundwater [78]. Nickel is derived from the dissolution of Ni-bearing silicates which are released to groundwater under morphological and geochemical conditions that do not favour the occurrence of Fe-hydroxides and other secondary minerals capable of adsorbing Ni [79]. Nickel exhibited statistically significant, moderate positive correlation coefficients with Mg2+, EC, and Si, further highlighting its geogenic origin. The two natural springs with high concentrations of Ni (“Potistis”‒W13, “Elafakia”‒W14), also exhibit high mean concentrations of dissolved Si, and are of Mg-Ca-HCO3 and Ca-Mg-HCO3 water type. A similar case of high concentrations of Cr and Ni in Mg-HCO3 groundwater has been reported by Margiotta et al. [40]. Unlike the spring waters, Cr in the irrigational wells in the lowland of the Sarigkiol Basin was reported to strongly correlate with NO3 and P, indicating the synergistic role of the agricultural activities [45].
The statistically significant, very strong positive correlation coefficients of Cr with Cr(VI) (Spearman’s rank correlation coefficient rs = 0.986) was further proven by their linear regression with a strong linear relationship (coefficient of determination R2 = 0.99, Figure 18b). In the analysed water samples, the Cr(VI)/Cr ratio ranges from 20% up to 100%. Specifically: (a) 62.5–90.3% in the spring “Mouratidis”‒S2, (b) 80–99% in the spring S1- Agios Dimitrios area, (c) 20–62.5% in the springs at the Vazelonas area, (d) 39–100% in the springs at the Agio Pnevma area, (e) 40–99% in the spring “Elafakia”‒W14, (f) 76–100% in the spring “Potistis”‒W13). The fluctuation in the Cr(VI)/Cr ratio depends on the prevailing geochemical conditions (redox reactions, pH), the presence of iron or manganese oxides, and competing anions in each area, suggesting that various processes take place [78]. Hexavalent chromium is the principal form of Cr in the natural water springs (“Potistis”‒W13, S1‒Agios Dimitrios area, “Mouratidis”‒S2); several factors contribute to the high Cr(VI)/Cr ratio. Specifically, the geological environment, which is enriched in Ca and Mg-bearing minerals, enhances Cr(VI) to form complexes with Mg and Ca and inhibits Cr(VI) reduction [71]. The presence of manganese oxides enhances the Cr(III) oxidation to Cr(VI) in ultramafic rocks, soils, and unsaturated zone releasing Cr(VI) to groundwater (Equation (8)) [24].
Cr ( III ) + 1.5 δ MnO 2 ( s ) + H 2 O HCrO 4 + 1.5 Mn 2 + + H +
The pronounced role of minerals in the concentrations of Cr in natural springs was investigated via SIs of selected mineral phases present in the study area (Figure 19). The percentage distribution of the SIs for the selected mineral phases is given in Table 6 for all collected water samples.
The water samples from the natural springs in which concentrations of Cr exceeded 50 μg/L (“Potistis”‒W13 and “Elafakia”‒W14) were oversaturated in: (a) the carbonate mineral calcite (100%) and (b) the oxide minerals chromite (100%) and magnetite (100%). On the other hand, they are undersaturated in: (a) the serpentine group minerals lizardite (100%) and chrysotile (100%), (b) the pyroxenes enstatite (100%) and diopside (100%), and the amphibole tremolite (100%) and (c) the olivine (100%).
In general, the mineralogical phases that appear undersaturated tend to dissolve in water. The dissolution reactions contribute major, minor elements and PTEs to the groundwater. Chromium-bearing silicate minerals (serpentine, amphibole, pyroxene, chlorite, talc) occurred mostly undersaturated in the water samples, whereas Cr-rich oxides (chromite and Cr-magnetite) were oversaturated in the water samples. Therefore, silicate minerals are the principal geogenic contributors of Cr and other major/minor elements (e.g., Mg2+, Ca2+, HCO3, Si) and PTEs (e.g., As, Ni) to the spring waters of western Vermio Mt.

4.3. NBLs of Cr in the Ultramafic Environment of Vermio Mt.

The geochemical characteristics of the natural springs, the geological environment, and the water–ultramafic rock interaction are reflected in springs’ water quality. Chromium constitutes the principal environmental component in groundwater of the Sarigkiol Basin, originating primarily from geogenic and incidentally from anthropogenic sources [45]. This paper aims to assess the NBLs of Cr, which is of great interest in the catchment scale of the Sarigkiol Basin.
Based on the above-discussed hydrogeochemical data (e.g., pH, DO, Eh, Mg2+, Si, Cr, alkalinity, etc.), the most representative natural springs, which flow through and interact with ultramafic rocks, are the S1-Agios Dimitrios area, “Mouratidis”‒S2, “Potistis”‒W13, and “Elafakia”‒W14.
Take into consideration the modified methodology for assessing NBLs, the PS method was applied to create the new dataset. All samples from the natural springs satisfied the two criteria (ORP > 100 mV, DO > 3 mg/L and NO3 < 10 mg/L). Regarding the third criterion, the available time series of measurements, two natural springs, those of “Potistis”‒W13 and “Elafakia”‒W14, sufficiently satisfied this criterion. The resulting population was examined for the normality of the dataset with the Shapiro–Wilks test, a method proposed to be appropriate for a sample size less than 50 [80]. Although the number of the sampling sites is limited (2), they are considered representative because of the available time-series measurements, their hydrogeochemical characteristics, and the elevated concentrations of Cr, Si, Ni, and Mg2+.
The normality test of Shapiro–Wilks showed that there was no normally distributed population of the samples, either for Cr or Cr(VI) (p < 0.05). The outliers were identified via Box plots to exclude these measurements in the next step until the total elimination of the outliers, and then the population of the remaining data was rechecked for its normal distribution (Figure 20 and Figure 21). The last datasets of each parameter without outliers were double-checked for their normality with Q–Q plots and the Shapiro–Wilks test. In the spring “Elafakia”‒W14, the NBLs of Cr constitutes the 95th percentile (57.24 μg/L) of the population as it was not normally distributed (Figure 20) [9]. On the other hand, Cr(VI) dataset was normally distributed, and based on the methodology, the NBL was defined to be 51.20 μg/L (NBL = the maximum value of the normally distributed dataset). Similarly, in the spring “Potistis”‒W13, the NBLs of Cr is equal to the max value (130 μg/L) of the normally distributed dataset while the NBL for Cr(VI) was calculated at 100 μg/L (NBL = the 95th percentile of the dataset as non-normally distributed) (Figure 21).
The estimated NBL is higher than the REF value (NBL/REF > 1) for both natural springs of Vermio. According to this suggestion, the TVs are determined as NBL for Cr in the natural springs of “Elafakia”‒W14 and “Potistis”‒W13. The spring “Potistis”‒W13 is considered to be the most suitable one to define the NBLs in the ultramafic environment of western Vermio Mt., because its water type (Mg-Ca-HCO3) indicated a strong influence by ultramafic rocks, whereas the spring “Elafakia”‒W14, with a Ca-Mg-HCO3 water type, was mainly influenced by carbonate formations. In Figure 22, the flow chart describes the NBLs and TVs assessment for Cr and Cr(VI), for the natural spring “Potistis”‒W13.
The assessment of NBLs in ultramafic springs is a challenging modern methodology based on the continuous monitoring of hydrochemical parameters. Nevertheless, it is essential to mention that, since the environmental systems are complex and multicomponent, NBLs should not be treated as the absolute value above which a parameter is of anthropogenic origin; instead, NBLs constitute the minimum target value and the guide for investigating elevated groundwater concentrations of a hydrogeochemical parameter and elucidating the influence of anthropogenic factors in a study area on a larger scale (e.g., at a catchment scale). The continued monitoring of water quality parameters is likely to provide higher concentrations of the specific parameters in the future, and therefore subsequent recalculation of NBLs may lead to higher NBLs in the study area.
The high potential leaching of Cr in Vermio Mt., as derived from the above-mentioned calculated NBLs, is imprinted in the lowland of the Sarigkiol Basin [45]. The surface runoff and the discharge of springs enriched in Cr follow various flow paths via torrents or streams through the weathered mantle of ultramafic rocks in western Vermio Mt. and leach into the lowland of the Sarigkiol Basin (Figure 23). Due to the hydraulic connection between the western Vermio Mt. and the eastern part of the lowland of the Sarigkiol Basin, the defined NBLs apply to the latter, supporting the dominance of the geogenic factor in the high groundwater concentration of Cr in the Sarigkiol Basin.

5. Conclusions

The ultramafic-dominated environment of western Vermio Mt. was fingerprinted on the groundwater chemistry and specifically on the elevated concentrations of Mg2+, Si, Ni, and Cr in natural spring waters. Chromium was recognized as the principal environmental parameter in the natural spring waters of western Vermio Mt; in 42% of the studied spring water samples, the concentrations of Cr exceeded the WHO guideline value of 50 μg/L for drinking water. The geogenic origin of Cr in groundwater is recorded in the very strong positive correlation coefficients of Cr with Si, the strong positive correlation coefficients with Mg2+, EC, and Ni, and the moderate positive correlation coefficients with HCO3 and alkalinity.
The main factors that determined the concentration of Cr in the studied spring waters were:
(a)
the time response of the aquifers systems to precipitations; direct infiltration on the geological formation of the aquifer results in immediate recharge of it. As a consequence, quick contaminant dilution takes place and fluctuations in Cr concentrations are observed depending on the time response of the aquifers to precipitation,
(b)
the water–rock contact time; the longer the water–rock contact time is, the higher the Cr leaching is,
(c)
the flow path of groundwater; a flow path through weathered ultramafic rocks results in the enrichment of groundwater in Cr,
(d)
the degree of the serpentinisation of ultramafic rocks; the more serpentinised the ultramafic rocks are, the higher their leaching potential in Cr is, and
(e)
the prevailing geochemical processes that favor the oxidation of Cr(III) to the soluble and mobile Cr(VI), such as alkaline pH, oxidative environment, presence of manganese oxides.
The absence of anthropic/anthropogenic activities in western Vermio Mt., the sufficient time-series data, and the hydrochemical characteristics of the studied springs allowed the assessment of NBLs of Cr by applying a multi-method approach. Considering the hydrogeological, hydrochemical, and hydrological data in western Vermio Mt. and applying the PS method, the spring “Potistis”‒W13 was selected as the most representative one to define the NBLs of Cr in the area. The applied methodology is fully harmonized with the GDD and the Water Framework Directive (WFD, 2000/60/EC) [81]. The NBL of Cr was defined at 130 μg/L, and that of Cr(VI) at 100 μg/L. Based on the NBLs of Cr, TVs for Cr at a catchment scale, i.e., the Sarigkiol Basin, were defined to be equal to the NBLs. Conclusively, the ultramafic environment in western Vermio Mt. presents a high geochemical potential to dissolve and mobilize geogenic Cr.
This first systematic study of the natural springs of western Vermio Mt. provides important hydrogeochemical data for the geogenic footprint of a natural ultramafic environment on the groundwater quality. The proposed methodology could be implemented in any catchment scale aiming to distinguish between geogenic and anthropogenic groundwater deterioration and to establish new TVs considering the NBLs.

Author Contributions

Conceptualization, E.V. and M.P.; Data curation, E.V., M.P. and P.P.; Formal analysis, E.V., P.P. and M.P.; Investigation, E.V., P.P., M.P. and D.D.; Methodology, E.V., P.P. and M.P.; Project administration, M.P.; Supervision, M.P.; Visualization, E.V., P.P. and M.P.; Writing—original draft, E.V.; Writing—review and editing, E.V., P.P., M.P. and D.D. All authors have read and agreed to the summited version of the manuscript.

Funding

This research was carried out within the framework of the Research Project NTUA 623748.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Some data is contained within the article. All hydrogeochemical data may be available for collaborative research projects by specific agreements. For information, contact [email protected].

Acknowledgments

Special thanks to G. Stamatis, Emeritus Professor in Hydrogeology (Agricultural University of Athens) and A. Dimitriadis, from Agios Dimitrios Power Plant of Public Power Corporation for their contribution to the fieldwork. A. Stamos in EAGME is acknowledged for valuable discussion on the operation mechanisms of the springs in western Vermio Mt. We would like to thank the two anonymous reviewers for their constructive com-ments and suggestions that significantly improved the quality of the paper. Special thanks are ex-pressed to the editors for their careful editorial handling.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. European Union. EU Groundwater Directive 2006/118/EC. Off. J. Eur. Union 2006, L 372, 19–31. [Google Scholar]
  2. Amiri, V.; Nakhaei, M.; Lak, R.; Li, P. An integrated statistical-graphical approach for the appraisal of the natural background levels of some major ions and potentially toxic elements in the groundwater of Urmia aquifer, Iran. Environ. Earth Sci. 2021, 80, 432. [Google Scholar] [CrossRef]
  3. Edmunds, W.; Shand, P.; Hart, P.; Ward, R. The natural (baseline) quality of groundwater: A UK pilot study. Sci. Total. Environ. 2003, 310, 25–35. [Google Scholar] [CrossRef] [Green Version]
  4. Gemitzi, A. Evaluating the anthropogenic impacts on groundwaters: A methodology based on the determination of natural background levels and threshold values. Environ. Earth Sci. 2012, 67, 2223–2237. [Google Scholar] [CrossRef]
  5. Preziosi, E.; Giuliano, G.; Vivona, R. Natural background levels and threshold values derivation for naturally As, V and F rich groundwater bodies: A methodological case study in Central Italy. Environ. Earth Sci. 2009, 61, 885–897. [Google Scholar] [CrossRef]
  6. Urresti-Estala, B.; Carrasco-Cantos, F.; Pérez, I.V.; Gavilán, P.J. Determination of background levels on water quality of groundwater bodies: A methodological proposal applied to a Mediterranean River basin (Guadalhorce River, Málaga, Southern Spain). J. Environ. Manag. 2013, 117, 121–130. [Google Scholar] [CrossRef] [PubMed]
  7. Sacchi, E.; Bergamini, M.; Lazzari, E.; Musacchio, A.; Mor, J.-R.; Pugliaro, E. Natural Background Levels of Potentially Toxic Elements in Groundwater from a Former Asbestos Mine in Serpentinite (Balangero, North Italy). Water 2021, 13, 735. [Google Scholar] [CrossRef]
  8. Libera, N.D.; Fabbri, P.; Mason, L.; Piccinini, L.; Pola, M. A local natural background level concept to improve the natural background level: A case study on the drainage basin of the Venetian Lagoon in Northeastern Italy. Environ. Earth Sci. 2018, 77, 487. [Google Scholar] [CrossRef]
  9. Parrone, D.; Ghergo, S.; Preziosi, E. A multi-method approach for the assessment of natural background levels in groundwater. Sci. Total. Environ. 2018, 659, 884–894. [Google Scholar] [CrossRef] [PubMed]
  10. Preziosi, E.; Parrone, D.; del Bon, A.; Ghergo, S. Natural background level assessment in groundwaters: Probability plot versus pre-selection method. J. Geochem. Explor. 2014, 143, 43–53. [Google Scholar] [CrossRef]
  11. Muller, D.; Blum, A.; Hart, A.; Hookey, J.; Kunkel, R.; Scheidleder, A.; Tomlin, C.; Wendland, F. D18: Final Proposal for a Methodology to Set Up Groundwater Threshold Values in Europe; Background Criteria for the Identification of Groundwater Thresholds; Bridge Publications: Vienna, Austria, 2006. [Google Scholar]
  12. Wendland, F.; Berthold, G.; Blum, A.; Elsass, P.; Fritsche, J.-G.; Kunkel, R.; Wolter, R. Derivation of natural background levels and threshold values for groundwater bodies in the Upper Rhine Valley (France, Switzerland and Germany). Desalination 2008, 226, 160–168. [Google Scholar] [CrossRef]
  13. Hinsby, K.; de Melo, M.T.C.; Dahl, M. European case studies supporting the derivation of natural background levels and groundwater threshold values for the protection of dependent ecosystems and human health. Sci. Total Environ. 2008, 401, 1–20. [Google Scholar] [CrossRef]
  14. Ducci, D.; de Melo, M.T.C.; Preziosi, E.; Sellerino, M.; Parrone, D.; Ribeiro, L. Combining natural background levels (NBLs) assessment with indicator kriging analysis to improve groundwater quality data interpretation and management. Sci. Total Environ. 2016, 569–570, 569–584. [Google Scholar] [CrossRef] [PubMed]
  15. Biddau, R.; Cidu, R.; Lorrai, M.; Mulas, M. Assessing background values of chloride, sulfate and fluoride in groundwater: A geochemical-statistical approach at a regional scale. J. Geochem. Explor. 2017, 181, 243–255. [Google Scholar] [CrossRef]
  16. Masciale, R.; Amalfitano, S.; Frollini, E.; Ghergo, S.; Melita, M.; Parrone, D.; Preziosi, E.; Vurro, M.; Zoppini, A.; Passarella, G. Assessing Natural Background Levels in the Groundwater Bodies of the Apulia Region (Southern Italy). Water 2021, 13, 958. [Google Scholar] [CrossRef]
  17. Filippini, M.; Zanotti, C.; Bonomi, T.; Sacchetti, V.; Amorosi, A.; Dinelli, E.; Rotiroti, M. Deriving Natural Background Levels of Arsenic at the Meso-Scale Using Site-Specific Datasets: An Unorthodox Method. Water 2021, 13, 452. [Google Scholar] [CrossRef]
  18. Oze, C.; Fendorf, S.; Bird, D.K.; Coleman, R.G. Chromium geochemistry in serpentinized ultramafic rocks and serpentine soils from the Franciscan complex of California. Am. J. Sci. 2004, 304, 67–101. [Google Scholar] [CrossRef]
  19. Tashakor, M.; Modabberi, S.; van der Ent, A.; Echevarria, G. Impacts of ultramafic outcrops in Peninsular Malaysia and Sabah on soil and water quality. Environ. Monit. Assess. 2018, 190, 333. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Kelepertzis, E.; Galanos, E.; Mitsis, I. Origin, mineral speciation and geochemical baseline mapping of Ni and Cr in agricultural topsoils of Thiva Valley (central Greece). J. Geochem. Explor. 2013, 125, 56–68. [Google Scholar] [CrossRef]
  21. Ryan, P.C.; Kim, J.; Wall, A.J.; Moen, J.C.; Corenthal, L.G.; Chow, D.R.; Sullivan, C.M.; Bright, K.S. Ultramafic-derived arsenic in a fractured bedrock aquifer. Appl. Geochem. 2011, 26, 444–457. [Google Scholar] [CrossRef]
  22. Nriagu, J.; Nieboer, E. Chromium in the Natural and Human Environments; John Wiley & Sons: New York, NY, USA, 1988; Volume 20. [Google Scholar]
  23. Rai, D.; Sass, B.M.; Moore, D.A. Chromium(III) hydrolysis constants and solubility of chromium(III) hydroxide. Inorg. Chem. 1987, 26, 345–349. [Google Scholar] [CrossRef]
  24. Sperling, M.; Xu, S.; Welz, B. Determination of chromium(III) and chromium(VI) in water using flow injection on-line pre-concentration with selective adsorption on activated alumina and flame atomic absorption spectrometric detection. Anal. Chem. 1992, 64, 3101–3108. [Google Scholar] [CrossRef]
  25. Kotaś, J.; Stasicka, Z. Chromium occurrence in the environment and methods of its speciation. Environ. Pollut. 2000, 107, 263–283. [Google Scholar] [CrossRef]
  26. Berna, E.C.; Johnson, T.M.; Makdisi, R.S.; Basu, A. Cr Stable Isotopes As Indicators of Cr(VI) Reduction in Groundwater: A Detailed Time-Series Study of a Point-Source Plume. Environ. Sci. Technol. 2009, 44, 1043–1048. [Google Scholar] [CrossRef] [PubMed]
  27. Johnson, C.; Xyla, A.G. The oxidation of chromium(III) to chromium(VI) on the surface of manganite (γ-MnOOH). Geochim. Cosmochim. Acta 1991, 55, 2861–2866. [Google Scholar] [CrossRef]
  28. Fendorf, S.E.; Fendorf, M.; Sparks, D.L.; Gronsky, R. Inhibitory mechanisms of Cr(III) oxidation by δ-MnO2. J. Colloid Interface Sci. 1992, 153, 37–54. [Google Scholar] [CrossRef]
  29. Fantoni, D.; Brozzo, G.; Canepa, M.; Cipolli, F.; Marini, L.; Ottonello, G.; Zuccolini, M. Natural hexavalent chromium in groundwaters interacting with ophiolitic rocks. Environ. Earth Sci. 2002, 42, 871–882. [Google Scholar] [CrossRef]
  30. Tziritis, E.; Kelepertzis, E.; Korres, G.; Perivolaris, D.; Repani, S. Hexavalent Chromium Contamination in Groundwaters of Thiva Basin, Central Greece. Bull. Environ. Contam. Toxicol. 2012, 89, 1073–1077. [Google Scholar] [CrossRef] [PubMed]
  31. Dermatas, D.; Mpouras, T.; Chrysochoou, M.; Panagiotakis, I.; Vatseris, C.; Linardos, N.; Theologou, E.; Boboti, N.; Xenidis, A.; Papassiopi, N.; et al. Origin and concentration profile of chromium in a Greek aquifer. J. Hazard. Mater. 2015, 281, 35–46. [Google Scholar] [CrossRef] [PubMed]
  32. Hausladen, D.M.; Alexander-Ozinskas, A.; McClain, C.N.; Fendorf, S. Hexavalent Chromium Sources and Distribution in California Groundwater. Environ. Sci. Technol. 2018, 52, 8242–8251. [Google Scholar] [CrossRef]
  33. Papazotos, P.; Vasileiou, E.; Perraki, M. Elevated groundwater concentrations of arsenic and chromium in ultramafic envi-ronments controlled by seawater intrusion, the nitrogen cycle, and anthropogenic activities: The case of the Gerania Mountains, NE Peloponnese, Greece. Appl. Geochem. 2020, 121, 104697. [Google Scholar] [CrossRef]
  34. Coyte, R.; McKinley, K.; Jiang, S.; Karr, J.; Dwyer, G.S.; Keyworth, A.J.; Davis, C.C.; Kondash, A.J.; Vengosh, A. Occurrence and distribution of hexavalent chromium in groundwater from North Carolina, USA. Sci. Total. Environ. 2019, 711, 135135. [Google Scholar] [CrossRef]
  35. Perraki, M.; Vasileiou, E.; Bartzas, G. Tracing the origin of chromium in groundwater: Current and new perspectives. Curr. Opin. Environ. Sci. Heal. 2021, 22, 100267. [Google Scholar] [CrossRef]
  36. Vithanage, M.; Kumarathilaka, P.; Oze, C.; Karunatilake, S.; Seneviratne, M.; Hseu, Z.-Y.; Gunarathne, V.; Dassanayake, M.; Ok, Y.S.; Rinklebe, J. Occurrence and cycling of trace elements in ultramafic soils and their impacts on human health: A critical review. Environ. Int. 2019, 131, 104974. [Google Scholar] [CrossRef]
  37. Liang, J.; Huang, X.; Yan, J.; Li, Y.; Zhao, Z.; Liu, Y.; Ye, J.; Wei, Y. A review of the formation of Cr(VI) via Cr(III) oxidation in soils and groundwater. Sci. Total Environ. 2021, 774, 145762. [Google Scholar] [CrossRef]
  38. Papazotos, P.; Vasileiou, E.; Perraki, M. The synergistic role of agricultural activities in groundwater quality in ultramafic environments: The case of the Psachna basin, central Euboea, Greece. Environ. Monit. Assess. 2019, 191, 317. [Google Scholar] [CrossRef]
  39. Oze, C.; Bird, D.K.; Fendorf, S. Genesis of hexavalent chromium from natural sources in soil and groundwater. Proc. Natl. Acad. Sci. USA 2007, 104, 6544–6549. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Margiotta, S.; Mongelli, G.; Summa, V.; Paternoster, M.; Fiore, S. Trace element distribution and Cr(VI) speciation in Ca-HCO3 and Mg-HCO3 spring waters from the northern sector of the Pollino massif, Southern Italy. J. Geochem. Explor. 2012, 115, 1–12. [Google Scholar] [CrossRef]
  41. Remoundaki, E.; Vasileiou, E.; Philippou, A.; Perraki, M.; Kousi, P.; Hatzikioseyian, A.; Stamatis, G. Groundwater Deteriora-tion: The Simultaneous Effects of Intense Agricultural Activity and Heavy Metals in Soil. Procedia Eng. 2016, 162, 545–552. [Google Scholar] [CrossRef] [Green Version]
  42. Kaprara, E.; Kazakis, N.; Simeonidis, K.; Coles, S.; Zouboulis, A.; Samaras, P.; Mitrakas, M. Occurrence of Cr(VI) in drinking water of Greece and relation to the geological background. J. Hazard. Mater. 2014, 281, 2–11. [Google Scholar] [CrossRef]
  43. Elango, L.; Kannan, R. Chapter 11: Rock–Water Interaction and Its Control on Chemical Composition of Groundwater; Elsevier: Amsterdam, The Netherlands, 2007; pp. 229–243. [Google Scholar] [CrossRef]
  44. Sharif, M.; Davis, R.; Steele, K.; Kim, B.; Kresse, T.; Fazio, J. Inverse geochemical modeling of groundwater evolution with emphasis on arsenic in the Mississippi River Valley alluvial aquifer, Arkansas (USA). J. Hydrol. 2008, 350, 41–55. [Google Scholar] [CrossRef]
  45. Vasileiou, E.; Papazotos, P.; Dimitrakopoulos, D.; Perraki, M. Expounding the origin of chromium in groundwater of the Sarigkiol Basin, Western Macedonia, Greece: A cohesive statistical approach and hydrochemical study. Environ. Monit. Assess. 2019, 191, 509. [Google Scholar] [CrossRef] [PubMed]
  46. Stamos, A.; Samiotis, G.; Tsioptsias, C.; Amanatidou, E. Natural presence of hexavalent chromium in spring waters of South-West Mountain Vermion, Greece. In Proceedings of the 16th International Conference on Environmental Science and Technology (CEST 2019), Rhodes, Greece, 4–7 September 2019; p. 4. [Google Scholar]
  47. Institute of Geology and Mineral Exploration of Greece. Geological Maps of Greece, Sheet: Kozani; Scale 1:50.000, Department of Geological Maps; Institute of Geology and Mineral Exploration of Greece: Athens, Greece, 1980. [Google Scholar]
  48. Perraki, M. Mineralogical, Petrological and Geochemical Study of Heavy Minerals with Emphasis on Chromium in the Geological Formations (Ultrabasic Rocks, Lignite, Clay Formations) and the Coal-Fired Products (Fly Ash) and the Quality of Surficial and Underground Aquifers of the Sarigkiol Basin (NW Greece); Technical Report; National Technical University of Athesn: Athens, Greece, 2016. [Google Scholar]
  49. Nematollahi, M.J.; Ebrahimi, P.; Razmara, M.; Ghasemi, A. Hydrogeochemical investigations and groundwater quality as-sessment of Torbat-Zaveh plain, Khorasan Razavi, Iran. Environ. Monit. Assess. 2015, 188, 1–21. [Google Scholar] [CrossRef]
  50. Esmaeili, A.; Moore, F. Hydrogeochemical assessment of groundwater in Isfahan province, Iran. Environ. Earth Sci. 2011, 67, 107–120. [Google Scholar] [CrossRef]
  51. Spearman, C. The Proof and Measurement of Association between Two Things. Am. J. Psychol. 1904, 15, 72. [Google Scholar] [CrossRef]
  52. Wuensch, K.L.; Evans, J.D. Straightforward Statistics for the Behavioral Sciences. J. Am. Stat. Assoc. 1996, 91, 1750. [Google Scholar] [CrossRef]
  53. Gauthier, T. Detecting Trends Using Spearman’s Rank Correlation Coefficient. Environ. Forensics 2001, 2, 359–362. [Google Scholar] [CrossRef]
  54. Wilk, M.B.; Gnanadesikan, R. Probability plotting methods for the analysis for the analysis of data. Biometrika 1968, 55, 1–17. [Google Scholar] [CrossRef]
  55. Parkhurst, D.L.; Appelo, C.A.J. User’s Guide to PHREEQC (Version 2): A Computer Program for Speciation, Batch-Reaction, One-Dimensional Transport, and Inverse Geochemical Calculations; U.S. Geological Survey, Water Resources Investigations Report 99-4259; United States Geological Survey (USGS): Washington, DC, USA, 1999. [Google Scholar]
  56. Merkel, B.J.; Planer-Friedrich, B.; Nordstrom, D.K. Groundwater Geochemistry: A Practical Guide to Modeling of Natural and Contaminated Aquatic Systems; Springer: Berlin, Germany, 2005. [Google Scholar]
  57. Zhang, F.; Jin, Z.; Yu, J.; Zhou, Y.; Zhou, L. Hydrogeochemical processes between surface and groundwaters on the north-eastern Chinese Loess Plateau: Implications for water chemistry and environmental evolutions in semi-arid regions. J. Geochem. Explor. 2015, 159, 115–128. [Google Scholar] [CrossRef]
  58. Christoforidou, P.; Panagopoulos, A.; Voudouris, K. Towards A New Procedure To Set Up Groundwater Threshold Values In Accordance With The Previsions Of The Ec Directive 2006/118: A Case Study From Achaia And Corinthia (Greece). Bull. Geol. Soc. Greece 2017, 43, 1678. [Google Scholar] [CrossRef] [Green Version]
  59. Molinari, A.; Guadagnini, L.; Marcaccio, M.; Guadagnini, A. Natural background levels and threshold values of chemical species in three large-scale groundwater bodies in Northern Italy. Sci. Total Environ. 2012, 425, 9–19. [Google Scholar] [CrossRef] [PubMed]
  60. Chidichimo, F.; de Biase, M.; Straface, S. Groundwater pollution assessment in landfill areas: Is it only about the leachate? Waste Manag. 2019, 102, 655–666. [Google Scholar] [CrossRef] [PubMed]
  61. Parrone, D.; Frollini, E.; Preziosi, E.; Ghergo, S. eNaBLe, an On-Line Tool to Evaluate Natural Background Levels in Groundwater Bodies. Water 2020, 13, 74. [Google Scholar] [CrossRef]
  62. Masetti, M.; Poli, S.; Sterlacchini, S.; Beretta, G.P.; Facchi, A. Spatial and statistical assessment of factors influencing nitrate contamination in groundwater. J. Environ. Manag. 2008, 86, 272–281. [Google Scholar] [CrossRef] [PubMed]
  63. Menció, A.; Mas-Pla, J.; Otero, N.; Regàs, O.; Boy-Roura, M.; Puig, R.; Bach, J.; Domènech, C.; Zamorano, M.; Brusi, D.; et al. Nitrate pollution of groundwater; all right, but nothing else? Sci. Total Environ. 2016, 539, 241–251. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Soltan, M.E. Evaluation Of Ground Water Quality In Dakhla Oasis (Egyptian Western Desert). Environ. Monit. Assess. 1999, 57, 157–168. [Google Scholar] [CrossRef]
  65. Singh, U.V.; Abhishek, A.; Singh, K.P.; Dhakate, R.; Singh, N.P. Groundwater quality appraisal and its hydrochemical char-acterization in Ghaziabad (a region of indo-gangetic plain), Uttar Pradesh, India. Appl. Water Sci. 2013, 4, 145–157. [Google Scholar] [CrossRef] [Green Version]
  66. World Health Organization. Guidelines for Drinking-Water Quality, 4th ed.; World Health Organization: Geneva, Switzerland, 2011. [Google Scholar]
  67. Marghade, D.; Malpe, D.B.; Zade, A.B. Geochemical characterization of groundwater from northeastern part of Nagpur urban, Central India. Environ. Earth Sci. 2010, 62, 1419–1430. [Google Scholar] [CrossRef]
  68. Lelli, M.; Grassi, S.; Amadori, M.; Franceschini, F. Natural Cr(VI) contamination of groundwater in the Cecina coastal area and its inner sectors (Tuscany, Italy). Environ. Earth Sci. 2013, 71, 3907–3919. [Google Scholar] [CrossRef]
  69. Barnes, I.; O’neil, J.R. The relationship between fluids in some fresh alpine-type ultramafics and possible modern ser-pen-tinization, Western United States. Bull. Geol. Soc. Am. 1969, 80, 1947–1960. [Google Scholar] [CrossRef]
  70. Cipolli, F.; Gambardella, B.; Marini, L.; Ottonello, G.; Zuccolini, M.V. Geochemistry of high-pH waters from serpentinites of the Gruppo di Voltri (Genova, Italy) and reaction path modeling of CO2 sequestration in serpentinite aquifers. Appl. Geochem. 2004, 19, 787–802. [Google Scholar] [CrossRef]
  71. Marques, J.M.; Carreira, P.M.; Carvalho, M.D.R.; Matias, M.J.; Goff, F.E.; Basto, M.J.; Graça, R.C.; Aires-Barros, L.; Rocha, L. Origins of high pH mineral waters from ultramafic rocks, Central Portugal. Appl. Geochem. 2008, 23, 3278–3289. [Google Scholar] [CrossRef]
  72. Richard, F.C.; Bourg, A.C. Aqueous geochemistry of chromium: A review. Water Res. 1991, 25, 807–816. [Google Scholar] [CrossRef]
  73. Zhang, B.; Zhao, D.; Zhou, P.; Qu, S.; Liao, F.; Guangcai, W. Hydrochemical Characteristics of Groundwater and Dominant Water—Rock Interactions in the Delingha. Water 2020, 12, 836. [Google Scholar] [CrossRef] [Green Version]
  74. Redwan, M.; Moneim, A.A.A.; Amra, M.A. Effect of water–rock interaction processes on the hydrogeochemistry of ground-water west of Sohag area, Egypt. Arab. J. Geosci. 2016, 9, 111. [Google Scholar] [CrossRef]
  75. Jalali, M.; Khanlari, Z.V. Cadmium Availability in Calcareous Soils of Agricultural Lands in Hamadan, Western Iran. Soil Sediment. Contam. Int. J. 2008, 17, 256–268. [Google Scholar] [CrossRef]
  76. Gibbs, R.J. Mechanisms Controlling World Water Chemistry. Science 1970, 170, 1088–1090. [Google Scholar] [CrossRef] [PubMed]
  77. Marandi, A.; Shand, P. Groundwater chemistry and the Gibbs Diagram. Appl. Geochem. 2018, 97, 209–212. [Google Scholar] [CrossRef]
  78. McClain, C.; Maher, K. Chromium fluxes and speciation in ultramafic catchments and global rivers. Chem. Geol. 2016, 426, 135–157. [Google Scholar] [CrossRef] [Green Version]
  79. Giammetta, R.; Telesca, A.; Mongelli, G. Serpentinites-water interaction in the S. Severino area, Lucanian Apennines, Southern Italy. GeoActa 2004, 3, 25–33. [Google Scholar]
  80. Hanusz, Z.; Tarasińska, J. Normalization of the Kolmogorov–Smirnov and Shapiro–Wilk tests of normality. Biom. Lett. 2015, 52, 85–93. [Google Scholar] [CrossRef] [Green Version]
  81. EUROPA. European Commission Water Framework Directive 2000/60/EC. Off. J. Eur. Communities 2000, L 327, 1–73. [Google Scholar]
Figure 1. A simplified geological map of the western Vermio Mt.; the natural springs studied herein are marked.
Figure 1. A simplified geological map of the western Vermio Mt.; the natural springs studied herein are marked.
Water 13 02809 g001
Figure 2. A simplified geological section of the western Vermio Mt.
Figure 2. A simplified geological section of the western Vermio Mt.
Water 13 02809 g002
Figure 3. Simplified hydrogeological section of the natural spring “Potistis”‒W13 in western Vermio Mt.
Figure 3. Simplified hydrogeological section of the natural spring “Potistis”‒W13 in western Vermio Mt.
Water 13 02809 g003
Figure 4. Flow chart of the modified conceptual model for the assessment NBLs of the target parameter.
Figure 4. Flow chart of the modified conceptual model for the assessment NBLs of the target parameter.
Water 13 02809 g004
Figure 5. Flow chart for the assessment TVs for the target parameter.
Figure 5. Flow chart for the assessment TVs for the target parameter.
Water 13 02809 g005
Figure 6. Box plots for the physical parameters of the natural springs of western Vermio Mt.
Figure 6. Box plots for the physical parameters of the natural springs of western Vermio Mt.
Water 13 02809 g006
Figure 7. Box plots for the major ions and Si of the natural springs of western Vermio Mt.
Figure 7. Box plots for the major ions and Si of the natural springs of western Vermio Mt.
Water 13 02809 g007
Figure 8. Concentrations of Cr and Cr(VI) in the natural springs of western Vermio Mt.
Figure 8. Concentrations of Cr and Cr(VI) in the natural springs of western Vermio Mt.
Water 13 02809 g008
Figure 9. Concentrations of Cr and Cr(VI) in the natural spring “Elafakia”‒W14 of western Vermio Mt.
Figure 9. Concentrations of Cr and Cr(VI) in the natural spring “Elafakia”‒W14 of western Vermio Mt.
Water 13 02809 g009
Figure 10. Concentrations of Cr and Cr(VI) in the natural spring “Potistis”‒W13 of western Vermio Mt.
Figure 10. Concentrations of Cr and Cr(VI) in the natural spring “Potistis”‒W13 of western Vermio Mt.
Water 13 02809 g010
Figure 11. Statistical analyses of concentrations of Cr and Cr(VI) of the natural springs in western Vermio Mt.
Figure 11. Statistical analyses of concentrations of Cr and Cr(VI) of the natural springs in western Vermio Mt.
Water 13 02809 g011
Figure 12. Piper diagram of major ion chemistry for the natural spring samples.
Figure 12. Piper diagram of major ion chemistry for the natural spring samples.
Water 13 02809 g012
Figure 13. Bivariate plots of: (a) Ca vs. Mg and (b) (Ca+Mg) vs. HCO3 for the natural springs of western Vermio Mt.
Figure 13. Bivariate plots of: (a) Ca vs. Mg and (b) (Ca+Mg) vs. HCO3 for the natural springs of western Vermio Mt.
Water 13 02809 g013
Figure 14. Gibbs diagrams of the natural springs in western Vermio Mt.
Figure 14. Gibbs diagrams of the natural springs in western Vermio Mt.
Water 13 02809 g014
Figure 15. Average concentrations of Cr vs. water types of the natural springs in western Vermio Μt.
Figure 15. Average concentrations of Cr vs. water types of the natural springs in western Vermio Μt.
Water 13 02809 g015
Figure 16. Fluctuation diagrams of rainfall and Cr concentrations of the natural spring “Potistis”‒W13 in western Vermio Mt.
Figure 16. Fluctuation diagrams of rainfall and Cr concentrations of the natural spring “Potistis”‒W13 in western Vermio Mt.
Water 13 02809 g016
Figure 17. Plot of Cr (μg/L) vs. discharge (L/h) in the natural spring “Potistis”‒W13 in western Vermio Mt.
Figure 17. Plot of Cr (μg/L) vs. discharge (L/h) in the natural spring “Potistis”‒W13 in western Vermio Mt.
Water 13 02809 g017
Figure 18. Plots of: (a) Mg/Si vs. Cr and (b) Cr(VI) vs. Cr (regression model) of the natural springs in western Vermio Mt.
Figure 18. Plots of: (a) Mg/Si vs. Cr and (b) Cr(VI) vs. Cr (regression model) of the natural springs in western Vermio Mt.
Water 13 02809 g018
Figure 19. Plots of saturation indices vs. Cr of the natural springs in western Vermio Mt.
Figure 19. Plots of saturation indices vs. Cr of the natural springs in western Vermio Mt.
Water 13 02809 g019
Figure 20. Box plots and normality tests for the natural spring “Elafakia”‒W14 in western Vermio Mt.
Figure 20. Box plots and normality tests for the natural spring “Elafakia”‒W14 in western Vermio Mt.
Water 13 02809 g020
Figure 21. Box plots and normality tests for the natural spring “Potistis”‒W13 in western Vermio Mt.
Figure 21. Box plots and normality tests for the natural spring “Potistis”‒W13 in western Vermio Mt.
Water 13 02809 g021
Figure 22. Flow chart of the modified conceptual model for assessing NBLs and TVs of Cr and Cr(VI) of the natural spring “Potistis”‒W13, in western Vermio Mt.
Figure 22. Flow chart of the modified conceptual model for assessing NBLs and TVs of Cr and Cr(VI) of the natural spring “Potistis”‒W13, in western Vermio Mt.
Water 13 02809 g022
Figure 23. The flow paths and the natural recharge from western Vermio Mt. towards the eastern part of the Sarigkiol Basin (Google Earth image, 2021).
Figure 23. The flow paths and the natural recharge from western Vermio Mt. towards the eastern part of the Sarigkiol Basin (Google Earth image, 2021).
Water 13 02809 g023
Table 1. Maximum, minimum, and median values of physical and chemical parameters of natural springs in western Vermio Mt.
Table 1. Maximum, minimum, and median values of physical and chemical parameters of natural springs in western Vermio Mt.
ParameterUnitQLDLThe Agio Pnevma AreaThe Agios Panteleimonas Area“Mouratidis”The Agios Dimitrios AreaThe Vazelona Area
MaxMinMedianMaxMinMedianMaxMinMedianMaxMinMedianMaxMinMedian
pH---8.47.67.97.947.377.388.47.88.18.58.18.37.77.77.7
DOmg/L--9.68.48.99.37.748.539.27.88.88.98.18.59.58.38.9
T°C--16.010.814.318.013.113.625.69.813.824.1816.115.815.015.4
TDSmg/L--561.3150.5226.0397.69361.25377.36398294.4369.8377319.3348.2462.2383.5422.9
ECμS/cm10-593.0293.0460.0505.0448.0494.0520405426.8456446451484.0389.0436.5
EhmV--320.0110.0160.0389.7303.0387.0409303.42346.2377340.7358.9301.0297.0299.0
Ca2+mg/L0.20.05119.024.043.894.941.693.260.454.255.249.546.648.1104.098.2101.1
Mg2+mg/L1.00.324.93.78.438.83.1311.434.921.321.431.830.931.413.53.18.3
Na+mg/L5.00.51.41.41.42.5BDL1.251.7BDL1.41.2BDL0.62.11.01.6
K+mg/L0.20.0533.00.40.61.630.591.2210.30.31.81.51.710.81.46.1
NO3mg/L5.01BDLBDLBDLBDLBDLBDL9.18.38.7BDLBDLBDL1.0BDLBDL
Clmg/L5.0131.02.07.5BDLBDLBDL5BDLBDLBDLBDLBDL12.01.06.5
SO42−mg/L10.0231.010.020.519.012.019.020161622BDL1121.013.017.0
HCO3mg/L10.02387.0119.0167.0276.0250.0271.0277192258271240255.5304.0259.0281.5
Alμg/LDL18.02.03.01.01.01.0154624313173.02.02.5
Asμg/LDL0.50.70.60.75.40.61.56.11.41.749.128.8391.80.51.2
Bμg/LDL59.05.06.012.012.012.0117910101019.015.017.0
Baμg/LDL0.0511.61.92.68.695.725.8113.866.916.114.415.27.06.46.7
Brμg/LDL525.08.013.014.011.012.01613141591215.013.014.0
Crμg/LDL0.147.80.53.818.01.51.938.31020.416.610.913.80.80.50.7
Cr(VI)μg/LDL0.136.70.51.87.21.01.033.971616.58.712.60.50.10.3
Cuμg/LDL0.12.80.91.21.20.60.720.71.41.31.11.22.81.12.0
Liμg/LDL0.10.70.30.64.30.10.23.40.50.64.73.64.20.20.10.2
Mnμg/LDL0.055.80.20.80.430.280.291.30.60.81.40.40.90.40.30.4
Niμg/LDL0.21.71.01.41.30.81.057.70.54.15.83.64.71.00.70.9
Pμg/LDL1011712.057.031.015.023.039373945132928.018.023.0
Siμg/LDL4014,3272231429624,8753250344121,350878610,69719,30717,24518,276349534673481
Srμg/LDL0.0184.044.559.879.0341.7670.1562.957.558.465.260.762.983.269.776.5
Uμg/LDL0.020.2BDL0.17.20.40.70.30.30.30.40.30.30.70.20.4
Vμg/LDL0.20.70.20.30.080.080.084.61.31.54.74.34.50.80.50.7
Znμg/LDL0.512.77.510.44.33.43.914.891311.97.39.633.113.623.4
BDL: Below the detection limit. DL: Detection limit. QL: Quantification limit.
Table 2. Maximum, minimum, and median values of physical and chemical parameters of the natural springs in western Vermio Mt.
Table 2. Maximum, minimum, and median values of physical and chemical parameters of the natural springs in western Vermio Mt.
ParameterUnitQLDLPotistisElafakia
MaxMinMedianMaxMinMedian
pH---8.37.37.98.37.37.7
DOmg/L--9.68.59.011.68.69.2
T°C--15.26.212.320.55.812.5
TDSmg/L--528.9386.3481.8522.5366.4458.0
ECμS/cm10-620.0374.8574.5718.0357.7546.0
EhmV--359.890.0325.5377.6194.9309.8
Ca2+mg/L0.20.0556.128.635.691.451.776.9
Mg2+mg/L1.00.373.334.261.736.624.830.3
Na+mg/L5.00.51.21.21.22.5BDL2.0
K+mg/L0.20.050.50.10.33.10.60.7
NO3mg/L5.01BDLBDLBDL1.0BDLBDL
Clmg/L5.011.0DL1.08.0BDL2.0
SO42−mg/L10.02BDLBDLBDL128.016.023.0
HCO3mg/L10.02409.0298.0382.0369.0250.0318.0
Alμg/LDL1111.01.03.013.01.02.0
Asμg/LDL0.5BDLBDLBDL1.70.91.5
Bμg/LDL513.08.010.527.06.014.0
Baμg/LDL0.0510.44.24.816.210.415.0
Brμg/LDL518.010.014.028.018.021.0
Crμg/LDL0.1131.539.0103.957.426.047.5
Cr(VI)μg/LDL0.1100.039.090.051.218.041.0
Cuμg/LDL0.14.00.30.56.40.61.3
Liμg/LDL0.11.00.70.81.20.91.1
Mnμg/LDL0.050.80.10.40.80.20.3
Niμg/LDL0.238.22.75.5314.05.17.0
Pμg/LDL1078.013.038.0147.018.032.0
Siμg/LDL4022,39415,12020,66317,71713,24814,685
Srμg/LDL0.0140.332.035.290.458.374.0
Uμg/LDL0.020.1BDLBDL1.10.60.8
Vμg/LDL0.21.60.20.41.20.71.0
Znμg/LDL0.537.34.15.174.74.68.7
BDL: Below the detection limit. DL: Detection limit. QL: Quantification limit.
Table 3. Geographical coordinates of the water sampling sites and the results of Cr-Cr(VI) for the natural springs in western Vermio Mt.
Table 3. Geographical coordinates of the water sampling sites and the results of Cr-Cr(VI) for the natural springs in western Vermio Mt.
Sample IDLatitudeLongitudeSampling PointCr (μg/L)Cr(VI) (μg/L)Sample IDLatitudeLongitudeSampling PointCr (μg/L)Cr(VI) (μg/L)
W13_06_201840°27′103″21°57′639″Potistis41.639.2W14_9a_201840°25′854″21°56′878″Elafakia46.043.0
W13_07_201840°27′103″21°57′639″Potistis39.039.0W14_9b_201840°25′854″21°56′878″Elafakia56.140.0
W13_09_2018a40°27′103″21°57′639″Potistis111.590.0W14_9c_201840°25′854″21°56′878″Elafakia56.340.0
W13_09_2018b40°27′103″21°57′639″Potistis109.590.0W14_9d_201840°25′854″21°56′878″Elafakia56.040.0
W13_09_2018c40°27′103″21°57′639″Potistis108.790.0W14_9e_201840°25′854″21°56′878″Elafakia54.541.0
W13_09_2018d40°27′103″21°57′639″Potistis112.690.0W14_05_201940°25′854″21°56′878″Elafakia42.540.0
W13_10_2018a40°27′103″21°57′639″Potistis131.5100.0W14_08_201940°25′854″21°56′878″Elafakia44.833.0
W13_10_2018b40°27′103″21°57′639″Potistis130.2100.0W14_11_201940°25′854″21°56′878″Elafakia46.341.0
W13_10_2018c40°27′103″21°57′639″Potistis111.990.0W14_02_202040°25′854″21°56′878″Elafakia47.532.0
W13_10_2018d40°27′103″21°57′639″Potistis127.8100.0W14_07_202040°25′854″21°56′878″Elafakia46.618.0
W13_10_2018e40°27′103″21°57′639″Potistis110.790.0W14_09_202040°25′854″21°56′878″Elafakia48.033.0
W13_11_201840°27′103″21°57′639″Potistis127.5100.0S10_11_201440°26′854″21°58′711″Agio Pnevma47.836.7
W13_04_201940°27′103″21°57′639″Potistis89.289.0S15_06_201840°26′689″21°58′801″Agio Pnevma15.915.2
W13_05_201940°27′103″21°57′639″Potistis92.890.0S10_06_201840°26′854″21°58′711″Agio Pnevma18.515.6
W13_06_201940°27′103″21°57′639″Potistis98.189.0S16_07_201840°27′338″21°58′750″Agio Pnevma2.41.8
W13_08_201940°27′103″21°57′639″Potistis103.388.0S17_07_201840°26′856″21°58′757″Agio Pnevma3.81.5
W13_10_201940°27′103″21°57′639″Potistis103.592.0S18_07_201840°27′871″21°58′475″Agio Pnevma2.21.6
W13_11_201940°27′103″21°57′639″Potistis105.899.0S19_07_201840°27′095″21°58′711″Agio Pnevma0.50.5
W13_02_202040°27′103″21°57′639″Potistis95.787.0S2_03_201440°25′789″21°56′216″Mouratidis38.333.9
W13_06_202040°27′103″21°57′639″Potistis99.097.0S2_09_201640°25′789″21°56′216″Mouratidis17.013.0
W13_07_202040°27′103″21°57′639″Potistis99.682.0S2_02_201740°25′789″21°56′216″Mouratidis10.07.00
W13_09_202040°27′103″21°57′639″Potistis103.988.0S2_04_201740°25′789″21°56′216″Mouratidis15.012.0
W13_10_202040°27′103″21°57′639″Potistis103.592.0S2_05_201740°25′789″21°56′216″Mouratidis28.023.0
W14_11_201440°25′854″21°56′878″Elafakia53.551.1S2_06_201740°25′789″21°56′216″Mouratidis25.022.0
W14_07_201440°25′854″21°56′878″Elafakia57.451.2S2_07_201740°25′789″21°56′216″Mouratidis26.019.0
W14_12_201540°25′854″21°56′878″Elafakia52.549.1S2_08_2017a40°25′789″21°56′216″Mouratidis12.08.00
W14_09_201640°25′854″21°56′878″Elafakia26.023.0S2_08_2017b40°25′789″21°56′216″Mouratidis23.721.4
W14_04_201740°25′854″21°56′878″Elafakia42.041.0S2_09_201740°25′789″21°56′216″Mouratidis16.010.0
W14_05_201740°25′854″21°56′878″Elafakia47.347.0S1_03_201440°25′224″21°55′889″Agios Dimitrios16.616.5
W14_06_201740°25′854″21°56′878″Elafakia49.047.0S1_08_201740°25′224″21°55′889″Agios Dimitrios10.98.70
W14_08_201740°25′854″21°56′878″Elafakia49.246.8S13_06_201740°25′810″21°56′900″Agios Panteleimonas1.501.00
W14_10_201740°25′854″21°56′878″Elafakia45.435.2S14_06_201740°25′801″21°57′099″Agios Panteleimonas1.901.00
W14_07_201840°25′854″21°56′878″Elafakia45.040.0W21_08_201940°25′842″21°56′856″Agios Panteleimonas18.07.2
W14_08_201840°25′854″21°56′878″Elafakia49.242.0S5_07_201440°25′713″21°56′725″Vazelonas0.500.10
W14_10_201840°25′854″21°56′878″Elafakia45.441.0S6_07_201440°25′840″21°56′882″Vazelonas0.800.50
Table 4. The abundance of major ions among the natural springs in western Vermio Mt.
Table 4. The abundance of major ions among the natural springs in western Vermio Mt.
Area/Sampling SiteSample IDCations OrderAnions Order
1Agios Dimitrios areaS1Ca2+ > Mg2+ > K+ > Na+HCO3 > SO42− > Cl > NO3
2ElafakiaW14Ca2+ > Mg2+ > Na+ > K+HCO3 > SO42− > Cl > NO3
3Agios Panteleimonas areaS2, S13, S14Ca2+ > Mg2+ > Na+ > K+HCO3 > SO42−
4PotistisW13Mg2+ > Ca2+ > K+HCO3 > SO42− > Cl > NO3
5Agio Pnevma areaS18, S16, S19, S10, S17, S15Ca2+ > Mg2+ > K+ > Na+HCO3 > SO42− > Cl
6MouratidisS2Ca2+ > Mg2+ > Na+ > K+HCO3 > SO42− > NO3 > Cl
7Vazelona areaS5, S6Ca2+ > Mg2+ > K+ > Na+HCO3 > SO42− > Cl
Table 5. Spearman’s rank correlation matrix of selected major and trace elements of the natural springs of western Vermio Mt.
Table 5. Spearman’s rank correlation matrix of selected major and trace elements of the natural springs of western Vermio Mt.
ParameterpHDOECEhCa2+Mg2+Na+K+NO3ClSO42−HCO3AlAsBBaBrCrCr(VI)LiMnNiPSiSrUVZnAlkalinity
pH1
DO−0.111
EC−0.150.211
Eh−0.110.060.191
Ca2+−0.180.10−0.200.151
Mg2+0.060.050.68 **0.32−0.42 *1
Na+−0.17−0.06−0.16−0.150.49 **−0.181
K+0.12−0.27−0.28−0.060.51**−0.511 **0.131
NO30.18−0.04−0.41 *0.130.31−0.180.47 **0.031
Cl−0.04−0.15−0.30−0.430.29−0.461 **0.37 *0.39 *0.191
SO42−−0.06−0.070.09−0.080.67 **−0.210.52 **0.49 **0.200.251
HCO3−0.110.040.62 **0.140.000.68 **−0.15−0.31−0.22−0.140.011
Al0.170.10−0.40 *−0.24−0.07−0.080.09−0.030.180.28−0.14−0.191
As0.17−0.20−0.100.190.42 *0.000.51**0.49 **0.33−0.040.60 **−0.250.091
B−0.110.310.39 *0.250.190.37 *0.040.000.20−0.190.250.41 *−0.140.181
Ba−0.010.040.180.040.47 **0.140.260.41 *0.12−0.100.63 **0.180.060.67 **0.45 **1
Br−0.170.100.28−0.220.33−0.020.320.210.040.310.70 **0.21−0.050.200.260.45 **1
Cr−0.110.320.71 **0.19−0.300.76 **−0.09−0.55 **−0.19−0.33−0.130.55 **−0.12−0.160.38 *0.070.171
Cr(VI)−0.140.320.73 **0.18−0.340.80 **−0.10−0.6 **−0.19−0.35 *−0.170.59 **−0.09−0.190.36 *0.060.170.98 **1
Li0.170.040.42 *0.13−0.170.611 **−0.080.05−0.15−0.320.230.220.110.48 **0.39 *0.63 **0.290.47 **0.48 **1
Mn0.17−0.26−0.30−0.44 *−0.15−0.19−0.130.040.050.06−0.07−0.110.30−0.04−0.43−0.05−0.09−0.41 *−0.32−0.141
Ni−0.270.300.60 **0.07−0.130.60**−0.10−0.25−0.25−0.300.010.51 **0.100.040.42 *0.45 **0.37*0.61 **0.67 **0.64 **−0.101
P−0.09−0.07−0.27−0.41 *−0.08−0.28−0.15−0.02−0.030.10−0.070.000.25−0.23−0.30−0.13−0.04−0.36 *−0.26−0.240.75 **−0.041
Si0.10−0.070.58 **0.34−0.43 *0.91 **−0.20−0.34−0.11−0.43 *−0.190.48 **−0.070.090.340.17−0.090.75 **0.76 **0.71 **−0.220.53 **−0.331
Sr−0.130.01−0.26−0.180.827 **−0.550.61 **0.61 **0.210.38 *0.74 **−0.230.060.53 **0.120.56 **0.47 **−0.37 *−0.43−0.04−0.09−0.13−0.07−0.531
U−0.140.12−0.03−0.090.676 **−0.270.59 **0.462 **0.280.150.82 **−0.150.050.72 **0.330.78 **0.58 **−0.16−0.180.29−0.050.18−0.03−0.280.85 **1
V0.100.02−0.080.080.150.160.300.160.33−0.110.41 *−0.130.250.70 **0.210.69 **0.19−0.10−0.090.59 **0.070.26−0.070.220.240.54 **1
Zn−0.01−0.10−0.50 **−0.230.22−0.330.240.180.52 **0.38 *0.04−0.220.320.09−0.050.09−0.04−0.23−0.22−0.180.34−0.160.41 *−0.240.230.210.091
Alkalinity−0.110.040.62 **0.140.000.68 **−0.15−0.31−0.22−0.140.021 **−0.19−0.240.41 *0.190.220.55 **0.59 **0.23−0.110.50 **0.000.48 **−0.22−0.14−0.13−0.221
* Correlation is significant at the 0.05 level (2-tailed). ** Correlation is significant at the 0.01 level (2-tailed).
Table 6. The percentage distribution of the saturation state with respect to the selected mineral phases in the collected water samples.
Table 6. The percentage distribution of the saturation state with respect to the selected mineral phases in the collected water samples.
Mineral PhaseOversaturatedUndersaturated
SI > 0 (%)SI < 0 (%)
Calcite63.836.2
Dolomite58.042.0
Magnesite29.071.0
Talc49.350.7
Chlorite30.469.6
Tremolite14.585.5
Enstatite0.00100
Diopside0.00100
Pyrolusite0.00100
Chromite1000.00
Magnetite1000.00
Chrysotile14.585.5
Lizardite14.585.5
Olivine (Forsterite)0.00100
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vasileiou, E.; Papazotos, P.; Dimitrakopoulos, D.; Perraki, M. Hydrogeochemical Processes and Natural Background Levels of Chromium in an Ultramafic Environment. The Case Study of Vermio Mountain, Western Macedonia, Greece. Water 2021, 13, 2809. https://doi.org/10.3390/w13202809

AMA Style

Vasileiou E, Papazotos P, Dimitrakopoulos D, Perraki M. Hydrogeochemical Processes and Natural Background Levels of Chromium in an Ultramafic Environment. The Case Study of Vermio Mountain, Western Macedonia, Greece. Water. 2021; 13(20):2809. https://doi.org/10.3390/w13202809

Chicago/Turabian Style

Vasileiou, Eleni, Panagiotis Papazotos, Dimitrios Dimitrakopoulos, and Maria Perraki. 2021. "Hydrogeochemical Processes and Natural Background Levels of Chromium in an Ultramafic Environment. The Case Study of Vermio Mountain, Western Macedonia, Greece" Water 13, no. 20: 2809. https://doi.org/10.3390/w13202809

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop