Next Article in Journal
Piezomagnetic Anomalies Associated with the 2021 MW 7.3 Maduo (China) Earthquake
Next Article in Special Issue
Investigation of Petroleum Hydrocarbon Fingerprints of Water and Sediment Samples of the Nestos River Estuary in Northern Greece
Previous Article in Journal
Modeling Flexural and Compressive Strengths Behaviour of Cement-Grouted Sands Modified with Water Reducer Polymer
Previous Article in Special Issue
Ecological Conditions of the Lower Dniester and Some Indicators for Assessment of the Hydropower Impact
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatiotemporal Variation and Ecological Risk Assessment of Heavy Metals in Industrialized Urban River Sediments: Fengshan River in Southern Taiwan as a Case Study

1
Institute of Environmental Engineering, National Sun Yat-sen University, Kaohsiung 80424, Taiwan
2
Department of Marine Environmental Engineering, National Kaohsiung University of Science andTechnology, Kaohsiung 81157, Taiwan
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2022, 12(3), 1013; https://doi.org/10.3390/app12031013
Submission received: 7 January 2022 / Revised: 17 January 2022 / Accepted: 17 January 2022 / Published: 19 January 2022
(This article belongs to the Special Issue Monitoring and Analysis of Environmental Pollution)

Abstract

:
The sediment pollution index acts as a useful indicator for assessing anthropogenic pollution within river drainage basins. An industrialized urban river, Fengshan River in Kaohsiung City, southern Taiwan has been suffering heavy metal pollution from surrounding factories. In this study, spatial and seasonal variations in heavy metals in sediments from seven sampling sites of Fengshan River were determined to assess sediment pollution status and potential ecological risk using multiple sediment pollution indices. Results showed that the heavy metal concentrations displayed large spatial variations. Severe contamination of heavy metals, especially for Cr, Hg, and Zn in the lower reaches of Fengshan River, may attribute to wastewater discharges from leather processing and metal finishing factories along the river drainage basin. An increase in metal concentrations from upstream to downstream indicated that heavy metals tend to accumulate in tidal reaches, probably as a result of the flocculation effect. Frequent heavy rainfall in the wet season can enhance surface runoff to discharge metal pollutants from non-point sources (scattered factories) into the river. Assessment of multiple pollution indices showed moderately polluted (mCd = 3.9, PLI = 2.6) and considerable ecological risk (RI = 540, mERMQ = 0.55), indicating Fengshan River sediments, particularly in the lower reaches, are considered toxic and can cause adverse effects to benthic organisms. Organic matters showed a good correlation with heavy metals, which play an important role in the spatiotemporal variations in heavy metal pollutants in the Fengshan River sediments. This study can provide valuable information for river pollution remediation, and urban planning and management.

1. Introduction

Heavy metal contamination in aquatic environments has attracted massive attention due to its toxicity, persistence, and accumulation, causing potential adverse effects on the ecosystems [1]. Those heavy metal pollutants such as cadmium (Cd), chromium (Cr), copper (Cu), mercury (Hg), and lead (Pb) can be discharged into a river system through anthropogenic sources through industrialization and distributed between aqueous and solid phases [2]. In river systems, suspended solids settling may cause the transfer of heavy metal pollutants from water to sediments [3]. The heavy metal pollutants in contaminated river habitats may enter into the food chain and pose a threat to human health [4]. Depending on the difference in particle size distribution, mineralogical composition, and organic matter content, river sediments may create anomalously heavy metals contamination depending on the river hydrological processes (e.g., weathering, erosion, and transport conditions) [5]. Moreover, heavy metal pollutants are not fixed permanently in the sediments, which may be desorbed/dissolved and released into the water column through the remobilization processes, negatively affecting water quality [6]. Thus, sediments act both as sinks and sources for heavy metals and strongly affect their distribution in river environments, playing an essential role in determining the pollution pattern of river systems [7,8]. Therefore, river sediments are valuable indicators for assessing anthropogenic pollution within river watersheds.
However, direct determination of heavy metal concentrations in sediments may not be suitable to assess the heavy metal pollution status and potential risks to exposed aquatic organisms; as heavy metals may be of natural or anthropogenic origin [9,10]. Therefore, the assessment of the sediment pollution level depends on sediment quality guidelines (SQGs) derived from numerous data of biological responses from aquatic organisms, including threshold-effects-level (TEL), probable-effects-level (PEL), effects-range-low (ERL), and effects-range-median (ERM) guidelines [11,12]. The SQGs can provide quick information on identifying the sediment quality and probability of biological toxicity. Based on the SQGs, many derived ecological and chemical approaches, such as mean ERM quotient (mERMQ) and the contamination severity index (CSI), have been used to evaluate the ecological risk of toxic metals for aquatic lives and the potential adverse effects on infaunal animals [13,14]. Furthermore, several contamination indices such as contamination factor (CF), geo-accumulation index (Igeo), and pollution load index (PLI) have also been widely used to distinguish anthropogenic pollutants from the natural background [15,16].
Therefore, whether it is a large natural river or a small river in an urban area, the heavy metal pollution level can be assessed using the sediment quality guidelines and the sediment pollution indices. Even though the actual ecological impacts of river sediment may vary with the water environment condition and sediment characteristics, which need to be confirmed by laboratory or on-site toxicity experiments, those sediment pollution indices can provide a simple and integrated judgement of river contamination levels and potential ecological risks of heavy metal pollution.

Case Study

Kaohsiung City is an important industrial city in Taiwan with the third-largest population of over 2.7 million. So far, over 55% of domestic wastewaters in Kaohsiung City are still directly discharged to rivers without treatment due to the uncompleted construction and household connection of sewage systems. Moreover, a lot of industries and factories, such as electronic industries, chemical manufacturing factories, metal processing plants, foundries, paint and dye industries, motor vehicle plating and finishing plants, and paper and board mills, are located in or around Kaohsiung City [13]. Those industries may discharge a considerable amount of untreated and minimally treated wastewaters into the adjacent rivers. Fengshan River, one of the most important rivers flowing through the downtown of Kaohsiung City, receives domestic, industrial, agriculture, and husbandry wastewater from the residential areas, factories, and stock farming in the watershed. In terms of the total amount of pollution, industrial wastewater accounts for the most contribution. The characteristics of heavy metal contamination strongly depend on the industries along the river. The situation leads to the need to evaluate the heavy metal contamination in the Fengshan River sediments and assess the potential environmental risk. Despite this great need, no studies have studied the heavy metal contamination in Fengshan River sediments and its associated ecological impacts.
Therefore, this study investigated spatial and temporal variations in heavy metals in the surface sediments of Fengshan River to assess the degree of heavy metal contamination and the potential ecological risk using multiple sediment pollution indices. The relationship between sediment properties and heavy metals, and the major factors controlling the spatiotemporal variations in sediment heavy metal concentrations were also determined.

2. Materials and Methods

This study focused on an industrialized urban river, Fengshan River, in Kaohsiung City southern Taiwan, and assessed the river sediment pollution status using a multivariate approach. The sediment pollution indices applied in this study includes the contamination factor (CF), geo-accumulation index (Igeo), modified degree of contamination (mCd), pollution load index (PLI), potential ecological risk index (RI), mean ERM quotient (mERMQ), sum of the toxic units (ΣTU), and contamination severity index (CSI), which were widely used in sediment heavy metal pollution assessment [13,14,15,16,17].

2.1. Study Area

Fengshan River originates from the Jiuqutang Mountain area and flows through Dashu, Daliao, Niaosong, and Fengshan District of Kaohsiung City. As the river flows southwest into Qianzhen District, it is renamed Qianzhen River and finally injected into Kaohsiung Port. Fengshan River (including the Qianzhen River) is approximately 20 km, with a drainage area of ~53.85 km2. The primary pollutant sources in the Fengshan River drainage basin include 43% of industrial wastewater, 40% of domestic sewages, 15% of husbandry wastewater, and 2% agricultural wastewater [18]. Industrial wastewater can be regarded as the primary source of pollution in the river due to the distribution of many small factories along and within the river basin, such as metal finishing, basic metal manufacturing, electroplating, printed circuit board manufacturing, and leather processing, especially the pollution from leather processing and metal finishing factories in the upstream area [19]. In the past, while the pollution prevention and control legislation were not complete and the sewage system was not popularized, most of the wastewater from leather processing was untreated or limitedly treated and discharged directly into Fengshan River, extremely polluting the river and causing environmental risk.

2.2. Sample Collection and Analysis

In this study, surface sediment samples (0–15 cm) from seven sampling sites (Figure 1) were collected from the main channel and river mouth of Fengshan River using a Van Veen dredge in the dry season (December 2015) and the wet season (May 2016). The sampling sites are situated upstream (F1, F2), midstream (F3, F4, F5), and downstream (F6, F7). Among them, site F3 is located nearby the effluent outfall of the Fengshan Water Resources Center. Sites F6 and F7 are situated at the tidal river reach, of which site F7 is located at the river mouth. Sites F1–F6 were sampled at the center of the river channel on bridge, and site F7 was sampled on a boat. According to Taiwan Central Weather Bureau [20], the accumulated precipitation in Kaohsiung City were 94 mm and 359 mm on October–December 2015 (dry season) and March–May 2016 (wet season), respectively.
After the sediment samples were air-dried, ground, and homogenized, approximately 1 g of dry sediments was digested in a Teflon digestion vessel with nitric acid (HNO3), hydrogen peroxide (H2O2), and hydrochloric acid (HCl) by referring to protocols of the standard method of the Taiwan Environmental Protection Administration (NIEA M353.02C). Firstly, the sediment sample was digested with 15 mL HNO3 at 95 ± 5 °C (not boiling) for 20 min. Then, the sediment sample was digested with 10 mL H2O2 at 95 ± 5 °C (not boiling) for 10 min. Finally, the sediment sample was digested with 10 mL HCl at 95 ± 5 °C (not boiling) for 15 min. After digestion, the digested solution was diluted to 50 mL with deionized water (~18 MΩcm). The concentrations of Cr, Cu, Ni, Pb, Zn, As, and Hg in the digested solution were analyzed by an inductively coupled plasma atomic emission spectrometer (Thermo Scientific iCAP 7400, Waltham, MA, USA) based on the standard method of the Taiwan Environmental Protection Administration (NIEA M104.02C). The standard reference material of sediment (NIST SRM 1646a, Gaithersburg, MD, USA) was analyzed using the same digestion and measurement procedures to validate the results.

2.3. Assessment of Sediment Contamination Level and Ecological Risk

2.3.1. Contamination Factor (CF)

The contamination factor (CF) was calculated using the following equation:
CF = M s a m p l e M b a c k g r o u n d
where Msample and Mbackground are heavy metal concentrations of the studied sediment sample and the background, respectively. The Mbackground was acquired from the continental upper crust average, i.e., 92, 28, 47, 17, 67, 4.8, and 0.05 mg/kg for Cr, Cu, Ni, Pb, Zn, As, and Hg, respectively [21]. Based on the CF value, the sediment heavy metal contamination status can be interpreted as very high (CF > 6), considerable (3–6), moderate (1–3), and low contamination (<1) [22].

2.3.2. Geo-Accumulation Index (Igeo)

The geo-accumulation index (Igeo) was calculated using the following equation:
I geo = log 2 M s a m p l e 1.5 M b a c k g r o u n d
where Msample and Mbackground are heavy metal concentrations of the studied sediment sample and the background, respectively. Based on the Igeo value, the sediment heavy metal contamination status can be classified as extremely polluted (Igeo > 5), strongly to extremely polluted (4–5), strongly polluted (3–4), moderately to strongly polluted (2–3), moderately polluted (1–2), unpolluted to moderately polluted (0–1), and unpolluted (<0) [23].

2.3.3. Modified Degree of Contamination (mCd)

The modified degree of contamination (mCd) was calculated using the following equation:
mCd = i = 1 n C F i n
where CFi and n are the contamination factor of each heavy metal obtained using the Equation (2) and the number of heavy metals studied, respectively. Based on the mCd value, the sediment heavy metal contamination level can be categorized into ultra-high degree (mCd > 32), extremely high degree (16–32), very high degree (8–16), high degree (4–8), moderate degree (2–4), low degree (1.5–2), and nil to very low degree of contamination (<1.5) [24].

2.3.4. Pollution Load Index (PLI)

The pollution load index (PLI) was calculated using the following equation:
PLI = C F 1 × C F 2 × C F 3 × × C F n n
where CFi and n are the contamination factor of each heavy metal, and the number of heavy metals studied, respectively. The sediment heavy metal contamination level can be ranked based on the PLI value and classified as strongly polluted (>3), moderately (2–3), and slightly (1–2). The simple interpretation of PLI value was no pollution for PLI < 1 and polluted for PLI > 1 [25].

2.3.5. Potential Ecological Risk Index (RI)

The potential ecological risk index (RI) was calculated using the following equation:
RI = i = 1 n ( C F i × T r i )
where CFi and Tri are the contamination factor and biological toxicity factor of each heavy metal, respectively. The Tri value was suggested by Håkanson [22] as 2, 5, 5, 5, 1, 10, and 40 for Cr, Cu, Ni, Pb, Zn, As, and Hg, respectively. Based on the RI value, the ecological risk of the studied sediment can be classified into very high (RI > 600), considerable (300–600), moderate (150–300), and low ecological risk (<150) [22].

2.3.6. Mean ERM Quotient (mERMQ)

The mean ERM quotient (mERMQ) was calculated using the following equation:
mERM Q = i = 1 n M s a m p l e / M E R M i n
where Msample and MERM are heavy metal concentrations of the studied sediment sample and the effect-range-median (ERM) criteria value, respectively, and n is the number of heavy metals studied. The MERM are 145, 390, 50, 110, 270, 85, and 1.3 mg/kg for Cr, Cu, Ni, Pb, Zn, As, and Hg, respectively [26]. Based on the mERMQ value, the toxicity probability of the studied sediment for benthic organisms can be categorized into high toxic (76% probability of being toxic for mERMQ > 1.5), considerable toxic (49%, 0.5–1.5), moderate toxic (21%, 0.1–0.5), and slightly toxic (9%, <0.1) [27].

2.3.7. Sum of the Toxic Units (ΣTU)

The sum of the toxic units (ΣTU) was calculated using the following equation:
Σ TU = i = 1 n M s a m p l e M P E L i
where Msample and MPEL are heavy metal concentrations of the studied sediment sample, and the probable-effects-level (PEL) criteria value, respectively. According to MacDonald et al. [26], the MPEL are 90, 197, 36, 91.3, 315, 17, and 0.486 mg/kg for Cr, Cu, Ni, Pb, Zn, As, and Hg, respectively. The potential acute toxicity of the studied sediment for benthic organisms can be ranked based on the ΣTU value, of which high ΣTU means high acute toxicity and vice versa [28].

2.3.8. Contamination Severity Index (CSI)

The contamination severity index (CSI) was calculated using the following equation:
CSI = i = 1 n W i M s a m p l e M E R L 1 / 2 + M s a m p l e M E R M 2 i
where Msample, MERL, and MERM are heavy metal concentrations of the studied sediment sample, the effect-range-low (ERL) and the effect-range-median (ERM) criteria values, respectively. The Wi value was calculated by the eigenvalue and loading value obtained from principal component analysis/factor analysis (PCA/FA) based on the following Equation (9), which the only factor attributed to anthropogenic pollution source was computed [29].
W i = (   Loading   Value i × Eigenvalue   ) i n   (   Loading   Value i × Eigenvalue   )
Based on the CSI value, the severity levels of ecological risk for sediment contamination can be classified into ultra-high (CSI > 5), very high (4–5), high (3–4), moderate to high (2.5–3), moderate (2–2.5), low to moderate (1.5–2), low (1–1.5), and very low severity of contamination (0.5–1), respectively [29].

2.4. Multivariate Statistical Analyses

The multivariate statistical analyses, including the Spearman correlation coefficient matrix and PCA/FA, were performed on the whole data set (13 variables) using IBM SPSS Statistics software version 20. The grain size composition (clay, silt, sand), organic matter (OM), total nitrogen (TN), and total phosphorus (TP) contents of the Fengshan River sediments in same study sites were obtained from our previous study [18].

3. Results and Discussion

3.1. Spatial and Seasonal Variations in Heavy Metals

Table 1 shows the concentrations of heavy metals in Fengshan River sediments. The spatial and seasonal variations in trace metals concentrations in the studied sediments are shown in Figure 2. The distribution of heavy metals in the sediments showed large variations with a range of 15.8–352.7, 17.3–175.3, 17.3–64.3, 11.6–72.8, 119–533, 11.2–44.1, and 0.18–1.44 mg/kg for Cr, Cu, Ni, Pb, Zn, As, and Hg, respectively. The average concentrations of those heavy metals (Cr, Cu, Ni, Pb, Zn, As, and Hg) were 118.2, 78.4, 39.1, 33.0, 277.3, 24.9, and 0.6 mg/kg, respectively. In particular, the concentrations of Cr, Cu, Ni, Pb, and Zn displayed a similar trend of spatial variations. During the dry season, Cr, Cu, Ni, Pb, and Zn all showed the highest concentrations at site F6, which were 294.9, 142.6, 57.1, 49.9, and 501 mg/kg, respectively. The second-highest concentration was observed at the river mouth (F7), which was 234.6, 117.7, 39.1, 33.5, and 343 mg/kg for Cr, Cu, Ni, Pb, and Zn, respectively. Those metal concentrations at other stations (F1–F4) were relatively low and gradually increased from upstream to downstream during the dry season. During the wet season, the highest concentration was found at site F6, which was 352.7, 175.3, 64.3, 72.8, 533 mg/kg for Cr, Cu, Ni, Pb, and Zn, respectively. The second-highest concentration was observed at the river mouth (F7) except for Ni, which was 292.6, 148.0, 49.5, and 431 mg/kg for Cr, Cu, Pb, and Zn, respectively. Similarly, other stations also showed an increasing trend from upstream to downstream during the wet season. The mercury concentration in the dry season rose from the upstream to the maximum of 1.44 mg/kg at site F5 and then decreased toward site F7, whereas the highest concentration was found at site F6 of 1.30 mg/kg in the wet season. There was no significant variation in the concentration of arsenic along the river in both dry and wet seasons.
Overall, the spatial distribution of heavy metals in the Fengshan River sediments showed a relatively high concentration in the lower reaches (F6, F7). The maximum was frequently found at site F6 where the boundary of tidal river reach is located, probably as a result of the flocculation effect [30,31]. The salinity of river water serves as an important controlling factor for the partitioning of heavy metals to sediments, particularly the flocculation process in the lower reaches and estuary of rivers affected the distribution of heavy metals in the river sediments [31]. The concentration at the river mouth (F7) slightly decreased compared with site F6, probably because of resuspension and transportation of contaminated sediments seaward.
The seasonal variation in precipitation plays an important role in the fate of anthropogenic pollutants in the aquatic environment, especially for the river basin with distinct wet and dry seasons, such as the Fengshan River in Kaohsiung City, Taiwan. As shown in Figure 2, most heavy metals showed relatively high concentrations in the sediments during the wet season compared with those collected during the dry season. During the wet season, frequent and heavy rainfall can induce the discharge of pollutants from non-point sources to the Fengshan River through surface runoff. Moreover, the river sediment in the upper reaches may be resuspended by a strong river stream when a heavy rainfall occurs, causing sediments to scour from upstream to downstream and estuaries. Some studies have also shown similar findings that the higher concentration of heavy metals in river sediments during the rainy season resulted from heavy rainfall-induced water erosion and runoff of heavy metals-containing materials from the upper catchment [32,33].
The comparison of the heavy metals (Cr, Cu, Ni, Pb, Zn, As, and Hg) concentrations in Fengshan River sediments and the sediment quality guidelines are shown in Figure 2 and Table 1. The heavy metal concentrations in the studied sediments were mostly higher than the criteria value of ERL. As exceeded the ERL at two sites (F2 and F4), Pb exceeded the ERL at three sites (F4, F6, and F7), both Cr and Cu exceeded the ERL at four sites (F2, F5, F6, and F7 for Cr, and F2, F4, F6, and F7 for Cu), Ni exceed the ERL at six sites (except F1), and both Hg and Zn exceeded the ERL at all the sampling sites. Moreover, Hg, Cr, Ni, Zn showed concentrations exceeding the criteria value of ERM, where Hg exceeded the ERM at one site (F5), Cr exceeded the ERM at two sites (F6 and F7), and both Ni and Zn exceeded the ERM at three sites (F4, F5, F6 for Ni and F4, F6, F7 for Zn). These findings indicated that the concentrations of Hg, Cr, Ni, and Zn were severely contaminated in Fengshan River sediments, leading to the frequent occurrence of adverse biological effects. Those kinds of heavy metals pollution may attribute to wastewater discharges from leather processing and metal finishing factories within the river watershed, particularly along the river upstream and tributaries [19,34].

3.2. Sediment Pollution Status

The sediment heavy metal contamination levels in this study were assessed using the contamination factor (CF), geo-accumulation index (Igeo), modified degree of contamination (mCd), and pollution load index (PLI). The CF and Igeo are often applied to differentiate the metal source between natural terrestrial and anthropogenic pollution in sediments. Both the mCd and the PLI were comparative indices for assessing overall heavy metals contamination level in a simple way.
As shown in Figure 3a, As, Hg, Zn, and Cu have CF values greater than 6, which denotes “very high contamination”. The CF value for Hg shows the highest among all metals and reveals the highest number of sites with CF > 6 among all studied metals. The second-highest occurrence of “very high contamination” is As, which appeared at the F2 and F4 both in the dry and wet seasons. Both Zn and Cu have CF > 6 at site F6, which appeared in the wet season for Zn and in each of the wet and dry seasons for Cu. The CF value of “considerable contamination” mainly came from As, Hg, Zn, and Cu. Other than the sites showing “very high contamination”, the CF value for Hg at other sites showed “considered contamination”. The CF value of “considered contamination” for As were shown at sites F1, F3, and F5 both in the dry and wet seasons, while those CF values for Zn were shown at sites F3–F5 and F7 in the dry season and sites F4–F5 in the wet season, and those for Cu were shown at site F7 in the dry season and site F4 and F7 in the wet season. In addition, the CF value with “considered contamination” for Cr also had one site (F6) in the dry season and two sites (F6 and F7) in the wet season, and those for Pb only appeared at site F7 in the wet season. In contrast, none of the sites showed “considered contamination” in Ni. Overall, the sediment contamination degree of heavy metals is ranked in the following order of Hg > As > Zn > Cu > Pb > Cr > Ni.
The geo-accumulation index (Igeo) results are shown in Figure 3b. The Igeo classes of the Fengshan River sediments were ranked from “moderately polluted” to “strongly to extremely polluted” for Hg, “unpolluted to moderately polluted” to “moderately to strongly polluted” for As and Zn, “unpolluted” to “moderately to strongly polluted” for Cu, “unpolluted” to “moderately polluted” for Pb and Cr, and “unpolluted” for Ni. Similar to the result of contamination factors (CF), the contamination degree of heavy metal in the Fengshan River based on the Igeo value followed the order of Hg > As > Zn > Cu > Pb > Cr > Ni.
The mCd is an overall assessment of the degree of contamination for all studied metals. According to the mCd assessment (Figure 3c), the mCd values classified as “high degree contamination” were shown at sites F4 and F6 in both the dry and wet seasons, at site F5 in the dry season and site F7 in the wet season. Both sites F2 and F3 in the dry and wet seasons, site F5 in the wet season, and site F7 in the dry season showed “moderate degree contamination”. In contrast, site F1 only indicated “low degree contamination” in both dry and wet seasons.
The PLI values ranged from 1.1 to 5.07 (Figure 3d), indicating that heavy metal pollution existed (PLI > 1) in Fengshan River sediments. Most sites in the Fengshan River showed “moderate degree polluted” to “strongly polluted” except sites F1 and F2 were “slightly polluted”, which further indicated that the industrial effluent discharge might be responsible for the accumulation of heavy metals in sediments of the Fengshan River. Overall, the degree of heavy metal contamination in Fengshan River can be classified as “moderate degree polluted,” mainly attributed to the industrial effluent discharge. The degree of heavy metal contamination showed an increasing trend from the upper reaches to the lower reaches. It slightly decreased in the river mouth, probably because of resuspension and transportation of contaminated sediments seaward.
A comparison between the heavy metal concentrations of sediments in Fengshan River and those from other regions is shown in Table 2. The concentrations of heavy metals in the Fengshan River sediments were mostly higher than those of river sediments from other countries, especially for Hg, Cr, Cu, and Zn. It highlighted the seriousness of heavy metal pollution in the Fengshan drainage area, which needs to pay more attention to this critical issue by the government and related organizations.

3.3. Potential Biological Risk

The potential biological risk of the Fengshan River sediments for benthic organisms was assessed by the potential ecological index (RI), mean ERM quotient (mERMQ), the sum of the toxic units (ΣTU), and contamination severity index (CSI). As shown in Figure 4a, the sediments in Fengshan River exhibited RI values ranging between 199 and 1233 with an average of 542, showing a considerable ecological risk according to the classification reported by Håkanson [22]. In particular, sites F4, F5, and F6 showed “very high ecological risk” with an average RI value up to 670. The mERMQ values varied from 0.24 to 1.14 with an average of 0.55, indicating that the contamination level of heavy metals in Fengshan River sediments may be considerably toxic (49% probability of being toxic) to benthic organisms (Figure 4b). A relatively high mERMQ value was found at site F6 (mERMQ = 1.14 in the wet season and mERMQ = 0.97 in the dry season), followed by site F7 (mERMQ = 0.82 in the wet season and mERMQ = 0.67 in the dry season). It is indicated that the sediments in the lower reaches and the river month were highly toxic (76% probability of being toxic), which may be harmful to the growth of local benthic organisms. Similarly, the ΣTU values ranged between 3.2 and 12.6, with an average of 6.7. The highest ΣTU value was observed at site F6 (ΣTU = 12.6 in the wet season and ΣTU = 10.5 in the dry season), indicating a high risk of potential acute toxicity effects of the river sediments on the benthic community, especially for those in the lower reaches (Figure 4c).
In addition, the CSI value was calculated by the eigenvalues of 6.387 and loading values of 0.776, 0.910, 0.881, 0.915, 0.819, 0.175, and 0.711 for Cr, Cu, Ni, Pb, Zn, As, and Hg, respectively. The eigenvalues and loading values adopted the VF1 obtained from PCA/FA, which can be interpreted as anthropogenic sources. As shown in Figure 4d, the CSI values ranged between 0.86 and 3.88, with a mean of 1.8. A similar trend to the mean ERM quotient was observed, which showed the highest value at site F6 (CSI = 3.9 in the wet season and CSI = 3.2 in the dry season), followed by site F7 (CSI = 2.7 in the wet season and CSI = 2.2 in the dry season). This indicated moderate to high severity of heavy metal contamination in the lower reaches and the river mouth.

3.4. Relationship among Sediment Properties and Heavy Metals

As shown in Table 3, significantly positive relationships between heavy metals (Cr, Cu, Ni, Pb) and organic matters (OMs) were observed in Fengshan River sediments, showing a correlation coefficient of 0.56–0.65 at a confidence level of 95%. It is indicated that OM is an important factor affecting the sediment heavy metal distribution in Fengshan River. Moreover, the sediment metal concentrations were also well correlated with total nitrogen (TN) and total phosphorus (TP) contents. Due to TN and TP being mainly derived from organic compounds, a good correlation of heavy metals with TN and TP strongly indicated that the organic compounds play an essential role in the spatiotemporal variations in heavy metals in the Fengshan River sediments. However, relatively low correlation coefficient values were observed for heavy metals with grain size (either sand, silt, or clay), indicating that the influence of grain size composition on the sediment heavy metal distribution was less critical than the organic matter.
In general, grain size and organic matter serve as two primary factors to affect the variations in heavy metal concentrations in sediments. The grain size effect has been documented in previous studies [45,46]. Muddy sediments, mainly consisting of aluminosilicate minerals (i.e., clay minerals), have high cation exchange capacity and large specific surface area [47], which are prone to adsorb anthropogenic pollutants (e.g., heavy metals). Consequently, heavy metals are accumulated in the muddy sediments [48]. On the contrary, sandy sediments commonly consist of coarse-grained biogenic carbonates and quartz. The carbonate or quartz sands with relatively low cation exchange capacity and specific surface area are not easy to adsorb heavy metal pollutants compared with clay minerals [49].
On the other hand, high concentrations of heavy metals are often associated with high organic matter content in sediments [17,50]. Most organic matters have large specific surface areas with negatively-charged sites, leading to the adsorption of anthropogenic heavy metals [51]. Chen et al. [17] reported that the influence of organic matter content on the distribution of heavy metals in sediments was more significant than the effect of sediment grain size, especially for the aquatic environments with rich organic matter (e.g., estuaries and harbors).
Furthermore, good relationships among the heavy metals (Hg, Cr, Cu, Ni, Pb, and Zn) in the studied sediments were observed with a correlation coefficient ranging between 0.53 and 0.94 at a confidence level of 95% (p < 0.05). This indicated that those heavy metals were probably derived from similar pollution sources, which may attribute to industrial wastewater discharges from leather processing and metal finishing factories within the river basin.
According to the results of PCA/FA, three factors with eigenvalues > 1 were extracted from our data set. The corresponding VFs, variable loadings, eigenvalues, and the explained variance are shown in Table 4. The VF1 displays strong positive loadings (>0.70) on OM, TN, TP, Hg, Cr, Cu, Ni, Pb, and Zn, which accounts for 49% of the total variance. This is consistent with the results of the Spearman correlation analysis, where good correlation among OM, TN, TP, and heavy metals (Cr, Cu, Ni, Pb) were observed (Table 2). Therefore, the VF1 can be interpreted as anthropogenic pollution sources. In contrast, the VF2 displays strong positive loadings on grain size composition (i.e., clay, silt, and sand) with 24.1% of the total variance, which can be interpreted as lithogenic sources (i.e., natural origin). The difference of PCA/FA factors obtained from grain size and heavy metals indicated grain size appears to be not a major controlling factor for the heavy metals’ spatiotemporal variations in the river sediments. This is consistent with the results of the Spearman correlation analysis. The VF3 displays strong positive loadings on As only with 13.6% of the total variance, probably indicating the source of As was different from the other metals (Hg, Cr, Cu, Ni, Pb, and Zn).

4. Conclusions

The results of this study clearly showed that Fengshan River sediments were extremely contaminated with Hg, Cr, Ni, and Zn. Those heavy metal pollutants were mainly attributed to the leather processing and metal finishing factories around the river upstream. The contaminated river sediments pose a considerable ecological risk based on the assessment of multiple pollution indices, especially for the sediments in the lower tidal reaches and river mouth, which may cause adverse effects on the benthic organisms. Even though the ecological impacts of river sediments on the local benthic organisms may depend on the water environment condition and sediment characteristics, the results of this study can provide a warning to the government, who should pay attention to river remediation and improve the current pollution status. Further laboratory and on-site toxicity experiments are needed to be conducted to understand the actual adverse effects on local organism in Fengshan River. This study can provide valuable information for river pollution remediation and urban planning and management.

Author Contributions

Conceptualization, C.-D.D. and C.-M.K.; methodology, C.-F.C.; validation, K.-N.L. and Y.-C.L.; formal analysis, K.-N.L.; investigation, Y.-C.L.; data curation, C.-W.C.; writing—original draft preparation, K.-N.L.; writing—review and editing, C.-W.C., C.-M.K. and C.-D.D.; visualization, K.-N.L.; supervision, C.-M.K.; project administration, C.-M.K.; funding acquisition, C.-D.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Lin, Q.; Liu, E.; Zhang, E.; Li, K.; Shen, J. Spatial distribution, contamination and ecological risk assessment of heavy metals in surface sediments of Erhai Lake, a large eutrophic plateau lake in southwest China. Catena 2016, 145, 193–203. [Google Scholar] [CrossRef]
  2. Sin, S.N.; Chua, H.; Lo, W.; Ng, L.M. Assessment of heavy metal cations in sediments of Shing Mun River, Hong Kong. Environ. Int. 2001, 26, 297–301. [Google Scholar] [CrossRef]
  3. Ghrefat, H.; Yusuf, N. Assessing Mn, Fe, Cu, Zn, and Cd pollution in bottom sediments of Wadi Al-Arab Dam, Jordan. Chemosphere 2006, 65, 2114–2121. [Google Scholar] [CrossRef] [PubMed]
  4. Deniseger, J.; Erickson, L.J.; Austin, A.; Roch, M.; Clark, M.J.R. The effects of decreasing heavy metal concentrations on the biota of Buttle Lake, Vancouver Island, British Columbia. Water Res. 1990, 24, 403–416. [Google Scholar] [CrossRef]
  5. Reimann, C.; de Caritat, P. Distinguishing between natural and anthropogenic sources for elements in the environment: Regional geochemical surveys versus enrichment factors. Sci. Total Environ. 2005, 337, 91–107. [Google Scholar] [CrossRef]
  6. Simpson, S.L.; Spadaro, D.A. Bioavailability and chronic toxicity of metal sulfide minerals to benthic marine invertebrates: Implications for deep sea exploration, mining and tailings disposal. Environ. Sci. Technol. 2016, 50, 4061–4070. [Google Scholar] [CrossRef]
  7. Maanan, M.; Saddik, M.; Maanan, M.; Chaibi, M.; Assobhei, O.; Zourarah, B. Environmental and ecological risk assessment of heavy metals in sediments of Nador lagoon, Morocco. Ecol. Indic. 2015, 48, 616–626. [Google Scholar] [CrossRef]
  8. Singh, K.P.; Mohan, D.; Singh, V.K.; Malik, A. Studies on distribution and fractionation of heavy metals in Gomti river sediments-a tributary of the Ganges. J. Hydrol. 2005, 312, 14–27. [Google Scholar] [CrossRef]
  9. Chapman, P.M.; Wang, F.; Janssen, C.; Persoone, G.; Allen, H.E. Ecotoxicology of metals in aquatic sediments: Binding and release, bioavailability, risk assessment, and remediation. Can. J. Fish. Aquat. Sci. 1998, 55, 2221–2243. [Google Scholar] [CrossRef]
  10. McCauley, D.J.; DeGraeve, G.M.; Linton, T.K. Sediment quality guidelines and assessment: Overview and research needs. Environ. Sci. Policy 2000, 3, S133–S144. [Google Scholar] [CrossRef]
  11. Long, E.R.; MacDonald, D.D.; Smith, S.L.; Calder, F.D. Incidence of adverse biological effects within ranges of chemical concentrations in marine and estuarine sediments. Environ. Manag. 1995, 19, 81–97. [Google Scholar] [CrossRef]
  12. MacDonald, D.D.; Carr, R.S.; Calder, F.D.; Long, E.R.; Ingersoll, C.G. Development and evaluation of sediment quality guidelines for Florida coastal waters. Ecotoxicology 1996, 5, 253–278. [Google Scholar] [CrossRef]
  13. Chen, C.W.; Lim, Y.C.; Chen, C.F.; Ju, Y.R.; Dong, C.D. Spatial distribution and ecological risk assessment of sediment metals in a highly industrialized coastal zone southwestern Taiwan. Environ. Sci. Pollut. Res. 2019, 26, 14717–14731. [Google Scholar] [CrossRef]
  14. Benson, N.U.; Adedapo, A.E.; Fred-Ahmadu, O.H.; Williams, A.B.; Udosen, E.D.; Ayejuyo, O.O.; Olajire, A.A. New ecological risk indices for evaluating heavy metals contamination in aquatic sediment: A case study of the Gulf of Guinea. Reg. Stud. Mar. Sci. 2018, 18, 44–56. [Google Scholar] [CrossRef]
  15. Kim, B.S.M.; Angeli, J.L.F.; Ferreira, P.A.L.; de Mahiques, M.M.; Figueira, R.C.L. A multivariate approach and sediment quality index evaluation applied to Baixada Santista, Southeastern Brazil. Mar. Pollut. Bull. 2019, 143, 72–80. [Google Scholar] [CrossRef] [PubMed]
  16. Ali, A.E.; Strezov, V.; Davies, P.J.; Wright, I. River sediment quality assessment using sediment quality indices for the Sydney basin, Australia affected by coal and coal seam gas mining. Sci. Total Environ. 2018, 616, 695–702. [Google Scholar] [CrossRef] [PubMed]
  17. Bhuiyan, M.A.H.; Dampare, S.B.; Islam, M.A.; Suzuki, S. Source apportionment and pollution evaluation of heavy metals in water and sediments of Buriganga River, Bangladesh, using multivariate analysis and pollution evaluation indices. Environ. Monit. Assess. 2015, 187, 4075. [Google Scholar] [CrossRef] [PubMed]
  18. Lin, K.N.; Chen, C.W.; Chen, C.F.; Lim, Y.C.; Kao, C.M.; Dong, C.D. Seasonal variation of phthalate esters in urban river sediments: A case study of Fengshan River system in Taiwan. Sustainability 2022, 14, 347. [Google Scholar] [CrossRef]
  19. Zou, E.T. Implementation of Environmental Policy: A Case Study of Water Pollution Control of Kaohsiung County’s Fengshan River. Master’s Thesis, I-Shou University, Kaohsiung, Taiwan, 2005. (In Chinese). [Google Scholar]
  20. Taiwan Central Weather Bureau. 2021. Available online: https://www.cwb.gov.tw/eng/ (accessed on 4 January 2022).
  21. Rudnick, R.L.; Gao, S. Composition of the continental crust. In The Crust, Treatise on Geochemistry; Rudnick, R.L., Ed.; Elsevier: Amsterdam, The Netherlands, 2003; Volume 3, pp. 1–64. [Google Scholar] [CrossRef]
  22. Håkanson, L. An ecological risk index for aquatic pollution control: A sedimentological approach. Water Res. 1980, 14, 975–1001. [Google Scholar] [CrossRef]
  23. Müller, G. Schwermetalle in den Sedimenten des Rheins-VeraÈnderungenseit. Umschau 1979, 79, 778–783. [Google Scholar]
  24. Abrahim, G.; Parker, R. Assessment of heavy metal enrichment factors and the degree of contamination in marine sediments from Tamaki Estuary, Auckland, New Zealand. Environ. Monit. Assess. 2008, 136, 227–238. [Google Scholar] [CrossRef]
  25. Tomlinson, D.; Wilson, J.; Harris, C.; Jeffrey, D. Problems in the assessment of heavy-metal levels in estuaries and the formation of a pollution index. Helgol. Meeresunters. 1980, 33, 566–575. [Google Scholar] [CrossRef] [Green Version]
  26. MacDonald, D.D.; Ingersoll, C.G.; Berger, T.A. Development and evaluation of consensus-based sediment quality guidelines for freshwater ecosystems. Arch. Environ. Contam. Toxicol. 2000, 39, 20–31. [Google Scholar] [CrossRef]
  27. Long, E.R.; MacDonald, D.D.; Severn, C.G.; Hong, C.B. Classifying probabilities of acute toxicity in marine sediments with empirically derived sediment quality guidelines. Environ. Toxicol. Chem. 2009, 19, 2598–2601. [Google Scholar] [CrossRef]
  28. Pedersen, F.; Bjørnestad, E.; Anderson, H.V.; Kjølholt, J.; Poll, C. Characterization of sediments from Copenhagen Harbour by use of biotests. Water Sci. Technol. 1998, 37, 233–240. [Google Scholar] [CrossRef]
  29. Pejman, A.; Bidhendi, G.N.; Ardestani, M.; Saeedi, M.; Baghvand, A. A new index for assessing heavy metals contamination in sediment: A case study. Ecol. Indic. 2015, 58, 365–373. [Google Scholar] [CrossRef]
  30. George, D.A.; Hill, P.S.; Milligan, T.G. Flocculation, heavy metals (Cu, Pb, Zn) and the sand–mud transition on the Adriatic continental shelf, Italy. Cont. Shelf Res. 2007, 27, 475–488. [Google Scholar] [CrossRef]
  31. Schindler, R.J.; Comber, S.D.W.; Manning, A.J. Metal pollutant pathways in cohesive coastal catchments: Influence of flocculation and biopolymers on partitioning and flux. Sci. Total Environ. 2021, 795, 148800. [Google Scholar] [CrossRef] [PubMed]
  32. Karthikeyan, R.; Vijayalakshmi, S.; Balasubramanian, T. Seasonal distribution of heavy metals in the sediments from Uppanar estuary (East coast of India). Int. J. Aquat. Biol. 2004, 19, 119–122. [Google Scholar]
  33. Nguyen, B.T.; Do, D.D.; Nguyen, T.X.; Nguyen, V.N.; Nguyen, D.T.P.; Nguyen, M.H.; Truong, H.T.T.; Dong, H.P.; Le, A.H.; Bach, Q.V. Seasonal, spatial variation, and pollution sources of heavy metals in the sediment of the Saigon River, Vietnam. Environ. Pollut. 2020, 256, 113412. [Google Scholar] [CrossRef] [PubMed]
  34. Dixit, S.; Yadav, A.; Dwivedi, P.D.; Das, M. Toxic hazards of leather industry and technologies to combat threat: A review. J. Clean. Prod. 2015, 87, 39–49. [Google Scholar] [CrossRef]
  35. Huang, K.M.; Lin, S. Consequences and implication of heavy metal spatial variations in sediments of the Keelung River drainage basin, Taiwan. Chemosphere 2003, 53, 1113–1121. [Google Scholar] [CrossRef]
  36. Li, R.; Tang, X.; Guo, W.; Lin, L.; Zhao, L.; Hu, Y.; Liu, M. Spatiotemporal distribution dynamics of heavy metals in water, sediment, and zoobenthos in mainstream sections of the middle and lower Changjiang River. Sci. Total Environ. 2020, 714, 136779. [Google Scholar] [CrossRef] [PubMed]
  37. Liu, C.; Xu, J.; Liu, C.; Zhang, P.; Dai, M. Heavy metals in the surface sediments in Lanzhou Reach of Yellow River, China. Bull. Environ. Contam. Toxicol. 2009, 82, 26–30. [Google Scholar] [CrossRef] [PubMed]
  38. Zhao, G.; Ye, S.; Yuan, H.; Ding, X.; Wang, J.; Laws, E.A. Surface sediment properties and heavy metal contamination assessment in river sediments of the Pearl River Delta, China. Mar. Pollut. Bull. 2018, 136, 300–308. [Google Scholar] [CrossRef] [PubMed]
  39. Gupta, A.; Rai, D.K.; Pandey, R.S.; Sharma, B. Analysis of some heavy metals in the riverine water, sediments and fish from river Ganges at Allahabad. Environ. Monit. Assess. 2009, 157, 449–458. [Google Scholar] [CrossRef] [PubMed]
  40. Islam, M.S.; Ahmed, M.K.; Raknuzzaman, M.; Habibullah-Al-Mamun, M.; Islam, M.K. Heavy metal pollution in surface water and sediment: A preliminary assessment of an urban river in a developing country. Ecol. Indic. 2015, 48, 282–291. [Google Scholar] [CrossRef]
  41. Nazeer, S.; Hashmi, M.Z.; Malik, R.N. Heavy metals distribution, risk assessment and water quality characterization by water quality index of the River Soan, Pakistan. Ecol. Indic. 2014, 43, 262–270. [Google Scholar] [CrossRef]
  42. El-Amier, Y.A.; Abd El-Gawad, A.M. Assessing the sediment pollution using heavy metals indices in the Nile River Branches in Egypt. J. Environ. Sci. Pollut. Res. 2016, 2, 107–112. [Google Scholar]
  43. Marengo, E.; Gennaro, M.C.; Robotti, E.; Rossanigo, P.; Rinaudo, C.; Roz-Gastaldi, M. Investigation of anthropic effects connected with metal ions concentration, organic matter and grain size in Bormida river sediments. Anal. Chim. Acta 2006, 560, 172–183. [Google Scholar] [CrossRef]
  44. Duodu, G.O.; Goonetilleke, A.; Ayoko, G.A. Comparison of pollution indices for the assessment of heavy metal in Brisbane River sediment. Environ. Pollut. 2016, 219, 1077–1091. [Google Scholar] [CrossRef]
  45. Bábek, O.; Grygar, T.M.; Faměra, M.; Hron, K.; Nováková, T.; Sedláček, J. Geochemical background in polluted river sediments: How to separate the effects of sediment provenance and grain size with statistical rigour? Catena 2015, 135, 240–253. [Google Scholar] [CrossRef]
  46. Sun, X.; Fan, D.; Liu, M.; Tian, Y.; Pang, Y.; Liao, H. Source identification, geochemical normalization and influence factors of heavy metals in Yangtze River Estuary sediment. Environ. Pollut. 2018, 241, 938–949. [Google Scholar] [CrossRef] [PubMed]
  47. Bergaya, F.; Lagaly, G. General introduction: Clays, clay minerals, and clay science. In Handbook of Clay Science, Developments in Clay Science; Bergaya, F., Theng, B.K.G., Lagaly, G., Eds.; Elsevier: Amsterdam, The Netherlands, 2006; Volume 1, pp. 1–18. [Google Scholar]
  48. Tansel, B.; Rafiuddin, S. Heavy metal content in relation to particle size and organic content of surficial sediments in Miami River and transport potential. Int. J. Sediment Res. 2016, 31, 324–329. [Google Scholar] [CrossRef]
  49. Lin, S.; Hsieh, I.J.; Huang, K.M.; Wang, C.H. Influence of the Yangtze River and grain size on the spatial variations of heavy metals and organic carbon in the East China Sea continental shelf sediments. Chem. Geol. 2002, 182, 377–394. [Google Scholar] [CrossRef]
  50. Wang, S.; Cao, Z.; Lan, D.; Zheng, Z.; Li, G. Concentration distribution and assessment of several heavy metals in sediments of west-four Pearl River Estuary. Environ. Geol. 2008, 55, 963–975. [Google Scholar] [CrossRef]
  51. Lin, J.G.; Chen, S.Y. The relationship between adsorption of heavy metal and organic matter in river sediments. Environ. Int. 1998, 24, 345–352. [Google Scholar] [CrossRef]
Figure 1. Sampling sites along Fengshan River in Kaohsiung City, Taiwan.
Figure 1. Sampling sites along Fengshan River in Kaohsiung City, Taiwan.
Applsci 12 01013 g001
Figure 2. Spatial and seasonal variations in heavy metals (As, Hg, Cr, Cu, Ni, Pb, Zn) and organic matter contents in the sediments of Fengshan River. The horizontal dot lines and dash lines represent the ERL and ERM values, respectively.
Figure 2. Spatial and seasonal variations in heavy metals (As, Hg, Cr, Cu, Ni, Pb, Zn) and organic matter contents in the sediments of Fengshan River. The horizontal dot lines and dash lines represent the ERL and ERM values, respectively.
Applsci 12 01013 g002
Figure 3. (a) Contamination factor (CF), (b) geo-accumulation index (Igeo), (c) modified degree of contamination (mCd), and (d) pollution load index (PLI) for heavy metals (As, Hg, Cr, Cu, Ni, Pb, Zn) in the sediments of Fengshan River. The CF is classified as low (<1), moderate (1–3), considerable (3–6), and very high contamination (>6); the Igeo is classified as unpolluted (<0), unpolluted to moderately (0–1), moderately (1–2), moderately to strongly (2–3), strongly (3–4), and strongly to extremely polluted (4–5); the mCd is classified as nil to very low (<1.5), low (1.5–2), moderate (2–4), and high degree of contamination (4–8); the PLI is classified as unpolluted (<1), slightly (1–2), moderately (2–3), and strongly polluted (>3).
Figure 3. (a) Contamination factor (CF), (b) geo-accumulation index (Igeo), (c) modified degree of contamination (mCd), and (d) pollution load index (PLI) for heavy metals (As, Hg, Cr, Cu, Ni, Pb, Zn) in the sediments of Fengshan River. The CF is classified as low (<1), moderate (1–3), considerable (3–6), and very high contamination (>6); the Igeo is classified as unpolluted (<0), unpolluted to moderately (0–1), moderately (1–2), moderately to strongly (2–3), strongly (3–4), and strongly to extremely polluted (4–5); the mCd is classified as nil to very low (<1.5), low (1.5–2), moderate (2–4), and high degree of contamination (4–8); the PLI is classified as unpolluted (<1), slightly (1–2), moderately (2–3), and strongly polluted (>3).
Applsci 12 01013 g003
Figure 4. (a) Potential ecological risk index (RI), (b) mean ERM quotient (mERMQ), (c) the sum of the toxic units (ΣTU), and (d) contamination severity index (CSI) for heavy metals (As, Hg, Cr, Cu, Ni, Pb, Zn) in the sediments of Fengshan River. The RI is classified into low (<150), moderate (150–300), considerable (300–600), and very high ecological risk (>600); the mERMQ is classified into 9% (<0.1), 21% (0.1–0.5), and 49% probability of being toxic (0.5–1.5); the CSI is classified into uncontaminated (<0.5), very low (0.5–1), low (1–1,5), low to moderate (1.5–2), moderate (2–2.5), moderate to high (2.5–3), and high severity of contamination (3–4).
Figure 4. (a) Potential ecological risk index (RI), (b) mean ERM quotient (mERMQ), (c) the sum of the toxic units (ΣTU), and (d) contamination severity index (CSI) for heavy metals (As, Hg, Cr, Cu, Ni, Pb, Zn) in the sediments of Fengshan River. The RI is classified into low (<150), moderate (150–300), considerable (300–600), and very high ecological risk (>600); the mERMQ is classified into 9% (<0.1), 21% (0.1–0.5), and 49% probability of being toxic (0.5–1.5); the CSI is classified into uncontaminated (<0.5), very low (0.5–1), low (1–1,5), low to moderate (1.5–2), moderate (2–2.5), moderate to high (2.5–3), and high severity of contamination (3–4).
Applsci 12 01013 g004
Table 1. The concentrations (mg/kg) of heavy metals (As, Hg, Cr, Cu, Ni, Pb, Zn) in the sediments of Fengshan River in dry and wet season.
Table 1. The concentrations (mg/kg) of heavy metals (As, Hg, Cr, Cu, Ni, Pb, Zn) in the sediments of Fengshan River in dry and wet season.
StationAsHgCrCuNiPbZn
Dry season (December 2015)
F124.30.2115.817.317.319.6148
F237.80.3958.732.726.722.7200
F328.50.3628.618.624.018.6251
F437.70.7231.977.831.931.3353
F525.21.4491.364.634.026.0218
F613.10.9929514357.149.9501
F712.60.3923511839.133.5343
Wet season (June 2016)
F119.50.1826.031.227.211.6119
F241.10.3381.979.437.728.0144
F318.40.2946.648.132.529.3139
F444.10.6646.491.662.144.0292
F521.10.3153.153.251.824.9211
F614.31.3035317564.372.8533
F711.20.4029314841.149.5431
ERL330.1580703035120
ERM851.314539050110270
TEL5.90.17437.335.71835123
PEL170.486901973691.3315
SEL33211011075250820
UCC4.80.059228471767
Criteria values for effect range low (ERL), effect range median (ERM), threshold effect level (TEL), probable effect level (PEL), and severe effect level (SEL) of sediment quality guideline for heavy metals in freshwater ecosystem [26]. Average contents of the upper continental crust (UCC) [21].
Table 2. Comparison among the heavy metal concentrations (mg/kg) of sediments in Fengshan River and those from other regions.
Table 2. Comparison among the heavy metal concentrations (mg/kg) of sediments in Fengshan River and those from other regions.
StationAsHgCrCuNiPbZn
Fengshan River, Taiwan11–440.18–1.4416–35317–17517–6412–73119–533
Keelung River, Taiwan [35]nanana3.5–120na9.3–20041–390
Changjiang River, China [36]9.0–340.01–0.09na14–3726–3417–4272–131
Yellow River, China [37]14–48na41–12830–102na26–7890–202
Xijiang River, China [38]4.9–600.01–0.486–1196–120na11–15528–479
River Ganges, India [39]nana1.8–6.41.0–4.4na4.3–8.410–20
Korotoa River, Bangladesh [40]2.6–52na55–18335–11837–16336–83na
River Soan, Pakistan [41]nana2.9–194.8–5919–357.5–787.3–189
Nile River Branches, Egypt [42]nanana12–681.0–3.715–3030–80
Bormida River, Italy [43]9.0–590.09–0.9663–39216–7728–23113–7748–228
Brisbane River, Australia [44]8.9–131–282–33220–11020–3425–126142–257
na: not available.
Table 3. Spearman correlation matrix for the sediment grain size, organic matter (OM), total nitrogen (TN), total phosphorus (TP), and heavy metals contents (As, Hg, Cr, Cu, Ni, Pb, Zn) in the sediments of Fengshan River.
Table 3. Spearman correlation matrix for the sediment grain size, organic matter (OM), total nitrogen (TN), total phosphorus (TP), and heavy metals contents (As, Hg, Cr, Cu, Ni, Pb, Zn) in the sediments of Fengshan River.
Clay aSilt aSand aOM aTN aTP aAsHgCrCuNiPb
Silt0.911 **
Sand−0.911 **−1.00 **
OM0.4220.534 *−0.534 *
TN0.4770.697 **−0.697 **0.763 **
TP0.1630.433−0.4330.688 **0.842 **
As−0.255−0.3630.363−0.051−0.407−0.538 *
Hg0.2180.332−0.3320.4970.673 **0.376−0.037
Cr0.3430.534 *−0.534 *0.622 *0.815 **0.749 **−0.4950.651 *
Cu0.1650.385−0.3850.648 *0.829 **0.807 **−0.4020.682 **0.842 **
Ni−0.0730.152−0.1520.635 *0.657 *0.741 **−0.3010.535 *0.714 **0.868 **
Pb0.1340.297−0.2970.556 *0.771 **0.714 **−0.3930.680 **0.754 **0.938 **0.820 **
Zn0.2440.371−0.3710.3140.631 *0.486−0.3580.801 **0.622 *0.763 **0.600 *0.780 **
a Grain size, OM, TN, and TP data obtained from Lin et al. [18]. ** significant at p < 0.01, expressed in bold; * significant at 0.01 < p < 0.05; n = 14.
Table 4. Corresponding variable loadings, eigenvalues, and explained variance obtained from PCA/FA for the entire data set.
Table 4. Corresponding variable loadings, eigenvalues, and explained variance obtained from PCA/FA for the entire data set.
VF1VF2VF3
Clay a−0.0960.938−0.043
Silt a0.2750.9170.232
Sand a−0.217−0.951−0.191
OM a0.7620.317−0.331
TN a0.8650.3650.292
TP a0.8230.1650.315
As−0.175−0.192−0.901
Hg0.7110.111−0.187
Cr0.7760.3420.507
Cu0.9100.1220.319
Ni0.881−0.2010.031
Pb0.9150.0420.229
Zn0.8190.1410.352
Eigenvalue6.3873.1301.768
% of Total variance49.1%24.1%13.6%
a Grain size, OM, TN, and TP data obtained from Lin et al. [18]. Bold values indicate strong loadings (>0.7).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lin, K.-N.; Lim, Y.-C.; Chen, C.-W.; Chen, C.-F.; Kao, C.-M.; Dong, C.-D. Spatiotemporal Variation and Ecological Risk Assessment of Heavy Metals in Industrialized Urban River Sediments: Fengshan River in Southern Taiwan as a Case Study. Appl. Sci. 2022, 12, 1013. https://doi.org/10.3390/app12031013

AMA Style

Lin K-N, Lim Y-C, Chen C-W, Chen C-F, Kao C-M, Dong C-D. Spatiotemporal Variation and Ecological Risk Assessment of Heavy Metals in Industrialized Urban River Sediments: Fengshan River in Southern Taiwan as a Case Study. Applied Sciences. 2022; 12(3):1013. https://doi.org/10.3390/app12031013

Chicago/Turabian Style

Lin, Kuan-Nan, Yee-Cheng Lim, Chiu-Wen Chen, Chih-Feng Chen, Chih-Ming Kao, and Cheng-Di Dong. 2022. "Spatiotemporal Variation and Ecological Risk Assessment of Heavy Metals in Industrialized Urban River Sediments: Fengshan River in Southern Taiwan as a Case Study" Applied Sciences 12, no. 3: 1013. https://doi.org/10.3390/app12031013

APA Style

Lin, K. -N., Lim, Y. -C., Chen, C. -W., Chen, C. -F., Kao, C. -M., & Dong, C. -D. (2022). Spatiotemporal Variation and Ecological Risk Assessment of Heavy Metals in Industrialized Urban River Sediments: Fengshan River in Southern Taiwan as a Case Study. Applied Sciences, 12(3), 1013. https://doi.org/10.3390/app12031013

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop