Introduction
Tailings are residual materials of various procedures in metal extraction from different ores, or coal washing processes. Usually milling and hydrometallurgical processes result in a huge volume of residual slurry which may contain heavy metals and many other toxic materials at concentrations higher than environmental standards [1]. Moreover, mining wastes may comprise specific chemical additives, although the concentration levels are generally of no concern [2]. Tailings disposal is an important issue in saving the environment, especially in the case of low grade deposits where the volume of tailing materials is considerable. The size of solid particles of tailing depends on the ore nature and its dressing procedures. As a case in point, tailings in heavy media processes are relatively coarse, whereas those resulted from flotation processes are fine. Usually, in mineral processing plants, tailings are concentrated into high solid percent pulp using thickeners and piled afterwards still having considerable amount of water [3]. Apart from aesthetic imperfections which are caused by stockpiling of tailings, leakage of toxic materials such as reagents and heavy metals may pose serious threats against the environment [4].
EIA is a process aimed to identify, predict, evaluate, and balance the biophysical, social, and other impacts prior to making basic decisions [5]. In fact, it is a tool in environmental administration used to assess the effects of project activities on the environment with an avoidance approach [6]. Regarding the importance and necessity of having full recognition of the area and its environmental status, accurate perception of the impacts caused by project activities, and the need for presentation and classification of the impacts to better demonstrate the results to the decision makers, various techniques are proposed by the researchers [7]. These techniques includes: Ad-Hoc, check lists, matrices, GIS mapping, and media methods [8].
Simple problems having few criteria and options for decision making may be solved with no need to specific methods; however, when the number of criteria and options increase, systematic methods are used to solve the problem and make the proper decisions [9]. Using these techniques help structuring the values and imaginations of decision makers.
One efficient tool in solving multi-objective/multi-criteria problems is multi-criteria decision analysis which is a model of decision making that reasonably optimizes the problem solving using multiple criteria (sometimes heterogeneous [10]. Multi-criteria decision analysis may be performed either by multi-objective decision making or multi-criteria decision making [11]. In multi-criteria decision making of problems, options are prioritized according to various criteria. A decision making problem can be organized in form of classic multi-criteria decision making techniques [12]. These techniques include AHP, SAW, TOPSIS, ANP, and KIKOR that started to develop in 1980s.
In the current study, an attempt was made to identify the best site for dumping tailing in the case study area using TOPSIS method. Additionally, the significant of environmental impact has been assessed using MCDM based on SAW techniques.
Materials and Methods
In the present study the SAW method was used to assess environmental impact assessment of mining activities in the study area, and TOPSIS method was applied to select the best location for the dumping site. Thus, the two mentioned method was described in this section. It should be noted that the questionnaire (n=15) was used to to obtain the necessary data and weight the criteria.
SAW method
The method has been first proposed by Hwang and Yoon in 1981 as a weighted linear combination method. In this method after de-scaling of the decision matrix, the weighted de-scaled decision matrix is obtained by applying weight coefficients of criteria; accordingly the score of each option is calculated [13]. In a multi-criteria decision making problem, if n criteria and m options are present, the decision matrix is as follows:
Where
is the operation of option i
in relation to criterionj
.
In order to de-scale the decision matrix, R matrix is defined as
Where the elements are calculated as
Regarding the significance coefficient of various criteria in decision making, criteria weight vector is defined as
and the best option is selected by
TOPSISI Method
In the present study the TOPSIS method was used to select the best location for the dumping site. In this regards, a negative ideal solution maximizes the cost criteria or attributes and minimizes the benefit criteria or attributes, whereas a positive ideal solution maximizes the benefit criteria or attributes and minimizes the cost criteria or attributes. The TOPSIS method is explained in a succession of six steps as follows:
Step 1: Calculate the normalized decision matrix. The normalized value
is calculated as follows:
i =1, 2, ..., m and j = 1, 2, ..., n.
where
is the weight of the
criterion or attribute and
Step 3: Determine the ideal (
) and negative ideal (
) solutions.
Step 4: Calculate the separation measures using the m-dimensional Euclidean distance. The separation measures of each alternative from the positive ideal solution and the negative ideal solution, respectively, are as follows:
Step 5: Calculate the relative closeness to the ideal solution. The relative closeness of the alternative
with respect to
is defined as follows:
Step 6: Rank the preference order
Case study area
The processing plant of refractory gold ore in As-Sbsulfurcontaining deposits in Zarshoran Gold Mine, Takab, WestAzarbaijan was selected as the case study area.
Result and Discussion
To choose the best siting dump, 9 locations and 10 criteria were defined. The questionnaire was applied to provide quantitative data to compare the sites locations and criteria. The dumping site distance from the mine is the first criteria that should be weighted in the questionnaire. Moreover, the distance from the residential area and sensitive ecosystems are the second and third criteria. The vulnerability to flood and earthquake estimate the vulnerability of the selected sites (A1 to A9) to flood or earthquake based on the expert opinions and distance from the faults. The energy consumption was estimated through the access to water and electricity based on the expert opinions. The ease of access to the local employees is categorized ad supply of human resource that weighted based on the expert opinions. The results of the applied MCDM methods are more described as follows.
Prioritizing locational options in order to dumping tailings using TOPSIS technique
This stage consisted of the following steps:
(1). First, 9 options (A1-A9) were selected as suggested options for the unit establishment.
(2). Important technical and locational sections were determined surveying the state of the art, and a specific score was assigned to each [14].
(3). A specific score (within 0-10) was assigned to each selected point in each notified item (Table 1).
Table 1: scoring to the locational options in each criterion |
decision matrix |
distancefrom the mine |
closenessto residential areas |
impacton surrounding ecosystem |
vulnerabilityto flood and earthquake |
access towater and electricity |
access tovehicles |
unitsecurity |
supplyof human resources |
closenessto agricultural fields |
propertopographical situations |
accessroads |
locational option
impact type |
negative |
negative |
negative |
negative |
Positive |
positive |
positive |
positive |
negative |
positive |
positive |
A9 |
5 |
0 |
2 |
0 |
10 |
10 |
10 |
10 |
2 |
5 |
10 |
A8 |
4 |
2 |
3 |
2 |
5 |
6 |
7 |
10 |
1 |
1 |
5 |
A7 |
2 |
8 |
8 |
7 |
8 |
10 |
10 |
10 |
8 |
5 |
10 |
A6 |
5 |
0 |
2 |
2 |
4 |
5 |
5 |
5 |
1 |
2 |
10 |
A5 |
2 |
0 |
3 |
2 |
5 |
3 |
4 |
3 |
8 |
5 |
3 |
A4 |
0 |
0 |
2 |
5 |
3 |
2 |
3 |
3 |
0 |
0 |
2 |
A3 |
4 |
5 |
5 |
5 |
10 |
10 |
8 |
10 |
5 |
3 |
10 |
A2 |
4 |
0 |
10 |
10 |
10 |
10 |
5 |
10 |
10 |
2 |
10 |
A1 |
5 |
0 |
10 |
10 |
10 |
10 |
5 |
10 |
10 |
2 |
10 |
(4). Next, the input data should be de-scaled using the vector method (Table 2).
Table 2: The normalized data of locational options. |
locational option
i |
distance from the mine |
closeness to residential areas |
impact on surrounding ecosystem |
vulnerability to flood and earthquake |
access to water and electricity |
access to vehicles |
unit security |
supply of human resources |
closeness to agricultural fields |
proper topographical situations |
access roads |
A9 |
0.44 |
0 |
0.11 |
0 |
0.43 |
0.42 |
0.49 |
0.39 |
0.11 |
0.51 |
0.4 |
A8 |
0.35 |
0.21 |
0.17 |
0.11 |
0.21 |
0.25 |
0.34 |
0.39 |
0.05 |
0.1 |
0.2 |
A7 |
0.17 |
0.83 |
0.45 |
0.4 |
0.34 |
0.42 |
0.49 |
0.39 |
0.42 |
0.51 |
0.4 |
A6 |
0.44 |
0 |
0.11 |
0.11 |
0.17 |
0.21 |
0.25 |
0.2 |
0.05 |
0.2 |
0.4 |
A5 |
0.17 |
0 |
0.17 |
0.11 |
0.21 |
0.13 |
0.2 |
0.12 |
0.42 |
0.51 |
0.12 |
A4 |
0 |
0 |
0.11 |
0.28 |
0.13 |
0.08 |
0.15 |
0.12 |
0 |
0 |
0.08 |
A3 |
0.35 |
0.52 |
0.28 |
0.28 |
0.43 |
0.42 |
0.39 |
0.39 |
0.26 |
0.3 |
0.4 |
A2 |
0.35 |
0 |
0.56 |
0.57 |
0.43 |
0.42 |
0.25 |
0.39 |
0.53 |
0.2 |
0.4 |
A1 |
0.44 |
0 |
0.56 |
0.57 |
0.43 |
0.42 |
0.25 |
0.39 |
0.53 |
0.2 |
0.4 |
(5). The criterion weights were normalized by vector method as they are listed in table 3.
Table 3: Normalized values of input weights for locational options |
|
distance from the mine |
closeness to residential areas |
impact on surrounding ecosystem |
vulnerability to flood and earthquake |
access to water and electricity |
access to vehicles |
unit security |
supply of human resources |
closeness to agricultural fields |
proper topographical situations |
access roads |
Criteria weight |
0.081 |
0.1162 |
0.116 |
0.1162 |
0.058 |
0.0581 |
0.0465 |
0.0581 |
0.1162 |
0.1162 |
0.1162 |
(6). The weighted normalized decision matrix may be constructed at this stage (Table 4).
Table 4: Weighted normalized data of the locational options. |
locational option
impact type |
negative |
negative |
negative |
negative |
Positive |
positive |
positive |
positive |
negative |
positive |
positive |
A9 |
0.04 |
0 |
0.01 |
0 |
0.02 |
0.02 |
0.02 |
0.02 |
0.01 |
0.06 |
0.05 |
A8 |
0.03 |
0.02 |
0.02 |
0.01 |
0.01 |
0.01 |
0.02 |
0.02 |
0.01 |
0.01 |
0.02 |
A7 |
0.01 |
0.1 |
0.05 |
0.05 |
0.02 |
0.02 |
0.02 |
0.02 |
0.05 |
0.06 |
0.05 |
A6 |
0.04 |
0 |
0.01 |
0.01 |
0.01 |
0.01 |
0.01 |
0.01 |
0.01 |
0.02 |
0.05 |
A5 |
0.01 |
0 |
0.02 |
0.01 |
0.01 |
0.01 |
0.01 |
0.01 |
0.05 |
0.06 |
0.01 |
A4 |
0 |
0 |
0.01 |
0.03 |
0.01 |
0 |
0.01 |
0.01 |
0 |
0 |
0.01 |
A3 |
0.03 |
0.06 |
0.03 |
0.03 |
0.02 |
0.02 |
0.02 |
0.02 |
0.03 |
0.04 |
0.05 |
A2 |
0.03 |
0 |
0.07 |
0.07 |
0.02 |
0.02 |
0.01 |
0.02 |
0.06 |
0.02 |
0.05 |
A1 |
0.04 |
0 |
0.07 |
0.07 |
0.02 |
0.02 |
0.01 |
0.02 |
0.06 |
0.02 |
0.05 |
(7). The maximum and minimum values in each column were determined (Table 5).
Table 5: The maximum and minimum values of locational options |
|
distance from the mine |
closeness to residential areas |
impact on surrounding ecosystem |
vulnerability to flood and earthquake |
access to water and electricity |
access to vehicles |
unit security |
supply of human resources |
closeness to agricultural fields |
proper topographical situations |
access roads |
Maximum |
0.03552 |
0.096 |
0.06506 |
0.066 |
0.0249332 |
0.0243 |
0.022881 |
0.022912 |
0.0613 |
0.059 |
0.046004 |
Minimum |
0 |
0 |
0.01301 |
0 |
0.0074799 |
0.0049 |
0.006864 |
0.006874 |
0 |
0 |
0.009201 |
(8). Next, the ideal and non-ideal values in each column were determined (Table 6).
Table 6: the ideal and non-ideal values |
|
C01 |
C02 |
C03 |
C04 |
C05 |
C06 |
C07 |
C08 |
C09 |
C10 |
C11 |
ideal |
0 |
0 |
0.01301 |
0 |
0.0249332 |
0.0243 |
0.022881 |
0.022912 |
0 |
0.059 |
0.046004 |
Non-ideal |
0.03552 |
0.096 |
0.06506 |
0.066 |
0.0074799 |
0.0049 |
0.006864 |
0.006874 |
0.0613 |
0 |
0.009201 |
(9). A matrix is constructed for the distance from ideal and non-ideal values (Table 7).
Table 7: Matrix for distance from ideal and non-ideal values |
|
A9 |
A8 |
A7 |
A6 |
A5 |
A4 |
A3 |
A2 |
A1 |
Distance from ideal for each option |
0.0376 |
0.0685 |
0.1248 |
0.0580 |
0.0689 |
0.0843 |
0.0861 |
0.1140 |
0.1160 |
Sum of distance from non-ideal |
0.0248 |
0.0140 |
0.0070 |
0.0198 |
0.0183 |
0.0181 |
0.0081 |
0.0122 |
0.0122 |
(10). Then the similarity index was constructed for each option as it is noted in table 8.
Table 8: Similarity index |
Similarity index for each option |
0.8073 |
0.6330 |
0.4021 |
0.7085 |
0.6624 |
0.6147 |
0.5112 |
0.4922 |
0.4874 |
|
A9 |
A8 |
A7 |
A6 |
A5 |
A4 |
A3 |
A2 |
A1 |
(11). Finally, prioritizing of each option was performed according to the value of options (Table 9).
Table 9: prioritizing of options |
A9 |
A6 |
A5 |
A8 |
A4 |
A3 |
A2 |
A1 |
A7 |
0.8073 |
0.7085 |
0.6624 |
0.6330 |
0.6147 |
0.5112 |
0.4922 |
0.4874 |
0.4021 |
As can be seen in Table 9, A9>A6>A5>A8>A4>A3>A2>A1>A7 with A9 as the most ideal and A7as the least proper options.
Determination the significance of the environmental impacts in zarshoran gold mine
Determination the significance of the environmental impacts has been always a challenging issue in EIA process. In this regard, 40 negative impacts of gold ore processing were extracted (Table 10).
Table 10: Project activity- environmental parameters. |
|
Project activities–environmental factors |
A1 |
Asphalting and widening the access road to the site-soil erosion |
A2 |
Asphalting and widening the access road to the site-air pollution |
A3 |
Asphalting and widening the access road to the site-sound pollution |
A4 |
Asphalting and widening the access road to the site-plants |
A5 |
Asphalting and widening the access road to the site-possessions |
A6 |
Soil excavation and embankment-soil erosion |
A7 |
Soil excavation and embankment-changes in ground morphology |
A8 |
Soil excavation and embankment- quantity of ground and surface water |
A9 |
Soil excavation and embankment-quality of ground and surface water |
A10 |
Soil excavation and embankment-air pollution |
A11 |
Soil excavation and embankment-plants |
A12 |
Soil excavation and embankment-ecosystem habitats |
A13 |
Soil excavation and embankment-social acceptance |
A14 |
Transportation-sound pollution |
A15 |
Construction of gable frames- changes in ground morphology |
A16 |
Construction of gable frames-ecosystem habitats |
A17 |
Establishment of tailings dump, secondary dump, and complementary dump- soil erosion |
A18 |
Establishment of tailings dump, secondary dump, and complementary dump- quantity of ground and surface water |
A19 |
Establishment of tailings dump, secondary dump, and complementary dump- quality of ground and surface water |
A20 |
Establishment of tailings dump, secondary dump, and complementary dump- vegetation cover |
A21 |
subsurface utilities- changes in ground morphology |
A22 |
subsurface utilities- soil contamination |
A23 |
Device installation- sound pollution |
A24 |
Worker's labor-fauna |
A25 |
Landscaping- changes in ground morphology |
A26 |
Extraction of gold- quantity of ground and surface water |
A27 |
Extraction of gold- soil contamination |
A28 |
Transportation of ROM to the pilot site- Safety |
A29 |
Collection and dump of input soil- quantity of ground and surface water |
A30 |
Collection and dump of input soil- soil contamination |
A31 |
Collection and dump of input soil- ecosystem habitats |
A32 |
Collection and dump of input soil- diseases |
A33 |
Activity of crushing unit- quantity of ground and surface water |
A34 |
Activity of crushing unit- sound pollution |
A35 |
Activity of crushing unit- air pollution |
A36 |
Crushed rocks dump-soil contamination |
A37 |
Crushed rocks dump- safety |
A38 |
Activity of grinding unit up to second preparation tank-air pollution |
A39 |
Activity of grinding unit up to second preparation tank-sound pollution |
A40 |
Activity of cyanidation unit up to carbon recovery and providing the product-safety and security |
The selected environmental impacts from gold ore processing is derived from Iranian Leopold matrix. The classification of this matrix has been done with numbers -1 to -5 into 4 classes (Table 11).
Table 11: Classification of Leopold matrix |
-3 |
High |
-1 |
low |
-4, -5 |
Very high |
-2 |
medium |
Table 12 carries the available criteria from national and international references. These include impact nature, magnitude, spatial extent, and duration as main criteria [15-19] and probability of occurrence, ease of implementing mitigation measures as complementary criteria.
Table 12: Assessing the environmental criteria in international references. |
Evaluation
criteria
researchers |
(Antunes et al.2001) |
(Bojórquez-Tapia et al., 1998) |
(Clark et al., 1983 ) |
(Duinker and Beanlands, 1986) |
(Gómez-Orea, 1999) |
(Lawrence, 2003) |
Repetition of criterion in different methods |
magnitude |
* |
* |
|
* |
* |
* |
5 |
spatial extent |
* |
* |
|
* |
* |
* |
5 |
duration |
|
* |
* |
* |
* |
* |
5 |
synergism |
|
* |
|
|
* |
|
2 |
cumulative impacts |
|
* |
|
|
* |
* |
3 |
conflict |
|
* |
|
* |
|
|
2 |
mitigation measure |
|
* |
|
|
* |
* |
3 |
sensitivity of resources |
* |
|
|
|
|
|
1 |
time framework |
* |
|
|
|
|
|
1 |
Vulnerable population |
* |
|
|
|
|
|
1 |
Positive/negative |
|
|
* |
|
* |
|
2 |
Reversibility |
|
|
* |
|
* |
* |
3 |
Direct and indirect |
|
|
* |
|
* |
* |
3 |
probability of occurrence |
|
|
|
* |
|
|
1 |
Ensure the prediction of impact |
|
|
|
* |
|
|
1 |
existence of compatible values |
|
|
|
* |
|
|
1 |
Being periodic |
|
|
|
|
* |
|
1 |
repetition |
|
|
|
|
|
* |
1 |
people and official priorities |
|
|
|
|
|
* |
1 |
level of risk and uncertainty |
|
|
|
|
|
* |
1 |
Results of sensitivity analysis for SAW model is delivered in figure 1.
Figure 1: Sensitivity analysis for SAW model.
After the options were outranked, the values are classified into classes: 1-10 having very high impacts (VH), 11-20 having high impacts (H), 21-30 having medium impacts (M), 31-40 having low impacts (L). Results of this model is given in table 13.
Table 13: Classification of TOPSIS results according to the suggested model. |
Classified values |
Numerical values |
options |
Classified values |
Numerical values |
options |
VH |
0.52 |
A21 |
M |
0.83 |
A01 |
H |
0.69 |
A22 |
L |
0.93 |
A02 |
L |
0.93 |
A23 |
L |
0.94 |
A03 |
H |
0.64 |
A24 |
M |
0.86 |
A04 |
L |
0.97 |
A25 |
M |
0.86 |
A05 |
M |
0.77 |
A26 |
M |
0.79 |
A06 |
VH |
0.2 |
A27 |
M |
0.85 |
A07 |
M |
0.77 |
A28 |
L |
0.89 |
A08 |
H |
0.65 |
A29 |
H |
0.72 |
A09 |
H |
0.57 |
A30 |
L |
0.88 |
A10 |
VH |
0.23 |
A31 |
L |
0.86 |
A11 |
M |
0.75 |
A32 |
VH |
0.29 |
A12 |
VH |
0.21 |
A33 |
VH |
0.39 |
A13 |
M |
0.77 |
A34 |
M |
0.81 |
A14 |
H |
0.67 |
A35 |
L |
0.97 |
A15 |
H |
0.56 |
A36 |
VH |
0.34 |
A16 |
VH |
0.28 |
A37 |
L |
0.94 |
A17 |
H |
0.62 |
A38 |
L |
0.99 |
A18 |
H |
0.63 |
A39 |
VH |
0.25 |
A19 |
VH |
0.05 |
A40 |
H |
0.7 |
A20 |
Conclusions
The achievements of current study could be summarized as follows:
In order to locate the dumping site for tailings, prioritizing of locational options was performed using TOPSIS technique. Prior to this, the decision matrix is constructed after the criteria and their weights are determined.
Determination of environmental impacts significance is performed by specifying the value of each option in each criterion and determination of weight for each criterion before constructing a matrix providing the raw data for decision making techniques.
Subsequently, SAW, TOPSIS, and ELECTRE-TRI methods were applied to classify the options. Sensitivity of each method was analysed and revealed that sensitivity of TOPSIS is maximum (20%) and ELECTRE-TRI has the minimum sensitivity (5%). In other words, ELECTRE-TRI has higher potential ability to determine the environmental impacts significance.