Heavy Metals in Bottom Sediments of Reservoirs in the Lowland Area of Western Poland: Concentrations, Distribution, Sources and Ecological Risk

The paper presents the results of a study of heavy metals (HMs) concentrations in six retention reservoirs located in the lowland area of western Poland. The objectives of this study were to analyze the Cd, Cr, Cu, Ni, Pb and Zn concentrations, assess contamination and ecological risk, analyze the spatial variability of HM concentrations and identify potential sources and factors determining the concentration and spatial distribution. The bottom sediment pollution by HMs was assessed on the basis of the index of geo-accumulation (Igeo), enrichment factor (EF), pollution load index (PLI) and metal pollution index (MPI). To assess the ecological risk associated with multiple HMs, the mean probable effect concentration (PEC) quotient (Qm-PEC) and the toxic risk index (TRI) were used. In order to determine the similarities and differences between sampling sites in regard to the HM concentration, cluster analysis (CA) was applied. Principal component analysis (PCA) was performed to assess the impact of grain size, total organic matter (TOM) content and sampling site location on HM spatial distribution. Additionally, PCA was used to assess the impact of catchment, reservoir characteristics and hydrological conditions. The values of Igeo, EF, MPI and PLI show that Cd, Cr, Cu, Ni and Pb mainly originate from geogenic sources. In contrast, Zn concentrations come from point sources related to agriculture. The mean PEC quotient (Qm-PEC) and TRI value show that the greatest ecological risk occurred at the inlet to the reservoir and near the dam. The analysis showed that the HMs concentration depends on silt and sand content. However, the Pb, Cu, Cd and Zn concentrations are associated with TOM as well. The relationship between individual HMs and silt was stronger than with TOM. The PCA results indicate that HMs with the exception of Zn originate from geogenic sources—weathering of rock material. However, the Ni concentration may additionally depend on road traffic. The results show that a reservoir with more frequent water exchange has higher HMs concentrations, whereas the Zn concentration in bottom sediments is associated with agricultural point sources.


Introduction
Reservoirs play important roles in water supply, irrigation, power generation and flood control. However, dam construction changes the hydrological regime and sediment transport [1][2][3]. The suspended sediments are stored in the reservoir and with them pollutants including heavy metals (HMs) [4]. Bottom sediment plays an important role in monitoring the aquatic environment [5], especially evaluating the contamination levels [6,7], and ecological risk assessment [8].
Wiatkowski [86] and Kasperek and Wiatkowski [87] show that the heavy metal concentrations in reservoir bottom sediments vary across Poland. Concentrations of heavy metals in bottom sediments depend on the catchment land use and the occurrence of pollution point sources [25,88,89]. Therefore, each reservoir should be considered individually in order to expose the factors affecting the bottom sediment pollution. In this study, six reservoirs located in the lowland area of western Poland were selected for analysis. Their selection was influenced by the time of their construction, different morphometric parameters, catchment land use and hydrological conditions. The objectives of this study are to: (1) analyze the concentrations of HMs in bottom sediments, (2) show the spatial variability of heavy metal concentrations, (3) analyze the bottom sediment contamination, (4) assess ecological risk, (5) identify potential sources and factors determining the content and spatial distribution of HMs in the reservoirs.

Study Area and Sampling Sites
The research was carried out on the example of six reservoirs located in lowland areas in western Poland ( Figure 1). Wiatkowski [86] and Kasperek and Wiatkowski [87] show that the heavy metal concentrations in reservoir bottom sediments vary across Poland. Concentrations of heavy metals in bottom sediments depend on the catchment land use and the occurrence of pollution point sources [25,88,89]. Therefore, each reservoir should be considered individually in order to expose the factors affecting the bottom sediment pollution. In this study, six reservoirs located in the lowland area of western Poland were selected for analysis. Their selection was influenced by the time of their construction, different morphometric parameters, catchment land use and hydrological conditions. The objectives of this study are to: (1) analyze the concentrations of HMs in bottom sediments, (2) show the spatial variability of heavy metal concentrations, (3) analyze the bottom sediment contamination, (4) assess ecological risk, (5) identify potential sources and factors determining the content and spatial distribution of HMs in the reservoirs.

Study Area and Sampling Sites
The research was carried out on the example of six reservoirs located in lowland areas in western Poland ( Figure 1). The reservoirs are located in the areas most vulnerable to drought in Poland [90]. The main functions of the reservoirs are to provide water for irrigation purposes. The reservoirs were built between 1967 and 2013 (Table 1). Their surface area ranges from 33 to 91 ha and their volume from 0.3 to 2.4 million m 3 . They are shallow reservoirs with mean depths ranging from 1.3 to 2.1 m. The greatest depths from 3.5 to 5.5 m are found near the dam. The reservoirs Jeżewo, Jutrosin, Rydzyna, Środa and Września have an elongated shape and their shoreline development ranges from 1.6 to 2.9. The reservoirs are located in the areas most vulnerable to drought in Poland [90]. The main functions of the reservoirs are to provide water for irrigation purposes. The reservoirs were built between 1967 and 2013 (Table 1). Their surface area ranges from 33 to 91 ha and their volume from 0.3 to 2.4 million m 3 . They are shallow reservoirs with mean depths ranging from 1.3 to 2.1 m. The greatest depths from 3.5 to 5.5 m are found near the dam. The reservoirs Jeżewo, Jutrosin, Rydzyna,Środa and Września have an elongated shape and their shoreline development ranges from 1.6 to 2.9. However, the Pakosław reservoir has a rounded shape and the shoreline development is 0.8. In the catchment land use structure, agriculture area dominates, its percentage ranging from 68 to 81%. The proportion of urban and industrial areas is small and ranges from 2 to 7% in theŚroda, Rydzyna and Jutrosin catchments, respectively. The stream and road density in the catchment range from 1.2 to 2.3 km·km −2 and from 0.9 to 1.1 km·km −2 respectively. The number of river and road crossings where HM inflow may occur ranges from 0.5 to 0.8 no.·km −2 .

Sample Collection and Preparation
The sampling and field surveys took place during September 2017. The total 28 surface sediment samples were collected from six reservoirs ( Figure 1). The sampling locations were recorded using a hand-held Global Positioning System (GPS) device (GARMIN OREGON 300). Due to the elongated shape of the Jeżewo, Jutrosin, Radzyny,Środa and Września reservoirs, five samples were taken along the potential flow path. In the case of the Pakosław reservoir, with a rounded shape, three samples were taken. The top 0-5 cm depth of sediment samples was manually collected by means of Nurek and Czapla devices and carried within zip-mouthed PVC packages. In the laboratory, the bottom sediment samples were subjected to grain size and total organic matter (TOM) analysis. The analysis was carried out in order to examine the relationship between bottom sediment characteristics and total HM concentration. In the first step the bottom sediment samples were air dried at room temperature and sieved through a 2 mm nylon sieve to remove coarse debris. Next the samples were homogenized and passed through 2 mm, 1 mm, 0.5 mm, 0.25 mm, 0.10 mm and 0.05 mm sieves. The analysis was done by means of an analytical sieve shaker (Model, AS 200). The grain size fraction less than 0.05 mm was determined according to the Casagrande aerometric method modified by Prószyński [91]. The sand (0.05-2.0 mm), silt (0.002-0.05 mm), and clay (<0.002 mm) fractions were determined. The total organic matter (TOM) was measured by the Tiurin method with wet oxidation, followed by ferrous ammonium sulfate titration [91]. Finally, the samples were prepared for chemical analysis. The bottom sediment samples were extracted with hydrochloric acid (Merck, Darmstadt, Germany) at a 1:4 ratio at 95 • C ± 5 • C in a Mars 5 Xpress microwave digestion system (CEM, Matthews, NC, USA).

Chemical Analysis
The concentrations of Cd, Cr, Cu, Pb, Ni, and Zn in the bottom sediment samples with a fraction of <2 mm were determined by inductively coupled plasma mass spectrometry (ICP-QQQ), model 8800 Triple Quad (Agilent Technologies Tokyo, Japan). The analytical procedure is described in previous papers [16,73,75]. The reagents used were analytically pure, and the water was purified to the resistivity of 18.2 MΩ·cm (at 25 • C) in a Direct-Q 3 Ultrapure Water Systems apparatus (Millipore Molsheim, France). During the determinations by the ICP-QQQ techniques, standard solutions produced by VHG Labs, Inc. (Manchester, UK) were used. We measured the Certified Reference Material CRM no. LGC 6187 used for River sediments (Manchester, UK), and high compliance with reference values was found.

Quantification of Metal Pollution
The assessment of bottom sediment pollution by HMs was performed on the basis of four indices: the index of geo-accumulation (I geo ), the enrichment factor (EF) index, the pollution load index (PLI), and the metal pollution index (MPI). I geo and EF are single indices, while PLI and MPI are integrated indices.

Index of Geo-Accumulation
The index of geo-accumulation I geo was determined by the equation proposed by Müller [92] as follows (1): I geo is a single factor pollution index where C n and B n are the concentration of the selected HMs in bottom sediments and the geochemical background value, respectively (mg.kg −1 ). The factor 1.5 is used for the possible variability of the background data due to lithological conditions. Values of I geo were used to classify sediment samples into 7 purity classes: class 0 (I geo ≤ 0)-uncontaminated, class 1 (0 < I geo ≤ 1)-uncontaminated to moderately contaminated, class 2 (1 < I geo ≤ 2)-moderately contaminated, class 3 (2 < I geo ≤ 3)-moderately to heavily contaminated, class 4 (3 < I geo ≤ 4)-highly contaminated, class 5 (4 < I geo ≤ 5)-heavily to extremely contaminated, class 6 (5 < I geo )-extremely contaminated [64].

Enrichment Factor
The enrichment factor (EF) is an efficient index for quantification of the human impact on the concentration of a given metal [93]. The EF is based on a ratio of measured metal concentration to the metal geochemical background value. To account for natural heavy metal concentrations, EF is normalized to sediment Al or Fe content. In this study the normalization with Fe was selected, as proposed by Deely and Fergusson [94]. The EF was calculated as follows [95] (2): where (C n /C Fe ) is the ratio of the concentration of heavy metal (C n ) to iron (C Fe ) in the sediment sample and (B n /B Fe ) is the same ratio for the geochemical background [96,97]. Values of EF were used to assess the pollution of bottom sediment samples into the following classes: 0 (EF ≤ 1) no enrichment; 1 (1 < EF ≤ 3) is minor enrichment; 2 (3 < EF ≤ 5) is moderate enrichment; 3 (5 < EF ≤ 10) is moderately severe enrichment; 4 (10 < EF ≤ 25) is severe enrichment; 5 (25 < EF ≤ 50) is very severe enrichment; and 6 (EF > 50) is extremely severe enrichment [98,99].

Pollution Load Index
The pollution load index (PLI) is an experimental formula developed by Tomlinson et al. [100] (3): where C n /B n is the ratio of selected elements to their background values and i is the number of elements. The empirical index provides simple comparisons of site average heavy metal pollution. PLI value = 0 denotes perfection, PLI < 1 denotes no pollution, and PLI > 1 indicates pollution [101].

Metal Pollution Index
The metal pollution index (MPI) was also used to compare the total metal content between different sampling sites within reservoirs. The MPI was calculated as follows (4): where C n are concentrations of HMs in the sample and i is the number of analyzed metals [96,102]. The background heavy metal concentrations in the study area were provided by Bojakowska and Sokołowska [103]. The geochemical background values for bottom sediments were as follows 1.0 for Cd, 10.0 for Cu, 10.0 for Cr, 10.0 for Ni, 25.0 for Pb, 100.0 for Zn and 10,000.0 mg · kg −1 for Fe.

Ecological Risk Assessment
A consensus-based sediment quality guideline (SQG) was introduced by Macdonald et al. [104] to predict the toxicity of a sediment samples. For each chemical contaminant there are two consensus-based values-threshold effect concentration (TEC) and probable effect concentration (PEC). The concentration of chemical contaminant below the TEC means that undesirable effects are not expected to occur, while above the PEC adverse effects are expected to occur more often than not [39,105,106]. To assess the effects of multiple HMs, the mean PEC quotient (Q m-PEC ) was proposed, which can be calculated using Equation (5): where n is the number of HMs, Cn is the measured concentration of a heavy metal and PEC n is the corresponding PEC. The PEC benchmark values for Cd, Cr, Cu, Ni, Pb and Zn are 4.98, 111, 149, 48.6, 128, and 459 mg·kg −1 , respectively. The categories for Q m-PEC are classified as: not toxic (Q m-PEC < 0.5) and toxic (Q m-PEC > 0.5). The toxic risk index (TRI) was developed by Zhang et al. [107] to assess the toxic risk for certain HMs in sediment samples. The TRI is based on threshold (TEL) and probable (PEL) level effects [108]; it can be calculated using Equation (6): To assess integrated risk of multiple HMs, the TRI can be calculated using Equation (7): where n is the number of HMs, C i is the content of each heavy metal in the sediment sample, and TRI i is the toxic risk index of each heavy metal. Five categories of the TRI are classified as: no toxic risk (TRI ≤ 5), low (5 < TRI ≤ 10), moderate (10 < TRI ≤ 15), considerable (15 < TRI ≤ 20), and very high (TRI > 20) [43].

Distribution and Source Identification
In order to determine the similarities and differences between sampling sites in regard to the HM concentration, cluster analysis (CA) was applied. The CA was carried out using the Ward method with square Euclidean distance as a measure of similarity between clusters. Groups and subgroups were distinguished using the cut-off criteria of 66% and 25% respectively [109]. For the purpose of comparison of HM concentrations between sample sites, the non-parametric Kruskal-Wallis test (K-W, p ≤ 0.05) followed by Dunn's test as a post-hoc procedure was used. The K-W test was used for verification of the hypothesis of significance of differences between concentrations of HMs in distinguished groups.
Subsequently, to choose the appropriate ordination methods for future analysis, detrended component analysis (DCA) was conducted. The obtained gradients of the first axis of the DCA were shorter than 3.0 SD, which indicates that the HMs were in linear distribution. Glińska-Lewczuk et al. [110] recommended in such cases principal component analysis (PCA). PCA was employed to identify a potential source of HMs and factors that have the greatest impact on their concentration and spatial distribution. The PCA analysis was done in two variants. First the impact of grain size, total organic matter content and sampling site location on HMs' spatial distribution was evaluated. The percentage content of sand (SA), silt (SI), clay (CL) and total organic matter (TOM) was calculated for each sample. In addition, the distance between sampling site location with respect to inlet (ID) and outlet (OD) of the reservoir was calculated. Second, the impact of other factors was analyzed. The following factors were distinguished: catchment characteristics (catchment area-CA, mean slope-MS, mean elevation-ME, urban and industrial-UI, agriculture-A, forest-F, stream density-DD, road density-RD, number of road and river crossings-RRC), reservoir characteristics (age-AG, area-AR, mean depth-MD, volume-VO, shoreline development-SD), and hydrological conditions (retention time-RT). The above factors were referred to the mean HM concentration in the samples located in the reservoir (J, Ju, P, R, S, W). Before starting the PCA analysis the factors were tested in the context of the occurrence of outliers. The analysis was performed using the two-sided Grubbs test at the significance level of 0.05. Subsequently the HM concentration and factors that may have an impact on HM content and distribution were log-transformed to obtain a normal distribution. In order to avoid misclassification due to differences in data dimensionality, multivariate statistical techniques were applied on standardized data through z-scale transformation [73,74]. The significant principal components were selected based on a Kaiser criterion with eigenvalues higher than 1. The correlation between principal components and analyzed data was classified according to the Liu et al. [111] criterion where values >0.75, 0.75-0.50, and 0.50-0.30 show a strong, moderate and weak relationship, respectively. CA and PCA were conducted using Statistica 13.1 and Canoco 5.0 respectively.

Heavy Metal Content
The characteristic concentrations of HMs (Cd, Cr, Cu, Ni, Pb and Zn) in bottom sediments from 28 sites are shown in Table 2. Generally, the mean concentrations of HMs in reservoirs follow a descending order of Zn > Pb > Cu > Cr > Ni > Cd. The Zn average concentration in reservoirs' bottom sediments was 452.2 mg·kg −1 . The lowest Zn concentrations occurred in the Jutrosin and the highest in the Jeżewo reservoir, 23.1 and 903.7 mg·kg −1 respectively. It was observed that in the Jeżewo, Rydzyna,Środa and Września reservoirs, Zn concentrations higher than 1000 mg·kg −1 were observed at the inlet to the reservoir or near the dam. Pb concentrations at the highest level were in the Jeżewo reservoir (average 17.6 mg·kg −1 ), while the lowest values were observed in Pakosław (average 2.6 mg·kg −1 ). A similar tendency was observed in the case of Cu, Cr and Ni concentrations. The average concentrations of Cu, Cr and Ni in the Pakosław reservoir were at a similar level and was on average 2.0 mg·kg −1 . In the Jeżewo reservoir, the average concentrations of these elements in bottom sediments were 3 to 5 times higher. The Cd concentrations in the reservoirs were the lowest and ranged from 0.01 to 0.72 mg·kg −1 . In the Września reservoir, the average concentration of Cd was about seven times higher than in the Pakosław reservoir. It was observed that Zn concentrations are characterized by the highest variability between the sampling sites and reservoirs; however, Cr concentrations have the lowest variability. The non-parametric Kruskal-Wallis test showed significant differences in Cd, Cu, Zn and Pb concentrations between bottom sediments at individual reservoirs. Differences were not significant only for Cr and Ni concentrations.

Bottom Sediment Pollution Assessment
The I geo was applied to assess the HM contamination level in bottom sediments. The I geo values for Cd, Cr, Pb and Ni were below zero and could be classified as unpolluted to moderately polluted. The I geo values for Cu were also lower than 0 except for two sampling sites, J1 and W3. At these points, I geo values were slightly higher than 0, which indicates a small local contamination of bottom sediments with Cu. The largest contamination of Zn was noted in the reservoirs. I geo values for Zn lower than 0 occurred only for the Jutrosin reservoir. In 23 measurement sampling sites located in other reservoirs, I geo values higher than 0 were recorded fourteen times. The sediment bottom at sample sites R2, R3 and S2 demonstrates weak contamination. The values of I geo were slightly higher than 0. The sediments in sampling sites P5, R5 and S4 demonstrate moderate contamination with the value of I geo higher than 1. The I geo for sampling sites J1, J2, J3, J5, R1, S5, W1 and W5 was between 2 and 4, which demonstrated moderate to severe pollution of bottom sediments (Figure 2a).

Bottom Sediment Pollution Assessment
The Igeo was applied to assess the HM contamination level in bottom sediments. The Igeo values for Cd, Cr, Pb and Ni were below zero and could be classified as unpolluted to moderately polluted. The Igeo values for Cu were also lower than 0 except for two sampling sites, J1 and W3. At these points, Igeo values were slightly higher than 0, which indicates a small local contamination of bottom sediments with Cu. The largest contamination of Zn was noted in the reservoirs. Igeo values for Zn lower than 0 occurred only for the Jutrosin reservoir. In 23 measurement sampling sites located in other reservoirs, Igeo values higher than 0 were recorded fourteen times. The sediment bottom at sample sites R2, R3 and S2 demonstrates weak contamination. The values of Igeo were slightly higher than 0. The sediments in sampling sites P5, R5 and S4 demonstrate moderate contamination with the value of Igeo higher than 1. The Igeo for sampling sites J1, J2, J3, J5, R1, S5, W1 and W5 was between 2 and 4, which demonstrated moderate to severe pollution of bottom sediments (Figure 2a). The enrichment factor (EF) was applied to assess the sources of HMs. The results of the EF calculations for Cd, Cr, Cu, Pb and Ni showed that the values generally did not exceed a value of 1 (Figure 2b). EF values ranging from 1 to 3 for Pb occurred once, for Cd and Cr twice, for Ni four times and for Cu nine times. In most cases, the values were less than 1.5, which may suggest that the HMs may come from natural weathering processes. In the case of Cr, at points S2 and W5 and for Cu at points S5 and W5, the EF values were higher than 1.5. However, an EF value higher than 1.5 indicates that HMs were delivered by other sources, such as point and non-point pollution [112,113]. The enrichment factor (EF) was applied to assess the sources of HMs. The results of the EF calculations for Cd, Cr, Cu, Pb and Ni showed that the values generally did not exceed a value of 1 ( Figure 2b). EF values ranging from 1 to 3 for Pb occurred once, for Cd and Cr twice, for Ni four times and for Cu nine times. In most cases, the values were less than 1.5, which may suggest that the HMs may come from natural weathering processes. In the case of Cr, at points S2 and W5 and for Cu at points S5 and W5, the EF values were higher than 1.5. However, an EF value higher than 1.5 indicates that HMs were delivered by other sources, such as point and non-point pollution [112,113]. The EF value of Zn was highest, and it ranged from 0.22 to 43.52. In fourteen sampling sites the EF values were higher than 3. Only in the Jutrosin reservoir were the EF values lower than 1. EF values for Zn higher than 10 were observed in the Jeżewo, Pakosław, Rydzyna,Środa and Września reservoirs, generally at sampling sites located at the inlet to the reservoir and at the dam. Two global indices of PLI and MPI were used to assess the bottom sediments' contamination with HMs. PLI values ranged from 0.06 to 1.54. PLI values higher than 1 were noted four times at points J1, J5, W1 and W3 (Figure 3).

Ecological Risk Assessment
The toxic risk index was used to assess the toxic risk for HMs in bottom sediment samples. The values of the TRI index ranged from 0.17 to 13.56. In general, the TRI values were lower than 5 in 20 out of 28 sampling sites, which indicates a lack of ecological risk associated with HM concentrations in bottom sediments. Values lower than 5 occurred at all sampling sites located in Jutrosin reservoirs. TRI values ranging from 5 to 10 occurred six times (J1, J2, J3, R1, S5 and W5), which indicated low ecological risk. Only in sampling sites J5 and W1 were the TRI values higher than 10, which indicated moderate ecological risk. The highest contribution in the TRI index were Zn concentrations (Figure 4a), on average 59%. The average TRI value for the Jutrosin reservoir was the lowest at 0.64, while more than ten times higher values were found in the Jeżewo reservoir. The reservoirs according to TRI values can be arranged in ascending order Jutrosin < Pakosław < Środa < Rydzyna < Września < Jeżewo.
For the toxicity assessment of combined HMs, a consensus-based sediment quality guideline (SQG) was used. The Qm-PEC values ranged from 0.07 to 4.88. The average Qm-PEC values of the reservoirs ranked as follows: Jeżewo > Września > Rydzyna > Środa > Pakosław > Jutrosin. Qm-PEC values higher than 0.5 were noted up to seventeen times, which may indicate that HMs can be toxic to certain sediment-dwelling organisms [104]. Values higher than 0.5 occurred at all sampling sites located in the Jeżewo reservoir, four times in the Rydzyna reservoir and three times in the Środa and Września reservoirs (Figure 4b). The contribution of Zn was the highest in the Qm-PEC index of the analyzed sampling sites. However, at the sampling sites located in the Jutrosin reservoir and points W2, W3 and W4 (Września Reservoir), R4 (Rydzyna Reservoir) and S3 (Środa Reservoir), the

Ecological Risk Assessment
The toxic risk index was used to assess the toxic risk for HMs in bottom sediment samples. The values of the TRI index ranged from 0.17 to 13.56. In general, the TRI values were lower than 5 in 20 out of 28 sampling sites, which indicates a lack of ecological risk associated with HM concentrations in bottom sediments. Values lower than 5 occurred at all sampling sites located in Jutrosin reservoirs. TRI values ranging from 5 to 10 occurred six times (J1, J2, J3, R1, S5 and W5), which indicated low ecological risk. Only in sampling sites J5 and W1 were the TRI values higher than 10, which indicated moderate ecological risk. The highest contribution in the TRI index were Zn concentrations (Figure 4a), on average 59%. The average TRI value for the Jutrosin reservoir was the lowest at 0.64, while more than ten times higher values were found in the Jeżewo reservoir. The reservoirs according to TRI values can be arranged in ascending order Jutrosin < Pakosław <Środa < Rydzyna < Września < Jeżewo.
For the toxicity assessment of combined HMs, a consensus-based sediment quality guideline (SQG) was used. The Q m-PEC values ranged from 0.07 to 4.88. The average Q m-PEC values of the reservoirs ranked as follows: Jeżewo > Września > Rydzyna >Środa > Pakosław > Jutrosin. Q m-PEC values higher than 0.5 were noted up to seventeen times, which may indicate that HMs can be toxic to certain sediment-dwelling organisms [104]. Values higher than 0.5 occurred at all sampling sites located in the Jeżewo reservoir, four times in the Rydzyna reservoir and three times in theŚroda and Września reservoirs (Figure 4b). The contribution of Zn was the highest in the Q m-PEC index of the analyzed sampling sites. However, at the sampling sites located in the Jutrosin reservoir and points W2, W3 and W4 (Września Reservoir), R4 (Rydzyna Reservoir) and S3 (Środa Reservoir), the percentages of HMs were at a similar level. It was observed that in the reservoirs, the greatest ecological risk associated with the pollution of HM deposits occurred at the inlet to the reservoir and near the dam.

Distribution and Source Identification
The CA analysis allowed the sampling sites to be divided into two groups, A and B, depending on the HM concentrations (Figure 5a). A total of 19 sampling sites were classified in group A and 9 in group B. The concentrations of all HMs in group A were lower than those observed in group B. In group B there are four measurement points located in the Jeżewo reservoir, three in the Września reservoir and two in the Środa reservoir. Group B was characterized by low variability of concentrations of individual HMs. Higher variability occurred in group A. Therefore, this group was divided into two subgroups: A1 and A2. Four sampling sites-S1, P1, Ju1 and P3-were classified in subgroup A1, in which the lowest concentrations of Cr, Ni, Cu and Pb were observed. However, in the case of Cd and Zn concentrations the differences between subgroups A1 and A2 were not statistically significant. Another division of the sampling sites into groups was obtained depending on the sediments' bottom texture and TOM (Figure 5b). Samples were divided into two groups, I and II, which included 3 and 25 sampling sites, respectively. Sampling sites in group I were characterized by the highest content of TOM and silt fraction, on average 8.5% and 49.7%. At these sampling sites, the lowest sand fraction was observed, on average 47.0%. Within group II, 3 subgroups-IIa, IIb and IIc-were separated, into which 9, 3 and 13 sampling sites were classified respectively. In subgroup IIa there were sampling sites with the highest sand content (mean 92.9%) and the lowest content of silt and clay (average 6.7% and 0.4%). Also, the content of TOM was at the lowest level: 0.4%. Subgroup IIb was characterized by the highest content of clay (average 7.7%) and high content of sand (84.7%). In subgroup IIc the clay content was at a low level (on average 1.2%) and silt content at a high level (on average 21.8%). Dunn's test showed that in groups I and II there are significant differences between sand, silt, clay and TOM content. Differences in all analyzed parameters were also present between subgroups IIa and IIc. However, between subgroups IIa and IIb there were only differences in the content of sand and clay and between subgroups IIb and IIc in the content of silt and clay. It was observed that the sample sites with the highest content of HMs (group B- Figure 5a) belonged to group I and subgroup IIc (Figure 5b). In these groups, high content of silt and TOM in group I was recorded.

Distribution and Source Identification
The CA analysis allowed the sampling sites to be divided into two groups, A and B, depending on the HM concentrations (Figure 5a). A total of 19 sampling sites were classified in group A and 9 in group B. The concentrations of all HMs in group A were lower than those observed in group B. In group B there are four measurement points located in the Jeżewo reservoir, three in the Września reservoir and two in theŚroda reservoir. Group B was characterized by low variability of concentrations of individual HMs. Higher variability occurred in group A. Therefore, this group was divided into two subgroups: A1 and A2. Four sampling sites-S1, P1, Ju1 and P3-were classified in subgroup A1, in which the lowest concentrations of Cr, Ni, Cu and Pb were observed. However, in the case of Cd and Zn concentrations the differences between subgroups A1 and A2 were not statistically significant. Another division of the sampling sites into groups was obtained depending on the sediments' bottom texture and TOM (Figure 5b). Samples were divided into two groups, I and II, which included 3 and 25 sampling sites, respectively. Sampling sites in group I were characterized by the highest content of TOM and silt fraction, on average 8.5% and 49.7%. At these sampling sites, the lowest sand fraction was observed, on average 47.0%. Within group II, 3 subgroups-IIa, IIb and IIc-were separated, into which 9, 3 and 13 sampling sites were classified respectively. In subgroup IIa there were sampling sites with the highest sand content (mean 92.9%) and the lowest content of silt and clay (average 6.7% and 0.4%). Also, the content of TOM was at the lowest level: 0.4%. Subgroup IIb was characterized by the highest content of clay (average 7.7%) and high content of sand (84.7%). In subgroup IIc the clay content was at a low level (on average 1.2%) and silt content at a high level (on average 21.8%). Dunn's test showed that in groups I and II there are significant differences between sand, silt, clay and TOM content. Differences in all analyzed parameters were also present between subgroups IIa and IIc. However, between subgroups IIa and IIb there were only differences in the content of sand and clay and between subgroups IIb and IIc in the content of silt and clay. It was observed that the sample sites with the highest content of HMs (group B- Figure 5a) belonged to group I and subgroup IIc (Figure 5b). In these groups, high content of silt and TOM in group I was recorded.

Relationships between HM Concentrations and Sediment Texture and Organic Matter Content
In the first stage, the aim of the PCA analysis was to identify factors that determine the concentration of HMs in the bottom sediments. The PCA analysis allowed two significant factors to be distinguished, PC1 and PC2, whose eigenvalues were higher than 1. The PC1 explains as much as 75.9% of the variance. The PC1 was strongly negatively correlated with Cd, Cr, Cu, Ni and Pb concentrations ( Figure 6).

Relationships between HM Concentrations and Sediment Texture and Organic Matter Content
In the first stage, the aim of the PCA analysis was to identify factors that determine the concentration of HMs in the bottom sediments. The PCA analysis allowed two significant factors to be distinguished, PC1 and PC2, whose eigenvalues were higher than 1. The PC1 explains as much as 75.9% of the variance. The PC1 was strongly negatively correlated with Cd, Cr, Cu, Ni and Pb concentrations ( Figure 6).

Relationships between HM Concentrations and Sediment Texture and Organic Matter Content
In the first stage, the aim of the PCA analysis was to identify factors that determine the concentration of HMs in the bottom sediments. The PCA analysis allowed two significant factors to be distinguished, PC1 and PC2, whose eigenvalues were higher than 1. The PC1 explains as much as 75.9% of the variance. The PC1 was strongly negatively correlated with Cd, Cr, Cu, Ni and Pb concentrations ( Figure 6). PC2, explaining about 11.7% of the variance, was strongly negatively correlated with Zn concentrations. The projection of sampling sites perpendicular to individual HM vectors or their extension beyond the origin of the coordinate system allows estimation of their content in bottom sediments. For example, the highest Zn values occur at sampling sites W1, J5, J1, S5, W5 and R1 and the lowest at sampling sites Ju3, Ju2, Ju1 and P1. The PC1 factor was moderately positively correlated with sand content and negatively with silt and TOM content. The angle created by the vectors of individual HMs informs about the relationship between them. The smaller the angle is, the stronger is the relationship; however, the angle value of 90 • means no correlation. The analysis showed that the concentrations of HMs in bottom sediments are positively correlated with the silt content and negatively with the sand content (except Zn), while no correlations were observed with clay. In addition, the studies showed a positive correlation between Pb, Cu, Cd and Zn concentrations and TOM content. However, the relationship between individual HMs and silt was stronger than with TOM. A positive correlation was observed between the sampling site distance from the inlet to the reservoir (ID) and the silt and TOM content and a negative correlation with the sand content. This indicates grain size segregation in reservoirs. Thicker, sandy fractions are deposited at the inlet of the reservoirs, while the silt and TOM fractions are deposited near the dam.
The second objective of the PCA was the identification of external factors that may affect the concentrations of HMs in reservoirs. The PCA analysis made it possible to distinguish two factors, PC1 and PC2, which explain 79.9 and 14.5% of variance respectively. The average concentrations of Cd, Cr, Cu, Ni and Pb were strongly negatively correlated with PC1, whereas Zn concentrations were strongly correlated with PC2 ( Figure 7). Among the other factors that may have an impact on the content of HMs in bottom sediments, PC1 was strongly positively correlated with mean catchment slope and strongly negatively correlated with reservoir shoreline development. In addition, a positive correlation was observed with the catchment area and the water retention time in the reservoir. Only the share of agricultural areas was negatively correlated with PC2. The results indicate that HMs with the exception of Zn have the same source. Low concentrations in bottom sediments indicate that they originate from geogenic sources-weathering of rock material. The positive correlation with reservoir shoreline development indicates a possible impact of adjacent areas. The content of Ni in bottom sediments may additionally result from road traffic, which confirms the correlation with number of road and river crossings. The negative correlation of Cd, Cr, Cu, Pb and Ni with retention time indicates that in the reservoir with more frequent water exchange, HM concentrations are higher, which is related to the suspended sediment deposition, whereas the Zn in bottom sediments is correlated with agricultural land use. High concentrations of Zn in bottom sediments may result from the existence of anthropogenic sources and a potential supply with domestic wastewater. The data provided by the Central Statistical Office in Poland indicate that the number of people using the water supply network in the reservoir catchment varied from 93 to 98%. On the other hand, the number of people using the sewerage network varied from 36 to 65% for the Rydzyna and Września catchment respectively. The ratio of the length of the water supply to the length of the sewage system varied from 0.21 to 0.35. The above data indicate possible supply of Zn with domestic wastewater. The second objective of the PCA was the identification of external factors that may affect the concentrations of HMs in reservoirs. The PCA analysis made it possible to distinguish two factors, PC1 and PC2, which explain 79.9 and 14.5% of variance respectively. The average concentrations of Cd, Cr, Cu, Ni and Pb were strongly negatively correlated with PC1, whereas Zn concentrations were strongly correlated with PC2 ( Figure 7). Among the other factors that may have an impact on the content of HMs in bottom sediments, PC1 was strongly positively correlated with mean catchment slope and strongly negatively correlated with reservoir shoreline development. In addition, a positive correlation was observed with the catchment area and the water retention time in the reservoir. Only the share of agricultural areas was negatively correlated with PC2. The results indicate that HMs with the exception of Zn have the same source. Low concentrations in bottom sediments indicate that they originate from geogenic sources-weathering of rock material. The positive correlation with reservoir shoreline development indicates a possible impact of adjacent areas. The content of Ni in bottom sediments may additionally result from road traffic, which confirms the correlation with number of road and river crossings. The negative correlation of Cd, Cr, Cu, Pb and Ni with retention time indicates that in the reservoir with more frequent water exchange, HM concentrations are higher, which is related to the suspended sediment deposition, whereas the Zn in bottom sediments is correlated with agricultural land use. High concentrations of Zn in bottom sediments may result from the existence of anthropogenic sources and a potential supply with domestic wastewater. The data provided by the Central Statistical Office in Poland indicate that the number of people using the water supply network in the reservoir catchment varied from 93 to 98%. On the other hand, the number of people using the sewerage network varied from 36 to 65% for the Rydzyna and Września catchment respectively. The ratio of the length of the water supply to the length of the sewage system varied from 0.21 to 0.35. The above data indicate possible supply of Zn with domestic wastewater.

Discussion
The HM content in bottom sediments of retention reservoirs and relations between them result mainly from the impact of anthropogenic sources [14] as well as sediment texture and organic matter content [26]. Meteorological, hydrological and geological conditions of the catchment affect the

Discussion
The HM content in bottom sediments of retention reservoirs and relations between them result mainly from the impact of anthropogenic sources [14] as well as sediment texture and organic matter content [26]. Meteorological, hydrological and geological conditions of the catchment affect the content of sand, silt, clay and TOM in water. These factors influence the transport of HMs in the river-reservoir system. In the world, bottom sediments of retention reservoirs are characterized by relatively high variability of HM concentration. The analyzed reservoirs are located in agricultural catchment; therefore, HM concentrations in bottom sediments are at a lower level. Only Zn concentrations are at a higher level with respect to those noted by other authors in Polish reservoirs [25]. Probably Zn content may result from potential supply with domestic wastewater. The results presented in this paper confirmed the results obtained by Haziak et al. [114] which indicated that Zn values exceeding the geochemical background were associated with anthropogenic activities (run-off from farmland and domestic wastewater).
High values of the I geo , EF, PLI and MPI indices at individual sampling sites of the reservoirs may suggest their origin of HMs from anthropogenic sources. High values at the inlet to the reservoir result from their inflow from point sources. Bing et al. [43] suggested that local anthropogenic activities increase contamination of specific HMs in the bottom sediments, especially of upper regions. Saleem et al. [55] observed relatively high content of HMs at sites which were adjacent to urban and semi-urban areas. The HMs may also come from the discharges of untreated urban/industrial wastes, agricultural runoffs and automobile emissions. High HM concentration near the dam is a consequence of changes in hydrological conditions in reservoirs during floods and high flow periods. HMs absorbed on silt and clay and TOM are transported and deposited near the dam. Also, Bing et al. [43] and Frémion et al. [17] reported relatively high values of HM concentrations close to the dam. Sojka et al. [25] noted that near the outlet pipe from the bottom sediments there can be discharged silt, clay and TOM, which are responsible for the transport of HM downstream from the reservoir.
The results have shown that concentrations of Cd, Cr, Cu, Ni, and Pb were strongly positively correlated with the silt and TOM contents in bottom sediments. However, the relationship with silt was at a higher level. Palma et al. [52] obtained slightly different results which showed stronger relationships between HMs and silt and clay content. However, the relationship between HMs and silt was stronger. Farhat and Aly [28] suggest that organic matter was more critical than grain size in controlling HM distribution in sediment. Different results were obtained by Frémion et al. [23], who found no relationships between trace elements and organic fraction.
The content in the bottom sediments of HMs involves the possibility of ecological risks. It was observed that the ecological risk in the reservoirs occurred at the inlet and near the dam, where the highest Q m-PEC and TRI values were obtained. Slightly different results were obtained by Bing et al. [43], who found that the potential ecological risk level increased towards the dam. Results obtained by Zhao et al. [50] suggested that the concentration and potential ecological risk of HMs in sediments near the dam were higher compared to the upstream sampling sites. The probable reason for the very high HM content at the reservoir's inlet was the date of sampling during the very low flow period. The obtained values of Q m-PEC and TRI were the most strongly correlated with TOM. Our research confirmed the results obtained by Lin et al. [29] that the organic matter-bound fraction of heavy metal was a key mediator for ecological risk. Saleem et al. [55] demonstrated that organic matter may retain heavy elements in sediments and play an important role in bottom sediment quality assessment.
The present results showed that there is a strong correlation between Cd, Cr, Cu, Ni, and Pb. Suresh et al. [8] and Wang et al. [7] indicate that the HMs have common geochemical behaviors and originated from similar pollution sources when the correlation coefficient between them is higher. The absence of a correlation among the HMs suggests that the contents of these metals are not controlled by a single factor. Low concentrations of Cd, Cr, Cu, Ni, and Pb indicate their origin mainly from natural sources. In addition, during the research, the impact of the areas immediately adjacent to the reservoirs is marked; moreover, in the case of Ni the number of road and river crossings (RRC) plays an important role. Studies have shown that reservoirs with a shorter water retention time are more likely to supply HMs. In the case of Zn, it was shown that it may come from anthropogenic sources related to wastewater management in the reservoirs' catchment.

Conclusions
The conducted analyses led us to draw the following conclusions:

•
The content of HMs in bottom sediments is an individual feature of reservoirs.

•
The PCA analysis and the values of I geo , EF, MPI and PLI indices show that Cd, Cr, Cu, Ni and Pb in bottom sediments originate from geogenic sources -weathering of rock material. In contrast, Zn comes from anthropogenic sources. • High variability of Zn concentrations between sampling sites and reservoirs confirms their origin from anthropogenic sources related to wastewater management.

•
The PCA analysis indicates that the areas adjacent to reservoirs may have an impact on HM distribution, while the Ni concentration may additionally be affected by road traffic.

•
In reservoirs with higher frequency of water exchange, higher HM concentrations were observed.

•
The highest concentrations of HMs are observed at the inlet to the reservoir and near the dam, which causes the greatest ecological risk.

•
The CA and PCA analysis show that the concentrations and spatial distribution of HMs depend on silt content. The Pb, Cu, Cd and Zn concentrations in reservoirs' bottom sediments are associated with TOM content. However, the relation between individual HMs and silt was stronger than with TOM. • Segregation of sediments along the reservoirs and the TOM content affect the concentrations and spatial distribution of HMs.

Conflicts of Interest:
The authors declare no conflict of interest.