Next Article in Journal
The Impact of Decreased Atmospheric Pressure on Forced Aeration of Discharged Flow
Previous Article in Journal
Appraisal of Daily Temperature and Rainfall Events in the Context of Global Warming in South Australia
Previous Article in Special Issue
Multi-Tier Validation of a Macroscale Nitrogen Model for Groundwater Management in Watersheds Using Data from Different Monitoring Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Index-Based Groundwater Quality Assessment of Nestos River Deltaic Aquifer System, Northeastern Greece

by
George Kampas
1,
Andreas Panagopoulos
2,
Ioannis Gkiougkis
1,*,
Christos Pouliaris
3,
Fotios-Konstantinos Pliakas
1,
Vasiliki Kinigopoulou
2 and
Ioannis Diamantis
1
1
Laboratory of Engineering Geology and Groundwater Research, Department of Civil Engineering, Democritus University of Thrace, 67100 Xanthi, Greece
2
Soil and Water Resources Institute, Hellenic Agricultural Organization, 57400 Sindos, Greece
3
School of Mining and Metallurgical Engineering, National Technical University of Athens, Iroon Polytechniou 9 Str., Zografou Campus, 15773 Athens, Greece
*
Author to whom correspondence should be addressed.
Water 2024, 16(2), 352; https://doi.org/10.3390/w16020352
Submission received: 8 December 2023 / Revised: 12 January 2024 / Accepted: 16 January 2024 / Published: 21 January 2024

Abstract

:
The Nestos River delta is one of the most important and sensitive basins in Greece and Europe due to its ecosystem functions, combining intensive agricultural production with low-enthalpy geothermal energy and important ecotopes. High water quality is of paramount importance to the sustainability of the system. Systematic and continuous assessment of water quality needs to be carried out in a way that is easy and quick for decision makers and non-expert societal partners to comprehend. In this way, decisions may be made more rapidly, and involved water users may be sensitized to rational water use. To this end, this paper presents the assessment of groundwater quality in the Nestos River’s western delta with the use of Poseidon (PoS), a versatile, index-based method. Groundwater samples collected from 24 and 22 wells tapping the unconfined and the confined aquifers, respectively, in four time periods (May and October 2019 and 2020) were analyzed. Using the PoS index, groundwater samples were classified according to their quality status, highlighting the parameters driving quality degradation issues, thus assisting water managers in obtaining an overview of quality status and evolution through datasets that were often large. PoS index is applied in the study area for the first time and provides a groundwater quality assessment through a unique score representative of the overall water quality status regardless of processes (anthropogenic or natural) or any kind of pressures.

1. Introduction

Groundwater quality assessment is often founded on traditional methods based solely on the comparison of analytical parametric values or calculation of molar ratios [1,2,3]. These methods can be very helpful, but, in most cases, they do not provide a convenient supervisory correlation between the examined samples, whilst identification of trends often proves to be unclear.
A variety of graphical techniques (those of Piper, Durov, Wilcox, etc., and diagrams) that have been proposed in the literature are also widely used [4,5,6] along with numerous synthetic ratios such as the sodium adsorption ratio and the Revelle factor [7]. Multivariate statistical techniques, such as correlation analysis, clustering analysis (CA), factor analysis/principle component analysis (FA/PCA), etc., are widely applied to understand processes, correlations, and hydrogeochemical evolution trends in groundwater systems. Correlation analysis is carried out as a bivariate statistic to determine the mutual relationships and strength of association between pairs of variables through calculation of the linear Pearson’s correlation coefficient [8]. Cluster analysis is a multivariate method based on the categories of data according to predetermined selection criteria. In this method, data are placed in clusters that indicate high internal (within-cluster) homogeneity and high external (between-cluster) heterogeneity [9]. Statistical methods such as those mentioned resolve issues relating to the identification of correlation and clustering of parameters that may be linked to the hydrogeochemical identity of the studied water samples, whilst they may also prove invaluable in deciphering or, at least, proving the hydrochemical types identified with the use of parametric ratios and graphical representation in one of the numerous available diagrams. Isotopic analyses are a valuable tool that often provide a higher resolution with regards to water origin where other conventional methods, such as those discussed above, fail.
Even though, in most of cases, the appropriate combination of some or all of the above discussed methods may provide definitive results with regards to groundwater origin, dominant evolution mechanisms, and spatiotemporal trends [10,11,12], they are not easy to be comprehend for non-experts, and they fail to provide an on-the-fly overview of the quality status to decision makers. Undoubtedly, reliable and objective assessments need to be based on systematic monitoring and consideration of a broad spectrum of parameters, which in several cases can be puzzling to decision makers. On the other hand, methods such as those mentioned above do not often allow for direct comparisons between different systems exposed to variable pressures and characterized by diverse hydrogeochemical settings and hydrodynamic evolution mechanisms.
Index-based methods may provide an alternative to overcome this shortfall, as they offer a quick assessment of groundwater quality in the form of a simple arithmetic value. They were introduced in the 1960s and were already adopted as official tools by water resource monitoring agencies in the 1970s [13]. Depending on the method, a different number of physicochemical and, in some cases, microbiological parameters are assessed, and quality evaluation is usually based on the criteria set for a particular use. Their usability, amongst others, relies on the fact that they reduce the inherent complexity of water sample analysis to a simple arithmetic value which is easy to comprehend. Over 30 such methods have been reported by other researchers [14]; these methods use a widely varying number of parameters ranging from 3 to 72. Often, index-based approaches adopt specific water-use quality standards, thereby potentially masking the overall quality attributes of a sample. Whilst a review of the developed index-based methods is not within the scope of this paper, it is worth noting the USA and the Canadian Water Quality Index (CWQI) methods as milestones and amongst the most well-known and widely used ones. A thorough benchmark study amongst these methods concluded that the CWQI yields more stringent classifications, i.e., provides a more end use-sensitive result, compared to the USA methods [15]. As is also the case with the USA methods, the CWQI considers as baseline for its assessments the potability criteria and classifies assessed water samples in five grades on a scale of 0–100, with 0 representing the worst and 100 the best water quality [16]. The number of parameters that are suggested for use in the assessment varies from 8 to 20; however, less than 8 or more than 20 may be used [17]. Calculation of the CWQI is based on three factors, namely, the scope, which represents the percentage of parameters that do not meet the quality standards at least once over the considered period of evaluation; the frequency, which refers to the percentage of individual tests failing the standards; and the amplitude, which represents the amount by which the failed test values do not meet the standards [16].
A relatively new method, the Poseidon index (PoS), is an alternative, versatile, index-based method for the assessment and classification of water samples [10]. The concept of the method relies on the fact that each water sample has an environmental quality status that, however, is shaped as a result of geogenic and anthropogenic factors and, thus, may not be uniquely associated to specific maximum admissible concentrations assumed for a particular use. Therefore, the method seeks to examine and assess the environmental quality status of a water sample, accounting for the origin of each of the compounds that synthesizes it. Any number of parameters may be considered in the assessment, and each of them is weighted against its toxicity classification as this is referenced in the key document of ATSDR [18]. Specific use standards are being considered, and, normally, these are the most stringent potability standards and, therefore, are environmentally sensitive. However, these standards only serve as a basis for initiating the assessment, and by no means are they related to the actual classification of the water, as explained in the following paragraphs. Classification is based on a scale of 0 to over 2400, where 0 corresponds to no/minimum degradation, and over 2400 corresponds to a severe degradation of the water’s quality. This scale is organized in six color-coded classes and corresponding subclasses (a, b, c) for classes two to five, thus enabling a more detailed categorization of assessed qualities. Moreover, the method allows for the identification of the parameters that control the quality assessment, helping to trace and correlate quality issues to their origin.

2. Poseidon (PoS) Index

PoS is a relatively new index-based operational method developed to address the aforementioned issues in assessing groundwater quality. PoS differs from other methods in that it provides a general operational comparison between different sampling points, irrespective of the measured parameters, sampling times, and time deviations (comparison between different points in different catchment areas and sampling periods) [3,10]. It allows temporal comparisons between different sampling periods in the same basin and assessments between different groundwater systems in the same or different basins, under different pressures, and controlled by different mechanisms [11].
Compared to other quality assessment methods that are proposed in the literature, PoS presents some key differences [10]. The main differences are [19]: (a) the calculated PoS index does not change with the parameters used (number and nature) and (b) other methods focus either on (i) specific objectives (e.g., sea water intrusion, chemical status), (ii) specific processes (e.g., origin of water), or, finally, (iii) end use-related characterization (e.g., quality of drinking or irrigation water).
PoS provides quality assessment through a unique score representative of the overall water quality status regardless of processes (anthropogenic or natural), objectives, and characterization, as described in the foundation paper of the index [10], and has been successfully applied in a series of cases of complex systems’ quality assessments [12,19,20], in conjunction with statistical or more conventional hydrochemical analysis methods, or in comparatively assessing water quality characteristics amongst vastly different hydrogeological basins [11].
It should be noted that PoS index not only provides the environmental quality status of a groundwater sample but also identifies the chemical parameter(s) which causes the deviation from a “good” environmental status, thus narrowing down the likely sources of the problem. Thus, the proposed approach directly relates the outcome of the assessment to the concentrations of selected chemical parameters which may have adverse effects on health and the environment. PoS serves, therefore, as a robust and fully comprehensive means of reducing the inherent complexity in hydrochemical data sets to an easy-to-understand quality index, and, as such, it does not reflect on end user or different types of water regulatory frameworks. The following table, Table 1, provides the parameters used in this paper for calculating the PoS index and classifies these parameters based on their toxicity, individual score, and weighting factor.
The calculation of the PoS index is based on the use of the most usual inorganic substances commonly found in groundwater. These substances may originate from a variety of sources such as agriculture, sea intrusion, rock weathering, industry, and household waste. The method provides a wide list of parameters that, ideally, should all be considered. However, in several cases, not all listed parameters are available, whilst, in other cases, quality evaluation focuses on a particular objective to which not all of the listed parameters are relevant. Each parameter selected is classified according to its overall (synergistic/additive) health and environmental impact. The toxicity of each agent to humans plays a key role in this; their classification is based on the “Priority List of Hazardous Substances” [18] but was further elaborated by toxicological data obtained from various published sources [21,22] and well-established international organizations, such as the WHO [23], in order to ensure a comprehensive examination of critical aspects of water quality.
The classification of each parameter is based on its toxicity; this is also a key difference from most index-based methods, and the main objective is the relative comparison between concentrations of the selected parameters for the purposes of the PoS index. Based on the above, the selected parameters are classified, accordingly, into six categories (P-class), as shown in Table 1.
Depending on the overall incidence of toxicity in humans, a partial score (PS) is assigned to each P-class which reflects the degree of overall toxicity for humans. The reference values used were based on those applied in references [10,11,12]. The individual PS score follows an exponential approach, but the scale between them is carefully chosen to differentiate between the recognized P-class classes without overlapping the lower ones.
The weighting factor/Wf for each parameter is subsequently calculated, based on the partial score, according to Equation (1):
Wfi = P S i Σ P S i
The partial scores were selected through a procedure using a trial-and-error approach. The main purpose of this procedure was to obtain the best possible results so that the weighting factor (Wf) of a parameter belonging to the P-class (V) will be a rank higher than that of a parameter belonging to the IV class. Thus, for Cd or Pb, the Wf was calculated at 0.15, while for Cu or Mn, the Wf was calculated at 0.03. This is because, while the differences between categories I and II (non-toxic) and V and VI (highly and severely toxic), in terms of overall environmental impact, are insignificant, the effects between non-toxic and toxic categories must be clearly distinguished.
The reliability of the final weights was confirmed by multivariate techniques (object analysis) by calculating Cronbach’s alpha (a) coefficient, which represents the internal consistency of the values considered [24]. The coefficient (a) ranging between 0 and 1 and the highest values indicate greater consistency between the examined values, with a historical reference point of acceptable values greater than 0.7. The lower the value in column points, the higher the P-class. The points in the following table, Table 1, are based on the value of toxicity [18].
The impact of each parameter on the overall quality assessment of each sample is attributed through each quality contribution factors (Qf). The quality contribution factors are derived from the following Equation (2):
Qfi = [ ( C i W i ) M A C i ] × 10 3
where Qfi is the quality contribution factor of the parameter i; Ci is the concentration of the parameter i (same units as the parameter); Wi is the weighting factor of the ith parameter; and MACi is the maximum parametric value of the ith parameter (the maximum acceptable concentration specified in Directive 2020/2184 EC [25] for water intended for human consumption; this serves only as a reference basis and does not actually interfere in the final quality characterization).
The final PoS index is obtained by summing all the individual Qfis according to the following equation, Equation (3). The PoS index is a dimensionless number that can be used for the qualitative evaluation of a sample.
PoS = Σ Qfi
The contribution of each parameter to the overall qualitative status of the sample, as a percentage, is derived from the Qfi calculation ratio. The prevailing control parameters are those with a contribution rate greater than 20% [10].
The PoS index score calculated for each sample classifies it into a category according to a classification system, based on a six-base color palette, according to the following table, Table 2. This is carried out so that the calculated PoS index score, which is simply a number, is represented by a color to understand and interpret each sample’s score. In the PoS index, every one of the six classes ranges from 0 to t, from t to 2t, from 2t to 3t, etc. Classes 2-5 are divided into three subclasses (a, b and c), offering a better resolution between the tested samples. Each of the main categories (except the first and last) is divided into three equally wide subcategories, on a progressively increasing scale, in order to ensure better representativeness and differentiation of the results. Categories 1 and 6 need not be further differentiated as they represent boundary quality conditions. As a result, these classes can be considered as critical limits of higher and lower groundwater quality, respectively. The upper and lower limits of each one of the three subclasses are tabulated in Table 2.
The procedure for calculating a value corresponding to the threshold of the PoS index was carried out based on a reference sample and with concentrations reflecting the average typical concentrations expected in natural groundwater for the selected parameters. The value used to define the classes and subcategories of the PoS index is 150 [1] and can be used as the fundamental constant step for the classification of the PoS index, regardless of the parameters involved at one time. Its calculation was based on several trial-and-error tests which were performed with various combinations of parameters and which produced similar values. The higher the class of the PoS index, the lower the overall environmental quality of the water, expressed by the color hue of the PoS index.
The final stage of the process involves the evaluation through a comprehensive and simplified presentation which provides the key elements of the indicator. This process includes:
(i)
The main category of the PoS index (1 to 6) and its subcategory, if any;
(ii)
The dominant factor (or factors with a Qfi > 20%) which controls the quality status of the sample;
(iii)
A final unique color characterization, which makes the data visually clearer, followed by accurate information on the level of quality degradation.
It is noted that the term “quality degradation” refers to a comparative characterization among the samples examined under the PoS index and, in no case, does it imply any specific characterization under the relevant legislation and/or international standards.

3. Case Study

The study area is located at the western part of the Nestos River delta (Figure 1) at the eastern coastal region of the Kavala Prefecture, Northeastern Greece. It is delineated by the Thracian Sea coastline to the south and west, the riverbed of the Nestos River to the east, and the Rhodope geological massif to the north, thus including the full length of the Nestos River along its deltaic plain. The river delta belongs to the National Park of Eastern Macedonia and Thrace, which was founded in the year 2008. In this section of the Nestos River delta, a great number of settlements are present. The river visibly splits the delta in almost two equal parts, expanding from the settlements of Nea Karvali to the west, up to Toxotes to the north, and Avdira to the east.
The deltaic plain slopes gently to the south with elevations ranging from +40 m at the northern margin to 0 ± 1 m at the southern end near the coastline. Limited-size morphological depressions are identified from north to south; these give rise to secluded seasonally filled stagnant water bodies (lakes or swamps). These features, and especially the marshes, were most distinguished in the early development phases of the delta and primarily owe their existence to residues of older secondary river courses or to mountain-zone currents.
In the whole study area (210.73 km2), agriculture is the dominant land use (over 77% landcover). A factor that contributes to deep percolation of groundwater recharge is flood irrigation, still practiced in many parts of the plain. An aftermath of this practice appears to be the prevention of seawater intrusion advancing into the coastal aquifer system; however, at the same time, this practice potentially affects groundwater quality by leaching pollutants originating from anthropogenic activities into the saturated zone. Finally, in recent years, shallow wells are mainly used for irrigation activities, while the use of deep wells is very limited and constrained to a short period of time [26].

4. Geo-Hydrological Setting

The deltaic region consists of Holocene sediments with an overall thickness of some tens of meters deposited mainly by the Nestos River and its substreams. These sediments are formed of alternate sand, clay, and silt layers, thus reflecting a broad spectrum of formational and depositional conditions. These conditions led to a very diverse geological domain. Moreover, the presence of organic clay in some spots is due to delta marshes.
The study area is part of the wider Nestos River delta. This delta is one of the wetlands of utmost significance in Greece and Europe because of the size and the diversity of its ecology. This ecosystem is part of the Ramsar Convention on Wetlands of International Importance and the Natura 2000 network and is a special protection area of the European Union Directive on the Conservation of Wild Birds. It is included in the National Park of Eastern Macedonia and Thrace, which extends from the Nestos River delta to the lakes of Vistonida and Ismarida. Overall, the Nestos River delta consists of various individual wetland ecosystems. The greatest and most significant of these are (in alphabetically order): The Agiasmos Lagoon, Bassova Lagoon, Chaidou Lagoon, Eratinos Lagoon, Keramoti Lagoon, Kokala Lagoon, Magana Lagoon, and Monastiraki Lagoon.
The minimum and maximum annual rainfalls for the period 1985–2019 were 227.30 mm and 968.20 mm, respectively. Taking into account that the mean annual precipitation value was 496.80 mm over this period, the mean annual volume of rainfall received in the study area was 104.69 × 106 m3. The wettest and driest periods were November–December and July–September, respectively.
Two distinct units are present: (i) an unconfined aquifer and (ii) a confined aquifer system located in the Miocene sediment series. The former is primary recharged directly by precipitation, through percolation of the Nestos riverbed, and by lateral crossflows from the karst system located to the northern margin of the deltaic system. Lateral inflows contribute greatly to the recharge of the Miocene series aquifers. The unconfined aquifer, which is formed in the Quaternary deposits, is characterized by minor hydraulic head fluctuations.
The confined aquifer system is situated in the porous media of the delta (mainly sands and gravel which alternates with clays in vertical as well as in horizontal directions, forming a series of superimposed aquifers). Hence, under the unconfined aquifer, a series of confined aquifers exist up to a depth of at least 150 m. In the series of confined aquifer units, some artesian ones have been documented; these start from a shallow depth of about 15 m. In the Keramoti area, similar artesian units have been identified to a depth of approximately 120 m.

5. Materials and Methods

In the framework of this study, in situ measurements of groundwater temperature, pH, and electrical conductivity were performed at a monitoring network consisting of 24 and 22 wells in the unconfined and the confined aquifers (Figure 2), respectively, in four time periods (May and October 2019 and 2020). The monitoring periods correspond to the end of the wet and the dry seasons (May and October, respectively), thus are deemed representative of the extreme hydrological conditions in the region that are expected to reflect on the maximum potential water quality differentiations within a year. Groundwater samples were also collected from this network over the same periods and analyzed at the ELOT EN ISO 17025:2017-accredited laboratory of the Soil and Water Resources Institute of the Hellenic Agricultural Organization, Sindos, Greece, for the determination of 16 physicochemical and chemical parameters used in the calculation of PoS index: electrical conductivity (EC), NH4+, ΝO3, ΝO2, Cl, Na+, Ca2+, Mg2+, Mn2+, SO42−, and the dissolved forms of As, Pb, Cd, Cu, Fe, and B (Asdis, Pbdis, Cddis, Cudis, Fedis, and Bdis).
Samples were collected in 1 L polyethylene bottles thoroughly cleaned with mild glass detergent, subsequently washed with deionized water, and air-dried. Samples were collected from taps or valves at the wellhead of each sampling point following operation of the pump for 5 to 10 min in order to ensure stagnant water from the pump mains was discharged and a representative groundwater sample was collected. A visual check of turbidity, and stabilization of the pH and EC values measured at the discharged water preceded sample collection, following which, sampling bottles were washed twice with discharged water, filled to the top, tapped, capped, and registered in dedicated sampling logs. In order to speed up field operations, no in situ filtration or chemical preservation was performed. Samples were transported to the laboratory facilities on the same day of sampling and cooled in thermally insulated containers with ice packs. At the laboratory facilities, samples were immediately filtered and stored in refrigerators at 4 °C. The fraction of the sample used for heavy metal determinations was acidified with ultra-pure HNO3- to pH < 2 [27]. Determinations were concluded within 2 days of collection for the major ions and EC and less than 10 days for the heavy metals. A set of counter samples was preserved for a period of a month as reference in case odd results dictated repetition of determinations.
Electrical conductivity was measured by a CRISON GLP 32 Conductimeter according to the ELOT EN 27888:1993 method [28]. Potentiometric measurement of alkalinity, carbonates, and bicarbonates was carried out using volumetric determination, according to the ELOT EN ISO 9963-1:1996 method [29]. Nitrates, nitrites, and ammonium were determined colorimetrically using a Perkin Elmer Lambda 35 UV/VIS spectrophotometer according to Standard Methods for the Examination of Water and Wastewater Handbook, APHA, 23rd edition [27], while boron was also determined colorimetrically using a Perkin Elmer Lambda 35 UV/VIS spectrophotometer according to Standard Methods for the Examination of Water and Wastewater, 2017, 1-43: “Sample Storage and Preservation and Methods of Soil Analysis, Part 3: Chemical Methods”, 1996, and Chapter 21, “Azomethine-Hydrogen Method” [30]. Cl- was determined by the Mohr titration method [31]. Sodium was determined by Sherwood M410 Flame Photometer according to the ISO 9964-3:1993 method [32]. Calcium and magnesium were determined with a Perkin Elmer AAnalyst 700 Atomic Absorption Spectrometer according to Standard Methods for the Examination of Water and Wastewater, APHA, 3500, B, 23rd edition [27]. Trace elements and heavy metals were measured with a Perkin Elmer AAnalyst 700 Graphite Furnace Atomic Absorption Spectrometer according to the ELOT EN ISO 15586:2003 method [33], while mercury was determined with hydride generation atomic absorption spectrophotometry (HGAAS) [34].
Based on the results of the laboratory determinations, only negligible differences were noted between the two wet (May) and the two dry (October) sampling periods. On the contrary, as expected, changes were apparent between the dry and the wet periods. Hence, individual analysis and presentation of the results per sampling period would not add to the scope of the work, which comprises the identification and presentation of spatiotemporal variations in groundwater quality, focusing on possible discretization between the unconfined and the confined aquifers. Analysis of the results is, therefore, presented in the meaningful pattern of averaged values for the two wet (May 2019 and 2020) and the two dry (October 2019 and 2020) periods, as tabulated in Table 3, Table 4, Table 5 and Table 6 for the identified aquifers.
The LOD values for each determined parameter, following the adopted analytical methods, are: Na+-K+-Ca2+, 0.33 mg/L; Mg2+, 0.26 mg/L; NO2, 0.011 mg/L; NO3, 0.67 mg/L; NH4+, 0.033 mg/L; Cl, 0.5 mg/L; Bdis, 0.022 mg/L; SO42−, 1.80 mg/L; Mndis, 2.17 μg/L; Asdis, 0.33μg/L; Pbdis, 1.70 μg/L; Cddis, 0.1 μg/L; Cudis, 3.33 μg/L; Fedis, 3.33 μg/L. As clearly illustrated in the tabulated results of the hydrochemical analyses, a considerable number of determinations are rather low, being below the limit of detection (LOD).
Such data may not be subject to further arithmetic or statistical analysis in this form, and, in this form, should be disregarded from further evaluation. However, this approach would lead to an erroneous perception of missing data, whereas, in fact, the dataset clearly illustrates very low concentrations in specific parameters. Various practices are proposed in the literature to address the condition of very low determined concentrations, apart from discarding the data [35]. Suggested practices include use of the limit of quantification (LOQ) instead of the LOD and setting as arbitrary concentration the value of LOD or LOQ. Alternatively, any other fraction of LOQ or LOD are used, usually 50% [36,37]. In the course of the current study, 50% of the LOD value for each determined parameter and for the analytical method used was adopted.

6. Results and Discussion

PoS index is calculated for the above 24 and 22 wells from the unconfined and the confined aquifers, respectively. In the following tables, Table 7, Table 8, Table 9 and Table 10, the classification of the results are presented in accordance with the PoS index.
The following figures, Figure 3, Figure 4, Figure 5 and Figure 6, present the spatial distribution of the PoS Index for the unconfined and the confined aquifers.
Calculated PoS indices were used to study the spatiotemporal variations of the groundwater quality characteristics and depict the key parameters relating to the assessed conditions that could potentially cause adverse effects to humans and the natural environment. The following main observations can be made from the above tables (Table 7, Table 8, Table 9 and Table 10) and the spatial distribution depicted in Figure 3, Figure 4, Figure 5 and Figure 6.
  • Unconfined aquifer
In the wet period of May, of the 22 total sampling points, 13 fall in the high-level degradation category, and only one sampling point (no. 42) falls into the very high-quality degradation category. In October, however, only 6 wells of the 22 belong in the high-level quality degradation category. Only samples of one well (no. 34) steadily present a low–none degradation level for both periods. In both periods, the main degradation factors are As and Mn. Water quality based on PoS appears higher during the wet period (May) compared to the dry (October) period (Table 11).
  • Confined aquifer
The overall picture presented is that the majority of the tested samples belong to the low and moderate degradation class for both periods (May and October). In the wet period (May), two samples (sampling points 2 and 15) demonstrate an exceptionally high level of degradation, two samples (sampling points 6 and 14) belong to very high degradation category, and two samples (sampling points 7 and 14) demonstrate a high degradation of quality. In the dry period of October, two samples (sampling points 2 and 15) belong to a high-degradation category, and three samples (sampling points 6, 7, and 14) belong to a high-degradation category. As in the unconfined aquifer, the main degradation factors are As and Mn, and a similar trend of improved overall quality status during the dry period of October is noted, although it is not pronounced (Table 12).
As already stated, observing the distribution of the dominant factors (Qfs) among the six PoS index classes (Table 11 and Table 12), it is clear that As and Mn dominate (classes 2–6 for the two aquifers). Also, the other factors that affect the PoS index for classes 1 and 2 are NO3, NH4, and Na. Each substance, regardless of the P-class to which it belongs, may affect, either weakly or strongly, the overall groundwater quality status. Hence, PoS index is not biased towards toxic substances which may pose adverse effects to the environment (e.g., Pb, As, Hg) but also considers the major contribution of less toxic substances (e.g., Na, NO3, NH4, Mn) to the control of sample’s hydrogeochemistry and, consequently, to its overall environmental characteristics.
The cumulative contribution (%) of each of the parameters considered in the calculation of the PoS index for each of the periods (wet and dry) and the unconfined and confined aquifer is reported as derived from the following equation, Equation (4):
Contribution   ( % ) = [ Σ Q f i Σ Q f i + Σ Q f j + + Σ Q f n ] × 100
where ΣQfi is the total coefficient of the parameter i; ΣQfj is the total coefficient of the parameter j; and ΣQfn is the total coefficient of the parameter n.
The diagrams in Figure 7 and Figure 8 indicate that, for the unconfined aquifer, the main driving factor for groundwater quality degradation is Mn, at 68% for the wet period and approximately at 53% for the dry period. For the same aquifer, As is the second most important driving factor, at 20% and approximately 32% for the wet and dry seasons, respectively.
The diagrams in Figure 9 and Figure 10 indicate that, for the confined aquifer, the main driving factor for groundwater quality degradation is As, at about 56% for the wet and dry periods, respectively, followed by Mn, at about 32% and 28% for the wet and dry periods, respectively. The percentage of contribution for each of the other parameters for both aquifers and periods is less than 10%.
Analysis of the results (Table 11 and Table 12, Figure 7, Figure 8, Figure 9 and Figure 10) and their spatial distribution clearly suggest that the controlling factors of groundwater quality are As and Mn, which show a slight variability between the wet and the dry sampling periods and between the confined system and the unconfined aquifer. In both aquifers, the most degraded quality characteristics appear around the villages of Chrisochori and, secondarily, Chrisoupoli, which coincides well with the documented geothermal field of Nestos basin [38]. Hence, the elevated concentrations of these parameters may well be linked to the geological structure and the low enthalpy fields documented, as suggested in several cases of such fields across the world [39,40,41]. As already discussed, in applying the PoS methodology, a set of reference values needs to be considered, which, in our case, is the potability set of criteria, as practiced in many analogous applications. The origin of these parameters can clearly be traced back to natural geogenic processes. As such, the quality status of the assessed water samples may not be characterized as environmentally poor, but these samples indeed demonstrate impaired potability. Assuming the environmental quality characterization, the thresholds set in the implementation of the method for the particular parameters could have been shifted to higher values to account for the natural background concentrations due to the geothermal fields. Had this been the case, the implementation of the method would have not pinpointed these parameters at all. In fact, this is the rationale adopted by the EU in the framework of groundwater quality characterization, as the Water Framework Directive (2000/60/EEC) [42] focuses solely on the identification of human-induced elevated concentrations. However, this approach, which indeed seeks to assess deviations of groundwater quality from natural conditions, may fail to identify quality issues that are of essence in particular uses. As and Mn concentrations in the unconfined aquifer are indeed elevated but, in general, are lower than in the confined aquifer system, which could be argued to be the result of mixing of groundwater of the unconfined aquifer with surface water from the Nestos River that is mainly used for irrigation.
Even though nitrogen compounds are not the dominant controlling factors of groundwater quality in the assessed samples, they do contribute to its degradation. Especially in the unconfined aquifer, nitrogen compounds in the forms of NO3- and NO4+ are apparent, and rather scattered, as reported in recent research works. Nitrogen concentrations are higher in the wet season compared to the dry season, and this may be attributed mainly to the following reasons: (a) during the wet season, the residues of fertilizers in the unsaturated zone are being leached down to the saturated zone due to the recharge the system receives, (b) parts of the unsaturated zone become saturated as the groundwater levels elevate in response to recharge, (c) typically, basic fertilization is practiced in late winter–early spring time, thus, in conjunction with the practiced over-irrigation, nitrogen leaches into the saturated zone. The confined aquifer system also exhibits nitrogen concentrations, but much lower, influencing the quality of the system to a much lower degree. The existence, however, of nitrogen compounds in the confined aquifer proves that the two systems are not entirely secluded by an aquiclude, which is a typical hydrogeological feature reported in other alluvial basins. Moreover, this finding may support the hypothesis of direct hydraulic interaction of the two aquifer systems through the production wells installed in the basin which tap all aquifer units drilled. Due to this phenomenon and the construction characteristics of the wells, it has been estimated that the preferential flow paths that develop through the gravel packs of the production wells may increase the vertical vector of the hydraulic conductivity of the confining geological strata by 2–3 orders of magnitude per basin—about the same scale as the study area.
Quality issues appear to be more concentrated in the confined aquifer system, especially around the Chrisochori–Chrisoupoli area, as already mentioned, that coincides with the geothermal field. On the contrary, the unconfined aquifer exhibits more scattered quality degradation issues that cover a much wider area, and these are, as already mentioned, related to nitrogen compounds. Clearly, nitrogen is an indicator of agricultural activity, the main economic activity in the region, as Nestos is amongst the high agricultural productivity basins of the country.
In contrast to conventional hydrochemical analyses concluded to date for the basin, salinization-related parameters do not seem to control or even strongly influence groundwater quality. Electric conductivity, or other conventionally considered parameters’ concentrations, such as Na+, for example, may often lead to misleading results or mask the quality-controlling factors, as these are clearly revealed in the studied case.

7. Conclusions

Accurate and reliable quality assessment is of essence in effective management of groundwater resources. As quick decisions need to be made, often by non-experts. understanding the key quality issues in a conclusive yet easy to conceive manner is necessary. To this end, index-based quality assessment methods offer a valuable tool. Amongst the numerous such methods proposed in the literature, PoS is a relatively new one that builds on the philosophy of quality assessments that are easy to use and comprehend, filtering out details that are essential and preserving accuracy.
Since, as its background, the method references toxicity studies listed in the literature and the quality standards for water intended for human consumption, it needs to be applied with great conscience and consideration with regards to the term “quality degradation level”. Hence, using unanimously accepted reference levels of a specific water use, the method yields a classification system that, however, does not necessarily imply degraded water quality “senso stricta”, but a potential harmful effect to biotic end users. It therefore follows that a critical view is always needed in carefully interpreting the results, thus avoiding misinterpretations. To this end, the term quality degradation level originally adopted by the method may be substituted with “quality category” or “level/class”, thus avoiding connotations that may be misleading. Indeed, the merit of the method compared to others, such as the chemical quality classification and assessment schemes for groundwater in the context of the 2000/60 EC Water Framework Directive [43]; 2006/118 EC Directives [44]; and 1st revision of the Thracian Water District (EL 12) [45], “Characterization and evaluation of the status of the groundwater bodies of the Thracian Water District (EL 12), Nestos Delta Groundwater System (EL 1200060)”, is that it employs a standardized approach regardless of the natural background levels for each constituent in each considered groundwater system, enabling a direct comparison of quality characteristics.
The results of the PoS index provide a new perspective to groundwater quality assessment. The key points of the deduced results are in line with the results of typical hydrochemical aproaches such as isoparametric spatial distribution maps, statistical analyses, and linear diagrams. PoS visualizes results in discrete classes and subclasses of quality degradation, thus offering a detailed assessment which does not lack spatial and temporal resolution. Highlighting the key parameters that control and shape the quality characteristics of the examined samples, the method helps focusing on the potential origin of the findings, these being either geogenic or manmade. In the studied case, it was deduced that the key quality issues relate to geothermal fields (As, Mn) and agricultural activities (nitrogen compounds), ruling out seawater intrusion as a key mechanism. Seawater intrusion could mistakenly have been assumed to be a key mechanism if not carefully interpreted through the study of conventional hydrochemical analysis methods. Based on this assessment, pressure analysis for a given catchment may easily confirm the potential origin of the quality degradation. As no actual restriction occurs to the parameters considered for the assessment, the method is adequately versatile to be used for comparisons between different aquifers of the same system, and of systems of different catchments, over various periods of time. Of course, any index-based assessment needs to be accompanied by in-depth analysis and study of the identified quality issues, employing conventional hydrochemical studies, statistical approaches, and hydrochemical diagrams, amongst others, to reach conclusive results and detailed, deep understanding of the dominant hydrogeochemical evolution mechanisms. Such analyses are indispensable, comprehended by experts, and constitute the full justification for groundwater quality assessment in support of the easy-to-comprehend PoS results.

Author Contributions

Conceptualization, G.K., F.-K.P. and A.P.; Methodology, G.K., A.P., F.-K.P. and V.K.; Validation, A.P., I.G. and F.-K.P.; Writing – original draft, G.K., A.P., I.G. and C.P.; Writing – review & editing, G.K., A.P., I.G., C.P. and F.-K.P.; Supervision, A.P., I.G. and I.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The results of the paper are part of the research of the doctoral dissertation (last stage of PhD) of George Kampas, Department of Civil Engineering, School of Engineering, Democritus University of Thrace, Greece.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kim, J.; Kim, R.; Lee, J.; Chang, H. Hydrogeochemical characterization of major factors affecting the quality of shallow groundwater in the coastal area at Kimje in South Korea. Environ. Geol. 2003, 44, 189–478. [Google Scholar] [CrossRef]
  2. Zhu, G.F.; Li, Z.Z.; Su, Y.H.; Ma, J.Z.; Zhang, Y.Y. Hydrogeochemical and isotope evidence of groundwater and recharge in Minqin Basin, Northwest China. J. Hydrol. 2007, 333, 239–251. [Google Scholar] [CrossRef]
  3. Sandeep, K.G.; Tziritis, E.; Sudhir, K.S.; Jayant, K.T.; Abhay, K.S. Environmental monitoring of water resources with the use of PoS index: A case study from Subarnarekha River basin, India. Environ. Earth Sci. 2018, 77, 70. [Google Scholar] [CrossRef]
  4. Piper, A.M. A graphic procedure in the geochemical interpretation of water analyses. Eos Trans. Am. Geophys. Union 1994, 25, 914–923. [Google Scholar]
  5. Durov, S.A. Natural water and graphical representation of their composition. Dokl. Akad. Nauk USSR 1948, 59, 87–90. [Google Scholar]
  6. Wilcox, L.V. Classification and Use of Irrigation Waters; Circular No. 969; U.S. Department of Agriculture: Washington, DC, USA, 1955; 19p. [Google Scholar]
  7. Revelle, R. Criteria for recognition of sea water in ground waters. Trans. Am. Geophys. Union 1941, 22, 593–597. [Google Scholar]
  8. Sivakumar, S.; Chandrasekaran, A.; Ravisankar, R.; Ravikumar, S.M.; Jebakumar, J.P.P.; Vijayagopal, P. Rements of natural radioactivity and evaluation of radiation hazards in coastal sediments of east coast of Tamilnadu, India using statistical approach. J. Taibah Univ. Sci. 2014, 8, 375–384. [Google Scholar] [CrossRef]
  9. McKenna, J.E. An enhanced cluster analysis program with bootstrap significance testing for ecological community analysis. Environ. Model. Softw. 2003, 18, 205–220. [Google Scholar] [CrossRef]
  10. Tziritis, E.; Panagopoulos, A.; Arampatzis, G. Development of an operational index of water quality (PoS) as a versatile tool to assist groundwater resources management and strategic planning. J. Hydrol. 2014, 517, 339–350. [Google Scholar] [CrossRef]
  11. Tziritis, E.; Arampatzis, G.; Hatzigiannakis, E.; Panoras, G.; Panoras, A.; Panagopoulos, A. Quality characteristics and hydrogeochemistry of irrigation waters from three major olive groves in Greece. Desalination Water Treat. 2016, 57, 11582–11591. [Google Scholar] [CrossRef]
  12. Tziritis, E.; Skordas, K.; Kelepertsis, A. The use of hydrochemical analyses and multivariate statistics for the characterization of groundwater resourses in the complex aquifer system. A case study in Amyros River basin, Thessaly, central Greece. Environ. Earth Sci. 2016, 75, 1–11. [Google Scholar] [CrossRef]
  13. Salcedo-Sánchez, E.R.; Garrido Hoyos, S.E.; Esteller Alberich, M.V.; Martínez Morales, M. Application of water quality index to evaluate groundwater quality (temporal and spatial variation) of an intensively exploited aquifer (Puebla Valley, Mexico). Environ. Monit. Assess. 2016, 188, 1–20. [Google Scholar] [CrossRef]
  14. Fernández, N.; Solano, F. Índices de Calidad y de Contaminación del Agua: Technical Report; Universidad de Pamplona: Norte de Santander, Colombia, 2005. (In Spanish) [Google Scholar]
  15. Lumb, A.; Sharma, T.C.; Bibeault, J.-F.; Klawunn, P. A Comparative Study of USA and Canadian Water Quality Index Models. Water Qual. Expo. Health 2012, 3, 203–216. [Google Scholar] [CrossRef]
  16. CCME. Canadian Environmental Quality Guidelines for the Protection of Aquatic Life; CCME water quality index: Technical report, 1.0; CCME: Winnipeg, MB, Canada, 2001. [Google Scholar]
  17. Basha, U.I.; Rajasekhar, M.; Ghosh, S. Spatial assessment of groundwater quality using CCME-WQI and hydrochemical indices: A case study from Talupula Mandal, Ananthapuramu district, South India. Appl. Water Sci. 2022, 12, 168. [Google Scholar] [CrossRef]
  18. Agency for Toxic Substances and Disease Registry (ATSDR). Detailed Data Table for the 2011 Priority List of Hazardous Substances. Agency for Toxic Substances and Disease Registry. 2011. Available online: https://stacks.cdc.gov/view/cdc/26622 (accessed on 20 November 2011).
  19. Soltani, S.; Moghaddam, A.A.; Barzegar, R.; Kazemian, N.; Tziritis, E. Hydrogeochemistry and water quality of the Kordkandi-Duzduzan plain, NW Iran: Application of multivariate statistical analysis and PoS index. Environ. Monit. Assess. 2017, 189, 455. [Google Scholar] [CrossRef] [PubMed]
  20. Vrouchakis, I. Hydrogeological Research Using Modern Management Tools in the Context of Qualitative and Quantitative Monitoring of the Aquifer System of the Tyrnavos Subbasin (Thessaly). Ph.D. Thesis, Athens Agronomist University, Athens, Greece, 2022. (In Greek). [Google Scholar]
  21. Valko, M.M.H.C.M.; Morris, H.; Cronin, M.T.D. Metals, toxicity and oxidative stress. Curr. Med. Chem. 2005, 12, 1161–1208. [Google Scholar] [CrossRef] [PubMed]
  22. International Toxicity Estimates for Risk (ITER). Toxicology Data Network. International Toxicity Estimates for Risk. 2011. Available online: https://pubmed.ncbi.nlm.nih.gov/15760833/ (accessed on 30 November 2012).
  23. World Health Organization (WHO). Guidelines for Drinking-Water Quality, 4th ed.; WHO: Geneva, Switzerland, 2011.
  24. Allen, M.J.; Yen, W.M. Introduction to Measurement Theory, printing ed.; Waveland Press: Long Grove, IL, USA, 2002. [Google Scholar]
  25. Directive 2020/2184/EC; European Directive on the Quality of Water Intended for Human Consumption. Publication Office of the European Union: Luxembourg, 2020.
  26. Delimani, P. Geological Changes of the Coastline in the Thrace Region and Impact on the Land Use of the Coastal Zone. Ph.D. Thesis, Department of Civil Engineering, Democritus University of Thrace, Xanthi, Greece, 2000. (In Greek). [Google Scholar]
  27. American Public Health Association. APHA Standard Methods for The Examination of Water and Wastewater, 23rd ed.; American Public Health Association (APHA): Washington, DC, USA, 2017. [Google Scholar]
  28. Novotny, V. Water Quality: Prevention, Identification and Management of Diffuse Pollution; Van Nostrand-Reinhold Publishers: New York, NY, USA, 1994; ISBN 0442005598. [Google Scholar]
  29. DIN EN ISO 9963-1:1996-02; Wasserbeschaffenheit—Bestimmung der Alkalinität—Teil 1: Bestimmung der Gesamten und der Zusammengesetzten Alkalinität (ISO 9963-1:1994); Deutsche Fassung EN ISO 9963-1:1995. DIN Standards Committee Water Practice (NAW): Berlin, Germany, 1995; p. 24.
  30. Keeney, D.R.; Nelson, D.W. Methods of Soil Analysis; Page, A.L., Ed.; Agronomy Monographs; American Society of Agronomy, Soil Science Society of America: Madison, WI, USA, 1983; ISBN 9780891189770. [Google Scholar]
  31. Belcher, R.; Macdonald, A.M.G.; Parry, E. On mohr’s method for the determination of chlorides. Anal. Chim. Acta 1957, 16, 524–529. [Google Scholar] [CrossRef]
  32. ISO 9964-3:1993; Water Quality—Determination of Sodium and Potassium—Part 3: Determination of Sodium and Potassium by Flame Emission Spectrometry. ISO: Geneva, Switzerland, 1993.
  33. ISO 15586:2003; Water Quality—Determination of Trace Elements Using Atomic Absorption Spectrometry with Graphite Furnace [Online]. ISO: Geneva, Switzerland, 2003. Available online: https://www.iso.org/standard/38111.html (accessed on 13 May 2018).
  34. Mohammed, E.; Mohammed, T.; Mohammed, A. Optimization of instrument conditions for the analysis for mercury, arsenic, antimony and selenium by atomic absorption spectroscopy. MethodsX 2018, 5, 824–833. [Google Scholar] [CrossRef]
  35. Nehls, G.J.; Akland, G.G. Procedures for handling aerometric data. J. Air Pollut. Control Assoc. 1973, 23, 180–184. [Google Scholar] [CrossRef]
  36. Clarke, J.U. Evaluation of censored data methods to allow statistical comparisons among very small samples with below detection limit observations. Environ. Sci. Technol. 1998, 32, 177–183. [Google Scholar] [CrossRef]
  37. Hornung, R.W.; Reed, L.D. Estimation of average concentration in the presence of nondetectable values. Appl. Occup. Environ. Hyg. 1990, 5, 46–51. [Google Scholar] [CrossRef]
  38. Kolios, N.; Koutsino, S.; Arvanitis, A.; Rydakis, G. Geothermal Situation in Northeastern Greece. In Proceedings of the World Geothermal Congress, Antalya, Turkey, 24–29 April 2005. [Google Scholar]
  39. Davraz, A.; Nalbantçılar, M.T.; Varol, S.; Önden, İ. Hydrogeochemistry and reservoir characterization of the Konya geothermal elds, Central Anatolia/Turkey. Geochemistry 2022, 82, 125867. [Google Scholar] [CrossRef]
  40. Morales-Simfors, N.; Bundschuh, J. Arsenic-rich geothermal fluids as environmentally hazardous materials—A global assessment. Sci. Total Environ. 2022, 817, 152669. [Google Scholar] [CrossRef] [PubMed]
  41. Kaasalainen, H.; Stefánsson, A.; Giroud, N.; Arnórsson, S. The geochemistry of trace elements in geothermal fluids, Iceland. Appl. Geochem. 2015, 62, 207–223. [Google Scholar] [CrossRef]
  42. Directive 2000/60/EC; European Framework Directive, Guidance Document 18: Guidance on Ground Water Status and Trend Assessment. European Commission: Luxembourg, 2009.
  43. Directive 2000/60/EC; European Directive on the Establishing a Framework for Community Action in the Field of Water Policy. European Commission: Luxembourg, 2000.
  44. Directive 2006/118/EC; European Directive on the Protection of Groundwater against Pollution and Deterioration. European Commission: Luxembourg, 2006.
  45. 1st Review of the Thracian Water District (EL 12)—National Gazette of Greece No 4680B. 2017. Available online: https://wfdver.ypeka.gr/el/consultation-gr/1revision-consultation-gr/consultation-1revision-el12-gr/ (accessed on 13 May 2018).
Figure 1. Study area.
Figure 1. Study area.
Water 16 00352 g001
Figure 2. Groundwater monitoring network points for the unconfined and the confined aquifers.
Figure 2. Groundwater monitoring network points for the unconfined and the confined aquifers.
Water 16 00352 g002
Figure 3. Spatial distribution of PoS Index for the unconfined aquifer (average May 2019 and 2020 values).
Figure 3. Spatial distribution of PoS Index for the unconfined aquifer (average May 2019 and 2020 values).
Water 16 00352 g003
Figure 4. Spatial distribution of PoS Index for the unconfined aquifer (average October 2019 and 2020 values).
Figure 4. Spatial distribution of PoS Index for the unconfined aquifer (average October 2019 and 2020 values).
Water 16 00352 g004
Figure 5. Spatial distribution of PoS Index for the confined aquifer (average May 2019 and 2020 values).
Figure 5. Spatial distribution of PoS Index for the confined aquifer (average May 2019 and 2020 values).
Water 16 00352 g005
Figure 6. Spatial distribution of PoS Index for the confined aquifer (average October 2019 and 2020 values).
Figure 6. Spatial distribution of PoS Index for the confined aquifer (average October 2019 and 2020 values).
Water 16 00352 g006
Figure 7. Cumulative contribution of the presented parameters in applying the PoS index for all analyzed samples during the wet season (May 2019 and 2020) for the unconfined aquifer.
Figure 7. Cumulative contribution of the presented parameters in applying the PoS index for all analyzed samples during the wet season (May 2019 and 2020) for the unconfined aquifer.
Water 16 00352 g007
Figure 8. Cumulative contribution of the presented parameters in applying the PoS index for all analyzed samples during the dry season (October 2019 and 2020) for the unconfined aquifer.
Figure 8. Cumulative contribution of the presented parameters in applying the PoS index for all analyzed samples during the dry season (October 2019 and 2020) for the unconfined aquifer.
Water 16 00352 g008
Figure 9. Cumulative contribution of the presented parameters in applying the PoS index for all analyzed samples during the wet season (May 2019 and 2020) for the confined aquifer.
Figure 9. Cumulative contribution of the presented parameters in applying the PoS index for all analyzed samples during the wet season (May 2019 and 2020) for the confined aquifer.
Water 16 00352 g009
Figure 10. Cumulative contribution of the presented parameters in applying the PoS index for all analyzed samples during the dry season (October 2019 and 2020) for the confined aquifer.
Figure 10. Cumulative contribution of the presented parameters in applying the PoS index for all analyzed samples during the dry season (October 2019 and 2020) for the confined aquifer.
Water 16 00352 g010
Table 1. Parameters for calculating the PoS index and classification of these parameters based on their toxicity, individual score, and weighting factor.
Table 1. Parameters for calculating the PoS index and classification of these parameters based on their toxicity, individual score, and weighting factor.
Toxic ClassSubstancesPointsP-ClassPartial Score (PS)Weighting Factor
(Wf)
1As1VI100.19
2Pb10V80.16
Cd10V80.16
3NH4100IV50.10
NO2100IV50.10
4NO31000III30.06
B5000II1.50.03
Cu5000II1.50.03
Mn5000II1.50.03
Ca50,000I10.02
Cl50,000I10.02
Fe50,000I10.02
K50,000I10.02
Mg50,000I10.02
Na50,000I10.02
SO450,000I10.02
EC50,000I10.02
Table 2. Presentation of the PoS indicator (range of categories, subcategories, and qualitative degradation color rendering).
Table 2. Presentation of the PoS indicator (range of categories, subcategories, and qualitative degradation color rendering).
PoS
Class
RangeMinMaxSubclass RangeabcQuality Degradation LevelColor
1a ≤ t0150-<tNone–Low
2t < a < 2t15030050150200200250250300Low
32t < a < 4t300600100300400400500500600Moderate
44t < a < 8t6001200200600800800100010001200High
58t < a < 16t12002400400120016001600200020002400Very High
6a > 16t2400---Severe
Table 3. Analytical results of groundwater determinations for unconfined aquifer (average values for May 2019 and 2020).
Table 3. Analytical results of groundwater determinations for unconfined aquifer (average values for May 2019 and 2020).
Well
Code
AsdisPbdisCddisNH4+NO2NO3BdisCudisMndisCa2+ClFedisMg2+Na+SO42−EC
μg/Lμg/Lμg/Lmg/Lmg/Lmg/Lmg/Lμg/Lμg/Lmg/Lmg/Lμg/Lmg/Lmg/Lmg/LμS/cm
251.00<LOD<LOD<LOD0.047.360.04<LOD844.7084.4010.5059.1811.889.0522.58546
2625.35<LOD0.190.790.0717.84<LOD<LOD1016.1076.4814.886.2716.5316.7645.04549.5
276.43<LOD0.310.90<LOD<LOD<LOD<LOD939.68114.3010.9173.1416.3322.3578.04732.5
288.00<LOD0.210.060.4117.340.14<LOD784.54132.3857.088.4239.5075.88122.651118.5
29<LOD<LOD<LOD0.301.2411.860.08<LOD791.3581.0813.28156.3516.3823.7817.85635
30<LOD<LOD<LOD<LOD0.255.330.11<LOD1262.3580.0513.73128.3223.6038.8159.76722.5
31<LOD<LOD0.460.070.238.21<LOD<LOD1347.7588.459.8713.8011.8011.5232.40645.5
321.34<LOD0.230.080.259.57<LOD<LOD797.2681.237.646.488.908.163.52422
33<LOD4.28<LOD4.90<LOD30.83<LOD<LOD20.7785.4695.05302.7011.6381.9871.27837
34<LOD<LOD<LOD0.11<LOD37.27<LOD<LOD25.46126.0055.6015.4210.2534.4877.76848
356.58<LOD0.244.110.02<LOD0.28<LOD760.8864.65126.2926.7642.13127.625.971086.5
367.08<LOD<LOD0.520.085.64<LOD<LOD1394.2595.2833.4432.7218.5828.7265.04717
377.98<LOD0.44<LOD0.4313.17<LOD<LOD369.5455.856.56<LOD6.807.06<LOD335
388.99<LOD<LOD0.340.02<LOD0.05<LOD1346.00100.6326.0127.8420.9349.5778.04749
394.45<LOD<LOD<LOD<LOD<LOD<LOD<LOD22.67106.73399.1119.3029.10182.9318.881728
4016.88<LOD<LOD2.89<LOD1.220.39<LOD964.5084.15268.3842.1827.50364.82<LOD2310
416.23<LOD0.195.30<LOD5.110.5212.63199.8640.4960.82769.9015.20155.56<LOD986.5
4245.58<LOD<LOD1.660.02<LOD<LOD<LOD2345.6067.486.89102.028.4810.68<LOD416
4310.98<LOD0.311.37<LOD<LOD<LOD<LOD1277.3073.339.19135.1011.789.44<LOD454
443.38<LOD0.270.67<LOD<LOD0.10<LOD981.0086.8325.1482.7519.4865.7974.67860.5
4510.09<LOD<LOD0.780.0612.500.42<LOD966.90105.4055.6012.1338.45134.16132.181206
464.94<LOD0.190.610.2612.020.07<LOD1314.25116.4524.9922.3030.8543.48109.64816
Note: LOD: Limit Of detection.
Table 4. Analytical results of groundwater determinations for unconfined aquifer (average values for October 2019 and 2020).
Table 4. Analytical results of groundwater determinations for unconfined aquifer (average values for October 2019 and 2020).
Well
Code
AsdisPbdisCddisNH4+NO2NO3BdisCudisMndisCa2+ClFedisMg2+Na+SO42−EC
μg/Lμg/Lμg/Lmg/Lmg/Lmg/Lmg/Lμg/Lμg/Lmg/Lmg/Lμg/Lmg/Lmg/Lmg/LμS/cm
255.11<LOD<LOD0.17<LOD5.09<LOD<LOD266.0097.639.3111.6811.309.1026.31553
2623.00<LOD<LOD0.960.0614.760.08<LOD387.5079.4813.529.5814.4815.2144.85571.5
276.21<LOD<LOD0.74<LOD<LOD<LOD<LOD709.50109.609.2515.3913.9517.2951.22653
2812.74<LOD<LOD0.070.1534.410.13<LOD315.7680.4336.554.8522.2566.1187.22884
294.37<LOD<LOD0.380.036.570.05<LOD263.4678.5810.5614.8212.5320.0824.04541.5
302.94<LOD<LOD0.140.1311.210.10<LOD896.7598.9514.6229.9124.6541.7564.67756
313.05<LOD<LOD<LOD0.3112.62<LOD<LOD712.00121.2510.5911.0010.9010.3942.13673.5
322.29<LOD<LOD0.120.207.29<LOD<LOD300.0689.5811.568.128.888.9032.49516
331.73<LOD<LOD10.120.06<LOD0.04<LOD17.7233.6485.81201.2412.8857.4383.40728
343.10<LOD0.27<LOD<LOD31.82<LOD<LOD7.00141.2339.5830.429.3837.8565.95830
358.39<LOD<LOD1.28<LOD<LOD0.74<LOD403.8873.5322.5138.6731.9043.4968.04704.5
362.84<LOD<LOD0.370.066.14<LOD<LOD466.75100.7511.2210.6518.9029.4665.40675.5
379.91<LOD<LOD0.17<LOD6.30<LOD<LOD138.5346.855.3735.855.656.878.31317
389.19<LOD<LOD0.34<LOD<LOD0.10<LOD848.1391.5015.0110.7322.1086.8767.90771
394.38<LOD<LOD<LOD<LOD<LOD0.09<LOD9.87154.43261.2732.9741.15268.7638.672295
4021.89<LOD<LOD2.040.023.110.27<LOD404.7664.55268.3827.8720.28238.585.971617
418.20<LOD<LOD2.92<LOD1.150.28<LOD180.8148.7826.12187.5711.8866.45<LOD626.5
4223.26<LOD<LOD1.41<LOD<LOD0.06<LOD732.3163.606.68119.717.8511.11<LOD390.5
438.12<LOD<LOD1.60<LOD<LOD0.04<LOD500.7573.885.7249.5011.931.00<LOD469
444.14<LOD<LOD0.77<LOD<LOD0.14<LOD939.80110.1818.8665.0121.8572.4974.95892.5
456.95<LOD<LOD0.11<LOD6.890.46<LOD313.25100.58130.7918.0138.83132.4097.591235
464.03<LOD<LOD0.230.0916.180.04<LOD760.38116.4519.6638.2029.5342.98120.55802.5
Note: LOD: Limit Of detection.
Table 5. Analytical results of groundwater determinations for confined aquifer (average values for May 2019 and 2020).
Table 5. Analytical results of groundwater determinations for confined aquifer (average values for May 2019 and 2020).
Well
code
AsdisPbdisCddisNH4+NO2NO3-BdisCudisMndisCa2+ClFedisMg2+Na+SO42−EC
μg/Lμg/Lμg/Lmg/Lmg/Lmg/Lmg/Lμg/Lμg/Lmg/Lmg/Lμg/Lmg/Lmg/Lmg/LμS/cm
19.72<LOD<LOD0.12<LOD<LOD0.50<LOD12.95189.63164.25189.632.41253.1952.951234.5
299.03<LOD<LOD2.84<LOD<LOD0.20<LOD756.85226.4842.87226.4818.4075.26<LOD741.5
34.80<LOD<LOD0.77<LOD<LOD0.16<LOD662.6927.4623.6527.4616.7095.96<LOD705.0
41.40<LOD<LOD0.50<LOD<LOD<LOD<LOD650.7649.0752.9849.0710.5536.592.79637.0
515.34<LOD<LOD0.75<LOD<LOD<LOD5.55830.1187.558.6187.558.4817.45<LOD412.0
671.51<LOD<LOD1.90<LOD<LOD0.10<LOD1346.95129.018.45129.0117.5546.30<LOD511.0
730.82<LOD<LOD3.02<LOD<LOD<LOD<LOD814.256.3033.076.305.0822.68<LOD553.5
80.97<LOD<LOD0.09<LOD<LOD0.12<LOD177.75<LOD7.47<LOD18.3044.003.15446.5
929.00<LOD<LOD0.930.02<LOD0.18<LOD838.207.204.277.2015.7331.65<LOD436.5
1019.03<LOD<LOD0.09<LOD<LOD0.52<LOD6.07116.2197.27116.211.61232.3841.761035.0
1114.72<LOD0.190.08<LOD<LOD0.3813.4019.7070.89164.6670.894.00249.6734.581172.5
123.56<LOD0.19<LOD<LOD<LOD0.16<LOD22.4352.64102.6452.645.05184.1730.13910.0
1326.10<LOD<LOD0.74<LOD<LOD0.09<LOD481.71373.8546.75373.8510.4348.932.70700.5
1438.68<LOD0.183.61<LOD1.070.18<LOD611.16202.3598.87202.3518.58142.08<LOD987.5
1595.70<LOD0.234.590.38<LOD0.57<LOD694.8130.53107.1930.5332.40232.79<LOD1424.5
1611.63<LOD0.200.21<LOD<LOD0.39<LOD18.3517.17113.8517.172.25219.664.521004.0
173.98<LOD0.18<LOD0.035.99<LOD<LOD465.936.4113.136.4111.0318.209.40435.0
1814.62<LOD0.171.96<LOD<LOD0.12<LOD782.25307.947.20307.949.6035.843.79464.5
192.72<LOD<LOD0.17<LOD<LOD0.41<LOD110.5927.00456.3427.0010.38399.6128.492127.5
2016.19<LOD0.190.05<LOD<LOD0.24<LOD27.0021.6924.5021.691.14160.95<LOD683.5
2115.83<LOD0.42<LOD<LOD<LOD0.47<LOD19.03<LOD115.05<LOD2.08222.8915.951046.0
22< LOD<LOD0.21<LOD<LOD142.83<LOD<LOD10.0317.2445.3917.2421.0364.4851.791166.0
239.72<LOD0.200.76<LOD55.69<LOD<LOD40.8155.26133.5555.2614.08121.2771.491085.0
24< LOD<LOD0.34<LOD0.0228.720.19<LOD624.4468.30825.4168.3053.18321.27102.221915.5
Note: LOD: Limit Of detection.
Table 6. Analytical results of groundwater determinations for confined aquifer (average values for October 2019 and 2020).
Table 6. Analytical results of groundwater determinations for confined aquifer (average values for October 2019 and 2020).
Well
Code
AsdisPbdisCddisNH4+NO2-NO3-BdisCudisMndisCa2+Cl-FedisMg2+Na+SO42-EC
μg/Lμg/Lμg/Lmg/Lmg/Lmg/Lmg/Lμg/lμg/Lmg/Lmg/Lμg/Lmg/Lmg/Lmg/LμS/cm
18.34<LOD<LOD<LOD<LOD<LOD0.37<LOD8.2443.7393.3043.733.17188.9133.67977.0
263.52<LOD<LOD2.67<LOD<LOD0.21<LOD567.5077.3645.7677.3618.5579.05<LOD747.0
36.08<LOD<LOD0.75<LOD2.840.16<LOD499.2511.2421.7811.2415.3587.08<LOD679.0
42.44<LOD<LOD0.61<LOD<LOD<LOD<LOD295.9214.8551.3514.859.2544.257.06638.0
512.99<LOD<LOD1.23<LOD1.88<LOD<LOD400.245.037.645.037.5515.01<LOD425.5
62.42<LOD<LOD0.65<LOD1.460.05<LOD1006.2510.927.6710.9215.9348.22<LOD477.0
721.11<LOD<LOD2.83<LOD3.60<LOD<LOD314.5741.9132.2841.918.8524.86<LOD547.5
81.38<LOD<LOD0.08<LOD<LOD0.14<LOD70.163.507.613.5016.0341.26<LOD447.0
926.09<LOD<LOD1.01<LOD<LOD0.10<LOD656.7521.604.2521.6014.3036.29<LOD434.0
1016.85<LOD<LOD0.12<LOD<LOD0.50<LOD<LOD50.6792.7350.671.34213.7642.311024.0
1113.18<LOD<LOD0.06<LOD<LOD0.38<LOD13.7422.86159.8622.864.40223.7132.951199.5
123.14<LOD<LOD<LOD<LOD<LOD0.15<LOD<LOD20.89100.4920.894.40159.0928.58908.0
1312.37<LOD<LOD0.50<LOD<LOD0.04<LOD105.7035.4740.8635.479.4863.967.06704.0
1429.25<LOD<LOD1.93<LOD<LOD0.38<LOD425.38435.6094.67435.6016.30154.935.40998.5
1556.23<LOD<LOD5.11<LOD1.060.76<LOD323.1637.83233.5137.8338.10318.16<LOD1437.5
1611.78<LOD<LOD<LOD<LOD<LOD0.46<LOD12.86345.00103.05345.0017.80208.9931.401021.0
173.75<LOD<LOD0.050.032.75<LOD<LOD128.8324.966.4324.967.7311.4511.13362.7
1810.36<LOD<LOD1.83<LOD<LOD0.14<LOD421.8872.956.4872.958.9337.85<LOD464.5
192.02<LOD<LOD0.11<LOD<LOD0.41<LOD26.6620.68337.7020.689.60434.8845.672200.0
2014.912.60<LOD<LOD<LOD<LOD0.25<LOD6.1335.7720.4835.770.74147.41<LOD666.0
21
221.30<LOD0.21<LOD<LOD133.65<LOD<LOD13.8615.8552.9415.8520.6351.9995.491182.5
239.00<LOD<LOD<LOD<LOD57.49<LOD<LOD34.0030.27101.2530.2713.50129.4085.491141.5
241.50<LOD<LOD0.090.0926.200.20<LOD411.3823.30590.3323.3051.88374.84125.222820.0
Note: LOD: Limit Of detection.
Table 7. Groundwater classification according to the PoS index (average May 2019 and 2020 values) (unconfined aquifer).
Table 7. Groundwater classification according to the PoS index (average May 2019 and 2020 values) (unconfined aquifer).
Well CodePoS IndexSample ClassificationQuality Degradation Level
ClassColor
255463c/Mndis Moderate
2611524c/Asdis, Mndis High
277364a/Mndis High
287114a/Asdis, Mndis High
295413c/Mndis Moderate
307904a/Mndis High
318394b/Mndis High
325323c/Mndis Moderate
332362b/NH4+ Low
34981/NO3- None–Low
357224a/Mndis High
3610024c/Mndis High
374203b/Mndis Moderate
3810054c/Mndis High
391862a/Asdis Low
4010544c/Mndis High
414543b/Asdis, NH4+, Mndis Moderate
4223055c/Asdis, Mndis Very High
4310194c/Asdis, Mndis High
447044a/Mndis High
458654b/Asdis, Mndis High
469464b/Mndis High
Table 8. Groundwater classification according to the PoS index (average October 2019 and 2020 values) (unconfined aquifer).
Table 8. Groundwater classification according to the PoS index (average October 2019 and 2020 values) (unconfined aquifer).
Well CodePoS IndexSample Classification Quality Degradation Level
ClassColor
252862c/Asdis, Mndis Low
267364a/Asdis, Mndis High
275773c/Asdis, Mndis Moderate
285193c/Asdis, Mndis Moderate
292782c/Asdis, Mndis Low
306384a/Mndis High
315223c/Mndis Moderate
322552c/Mndis Low
332882c/NH4+ Low
341471/As, NO3- None–Low
354833b/Asdis, Mndis Moderate
363743a/Mndis Moderate
372972c/Asdis, Mndis Low
387214a/Asdis, Mndis High
391952a/Asdis Low
407854a/Asdis, Mndis High
413653a/Asdis, Mndis Moderate
429294b/Asdis, Mndis High
435013c/Asdis, Mndis Moderate
446914a/Mndis High
454073b/Asdis, Mndis Moderate
465943c/Mndis Moderate
Table 9. Groundwater classification according to the PoS index (average May 2019 and 2020 values) (confined aquifer).
Table 9. Groundwater classification according to the PoS index (average May 2019 and 2020 values) (confined aquifer).
Sample CodePoS IndexSample ClassificationQuality Degradation Level
ClassColor
12782c/Asdis Low
224686/Asdis Severe
35263c/Mndis Moderate
44483b/Mndis Moderate
58164b/Asdis, Mndis High
622425c/Asdis, Mndis Very High
711564c/Asdis, Mndis High
81542a/Mndis Low
910924c/Asdis, Mndis High
104403b/Asdis Moderate
113713a/Asdis Moderate
121341/Asdis None–Low
138524b/Asdis, Mndis High
1412435a/Asdis, Mndis Very High
1524436/Asdis Severe
162992a/Asdis Low
173803a/Asdis, Mndis Moderate
188244b/Asdis, Mndis High
192342b/Asdis, Mndis Low
203693a/Asdis Moderate
213873a/Asdis Moderate
222212b/NO3- Low
233513a/Asdis Moderate
245643c/Mndis Moderate
Table 10. Groundwater classification according to the PoS index (average October 2019 and 2020 values) (confined aquifer).
Table 10. Groundwater classification according to the PoS index (average October 2019 and 2020 values) (confined aquifer).
Sample CodePoS IndexSample ClassificationQuality Degradation Level
ClassColor
12182b/Asdis Low
216575b/Asdis Very high
34573b/Asdis, Mndis Moderate
42582c/Mndis Low
55263c/Asdis, Mndis Moderate
66684a/Mndis High
76754a/Asdis, Mndis High
8901/Asdis, Mndis None–low
99294b/Asdis, Mndis High
103873a/Asdis Moderate
113263a/Asdis Moderate
122962c/Asdis Low
133433a/Asdis Moderate
149344b/Asdis, Mndis High
1514845a/Asdis Very high
163183a/Asdis Moderate
171682a/Asdis, Mndis Low
185083c/Asdis, Mndis Moderate
191652a/Asdis, Na+ Low
203443a/Asdis Moderate
21
222472b/NO3- Low
233163a/Asdis Moderate
244583b/Mndis Moderate
Table 11. Comparative display of the two periods of the unconfined aquifer according to PoS index classes and identified dominant factors (Qfs).
Table 11. Comparative display of the two periods of the unconfined aquifer according to PoS index classes and identified dominant factors (Qfs).
Wet Period (May 2019 and 2020)
PoS Index Class123456
No. of samples1251310
% (percentage)4.559.0922.7359.094.550.00
Dominant factors (Qfs)NO3-NH4+MndisMndisMndis
AsdisAsdisAsdisAsdis
NH4+
Dry Period (October 2019 and 2020)
PoS Index class123456
No. of samples269500
% (percentage)9.0927.2740.9122.730.000.00
Dominant factors (Qfs)NO3-AsdisAsdisAsdis
AsdisMndisMndisMndis
NH4+
Table 12. Comparative display of the two periods of the confined aquifer according to PoS index classes and identified dominant factors (Qf).
Table 12. Comparative display of the two periods of the confined aquifer according to PoS index classes and identified dominant factors (Qf).
Wet Period (May 2019 and 2020)
PoS Index Class123456
No. of samples169322
% (percentage)4.1725.0037.5012.508.338.33
Dominant factors (Qfs)AsdisAsdisAsdisAsdisAsdisAsdis
MndisMndisMndisMndis
NH4+
NO3-
Dry Period (October 2019 and 2020)
PoS Index class123456
No. of samples1610420
% (percentage)4.3526.0943.4817.398.700.00
Dominant factors (Qfs)AsdisAsdisAsdisAsdisAsdis
MndisMndisMndisMndis
Na+
NO3-
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kampas, G.; Panagopoulos, A.; Gkiougkis, I.; Pouliaris, C.; Pliakas, F.-K.; Kinigopoulou, V.; Diamantis, I. Index-Based Groundwater Quality Assessment of Nestos River Deltaic Aquifer System, Northeastern Greece. Water 2024, 16, 352. https://doi.org/10.3390/w16020352

AMA Style

Kampas G, Panagopoulos A, Gkiougkis I, Pouliaris C, Pliakas F-K, Kinigopoulou V, Diamantis I. Index-Based Groundwater Quality Assessment of Nestos River Deltaic Aquifer System, Northeastern Greece. Water. 2024; 16(2):352. https://doi.org/10.3390/w16020352

Chicago/Turabian Style

Kampas, George, Andreas Panagopoulos, Ioannis Gkiougkis, Christos Pouliaris, Fotios-Konstantinos Pliakas, Vasiliki Kinigopoulou, and Ioannis Diamantis. 2024. "Index-Based Groundwater Quality Assessment of Nestos River Deltaic Aquifer System, Northeastern Greece" Water 16, no. 2: 352. https://doi.org/10.3390/w16020352

APA Style

Kampas, G., Panagopoulos, A., Gkiougkis, I., Pouliaris, C., Pliakas, F. -K., Kinigopoulou, V., & Diamantis, I. (2024). Index-Based Groundwater Quality Assessment of Nestos River Deltaic Aquifer System, Northeastern Greece. Water, 16(2), 352. https://doi.org/10.3390/w16020352

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