Received 9 November 2015; accepted 30 January 2016; published 4 February 2016
Barapukuria coal mine is the first developed coal mine in Bangladesh. Here underground mine water discharged through the drainage system and is usually used for irrigational purpose in the adjoining area. Pollution of water bodies and soil with mine drainage is a major environmental concern worldwide. As coal mine discharge is frequently acidic and contains high concentrations of metals and metalloids, it can create problems of groundwater, surface water and topsoil in the surrounding area of the mine. The low pH water dissolves and mobilizes metals from coal and residue deposits  , thereby adversely impacting on aquatic life and the surrounding vegetation  . Barapukuria coal mine has the potential of polluting groundwater, surface water and soil quality in the surrounding area of mine site with mine discharge water. Little concentration of heavy metal can create harm and long-term effect on the soil and so on to the plant. Varying levels of heavy metals tend to accumulate in coal mining areas. In addition to damage structure, function and quality of soils, this heavy metal contamination directly affects human health and development. The soil heavy metals in the environment are relatively stable and difficult to remove through natural processes  .
The present study is carried out to explore the degree & spatial distribution of metal contamination in agricultural soils in the vicinity of Barapukuria coal mine in northern part of Bangladesh. To measure the contamination level of heavy metals, enrichment factor and geoaccumulation index were used. Hakanson potential ecological risk index was done to assess the risk factor of heavy metals.
2. Materials and Methods
2.1. Study Region
The Barapukuria Coal Mine Company Limited (BCMCL) is located near Barapukuria village of Hamidpur union council under Parbatipur thana, Dinajpur district. The study area and its adjoining area is situated at the northern part of Bangladesh and the area lies between latitudes 25˚31'N to 25˚35'N and longitude 88˚57'E to 88˚59'E (shown in Figure 1) included in the survey of Bangladesh topographic sheet No. 78 c/14.
2.2. Soil Sample Collection and Processing
A total of nine topsoil (0 - 30 cm depth) samples were collected from the discharge canal at definite linear distance of about 200 m during the time of July, 2012. We have taken soil samples from nine different stations. Samples were collected by using a hand auger (a stainless steel screw) and a spade. Samples were taken in a clean polyethylene bags to avoid contaminations. Samples were well labeled with sample ID i.e., St-1 to St-9. Wet samples were air dried before placing in the plastic bags. Soil samples were analyzed in the laboratories of
2.3. Analysis and Determination of Heavy Metals
The selected sediment samples were crushed for 20 min in a planetary ball mill (PM-200, Restsch, Germany) to make powder. The powder samples were then pulverized in a pulverizer machine. The finely ground powder (<75 um) was then put in a porcelain crucible and dried at 100˚C in an oven overnight to remove moisture. The dried powder samples were mixed with binder (stearic acid: sample at a ratio of 1:10) and pulverized for two minutes. The resulting mixture was spooned into an aluminum cap (30 mm). The cap was sandwiched between two tungsten carbide pellets using a manual hydraulic press with 10 tons/sq. in. for 2 min and finally pressure was released slowly. The pellet was then ready for X-Ray analysis. The elements were determined by X-ray fluoresecence (XRF) in the Institute of Mining, Mineralogy and Metallurgy, BCSIR, Joypurhat following the standard procedures using Rigaku ZSX Primus equipped with an end window 4 KW Rh-anode X-ray tube. The trace elements were determined by using crystal LiFi at a 40 kv voltage with 60 mA current. Certified reference material SGR-1 was used in this work. Reference materials were firstly analyzed as unknown samples before analysis of selected sediments to ensure the standard deviation.
2.4. Determination of Enrichment Factor of Soil Samples
The extent of metal contamination compared to the background area was assessed using the enrichment factor,
Figure 1. Location map of study area  (compiled from WA, 1991).
EF   . Metal concentrations were normalized to the textural characteristics of soils with respect to Fe  . According to Aprile and Bouvy (2008), the enrichment factor (EF) for metal concentration in soils at all stations was calculated as Equation (1).
Iron was chosen as the element of normalization because of natural sources (98%) vastly dominates its input  . The crustal abundance data of Krauskopf and Bird  were used for EF values.
2.5. Determination of Geo-Accumulation Index (Igeo) of Soil Samples
The geo-accumulation index for the metal concentrations in nine stations in BCMCL area was calculated using Muller’s  expression and regional backgrounds available for this region (Muller’s (1979)  - ). Abrahim and Paker (2008)  method was as Equation (2).
The sediment pollution was investigated by following the concentration of nine studied trace metals such as Fe, Mn, Zn, Ni, Pb, Cu, Cd, Co and As. The factor 1.5 was introduced to minimize the effect of possible variations in sediments  . According to crustal abundance data of Krauskoph and Bird (1995)  , the reference samples were Fe: 50,000, Mn: 950, Zn: 70, Ni: 75, Cr: 100, Pb: 13, Cu: 55, Cd: 0.2, Co: 25 and As: 1.8 mg/kg.
2.6. Determination of Hakanson Potential Ecological Risk Index
The potential ecological risk index method was proposed by Hakanson from a sedimentary perspective to assess the characteristics and environmental behavior of heavy metal contaminants in soil sediments. In addition to considering the heavy metals in soil, the Hakanson method assesses their potential ecological and environmental effects with toxicology. The assessment, which is conducted using parallel and equivalent index classification, provides a quantitative method of directly isolating the extent of potential hazards. This method is able to reflect the effects of various contaminants and reveal the comprehensive influence of multiple contaminants in a particular environment. The method comprises a single contamination coefficient, a comprehensive contamination measure, the toxic response factor for heavy metals, and a potential ecological risk index. The specific calculating formulas are as follows:
Equation (3) is the single contamination coefficient
where is contamination coefficient of a particular heavy metal, is the measured data of soil heavy metals, and is the reference value.
The comprehensive contamination was measured by using Equation (4)
Equation (5) is the potential ecological risk index of a particular heavy metal
where is the toxic response factor.
The potential ecological risk index was calculated by using Equation (6).
The corresponding degrees of contamination and the grading standards for the levels of potential ecological risk in, , , and RI based on relevant studies  - are shown in Table 1.
Hakanson proposed the highest values of the heavy metals in sediments before the modern industrial era as reference values  . In this study background values were used as reference values to accurately reflect the pollution levels at BCMCL area. According to standardized toxic response factor proposed by Hakanson  Cd, Hg, As, Pb, Cr, Cu, Zn, and Ni have toxic response factors of 30, 40, 10, 5, 2, 22, 1, and 5, respectively.
3. Results and Discussion
The concentration of Fe was between 16,108 to 38,734 ppm, with an average value of 22,828.6 ppm. Fe concentration was found 38,734 ppm at St-6 and minimum was 16,108 ppm at St-2. The concentration of As was between 34.49 to 39.72, with an average value of 36.2, and the value of St-9 sample was above that of the background soil. The concentration of Cr was between 33.62 and 90.7, with an average value of 47.39 and value of
Table 1. Degree of contamination for particular heavy metals and the corresponding grading standards for potential ecological risk.
St-8 was below that of background soil but only St-1 had higher value (90.7) than that of its background. The Cu ranges from 40.1 to 51.49, with an average value of 43.75 and value of St-9 was above that of soil background. The concentration of Mn was between 117.15 and 900.52 with an average 261.79 and value of St-8 was below that of its background soil but St-7(900.52) was exception. The concentration of Zn was between 44.43 and 110.56 with an average value of 57.63 and the value of St-5 was above that of background soil. The concentration of Ni was between 8.51 and 44.58 with an average value of 22.34 and the value of St- 4 was above that of background soil. The concentration of Pb was between 22.12 and 24.81 with an average value of 22.55 and values of all stations were little higher than that of its background soil. The Cd concentration in the studied samples were in lower concentration at all stations (0.76 - 0.89). The soil distribution characteristics of heavy metals in the study area are shown in Table 2.
Table 2. Metal concentrations in sediments collected from the surrounding area of BCMCL, Dinajpur.
The heavy metal concentration in BCMCL area were in this order Fe > Ba > Mn > Zr > Sr > Cr > Cu > Zn > As > Ni > Mo > Pb > Co > Sb > Se > Sn > Cd. Radioactive elements such as Thorium & Urenium in this samples ranges from 4.12 to 5.14 & 4.15 to 5.28 respectively. Sulphur content of this studied area were very low (0.017 - 0.28) and Carbon content little high (6.13 - 20.93). The results were validated using two Certified Reference Materials (Soil-7 and SGR-1) as shown in Table 3. The experimental data was compared with certified data. Higher concentration of Ba (353.8) & Zr (133.8) was found in the analyzed soil as compare to SGR-1 values. In study area the concentration of Cu, Ba, Sb and Se were higher as compared to Soil-7, Certified Reference Materials UGCS/SGR-1 and IAEA/Soil-7  .
The geoaccumulation Index of some heavy metals of soil samples in the BCMCL area were shown in Table 4. values were interpreted with support of classification Table 5 of Abrahim and Paker  . Generally the consists of 7 grades or classes.
The classes for 10 studied heavy metals for each stations are listed in Table 6. The values revealed that the value of Ni, Fe, Mn, Cr, Co, Zn and Cu in all stations fell into class 0.This indicated that the sediments in these stations were uncontaminated by these metals. The values for Pb & Cd in all stations were fell in class 1. This indicated that sediments in these stations were uncontaminated to moderately contaminated by Pb and Cd. The values for As at all the stations fell in class 2 and this indicated that sediments in these stations were moderately contaminated by As. Some studied heavy metals (Ni, Mn, Cr, Co, Zn & Cu) occured no contamination in that area. But the studied area at all stations were suffering from moderately contamination with As. That means As-toxicity was in moderate condition in BCMCL area.
The results of enrichment factor for nine studied heavy metals for each station were listed in Table 7. There is
Table 3. Comparison of detected values of heavy metals with two certified reference materials.
Table 4. The geoaccumulation Index of some heavy metals of soil samples in the BCMCL area, Dinajpur.
Table 5. The degree of metal pollution of seven enrichment classes (Abrahim & Paker, 2008).
Table 6. The results of geoaccumulation index classes of heavy metals of BCMCL area, Dinajpur.
Table 7. The results of enrichment factor of heavy metals of sediment samples of BCMCL area.
no accepted pollution ranking system or categorization of degree of pollution on the enrichment ratio and/or factor methodology. Five contamination categories are recognized on the basis of enrichment factor: EF < 2 states deficiency to minimal enrichment, EF = 2 − 5 moderate enrichment, EF = 5 − 20 significant enrichment, EF = 20 − 40 very high enrichment and EF > 40 extremely high enrichment   . Enrichment factors of heavy metals were calculated for each soil sample relative to the background values of abundance of chemical elements in the continental earth crust, choosing Fe as the reference element.
According to the classification (Sezgin et al., 2003  ), Cd results (12.67 - 4.90) showed that eight stations have significant enrichment and St-6 have moderate enrichment. The heavy metals Zn, Cu, Mn and Ni showed deficiency to minimal enrichment. Mn (1.76 - 0.32) showed deficiency to minimal enrichment at every station. As ranges from 24.73 - 60.87 which have extremely high enrichment at eight stations and only St-6 (24.73) have very high enrichment. These results revealed that the soils in these stations were in extremely high enrichment condition by As metal. Co ranges from 1.45 - 0.84 which has deficiency to minimal enrichment. We saw Cr (1.99 - 0.51) have deficiency to minimal enrichment at all stations. Cu ranges from 2.35 - 0.98 that means six stations have (<2 values) deficiency to minimal enrichment and at stations 2, 4 & 9 (2.35, 2.01 & 2.23) have moderate enrichment. In case of Pb (5.24 - 2.12) seven stations have moderate enrichment and St-2 & 9 (5.24 & 5.18) have significant enrichment. Zn ranges from 2.93 - 0.98 that means seven stations have deficiency to minimal enrichment but station 2 and 7 showed exception (2.12 & 2.93) moderate enrichment. For Ni (range 1.14 - 0.20), every station have deficiency to minimal enrichment.
Hakanson potential ecological risk method: The single contamination indices () and the overall contamination () of heavy metals in BCMCL soils were shown in Table 8. The degree of soil contamination at the Barapukuria coal mine area ranged from 12.48 to 18.20, with an average value of 14.27, reflecting moderate contamination. The order of contamination of the eight heavy metals is Cu > As > Cd > Zn > Pb > Cr > Ni > Mn and the order of contamination at the sites is St-7 > St-5 > St-1 > St-6 > St-8 > St-4 > St-2 > St-3 > St-9. The average value of contamination coefficients was between 0.47 and 5.22, indicating a light level (followed by Table 1) of contamination. Cu had the largest contamination coefficient, especially in the St-7 (6.15), indicating heavy contamination. The contamination coefficient of As ranges from 3.64 to 4.19, that indicates heavy contamination at every stations. In case of Pb, the contamination coefficient at all stations was >1 that indicates light contamination. Comprehensive contamination () ranges 12.48 to 18.50, average value is 14.27 which represents moderate level of contamination and only St-7 has relatively high class of.
The potential ecological risk index of a particular heavy metal in the soil of BCMCL area was calculated according to Hakanson index method and its division standard and the results are shown in Table 9. The risk coefficients of As in most of soil samples were low. Among seven heavy metals, the level of contamination posed low potential (followed by Table 1) ecological risk, with an index of between 1.09 to 38.25 and only Cd has 61.21. The risk coefficients of As, Pb, Cr, Ni, Zn, Cd and Cu at every station was low, indicating that these seven heavy metals had a limited environmental impact on BCMCL area. The order of ecological risk associated with seven heavy metals is Cd > As > Cu > Pb > Ni > Cr > Zn. Cd has the highest potential ecological risk index and poses greatest damage to the BCMCL area’s soil. Zn plays an important role in plant nutrition & has a relatively lower toxicity response factor. Here Zn represents least ecological hazard.
The potential ecological risk index (RI) is a comprehensive measure of the degree of contamination and potential ecological hazard of total heavy metals. Here seven heavy metals were calculated because standard toxic
Table 8. Contamination indices () and the overall contamination of heavy metals in BCMCL soils.
Table 9. Ecological risk factor and potential ecological risk index (RI) of heavy metals in BCMCL area.
Figure 2. Soil contamination indices of heavy metals in the study area.
response factor proposed by Hakanson has no value for Mn. So, Manganese is not shown in Table 9. The potential ecological risk index of heavy metals found at BCMCL area ranged from 126.24 to 153.94, with an average value of 137.13, representing a low ecological risk (followed by Table 1) and only st-7 fell under moderate class of RI. Thus comprehensive treatment of heavy metals especially for Cd should be started to improve environmental protection at that site. The soil contamination indices or pollution index (PI) of heavy metals in study area are shown graphically in Figure 2.
The soil samples collected from surrounding area of Barapukuria Coal Mine Company Limited have been studied in this study. The concentration of heavy metal (Cd, As, Pb, Cr, Cu, Zn, Ni and Mn) and their contamination levels of that area were determined. The enrichment factor, geoaccumulation index and Hakanson potential ecological risk index were used to analyze the level of contamination and potential ecological toxicity of heavy metals in that area.
1) The enrichment factor (EF) results indicated that the soil of BCMCL area was in extremely high enrichment condition by As metal. Co has deficiency to minimal enrichment at that site. Pb has moderate enrichment in seven stations and significant enrichment in two stations.
2) According to geoaccumulation index () values, the sediment of BCMCL area was suffering from uncontaminated to moderate contamination with Pb and Cd. The studied areas at all stations were suffering from moderate contamination with As. That means As-toxicity was in moderate condition in BCMCL area.
3) According to Hakanson potential ecological risk index, the level of soil contamination (level) at BCMCL area reflects moderate contamination. The order of contamination coefficients () of eight heavy metals is as follows: Cu > As > Cd > Zn > Pb > Cr > Ni > Mn. The average value of contamination coefficients indicates a light level of contamination. The contamination coefficient of As and Cu indicates heavy contamination at that area. So, targeted pollution control and management measures are urgently required to prevent the increase of As & Cu contamination and to limit potential harm.
4) By considering potential ecological risk index of a particular heavy metal the order is as follows: Cd > As > Cu > Pb > Ni > Cr > Zn. Cd has the highest potential ecological risk index and poses greatest damage to this area. The potential ecological risk index (RI) and RI grade showed low ecological risk in this area.
5) The overall results showed that there is a moderate contamination at BCMCL area. Therefore the contamination control and comprehensive treatment of heavy metals should be strengthened, specifically for As, Cu and Cd to improve the ecological restoration and environmental protection at that site.
The Research and Development (R&D) project has been submitted to the BCSIR through Director, IMMM for financial year 2010-2011. The project was approved at 25th May, 2011. We would like to thank Chairman, BCSIR for approval of this R&D project and guideline to complete the project. We also thank to Director, IMMM, BCSIR, Joypurhat for his cooperation and inspiration. We are also grateful to the scientists, technicians and lab attendants for their support to complete the project.
 Cherry, D.S., Currie, R.J., Souek, D.J., Latimer, H.A. and Trent, G.C. (2001) An Integrative Assessment of a Watershed Impacted by Abandoned Mined Land Discharges. Environmental Pollution, 111, 377-388.
 Gergen, I. and Harmanescu, M. (2012) Application of Principal Component Analysis in the Pollution Assessment with Heavy Metals of Vegetable Food Chain in the Old Mining Areas. Chemistry Central Journal, 6, 156.
 Selvaraj, K., Ram Mohan, V. and Szefer, P. (2004) Evaluation of Metal Contamination in Coastal Sediments of the Bay of Bengal, India: Geochemical and Statistical Approaches. Marine Pollution Bulletin, 49, 174-185.
 Woitke, P., Wellmitz, J., Helm, D., Kube, P., Lepom, P. and Litheraty, P. (2003) Analysis and Assessment of Heavy Metal Pollution in Suspended Solids and Sediments of the River Danube. Chemosphere, 51, 633-642.
 Aprile, F.M. and Bouvy, M. (2008) Distribution and Enrichment of Heavy Metals in Sediments at the Tapacura River Basin, Northeastern Brazil. Brazilian Journal of Aquatic Sciences and Technology, 12, 1-8.
 Rubio, B., Nombela, M.A. and Vilas, F. (2000) Geochemistry of Major and Trace Elements in Sediments of the Ria de Vigo (NW Spain): An Assessment of Metal Pollution. Marine Pollution Bulletin, 40, 968-980.
 Sahu, K.C. and Bhosle, U. (1991) Heavy Metal Pollution around the Island City of Bmbay, India. Part 1: Quantification of Heavy Metal Pollution of Aquatic Sediments and Recognition of Environmental Discriminants. Chemical Geology, 91, 263-283.
 Subramanian, V. and Mohanachandran, G. (1990) Heavy Metals Distribution and Enrichment in Sediments of Southern East Coast of India. Marine Pollution Bulletin, 21, 324-330.
 Abrahim, G.M.S. and Paker, R.J. (2008) Assessment of Heavy Metal Enrichment Factors and the Degree of Contamination in Marine Sediments from Tamaki Estuary, Auckland, New Zealand. Environmental Monitoring and Assessment, 136, 227-238.
 Al-Haidarey, M.J.S., Hasan, F.M., Al-Kubaisey, A.R.A. and Douabul, A.A.Z. (2010) The Geoaccumulation Index of Some Heavy Metals in Al-Hawizeh Marsh, Iraq. E-Journal of Chemistry, 7, S157-S162.
 Dong, J., Bian, Z. and Wang, H. (2007) Comparison of Heavy Metal Contents between Different Reclaimed Soils and the Control Soil. Journal of China University of Mining and Technology, 3, 531-536.
 Lee, S., Lee, B., Kim, J., Kim, K. and Lee, J. (2006) Human Risk Assessment for Heavy Metals and as Contamination in the Abandoned Metal Mine Areas, Korea. Environmental Monitoring and Assessment, 119, 233-244.
 Qu, C., Sun, K., Wang, S., Huang, L. and Bi, J. (2012) Monte Carlo Simulation-Based Health Risk Assessment of Heavy Metal Soil Pollution: A Case Study in the Qixia Mining Area, China. Human and Ecological Risk Assessment, 18, 733-750.
 Duzgoren-Aydin, N.S., Wong, C.S.C., Aydin, A., Song, Z., You, M. and Li, X.D. (2006) Heavy Metal Contamination and Distribution in the Urban Environment of Guangzhou, SE China. Environmental Geochemistry and Health, 28, 375-391.
 Sezgin, N., Ozcan, H.K., Demir, G., Nemlioglu, S. and Bayat, C. (2003) Determination of Heavy Metal Concentrations in Street Dusts in Istanbul E-5 Highway. Environment International, 29, 979-985.