Determination of Changes in the Quality of Surface Water in the River—Reservoir System

: proﬁle combined with PCA and RDA allows more effective monitoring for the needs of water management in the reservoir catchment area. The analyzed WQI for water below the reservoir (bTR) was lower than that of the reservoir water (TR), which indicates an improvement in water after passing through the reservoir.


Introduction
Intensive land use management, such as that undertaken by agriculture, industrial, and urban activities, significantly impact catchment areas, and particularly river valleys. Furthermore, land use management changes the quality and quantity of surface waters [1,2]. Particularly important is the nutrient increase (such as nitrogen and phosphorus), which poses a threat to the aquatic ecosystems of the downstream sections of rivers [1][2][3][4]. The sources of nutrients are mainly runoff from agricultural areas or point pollution of domestic wastewater [5]. The increase in nutrients in the water may lead to factors responsible for the eutrophication process of the Turawa reservoir. Additionally, the WQI will assess how the water quality changes in the analyzed river-reservoir system.

Study Area and Sampling Points
The Turawa reservoir (50°43′25″ N 18°07′13″ E) is located on the Mała Panew river in the south-western part of Poland, in Central Europe, in the Upper Odra river basin ( Figure  1). This reservoir was built during 1933-1939. It is a multipurpose reservoir. Its main tasks include flood protection, water supply to the power plant, provision of water for shipping purposes on the Odra river, fishing, and recreation [33]. Table 1 presents the parameters of the Turawa reservoir. The bottom of the reservoir contains about 4 million m 3 of sapropelic mud and less-contaminated sandy sediments [40]. According to Latała and Wierzba [41], the bottom sediment is in the form of dark grey fine-grained mud consisting of mineral and organic matter. The organic matter content comprises from 20% to 25% of the total mass of the bottom sediment. The bottom sediment is involved in the circulation of matter in the Turawa reservoir [30] and, consequently, has a significant impact on water quality [11,12]. The important aspect in the water quality assessment of the water reservoir is the microbiological characteristic [42,43].  50 3 Surface area at Normal Pool Level in the reservoir km 2 18. 29 4 Surface area at Maximum Pool Level in the reservoir km 2 19.41 (Source: [33]). This reservoir was built during 1933-1939. It is a multipurpose reservoir. Its main tasks include flood protection, water supply to the power plant, provision of water for shipping purposes on the Odra river, fishing, and recreation [33]. Table 1 presents the parameters of the Turawa reservoir. The bottom of the reservoir contains about 4 million m 3 of sapropelic mud and less-contaminated sandy sediments [40]. According to Latała and Wierzba [41], the bottom sediment is in the form of dark grey fine-grained mud consisting of mineral and organic matter. The organic matter content comprises from 20% to 25% of the total mass of the bottom sediment. The bottom sediment is involved in the circulation of matter in the Turawa reservoir [30] and, consequently, has a significant impact on water quality [11,12]. The important aspect in the water quality assessment of the water reservoir is the microbiological characteristic [42,43]. Source: [33].
The catchment area of the reservoir is 1423 km 2 . It is dominated by a forest area, which accounts for 59%. Arable land covers an area of 23%. Urban buildings, meadows, and pastures cover 5% of the area. The remaining use covers 8% of the catchment area. The average gradient of the Mała Panew valley is 1.58‰. The average temperature is 8.4 • C [33]. The average annual rainfall during the period 1986-2019 was 652 mm.
The average sum of the rainfall during the analyzed period was 684 mm. The highest amount of precipitation in the period from June 2019 to June 2020 was recorded in June 2020 (219 mm), and the lowest in April 2020 (6 mm) ( Table 2). The standardized precipitation index (SPI) for the analyzed period, from June 2019 to June 2020, indicates a normal period (SPI = 0.3). Own work based on IMWM-NRI 2020 and [44].
The analysis of the SPI index showed that the normal period in 2019 occurred only in three months (September, October, and December), and in 2020 it was characteristic only for March. The dry period occurred in four months: June 2019, July 2019, November 2019, and January 2020 (Table 2). In addition, the wet period occurred in three months (August 2019, February 2020, and May 2020), and the very dry and extremely wet periods each occurred on one occasion, in April 2020 and June 2020, respectively (Table 2).

Physico-Chemical Analyses
The 17 water quality parameters were measured twice per month from June 2019 to June 2020 in the mentioned stations of the Mała Panew river and Turawa reservoir. Measurements were taken at 13 points of the Turawa reservoir and one point each below and above the reservoir (Figure 1). A total of 360 samples were analyzed. Selection of water uptake was associated with the variability of the analyzed parameters in a river-reservoir system [33]. Temperature (Tw) and reaction (pH) were measured in situ, whereas dissolved oxygen (DO), biochemical oxygen demand (BOD 5 ), chemical oxygen demand (COD Mn ), ammonium nitrogen (NH 4 -N), dissolved substances (DSS), total suspended solids (TSS), nitrate nitrogen (NO 3 -N), nitrite nitrogen (NO 2 -N), organic nitrogen, total nitrogen (TN), phosphates (PO 4 -P), total phosphorus (TP), Zn, Cu, and Cr were determined in the laboratory (Table 3). Tw and pH were measured with Ezodo 7200. The time of the day affects the dissolved oxygen level. Water samples were measured and collected in the morning hours. The following equipment was used in the spectrophotometric methods: Thermo Scientific, Evolution 220 UV-Visible Spectrophotometer (for the determination of: NH 4 -N, NO 2 -N, NO 3 -N, PO4-P, TP, Kjeldahl nitrogen), DKL 12 Automatic Digestion Unit (for the determination of: TP, Kjeldahl nitrogen), UDK 139 Semi-Automatic Kjeldahl Distillation Unit (for the determination of Kjeldahl nitrogen). Kjeldahl nitrogen was needed to calculate total nitrogen (TN). The following measurement wavelengths were used: λ NH4 = 655 nm, λ NO2 = 540 nm, λ NO3 = 410 nm, λ PO4 = 880 nm, λ TP = 880 nm. The cuvette test was not performed as the applicable standards did not require it. For AAS Thermo Scientific iCE3000 Series AAS was used. The method uses a flame in the conditions of compressed air and acetylene. The measuring wavelengths are as follows: λ Zn = 213.9 nm, λ Cu = 324.7 nm, λ Cr = 357.9 nm. A three-point calibration mode was used. A certified standard solution was used in calibration. The standard curve was checked with ICP multi-element standard solution IV. Each sample was measured three times. The mean was taken as the final value.

SPI Analyses
The source of the precipitation data for the period of 1986-2020 is the Institute of Meteorology and Water Management-National Research Institute (IMWM-NRI). These data have been processed.

Statistical Study
All analyzes were performed in R software version 3.6.3. For the analyzed water quality indicators, basic statistics were calculated, such as maximum and minimum value, mean, standard deviation, and coefficient of variation.
The normality of 17 physico-chemical parameters was tested using the Shapiro-Wilk test for each parameter, which is standard procedure in water quality data assessment [2,33]. The water quality parameters with non-normal distribution were logarithmically transformed using log (x + 1).
One-way analysis of variance (one-way ANOVA), with a significance level of 5%, was performed to compare differences in water quality parameters (dependent variable) at different sampling sites (independent variable). One-way ANOVA calculations were made using the "stats" package. For PCA and RDA, only the data were used for which oneway ANOVA showed significant differences. For PCA and RDA, the "vegan" package was used [29]. For RDA, site one-way ANOVA was the independent variable. The significance of the axis was checked using the one-way ANOVA. The assumption of the water quality profile of the indicator changes is that the points of the analyzed places are marked on the horizontal axis, and the concentrations of the parameters are marked on the vertical axis [11]. The plot contains two pieces of information: boxplot (information on the minimum value, lower quantile, median, upper quantile, and maximum value) and the water quality profile along the watercourse [31].

Water Quality Index
In this study, the method for determining the WQI was provided by [24,26,27]. The WQI was calculated using the following mathematical expression: where W i is the unit weight for the ith water quality parameter, i represents the number of the parameters taken into consideration, and q i represents the rating for the ith water quality parameter. The WQI has a value between 0 and 100 (it can exceed 100 if the area is heavily polluted); q i was calculated using the following equation: where V i represents the measured value of the ith parameter; V 0 represents the ideal value of that parameter (it is 0 for all the parameters, with the exception of DO, for which the value is 14.6 mg dm −3 ); S i represents the value legally accepted for the water category in which the water sample was included. S i depends on the number of indicators required in the legal regulations for determining water quality classes in a given country. S i was determined based on [45] for the following indicators: DO, BOD 5 , NH 4 -N, NO 3 -N, PO 4 -P, TN, and TP. W i represents a factor which was calculated using the following formula: where K is a relative constant calculated by the formula:

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 Sustainability 2021, 13, 3457 9 of 18 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 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 physicochemical 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.

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.  Figure 2). This may indicate that the pH has no significant influence on the higher concentration of PO4-P in the waters of the TR reservoir. At the aTR point, pH and PO4-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 NO2-N and NO3-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 NO3 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 PO4 were not concentrated around the PC1 axis, which indicates that they had a smaller impact on water quality at the analyzed 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.

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 NO3-N and NO2-N indicators followed the same trend and obtained the highest values at the same site (aTR). In addition, NO2-N and NO3-N concentrations were closely related to the proportions of aTR. At the bTR point, the highest values were obtained for PO4-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 NH4-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 NH4-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  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.

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 riverreservoir 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. Explanations: *-maximum value. 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 NO2-N and NH4-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 PO4-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, NO3-N, and NO2-N decrease after passing through the Turawa reservoir. After the waters passed through the Turawa reservoir, the median concentrations of TP and NH4-N, PO4-P increased (Figure 4). The water quality profile of the river-reservoir

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  (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 BOD 5 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.

Conclusions
The water quality profile showed that TN, NO 2 -N, and NO 3 -N decreased after passing through the reservoir, whereas the values of phosphorus compounds (TP and PO 4 -P) increased. Despite this, the use of PCA and RDA showed that only NO 2 -N and NO 3 -N contributed to the observed variability in water quality in the river-reservoir system. Moreover, these analyzes showed that pH and PO 4 -P did not have the same influence on the water quality in the reservoir as the nitrogen compounds. In addition, the use of RDA indicated that the analyzed points in the river-reservoir system were of significant importance, showing that the values of the NO 2 -N and NO 3 -N indicators obtained the highest values at the aTR point, PO 4 -P at the bTR point, and pH at the TR point.
The investigation also showed that during the summer season in the dry months, high temperatures and alkaline reaction may cause the release of nitrogen and phosphorus compounds from sediment, which indicates an increased concentration of TP, PO 4 -P, and N org in the waters at the TR point and TP, PO 4 -P, and NH 4 -N concentrations at the bTR point. This finding is important in the analysis of the physico-chemical processes at the catchment scale and could be overlooked in conventional long-term analyses.
The WQI of the analyzed points indicated that the water of the Mała Panew river improved its properties after passing through the reservoir. The S i parameter was used to determine the WQI based on the DO, BOD 5 , NH 4 -N, NO 3 -N, PO 4 -P, TN, and TP indices.
The main finding is that both parameters (NO 2 -N and NO 3 -N) could be responsible for the eutrophication process. It is also evident that the Turawa reservoir reduces the concentration of nitrate nitrogen and nitrite nitrogen compared to the concentration of these compounds flowing into the reservoir.
This study has several practical implications: (1) the short-term sampling period may affect the correlations between the water quality in the inflow, in the reservoir, and in the outflow of the reservoir, and (2) the water quality profile combined with multidimensional PCA and RDA exploration techniques allows more effective monitoring to be implemented and improves the allocation of resources for the needs of water management at the catchment scale.

Informed Consent Statement: Not applicable
Data Availability Statement: The data presented in this study are available on request from the corresponding author.