1. Introduction
Manufacturing heavy industries, which are located in the coastal area or near a harbor, are one of the major anthropogenic sources for heavy metal pollution in environment systems such as soil, sediment, and seawater (Yang et al., 2017;Gottesfeld et al., 2018;Kulkarni et al., 2018). Exposed metals released from the industrial process are a great concern about potential biological toxicity as well as metal contamination (Gough et al., 1994;Govil and Narayana, 1999;Acero et al., 2003;Krishna and Govil, 2007; Han et al., 2016).
The harmful heavy metals generated during the process for building or repairing the ship are transferred into sediment and soils through the atmosphere, or released into the coastal area along with the rain. The study on contamination of coastal sediment or toxicity of benthic organisms (shellfishes) induced by manufacturing process in the industrial complex area has been variously carried out (Pekey, 2006;Chen et al., 2007). Untreated wastewater discharged from the neighboring industrial areas (Chen et al., 2007) and metal dust (Pekey, 2006) were significant sources of metal contamination in sediments. However, it has been relatively insufficient to investigate heavy metal pollution in the residential area near the shipbuilding manufacturing industries. Since soil contamination due to anthropogenic processes has been reported to a large extent (Wang et al., 2003;Sharma et al., 2006, 2007;Gottesfeld et al., 2018;Liu et al., 2018), it should be investigate heavy metal contamination that can be attributed to the marine industry not only in the marine environment but also in the vicinity of the heavy industry.
Statistical analysis for metal contamination and biological toxicity is a useful technique to qualify degree of contamination (Li et al., 2017;Yang et al., 2017). There are many indexs such as enrichment factor, geo-accumulation index, potential ecological risk index (PERI), and contamination factor (CF) including spatial distribution, multivariate statistical analysis, principal component analysis (PCA) and correlation analysis. Many researchers have expressed the degree of contamination for heavy metals as numerical value instead of the concentration (Ghrefat and Yusuf, 2006;Li et al., 2017;Yang et al., 2017;Kulkarni et al., 2018), so that level of contamination is easily understood.
In this study, the assessment of heavy metal contamination in field soils was conducted by applying the geo-statistical index. To investigate the degree of contamination level, metal concentration in soils were estimated and measured using inductively coupled plasma-optical emission spectrometer (ICP-OES). And then, geo-statistical index such as enrichment factor (EF), geo-accumulation index, and potential ecological risk index (PERI) was calculated. Metal correlation using statistical method was investigated to identify the potential anthropogenic source.
2. Materials and methods
2.1 Sample collection and analytical procedures
We focused on manufacturing industry which develops large-scale engines, modules for plant facilities, and manufactures stationary floating facility and shipbuilding for various types. Eleven soil samples (5-20 cm depth) were collected using soil auger with a 300 m of manufacturing industries (latitudes of 35°683’N-35°5.715’N and longitudes of 129°2.994’E-129°3.124’E). The position of samples for each sample is summarized in Fig. 1. Collected samples were stored in a polyethylene bag at room temperature and air-dried for 75 hr. Prior to the analysis, the soil samples were sieved through a No. 100 sieve.
To estimate the metal concentration in a soil sample, the samples were analyzed using the proposed method by Jeon et al. (2017). Briefly, the soil was digested with hydrochloric acid (HCl) and nitric acid (HNO3) at 105°C for 2h. And separation for 2h. And separation of supernatant in digested solution was conducted using a centrifuge at 4000 rpm for 20 min, and the supernatant was diluted properly. Metal concentration was analyzed using inductively coupled plasma-optical emission spectrometer (ICP-OES). To ensure the quality of treatment method, the same standard reference proposed by Jeon et al. (2017) was analyzed and metal concentrations varied by less than ± 5 % of the values provided in the SRM. Details of the analytical quality assurance are summarized in Table 1. All samples were measured three times and average value was collected.
2.2 Enrichment factor (EF)
Many researchers have suggested qualifying the extent of heavy metal contamination in the soil. For the estimation of anthropogenic sources, the enrichment factor was often calculated in the soil (Ghrefat and Yusuf, 2006;Chen et al., 2007;Çevik et al., 2009). The metal concentration in soil was normalized to metal concentration of background value, and widely Al (Chen et al., 2007) and Fe (Ghrefat and Yusuf, 2006) for normalization were used. The upper continental crust value (Taylor and McLennan, 1995) was used to background concentration of metal (Upper crust value; Al: 8.04 wt%, Fe: 3.5 wt%, Pb: 20 g/μ g, Zn: 71 μg/g, Cd: 98 μg/g, Cu: 25 μg/g, Cr: 85 μg/g). The enrichment factor (EF) is defined as follows:
Where, (C/Fe or C/Al)sample (C/Fe or C/Al)background represents the relative toxic metal concentration in the sample and geochemical background.
The numerical results are indicated by different pollution level and EF is classified the five levels (Sutherland et al., 2000); EF < 2, minimal; 2 < EF < 5, moderate; 5 < EF < 20, significant; 20 < EF < 40, very high; and EF > 40, extremely high.
2.3 Index of geo-accumulation (IGEO)
IGEO enables and can be calculated by comparing background concentration as below equation.
Where, Cn is the estimated concentration of each element in soil, Bn is the background concentration as geochemical value. The constant 1.5 was considered to the natural variation in the contents of a given substance in the environment and to the effect of small anthropogenic source.
The geo-accumulation index is representative value to compare relative soil contamination and consists of seven grades (Muller, 1969;Yaqin et al., 2008) (Table 2). The second grade may be 10 fold greater than the background value, and sixth grade means extremely contamination with hundredfold greater than background (Teng et al., 2002).
2.4 Ecological toxicity
Biological toxicity can be expressed potential ecological risk index (PERI) proposed by Hakanson in 1980. PERI consists of three modules and is calculated as follow (Jiang et al., 2014); a degree of contamination (CD), toxic-response factor (TR) and potential ecological risk factor (ER).
Where, CDi is a concentration of heavy metals in each sample, CRi is the background value as base line, Cfi is the contamination level of a single element, ERi is the biological hazard of a single element, RI is the overall biological hazard, TRi is biological toxicity weight of a single element with 1for zinc (Zn), 2 for chromium (Cr), 30 for cadmium (Cd), 5 for copper (Cu) and lead (Pb). The meanings of the values calculated by the above equations are shown in Table 3.
3. Results and discussion
3.1 A degree of metal contamination in soils
The individual results obtained of each metal are presented in Table 4 and Fig. 2. The metal concentration had a range of 21.2-85.3 mg/kg, 18.7-44.2 mg/kg, 47.1-335.6 mg/kg, 34.7-338.0 mg/kg, and 0.2-0.8 mg/kg for Pb, Cr, Zn, Cu, and Cd. Overall, no high level of contamination in soil was observed comparing background values in upper crust (Pb: 20 mg/kg, Zn: 71 mg/kg, Cd: 98 mg/kg, Cu: 25 mg/kg, Cr: 85 mg/kg) (Taylor and McLennan, 1995) and government guidelines. But concern levels of pollutants in specific soil showed with giving rise to concern anthropogenic pollution. Especially, metal concentration in specific points such as S-4, S-5, and S-10 was significantly higher than in the other position and it could be assumed that metal particles scattered from heavy industry facility were deposited and accumulated into soil. S-4 and S-5 were classified into area near and far from heavy industry on the basis of roads, and metal concentration in their samples was extremely different according to sample position. The difference of linear distance from the manufacturing industry to the sample is only about 10 meters, but the difference of heavy metal concentration is more than double. It can be seen that the possibility of the dispersion and deposition of heavy metal particles increase as closer from manufacturing industry, and it could be assumed that heavy metal industry is one of the anthropogenic sources for soil contamination. And, characteristically, the concentration of heavy metals in S-3 and S-10 were higher than other samples such as S-2, S-6, S-7, S-8, and S-9, which indicated metal particle scattered from heavy industry could be deposited to the other side of roads. The linear distance from heavy industry is approximately 50-200 m, so it is possibility to disperse metal particles from anthropogenic sources.
Heavy metal contamination can be found in many industries area and disposal dumping, which can release high concentration of heavy metals (Govil and Narayana, 1999). These elements can be leached into ground water or sea water, and later the plant that absorbed heavy metal may affect human health (Gough et al., 1994;Acero et al., 2003;Krishna and Govil, 2007). Although metal contamination in this study was not severe, consideration and monitoring of potential pollution hazards in the future around manufacturing industries are needed.
3.2 Enrichment factor and geo-accumulation index
For the estimation of anthropogenic sources, Enrichment factor (EF) and geo-accumulation index were calculated based on metal concentration (Table 5 and Figs. 3-4). EF factor and geo-accumulation index effectively represent the degree of heavy metal present in the soil with comparing background value (Salomons and Förstner, 1984;Feng et al., 2004). An EF value of more than 5 indicated environmental concern about heavy metal contamination and also IGEO index more than 1. The EF value was not severe level in overall samples and has a range with 1.1-4.4, 0.3-4.9, 1.4-14.0, 0.2-0.5 for Pb, Zn, Cu, and Cr based on Al (Based on Fe: 2.6-10.4 for Pb, 0.3-5.0 for Zn, 1.5-14.4 for Cu, 0.2.-0.6 for Cr) (Fig. 3). Contamination of Cd was not observed unlike other metals (EF < 0.1 based on Al and Fe). The relative values of Pb and Cu were higher than those of other metals. Especially, pollution level for S-4 and S-5 were significant so soil contamination in soil induces by anthropogenic source could be estimated with comparing other samples.
Geo-accumulation index (IGEO) showed similar trend to the enrichment factor. IGEO value had a range with -0.3-1.0, -1.8-1.1, -0.1-2.2, -1.0-1.8 for Pb, Zn, Cu, and Cr (Fig. 4). In specific soil samples, moderately contamination was observed for Pb, Zn, and Cu not Cd. Especially, in S-3-S-5 and S-10, copper contamination was severe (IGEO: 1.6-2.2) and the value was more than two which indicated heavily contaminated condition. The IGEO values of Pb in most soils fell in group 1(uncontaminated to moderately contaminated) except S-4 and S-10 (in group 2) (moderately contaminated). For the Zn metal IGEO grade varied from group 0 (practically uncontaminated) to group 2 (moderately contaminated). And, for the Cu metal, overall IGEO grad was high from group 0 (practically uncontaminated) to group 3 (moderately to heavily contaminated condition). Pollution of Cu and Pb could be due to anthropogenic sources including manufacturing industries.
3.3 Metal correlation in soils
The concentration distribution of heavy metal in soil was similar in places where pollution level was high. Lead, zinc, and copper were commonly observed in contaminated soil, so the correlation was calculated using statistical methods to confirm the correlation of these metals. A relative correlation was expressed as a number from -1 to +1, with 1 indicating definitely positive and -1 indicating a negative relation. Also, to observe these results more clearly and easier, a heating map was used to express correlation degree using red, white and blue color (Table 6). There was a high correlation (p < 0.05) between Pb, Zn, and Cu in soils, and their correlation value was about 0.69 for Zn and 0.72 for Cu based on Pb. And the correlation of Cd between other metals was not significant.
The manufacturing industry in study area develops large-scale engines, modules for plant facilities, and manufactures stationary floating facility and shipbuilding for various types. Various metals are used in the related process, and the metals could be exposed and scattered to the atmosphere and then deposited. Al-Mg-Zn alloy was interested in shipbuilding industries because strength property was better than Al-Mg alloy (Dudzik and Czechowski, 2009). Metals pollution such as Cu, Zn, Cd, and Cr in samples was observed due to shipbuilding actives (Adamo et al., 2005) and also inorganic compounds oxides of Ti, Zn, Cu, Fe, Pb generally used in paint with trace elements (Jolly et al., 2006) for primers, anticorrosive and antifouling (EPA, 1997;EPA, 2000). Manufacturing industry, which had various processes, could be one of the important anthropogenic sources affecting soil contamination.
3.4 Ecological toxicity
The potential ecological risk index (PERI) for Pb, Zn, Cu, and Cd was calculated and presented in Table 7 and Fig. 5. In Fig. 5, the RI which indicated the overall biological hazard was only presented at each sample. The RI value was in medium grade (40 < RI < 80) for S-3, S-5, S-8, and S-10. And, strong biological risk (80 < RI < 160) was only observed in S-4. Among the various metals, ER indicated that the biological hazard of a single element of Cu was largest.
Leung et al. (2017) reported that the level of heavy metal such as Pb, Zn, and Cd in sediment extremely exceeded comparing reference value and their bioaccumulation was also a concern level. The metals released from anthropogenic sources were present in the groundwater, sea water, and sediment and human health could be affected by plant absorbed metals (Acero et al., 2003;Krishna and Govil, 2007). Although the pollution of metal was localized in a specific point, consideration of potential biological risks and steady monitoring are needed.
4. Conclusion
In this study, to estimate the effect of an anthropogenic source to soil contamination, statistical analysis methods such as the enrichment factor, geo-accumulation index, and ecological toxicity index were calculated. Various statistical methods were used because the pollution degree of metals can be estimated differently according to applied value. Metal contamination in soils was not severe but concern level of pollution for metals was observed in some points, especially Zn, Cu, and Pb. And also, strong biological risk (80 < RI < 160) and heavily contamination condition (2 < IGEO) were resulted in specific sample which was located in front of manufacturing industry. The degree of metal contamination by various statistical methods was not severe, but it indicates that it is necessary to consider potential pollution by the heavy manufacturing industry and monitor it continuously.
The metals generated during the process for building or repairing the ship are transferred into sediment and soils through the atmosphere, or released into the coastal area along with the rain. However, it has been relatively insufficient to investigate heavy metal pollution in residential area near the shipbuilding manufacturing industries. In this study, metal contamination caused by the marine industry can be observed not only in the coastal but also in the soil, and it can be used as basic data for prevention against heavy metal pollution. In future work, after additional sampling for field and potential anthropogenic sources, researches as source identification for contamination, transportation mechanism of metals, and characteristic of metals in soil is needed using a variety of analytical methods such as isotope analysis and X-ray technique.