Risk Factor Assessment of Coal Mine Drainage Water on Surrounding Agricultural Soil : A Case Study at Barapukuria in Bangladesh

This study is performed to know the contamination level of heavy metals on the surrounding agricultural soils at the site of Barapukuria Coal Mine Company Limited (BCMCL). Total concentrations of eight heavy metals (Cu, As, Cd, Zn, Pb, Cr, Ni and Mn) were measured in that soil. To calculate the contamination level, the enrichment factor (EF), geoaccumulation index (Igeo) and Hakanson potential ecological risk index were done. Igeo values indicated that the study areas at all stations were suffering from moderate contamination with As. The results of enrichment factor revealed that soils in all stations were in extremely high enrichment condition by As metal. Under Hakanson potential ecological risk index, the contamination coefficient ( f C ) and potential ecological risk index ( i r E ) were calculated. The degree of contamination at BCMCL area reflects moderate contamination. By considering the ecological risk index of a particular heavy metal, the order is as follows: Cd > As > Cu > Pb > Ni > Cr > Zn and Cd has the highest risk index and poses the greatest damage to the BCMCL area’s soil. EF, Igeo values and potential ecological risk index indicate that there is considerable As and Cd pollution, mainly caused by coal mine drainage water.


Introduction
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 [1], thereby adversely impacting on aquatic life and the surrounding vegetation [2].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 [3].
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.

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.

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 Institute of Mining, Mineralogy and Metallurgy (IMMM), Bangladesh Council of Scientific and Industrial Research (BCSIR), Joypurhat.

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.

Determination of Enrichment Factor of Soil Samples
The extent of metal contamination compared to the background area was assessed using the enrichment factor, EF [5] [6].Metal concentrations were normalized to the textural characteristics of soils with respect to Fe [7].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 [8].The crustal abundance data of Krauskopf and Bird [9] were used for EF values.

Determination of Geo-Accumulation Index (I geo ) of Soil Samples
The geo-accumulation index for the metal concentrations in nine stations in BCMCL area was calculated using Muller's [10] expression and regional backgrounds available for this region (Muller's (1979) [11]- [14]).Abrahim and Paker (2008) [15] method was as Equation (2).

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 C is contamination coefficient of a particular heavy metal, i sl C is the measured data of soil heavy metals, and i n C 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 i r T 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 [17]- [19] are shown in Table 1.
Hakanson proposed the highest values of the heavy metals in sediments before the modern industrial era as reference values [20].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 [21] Cd, Hg, As, Pb, Cr, Cu, Zn, and Ni have toxic response factors of 30, 40, 10, 5, 2, 22, 1, and 5, respectively.

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 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 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 [22].
The geoaccumulation Index of some heavy metals of soil samples in the BCMCL area were shown in Table 4. 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     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 [24] [25].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 [25]), 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 ( i f C ) and the overall conta- mination ( d C ) 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 The average value of contamination coefficients was between 0.47 and 5.22, indicating a light level (followed by The potential ecological risk index of a particular heavy metal ( ) i r E in the soil of BCMCL area was calcu- lated 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  response factor proposed by Hakanson has no value for Mn.So, Manganese is not shown in

Conclusions
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 ( geo I ) 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 ( d C level) at BCMCL area reflects moderate contamination.The order of contamination coefficients ( i f C ) 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 ( ) i r E 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.
with support of classification Table 5 of Abrahim and Paker [23].Generally the geo I consists of 7 grades or classes.The geo I classes for 10 studied heavy metals for each stations are listed in Table 6.The geo I 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 geo I 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 geo I

Figure 2 .
Figure 2. Soil contamination indices of heavy metals in the study area.

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)wasexception.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 Table2.

Table 2 .
Metal concentrations in sediments collected from the surrounding area of BCMCL, Dinajpur.

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 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.

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 ( d C ) 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 d C .

Table 8 .
Contamination indices ( i f C ) 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.

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 Table1) 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 Figure2.