*Review* **Trace Metal Contamination of Bottom Sediments: A Review of Assessment Measures and Geochemical Background Determination Methods**

**Nicole Nawrot 1, \* , Ewa Wojciechowska 1 , Muhammad Mohsin 2 , Suvi Kuittinen 2 , Ari Pappinen 2 and Shahabaldin Rezania 3**


**Abstract:** This paper provides an overview of different methods of assessing the trace metal (TM) contamination status of sediments affected by anthropogenic interference. The geochemical background determination methods are also described. A total of 25 papers covering rivers, lakes, and retention tanks sediments in areas subjected to anthropogenic pressure from the last three years (2019, 2020, and 2021) were analysed to support our examination of the assessment measures. Geochemical and ecotoxicological classifications are presented that may prove useful for sediment evaluation. Among the geochemical indices, several individual pollution indices (*CF*, *Igeo*, *EF*, *Pi* (*SPI*), *PTT*), complex pollution indices (*PLI*, *Cdeg*, *mCdeg*, *Pisum*, *PIAvg* , *PIaAvg* , *PIN*, *PIProd* , *PIapProd* , *PIvectorM*, *PINemerow*, *IntPI*, *MPI*), and geochemical classifications are compared. The ecotoxicological assessment includes an overview of Sediment Quality Guidelines (SQG) and classifications introduced nationally (as LAWA or modified LAWA). The ecotoxicological indices presented in this review cover individual (*ER i* ) and complex indices (*CSI*, *SPI*, *RAC*, *PERI*, *MERMQ*). Biomonitoring of contaminated sites based on plant bioindicators is extensively explored as an indirect method for evaluating pollution sites. The most commonly used indices in the reviewed papers were *Igeo*, *EF*, and CF. Many authors referred to ecotoxicological assessment via *SQG*. Moreover, *PERI*, which includes the toxic response index, was just as popular. The most recognised bioindicators include the *Phragmites* and *Salix* species. *Phragmites* can be considered for Fe, Cu, Cd, and Ni bioindication in sites, while *Salix* hybrid cultivars such as Klara may be considered for phytostabilisation and rhizofiltration due to higher Cu, Zn, and Ni accumulation in roots. *Vetiveria zizanoides* demonstrated resistance to As stress and feasibility for the remediation of As. Moreover, bioindicators offer a feasible tool for recovering valuable elements for the development of a circular economy (e.g., rare earth elements).

**Keywords:** sediments; trace metals; contamination assessment; pollution indices; bioindicators; geochemical background

#### **1. Introduction**

Nowadays, the contamination of sediments with trace metals (TMs) is a salient environmental issue [1,2]. TMs delivered to the environment from anthropogenic sources are mostly represented by these elements: Cd, Co, Cr, Cu, Fe, Hg, Mn, Mo, Ni, Pb, Sn, V, and Zn. Supplementation from external sources also includes other potentially toxic elements (PTEs) such as metalloids (As, Sb) and non-metals (Se) [3]. PTEs are not biodegradable, persistent, and can accumulate in the environment in large amounts. TM ions cause health hazards; just to mention a few, these include emphysema (Cd), anaemia (Cd, Zn, Pb), damage to the brain (Pb), kidneys, and bones (Cd), necrosis nephritis (Cr), DNA damage

**Citation:** Nawrot, N.;

Wojciechowska, E.; Mohsin, M.; Kuittinen, S.; Pappinen, A.; Rezania, S. Trace Metal Contamination of Bottom Sediments: A Review of Assessment Measures and Geochemical Background Determination Methods. *Minerals* **2021**, *11*, 872. https://doi.org/ 10.3390/min11080872

Academic Editors: Ana Romero-Freire and Hao Qiu

Received: 16 June 2021 Accepted: 10 August 2021 Published: 12 August 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

(Ni), eczema on the hands (Ni), risk of developing cancer (As), central nervous system irritation (Cu), and even death in humans [4]. Due to the toxic nature of PTEs, there is an obvious need to monitor their content. The most significant anthropogenic sources of PTEs are related to industry, urbanisation, transportation, agriculture, and mines [5–8] TMs adhered to particulate matter (PM; including PM10 and PM2.5) could be released into the surrounding environment causing contamination of soils, groundwaters, sediments, and surface waters. PTEs entering the aquatic ecosystem are associated with fine-grained fractions of suspended solids and sediments due to their large surface areas and high sorption capacities [9]. Throughout the hydrological cycle, 99% of PTEs are stored in sediments, which, therefore, act as the main sinks and carriers for contaminants in aquatic environments [10]. Contemporary studies [11–18] have revealed that stormwater run-off carried metals originating from various sources of everyday use associated with tyre wear, corrosion, roof run-off, and fuel combustion products. Zn, Pb, Cu, Cd, Ni, Ce, and As are typically observed in urban and industrial catchments [11,12].

Sediments are the ultimate repository for contaminants and provide valuable information about the environmental status of a catchment area [17]. According to USEPA [19], sediment site management and silting processes in reservoirs present a unique challenge because of the large number of difficult to control contaminant sources, dynamic phenomena influencing movement and sediment stability, dispersion of contaminants, diversified use of areas affecting the sediment environment, as well as difficulties involved in engineering cleanup works.

Sediment contamination may be evaluated in several ways. Geochemical, ecotoxicological, and bioindicator-related approaches can be distinguished. Extensive evaluation includes geochemical and ecotoxicological methods. The geochemical assessment relies on comparing TMs at particular sites to the geochemical background or baseline. The term "geochemical background" (or "natural background") has evolved. In exploration geochemistry, it refers to the non-anomalous concentration of elements in the rocks surrounding the mineral ore body. In environmental sciences, it is used to detect anthropogenic impact on the environment [20]. Establishing the geochemical background raises research concerns because it may determine the overall pollution status of an examined site. Ecotoxicological approaches are based on assessing the effects of metal toxicity on living organisms and determining certain dose limits. Bioindication is an indirect method of verifying the quality of the environment—i.e., by using aquatic plants and hydrophytes able to accumulate a large number of metals, making them suitable for environmental biomonitoring. The distribution and behaviour of aquatic plants are frequently correlated with water and sediment quality [21]. Macrophytes are able to accumulate pollutants at a higher level, irrespective of their content in the surrounding environment [22]. Moreover, the analyses of plant tissues provide time-integrated information about the quality of an examined site, even if contaminant inputs are discontinuous and quickly diluted in water [23]. One of the species widely recognised as featuring correlations between concentration of metals in sediment/soil and roots is *Phragmites australis* [15]. Hydrophytes uptake PTEs through their roots and rhizomes, stems, and leaves [24]. Plant species with a strong ability to accumulate TMs in the aboveground tissues are good candidates for phytoextraction [25].

The most simple way to evaluate TM contamination of sediments is based on total metal content. This approach reflects the geological origin and anthropogenic influence, although it is ineffective in terms of documenting mobility or bioavailability [26]. According to Ure and Davidson [27], a more accurate method taking into account these two features is to measure the "pseudo total" metal content in sediments. Pseudo total content presents the maximum potential soluble or mobile metal content, which can be identified with maximum potential hazard. Sequential extraction (or fractionation) distinguishes several factions—mobile, conditionally mobile, and immobile (depending on the applied method). Therefore, using this approach, potential metal contamination hazards can be more accurately assessed and explored. Combining fractionation with different contamination

assessment approaches (geochemical and ecotoxicological) may also offer valuable insight into sediment quality.

Sediment contamination is a well-known problem across the world; however, monitoring programs based on different PTE contamination assessment methods as well as bioindicators have to be implemented to ensure that the assessment is credible and reliable. Various studies use different approaches to evaluate the quality of sediments, which impairs the comparison of results from diverse studies. Even the definition of background values, which is a fundamental step in geochemical assessment, is often controversial. Therefore, there are four aims of this review: (1) to summarise the levels of TM content and methods used to evaluate the contamination status of bottom sediments related to surface run-off receivers from the last 3 years, (2) to assemble and critically evaluate the geochemical and ecotoxicological classifications and contamination indices used to evaluate the PTE-polluted sediment and the subsequent risk, (3) to explain the difference between geochemical background and baseline; and (4) to specify the bioindicators of TMs that have been extensively researched at different sites across the world.

#### **2. Materials and Methods**

A literature review was conducted via the Web of Science platform. In order to access a complete spectrum of data on PTEs in the bottom sediments of urban water bodies, 5 specific keywords were chosen: sediments, trace metals, contamination, pollution indices, and urban. The Web of Science Core Collection displayed a total of 519 papers, which are available in the database (starting from the year 1994). From the last 3 years (2019, 2020, 2021), there are 202 publications in total, and a further 145 results in the field of environmental science (as of May 2021). We used 25 of the latest publications from the last 3 years (2019, 2020, and 2021). To explore the last 10 years of updates regarding bioindication species, the search process was based on the following keywords: plant bioindicators, metal indicator, *Salix* or *Willow* as bioindicator, phytoextraction, and bioindicator.

#### **3. PTEs in Bottom Sediments of Urban Water Bodies**

Based on 25 of the latest publications [28–52] on the topic of bottom sediments in urban waterbodies, Table A1 in Appendix A presents (1) a list of TMs evaluated under each study, (2) their concentration ranges (in mg/kg d.w.), (3) a basic description of the site and land use, (4) indices, classifications, and an evaluation of TM bioavailability with the use of metal fractionation.

In most of the examined papers, the authors focused on five TMs: Cu, Zn, Pb, Ni, Cr, and a metalloid: As. All of the mentioned elements can be delivered by anthropogenic activities such as metalliferous mining, smelters, metallurgy, rolling, electronic, industry, the dye and paint industry, metal corrosion of materials, plastic industry, chemical industry, combustion of fossil fuels, wood industry, as well as tyre wear and lubricant oils [6]. The scope of analyses in many studies also included Fe, Mn, Co, Hg, as well as rare earth elements, i.e., La and Sc.

Most authors sampled the surface layer of sediments (0–2 cm, 0–5 cm, or 0–10 cm), where reactions occur at the water–sediment interface and where living organisms dwell. The sampling strategy is crucial for the successful collection of sediment samples with undisturbed stratification. With this aim, a grab sampler (for surface sediment extraction), hand (core extraction < 60 cm) and box corer (core extraction < 1 m) devices, or a multiple corer can be used to obtain proper samples. Sediment samples are collected to describe the condition of the sites. The specific reasons for sediment sampling include bioassays (the top 10 cm are usually collected because of biological activities), biosurvey (macroinvertebrate analyses), monitoring, examining non-point pollution, and pollutant sources. Monitoring trends and relationships in sediments requires information on current PTE burden in the uppermost sediment layer (e.g., 0–1 cm). Depending on the receiver (large/small body of water, river, stream, bay), the rate of deposition varies, which means that 1 cm of sediment builds up over different times. Therefore, the EU [53] has recommended sampling the

top layer of sediment at a depth of 1–5 cm; however, this depends on the deposition rate. This is also in line with data from the United States, which stated that more than 10% of the volume of bottom sediments (the upper 5 cm of surface sediments) are estimated to have been contaminated. The depth to which sediments are sampled should be one of the key monitoring objectives. Core sample collection is crucial for evaluating the PTE footprint. Sediments from depths greater than 15 cm below the sediment–water interface are usually collected to determine the spatial vertical properties and sedimentation history. Each sampling program requires a compromise between the number of samples, sampling frequency, measurement parameters, and cost effectiveness [54].

The most frequently used indices were Geoaccumulation Index (*Igeo*) (17 studies), Enrichment Factor (*EF*) (13 studies), Contamination Factor (*CF*) (10 studies), Pollution Load Index (*PLI*) (6 studies), Metallic Pollution Index (*MPI*) and Modified Contamination Degree (*mCd*) (each used in 3 studies), Nemerow Pollution Index (*PINemerow*) (2 studies), *MRI* and Pollution Index (*PI*) (each used in 1 study). Two studies undertook a sediment assessment using the German classification LAWA, while one study refers to Romanian legislation. Reference to the ecotoxicological assessment via Sediment Quality Guidelines (SQG) was carried out in eight studies, while Ecological Risk Index (*Er<sup>i</sup>* ), Risk Index (*RI*), and Potential Ecological Risk Index (*PERI*) (which account for metal toxicity) were used in 17 studies. Moreover, only two studies assessed health risk. Three papers considered the use of chemical sequential extraction for element analyses, which is important to assess the bioavailability and mobility of TMs. Only one study included a Risk Assessment Code (*RAC*) analysis, which takes into account the risk determination based on the metal binding to the exchangeable and carbonate fraction. Health risk—as Hazard Index (*HI*) and Hazard Quotient (*HQ*)—was calculated in two papers. The common use of indices enables the results to be compared between various studies. However, calculations of geochemical indices are strongly dependent on the adopted reference value (geochemical background).

Another important issue, addressed in three studies, is the fractionation analysis of TMs. Fractionation allows the environmental behaviour of trace elements to be assessed [55]. However, the procedure itself is highly labor-intensive and cost-prohibitive. Most works are based on assessing the total metal content of the samples; nevertheless, fractionation is an indispensable method when it comes to exploring the actual metal bioavailability and mobility.

The above characteristic offers an overview of the most often used indices for geochemical and ecotoxicological evaluation of sediments. The next part of this review presents a "conventional" division of tools for assessment of sediment quality, taking into account geochemical and ecotoxicological approaches.

#### **4. Assessment Techniques for Sediments Contaminated with HMs**

The most effective approach in assessing sediment contamination with PTEs relies on the use of several assessment tools that can be grouped according to geochemical and ecotoxicological evaluation, as shown in Figure 1. Another method for indirect assessment of TM sediment contamination is the use of bioindicators.

Geochemical evaluation is usually based on pollution indices (PIs), while geochemical classifications are usually implemented by national governments. Both approaches rely on the establishment of a geochemical background. Among ecotoxicological methods, descriptive classifications (usually based on the ecotoxicological effect on benthic organisms—Sediment Quality Guidelines (SQG) as well as those based on toxicological, chemical, and ecological tests) and ecotoxicological indices (EIs) are used. Since the preference for adopting a sediment assessment method is usually on the side of the lawmakers, many classification and indices have been introduced since 1980.

**Figure 1.** Division of the methods of assessment of sediment contamination with PTEs.

#### *4.1. Geochemical Evaluation*

4.1.1. Geochemical Background and Geochemical Baseline

There are two controversial issues within contemporary geochemistry and environmental biochemistry: the geochemical background and the geochemical baseline of an element or chemical compound. The geochemical baseline is easier to determine—it is commonly used for assessing the environmental status before embarking on a specific investment (e.g., constructing a power plants). In exploratory geochemistry, geochemical background means the natural concentration of an element in the environment, i.e., a non-anomalous concentration of elements that corresponds to Clarke's characteristic for their geochemical environment [20]. There are several methods to determine the accurate geochemical background. Once this value is determined, the pollution status can more precisely be assessed. These methods rely on statistical (indirect, theoretical) and empirical (direct, geochemical) determinations. An integrated approach combining both methods is also in use. Geochemical methods involve analysing samples unaffected by anthropogenic activities, which are also referred to as preindustrial samples. Usually, deep core samples and/or samples collected at a certain distance from pollution sources are used to establish geochemical background values. Statistical methods involve a determination and elimination of outliers, which are identified by standard deviation, regression analysis, or by using the Tukey boxplot method [56].

The geochemical background is not a constant value over time due to natural processes taking place at the interface of water and sediment. Additionally, the geochemical background value is local or regional rather than global. The adopted geochemical background affects the overall assessment environment quality. Sometimes, the reference level follows the shale standard (from 1961), i.e., the average content of elements in the earth's crust [57], although this approach is incorrect.

#### 4.1.2. Pollution Indices

According to Weissmannová and Pavlovský [6], pollution indices (PIs) can be divided into two categories including Individual (single) PIs (IPIs) and Complex (total) PIs (CPIs). In such a division, IPIs are calculated based on the contents of each metal separately in particular sediments. IPIs can be used to classify sediments into several classes according

to the degree of pollution presented by a single metal. CPIs are calculated as multielement indices based on IPIs. Moreover, CPIs can be composed of IPIs separately. PIs are of most use when they meet several criteria: (1) they are easy to calculate, (2) they make reference to the geochemical background, and (3) they are widely used, which makes it possible to compare sediments from different sites.

#### Individual Pollution Indices (IPIs)

The IPIs used for sediment assessment are presented in Table A2 in Appendix B. In this group, five indices were identified: Contamination Factor (*CF*) [58], Geoaccumulation Index (*Igeo*) [59], Enrichment Factor (*EF*) [60], Pollution Index (*P<sup>i</sup>* ) [15,54], and Threshold Pollution Index (*PIT*) [28]. Each of these indices compares the concentration of a given metal in a sample (*CmSample*) to its comparative concentration—represented by *CmPre* (preindustrial concentration of TM) or *CmBackground* (geochemical background concentration of TM) [40,58,61–64].

*CF* is originally defined as a ratio of mean TM content averaged from at least five samples (*CmSample*) and metal concentration in preindustrial sediments [58,61,62]. According to Dung et al. [20] the ratio can also refer to background TM content. This approach is included in *P<sup>i</sup>* . However, usually both indices (*CF* and *P<sup>i</sup>* ) are treated by authors interchangeably. The difference lies in the classification of the contamination levels; it also depends on whether the values of *CmPre* and *CmBackground* are consistent with each other. *CF* and *P<sup>i</sup>* are easy to calculate, although both omit natural variability, grain-size, and reference elements ratios. Moreover, they do not take into account metal fractionation even though the formulas of the equations themselves make it possible.

*Igeo* proposed by Müller [59] is a commonly used index, which uses a completely different approach from other IPIs, since it defines seven classes of TM sediment contamination. The TM content in a sample is referred to as *CmBackground*. The formula includes a constant value of 1.5, which is used to offset the natural TM fluctuations in the environment and detect minor anthropogenic impact [6].

*EF* is an effective tool for TM content comparison [13]. In comparison to other IPIs, *EF* requires the most extensive data input, such as geochemical background concentration of TM (*CmBackground*) and reference element concentration in a sample (*Cre f*) as well as in the reference environment (*Bre f*). The reference element has a low occurrence variability. The reference element—and this could be Sc, Mn, Ti, Al, Fe, and Ca—has a low occurrence variability. Frequently used reference metals are Mn, Al, and Fe [12]. Usually, Fe has a relatively high natural concentration in comparison to other metals and is therefore not expected to be enriched from anthropogenic sources [64]. The use of *EF* seeks to reduce the metal variability associated with the granulometric composition (mud/sand ratios). The *EF* could be also expressed as a percentage.

*PI<sup>T</sup>* is based on the single Pollution Index (*P<sup>i</sup>* ), which refers the TM content in a sample to the tolerance levels of metal concentration (*CTL*). *CTL* could be established by national guidelines or criteria regarding health-threatening TM content [63]. The determination of the threshold parameter (*CTL*) is freeform and can be adjusted to the state of the environment and defined individually by national regulations. Therefore, this index may be an effective tool in the hands of environmental managers.

All of the described IPIs omit the metal fractionation and mobility/bioavailability. IPIs are usually calculated in accordance with total metal content; however, they can be easily applied in calculations that include fractionation (as described in Section 4.3).

#### Complex Pollution Indices (CPIs)

A list of CPIs is presented in Appendix C. CPIs mostly rely on *CF* or *P<sup>i</sup>* values for TMs under analysis. Thus, the starting point for a comprehensive metal contamination assessment is IPI analysis. Table A3 in Appendix C presents calculation formulas with classifications and descriptions (pros and cons considered by the authors) for the following CPIs: Pollution Load Index (*PLI*) [62], Contamination Degree (*CDeg*) [58], Modified Con-

tamination Degree (*mCDeg*) [64], Sum of Pollution index (*PIsum*) [63], Average of Pollution Index (*PIAvg*) [63], Weighted average of Pollution Index (*PIwAvg*) [58,59,63], Background Enrichment Factor = New Pollution Index (*PIN*) [65], Product of Pollution Index (*PIProd*) [63], Weighted power product of Pollution Index (*PIwpProd*) [63], Vector modulus of Pollution Index (*PIvectorM*) [63], Nemerow Pollution Index (*PINemerow*) [66], Integrated Pollution Index (*IntPI*) [67], and Metallic Pollution Index (*MPI*) [68,69].

*PLI* is a very easy-to-use index based on *CF* values specified for a selected number of metals. Moreover, the deterioration of site quality could be straightforwardly specified because the unity value of PLI (PLI = 1) is the boundary differentiating the adverse impact on the quality of sediments from non-polluted sites [62]. The overarching merit is a comprehensive assessment of sediments, while the disadvantages are related to the shortcomings of the *CF* index.

*CDeg* and *mCDeg* are also based on *CF* values for analysed TMs. The difference lies in the classification of sediment contamination status. *CDeg* defines sediment contamination referring to four categories, while *mCDeg* distinguishes seven ranges of sediment quality [58,64].

The equations for calculating *PIsum* and *PIAvg* are similar to those for *CDeg* and *mCDeg*. The difference lies in the use of individual *Pi* (*Pi = SPI = PI*) factors related to *CmBackground* (not to *CmPre*) [63]. In order to evaluate the sediments with *PIsum*, it is necessary to specify the classification scale for each individual case. The defined scale should refer to the number of selected TMs used to assess contamination. The pollution evaluation with *PIAvg* is easy to apply due to the reference to unity (1). *PIwAvg* application when the TM weight is equal to unity (*w<sup>i</sup>* = 1) is a special form of *PIsum*.

*PIN* is based on *Pi* (*Pi = SPI = PI*) and requires the TM class to be determined (from 1 to 5 in reference to *Pi*). This very precise scale classifies sediments into five categories [68].

*PIProd* is a Cartesian product of *Pi* applied for a selected number of metals. The classification should be powered by n (count of TMs) to determine the contamination of sediments [63]*. PIwpProd* uses the weight of *Pi* in an equation, which is similar to those used for calculating *PIProd* [63]*. PIvectorM* was introduced by Gong et al. [63] and depends on *Pi* and the number of TMs selected for investigation. *PINemerow* determination emphasises the maximum value of a single Pollution Index (*Pi*) for all TMs; however, it does not include the weight of *Pi* [66]. *IntPI* is an easy-to-calculate index based on the *Pi* mean value. The classification based on *IntPI* divides sediment contamination into low, moderate, and strong [67,70].

*MPI* is equal to *PLI*. The general difference lies in *CmBackground* used in the equation for the *Pi* calculations [71].

In relation to CPIs, one major advantage is the possibility to comprehensibly evaluate the effects of several metals simultaneously. The authors of this review believe that the the *Pi* weight, which represents the impact of an individual metal on the overall sediment pollution rating, plays a significant role in the calculations. This approach is used by *PIwAvg* and *PIwpProd*. Moreover, it would seem to be more accurate to include the maximum *Pi* as a reference in the calculations of *PINemerow*. Despite the mentioned advantages of *PIwAvg*, *PIwpProd*, and *PINemerow*, it is still desirable to use simpler calculations, which are included by *PLI*, *CDeg*, *mCDeg*, *PIsum*, etc.—due to their simplicity and frequent use in different studies.

#### 4.1.3. Geochemical Classifications Referring to the Geochemical Background

Extensive environmental monitoring is required to develop geochemical classification and a database. Often, such studies end by mapping a given area to present geochemical data as accurately as possible [72]. Since determining the geochemical background is usually a fundamental step while assessing the environmental status of a given site, such classifications have a temporary value [20].

In Poland, the first attempts to classify lake and river sediments contaminated with TMs were made by Bojakowska and Sokołowska [73]. This so-called "geochemical classification of aquatic sediments" approach involves standardising samples, taking into account

a grain size smaller than 0.2 mm. Moreover, this classification applies to samples digested with hydrochloric acid, nitric acid, or aqua regia. The geochemical criteria distinguish four geochemical quality classes presented and described in Table 1. This classification was developed in 1998, so the geochemical background could have changed since then. However, it is still used in Poland in the absence of other regulations. In Table 1, the geochemical quality classes are collated with geochemical background values established using different approaches. Depending on the goals, TM contents in analysed sediments can be compared with preindustrial concentrations [58], trace metal content in world soils [74], or shale standard [57].

**Table 1.** The concentration of various trace metals (TMs) in mg/kg d.w. for each class of sediments—geochemical classification of aquatic sediments [57,58,73–75].


Bold—geochemical quality classes description according to Bojakowska and Sokołowska: Class I—uncontaminated sediments (boundary value of elements concentration determined as from 5 to 10 times higher than the geochemical background); Class II—moderately contaminated sediments (TMs content from 10 to 20 times higher than the geochemical background); Class III—contaminated sediments (TMs content from 20 to 100 times higher than the geochemical background); Class IV—highly contaminated sediments [73].

#### *4.2. Ecotoxicological Evaluation*

In order to assess the level of contamination that is still safe for living organisms, it is necessary to apply ecotoxicological criteria or use indicators that take into account the level of metal toxicity.

#### 4.2.1. Descriptive Classifications Based on Total TM Content

Internationally, a range of sediment quality guidelines (SQGs) have been proposed for contaminants. Originally, they provided the information on SQGs compared with a reference or background concentration of analysed substances [76]. After the 1980s, SQGs were developed to evaluate sediment quality where harm to aquatic organisms and the environment was observed [77]. In general, SQG approaches that use different criteria and factors are divided into two main categories [76]:


Both types of SQG predict adverse ecological effects caused by sediment contamination by the response of benthic organisms. SQGs feature two concentration thresholds: the first below which adverse effects are not expected to occur (threshold effect concentration— TEC) and the second above which adverse effects are expected to occur more often than not (probable effect concentration—PEC). The limitation of this approach is that it leaves significant uncertainty in the "grey" area of contaminant concentrations that lie between TEC and PEC. Therefore, it may be necessary to perform a site-specific analysis. Another relevant issue in the case of SQGs is that these descriptions apply to only one type of

contaminant, which may lead to inaccurate conclusions in the contamination description, as the sediments are a sink for various type of contaminants that can contribute to their overall toxicity. Table 2 outlines the most popular empirical SQGs used in practice. All of them are based on differently described terms of TEC and PEC (e.g., TEL/PEL, ERL/ERM, or LEL/SEL) [78–80]. Moreover, we presented the average TEC and PEC values calculated based on the literature (TEClit and PEClit).

**Table 2.** Sediment quality guidelines (SEQs) for TMs and metalloid in mg/kg d.w. that reflect threshold effect concentration (TEC) and probable effect concentration (PEC) [77–82].



TEC—Threshold-Effects Concentration; PEC—Probable-Effects Concentration; TEL—Threshold-Effects Level; PEL—Probable-Effects Level; ERL—Effects Range: Low; ERM—Effects Range: Median; MET—Minimal Effect Threshold; TET—Toxic Effect Threshold; SQGV—Sediment Quality Guideline Value; ANZECC/ARMCANZ [82]—currently used classification in Australia and New Zealand.

> In parallel with the intense work undertaken on SQGs in North America and Australia, several approaches and descriptive classifications and recommendations were also developed in Europe. It is worth noting that the authorities and legislators usually consider water quality as a priority, paying less attention to the state of sediments. It should be emphasised that contaminants accumulated for decades in the bottom sediments of water bodies could be released into the water column and cause its recontamination, thereby directly affecting water quality. Despite the large-scale processes that can affect the stability of sediments, there are different inceptors of secondary water recontamination, such as pH fluctuation, or the disturbance of sediment structure during dredging operations or boat movement.

> In Germany, the Länder-Arbeitsgemeinschaft Wasser classification (LAWA) [83] is used to assess Hg, Cd, Zn, Pb, Cu, Cr, and Ni content in water, suspended solids, and sediments. The LAWA classification is constrained to the methodology of digestion with the use of aqua regia. This approach does not allow silicates to dissolve; therefore, the term "pseudo total" concentration of TMs is used in the literature. However, in terms of the actual contamination status of sediments and their impact on the aquatic environment, the mobile and unstable forms of TMs are most relevant [75]. In the light of this, digestion by using aqua regia or diluted acids (e.g., hydrochloric acid, nitric acid) may determine the bioavailable fraction of TMs that take part in environmental processes, while it remains impossible to determine the content of TMs bound in minerals resistant to weathering processes. LAWA classification divides sediments into seven classes of contamination, depending on the TM contents in aquatic ecosystems. The boundary values for each class are presented in Table 3.


**Table 3.** The concentration of various trace metals (TMs) in mg/kg d.w. for each class of sediments—classification LAWA and modified LAWA [83,84].

\* Quality goal—Class II. Class I refers to uncontaminated sediment; Class I–II—sediments unpolluted or with very small anthropogenic interference; Class II corresponds to moderately polluted sediments, and it is a reference for the other levels of pollution; Class II–III—twice as much as Class II—defines moderate to significant contamination of sediments; Class III—four times more than Class II—defines significant pollution; Class III–IV—eight times more than Class II—refers to very strong pollution; Class IV—means more than eight times exceeding Class II—ultimate pollution. Comparative value is the average contents of TM expressed in mg/kg [83,84].

> Following the LAWA classification, Ahlf et al. [84] introduced guidelines for evaluating sediment quality based on integrated toxicological, chemical, and ecological approaches (hereafter referred to as modified LAWA). This approach is based on a series of three different tests that should be performed step-by-step, starting from an ecotoxicological assessment. If the ecotoxicological test confirms sediment contamination, the second step is introduced, which focuses on chemical analyses. The final assessment stage involves an additional and thorough analysis of factors that may harm the environment. The classification of TM sediment contamination applied by Ahlf et al. [84] (Table 3) combines temporary valuation standards of river mud and ATV standards concerning sewage and garbage. The final evaluation is based on the concentration limits for class II, which constitutes a point of reference.

> In accordance with the Water Framework Directive, water and sediment quality assessment considers the extent of deviation from reference conditions with no or with very minor anthropogenic influence. The general conclusion from analyses of sediments guidelines is that they have advantages as well as limitations for their application in sediment quality assessment [85].

#### 4.2.2. Ecotoxicological Indices (EIs)

EIs as indicators of ecotoxicology comprise numerous indices, such as Ecological Risk Factor (*ER*= *Er<sup>i</sup>* = Risk Index (*RI*)) [58]—the only IEI, Contamination Severity Index (*CSI*) [86], Sediment Pollution Index (*SPI*) [58,64], Risk Assessment Code (*RAC*) [35], Risk Index (*PERI*) [87,88], and the probability of Toxicity Index (*MERMQ*) [86]. The list of IEIs and CEIs are presented in Table A4 Appendix D with classification and a summary of strengths and weaknesses.

*ER* describes the ecological risk caused by TM. *T i r* in the *ER* formula represents the toxic-response factor for a given metal, and so this index ties the risk assessment with the metal content in a sample [6]. *SPI* also refers to *T i <sup>r</sup>* of TMs and additionally takes into account the average shale concentration of TMs. Moreover, *SPI* is tailored for the comprehensive evaluation of several elements at the same time. Another difference between *ER* and *SPI* lies in diverse points of reference—*ER* classification refers to ecological risk, and *SPI* refers to the sediment pollution level. *PERI* is a complex index based on a single

*ER*. The toxic response factors for TMs according to Håkanson [58] are as follows: Hg—40, Cd—30, As—10, Cu—5, Pb—5, Cr—2, and Zn—1. *T i r* in the limnic system was discussed and established by Håkanson according to the "abundance principle" of elements in different matrixes. The element with the highest mean concentration must be ranked as 1.0; according to this assumption, analysed elements were marked with the following abundance number: Zn = 1 < Cu = 3.4 < Pb = 13 < Cr = 110 < As = 140 < Cd = 230 < Hg = 1160. Note that the abundance number is not equal to the toxic factor. Håkanson corrected the abundance number to reflect the "sink-effect". This effect is related to different footprints caused by substances (e.g., Cr has the lowest sink-factor = 2, which implies that this element leaves the heaviest footprint in the sediment, while Hg has the highest sink-factor =320, which means that much more Hg can be found in water compared to sediments). Taking into account the dimension problem, Håkanson corrected the abundance numbers that correspond to *T i r* . Moreover, according to Kabata-Pendias and Pendias [58], the highest *T i r* corresponds to strong chemical and biological activity in the case of Hg and high solubility in acidic environments and mobility for Cd. Quite low *T i r* for Pb is related to its low mobility. *T i <sup>r</sup>* = 2 for Cr is related to the principle that only Cr(VI) is toxic; at the same time, Cr(VI) is easily reduced to poorly soluble Cr(III), which is generally difficult for plants and benthic organisms to absorb [58].

*CSI* requires the use of computed weight of each TM according to Pejman et al. [86]: (Cu—0.075, Zn—0.075, Cr—0.134, Ni—0.215, Pb—0.251, Cd—0.25). The weight of TMs may result from statistical analyses (e.g., Principal Component Analyses—PCAs and/or Factor Analyses—FAs). Moreover, *CSI* requires the application of ERL and ERM, and at the same time, it is highly accurate and could provide a precise evaluation of sediments.

*RAC* depends on exchangeable and carbonate fractions determination (determined with CH3COOH 0.11 M solution) [35]. The classification refers to the percentage of mobile fractions in relation to the total concentration of bound TMs in all fractions together. It is worth noting that *RAC* is the sole index that combines the speciation analyses of TMs. Therefore, it is very valuable for evaluating the mobility and bioavailability of metals.

*MERMQ* is also referred to in the literature as a combined effect of Toxic Metals (*mPECQ).* This complex EI is based on TM content in the sample and the second threshold level of effects (e.g., ERM for *MERMQ* and PEC for *mPECQ)* [86].

#### *4.3. Extended Sediments Quality Ranking*

### 4.3.1. Integrated Geoaccumulation Index—*Igeointegrated*

The approach described by Von Tümpling et al. in 2013 [89] relies on *Igeo* introduced by Müller [59], and this index is combined with metal speciation analysis based on BCR (Community Bureau of Reference) three-stage sequential extraction—this approach is referred to as "extended use of *Igeo* in combination with BCR fractionation". To compare the results obtained with *Igeointegrated* between fractions, an integrated *Igeo* was introduced (Σ*Igeointegrated*). The sum of *Igeointegrated* offers a general overview of the contamination degree of the considered elements in one fraction in comparison to Σ*Igeointegrated* of other fractions. The Σ*Igeointegrate*<sup>d</sup> equation is as follows:

$$\sum Igeo\_{integrated} = \sum\_{j=1}^{m} \sum\_{k=1}^{n} m\_{j,k}$$

where *m* is the quality class, *j* is the sampling point (SP), and *k* is the element. The geoaccumulation classes (described in Table A2 Appendix B) of all the examined TMs are added for each fraction. If there is no contamination or no accumulation of an investigated element, the resulting sum is 0. The maximum class is 6.

The *Igeointegrated* approach can completely change the result of the sediment assessment based on the total TM content performed by *Igeo*. The results reported by Von Tümpling et al. [89] for Klinke River in Germany indicated that F3, F2, and F1 had the greatest Σ*Igeointegrated*. At the same time, the "good chemical status" classification in accordance with total TM content was overturned by the sequential extraction method, which demonstrated that changes in element speciation can lead to higher TM mobility, which may increase the hazardous potential of sediments.

#### 4.3.2. Fractionation of TMs

In order to establish how TMs behave in the sediment–water or soil–water environments, the presence of physical and chemical forms of metals—i.e., fractionation, also called speciation—must be determined. Therefore, speciation is an identification of various physicochemical forms of elements in the analysed material and/or mutual proportions of these forms. This approach helps to understand the mobility of metals in the environment, although the sample processing is time consuming [90].

The analytical methods for the chemical section of TMs in the solid phase involve single reagent leaching, sequential extractions, and ion exchange resins [55]. In general, the processes involving the separation of different species are based on selective chemical reactions and/or extraction procedures (e.g., liquid–liquid, solid-phase extraction). The final result of element contents is presented in mg/kg and or mg/L (ppm—parts per million) or ng/kg and/or ng/L (ppb—parts per billion). Therefore, the overall process is complex, arduous, and susceptible to critical errors that are sometimes neglected. In the field of geochemistry, sequential extraction has proven its value. Nevertheless, accuracy at every stage of sample processing and analysis is crucial to guarantee an acceptable level for detecting (LOD—limit of detection) and quantifying the analytes [91]. TMs may be determined by commonly used instruments such as flame (FAAS), hydride generation (HGAAS), and electrometry (ETAAS) atomic absorption spectrometry, as well as inductively coupled plasma mass spectrometry (ICP-MS) and emission spectrometry (ICP-OES). In comparison with procedures in which total TM contents are determined in a direct manner, the extraction methods generate more uncertainties. This inconvenience is caused by difficulties in isolating compounds from substrates, upsetting the equilibrium between different chemical species, insufficient analytical sensitivity, or a lack of certified reference materials [92].

Many selective extraction methods have been developed, among others, by Tessier et al. [93], Boszke [94], Miller et al. [95], and Fiedler et al. [21]. One method that seeks to minimise errors during sample processing and analysis was proposed by the European Community Bureau of Reference—the BCR method [96].

The BCR three-step extraction procedure is similar to those proposed by Tessier et al. [90]. The key difference involves combining the first fraction of the procedure (instead of evaluating the exchangeable and carbonate fractions separately, the BCR procedure combines both); therefore, BCR distinguishes four factions (F1—exchangeable, F2—reducible, F3—oxidisable, and F4—residual). The description of subsequent steps of the BCR extraction procedure is presented in Table 4.


**Table 4.** The BCR three-steps sequential extraction scheme [90].

#### *4.4. Bioindicators*

In recent decades, the interest in using bioindicators as monitoring tools to examine environmental pollution caused by TMs has been growing steadily. Numerous plants are used as bioindicators to construct assumptions regarding environmental conditions. In this context, plants including macrophytes can be used as bioindicators for TM-related environmental pollution due to their ability to absorb metals from air, water, sediment, soil, and the food chain [97]. In general, bioindicators are characterised by the following abilities: they are easy to grow and cultivate, can withstand polluted environments, and can accumulate high pollutant concentration in their biomass; they are ecologically and economically important for society, and they should be abundant and widespread. In addition, bioindicators not only provide sediment, soil, and water remediation but are also a feasible tool for the recovery of valuable elements to develop a circular economy [98]. In practice, aquatic and related fauna, microbial systems, fungi, animals, and plants are popular targets of TM monitoring programs used as bioindicators to formulate conclusions regarding the environmental contamination status. Biological indicators provide information on the long-term impact of TM contamination [97]. There are many reasons to use a given species as bioindicators, while the key issue lies in the correlation of analysed toxic compounds/elements between species and contaminated matrix (as waters, sediments, or soils) [99].

In this review, different bioindicator plant species were identified from different studies based on their survival, highest metal accumulation, and biomass production in different polluted environments (soil and water) (Table 5). Moreover, metals such as Zn, Cd, Cu, and Ni were predominantly investigated in most regions. For instance, Salinitro et al. [100] depicted the linear relation of Ni in soils and in *Poa annua* shoots (R <sup>2</sup> = 0.78) and *Senecio vulgaris* (R<sup>2</sup> = 0.88). Therefore, these species have been demonstrated to be reliable indicators of both total and bioavailable Ni fractions in an anthropogenic environment. Extensive research on bioindicators in an estuarine and coastal environment was performed by Farias et al. [98]. In this study, *Ulva australis*, *Zostera muelleri,* and *Ruppia megacarpa* were selected for bioindicator identification due to their great abundance and the highest content of TMs inside their tissues. These species were verified as potential bioindicators of TM (As, Cu, Pb, and Zn) pollution in estuary sediments.

Some authors also focused on highly toxic elements such as Hg and As. These nonessential elements are uptaken passively by plants [24]. Many plant species presented tolerance to As contamination, as presented by Singh et al. [101] based on the example of *Vetiveria zizanoides*, which demonstrated resistance to As stress and feasibility for revegetation/remediation of As-contaminated soils. Although Hg in plants has no metabolic function, it could be easily translocated from roots and rhizomes to shoots [102].

High resistance to changes in salinity, fertility, textures, and pH is presented by an extensively distributed species Phragmites australis [103]. Ganjali et al. [104] demonstrated that Phragmites roots can be used as bioindicators of Fe, Cu, Cd, and Ni.

*Phragmites* and *Salix* species were studied widely for the phytoremediation of TMs in temperate and boreal climates. *Salix schwerinii* and *S. myrsinifolia* are known as bioindicators in boreal climates. In addition, hybrid cultivars of *Salix*, such as Klara (*Salix viminalis* × *S. schwerinii* × *S. dasyclados*) and/or Karin (*S. viminalis* × *S. schwerinii*) are in use in bioindication. Salam et al. [105] proved that Klara (*Salix viminalis* × *S. schwerinii* × *S. dasyclados*) accumulated high concentrations of Cu, Zn, and Ni in leaves, shoots, and roots from soils. At the same time, Klara may be considered for phytostabilisation and rhizofiltration due to the highest TM accumulation in the roots. *Phragmites* is an emergent plant worldwide that has demonstrated high TMs phytostabilisation properties [106,107]. In three reviewed papers, the authors focused on rare earth elements (lanthanide) such as Pr, Nd, Tb, Sm, Dy, and Er. Mikołajczyk et al. [108] have thoroughly documented the potential of *Artemisia vulgaris*, *Papaver rhoeas*, and *Taraxacum officinale* for phytoextraction of rare earth elements whose occurrence is increasing due to the development of modern technologies.


**Table 5.** A list of plant species tested as candidates for bioindication across the world.

Bioindicator use in biomonitoring programs may be an effective approach, although each time, it requires extended research to recognise how metals in a potential bioindicator correlate with the evaluated matrix (water, sediments, or soils). Species properties regarding metal uptake, the antagonistic and synergistic properties between the elements as well as whether the metal is essential or non-essential play a key role in evaluating bioindication technology [23]. Nevertheless, this is a useful assessment tool, offering the additional potential for recovering valuable elements, and so it is expected to develop relatively soon.

#### **5. Conclusions**

By applying geochemical and ecotoxicological assessment and classification indices, a numeric evaluation may be obtained regarding the extent of sediment contamination. An analysis of the results yielded by recent studies using the Web of Science Core Collection (based on five keywords: sediments, trace metals, contamination, pollution indices, urban) shows that the following occur most often: (1) the researchers studied metal concentration in the surface layer of sediments (0–2 cm, 0–5 cm, 0–10 cm); (2) the most frequently used geochemical indices were *Igeo*, *EF*, and *CF;* (3) ecotoxicological evaluation was mainly performed by using *SQG* and *PERI,* which includes the toxic response index. The abovepresented indices and classifications summarise the geochemical and ecotoxicological approaches of sediment assessment. The geoaccumulation indices help to determine whether the accumulation of TMs is due to natural processes or anthropogenic interference. Ecotoxicological assessment supports the overall contamination assessment of sediments and its harm for living organisms. The final results can be calculated for individual elements (*CF*, *Igeo*, *EF*, *Pi* (*SPI*), *PTT*) and by using complex pollution indices (*PLI*, *Cdeg*, *mCdeg*, *Pisum*, *PIAvg*, *PIaAvg*, *PIN*, *PIProd*, *PIapProd*, *PIvectorM*, *PINemerow*, *IntPI*, and *MPI*) for a larger set of TMs. Similar approaches can be applied for ecotoxicological indices (individual: *ER<sup>i</sup>* and complex: *CSI*, *SPI*, *RAC*, *PERI*, and *MERMQ*).

Geochemical indicators largely refer to the total metal content in the sample, although the inclusion of mobile and bioavailable fractions would appear to be indispensable for geochemical and ecotoxicological assessment. The use of *Igeointegrated* may help comprehensively verify the sediment quality. This approach includes contamination analysis, taking into account metal fractionation. A "good quality assessment" of sediments via *Igeo* may turn out to be worse after checking TM fractionation and referring to *Igeointegrated* (especially when a large proportion of TMs are bound in mobile or water-soluble fractions). Conversely, when TMs are mainly bound to an immobile and/or residual fraction, the *Igeo* assessment may be mendacious, resulting in an overestimation of the pollution level.

Bioindication is an indirect method of environmental assessment, which can support environmental quality assessment monitoring programs. Nowadays, the simultaneous use of different classifications and indices reflects a more accurate assessment of TM pollution. The most recognised bioindicators include *Phragmites* and *Salix* for TM phytoremediation, particularly Zn, Cd, Cu, and Ni in a boreal and temperate climate. Hybrid cultivars of *Salix* such as Klara (*Salix viminalis* × *S. schwerinii* × *S. dasyclados*) and Karin (*S. viminalis* × *S. schwerinii*) are also applied in biomonitoring and considered for phytostabilisation and rhizofiltration. *Vetiveria zizanoides* demonstrated resistance to As. Rare earth elements (lanthanides) such as Pr, Nd, Tb, Sm, Dy, and Er constitute a separate group of analysed metals investigated in plant research. This issue is of serious concern for the recovery of valuable elements in the development of a circular economy. *Artemisia vulgaris*, *Papaver rhoeas*, and *Taraxacum officinale* present lanthanide phytoextraction potential.

In order to be properly managed or remediated, sediments must be comprehensively evaluated. Monitoring programs for sediments in urban retention systems are crucial in this regard. A promising evaluation method is the use of bioindicators, which also have the potential to recover valuable elements. Therefore, future research should focus on developing this tool while evaluating the quality of sediments.

**Author Contributions:** Conceptualisation, N.N. and E.W.; methodology, N.N.; investigation, N.N., M.M., S.K.; resources, N.N.; data curation, N.N.; writing—original draft preparation, N.N. and E.W.; writing—review and editing, N.N. and E.W.; visualisation, N.N.; supervision, E.W., S.R., A.P.; funding acquisition, N.N. and E.W. All authors have read and agreed to the published version of the manuscript.

**Funding:** The authors acknowledge the support provided by InterPhD2 program funded by the European Union (Project No. POWR.03.02.00-IP.08–00-DOK/16).

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

#### **Appendix A**

**Table A1.** Comparison of trace elements investigated in bottom sediments of diverse urban water bodies (river, lakes, receivers, retention tanks, etc.) with concentration ranges in mg/kg d.w., indication of indices used for pollution evaluation, classifications, and verification of health risks and an indication of the use of metal fractionation.


**Table A1.** *Cont.*


**Table A1.** *Cont.*


**Table A1.** *Cont.*



**Table A1.** *Cont.*

Explanations: \* BDL—below detection limit; \*\* M—in monsoon; \*\*\* W—in winter; \*\*\*\* RCFG—rural commercial fishing group; \*\*\*\*\* RRG—rural reservoir group; \*\*\*\*\*\* UPG—urban park group; \*\*\*\*\*\*\* UCFG—urban commercial fishing group \* WL Wigry Lake in Poland is located in Wigry National Park (non urban area)—just to comparison with other studies.

#### **Appendix B**

**Table A2.** A list of Individual (single) Pollution Indices (IPIs) given in the literature with classification and a basic description of strengths and weaknesses.


**Table A2.** *Cont.*


**Table A2.** *Cont.*


\* Preindustrial concentration of TM according to Håkanson [58] [mg/kg]: Hg—0.25, Cd—1.0, As—15, Cu—50, Pb—70, Cr—90, Zn—175.

#### **Appendix C**

286

**Table A3.** A list of Complex Pollution Indices (CPIs) given in the literature with classification and basic description of strengths and weaknesses.


287

 **Index Formula Classification Description (Pros "+" and Cons "–") References**Degree of contamination (*Cdeg*) *<sup>C</sup>deg* = *n* ∑ *i*=1 *CFi CFi*Contamination Factor of i element *n*—the number of analysed TMs *<sup>C</sup>deg* < 8—low degree of contamination 8 ≤ *<sup>C</sup>deg* < 16—moderate degree of contamination 16 ≤ *<sup>C</sup>deg* < 32—considerable degree of contamination *<sup>C</sup>deg* ≥ 32—very high degree of contamination + allows for identifying the contamination in relation to several trace metals + easy to apply (easy to calculate) + precise scale + assesses a sum of contamination factors + used *CF* (includes *CmPre*) \* − not widely used − does not include geochemical background − the preindustrial reference value is necessary − omits the TMs availability and mobility in the environment − omits the natural variability process, grain size, and reference elements ratios [58] Modified contamination factor (*mCdeg*) *mCdeg* = 1*n* · *n*∑*i*=1 *CFi CFi*Contamination Factor of i element *n*—the number of analysed TMs *mCdeg* < 1.5—very low 1.5 ≤ *mCdeg* < 2—low 2 ≤ *mCdeg* < 4—moderate 4 ≤ *mCdeg* < 8—high 8 ≤ *mCdeg* < 16—very high 16 ≤ *mCdeg* < 32—extremely high *mCdeg* ≥ 32—ultra-high + allows for identifying the contamination in relation to several trace metals + easy to apply (easy to calculate) + precise scale + assesses a sum of contamination factors + used *CF* (includes *CmPre*) \* + widely used − does not include geochemical background − the preindustrial reference value is necessary − omits the TMs availability and mobility in the environment − omits the natural variability process, grain size, and reference elements ratios [64] Sum of Pollution Index (*PIsum*) *PIsum* = *n* ∑ *i*=1 *Pi Pi*—calculated value for Pollution Index *n*—the number of analysed TMs The classification for *Pi* can be used in *PIsum.* The values in *PIsum* should be multiplied by n (count of TMs): *PIsum* < <sup>1</sup>*n*—unpolluted, low level of pollution1*n*≤*PIsum*≤3*n*—moderatelypolluted+ defined as the sum of all determined contents of TMs, expressed as *Pi* + dedicated to combining all analysed TMs + based on *Pi* (includes *<sup>C</sup>mBackground*) − does not require the variation of natural processes − omits the TMs availability and mobility in theenvironment[63]

3*n* <

*PIsum*—heavily

 polluted the key is the choice of appropriate

does not include precise scale

−

− *<sup>C</sup>mBackground* value

**Table A3.***Cont.*

**Table A3.** *Cont.*


289


**Table A3.** *Cont.*

\* (includes *CmPre*)—with reference to the literature, the CPIs are related to the IPI, but it should be noted that in the IPI *CmPre* and *<sup>C</sup>mBackground* are two different concepts that are often confused/or taken interchangeably.

#### **Appendix D**


**Table A4.**A list of ecotoxicological indices (EIs) given in the literature with classification and basic description of strengths and weaknesses.

291

**Table A4.** *Cont.*


\* Toxic response factor for subsequent trace metals according to Håkanson [58]: Hg—40, Cd—30, As—10, Cu—5, Pb—5, Cr—2, Zn—1. \*\* Weight of each heavy metal according to Pejman et al. [86]: Cu—0.075, Zn—0.075, Cr—0.134, Ni—0.215, Pb—0.251, Cd—0.25.

### **References**


*Article*
