3.1. Summary Statistics of Water Quality Parameters
Water samples were measured and collected at the bTR, TR, and aTR points (
Figure 1). The one-way ANOVA results showed that among the 17 parameters of water quality, significant statistical differences were noted for pH, NO
3-N, NO
2-N, and PO
4-P (
Table 4).
A higher standard deviation (SD) suggested greater differentiation of individual indicators. The SDs in all points were the highest for DSS and TSS and the lowest for heavy metals (Cu and Cr), and NO
2-N (
Table 4). SDs of NO
2-N, NH
4-N, and PO
4-P were <1 mg∙dm
−3, and SDs of TN and DO, BOD
5, and COD
Mn were >1 mg∙dm
−3 at all points. The SDs of certain parameters differed between points. The SD of NO
3-N was less than 1 mg∙dm
−3 only at bTR and TR, whereas the SDs of N
org and TP were less than 1 mg∙dm
−3 only at the bTR and aTR points.
It was found that the coefficient of variation (CV) of Cu at bTR and aTR, in addition to BOD
5, COD
Mn, TSS, Cr, and TP at TR, exceeded 100%. The remaining physico-chemical parameters showed smaller temporal changes (
Table 4). The CV coefficient also made it possible to evaluate the homogeneity of the analyzed indicators. Two indicators, pH and DSS at all points, in addition to Zn at bTR and DO at aTR, were more concentrated around the mean value, which made this coefficient much lower than for other indicators.
An investigation of water quality at three points in the Turawa Reservoir, (W_res), and below (Out_res) and above (In_res) this reservoir, covered the years 1998–2009, 2011, 2014, and 2016, indicated that the lowest SD value was obtained for the biogenic factor, in this case NO
2, and the highest for TSS. These results were consistent with the present study. Wiatkowski and Wiatkowska [
33] showed that the CV for TSS, (W-res), NO
2 (In_res), and PO
4 (Out_res) were characterized by large temporal changes. The number of samples in this study was greater; hence, the CV coefficients did not exceed 130%, as in the case of the present study in which the temporal changes in parameters were more visible [
33].
Physico-chemical parameters, such as DSS, NO
3-N, NO
2-N, TN, and Zn, obtained the highest mean concentrations in the water flowing into the reservoir (aTR point) (
Table 4). The highest mean concentrations of Tw, pH, DO, BOD
5, COD
Mn, TSS, and TP were recorded in the water stored in the reservoir (TR) (
Table 4). This is consistent with the results of Wiatkowski and Wiatkowska [
33]. The highest mean concentrations of NO
3, NO
2, TN, and DSS and additionally NH
4, PO
4, and TP were recorded in the water flowing into the reservoir (In_res). Moreover, the highest mean concentrations of Tw, TSS, pH, and DO occurred in the water stored in the reservoir (W-res). The explanation could be that water reservoirs break the river continuum and disrupt the transport of organic matter and metals in the river system, and ultimately affect the chemical and biological properties of the river ecosystem [
46,
47].
The mean concentrations of trace elements in the Turawa reservoir can be arranged in descending order: Zn > Cu > Cr. Similar results for these concentrations were obtained by Siepak and Sojka [
32]. They studied trace elements in two reservoirs: Kowalskie and Stare Miasto. However, the order below (bTR) and above (aTR) the Turawa reservoir is different: Zn > Cr > Cu.
The lowest water temperature (Tw) in the reservoir occurred in January 2020 (3.6 °C). Tw in the reservoir above 24 °C occurred in the months of June 2019 and August 2019. The investigation carried out in the period from June 2019 to June 2020 covered the normal period (SPI = 0.30), although 2019 was classified as a dry year (SPI = −0.8). In addition, dry and extremely dry periods were recorded in five months: June 2019, July 2019, November 2019, January 2020, and April 2020. Monthly periods in which measurements were performed had an impact on changes in water quality parameters in the river–reservoir system. The pH value recorded in the reservoir (TR) was above neutral—the lowest pH value was recorded in July 2019 (7.5 pH) and the highest in August 2019 (8.6 pH). The conditions were similar to those found in the research conducted by Hu et al. [
12] and Wiatkowski and Czerniawska-Kusza [
30]. This may be the cause of phosphorus release from sapropel sediments of the Turawa reservoir, which indicates an increased concentration of TP and PO
4-P in the reservoir waters (TR) and at the point below the reservoir (bTR). According to Hu et al. [
12], as a result of the mineralization of sedimentary organic matter and of the long-term accumulation of pollutants, the concentrations of nutrients in sediments are many times greater than those in water. When the external environment changed due to the increase in air temperature in summer and autumn, TP and NH
4-N were released from the sediment into the water, which would lead to excessive algae bloom and water quality degradation. Thus, the sediments act as a source of water pollution. The water analyzed in the Turawa reservoir, due to rising temperatures (25–27 °C) and elevated pH (8.12–9.95 pH), was characterized by the release of phosphorus from deeper sediment layers into the surface water [
30]. Wang et al. [
11] reported that both water temperature (25.3 °C in July) and pH (8.00–8.73 pH) significantly influenced the rate of phosphorus release from the sediment, but the effect of water temperature was more significant. In addition, Bartoszek et al. [
15] noted that the water in the reservoir does not have to be alkaline. As further emphasized by Wang et al. [
11], acid–base conditions favored the release of phosphorus from sediments, whereas the neutral environment did not. The mean NH
4-N concentration had the lowest value in the reservoir water (TR, 0.17 mg∙dm
−3), and the highest value at bTR (0.32 mg∙dm
−3) (
Table 4). By comparison, the mean concentration of N
org had the lowest value in the inflow to the reservoir (1.34 mg∙dm
−3) and the highest in the reservoir (TR, 2.66 mg∙dm
−3) (
Table 4). At the point below the reservoir (bTR), the mean concentration of N
org was lower than in the reservoir (1.62 mg∙dm
−3), but still higher than at the inflow to the reservoir (
Table 4). This is in line with the research by Hu et al. [
12]—when Tw increased at alkaline pH, N
org was transformed into ammonium nitrogen in the sediments, which promoted NH
4+-N to be released from sediments. In addition, the reduced oxygen caused by the intensifying mineralization of the sediments hindered the nitrification reaction in the sediments.
Wei et al. [
47] report that the Tw in the reservoir is often high and leads to a decrease in the amount of dissolved oxygen (DO) and a decrease in the rate of pollutant decomposition. However, the research of the Turawa reservoir (TR) showed that the mean value of Tw in the reservoir (TR) was higher than Tw in the Mała Panew river at the inflow (aTR) and at the outflow from the reservoir (bTR) and DO concentration in the TR had a higher value than in the river water. Nevertheless, the dissolved oxygen consumption in the reservoir (TR) was higher than that in the river, which indicates the mean value of BOD
5 and COD
Mn (
Table 4). Both indicators—BOD
5 and COD
Mn—as clearly indicate anthropogenic pollution [
28]. In addition, Maavara et al. [
10] report that damming modifies the ecological function of river systems. This indicates that the reservoir takes up nutrients and thus limits nutrient transfer downstream to flood plains. As reported by Wiatkowski [
13], depending on the depth of the reservoir and water retention, time, and the incoming load, up to 90% of all of the influencing matter in the water reservoir may be periodically retained.
Additionally, Kanownik et al. [
16] showed that in the waters of the Wisłoka river during the period of 2004–2013, the maximum values of most of the analyzed nutrients were recorded in June (TKN = 5.90 mg∙dm
−3, TN = 7.40 mg∙dm
−3, PO
4-P = 0.10 mg∙dm
−3, and TP = 0.76 mg∙dm
−3) and July (NO
2-N = 0.13 mg∙dm
−3). The opposite situation was recorded in the mountain river of Raba, where the highest average NO
3-N (2.71 mg∙dm
−3) concentrations occurred above the Dobczyce reservoir in the winter season. In turn, the highest average PO
4-P concentration was recorded below the reservoir (0.05 mg∙dm
−3) and the lowest concentration in the reservoir (0.03 mg∙dm
−3) [
23]. The NH
4-N concentration was highest in the inflow to the reservoir (0.98 mg∙dm
−3) and lowest in the outflow from the reservoir (0.14 mg∙dm
−3), which, according to Kijowska-Strugała et al. [
23], was caused by the absorption by autotrophs of ammonia in the reservoir because it is the main source of nutrients for phytoplankton. The investigation of the Rzeszów reservoir conducted by Bartoszek et al. [
14] during the summer season of 2013 showed that the reaction is alkaline; the values of Tw and PO
4-P were similar to the values in the current study (in summer), whereas the concentrations of TN and TP were lower. Zeinalzadeh and Rezaei [
2] obtained different results, finding that changes in water quality above and below the Shahr Chai reservoir in Iran did not show a significant difference. Furthermore, Rinke et al. [
22] provided information on data logging, at intervals of 15 min, for the following physico-chemical parameters: pH, conductivity, Tw, dissolved organic carbon (DOC), oxygen, TSS, NO
3, and chlorophyll. They showed that observations of the inflow water to the reservoirs over several months documented the highly temporal fluctuations which were caused by changes in water discharge.
3.2. Principal Component Analysis
Two principal components were obtained with eigenvalues >1, which accounted for 77.26% of the total variance in the bTR point and 66.86% in the TR point (
Table 5).
Table 6 presents the strength of the correlation for individual parameters, and
Figure 2 illustrates the first two principal components for each of these points.
At the bTR point, the first principal component (PC1) correlated the most highly with pH, NO
2-N, and NO
3-N, and weaklier with PO
4-P (
Table 6,
Figure 2). The NO
2-N and NO
3-N components appeared to have the lowest concentrations after transit through the reservoir (
Table 4). The nitrification and denitrification processes, which occur in the aquatic ecosystem are linked to microbial quantity and function [
48]. In the study on Chen et al. [
49] in Lake Taihu the cyanobacterial bloom promoted the efficacy of nitrogen removal in aquatic ecosystem. Therefore, the nitrogen concentration could be lower in the reservoir with the high degree of euthrofization [
49,
50]. It can be assumed, that in the Turawa reservoir occur similar processes, which could be confirmed in further studies. The second principal component (PC2) correlated most highly with PO
4-P (
Table 6,
Figure 2). As reported by Hu et al. [
12], the increased concentration of PO
4-P may be affected by bottom sediment of the reservoir (
Table 4). At the TR point, parameters were more correlated with PC1 and PC2, than with bTR and aTR points. At the TR point, the PC1 correlated most highly with NO
2-N and NO
3-N, whereas PC2 correlated most highly with pH, and PC3 correlated most highly with PO
4-P. Both NO
2-N and NO
3-N had lower concentrations at the TR point than in the water flowing into the reservoir (aTR), hence their strong correlation (
Table 6,
Figure 2). By comparison, pH and PO
4-P were not correlated with each other at the TR point (
Table 6,
Figure 2). This may indicate that the pH has no significant influence on the higher concentration of PO
4-P in the waters of the TR reservoir. At the aTR point, pH and PO
4-P were negatively correlated with each other, which indicates that they did not have such an impact on the water quality in the reservoir as nitrogen compounds (
Table 6,
Figure 2). The correlation between NO
2-N and NO
3-N at all points may indicate that they had the greatest effect on water quality at the analyzed points (
Figure 2).
The PCA analysis showed similar results of the PC1 and PC2 factors between the three analyzed points to those found in the study of Wiatkowski and Wiatkowska [
33]. Several parameters that influenced the variability of the water quality at the aTR point were the same as those observed by other researchers. As reported by Wiatkowski and Wiatkowska [
33], the first principal component (PC1) correlated most highly with NO
3 below, in, and above the Turawa reservoir. This indicator, like other nutrients, had the greatest impact on the qualitative differentiation of water in the Turawa reservoir. In addition, at the inflow and outflow from the reservoir, indicators such as pH and PO
4 were not concentrated around the PC1 axis, which indicates that they had a smaller impact on water quality at the analyzed points.
In the PCA presented by Zeinalzadeh and Rezaei [
2], the first four PCs had eigenvalues close to one. The eigenvalues accounted for 78.6%, 74.7%, 70.8%, and 80% of the total variability of the information contained in the original data set in spring, summer, autumn, and winter, respectively. From their study, it appears that changes in water quality at points above the reservoir and in the reservoir showed no significant difference. The point below the reservoir showed considerably worse water quality only in spring due to the dynamic development of the inhabitants and their commercial activities. In another study, water samples from four dam lakes: Grodzisk Duży, Karwacz, Łoje, and Wykrot, located in north-eastern Poland, from 2007 to 2013 were analyzed [
28]. This investigation showed that the parameters of the highest positive eigenvectors associated with PC1 included BOD
5 and COD. Negative relationships were found for SO
4 and NO
3. Cl and EC were negatively related to PC2, whereas DO and pH were positively related. In addition, in an investigation of 10 Polish reservoirs, Mazur et al. [
36] showed that, in the case of physico-chemical factors, first principal component (PC1) correlated most highly with TOC and TP. They found a significant association between TN and PC2, and EC and DO with PC3.
Similar studies using the PCA method can be found in the literature. However, this approach does not analyze the impact of reservoirs on water quality and flow, but rather studies the river systems. They are important because planning to build a reservoir can be associated with adequate water quality in the river. PCA is widely used to determine how different river points contribute to the overall pollution load [
7,
35,
51,
52]. Nazeer et al. [
3] showed that the most influential parameters for PC1 were PO
4, EC, salinity, Ca, Mg, K, pH, and DO. They studied the waters of the Soan River in Pakistan. In another study, of the Chillán river and in its tributaries, due to point pollutants, the first principal component (PC1) correlated with NH
4, PO
4, EC, COD, and DO (negatively correlated), and BOD
5. PC2 was highly influenced by the levels of NO
3 and NO
2 (both positively correlated) and PC3 Tw and pH (both positively correlated) [
1]. Siepak and Sojka [
32] used PCA, cluster analysis (CA), and discriminant analysis (DA) to assess the temporal and spatial variability in the investigated reservoirs. In the Kowalskie reservoir, Siepak i Sojka [
32] found that PC1 was most influenced by Cu, Zn, and Cd, which may indicate the natural origin of these components in water. They also showed that in the Old Town reservoir, PC1 was strongly positively correlated with the concentrations of Cr, Ni, Cu, and Cd.
3.3. Redundancy Analysis
The proportion of the variance explained was 26.94%; the first axis was 6% and was significant according to the RDA model (formula = data ~ site) (
p = 0.001) (
Table 7). The highest parameter scores of RDA1 were assigned to pH, NO
3-N, and NO
2-N (
Figure 3,
Table 8). Similar values of the variants were obtained by Zhou et al. [
39], who found that RDA could explain 20.13% of microbial community changes, with the contributions of RDA1 and RDA2 being 11.94% and 8.19%, respectively. Zhou et al. [
39] performed RDA and proved that physical (temperature, DO) and chemical (CHl-a, TOC, TN) parameters are the key factors that influence changes in the microbial community in the reservoir water.
The performed redundancy analysis showed that the distribution of parameters at the analyzed points (sites) (
Figure 3) was a grouping variable (the ellipsoid was defined by the 95% confidence interval), thus indicating the highest values of physicochemical parameters at a given site. The values of the NO
3-N and NO
2-N indicators followed the same trend and obtained the highest values at the same site (aTR). In addition, NO
2-N and NO
3-N concentrations were closely related to the proportions of aTR. At the bTR point, the highest values were obtained for PO
4-P. pH was closely related to the proportions of TR. For similar purposes, RDA was used by Palmer et al. [
38], Nazeer et al. [
3], and Bo [
34]. Palmer et al. [
38] analyzed water from Ontario lakes (Canada) and showed that lake development and road maintenance were associated with higher concentrations of Na and Cl ions and decreased conductivity in lake waters. The second axis of the RDA was positively correlated with the volume of the lake, runoff, and distance from the road [
38]. Nazeer et al. [
3] indicated that RDA analysis can be used to identify sources of pollution in river ecosystems because high scores correspond to high influence of the parameter at the sites. Bo et al. [
34] presented results of RDA analysis covering land use and water quality parameters. They found that the concentrations of COD and NH
4-N were closely related to the proportions of urban green areas and impermeable land, and BOD was closely related to the concentration of COD and the share of forests in the catchment area. However, ten years later, forest area was no longer in the group of dominant areas, and only COD and NH
4-N were closely related to urban green areas and impermeable land.
Other applications of RDA were presented by many scientists [
29,
37]. For example, Kowalczewska-Madura et al. [
37] performed RDA for P sedimented in bottom sediment and P released from bottom sediment of the eutrophic Durowskie Lake located in West Poland. This analysis showed that in lake water P was most negatively correlated with DO. Their investigation showed that the release of P depends partly on the N content in the sediment [
37]. On the basis of RDA, Ortiz-Vera et al. [
29] showed that the variability of physico-chemical parameters depended on the seasons. In the dry season, DO, NO
3, and pH concentrations appeared to be predictors of community structure, whereas in the rainy season temperature was also a predictor. Moreover, they observed greater differences in the composition of fungi communities in the dry season than in the rainy season [
29].
RDA showed that the RDA1 axis was the most positively correlated with NO
2-N and NO
3-N, and the most negatively correlated with pH, whereas RDA2 was the most negatively correlated with PO
4-P (
Table 8).
The results of this study were similar to those of the research conducted by Stendera and Johnson [
53], which show that TP determined in streams and lakes was negatively strongly associated with pastures and arable land.
3.4. Water Quality Profile of the Indicator Changes
Changes in the concentrations of the water quality parameters (TN, TP, NO
3-N, NO
2-N, NH
4-N, and PO
4-P) of the Mała Panew river at all analyzed points (bTR, TR, and aTR) are shown in
Figure 4. The analysis of this profile showed that the highest values of mean concentrations of TN and NO
3-N were recorded at the aTR station and the lowest at the bTR station.
The highest values of the median TP were recorded at the TR station, and the lowest values at the aTR station (
Figure 4). In addition, the lowest values NO
2-N and NH
4-N were recorded at the TR station, and the highest at the aTR and bTR stations, respectively (
Figure 4). The highest value of the median PO
4-P concentration was found at the bTR station and the lowest at the aTR station (
Figure 4). The profiles showed clearly that the concentrations of TN, NO
3-N, and NO
2-N decrease after passing through the Turawa reservoir. After the waters passed through the Turawa reservoir, the median concentrations of TP and NH
4-N, PO
4-P increased (
Figure 4). The water quality profile of the river–reservoir system, including the monitoring of the physico-chemical parameters, shows the influence of the reservoir on water quality during the analyzed period. This was also confirmed by Wiatkowski et al. [
31] in the research on the water quality of the Stobrawa river. Their investigation was performed before the construction of the Kluczbork reservoir.
Many scientists also indicated that the creation of a source profile of various factors helps to identify contaminated areas in a catchment, which may be useful in assessing the quality of waters in the catchment [
3,
11,
54]. Yan et al. [
54] provided a spatial diagram of polluted and risky sensitive zones in the basin based on the pollution index. Other researchers provided descriptive statistics according to seasonal changes [
7]. Debels et al. [
1] presented selected descriptive statistics in tabular form according to measuring stations. In the present study, the results of profile analysis are in line with those of the RDA. However, in our opinion, profile analysis is more suitable for analysis of the parameters at the sampling points. By comparison, RDA enables the correlations between the sampling points and parameters to be shown simultaneously. Therefore, profile analysis can be used to support the interpretation of RDA.
3.5. Water Quality Index
The WQI was determined for the period from June 2019 to June 2020 for three points: bTR, TR, and aTR. The index reached its minimum at aTR (WQI = 88), which indicates very poor water quality (
Table 9). The waters of the Mała Panew river at the point below the Turawa reservoir (bTR) and in the reservoir (TR) were classified as very polluted. The index peaked at TR. The WQI value at this point was the highest (WQI = 194) (
Table 9). However, at the point bTR, the value was WQI = 139 (
Table 9).
This was confirmed by the analyses of the average concentrations of water quality parameters at these three research points. The highest concentrations of TP, DO, and BOD5 occurred in the waters of the Turawa reservoir, which may be caused by the release of pollutants from sediment. Despite this, the WQI indicates that the analyzed index for waters below the reservoir (bTR) was lower than the index for waters in the reservoir (TR).
The WQI values of the points at the dam (3.1, 3.13,
Figure 1) are as follows: WQI
3.1 = 153 i WQI
3.13 = 158 (
Table 1). Moreover, the WQI values at points above the reservoir are twice as high and amount to WQI
3.7 = 253 i WQI
3.8 = 221 (
Table 10). This indicates that the WQI values determined from the water of the reservoir at the inlet and outlet are consistent with the values presented in the points aTR and bTR. The highest values were determined in the following bays: WQI
3.2 = 658, WQI
3.3 = 250 i WQI
3.12 = 529 (
Table 10). WQI determined for sample points on the left and right edges range from 100–200 (
Figure 1,
Table 10).
The WQI analysis provided by Imneisi and Aydin [
26] and Iticescu et al. [
27] concerned research of seasonal changes at the analyzed places. Imneisi and Aydin [
26] showed that the WQI index was the worst in January (poor water quality) and the best in September (excellent water quality), and good water quality was found in other months. According to Iticescu et al. [
27], the highest WQI index was in summer 2015 (WQI above 80), and the lowest in winter 2014 (WQI below 50). Tomczyk et al. [
19] emphasized that, in the assessment of water quality in a river–reservoir system, the physico-chemical properties are important in assessing the degree of water pollution.