Floodplain soils contamination assessment using the sequential extraction method of heavy metals from past mining activities

Floodplains are among the most precious and threatened ecosystems in the world. The study deals with floodplain soil contamination caused by 8 heavy metals (HMs) (Cd, Co, Cr, Cu, Mo, Ni, Pb, Zn) originating and transported from old mine works along the Štiavnica River in Slovakia. We determined the total HMs content and the HM fractions using BCR sequential extraction method. We selected 12 alluvial sites (AS), two contaminated sites (CS), and one reference site (RS). The sampling points were located within the riparian zones (RZ), arable lands (AL), and grasslands (GL). We confirmed soil contamination by HMs and the related ecological risk by different factors. The contamination by HMs at many AS localities was similar or even higher than at CS localities. The highest contamination factor was calculated for Cu (39.8), followed by Pb (27.4), Zn (18.2), and Cd (7.2). The HMs partitioning in the different fractions at the CS and AS localities revealed that Cd, Zn, and Pb were mainly associated with the exchangeable and reducible fractions, while Cu was mainly associated with the oxidisable fraction. The soil properties were selectively correlated with the HM fractions. Based on the ANOVA results, the effect of different ecosystem types on HM fractions was revealed.

Floodplains are relatively flat land surfaces created when alluvial material carried by surface water is deposited. Rivers and floodplains are connected, exchanging water, sediments, biota, and nutrients in a shared natural ecosystem. In the natural state, floodplains are characterised by high biodiversity with the presence of riparian areas, natural or artificial grasslands, and agricultural land. Floodplains fulfil many ecosystem services, including regulating such as filtration. They cover 7% of the total European area and host 12% of the European population. In Slovakia, more than 25% of the population lives in floodplain areas 1 . Floodplains are among the most dynamic, productive, diverse, and precious ecosystems in the world 2,3 . On the other hand, floodplains are the most threatened ecosystems in the world. In degraded floodplains, the quality of ecosystem services is reduced. In Europe and North America, up to 90% of floodplains are cultivated and, therefore, functionally extinct. In countries with fast growing economies, there is a huge discharge of pollutants into floodplains that are so vulnerable and highly threatened 4 .
Today, many river floodplains are highly contaminated by HMs that exceed standards for soils and sediments and also the filtering capacity of alluvial soils. Pollution is a very prominent and lasting environmental problem around the world, despite the efforts made in recent decades to clean up the environment. Soil contamination in Europe is a great challenge and puts human and environmental health at risk. Soil contamination is mainly located near waste landfills and industrial activities that spread heavy metals (HMs) and other potentially toxic elements. Important anthropogenic sources of HMs are present in former mining activities, foundries, smelters, and other diffuse sources 5 . Historically, the mining industry is one of the oldest documented types of human activity. Mining waste rocks and mining tailings contain chemicals and residues characterised by a high content of HMs. Mining activities as the main source of HMs constitute a real concern around the world, and many researchers study their impact on soil pollution [6][7][8][9] . HMs contained in mining waste can be released from the mining field into the environment and can be dispersed by gravitational transport and fluvial pathways. Floods are directly involved as serious agents of contaminant dispersion 10,11 , resulting in sedimentation on lands, where contaminants can remain for decades and centuries until they are remobilised by surface or river bank erosion 12 .

Material and methods
Sampling sites and soil sampling. Sampling sites were selected in the floodplain along the Štiavnica River in central Slovakia. The riverbanks consist of fluvial deposits (sand and silt). The Štiavnica River is 55 km long, with a 443 km 2 large catchment area. In the upstream part of the basin in the Štiavnické vrchy Mounatins of volcanogenic origin, the relief is dominated by hills and valleys with forests and partly by permanent grasslands used as meadows. The bedrock consists mainly of andesites and rhyolites with scattered occurrences of conglomerates and shales. In the downstream part of the basin in the Podunajská pahorkatina Hills, the relief is dominated by a flat surface. The bedrock consists of clay-and sand-stones on slopes, sand and gravel in areas of floodplains. The land is used as arable land. Due to volcanogenic origin, the region of the Štiavnické vrchy Mountains is mineral-rich. The extraction of precious metals in the central part of the Štiavnické vrchy Mountains has a very long history. Large-scale mining started in the twelfth century. The discovery of precious metals encouraged rapid expansion. Silver ore prevailed among the mined metals. In the eighteenth century, the region became the biggest mining centre in the Habsburg Monarchy. In the period between 1790 and 1863, the mountains yielded 490 metric tons of silver and 11 tons of gold. At that time, there was no environmental legislation preventing the release of contaminated water and sediments into watercourses. At the end of the nineteenth century, mining fell in decay. This trend also continued in the twentieth century. Past mining activities left imprint on the environment, including contamination by HMs and other potentially toxic elements. Nowadays, the storage and handling of tailings is a major environmental issue in the region. There exist several old mine tailings, waste heaps, and dumps. The ore field is drained by the Štiavnica River.
To study the contamination of floodplain soils by HMs from past mining activities and transported by the Štiavnica River and its sediments, we selected a total of 15 study sites along the Štiavnica River, including 2 contaminated sites (2CS, 4CS), 1 reference site (7RS) and 8 alluvial sites (AS) located in the Štiavnické vrchy Mountains (1AS, 3AS, 5AS, 6AS, 8-11 AS), and 4 alluvial sites situated in the Podunajská pahorkatina Hills (12-15 AS) (Fig. 1).
The soil types on the floodplain were classified as Haplic Fluvisols according to the World Reference Base for Soil Resources. The sampling points were located on the river bank within the riparian zones (RZ) in the Štiavnické vrchy Mountains, and on arable land (AL) in the Podunajská pahorkatina Hills. In the case of a larger floodplain, we also collected samples distant from the water flow at the end of the floodplain, that were used as permanent grasslands (GL) in the Štiavnické vrchy Mountains, and as arable land in the Podunajská pahorkatina Hills (at 5 of 12 alluvial sites). Outside the river floodplain, we selected the two contaminated sites (CS) classified as Ekranic Technosols, which are poorly developed on metallurgical slag and shallow soils. As CS sites, we selected the area of the Maximilian shaft that is registered in the Register of environmental burdens of the Slovak republic as probable environmental burden (PEB), and the hard rock mine tailing Lintich that is registered as a confirmed environmental burden (CEB), where hazardous substances caused by mining activities can pose a risk to human health or to the environment, soil and groundwater. The Maximilian shaft area was located approximately 1 km and the Lintich mine tailing approximately 100 m from the Štiavnica River. The RS site represents a site with a ´normal´ level of contaminants, and arising from a combination of natural and diffuse pollution contributions. The RS site was used to determine the background values. The RS site was located www.nature.com/scientificreports/ 0.1 M Na 4 P 2 O 7 ·10H 2 O. Oxalate-extractable Al, Fe and Mn as an indicator of hydrous oxides were determined by the acid oxalate extraction method described by Pansu and Gautheyrou 20 with a (NH 4 ) 2 -oxalate solution adjusted to pH 3, shake for 4 h in the dark and centrifuge for 10 min at 10,000×g. The hydrological soil properties were determined using overall undisturbed soil samples using core-extracting tubes (Eijkelkamp Equipment for Soil Research, The Netherlands). Undisturbed soil samples in tubes were taken to the laboratory for hydrophysical analyses by the gravimetric method of particle density (PD), bulk density (BD), maximum capillary water capacity (MCWC), field water capacity (FWC), actual soil moisture (SM), porosity (as total-Pt, capillary-Pc, non-capillary-Pn, semi-capillary-Ps), and minimal air capacity (MAC) according to Novak. Particle size analysis was performed by pipette method using soil particle sedimentation and particle size fractions (sand, silt, clay) were classified according to the United States Department of Agriculture (USDA) system. Data analysis. The degree of the soil pollution and the related ecological risk we described with the help of contamination factor (CF), pollution load index (PLI), potential ecological risk of individual metals (Er), and potential ecological risk of the environment (RI).
The CF also designated as an enrichment factor was calculated according to the formula: where Ci is the measured content of HM in soil, and Cb is the background value of the HM contents of the reference site. According to Varol 21 , the CF < 1 means low contamination, 1-3 is moderate contamination, 3-6 is considerable contamination, and CF > 6 means very high contamination.
To identify heavy metal contamination as a total scale in soil, Tomlinson et al. 22 recommended an equation for calculating the pollution load index (PLI): where n is the number of metals. PLI > 1 indicates pollution and demonstrates dynamic deterioration of quality.
According to Hakanson's 23 methodology, the potential ecological risk (Er) of a given heavy metal is defined as follows: where Tr is the toxic response factor for a given substance, CF is the contamination factor. The Tr was suggested by Hakanson 23 for six metals Cd (30), Pb (5), Hg (40), Cu (5), Ni (5), Zn (1). The sum of individual potential risks (Er) is the potential ecological risk of the environment (RI) of the study site: The Er < 40 means low risk, 40-80 is moderate risk, 80-160 is considerable risk, 160-320 is high risk, and > 320 is very high risk. The RI < 150 means low risk, 150-300 is moderate risk, 300-600 is considerable risk, and > 600 is very high risk.
To identify the relationships between HMs content and soil parameters, statistical analyzes were conducted using the Pearson correlation coefficient. Hierarchical cluster analysis (HCA) was applied to group similar localities (spatial variability) using the Ward´s method and Euclidean distances. A two-way ANOVA analysis was performed to assess statistically significant differences in HM fractions between different types of ecosystems (riparian zone, alluvial grassland, arable land), and soil depths (0-10, 20-30 cm). The SPSS Statistics 19 software was used to perform statistical analysis of the data.

Results
Basic characteristics of soils. The basic statistical characteristics of the chemical and physical properties of the soils were tested in the laboratory and calculated at the CS, AS, and RS sites in two depths. Data are reported in Tables 1, 2. There are differences between the observed properties at specific sites (CS, AS, RS) and in different depths, but also between AS sites themselves, confirming their heterogeneity.
Soil contamination by heavy metals. The total HM content in the two soil depths of the Štiavnica River floodplain together with the Slovak, Finish and Canadian guideline values of the contaminated soil are presented in Table 3. The Finish legislation sets content levels for each hazardous element to identify soil contamination and remediation needs. It distinguishes the so-called ´threshold value´ and ´guideline value´. If the guideline value is exceeded, the area has a contamination level that presents ecological or health risks.
For all HMs, except Cr and Ni, we recorded at least one locality exceeding one of the limit values. In the case of Mo, the total content exceeded the limit value only at one CS locality. The total Co content slightly exceeded only the Slovak limit value at CS localities and several AS localities, and was below other national limit values. We recorded 4 HMs (Cd, Cu, Pb, Zn) whose total content in most localities significantly exceeded all limit values, including the Finnish guideline values for soils with decontamination and remediation needs. At AS localities, the situation in total HM content of both observed soil depths was similar. HM contamination was also observed across broader floodplains. www.nature.com/scientificreports/ Contamination and ecological risk assessment for heavy metals in soil. Contamination and ecological risk at the CS and AS localities were confirmed by calculated indices along the entire Štiavnica River floodplain including the downstream part of the Podunajská pahorkatina Hills (Table 4). In contaminated sites, the average CF showed a very high contamination by Zn and Cu; considerable contamination by Cd, Mo, Pb; and moderate contamination by Co, Cr, and Ni. At alluvial sites, the average CF showed very high contamination by Cu, Pb, and Zn; considerable contamination by Cd; moderate contamination by Cr, Mo, and Ni; and low contamination by Co. Very high contamination (CF > 6) was calculated for Pb, Zn, and Cu at most AS localities. Overall, the highest CF was calculated in the case of Cu (39.8 in the 15AS-AL locality), Pb (27.4 in the 1AS-RZ locality), Zn (18.2 in the 6AS-GL locality), Cd (7.2 in the 1AS-GL locality). Contamination by HMs was, in many AS localities, similar or even higher than in CS localities. Metal PLI values confirmed that the CS and AS sites were contaminated except for five AS localities, two of them are used as grasslands, two as arable land, and one is located in the riparian zone. The average Er showed a considerable risk of Cd at the AS and CS localities, and moderate risk of Cu and Pb at the AS localities. The RI values indicated a considerable ecological and environmental risk at 7 study sites located mainly in riparian forests, moderate risk at 11 study sites and only 2 study sites located in the riparian zone were at low risk.
Heavy metals fractionation. The sequential extraction procedure enables us to measure broader forms or phases.  Table 2. Soil physical properties in the soil depth of 0-10 and 20-30 cm in contaminated, alluvial and reference sites with basic statistical characteristics. AS alluvial site, CS contaminated site, RS reference site.

Depth (cm) Site
Content of size particles %
The content of HMs in the residual fraction was dominated at the RS localities (in 7 of 8 HMs). The study of the partitioning of HMs in the different fractions at the CS and AS localities revealed that Cd, Cu, Pb, and Zn were mainly associated with the first three fractions. The first exchangeable fraction was the highest in the case of Cd and Zn. Cu was mainly associated with the oxidisable fraction bound to organic matter and sulphides. The distribution of HMs in the different fractions, mainly in exchangeable ones, suggested that the mobility and bioavailability of the metals probably declined in the order of Cd, Zn, Cu, and Pb at AS localities, and Zn, Cd, Pb, and Cu at CS localities.
The  Table 4. Contamination factor (CF), pollution load index (PLI), potential ecological risk factor (Er), potential ecological risk of the environment (RI) for contaminated and alluvial sites in the soil depth of 0-10 cm. AS alluvial site, CS contaminated site, RZ riparian zone, AL arable land, GL grassland, PEB probable environmental burden, CEB confirmed environmental burden. The CF < 1 means low contamination, 1-3 is moderate contamination, 3-6 is considerable contamination and CF > 6 means very high contamination. The PLI > 1 indicates pollution. The Er < 40 means low risk, 40-80 is moderate risk, 80-160 is considerable risk, 160-320 is high risk, and > 320 is very high risk. The RI < 150 means low risk, 150-300 is moderate risk, 300-600 is considerable risk, and > 600 is very high risk.  The selected HM fractions of all HMs were significantly positively correlated with pH except for a significant negative correlation with the I Cu fraction. No statistically significant correlations were found between all HM fractions and Eh. The most numerous, mainly positive, correlations were found with the oxalate-extractable Fe and Mn in the case of Mo followed by Cd, Co, and Cr, while they did not occur in the case of Pb. Correlations with SOC occurred only with Cu (I and II fractions). Significantly positive correlations occurred between HA/ FA in the case of selected fractions of Cr, Mo, Ni, Pb, Zn and significantly negative only in the case of Cu (I and II fractions) ( Table 5).
The most numerous correlations were found with the different particle size fractions in the case of Ni followed by Cr, Cu, Mo, and Co, while they did not occur in the case of Pb. While the selected HM fractions were mostly positively and significantly correlated with the content of clay particles, a reversed situation occurred in the case of the content of sand particles. No statistically significant correlations were found between all HM fractions and PD and only one negative correlation was found between BD and the IV Cd fraction. The water capacity parameters (MCWC and FWC), only negatively correlated with the selected HM fractions of Co, Cr, Cu, Ni, and Pb. The Cd, Mo, and Zn fractions were not correlated with MCWC and FWC. The most numerous correlations were found with the different porosity categories and MAC in the case of selected fractions of Cu followed by Co, Ni, and Cd (Table 6).
In Fig. 3, the hierarchical dendrogram of the linkage between the localities is shown for the surface soil layer (depth 0-10 cm). The variables taken for this analysis are the same as the Pearson´s correlation analysis.
The HCA has highlighted three large clusters, A, B, and C. Cluster A includes sites located in the lower part of the river located in the Podunajská pahorkatina Hills in the broader floodplains predominantly used as arable land. Cluster B includes two contaminated sites located in the Štiavnické vrchy Mountains and one site located in the lower part of the river situated in the Podunajská pahorkatina Hills in the broader floodplains used as arable land. The other localities are grouped in cluster C. Two clusters of lower hierarchical order can be recognised in  www.nature.com/scientificreports/ C (C1 and C2). Cluster C1 includes three sites located in the middle part of the river located in the Štiavnické vrchy Mountains, of which one site is located within the riparian zone and two sites are used as grasslands. Cluster C2 includes predominantly sites located in the upper and middle part of the river located in the Štiavnické vrchy Mountains within the riparian zones or used as grasslands.
The ANOVA test was conducted to determine whether the differences in the content of the different HM fractions in the soils could be significant between the three ecosystem types (RZ, PG, AL) and the soil depths (0-10, 20-30 cm).
Based on the ANOVA results, significant differences emerged between the different ecosystem types for Cd, Cr, and Pb in fractions I and II, for Cr, Mo, Ni, and Pb of fraction IV, and for Co, Cr, Mo, Ni, and Pb of the total content. Within fraction III, the ecosystem types had an effect on all the observed HMs. The most significant differences between the different ecosystem types and the HM fractions were found in the case of Cr, followed www.nature.com/scientificreports/ by Pb, Cd, Mo, Ni Co, Cu, and Zn ( Table 7). The effect of soil depth on HM fractions was revealed only in the case of Cr(I) and Co(III).

Discussion
The high content of HMs in different fractions in floodplain soils along the entire Štiavnica River watercourse proved to be a consequence of past mining activities. Many regions in Europe suffer from soil contamination caused by past mining activities 24 . On the other hand, in many developing countries, soil HM contamination has also become a severe problem due to rapid industrialization over the past several decades 25 . Based on the findings, Cd, Cu, Pb and Zn are the most serious pollutants of the Štiavnica River floodplain. These metals are listed as priority control pollutants by the US Environmental Protection Agency due to their potential toxic, persistent, and irreversible characteristics. Therefore, their excessive accumulation in soils is of Table 6. Pearson´s correlation coefficient matrix of HMs in various fractions and soil physical properties. I exchangeable fraction, II reducible fraction, III oxidisable fraction, IV residual fraction, T total content, *r values shown in bold are significant at p < 0.05, **r values shown in bold are significant at p < 0.01, N.C no correlation.  27 for the assessment of agricultural soils contamination in the EU. Contamination factors confirmed that the contamination by HMs at many AS localities was similar or even higher than at CS localities. The mine tailings contain residues characterised by a high content of HMs. The tailings were usually deposited untreated and left without improper management, being so unstable and prone to wind and water erosion 28,29 . HMs were generated during active mining in the eighteenth and nineteenth centuries and later released and transported by watercourse, and subsequently dispersed during floods along the entire Štiavnica River floodplain. This fact was revealed by HCA showing grouping localities on the basis of their similarities. Many of the study sites were under considerable ecological risk to the environment. Nevertheless, some of the localities located in the downstream part of the Podunajská pahorkatina Hills are still used intensively in agriculture. The natural characteristics of the basin and the water course, including steepness and flow rate, have a significant impact on contamination, and they spread to the more distant surroundings. In the Štiavnické vrchy Mountains, the Štiavnica River flow is steep, narrow, and relatively shallow. In the downstream part of the basin in the Podunajská pahorkatina Hills, the Štiavnica River flow is milder, wider, and deeper. This situation supports heavy metal accumulation downstream, mainly in the part of the floodplain used for agriculture. In addition to the influence of anthropogenic activities, the higher content of clay particles could be the important environmental factor causing the absorption and adsorption of HM by soil. HM accumulated to toxic levels in agricultural land can potentially pose an environmental risk to other ecosystems and food production. Although the Štiavnica River is rather small and reflects the ecological conditions of its local environment and has an impact on the surrounding ecosystems and larger watercourses. Small rivers can have a big effect on larger water bodies and their surroundings 30,31 .
HM toxicity is related not only to the total content in the soil, but also to the distribution of its fractions. It is accepted that a large and relatively stable portion of HMs is present in the crystal lattice of minerals and the residual fraction. On the contrary, it is the mobile exchangeable fraction that can be used to evaluate the extent of environmental bioavailability. Among the HM fractions, the exchangeable fraction determines the most the actual environmental risk, and the reducible and oxidizable ones are the potentially bioavailable fractions. FES in floodplains is directly dependent on the capacity of soil structures to immobilize these mobile and/or potentially mobile and bioavailable fractions. Therefore, it is important to control the bioavailable fractions also for soil remediation 32 . In terms of bioavailability, various species of metals are biologically more available than others. In our study, the content of Cd, Pb, and Zn in all, and Cu in most localities exceeded the Slovak threshold values for mobile forms of HMs in agricultural soils. This suggests that these elements can be readily available to plants and soil organisms 33 . Most HMs are toxic to humans, even at low contents.
The content of HMs in the residual fraction dominated at the RS locality. HMs bound in the residual fraction are naturally present in the parent rock. HMs of anthropogenic input tend to reside commonly in the first three fractions 34 and this is what we observed at the CS and AS localities in the case of Cd, Cu, Pb, and Zn. In the case of Cd, the results indicated the relative importance of Fe-Mn oxide bound, easily soluble, and exchangeable fractions for Cd 35 . The usual distribution of Cd in contaminated soils is exchangeable > reducible > oxidisable > residual fraction 36,37 what, in our study, was more similar to the HM fractionation at the CS localities. In the case of Cu, we observed the great importance of organic matter bound, similar to what other authors have observed 26,38 . In the case of Pb, we found the relative importance of Fe and Mn oxide bound, in agreement with the findings of many authors [39][40][41] .
The correlation analysis allowed us to indicate the key factors that determine the binding of HMs in the soils and their mutual interactions. Many studies have dealt with the correlation relationships between HMs www.nature.com/scientificreports/ and physicochemical properties [42][43][44] . However, the results obtained are not unequivocal 45 but can still serve as a valuable source of information for the management of contaminated sites. Although soil pH is considered to affect the chemical composition of soil solution and also the bioavailability of heavy metals 46,47 , we found also no correlations between the I HM fraction and the pH in the case of Cd, Co, Ni and Pb. Furthermore, in our study, no correlations with Eh were found, although heavy metal behaviour is strongly dependent on redox gradients 46 . Fe and Mn can strongly adsorb divalent HM cations, and their positive correlation with HMs is more prevalent in the studies 48 similar to our study. Negative correlation is less common 49 .
In our study, we confirmed the importance of the content of clay and silt particles as well as SOC as sorbents with great absorptive capacity that support FES. This is consistent with the findings of other authors 50 who claimed that SOC and clay contents 51 have substantial influences on the variability of the HMs. The grain size is one of the essential factors influencing HM contents in soil 52 . The I fraction of Cu, Mo, and Zn is negatively correlated with the content of clay or silt particles, and the I fraction of Cu is negatively correlated with SOC. Thus, a higher clay and silt content can immobilize metals, such as Cu, Mo, and Zn, and a higher organic content Table 7. Effect of ecosystem type and soil depth on heavy metal fractions. The P-values highlighted in bold are highly statistically significant (P < 0.001) and statistically significant (P < 0.05). www.nature.com/scientificreports/ can immobilize metals such as Cu. This finding confirmed the importance of organic soil amendments used to immobilize HMs in soils by changing the speciation from initially highly bioavailable forms to much less bioavailable fractions associated with organic matter 53 . Organic matter is known to have significant immobilizing effects on HMs including Cu 54,55 . Therefore, the retention of HM induced by organic amendment is, on the whole, mainly attributed to an increase in surface charges 56 and to FES support. HM interactions with soil hydro-physical properties are not observed and evaluated so often. But these interactions exist and even poor physical properties such as lower soil porosity, and reduced water holding capacity may indirectly result from HM contamination 57 . Zhang et al. 44 similarly to our results, consider soil porosity as one of the dominant factors determining the distribution of HMs and their fractions in the micropore-dominated riparian soil.
We recorded significant differences between different ecosystem types for all HMs in different fractions but mainly for Cr and Pb, and no effect of soil depth on HM fractions was revealed, except for Cr (I) and Co (III). This indicates that different types of ecosystem and land-use can influence HMs and their fraction content in soils, as documented by many authors 58,59 .

Conclusions
Based on our findings, there is a significant degree of HM contamination, particularly Cd, Cu, Pb, and Zn in the floodplain soils around the Štiavnica River. This pollution originated mainly from past mining activities. The entire downstream transport of suspended solids was a source of HMs for distant soils and, furthermore, pollution caused serious contamination to distant agricultural land in Podunajská pahorkatina Hills.
We comprehensively analysed soil properties that affect HM fractionation and bioavailability. The distribution of HMs in the different fractions, mainly in exchangeable ones, suggested that the mobility and bioavailability of the metals probably declined to the order of Cd, Zn, Cu, Ni, and Pb. The dominance or significant portion of HM in the exchangeable fraction indicates a lower or higher FES capacity provided by floodplain soils, also determined by other soil properties. Soil properties were selectively correlated with the HM fractions. Based on the ANOVA results, the effect of different ecosystem types on the HM fractions was revealed. Thus, soil properties and ecosystem types confirmed a special role in the fate of heavy metals.
These results highlight the importance of considering the HM fractionation as well as the soil properties of different ecosystem types, when making strategic decisions to arrange anthropogenic activities and prevent risk to the environment and human health. There is an urgent need to preserve existing floodplains as natural resources and begin to restore them. This study could help managers adopt a more ecosystem-based approach to floodplain management by helping to reduce pollution. Furthermore, we should consider the adverse influences on ecosystems, living organisms, and public health.