Next Article in Journal
Numerical Modelling of a Dam-Regulated River Network for Balancing Water Supply and Ecological Flow Downstream
Next Article in Special Issue
Physics-Informed Neural Networks-Based Salinity Modeling in the Sacramento–San Joaquin Delta of California
Previous Article in Journal
When Climate Change and Overexploitation Meet in Volcanic Lakes: The Lesson from Lake Bracciano, Rome’s Strategic Reservoir
Previous Article in Special Issue
Data-Driven Parameter Prediction of Water Pumping Station
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Principal Component Analysis and the Water Quality Index—A Powerful Tool for Surface Water Quality Assessment: A Case Study on Struma River Catchment, Bulgaria

1
“Georgi Benkovski” Air Force Academy, 1 St. St. Cyril and Methodius Str., 5855 Dolna Mitropolia, Bulgaria
2
Bulgarian Academy of Sciences, National Institute of Geophysics, Geodesy and Geography, Acad. G. Bonchev Str., bl. 3, 1113 Sofia, Bulgaria
3
Chair of Water Supply, Sewerage, Water and Wastewater Treatment, University of Architecture, Civil Engineering and Geodesy, 1 Hr. Smirnenski Blvd., 1046 Sofia, Bulgaria
4
Chair of Analytical Chemistry, Faculty of Chemistry and Pharmacy, Sofia University “St. Kliment Ohridski”, 1 J. Bourchier Blvd., 1164 Sofia, Bulgaria
*
Author to whom correspondence should be addressed.
Water 2023, 15(10), 1961; https://doi.org/10.3390/w15101961
Submission received: 26 March 2023 / Revised: 13 May 2023 / Accepted: 16 May 2023 / Published: 22 May 2023

Abstract

:
The water quality assessment of the surface water bodies (SWBs) is one of the major tasks of environmental authorities dealing with water management. The present study proposes a water quality assessment scheme for the investigation of the surface waters’ physicochemical status changes and the identification of significant anthropogenic pressures. It is designed to extract valuable knowledge from the Water Frame Directive (WFD) mandatory monitoring datasets. The water quality assessment scheme is based on the Canadian Council of Ministers of the Environment water quality index (CCME-WQI), trend analysis of estimated WQI values, and Principal Component Analysis (PCA) using calculated excursions during the determination of WQI values. The combination of the abovementioned techniques preserves their benefits and additionally provides important information for water management by revealing the latent factors controlling water quality, taking into account the type of the SWB. The results enable the identification of the anthropogenic impact on SWBs and the type of the corresponding anthropogenic pressure, prioritization and monitoring restoration measures, and optimization of conducted monitoring programs to reflect significant anthropogenic pressures. The proposed simple and reliable assessment scheme is flexible to introducing additional water quality indicators (hydrological, biological, specific pollutants, etc.), which could lead to a more comprehensive surface water quality assessment.

1. Introduction

Everyday human life necessitates the use of one of the Earth’s most vital resources—water. This use often leads to adverse impacts on the flora and fauna of the planet since about 80% of the used water is returned to the environment untreated, leading to decreased water quality. As a result, a worldwide increase in freshwater scarcity is observed. To reflect this problem, an adaptation of the relevant legislation has been initiated at the start of the new millennium—WHO Guidelines [1], Millennium Development Goals [2] and Directives [3] with the sole goal of improving the management and utilization of natural resources on the one hand, and to reduce the human pollution, on the other.
In addition to the anthropogenic influence, deterioration of the physical, chemical, or biological water quality parameters in surface waters is also due to natural processes such as the wearing of rocks, run-off, and ion exchange, which often makes the water resource unsuitable for its intended use. Therefore, to understand the changes in water quality, it is mandatory to monitor the possible increasing or decreasing trends in water quality to successfully perform water management. To achieve this goal, the responsible stakeholders should execute appropriate monitoring programs, including all required parameters to figure out the natural and anthropogenic processes controlling water quality [4]. The monitoring conducted by the environmental authorities produces large and complex data sets including water quality parameters with different sampling frequencies. Furthermore, it involves the determination of many water quality indicators at sampling locations positioned in different types of surface water bodies (SWBs) during different seasons [5,6].
The water quality index (WQI) approach is a preferred tool for water resource management as it estimates surface water quality in a single value. This value estimates the water quality at a specific place and at a specific time, taking into account water quality indicators’ adopted threshold values [7,8]. This way of presentation is convenient for the assessment of temporal surface water quality changes by combining WQI with the Mann–Kendall test [9,10,11] or with other techniques such as Spearman correlation tests [12] and rescaled range analysis [13]. The main disadvantage of such an approach is that the calculation of the integral WQI is performed by dimensionality reduction in the monitoring data set without taking into account its structure.
The WQI concept was introduced long ago in Germany [14], and since then, numerous WQIs have been developed [15,16], such as the National Sanitation Foundation WQI (NSFWQI) [17] and Oregon WQI (OWQI) [18]. They are applied in many regions of the world [19,20]. Different modifications of WQI have been introduced by the researchers, and comparisons of the different modifications exist in the scientific literature [21,22]. The Canadian Council of Ministers of the Environment (CCME) WQI is currently one of the most widely used, not only in all the provinces of Canada but in many countries worldwide [23,24,25,26], as it is one of the main recommended indices for assessing the quality of drinking water by the United Nations Environment Programme Global Environment Monitoring System (GEMS)/Water Programme [27]. Although there are many WQIs available worldwide, none is universally accepted due to the difference in climate and land use. The main advantages of CCME-WQI facilitating its application are its simplicity, flexibility to input water quality indicators and thresholds, and tolerance to missing data [7]. Thus, the adoption of CCME-WQI to different national environmental legislations could be easily performed.
To reveal the relationships within the data set and extract knowledge concerning surface water quality, a variety of multivariate statistical methods are added to the WQI. The Cluster Analysis (CA) is used to detect the groups of similarities between sampling locations and/or water quality parameters [12]. The Principal Component Analysis (PCA) is used to reveal the “hidden” factors controlling surface water quality [28,29,30,31,32]. Numerous groups of studies use both CA and PCA as additions to the WQI to introduce the specific outputs of both techniques to water quality assessment [11,33,34,35,36,37,38,39]. To investigate the temporal behavior of water quality indicators at regions with pronounced seasonality, Discriminant Analysis (DA) is used [13,40,41,42,43,44]. In some studies, multivariate methods for the selection of important water quality indicators, such as redundancy discriminant analysis [45] and the Boruta algorithm [10], are used. Where the monitoring scheme includes apportioning studies for enough sampling locations and specific for the investigated region water quality indicators [40,44], water quality risk assessment by prediction models is carried out [10,46,47]. In the case of high-density monitoring networks, GIS techniques are applied for the visualization of spatial patterns [48,49].
The above-mentioned multivariate statistical methods reveal the relationships between sampling points and/or water quality indicators in the original monitoring data but do not take into consideration the specificity of the investigated catchments and environmental legislation. None of these multivariate approaches added to the WQI in the proposed schemes for water quality assessment takes into account the adopted environmental thresholds, which could bias the estimation for catchments including different surface water bodies (SWBs) with different adopted thresholds. Additionally, the multivariate analysis of original monitoring data could not give a real insight into the “structure” of composite WQI, which could reflect different water quality scenarios [6].
The aim of this study is to develop a simple and reliable scheme for water quality assessment of a catchment that includes different types of SWBs. The case study of the Struma River catchment in Bulgaria is explored. The assessment is performed by the following steps: (i) calculation of WQI (using CCME-WQI) for each sampling location and estimation of temporal changes in water quality by the Mann–Kendall test; (ii) application of PCA using calculated excursions during WQI calculation to reveal latent factors controlling water quality, and (iii) spatial analysis based on factor score plots obtained by the PCA.

2. Materials and Methods

2.1. Struma River Catchment

The Struma River is one of the longest rivers in Bulgaria and is second only to the Maritsa River in terms of catchment area (10,797 km2). The catchment basin has a predominantly mountainous and hilly terrain. The river springs at a height of 2180 m at the southern slope of Vitosha and leaves Bulgarian territory at the village of Kulata (62 m above sea level). The river catchment has an oblong shape and asymmetry. It is located in Southwestern Bulgaria and includes parts of Kraishteto, Vitosha, Osogovo-Belasishka mountain range, Rila and Pirin, as well as many valleys closed between them (Figure 1).
The Struma River basin is divided into three parts [50]. The upper course covers the western and southwestern slopes of Vitosha, the Pernik, Radomir, and Kyustendil valleys. The middle course of the river starts from the southern periphery of the Kyustendil Valley to the southern border of the country. The lower course of the river is entirely on the territory of Greece.
The Struma River receives most of its right tributaries in Kraishteto (Konska, Svetla, Treklyanska, and Dragovishtitsa), and larger right tributaries in the middle course are Eleshnitsa, Logodashka, Lebnitsa, and Strumeshnitsa. The left tributaries of the river are significantly more abundant because they are formed in the high Rila and Pirin mountains.
The discharge of the river varies significantly along the course. At Pernik, the amount of water is around 2.2 m3/s, while at the border, it reaches 76.1 m3/s. The high water is observed in the period March–May, and the low water in July–October.
In 2021, about 412,000 people lived in the river catchment area, and compared to 2011, the population decreased by about 50,000. The urban population predominated (73 to 27%). The main part lives in the large regional cities—Pernik, Blagoevgrad, and Kyustendil, as well as in several municipal centers—Dupnitsa, Sandanski, Petrich, and Radomir.
Blagoevgrad district has the smallest share of the population connected to wastewater treatment plants (WWTPs)—32.7%. The Kyustendil region has 62.1%, and the Pernik region—73.8% [51]. There are 5 WWTPs downstream, which discharge directly into the Struma River or in its tributaries—WWTP Batanovtsi, WWTP Radomir, WWTP Dupnitsa, WWTP Kyustendil, and WWTP Blagoevgrad. The agricultural land, including arable land and pastures, covers 29.2% of the area of the Struma catchment.

2.2. Sampling Locations

The codes and descriptions of the 30 monitoring locations are presented in Table 1. The letter in the code represents the river course (A—upper, B—middle). The sampling stations located on the Struma River are numbered from the beginning of the course (A1, B1, etc.), while for tributaries, an additional number is added to specify the number of inlets between two neighboring locations on the main river. (A 2-1, B 1-1, etc.). On the transborder river Strumeshnitsa, two sampling locations (B 6-2U and B 6-2U) are sampled. Additionally, the type of SWB where the corresponding sampling point is located is given.

2.3. Data Acquisition

In this study, data obtained from the mandatory monitoring program of the Ministry of Environment and Waters for the Struma River catchment, conducted at the sampling locations described in Table 1 between 2010 and 2021, is used. The mandatory monitoring data for surface water quality include the following physicochemical parameters: dissolved oxygen (DissO2), pH, electrical conductivity (EC), ammonia (NH4+), nitrates (NO3), nitrites (NO2), total nitrogen (N), phosphates (PO43−), total phosphorus (P), and biochemical oxygen demand after 5 days (BOD5). The selected physicochemical parameters are analyzed at least 4 times a year [52] and are an important integral part of the ecological status assessment of SWBs. All analyses are performed in accredited laboratories every month according to national legislation requirements [52].

2.4. Water Quality Index

The water quality integral indicator applied in this study is the composite WQI developed by the CCME [53]. Since 2006 it has been recommended by the United Nations Environment Programme (UNEP) for application in assessing surface water quality [54]. The CCME-WQI is based on three factors characterizing the anthropogenic impact on the water quality:
  • F1 (Scope) represents the percentage of water quality indicators not meeting the regulatory guideline values (“failed variables”) over the total number of variables included in the water quality assessment;
F 1 = N u m b e r   o f   f a i l e d   v a r i a b l e s T o t a l   n u m b e r   o f   v a r i a b l e s × 100
  • F2 (Frequency) represents the percentage of measurements in which a water quality indicator exceeds the guideline values (“failed tests”) over the total number of tests (measurements);
F 2 = N u m b e r   o f   f a i l e d   t e s t s T o t a l   n u m b e r   o f   t e s t s × 100
  • F3 (Amplitude) represents the extent of deviation of the “failed tests” values relative to the corresponding guideline values. The amplitude is calculated utilizing a three-step algorithm, at the beginning of which an assessment is made of the magnitude of the deviations (excursion) of the so-called “bad samples” relative to the corresponding maximum allowable concentrations:
e x c u r s i o n i = F a i l e d   T e s t   V a l u e i O b j e c t i v e j 1
or
e x c u r s i o n i = O b j e c t i v e j F a i l e d   T e s t   V a l u e i 1 ,
where Failed Test Valuei is the value of the “bad” sample, and Objectivej is the reference value of the maximum allowable concentration for the corresponding quality indicator. Equation (4) is applied to calculate the amplitude of indicators whose reference values of maximum allowable concentration decrease when moving to a worse category (e.g., dissolved oxygen content). The normalized sum of deviations NSE is then calculated:
N S E = i = 1 n e x c u r s i o n i N u m b e r   o f   t e s t s .
Finally, the amplitude (F3) is calculated using the formula:
F 3 = N S E 0.01 N S E + 0.01 .
The calculation of WQI is performed by aggregation of the obtained factors as follows:
W Q I = 100 F 1 2 + F 2 2 + F 3 2 1.732 .
The equation above gives a value between 0 and 100, with 0 being the “worst” and 100 indicating the “best” water quality. Within this range, the original CCME-WQI scale [53] ranks water quality into 5 categories: poor, marginal, fair, good, and excellent (Table 2). According to the Bulgarian legislation [52], harmonized with the European legislation through the WFD, the characterization of the river waters status based on physicochemical parameters is performed in 3 categories—moderate, good, and very good (Table 2) [6]. Currently, there is no national WQI in Bulgaria. Testing of various WQIs began in 2007, and the CCME-WQI was found to be the most flexible to be applied according to the requirements of national legislation [55,56].
The classical CCME-WQI approach [53] was harmonized with the adoption of the Bulgarian legislation reference ranges [52] using threshold values for “good” status in the subsequent calculations (see Table 3). For most of the investigated physicochemical indicators (except for pH and DissO2), the threshold value is the lower limit of the “good’ status range presented in Table 3. “Good” status based on the respective parameter is achieved when the measured value is lower than the threshold value (for example, less than 750 µS/cm for EC). When the measured value of the indicator is higher than this threshold, the status is classified as “moderate”. The excursions of these water quality parameters are calculated using Equation (3). For DissO2, the opposite approach was adopted; a “good” status is achieved when being higher than the threshold value for the respective range (6 mg/L), and Equation (4) is used for calculations of the excursions. The “good” status is only achieved for pH when the measured value falls within the reference range (Table 3). The proposed harmonization was already presented and validated for another Bulgarian river catchment [6]. Following the aforementioned scheme, annual WQIs are calculated for each sampling point during the monitoring period.

2.5. Statistical Data Treatment

2.5.1. Mann–Kendall Test

The Mann–Kendall test is a common and practical nonparametric method to determine the monotonic trend in different environmental studies producing time series [57]. Its main advantage is that there is no requirement for a specific distribution of the data, and missing data is allowed [58]. In this study, the variance correction version of the Mann–Kendall test is performed [59]. The ‘modifiedmk’ R package (version 1.6) is used for all trend analysis calculations.

2.5.2. Principal Component Analysis (PCA)

A well-documented statistical approach in surface water studies, the Principal Component Analysis (PCA) [60,61,62], is used for the analysis of the obtained monitoring data. The data interpretation is achieved through the latent factors (principal components) formed by original variables (the water quality indicators) by looking at the major sources of variance in the data [63]. The latent factors form a matrix composed of factor loadings (weights of the original variables) and factor scores (projections of the sampling locations on the principal component axis). Thus, data can be interpreted using a fewer number of principal components than the number of original variables extracting useful information about hidden relationships between original variables and revealing the impact of latent factors on different sampling locations. All calculations for PCA were performed using PLS Toolbox 9.0 (Eigenvector Research, Manson, WA, USA) of MATLAB R2021a (The Mathworks, Natick, MA, USA).

3. Results

3.1. Statistical Analysis

The basic statistics for the 10 water quality indicators is presented in Table 4. The comparison of the results obtained in this study (2010–2021) with previously published data for the period 1989–1998 outlines a decrease in measured water quality indicator values except for pH and DissO2 [64]. The comparison between monitoring results and guideline regulation values reveals three groups of indicators responsible for the moderate status of the surface waters in the Struma River catchment. The first group consists of water quality indicators with more than 50% of exceedings: P (69.4%), PO43− (68.4%), and N (51.2%). They are followed by NO2 (38.1%), BOD5 (36.6%), NO3 (25.8%), and NH4+ (23.7%). The pH (12.5%), DissO2 (12.1%), and EC (11.1%) have the smallest contribution to the moderate water status.

3.2. WQI

The results of the calculated water quality indices for all 30 monitoring stations along the Struma River catchment are presented in Figure 2. They demonstrate that in more than half of the points, the surface water status is “moderate”, in about 30%—“good”, and in only 16%—“very good”. In most parts of the catchment, the water quality does not meet the goals set in the Struma River Basin Management Plan to improve the status to “good”. The proportion of “moderate” status in the upper course (A) of the river is particularly high (almost ¾).
The annual trend of WQI for the investigated sampling locations is presented in Figure 3. The performed Mann–Kendall test outlines 15 out of 30 sampling locations with a significant rate of change in the WQI, while only B 1-1 (the Eleshnitsa river) has a decreasing trend during the investigated period.

3.3. PCA

The relationships between water quality indicators and/or sampling locations were examined by PCA. The input data set used for PCA consists of variables—water quality indicators and objects—and sampling locations. The water quality indicators are presented by the excursions, calculated according to Equations (3) and (4). Only for pH was the original data used. Prior to the analysis, the input data underwent autoscaling pretreatment. The PCA reveals that the first three principal components explain almost 70% of dataset variation. The number of latent factors is determined based on their eigenvalues and the internal model validation error. In the formation of the first principal component, explaining 42.68% of data variance, the nutrients (N, P), PO43−, NH4+, and BOD5 have a significant impact (Figure 4). This latent factor could be conditionally named “urban waste waters”, but to a smaller extent, could also be indicative of agricultural activity, which is not recognized as driving pressure for the Struma River catchment. The second principal component (16.11%) is formed predominantly by EC and NO3 and reflects mainly the industrial wastewaters in the Struma River catchment. The last principal component (11.03%) is positively correlated with pH and negatively correlated with DissO2. This latent factor could be related to the surface water pollution response at the investigated sampling locations. A high pH may inhibit the photosynthesis of algae, which will reduce the DissO2 [65].
The factor score plots resemble the impact of the abovementioned latent factors at the sampling locations in the Struma River catchment. The PC1 score plot outlines the more pronounced impact of urban wastewater on sampling locations in the upper course (A) located on the Struma River—A3 (the Struma River at Batanovtsi town), A4 (the Struma River before the Pchelina Dam) and A5 (the Pchelina Dam), and on tributaries from the middle course (B) B 0-3 (the Banshtitsa River), B 2-1 (the Blagoevgradska Bistrica River), B 5-1 (the Sandanska Bistrica River), and B 6-2D (the Strumeshnitsa River before flowing in the Struma River) (Figure 5).
The influence of urban wastewater is presented in Figure 6. The most influenced are sampling stations located on the upper course of the Struma River—A3, A4, and A5 followed by the tributary Razmetanitsa River in the middle course (B1-4).
The impact of surface water pollution response does not discriminate sampling locations, but locations influenced by urban wastewaters: A5 (the Pchelina Dam) and B1-4 (the Razmetanitsa River), which possess more scattered factor scores (Figure 7).

4. Discussion

The Struma River catchment, one of the longest in Bulgaria, has a significant impact on water quality, and as such, it has been explored in this study. Our results show that there is a significant difference between surface water quality in the upper and the middle course of the catchment due to the specifically located urban and industrial pressures.
The upstream sampling location of the Studena dam (A1) in the upper course is used as a drinking water supply source for the town of Pernik and has had a “good” or “very good” status for the entire period of investigation. The PCA analysis does not reveal any anthropogenic pressures, and the trend analysis results show a small increasing trend in WQI of 0.22 units per year. The following sampling locations situated along the Struma River (A2–A6) have moderate physicochemical status during the investigated period. The A2 (Struma River at Pernik) and A3 (Struma River at Batanovtsi town) are impacted by the urban and industrial wastewaters from the Pernik region, while the wastewater from the Radomir region affected the downstream sampling stations A4 (the Struma River before the Pchelina Dam), A5 (the Pchelina Dam), and A6 (the Struma River at Zabliano village). The Pchelina Dam was built in 1975 as a precipitator of the polluted waters from the communal and industrial pollution of the Struma River waters from the Pernik and Radomir regions. The complex nature of the anthropogenic pressure is confirmed by the previous studies on surface waters [66] and sediments of the Pchelina Dam [67]. The high variation of PC3 factor scores for A5 is potentially caused by the changes in the surface water quality and operation of the small hydroelectric station located near the Pchelina Dam. The WQIs of the sampling locations downstream of the Struma River are improved after passing through the reservoir—a conformation of the study in the period between 1989 and 1998 [68]. The sampling location of the Struma River at Razdavitsa village (A7) has had a good status for 6 out of 12 investigated years. The significant water quality improvement trend (0.43 WQI units per year) for A3 could be related to the improvement of the WWTP of Pernik and the related sewage network. The tributary Konska River (A2-1) has a pattern similar to the sampling location A2 (Struma near Pernik town), and its “moderate” status is due to urban wastewaters and small industries from Breznik town and surrounding villages. The tributary Arkata River (A3-1) has dominantly “good” and “very” good status with a significant positive water quality trend. The “moderate” status during 2017 and 2018 is caused by the small industries located in the tributary catchment (see Figure 7).
Based on the obtained results for the upper course of the Struma River catchment, the main mitigation measures for improving surface water quality should be directed to the more efficient treatment of industrial wastewaters from the regions of Pernik and Radomir, accompanied by the improvement of their sewage networks.
The sampling locations of the Struma River in the middle course (B1–B7) have predominantly “good” physicochemical status with increasing downstream water quality tendency, similar to the tendency of sampling stations in the upper course, located after the Pchelina Dam. The PCA results do not reveal any significant urban or industrial impact on the abovementioned sampling locations. There is a significant improvement in the water quality trend for the Struma River at Nevestino village (B1), while the most pronounced within the entire catchment of the Struma River is at Kresna town (B5) with 0.61 WQI units per year. The first two tributaries in the middle course—the Dragovishtitsa River (B0-1) and the Sovolyanska Bistrica River (B0-2)—have significant positive water quality trends. It should be mentioned that the Dragovishtitsa River is a transboundary river, and the WQI categories of B0-1, located close to the border, reflect the anthropogenic pressures on the territory of the Republic of Serbia. The Banshtitsa River (B0-3) has a “moderate” status during the whole investigated period, which is caused by the urban wastewater of Kyustendil town. Extreme exceedings (over 25 times) for PO43− and P are observed. This is in agreement with the results previously reported [69], where the water quality of the Banshtitsa River is classified as a WWTP effluent. The tributaries that flow into the Struma River between sampling locations B1 and B2 could be divided into three groups. The Eleshnitsa River (B1-1) and the Dyakovo Dam (B1-2) have “good” and “very good” statuses without evidence of any anthropogenic impact from the PCA analysis. The Dupnishka Bistrica River (B1-3) and the Razmetanitsa River (B1-4) have varying “moderate” to “good” water quality statuses during the period of investigation. The reasons for the estimated “moderate” status are urban wastewater of Dupnitsa town and industrial wastewater of the thermal power plant Bobov dol for B1-3 and B1-4, respectively. The variation in the surface water quality of the Razmetanitsa River is supported by the scattered surface water pollution response factor scores (PC3). The German River (B1-5) has a “moderate” status during the whole period of observation. The PCA results do not detect upstream urban and industrial pressures at sampling points B1-3 and B1-4. It seems that surface water processes that take place between sampling locations B1-4 and B1-5 “change” the relationships between the physicochemical parameters from that included in PC1 and PC2. This assumption is supported by the registered extreme exceedings (over 25 times) for PO43− and P in the German River and from the very high measured concentration of total organic carbon during 2021–9.48 mg/L. The introduction of additional water quality indicators specific to upstream pressures could overcome their “masking”. The next two downstream tributaries, influenced by urban wastewater, are the Blagoevgradska Bistrica River (B2-1) and the Sandanska Bistritsa River (B5-1). Both sampling locations have a significant increase in the WQI. This improvement is due to the putting into the operation of WWTP Blagoevgrad in 2012 and the improvement of the sewage network of Sandanski town for B2-1 and B5-1, respectively. The last two tributaries sampling locations where PCA results indicate urban wastewater impacts are B6-2U and B6-2D, located on the transboundary Strumeshnitsa River. The “moderate” status of the surface waters coming from the Republic of North Macedonia (B6-2U), which is characterized by PO43− and P exceedings, does not undergo significant change at B6-2D.
In contrast with the upper course, the main anthropogenic pressures in the middle course of the Struma River have urban origins. Since the introduced mitigation measures for the treatment and collection of urban waters from Blagoevgrad and Sandanski town led to the improvement of surface water quality, the reconstruction of the sewage network of Kyustendil town is necessary to be planned in further river catchment management.
In general, along the course of the main river, the anthropogenically impacted tributaries do not lead to the deterioration of the water quality of the Struma River. The PC1 and PC2 factor scores have similar values for the first Studena Dam (A1) and the last Struma River at Topolnitsa village (B7) sampling points in the Struma River catchment.

5. Conclusions

The present study introduces an applicable and reliable scheme for water quality assessment of catchments, including different types of SWBs using WFD mandatory monitoring datasets. The assessment is performed by combining CCME-WQI, the Mann–Kendall test, and PCA, exploring in this way their advantages and complementarity.
The proposed approach provides important knowledge to the environmental authorities for (i) identification of the anthropogenic impact on SWBs and the type of the corresponding anthropogenic pressure, (ii) prioritization and monitoring of mitigation measures, and (iii) optimization of conducted monitoring programs to reflect significant anthropogenic pressures. The obtained knowledge could be used for future river catchment management. Additionally, the proposed water quality assessment scheme could be applied to monitor the effectiveness of the implemented mitigation measures.
The assessment scheme is flexible to introducing additional water quality indicators (priority pollutants, hydrological, and biological), which could lead to a more comprehensive surface water quality assessment. The opportunity for introducing special pollutants for a particular river basin opens the door for the proposed approach to be implemented not only in another river catchment but to be used on the river basin management level.

Author Contributions

Conceptualization, M.V. and S.T.; methodology, I.B., M.V., T.V. and S.T.; software, M.V. and S.T.; validation, I.B., M.V., T.V. and S.T.; resources, M.V. and T.V.; data curation, I.B. and M.V.; writing—review and editing, I.B., M.V., T.V. and S.T.; visualization, S.T.; funding acquisition, M.V., T.V. and S.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science Program “Environmental Protection and Reduction of Risks of Adverse Events and Natural Disasters”, approved by the Resolution of the Council of Ministers No. 577/17.08.2018 and supported by the Ministry of Education and Science (MES) of Bulgaria (Agreement No. Д01-271/09.12.2022)” and the APC was partially funded by T.V. and S.T.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

This work has been carried out in the framework of the National Science Program “Environmental Protection and Reduction of Risks of Adverse Events and Natural Disasters”, approved by the Resolution of the Council of Ministers No. 577/17.08.2018 and supported by the Ministry of Education and Science (MES) of Bulgaria (Agreement No. Д01-271/09.12.2022).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Victor, R.; Kotter, R.; O′Brien, G.; Mitropoulos, M.; Panayi, G. WHO Guidelines for the safe use of wastewater, excreta and greywater—Volume 1: Policy and regulatory aspects. Int. J. Environ. Stud. 2006, 65, 157–176. [Google Scholar] [CrossRef]
  2. The United Nation. The Millennium Development Goals Report; The United Nation: New York, NY, USA, 2015. [Google Scholar]
  3. European Parliament, Council of the European Union. Directive 2000/60/EC of the European Parliament and of the Council establishing a framework for Community action in the field of water policy. Off. J. Eur. Union 2000, 327, 1–73. [Google Scholar]
  4. European Commission, Directorate-General for Environment, Publications Office. Monitoring under the Water Framework Directive. Guidance Document No 7; European Commission, Directorate-General for Environment, Publications Office: Luxembourg, 2012. [Google Scholar]
  5. Gurjar, S.; Tare, V. Spatial-temporal assessment of water quality and assimilative capacity of river Ramganga, a tributary of Ganga using multivariate analysis and QUEL2K. J. Clean. Prod. 2019, 222, 550–564. [Google Scholar] [CrossRef]
  6. Yotova, G.; Varbanov, M.; Tcherkezova, E.; Tsakovski, S. Water quality assessment of a river catchment by the composite water quality index and self-organizing maps. Ecol. Indic. 2021, 120, 106872. [Google Scholar] [CrossRef]
  7. Tyagi, S.; Sharma, B.; Singh, P.; Dobhal, R. Water Quality Assessment in Terms of Water Quality Index. Am. J. Water Res. 2013, 1, 34–38. [Google Scholar] [CrossRef]
  8. Poonam, T.; Tanushree, B.; Sukalyan, C. Water Quality Indices—Important Tols for Water Quality Assessment: A Review. Int. J. Adv. Chem. 2015, 1, 15–29. [Google Scholar]
  9. Hashim, M.; Nayan, N.; Setyowati, D.; Said, Z.; Mahat, H.; Saleh, Y. Analysis of Water Quality Trends Using the Mann-Kendall Test and Sen’s Estimator of Slope in a Tropical River Basin. Pollution 2021, 7, 933–942. [Google Scholar] [CrossRef]
  10. Wang, X.; Liu, X.; Wang, L.; Yang, J.; Wan, X.; Liang, T. A holistic assessment of spatiotemporal variation, driving factors, and risks influencing river water quality in the northeastern Qinghai-Tibet Plateau. Sci. Total Environ. 2022, 851, 157942. [Google Scholar] [CrossRef]
  11. Chong, L.; Li, B.; Sun, Z.; Hu, C.; Meng, X.; Gao, J. Temporal and spatial variation in water quality in the Yangtze Estuary from 2012 to 2018. Environ. Sci. Pollut. Res. 2022, 29, 76235–76250. [Google Scholar] [CrossRef]
  12. Fraga, M.D.S.; da Silva, D.D.; Reis, G.B.; Guedes, H.A.S.; Elesbon, A.A.A. Temporal and spatial trend analysis of surface water quality in the Doce River basin, Minas Gerais, Brazil. Environ. Dev. Sustain. 2021, 23, 12124–12150. [Google Scholar] [CrossRef]
  13. Chong, L.; Zhong, J.; Sun, Z.; Hu, C. Temporal variations and trends prediction of water quality during 2010–2019 in the middle Yangtze River, China. Environ. Sci. Pollut. Res. 2023, 30, 28745–28758. [Google Scholar] [CrossRef] [PubMed]
  14. Horton, R.K. An index number system for rating water quality. J. Water Pollut. Control Fed. 1965, 37, 300–306. [Google Scholar]
  15. Bharti, N.; Katyal, D. Water Quality Indices Used for Surface Water Vulnerability Assessment. Int. J. Environ. Sci. 2011, 2, 154–173. [Google Scholar]
  16. Abbasi, T.; Abbasi, S. Water-Quality Indices; Elsevier: Amsterdam, The Netherlands, 2012; pp. 353–356. [Google Scholar]
  17. Brown, R.; McClelland, N.; Deininger, R.; Tozer, R.G. A water quality index—Do we dare? Water Sew. Works 1970, 117, 339–343. [Google Scholar]
  18. Dunnette, D.A. A geographically variable water quality index used in Oregon. J. Water Pollut. Control Fed. 1979, 51, 53–61. [Google Scholar]
  19. Dezfooli, D.; Hosseini-Moghari, S.-M.; Ebrahimi, K.; Araghinejad, S. Classification of water quality status based on minimum quality parameters: Application of machine learning techniques. Model. Earth Syst. Environ. 2018, 4, 311–324. [Google Scholar] [CrossRef]
  20. Noori, R.; Berndtsson, R.; Hosseinzadeh, M.; Adamowski, J.F.; Abyaneh, M.R. A critical review on the application of the National Sanitation Foundation Water Quality Index. Environ. Pollut. 2019, 244, 575–587. [Google Scholar] [CrossRef]
  21. Wong, Y.J.; Shimizu, Y.; He, K.; Sulaiman, N.M.N. Comparison among different ASEAN water Qual-Ity indices for the assessment of the spatial variation of surface water quality in the Selangor river basin, Malaysia. Environ. Monit. Assess. 2020, 192, 644. [Google Scholar] [CrossRef]
  22. Gikas, G.D.; Lergios, D.; Tsihrintzis, V.A. Comparative Assessment of the Application of Four Water Quality Indices (WQIs) in Three Ephemeral Rivers in Greece. Water 2023, 15, 1443. [Google Scholar] [CrossRef]
  23. Panagopoulos, Y.; Alexakis, D.E.; Skoulikidis, N.T.; Laschou, S.; Papadopoulos, A.; Dimitriou, E. Implementing the CCME Water Quality Index for the Evaluation of the Physicochemical Quality of Greek Rivers. Water 2022, 14, 2738. [Google Scholar] [CrossRef]
  24. Hu, L.; Chen, L.; Li, Q.; Zou, K.; Li, J.; Ye, H. Water quality analysis using the CCME-WQI method with time series analysis in a water supply reservoir. Water Suppl. 2022, 22, 6281–6295. [Google Scholar] [CrossRef]
  25. Dao, V.; Urban, W.; Hazra, S.B. Introducing the modification of Canadian water quality index. Groundw. Sustain. Dev. 2020, 11, 100457. [Google Scholar] [CrossRef]
  26. Bilgin, A. Evaluation of surface water quality by using Canadian Council of Ministers of the Environment Water Quality Index (CCME WQI) method and discriminant analysis method: A case study Coruh River Basin. Environ. Monit. Assess. 2018, 190, 554. [Google Scholar] [CrossRef]
  27. UNEP—United Nations Environment Programme. Global Drinking Water Quality Index Development and Sensitivity Analysis Report; UNEP: Toronto, ON, Canada, 2007; p. 60. [Google Scholar]
  28. Kumar, R.; Dutt, V.; Raina, A.; Sharma, N. Spatial water quality assessment of a mountain stream in northwestern India using multivariate statistical techniques. Environ. Monit. Assess. 2022, 194, 785. [Google Scholar] [CrossRef] [PubMed]
  29. Subramanian, A.; Baskar, S. Water quality assessment of Noyyal river using water quality index (WQI) and multivariate statistical techniques. Water Sci. 2022, 36, 85–98. [Google Scholar] [CrossRef]
  30. Ustaoğlu, F.; Taş, B.; Tepe, Y.; Topaldemir, H. Comprehensive assessment of water quality and associated health risk by using physicochemical quality indices and multivariate analysis in Terme River, Turkey. Environ. Sci. Pollut. Res. 2021, 28, 62736–62754. [Google Scholar] [CrossRef] [PubMed]
  31. Franco, H.; Custodio, M.; Peñaloza, R.; De La Cruz, H. Application of Multivariate Statistical Methods and Water Quality Index for the Evaluation of Surface Water Quality in the Cunas River Basin, Peru. Asian J. Water Environ. Pollut. 2021, 18, 19–27. [Google Scholar] [CrossRef]
  32. Panjgotra, S.; Rishi, M.; Awasthi, A. Water Quality Assessment of River Tawi, Jammu Using Water Quality Index and Multivariate Statistical Techniques. Water Resour. 2022, 49, 1059–1069. [Google Scholar] [CrossRef]
  33. Roy, B.; Manna, A. The Status of Surface Water in West Tripura District, India: An Approach by Using Water Quality Index and Multivariate Statistical Technique. Asian J. Water Environ. Pollut. 2021, 18, 27–36. [Google Scholar] [CrossRef]
  34. Gupta, A.K.; Kumar, A.; Maurya, U.K.; Singh, D.; Islam, S.; Rathore, A.C.; Kumar, P.; Singh, R.; Madhu, M. Comprehensive spatio-temporal benchmarking of surface water quality of Hindon River, a tributary of river Yamuna, India: Adopting multivariate statistical approach. Environ. Sci. Pollut. Res. 2023, in press. [Google Scholar] [CrossRef]
  35. Aydin, H.; Ustaoğlu, F.; Tepe, Y.; Soylu, E.N. Assessment of water quality of streams in northeast Turkey by water quality index and multiple statistical methods. Environ. Forensics 2021, 22, 270–287. [Google Scholar] [CrossRef]
  36. Zhang, Z.-M.; Zhang, F.; Du, J.-L.; Chen, D.-C. Surface Water Quality Assessment and Contamination Source Identification Using Multivariate Statistical Techniques: A Case Study of the Nanxi River in the Taihu Watershed, China. Water 2022, 14, 778. [Google Scholar] [CrossRef]
  37. Das, S.; Sarkar, R. Monitoring and evaluating the spatiotemporal variations of the water quality of a stretch of the Bhagirathi-Hugli River, West Bengal, India, using geospatial technology and integrated statistical methods. Environ. Sci. Pollut. Res. 2021, 28, 15853–15869. [Google Scholar] [CrossRef] [PubMed]
  38. Alam, R.; Ahmed, Z.; Seefat, S.M.; Nahin, K.T.K. Assessment of surface water quality around a landfill using multivariate statistical method, Sylhet, Bangladesh. Environ. Nanotechnol. Monit. Manag. 2021, 15, 100422. [Google Scholar] [CrossRef]
  39. Dimri, D.; Daverey, A.; Kumar, A.; Sharma, A. Monitoring water quality of River Ganga using multivariate techniques and WQI (Water Quality Index) in Western Himalayan region of Uttarakhand, India. Environ. Nanotechnol. Monit. Manag. 2021, 15, 100375. [Google Scholar] [CrossRef]
  40. Varol, M.; Karakaya, G.; Alpaslan, K. Water quality assessment of the Karasu River (Turkey) using various indices, multivariate statistics and APCS-MLR model. Chemosphere 2022, 308, 136415. [Google Scholar] [CrossRef] [PubMed]
  41. Giao, N.; Dan, T.; Ni, D.; Anh, P.; Nhien, H. Spatiotemporal Variations in Physicochemical and Biological Properties of Surface Water Using Statistical Analyses in Vinh Long Province, Vietnam. Water 2022, 14, 2200. [Google Scholar] [CrossRef]
  42. Giao, N. Analysis of Surface Water Quality using Multivariate Statistical Approaches: A case study in Ca Mau Peninsula, Vietnam. Pollution 2022, 8, 463–477. [Google Scholar] [CrossRef]
  43. Ofosu, S.; Adjei, K.; Odai, S. Assessment of the quality of the Densu river using multicriterial analysis and water quality index. Appl. Water Sci. 2021, 11, 183. [Google Scholar] [CrossRef]
  44. Mamun, M.; An, K.-G. Application of Multivariate Statistical Techniques and Water Quality Index for the Assessment of Water Quality and Apportionment of Pollution Sources in the Yeongsan River, South Korea. Int. J. Environ. Res. Public Health 2021, 18, 8268. [Google Scholar] [CrossRef]
  45. Gruss, Ł.; Wiatkowski, M.; Pulikowski, K.; Kłos, A. Determination of Changes in the Quality of Surface Water in the River—Reservoir System. Sustainability 2021, 13, 3457. [Google Scholar] [CrossRef]
  46. Elsayed, S.; Hussein, H.; Moghanm, F.; Khedher, K.; Eid, E.; Gad, M. Application of Irrigation Water Quality Indices and Multivariate Statistical Techniques for Surface Water Quality Assessments in the Northern Nile Delta, Egypt. Water 2020, 12, 3300. [Google Scholar] [CrossRef]
  47. Wong, Y.J.; Shimizu, Y.; Kamiya, A.; Maneechot, L.; Bharambe, K.P.; Fong, C.S.; Sulaiman, N.M.N. Application of artificial intelligence methods for monsoonal river classification in Selangor river basin, Malaysia. Environ. Monit. Assess. 2021, 193, 438. [Google Scholar] [CrossRef] [PubMed]
  48. Azhari, H.E.; Cherif, E.K.; Sarti, O.; Azzirgue, E.M.; Dakak, H.; Yachou, H.; Esteves da Silva, J.C.G.; Salmoun, F. Assessment of Surface Water Quality Using the Water Quality Index (IWQ), Multivariate Statistical Analysis (MSA) and Geographic Information System (GIS) in Oued Laou Mediterranean Watershed, Morocco. Water 2023, 15, 130. [Google Scholar] [CrossRef]
  49. Thuan, N.C. Assessment of Surface Water Quality in the Hau Giang Province Using Geographical Information System and Statistical Approaches. J. Ecol. Eng. 2022, 23, 265–276. [Google Scholar] [CrossRef]
  50. FORCOM. Geography of Bulgaria: Physical and Socio-Economic Geography; FORCOM: Sofia, Bulgaria, 2002; p. 760. [Google Scholar]
  51. Regions, Districts and Municipalities in the Republic of Bulgaria. Electronic Edition, NSI. Available online: https://www.nsi.bg/sites/default/files/files/publications/ROO_2021.zip (accessed on 17 March 2023).
  52. Ordinance N-4/2012 for Characterisation of Surface Waters. D V 2013, 22, 9–46. Available online: http://eea.government.bg/bg/legislation/water/Naredba13.pdf (accessed on 17 March 2023). (In Bulgarian).
  53. Neary, B.; Cash, K.; Hébert, S.; Khan, H.; Saffran, K.; Swain, L.; Williamson, D.; Wright., R. Canadian Water Quality Guidelines for the Protection of Aquatic Life; CCME Water Quality Index 1.0 Technical Report; Canadian Council of Ministers of the Environment: Winnipeg, MB, Canada, 2001. [Google Scholar]
  54. Radeva, K.; Seymenov, K. Assessment of Physicochemical Properties and Water Quality of the Lom River (NW Bulgaria). In Smart Geography. Key Challenges in Geography; Springer: Berlin/Heidelberg, Germany, 2020. [Google Scholar] [CrossRef]
  55. Varbanov, M.; Gartsiyanova, K. Index assessment of the water quality—A case study of Bulgarian rivers. Geogr. Tour. 2017, 5, 41–49. [Google Scholar] [CrossRef]
  56. Rickwood, C.; Carr, G. Development and sensitivity analysis of a global drinking water quality index. Environ. Monit. Assess. 2009, 156, 73–90. [Google Scholar] [CrossRef]
  57. Kisi, O.; Ay, M. Comparison of Mann–Kendall and innovative trend method for water quality parameters of the Kizilirmak River, Turkey. J. Hydrol. 2014, 513, 362–375. [Google Scholar] [CrossRef]
  58. Güçlü, Y. Multiple Şen-innovative trend analyses and partial Mann-Kendall test. J. Hydrol. 2018, 566, 685–704. [Google Scholar] [CrossRef]
  59. Yue, S.; Wang, C. The Mann-Kendall Test Modified by Effective Sample Size to Detect Trend in Serially Correlated Hydrological Series. Water Resour. Manag. 2004, 18, 201–218. [Google Scholar] [CrossRef]
  60. Platikanov, S.; Rodriguez-Mozaz, S.; Huerta, B.; Barceló, D.; Cros, J.; Batle, M.; Poch, G.; Tauler, R. Chemometrics quality assessment of wastewater treatment plant effluents using physicochemical parameters and UV absorption measurements. J. Environ. Manag. 2014, 140, 33–44. [Google Scholar] [CrossRef]
  61. Navarro, A.; Tauler, R.; Lacorte, S.; Barceló, D. Occurrence and transport of pesticides and alkylphenols in water samples along the Ebro River Basin. J. Hydrol. 2010, 383, 18–29. [Google Scholar] [CrossRef]
  62. Acquavita, A.; Aleffi, I.F.; Benci, C.; Bettoso, N.; Crevatin, E.; Milani, L.; Tamberlich, F.; Toniatti, L.; Barbieri, P.; Licen, S.; et al. Annual characterization of the nutrients and trophic state in a Mediterranean coastal lagoon: The Marano and Grado Lagoon (northern Adriatic Sea). Reg. Stud. Mar. Sci. 2015, 2, 132–144. [Google Scholar] [CrossRef]
  63. Jolliffe, I. Principal Component Analysis, 2nd ed.; Springer: New York, NY, USA, 2002. [Google Scholar]
  64. Mihailov, G.; Simeonov, V.; Nikolov, N.; Mirinchev, G. Water Quality Environmetric Study of the Struma River Basin, Bulgaria, Part I: Water quality Long-Term Trends (1989–1998). Toxicol. Environ. Chem. 2002, 83, 1–12. [Google Scholar] [CrossRef]
  65. Zang, C.; Huang, S.; Wu, M.; Du, S.; Scholz, M.; Gao, F.; Lin, C.; Guo, Y.; Dong, Y. Comparison of Relationships Between pH, Dissolved Oxygen and Chlorophyll a for Aquaculture and Non-aquaculture Waters. Water Air Soil Pollut. 2011, 219, 157–174. [Google Scholar] [CrossRef]
  66. Gartsiyanova, K. Assessment of the water quality of the “Pchelina” Reservoir. Probl. Geogr. 2017, 1, 62–71. [Google Scholar]
  67. Venelinov, T.; Mihaylova, V.; Peycheva, R.; Todorov, M.; Yotova, G.; Todorov, B.; Lyubomirova, V.; Tsakovski, S. Sediment Assessment of the Pchelina Reservoir, Bulgaria. Molecules 2021, 26, 7517. [Google Scholar] [CrossRef]
  68. Astel, A.; Tsakovski, S.; Barbieri, P.; Simeonov, V. Comparison of self-organizing maps classification approach with cluster and principal components analysis for large environmental data sets. Water Res. 2007, 41, 4566–4578. [Google Scholar] [CrossRef] [PubMed]
  69. Yotova, G.; Lazarova, S.; Kudłak, B.; Zlateva, B.; Mihaylova, V.; Wieczerzak, M.; Venelinov, T.; Tsakovski, S. Assessment of the Bulgarian Wastewater Treatment Plants’ Impact on the Receiving Water Bodies. Molecules 2019, 24, 2274. [Google Scholar] [CrossRef]
Figure 1. Sampling map of the Struma River catchment, Bulgaria.
Figure 1. Sampling map of the Struma River catchment, Bulgaria.
Water 15 01961 g001
Figure 2. The annual WQI values for all sampling locations during the period of observation.
Figure 2. The annual WQI values for all sampling locations during the period of observation.
Water 15 01961 g002
Figure 3. Annual change in WQI for investigated sampling locations.
Figure 3. Annual change in WQI for investigated sampling locations.
Water 15 01961 g003
Figure 4. Factor loading plots for the first 3 latent factors.
Figure 4. Factor loading plots for the first 3 latent factors.
Water 15 01961 g004
Figure 5. Factor score plot for PC1.
Figure 5. Factor score plot for PC1.
Water 15 01961 g005
Figure 6. Factor score plot for PC2.
Figure 6. Factor score plot for PC2.
Water 15 01961 g006
Figure 7. Factor score plot for PC3.
Figure 7. Factor score plot for PC3.
Water 15 01961 g007
Table 1. Sampling locations and water bodies’ types.
Table 1. Sampling locations and water bodies’ types.
CodeDescriptionType *
A1Studena Dam (surface sample)L3
A2Struma River near Pernik town, after the Church quarter, before the mouth of the Rudarska RiverR5
A2-1Konska River before flowing into the Struma RiverR5
A3Struma River on the bridge near Batanovtsi town, after the Wastewater Treatment PlantR5
A3-1The Arkata River before flowing into the Struma RiverR13
A4Struma River at the bridge near Priboy village, after the confluence of the Arkata River, before the Pchelina DamR5
A5Pchelina Dam (surface sample)L13
A6Struma River after the Pchelina dam, near Zabliano villageR5
A7Struma River near Razdavitsa villageR5
B0-1Dragovishtitsa River near the borderR3
B0-2Sovolyanska Bistrica River before its mouthR5
B0-3Banshtitsa River after Kyustendil town, before flowing into the Struma RiverR5
B1Struma River near the village of NevestinoR5
B1-1Eleshnitsa river before its mouth, near Chetirtsi villageR3
B1-2Dyakovo Dam (surface sample)L13
B1-3Dupnishka Bistrica River after the village of Bistrica near the town of Dupnitsa before its mouthR5
B1-4Razmetanitsa River before its mouthR13
B1-5German River before flowing into the Struma River, bridge Boboshevo townR5
B2Struma River before Blagoevgrad townR5
B2-1Blagoevgradska Bistrica River after Blagoevgrad town before flowing into the Struma RiverR5
B3Struma River after Blagoevgrad townR5
B4Struma River after Simitli town, road. bridge on E79 in the city of OrlovetsR5
B4-1Brezhanska River before its mouth, Poleto villageR5
B5Struma River before Kresna townR5
B5-1Sandanska Bistrica River before its mouthR3
B6Struma River after the confluence of the Sandanska Bistrica River, after Sandanski townR5
B6-1Lebnitsa River before its mouthR3
B6-2UStrumeshnitsa River near the border (bridge to the village of Gabrene)R5
B6-2DStrumeshnitsa River before its mouth, bridge to Mitinovo villageR5
B7Struma River near the border with Greece (bridge to Topolnitsa village)R5
Notes: * Types of water bodies: R3—mountain river types; R5—semi-mountainous river types; R13—plain river types; L3, L13—types of “lakes” with oligotrophic conditions.
Table 2. Modified WQI categorization scheme for Bulgarian legislation compliance.
Table 2. Modified WQI categorization scheme for Bulgarian legislation compliance.
Status [52]WQINotesOriginal CCME-WQI Scale [53]
Very good80–100Clean and conditionally clean waters—the water quality meets the reference values for “Good” statusExcellent (95–100)
Good (80–94)
Good65–79Weakly polluted waters—the water quality randomly deteriorates from the reference values for “Good” statusFair (65–79)
Moderate0–64Polluted waters—the water quality does not meet the reference values for “Good” statusMarginal (45–64)
Poor (0–44)
Table 3. National Bulgarian guidelines for characterizing the ecological status of surface waters based on physicochemical indicators for quality [52].
Table 3. National Bulgarian guidelines for characterizing the ecological status of surface waters based on physicochemical indicators for quality [52].
Type *StatusDissO2,
mg/L
pHEC
µS/cm
NH4+
mg/L
NO3
mg/L
NO2
mg/L
N
mg/L
PO43−
mg/L
P
mg/L
BOD
mg/L O2
R3Very good10.5–8.0650<0.04<0.2<0.01<0.2<0.01<0.012<1
Good8.0–6.06.5–8.57500.04–0.40.2–0.50.01–0.0250.2–0.80.01–0.020.012–0.031–2.5
Moderate<6.0>750>0.4>0.5>0.025>0.8>0.02>0.03>2.5
R5Very good10.5–8.0700<0.04<0.5<0.01<0.5<0.02<0.025<1.2
Good8.0–6.06.5–8.57500.04–0.40.5–1.50.01–0.030.5–1.50.02–0.040.025–0.0751.2–3
Moderate<6.0>750>0.4>1.5>0.03>1.5>0.04>0.075>3
R13Very good9.0–7.0700<0.10<0.7<0.03<0.7<0.07<0.15<2
Good7.0–6.06.5–8.57500.10–0.30.7–20.03–0.060.7–2.50.07–0.150.15–0.32–4
Moderate<6.0>750>0.3>2>0.06>2.5>0.15>0.3>4
L3, L13Very good10.5–8.0650<0.03<0.2<0.01<0.20.007–0.0125<0.0125<1
Good8.0–6.06.5–8.77500.03–0.080.2–0.50.01–0.0250.2–0.80.0125–0.040.0125–0.041–2.5
Moderate<6.0>750>0.08>0.5>0.025>0.8>0.04>0.04>2.5
Note: * Types of water bodies as in Table 1.
Table 4. Basic statistics of the water quality indicators.
Table 4. Basic statistics of the water quality indicators.
IndicatorUnitnMeanMedianMinimumMaximumStandard
Deviation
DissO2mg/L13008.308.302.3915.992.12
pH13018.088.126.2411.390.45
ECμS/cm1300467.16424.0025.001883.00270.74
NH4+mg/L12840.450.100.0016.801.12
NO3mg/L12811.150.850.0010.351.19
NO2mg/L12720.050.020.002.700.10
Nmg/L12722.201.560.0039.002.36
PO43−mg/L12850.140.080.004.900.24
Pmg/L12770.220.120.005.300.36
BOD5mgO2/L12983.502.490.2579.004.65
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

Benkov, I.; Varbanov, M.; Venelinov, T.; Tsakovski, S. Principal Component Analysis and the Water Quality Index—A Powerful Tool for Surface Water Quality Assessment: A Case Study on Struma River Catchment, Bulgaria. Water 2023, 15, 1961. https://doi.org/10.3390/w15101961

AMA Style

Benkov I, Varbanov M, Venelinov T, Tsakovski S. Principal Component Analysis and the Water Quality Index—A Powerful Tool for Surface Water Quality Assessment: A Case Study on Struma River Catchment, Bulgaria. Water. 2023; 15(10):1961. https://doi.org/10.3390/w15101961

Chicago/Turabian Style

Benkov, Ivan, Marian Varbanov, Tony Venelinov, and Stefan Tsakovski. 2023. "Principal Component Analysis and the Water Quality Index—A Powerful Tool for Surface Water Quality Assessment: A Case Study on Struma River Catchment, Bulgaria" Water 15, no. 10: 1961. https://doi.org/10.3390/w15101961

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