The Effects of Land Use on Concentrations of Nutrients and Selected Metals in Bottom Sediments and the Risk Assessment for Rivers of the Warta River Catchment, Poland

Changes in the environment, aiming at agricultural intensification, progressive urbanisation and other forms of anthropopression, may cause an increase in soil erosion and a resulting increase in the pollution inflow to surface water. At the same time, this results in increased nutrient pollution of bottom sediments. In this study, the concentrations of total nitrogen (TN), total phosphorus (TP), total organic carbon (TOC), calcium (Ca), iron (Fe) and potassium (K) were analysed using bottom sediment samples collected at 39 sites located along the entire length of the Warta River and its tributaries. Agricultural use of land adjacent to rivers was found to significantly degrade sediment quality, while anthropogenic land use (as defined by Corine Land Cover classification—CLC), unlike previous studies, reduces the pollution loads in the bottom sediments. Forest use also contributes to the reduction of the pollution load in sediments. It was found that the significance of the relationship between pollutant concentrations and land use depends on the length of the river–land interface. According to the analyses, the level of correlation between the analysed constituents depends on the use of land adjacent to rivers. The impact of agricultural land use has the strongest effect in the 1 km zone and 5 km in the case of anthropogenic land use. The results showed that the variability of total phosphorus TP concentrations is strongly correlated with the variability of iron concentrations. SPI values indicate that the risk to sediment quality is low due to TOC and Fe concentrations. In contrast, the risk of sediment pollution by TN and TP shows greater differentiation. Although the risk is negligible for 40% of the samples, at the same time, for 33% of the samples, a very high risk of pollution with both TN and TP was found.


Introduction
River sediments are an important part of the cycling of materials in the aquatic ecosystem. The understanding of the pollution processes of river bottom sediments is of great importance because it can be an indicator of the ecological health of waters [1]. Sediments can accumulate pollutants and act as a buffer to absorb and release pollutants into the aquatic environment [2]. Loads of constituents such as carbon, nitrogen and phosphorus in bottom sediments and their interrelationships enable identification of the sources of organic matter and its transformations. The ratios of these constituents are affected by several environmental factors such as climate [3,4], terrestrial inflows, morphometry and use of adjacent areas [5] or the mineralogical composition of sediments [6]. Biochemical and biological transformations taking place in sediments and at the water-sediment interface [7,8], as well as the hydrodynamic processes of sediment transport in riverbeds [9], are also of great importance for concentrations and mutual proportions of constituents that can be found in bottom sediments.
The land-use pattern of the area located in the catchment is one of the more important factors affecting the concentrations of nutrients, as well as other pollutants, in river bottom

Site Description
The catchment area of the Warta River is 54,519 km 2 . In terms of length, 808 km, it is the third-longest river in Poland ( Figure 1). The tributaries of the Warta River, which were also analysed in the presented study, are the Prosna River; the Bobrowski Canal; the Mosiński Canal; the Wełna, Obra and Noteć rivers as well as the additional tributaries of the Noteć River: the Gwda and Drawa rivers ( Figure 1). The source of the Warta River is 380 m above sea level, and its mouth to the Oder is 12 m above sea level. The catchment area of the Warta River varies in terrain. The upper-southern part of the catchment is located in the Kraków-Częstochowa Upland, passing through the Central Polish lowlands to the northern lake districts. The tributaries of the Warta River range from 11 km (the Bobrowski Canal) to 391 km (the Noteć River). Their areas range from 28 km 2 in the case of the Bobrowski Canal to 17,319 km 2 in the case of the Noteć River (Table 1). Almost along the entire length of the Warta river as well as its tributaries drains a post-glacial landscape formed during the retreat of the Vistulian ice-sheet. The layer of glacial sediments with a thickness of more than 200 m consists of tills, glacial sands and gravels and outwash sands and gravels [35]. In the Warta River catchment area, sandy soils make up 45% of the land area, clay soils 41% and organic soils and alluvial soils 14% of the land area [36]. According to the European Soil Database v2.0 m, the main soils are Podzol-32%, Luvisols-30%, Fluvisols-12% and Arenosols-12%.
According to the CLC classification, in the whole catchment area of the Warta River, agricultural land use makes up 60.2% of the catchment area, forest and semi-natural land use make up 32.5%, anthropogenic land use makes up 5.6%, surface water makes up 1.4% and the rest is covered by wetlands and peat moors (Table 1, Figure 2).  The highest proportion of agricultural land has the catchment area of the Mosiński Canal-approximately 78%. This catchment also has the lowest proportion of forest land and semi-natural land in a total catchment area-less than 16%. A slightly lower, compared to the aforementioned catchment, proportion of agricultural land in the total catchment area, 72%, is found in the catchment areas of the Prosna and Wełna rivers. On the other hand, the catchment area with the highest proportion of forest land and semi-natural forest land is the Drawa River catchment located in the northern part of the Warta catchment, where this proportion is 61.8%. This catchment also has the highest proportion of the open water area-3.9% of the total catchment area. The proportion of wetlands in all catchments is very small, reaching a maximum of 0.3% for the Mosiński Canal catchment.
An important source of bottom sediment contamination is the material brought into rivers as a result of erosion processes. One of the indicators that identify areas sensitive to erosion is the Topographic Wetness Index TWI. The TWI takes into account the upslope area and its slope, allowing the calculation of the steady-state wetness and runoff across the analysed area. Calculated TWIs show values in the range of 0. 9-14.4 and are shown in Figure 3. Geochemical conditions in the Warta River catchment, which describe the phosphorus, calcium and iron concentrations found in a 20 cm topsoil, are shown in Figure 4. The Warta River catchment has the lowest annual precipitation in Poland [37]. Climatologically, according to the Köppen-Geiger's classification, the catchment area is located in the oceanic climate (Cfb) that has warm summers and milder winters [38,39]. The average annual air temperature for the WRC is approximately 8 • C, while the annual precipitation ranges from 520 mm in the north-eastern part of the catchment to 675 mm in the southern part [40]. The mean annual runoff is 3.86 dm 3 ·s −1 ·km −2 for the Warta catchment and it shows considerable spatial variability [41]. The lowest runoff values, below 2 dm 3 ·s −1 ·km −2 , are observed in the upper catchment of the Noteć River [42].

Materials
The analyses used data concerning the measurements of N, P, K, Ca, Fe and TOC concentrations in bottom sediments of the Warta River and its tributaries, which were obtained as part of the State Monitoring of the Natural Environment programme. Bottom sediment samples were collected from 39 sites in 2016 ( Figure 1). Sampling sites were located at the borders of the catchment area, at the mouths of tributaries and below cities where large industrial plants are located. For chemical analyses, a 5 cm top layer of sediments was collected from 4-5 sites over a 50 m distance for each location using a van Veen grab sampler. Samples for each site were mixed and rubbed through a nylon sieve with 2 mm mesh. Laboratory measurements were performed for Ca, K, Fe and P using inductively coupled plasma atomic emission spectrometry (ICP-AES) according to PN The boundaries of the catchment area of the Warta River, its course and the course of its analysed tributaries were determined based on a Computer Map of the Hydrographic Division of Poland, which contains full hydrographic data of Poland in vector format. The land-use structure in the catchment area was made based on the CLC database provided by the Chief Inspectorate for Environmental Protection. The landform was described using a Digital Elevation Model (DEM) on a grid of at least 100 m, which was obtained from the Central Office of Geodesy and Cartography. The Ca, Fe and P concentrations found in the 20 cm topsoil were estimated using the Geochemical Atlas of Poland at a scale of 1:500,000, which was obtained from PGI-PIB (Polish Geological Institute-National Research Institute).

Overall Statistics
To evaluate the variability of N, P, K, Ca, Fe and TOC concentrations in bottom sediments of the analysed rivers, the mean values of the concentrations, the median and their maximum and minimum values were calculated. Concentration distributions were determined using the Shapiro-Wilk test. Outliers were determined using the Interquartile Range (IQR) test. The values that fall below or above 1.5-IQR were considered as outliers of the variables: where Up-concentrations for values that fall above the 75th quartile (Q 3 ) and Lowconcentrations for values that fall below the 25th quartile (Q 1 ). Correlations were determined using Spearman's rank coefficients. All statistical calculations and visualisation of results were performed using R version 4.0.5.

Risk Assessment Method
An ecological risk assessment of nutrients in river bottom sediments was performed using a method based on the Single Pollution Index (SPI) [43].
The SPI i is calculated based on the formula: where C i is a measured concentration of the evaluated i factor and C s is a standard concentration of the evaluated i factor. Calculations were performed for TN, TP, TOC and Fe. Values of 550 mg-kg −1 , 600 mg-kg −1 , 1% and 2% were used as standard values for TN, TP, TOC and Fe, respectively, which are derived from safe nutrient concentration limits that can be found in the Sediment Quality Guidelines [25,43]. Based on the P i , four hazard classes can be distinguished as shown in Table 2 [43].
Seriously polluted C:P:N stoichiometry enables the description of potential sources of pollution as well as geochemical and biochemical processes occurring in the environment [31,44]. TOC/TN is one of more frequently used ratios, which identifies potential sources of organic matter [43]. TOC/TP reflects to some extent the rate of conversion of organic carbon and phosphorus compounds [45]. In contrast, TN/TP reflects the dynamics of accumulation, deposition and release of nitrogen and phosphorus in water [46,47]. TOC/TN > 10 indicates primarily a terrestrial source of organic matter. In contrast, TOC/TN < 10 indicates an aquatic source of organic matter. TOC/TN~10 indicates that the organic matter found in sediments is of both terrestrial and aquatic origin [48,49].

Analyses of Spatial Variability of Concentrations of Selected Elements Found in Bottom Sediments
The principal component analysis (PCA) and cluster analysis (CA) were used for the differentiation of groups of sampling sites exhibiting similar variability and similar sources. Due to concentrations of analysed substances in bottom sediments, the cluster analysis of sampling sites was performed using k-medoids [50,51]. K-medoids is a technique that is less sensitive to the effect of noise and outliers in analysed data, compared to the k-means. The number of cluster groups was found using the silhouette method [52].
Land cover was defined for river-adjacent zones of a width of 500 m, in % of the total area of the zone. The adopted zone lengths were 1, 2, 5 and 10 km from the sampling sites ( Figure 5). In zones defined as above, GBI's content index for P, Ca and Fe was also calculated based on the formula: where GBC i means the topsoil content class for each element, while A i is the area per class in an analysed zone. Spatial variability data concerning the parameters were determined using QGIS 3.18 whereas statistical calculations were performed using R 4.0.

General Characteristics of Concentration Variability
The results of measurements made in 2016 for TN, TP, Ca, Fe and TOC concentrations found in the bottom sediments of the Warta River and its tributaries are shown in Table A1 The average concentrations of the analysed elements ranged from 13,182 mg-kg −1 for Fe through 9320 mg-kg −1 for Ca, 2169 mg-kg −1 for N to 1130 mg-kg −1 for P and 3.76% for TOC ( Table 3). The median values are 8269 mg-kg −1 , 11,078 mg-kg −1 , 998 mg-kg −1 , 464 mg-kg −1 and 2.18% for Fe, Ca, N, P and TOC, respectively. The TOC values measured in 1998-2000 are similar to those measured in bottom sediments of the Oder, which is a receiving body for the water of the Warta River, when they ranged from 0.2% d.m. to 18.2% d.m., averaging 4.9% d.m [53]. Similar values for P, Ca, Fe and TOC concentrations were observed by House and Denison [31,54] in seven rivers in southern England. Their reported data concerning concentrations were 40-47,000 mg-kg −1 for Fe, 40-19,000 mg-kg −1 for Ca, 10-2000 mg-kg −1 for P and 0.6-19% d.m. for TOC. This may indicate similar sources of sediment pollution and similar land use of catchments. The IQR values shown in Table 3 were 2685 mg-kg −1 , 1084 mg-kg −1 , 9934 mg-kg −1 , 13,821 mg-kg −1 and 3.18% for TN, TP, Ca, Fe and TOC, respectively. The IQR test showed several outliers involving N concentrations in the Warta River bottom sediments at site no. 4 and the Noteć River sediments at site no. 82. For phosphorus, concentration outliers were found in the Warta River sediments at sites no. 4, 19, 20 and in the Noteć River sediments, at site no. 82. For Ca, outliers were observed in the Mosiński Canal at sites no. 51 and 53. For Fe, outliers were found in the Noteć River at site no. 82. Outlier TOC values were observed in the Warta River sediments at site no. 4 and in the Noteć River bottom sediments at sites no. 82, 83 and 84. The outliers were most frequently observed at site no. 82 (Noteć River)-four times and three times higher/lower at sites no. 4 and 23 (Warta River).
The analysis of distributions of the analysed concentrations showed that none of the concentration distributions is a normal distribution ( Figure 6). In contrast, for all p > 0.05 parameters, all distributions are similar to a log-normal distribution.
All analysed concentrations are strongly positively correlated with one another. Spearman's correlation coefficients range from 0.50 for the correlation between TOC and Ca to 0.93 for the correlation between TP and Fe. Almost all correlation coefficients are significant at the p < 0.0001 level (Table 4). Only the TOC-Ca correlation is significant at the p < 0.01 level, whereas the TOC-K correlation is significant at the p < 0.001 level. The strongest correlation is between Fe concentrations and TP, Ca and K concentrations, for which the correlation coefficients are 0.93, 0.87 and 0.86, respectively. In contrast, the correlations between TOC and Ca/K are the least correlated, but still significant at the p < 0.01 level.  Inorganic constituents that can combine with or adsorb phosphates are one of the factors affecting TP concentrations in bottom sediments [55,56]. Immobilising components include metals, such as iron, which combine with phosphorus, forming crystalline metal phosphates (e.g., strengite or vivianite), or they can adsorb phosphorus on the oxide/hydroxide layer. In the case of the presence of large amounts of calcium, phosphate precipitation may occur in the form of hydroxylapatite (HA) or various calcium phosphates [55,57,58]. Furthermore, iron (III) plays an important role in the phosphate sorption by humic elements, forming iron (III)-humus complexes with phosphates [55,59].

C:N:P Stoichiometry
The calculated characteristic values of C:N, C:P and N:P ratios indicate their high variability as shown in Table 5 and Figure 7. C:N values range from 4.51 to nearly 112, with a median value of 22.98. The highest C:N values of 112, 81, 68 and 61 were observed at sites no. 11, 52, 3 and 21, respectively. A total of 33 sites had C:N > 10, two sites had ratios similar to 10, and only four sites had C:N < 10 ( Figure 7). This indicates that the organic matter found in bottom sediments derives mainly from runoff from surrounding areas. Similarly, high C:N values were also obtained by Lu [60], indicating that exogenous sources of pollution are the cause. On the other hand, the lowest C:N values have sediments at sites located in the vicinity of the aforementioned site no. 52, i.e., sites no. 51 and 53-4.5 mg-kg −1 and 6.9 mg-kg −1 , respectively. The bottom sediments at these two sites have high nitrogen concentrations as well as very high P, Ca and Fe concentrations (Table A1).  N:P ratios in the bottom sediments of rivers of the Warta River catchment ranged from 0.28 to 15.37, with a mean value of 2.32. According to Chen et al. [45], low N:P values correspond to eutrophic and mesotrophic conditions with the predominant supply of nutrients derived from external sources. The nutrients from such sources are complex and have low N:P ratios. In contrast, high TN/TP values can be observed under oligotrophic conditions when natural nutrient supply sources have high N:P ratios. Knösche [55] reports that the TN:TP ratio in river sediments decreases significantly from low-density Paleopotamal sediments (median TN:TP = 23.7) to high-density Eupotamal sediments (median 8.5). He also observed a similar correlation for the TOC:TP ratio. At the same time, Ostendorp [61] indicated that N:P:OC ratios could be explained by the variability of organic content rather than by eutrophication processes in each case.

Single Pollution Index (SPI)
The pollution status of bottom sediments of the rivers of the Warta River catchment shows great spatial variation, as shown in Figure 8. SPI values for nitrogen range from 0.13 to 29.5. The largest proportion of sediment samples, 41%, can be classified as Class I (Table 6). Class II includes 13% of the samples. Sampling sites with a low risk of nitrogen pollution are located mainly in the upper part of the Warta River (sites no. 1-14, excluding sites no. 4 and 12) and in its lower reaches (sites no. 18-24, excluding site no. 20). On the other hand, Class III includes 13% of sediment samples and Class IV, with the highest pollution, includes as many as 33% of samples. These sites are located primarily in the middle course of the Warta River and in its tributaries.   I  II  III  IV   TN  41  13  13  33  TP  39  15  10  36  TOC  23  54  15  8  Fe  56  23  18  3 SPI values for sediments polluted with phosphorus ranged from 0.08 to 14.8. Classification of sediment samples classifies 39% of samples into Class I, 15% into Class II, 10% into Class III and 36% into Class IV. These values are similar to SPI for TN. As with nitrogen, the sampling sites classified as Class I and II are located in the upper and lower reaches of the Warta River. The SPI values for TOC allow 23% of sediment samples to be classified as Class I and 54% as Class II. On the other hand, 15% and 8% of sediment samples should be classified into Class III and IV, respectively. This indicates a much lower risk of bottom sediment pollution by TOC than by TN or TP. The lowest pollution risk to bottom sediments was found for Fe. SPI values ranging from 0.08 to 3.4 enabled 56%, 23%, 18% and 3% of the samples to be classified as classes II-V, respectively.
The obtained SPI values indicate a lower pollution risk to bottom sediments of the Warta River than those obtained by He et al. [43] for Chinese rivers.

Spatial Variation of Pollutant Concentrations
The clustering analysis of TN, TP, TOC, Ca, Fe and K concentrations in bottom sediments using the Partitioning Around Medoids (PAM) identified three groups of bottom sediment sampling sites (Figure 9). A simultaneous principal component analysis (PCA) found that the first two principal components, for which the eigenvalues are greater than 1, explain approximately 87% of the total variance (Table 7).  The land-use structure, as determined by the CLC classification, shows varying effects on the levels of constituents found in river sediments of the Warta River catchment. Anthropogenic lands (C_1) and forests and semi-natural ecosystems (C_3) show a negative correlation, while agricultural lands (C_2) and lands covered by water (C_5) show a positive correlation with concentrations of all elements; however, this relationship is not statistically significant in most cases. It should be noted that previous studies highlighted the negative impact of urbanised areas on the quality of bottom sediments in rivers. This was mainly due to the impact of direct discharge of wastewater into rivers or increasing surface runoff from areas with impervious surfaces [62,63]. However, the following phenomena that have been marked in recent years, such as the connection of most wastewater sources to wastewater treatment plants, the increase of areas that enable rainwater infiltration or, for example, measures to retain rainwater runoff from impervious surfaces, may reduce or neutralise the negative impact of urbanisation on the quality of surface water and bottom sediments. The negative impact of agricultural areas is mainly caused by an increased input of nutrients [12], inorganic suspended solids [64] and organic matter [65] from these areas to surface waters.
Factor loadings for the first two principal components are shown in Table 8. According to the classification presented by Liu et al. [66], factor loadings were classified as "strong," "average," and "weak" for values of >0.75, 0.75-0.50, and 0.50-0.30, respectively. PC1, which explains 67.4% of the total variance, has strong positive Fe, TN, TP and TOC loadings, while also showing average K and Ca loadings.  Table 9 shows the Spearman's correlation coefficient values of TN/TP/K/Ca/Fe/TOC concentrations in bottom sediments with land cover structure in the zones of a width of 1000 m, adjacent to the river on a length of 1, 2, 5 and 10 km, and the calculated Ca, Fe and P content index of the topsoil.
The significance of the land-use impact depends not only on the land-cover structure itself, but also on the length of the contact zone of land adjacent to the watercourse. A significant (p < 0.05) effect of agricultural land (C_2) on the concentration levels of all analysed elements was observed for the 1 km buffer area. For longer buffer lengths, excluding K and TOC concentrations in the 2 km buffer area, that effect was no longer statistically significant. In contrast, the proportion of anthropogenic land (C_1) is significantly linked to P, K, Fe and TOC concentrations in the 5 km buffer area. Additionally, the proportion of the area covered by water (C_5) shows the strongest relationships with P, K, Ca and TOC concentrations in the 5 km buffer area. At the same time, the proportion of the land covered by water is the only one that is significantly correlated with phosphorus concentrations for all buffer lengths. Except for the effects of forest land proportion (C_3) on Fe concentrations in the 10 km buffer area, there was no significant impact of that land-use category on concentrations of analysed substances found in bottom sediments.
The results indicate a significant positive relationship between the geochemical background level of phosphorus in the 2 km buffer area and N, P, Ca, Fe and TOC concentrations in the bottom sediments analysed. In contrast, the geochemical background levels of Ca and Fe do not show such a relationship.

The Relationships between Water Quality Parameters and Land Use
The analysis of land use as identified based on the CLC data in individual buffers enabled four land-cover groups to be distinguished. The clustering results and the distribution of each group are shown in Figure 10. The first group includes areas dominated by anthropogenic lands, with a low proportion of forests. The second group includes buffer areas that are predominantly forested. The third group includes buffer areas dominated by agricultural lands. The fourth group consists of buffer areas where the proportion of lands covered by water is at least 10%, regardless of other land-use categories.  The Spearman's correlation coefficients between parameters describing constituent concentrations in bottom sediments for defined land-use groups are shown in Table A2. In group 1, where anthropogenic land use is the main factor, the strongest relationship is between TP and Fe, TN and TOC concentrations. However, the relationships between the elements show a higher correlation in the case of buffers that were identified based on a significant proportion of forest land. The highest correlation is between TP and TN, TN and TOC and TP and Fe. Unlike cluster 1, there is also a strong correlation between TOC and Fe in cluster 2.
The highest levels of correlation between bottom sediment constituents occur when the river is adjacent to agricultural land (cluster 3). In this case, all correlation coefficients are statistically significant, at least at the p < 0.05 level. The highest correlation in cluster 3 is between Ca and K, TP and Fe and Fe and K. The lands included in cluster 4 have surface water within the buffer area, the proportion of which is between 10% and 20% (highest proportion) of the total buffer area. In the case of sediment samples collected at sites assigned to this cluster, the highest correlation occurred between Ca and K concentrations, where r = 1.0. Again, very high levels of correlation occurred between Fe and TP and Ca and K.
It should be noted that regardless of the influence of the environment, phosphorus concentrations show the highest correlation with iron concentrations and slightly lower correlation with calcium concentrations, which is due to the already mentioned mutual interactions between these constituents, leading to precipitation or release of phosphorus [67][68][69]. The strong correlations between TP and TN concentrations for all use categories of the river-adjacent zones, except for the relatively low correlation for agricultural use, may indicate similar sources of nutrients in these areas. In agricultural lands, however, sources of pollution are of different nature-diffuse or point. Nutrients may derive from both field crops and animal husbandry, and they have different nutrient ratios, as also indicated by Lee [70] and Jones [10]. At the same time, bottom sediments collected from sites included in cluster 3 have the highest average TN, TOC, Ca and Fe concentrations and also high TP concentrations (Table 10). In contrast, the concentration values for clusters 1 and 2 (i.e., river-adjacent zones with anthropogenic and forest land uses) are significantly lower. Average sediment constituent concentrations for cluster 3 compared to concentrations for clusters 1 and 2 are 4.3-5.8 times greater for nitrogen, 2.3-4.0 times greater for phosphorus, 3.4-4.2 times greater for TOC, 4.3-5.8 times greater for calcium and 1.9-2.8 times greater for iron. The intermediate values of constituent concentrations were found in the sediments for cluster 3 determined on the basis of the proportion of surface water in the river-adjacent zone. As this proportion is between 10-20%, the remaining area of the land may have very different uses and, consequently, average out the effects of different land uses on the bottom sediment composition.

Conclusions
Based on the results of analyses concerning spatial variability of N, P, TOC, Ca, Fe and K concentrations in the bottom sediments of the Wart River and its tributaries, the following conclusions can be made:

•
The dominance of agricultural activities in the zone of land adjacent to the river was found to be a significant factor increasing the pollutant content found in the bottom sediments of the analysed rivers. On the other hand, artificial and wooded areas reduce the concentrations of analysed substances in bottom sediments. The significance of this impact also depends on the length of the contact zone between the river and the land adjacent to it.

•
The results of C:N:P ratios indicate that the organic matter found in sediments has mainly a terrestrial source. These results also indicate that there is a threat of eutrophication to a substantial number of the rivers.

•
Based on the analysis of SPI values, it can be concluded that the risk resulting from TOC and Fe is relatively low or absent in almost 80% of the sampling sites, and it is medium to high in the remaining sampling sites. On the other hand, SPI values for TN and TP indicate that although the risk is absent for 40% of the samples, at the same time, more than a third of the samples are severely polluted. • Future analyses will require consideration of the temporal variability of pollutant concentrations. Numerous studies indicate that sediment parameters vary not only over annual periods but also over multi-year periods [32]. This is due to the influence of river processes on sediment composition, in addition to the inflow of substances from adjacent areas, due to different water flow rates and the resulting migration of pollutants.

Funding:
The publication was co-financed within the framework of the Ministry of Science and Higher Education programme as "Regional Initiative Excellence" in years 2019-2022, Project No. 005/RID/2018/19.

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

Conflicts of Interest:
The author declares no conflict of interest.