Weighting Criteria and Prioritizing of Heat stress indices in surface mining using a Delphi Technique and Fuzzy AHP-TOPSIS Method

Background Heat stress as a physical harmful agent can increase the risk of health and safety problems in different workplaces such as mining. Although there are different indices to assess the heat stress imposed on workers, choosing the best index for a specific workplace is so important. Since various criteria affect an index applicability, extracting the most effective ones and determining their weights help to prioritize the existing indices and select the optimal index. Methods In order to achieve this aim, present study compared some heat stress indices using effective methods. The viewpoints of occupational health experts and the qualitative Delphi methods were used to extract the most important criteria. Then, the weights of 11 selected criteria were determined by Fuzzy Analytic Hierarchy Process. Finally, fuzzy TOPSIS technique was applied for choosing the most suitable heat stress index. Results According to result, simplicity, reliability, being low cost, and comprehensiveness were the most determinative criteria for a heat stress index. Based on these criteria and their weights, the existing indices were prioritized. Eventually, wet bulb glob temperature appropriated the first priority and it was proposed as an applicable index for evaluating the heat stress at outdoor hot environments such as surface mines. Conclusions The use of these strong methods allows introducing the most simple, precise, and applicable tool for evaluation the heat stress in hot environments. It seems that WBGT acts as an appropriate index for assessing the heat stress in mining activities at outdoors.


Background
Most workplaces depending on specific jobs and tasks may include the chemical, physical, biological, and ergonomic hazards that they cause a variety of effects ranging from minor injuries to serious diseases and even death. Heat is one of the most important physical hazards in many workplaces. Heat, as an energy source, exist in many industrial processes and workers who expose to it may be at risk of health problems and disorders [1]. Therfore, this adverse agent is an undeniable fact and considered as a major health problem around the world, especially in developing countries [2].
Heat stress is a combination of internal and external factors leading to heat related illnesses. The internal factors include internal body temperature, metabolic rate, and physiological adaptations to hot environments. The air temperature, thermal radiation, air velocity, and humidity contribute as the external factors [3]. Short-term exposure to extreme heat (acute exposure) can lead to rise the core body temperature, which it may directly cause heat related illnesses such as mild rash, cramps, heat exhaustion and heat stroke. It is reported that long-term chronic exposure to heat leads to chronic kidney diseases [4], cardiovascular diseases, and mental health problems [5].
In some industries, heat stress is a common problem by which workers continuously exposed to. These conditions are found in indoor workplaces such as ceramics, foundry, and glass and outdoor works like mining, agriculture and construction industries [6].
Mining activities lead to workersʼ exposure to undesirable environmental conditions. There are various causes for heat stress in mines. The sources of heat in surface mines include sun, electrical tools, machines, and mechanical processes. The efficiency of diesel engines typically is 33% and two-thirds of input energy is lost as heat, resulting in irreparable damages and disorders [7]. According to a study in New Zealand and India, miners are highly vulnerable to heat related illnesses than other workers, especially workers of indoor workplaces [6]. In 2013, an investigation assessing the sings and symptoms of heat-related illnesses showed 87% of workers in surface mines have experienced such symptoms, so that over 90% of them have occurred in warm seasons and over 80% happened more than once. This frequency of occurrences demonstrates the experiencing of heat-related symptoms is a common problem among workers of surface mines. High body temperature, headache, fatigue, and muscles cramps were the most common symptoms [8].
In Iran there are a lot of jobs and tasks with heat-related problems due to the weather conditions and the hot nature of most mining activities. According to the Statistical Centre of Iran in 2012, the total miners who worked in 6206 surface mines was 84528 [7]. Unfortunately, there is no statistics about miners exposed to heat stress and also no comprehensive study has been done to assess the heat stress and determine an index for mine workers. Therefore, it is necessary to determine an appropriate index to present the heat permissible limit due to geographical, environmental, and personal conditions.
There is a need for using the heat indices to assess the heat stress, protect the workers against extreme heat, and control the workplace heat stress. Based on literatures, since 1905 many efforts have been made to measure the levels of heat stress in workplaces and estimate the heat strain. The aim of presenting an suitable index is to specify the relationship between environmental parameters, clothing and activity as a number. Due to the complexity of variables related to heat stress in real working environment, determining a comprehensive heat stress index for each variable is maybe not possible. Differences in the physical work demands, health status, heat tolerance, and heat sources, as well as involvment of the heat-related mechanisms, job rotation, and environmental temperature may affect the risk levels [9]. Consequently, there have been over 60 heat stress indices indicating the scientific basis to establish the safety standards and limits for workers in hot environments. Heat indices are divided into 3 groups: logic indices based on heat balance equation, the experimental indices on the basis of objective and subjective measurement of heat strain, and direct indices based on the measurement of environmental variables [9].
The history of introducing and applying the heat stress indices is related to over a hundred years ago. The wet-bulb temperature index was firstly proposed by Haldane in 1905 as a measure to express the heat stress. After that, many indices were developed and applied around the world. Some of them include wet bulb glob temperature (WBGT), Effective temperature index, the Oxford index, 4-h Sweat Rate (P4SR), The new Universal Thermal Climate Index (UTCI). But none has been widely accepted [10].
This study was aimed to determine a heat stress index among the existing ones for sarface mines using a set of criteria. Considering the lack of certain and clear criteria for choosing a heat stress index, at first it was necessary to get the viewpoints of occupational health experts and use a qualitative method (Delphi) for defining some criteria to select an index.
Delphi is a method for acquiring the group knowledge. It has a structural process to predict and make decisions by a series of rounds, gather the information and eventually, achieve a group consensus [11]. Delphi is a systematic process to extract the experts opinions on a specified topic through a series of questionnaires. The anonymity of respondents and controlled feedback are essential [12]. Delphi is a simple, easy to use, and low cost tool that can be applied for gaining judgments on complex matters in lack of precise information.
The other advantages of the Delphi Method include high flexibility, elimination of the geographical impediments, providing open discussions, and focusing attention on the relevant issues [13]. Detailed planning and considering the effective components of Delphi study should be taken into account. Then, the weight of each criterion was determined by Fuzzy Analytic Hierarchy Process (AHP). Finally, using the fuzzy TOPSIS technique, a heat stress index was selected among the existing indices.

Methods
This study is a qualitative practical research that was designed and implemented in 11 steps. Due to the fact that there were no common criteria for selecting a heat stress index, the Delphi technique was used.
This study was accomplished in following steps: 1. The formation of implementing team and monitoring the Delphi process 2. Selecting the experts and participants 3. Adjusting the questionnaire for first round 4. Editing the questionnaire grammatically (deductive and remove ambiguities) 5. Sending the questionnaire for experts 6. Analyzing the obtained responses in the first round 7. Preparing the second round questionnaire considering the required revisions 8. Sending the second questionnaire to the same experts 9. Analyzing the results of second questionnaire 10.Determining the relative weights of each criterion using the fuzzy AHP. 11.Choosing a heat stress index among the existing ones in the study using the fuzzy TOPSIS method.
A team including experts and consultants was formed and the aims of study were reviewed and revised. The number of 30 experts and panel members was selected from occupational health professionals of medical universities, administrative experts working in the department of environmental health in ministry of health, post-graduate students who were working on a thesis related to heat stress. According to studies, if the group is homogeneous, a sample consistent of ten to fifteen people will be sufficient for Delphi [14]. Therefore, selecting 30 participants in this study ensures the validity of findings in case of probable loss of respondants.
As a pilot test, an open question was designed and sent to 10 experts in order to survey their perceptions and response. Given the results of the pilot study, the questionnaire was sent to 30 mentioned participants. Then, the responses of first phase of study were analyzed, criteria extracted and prioritized. The criteria with less scores were also investigated. During several meetings with professors, some overlapped criteria with the same concept were merged. At the end of meetings, the criteria with the highest scores were extracted and used for prioritizing and pairwise comparing in the next step. The second questionnaire was designed and sent to the 30 previous members for pair-wise comparison between the criteria and finally determinig the weight of each criterion using Fuzzy AHP. For prioritizing the indices existing in the study using fuzzy TOPSIS, a questionnaire including m alternatives and n criteria was prepared and sent to the same experts. They were asked to assign a score for each alternative using the linguistic scale in order to select the most appropriate heat stress index for surface mines. Indices in this study were used with respect to their prevalence and validity, strong correlations with physiologic indices and other valid ones, and also their applicability for outdoor environments. Figure 1 indicates the conceptual model of the study.

Fuzzy analytic hierarchy process
The Analytic Hierarchy Process (AHP) was proposed first time by Thomas. El. Saaty. This technique combines the opinions of all experts and provides the decision-making in a simple way. Then, a scale is used to pair-wise comparisons of criteria and assess their priorities [15]. This technique assays the quantitative and qualitative indices efficiently [16]. The advantages of this method include formulating the problem in question, improving the consistency of judgments, handling and solving the various problems, obtaining the opinions of members for making decision, aggregating the judgments of experts to determine the best alternative, and priotitizing through the pair-wise comparisons of criteria. The other advantage of AHP is the use of qualitative criteria for decisionmaking and expressing the results quantitatively by mathematical techniques. Using the quantitave results leads to intelligible data and accurate judgments [17]. the use of expert opinions based on the Delphi method ensures the acceptable results in the decision-making process. However, there are some disadvantages in the use of this method. For instance, AHP provides judgment based on the certain numbers without the reflecting of experience and knowledge of experts. To overcome this weakness, the combination of AHP and fuzzy logic and the use of fuzzy numbers are recommended. This approach can help resolve the weakness of criteria weighting [18]. After determining the criteria by Delphi technique, the following steps were performed to specify the weights of criteria using fuzzy AHP: Determining the linguistic terms for pair-wise comparison the criteria [19] using Table 1. Using the triangular fuzzy numbers to form the pairwise comparison matrix.
In this study, to remove ambiguity of uncertainties in decision-making, the triangular fuzzy numbers presented in Table 1 were used in all steps for pair-wise comparisons in AHP. A triangular fuzzy number is composed of three factors, as Ã = (l + m + u), and its membership function is illustrated as follow: Where, m is the maximum grade of membership function, l and u are the lower and upper bounds to allow more reasonable evaluation. A triangular fuzzy number (m, l, u) and the membership function μ (X) is expressed in Fig. 2. Matrix 2 shows the triangular membership function for linguistic terms [20].
2. Calculating the fuzzy weight of each elements.
3. The use of best non fuzzy performance (BNP) 4. Calculating the inconsistency rate of matrix The subjective comparisons in this method leads to some inconsistencies in AHP. The consistency ratio (CR) can be calculated as follow. If the CR is unacceptable, pair-wise comparison should be reconsidered.
The consistency index (CI) presents the deviation from the consistency and is calculated as follow: where, n is the size of pair-wise comparison matrix, λ max is the maximal value of comparison matrix, and RI is the random consistency index obtained randomly from related tables [15]. If CR < 0.1, the comparisons are acceptable. The CR > 0.1 indicates the inconsistent judgments and pair-wise comparison should be revised [22]. The TOPSIS is a multi criteria decision-making method introduced firstly by Hwang and Yoon. The basic concept of TOPSIS is providing an ideal solution to maximize the benefit criteria and minimize the cost criteria. Briefly, the positive ideal solution includes the best available criteria while the negative ideal solution is consistent of the worst values of criteria. Therefore, optimized alternative has the least distance from the positive ideal solution and the farthest distance from the negative one. The TOPSIS algorithm is a strong compensatory multi-criteria decisionmaking technique to prioritize the alternatives with regard to the ideal solution. It has a little sensitivity to type of weighting technique and its responses don't significantly change. In addition, an advantage of this model is its ability for fast identifying the best alternative. In this algorithm, it is assumed that each index and criterion in decision-making increases or decreases uniformly [23].
It is often difficult for decision-makers to present an precise value for criteria due to some errors in this respect. In this case, it is recommended the use of fuzzy numbers for assessments. Therfore, in this study the fuzzy TOPSIS was applied for prioritization. The steps of fuzzy TOPSIS are as follows: Step 1. Constructing the initial comparison matrix. The fuzzy values in this matrix have been defined in the scale of membership functions in Table 2.
Step 2. Determining the positive (benefits) and negative (costs) ideal solutions as A+ and A-respectively.
where J 1 and J 2 respectively present the criteria benefit and cost.
Step 3. Calculating the distance of each alternative from negative and positive ideal solutions by Euclidian (n-dimension) method.
Where di + and d − represents the distance of alternative from positive and negative ideals. Also, distance (d) between two fuzzy numbers (Pij, Pi + ) is calculated by following equation: Step 4. Calculating the relative closeness of each alternative to the best solution as follow: Step 5. Prioritizing the alternatives using the relative closeness. Each alternative which has the largest closeness is the best for choosing.

Results
The aim of this study was to determine the effective criteria for choosing a heat stress index by using the Delphi model and specify those criteria weights by fuzzy AHP. Using 30 completed questionnaires by experts and reviewing their viewpoints, 30 criteria were extracted. In another meeting, due to the same concepts of some criteria and the lack of relevance of some of them, the effective criteria to select the index decreased to 14. Finally, 11 criteria were selected considering the percentage of experts responses. They can be seen in Table 3. According to the table, the Simplicity criterion by 76.6% appropriated the highest percentage of responses.
As shown in Table 2, in this study 11 criteria were considered with regard to problem in question. Then, the selected criteria with high agreement were sent back to experts for pair-wise comparisons. Finally, the weight of each criterion was analyzed by fuzzy AHP considering the importance of criteria toward the interested aims. The consistency rate (CR) of this study was 0.083, lower than 0.1, that is acceptable. According to Table 4, the criterion of "being standard" had the highest relative weight (0.141). The relative weights (BNP) of other criteria can be seen in Table 3.
Results related to the closeness factor (the distanse from positive and negative ideal solutions) and prioritization of heat stress indices using fuzzy TOPSIS technique are seen in Table 5. According to findings, the first priority for surface mines was the WBGT index with closeness factor of Medium good (MG) (5,7,9) Good (G) (7,9,10) Very good (VG) (9,10,10) 0.730. The normalized fuzzy decision-making matrix has beeb presented in Table 6. Figure 3 indicates the prioritizing (ranking) of heat stress indices. All criteria of this study were positive.

Discussion
Heat stress is a hazardous physical agent in the workplaces. Nowadays, due to the technological advances, climate change, and consequently the increasing of global temperature, the heat stress is considered as a adverse factor in many workplaces. Exposure to heat causes heat syncope, heat exhaustion, heat stroke, confusion, poor concentration, and also exacerbation of some diseases, as well as reducing the production and increasing the indirect costs. Particularly, long-term exposure to hot conditions can result in so much stress on individuals [24]. The heat stress indices are used to assess this hazardous agent. These indices summarize the existing environmental conditions and express it as a number. In spite of many indices for assessment of heat stress and their advantages and disadvantages, there is no acceptable global comprehensive index for that. Therefore, it is necessary determining effective criteria to choose the most appropriate and precise heat stress index among existing ones [10]. In this study, the effective criteria for selecting a heat stress index were investigated. A group of the most important criteria, determined using the Delphi method and experts viewpoints, were simplicity, reliability, comprehensiveness, being direct reading, precision, strong correlation with the physiological strain indices, and so on. At first, the experts were asked to list the criteria affected on the selecting a heat stress index in order to extract the significant criteria. After some meetings with experts and omitting a number of criteria, 11 of 30 criteria were selected ( Table 2). The extracted criteria were next sent to experts to determine the importance of each of them by pair-wise comparisons. In the other word, using the pair-wise comparisons and fuzzy AHP was resulted in specifying the real weights of criteria. Among the obtained responses, the simplicity criterion had the highest frequency. Simplicity of an index means the calculations are feasible without complex computing and multiple devices. Although this criterion appropriated the highest percentage of responses (76.6%), its relative weight was 0.041 (the least weight after the criterion of "being lowcost") and had the less importance in comparison with other criteria. The next criterion with the most frequency was the reliability of index. Reliability refers to the fact that if the measurements are repeated several times in the same weather conditions, the similar results will obtain. The percentage of responses was 53.3%. The percentage for "being low-cost" the heat stress index was the same as the reliability (53.3%). The criteria of comprehensiveness, being direct  reading, and precision were located in the next priorities with the percentage of 50, 50, and 46.6%, respectively. The Fuzzy analytic hierarchy process (AHP) showed that according to the opinions of experts, the maximum weight was related to "being standard" of index. This criterion means that if the index is acceptable by international organizations or not? For instance, The International Organization for Standardization (ISO) publish the standards in many fields [10]. The criteria of precision, reliability, strong correlation with the physiological strain indices, and being unaffected by other factors had the higher weights, respectively (the priorities of secend to fifth). Although there is no study to determine the effective criterion for choosing a heat stress index, the National Institute for Occupational Safety and Health (NIOSH) has proposed some criteria for using an index. They are as follow: 1. It should be applicable and precise when using; 2. It should be consisted of all important factors (environmental, metabolism, and clothing ones); 3. Its measurements and calculations should be simple; 4. The means and methods used to measure the workers' exposure should not interfere with their activities; 5. The index should be applicable in the wide range of environmental and metabolic conditions [25].
It is worth mentioning that, the NIOSH recommendations to choose an appropriate index is aligned with the results of this study. So it can be concluded that to have a reliable and accurate measurements, the first step is choosing a suitable and precise index. To select an index based on the viewpoints of experts, its being standard is the most important criterion. The precision and reliability of index, and its strong correlation with the physiological strain indices should be considered. In this study for choosing an index among existing ones, the fuzzy TOPSIS technique was used. The TOPSIS algorithm is one of the most reliable scientific and management methods for decision-making. Using this technique, with regards to all effective aspects, it is possible assessing the criteria types, priorities, and weights toward each other, evaluating the indices relative to the criteria, and also prioritizing the indices in a resonable way. The results of closeness factor and indices priorities can be seen in Table 5. The results showed that for surface mines the WBGT index by the closeness factor of 0.730 appropriated the first priority among the existing indices. There are some important reasons for this result. This index by using the natural wet bulb temperature, globe and air temperatures can express the thermal conditions numerically. In addition to simplicity, WBGT has the high efficiency for assessing the thermal conditions. It is able to evalute the effects of heat during a period of time in which one performs an individual activity. It should be noted that the WBGT is the most widely used index around the world. By this index, there is an acceptable relationship with the physiological parameters at high temperatures [26]. The World Health Organization (WHO), International Organization for Standardization (ISO), and the National Institute of Occupational Safety and Health (NIOSH) have also recommended the WBGT as an index for assessing the heat stress. The Standard ISO-7247 has been introduced by ISO for measurement and evaluation of heat stress [10]. Alfano in 2014 conducted a review study on revision the WBGT index after 60 years of use.
It is said that, the WBGT index is simple to understand and use, and also it has suitable validity to organize the work-rest cycle. The fact that, there are some limitations related to WBGT and they would be preferable to consider for protecting the workers who work in hot environments [27]. According to Parsons (2006), the WBGT index, accepted by ISO, has adequate validity, releability, and applicability to control the heat  Table 6 The normalized fuzzy decision-making matrix using triangular fuzzy number (C1-C11: The criteria number 1 to 11)  Table 5, among the heat stress indices, the Predicted Heat Strain index (PHS) had the second priority with the closeness factor of 0.625. The PHS, the resultant of modifying the required sweat rate (SWreq) index, was proposed by ISO-7933 standard in 2004 [28]. This index was developed by a team of European researchers [29]. The PHS predicts the sweat rate and the internal core temperature for workers in hot environments. Although the PHS is not able to predict the response of each worker, it can predict the heat stress conditions by which the body core temperature increases. Consequently, it can estimate the maximum allowable exposure duration for work in hot workplaces. Because of complex calculations, this index is computed by a software. The required factors include air temperature, wet temperature, radiant temperature, relative humidity, air velocity, metabolic rate, partial vapor pressure, and thermal resistance of clothing [30]. Based on the fuzzy AHP and experts' viewpoints, among the criteria the highest weight was obtained for the "being standard" of predicted heat strain index. Being standard means that an index has been accepted by international organizations. It is consistent with the results of this study, because the WBGT and PHS indices have been recommended by ISO [27].
Holmer in 2010 conducted an investigation using internationally proposed methods including WBGT-ISO 7243 and PHS-ISO 7933. The results obtained from the comparative evaluation of heat stress at the same conditions with PHS and WBGT showed the WBGT index has more conservative assessment approach compared to PHS, resulting in much shorter working time. It seems that, WBGT overestimates the heat stress. This is aligned with findings related to some developing countries in which work is done under climatic conditions. In this case, considerable limitations are necessary when following the WBGT recommendations. Although, considering the data from warm countries, the predicted physiological strain by PHS may be more appropriate, the conservative philosophy of WBGT can be valuable in regard to safety [31]. Among the indices, the third priority was related to Thermal Work Limit (TWL) index with closeness factor of 0.551. It was proposed by Bates and Brake in 1997. The allowable limits of this index were firstly provided for underground industries (2002), and then these limits were proved for surface workplaces by Miller and Bates [32]. The TWL is a rational heat stress index that offers maximum tolerable rate of metabolism. So that, acclimatized person with suitable hydration status can work in the specific hot environment, provided that the body core temperature is less than 38.5°C and the sweat rate is less than 1.2 kg/h. This index is based on the experimental studies of human heat transfer and humidity and heat tranfer equations through clothing. The parameters of clothing are variable and this protocol can be used for unacclimatized people [33]. As an advantage of this index, there is no need for estimating the actual metabolic rate. The aim of the index is the calculating the maximum metabolic rate (W/m 2 ) that it can be continuously tolerated in hot environments. The maximum level of work in hot environments is calculated by measuring dry, wet, and radiation temperatures, wind speed, and atmospheric pressure, and also by considering the adaptability status and type of clothing. Also the safe work duration and guidelines for work-rest cycle in hot environments can be determined by this index. It is able to calculate the sweat rate and determine the needs for fluid intak for preventing dehydration. This index has taken the code of practice from the Abu Dhabi management system of health safety and environment to manage the heat stress [33]. Table 5 illustrates the forth priority of Universal Thermal Climate Index (UTCI) index with closeness factor of 0.543. As a new heat stress index for outdoor hot environments, the UTCI was developed in 2011 by international institute of biometeorology. It was extended by COST (a European Union program promoting Cooperation in Science and Technology) Action 730 [34]. It was derived from the scientific advances in the fields of human thermal physiology, biophysic, and heat transfer theories. It was proposed for developing appropriate standards to assess the heat stress in outdoor environments. Compared to other indices, it is sensitive to small changes in air temperature, wind, radiation, and humidity and it can be used in various climatic conditions [35]. Therefore, the UTCI index provides an appropriate assessment for impacts of climate changes on human. This index categorizes the thermal stresses in 10 classes, from extreme cold stress to extreme heat stress, as an equivalent temperature between −50 and +50 [36].
Among the indices, Environmental Stress Index (ESI) with the closeness coefficient of 0.517 had the fifth priority. Moran et al. (2001) proposed ESI based on the Fig. 3 Ranking the heat stress indices using the fuzzy TOPSIS measurement of air temperature, relative humidity (RH), and solar radiation (SR). These three variables are commonly used due to the quick response and simple measurement. Also, they can be determined using commercially available fast-reading sensors [37]. In addition to a high correlation between the ESI and WBGT index, the ESI correlation with the physiological parameters (core temperature, heart rate, and sweat rate) was considerable [38]. Following the assessment under various thermal conditions, ESI was proposed as a good alternative index for WBGT at the same conditions [39].
According to the results of studies conducted in mines of different countries, the WBGT, TWL, and effective temperature indices have been used seperately or in combination with each other, that it is consistent with selected indices in present study [40]. A study by Naghib indicated good correlation between WBGT and physiological parameters in surface mines, while this relationship was weak for undergroud mines [41].
To determine the optimal heat stress indices, Omidvari et al. evaluated the SWreq, WBGT, P4SR, and CET in surface mines located in Kermanshah province of Iran. In this respect, the criteria of heart rate and skin and core temperatures were also considered. The obtained results recommended WBGT and SWreq as the most appropriate indices to assess the heat stress [42]. Golbabaei et al. compared the indices DI, WBGT and SWreq in an outdoor industry with warm and humid conditions. Their investigation illustrated the WBGT index had the highest correlation with heart rate and the most optimal value compared to other indices. It was recommended for fast assessment of thermal stress in such thermal conditions, WBGT is often more valid index than the other indices [43].  compared different indices of heat stress to determine the best one. It was found that although there is considerable difference between the allowable limits of various indices, WBGT is one of the most appropriate indices to evaluate the thermal stress [44].
Since the activities in surface mines are performed under the sunlights, as a heat source, it can be concluded that in such conditions the WBGT index is more suitable than other indices and indicates the actual heat stress imposed by the environment. A study by Brief found that the WBGT index acts more accuratly than other indices in presence of radiation energy in the environment and it can indicated the actual thermal conditions. These findings are in line with results obtained in the present study.

Conclusion
In conclusion, this study demonstrates that WBGT is often the most appropriate index to evaluate the heat stress in surface mines. Some criteria may affect prioritizing the heat stress indices and choosing a suitable index among them. In this respect, the most important criteria can be extracted using experts' viewpoints and Delphi method. Simplicity, reliability, being low cost, and comprehensiveness are determinative criteria for applicability of an index. Determining the weights of criteria and prioritizing the indices are performed by the fuzzy AHP and TOPSIS techniques, respectively. The use of these strong methods allows introducing the most simple, precise, and applicable tool for evaluation the heat stress in hot environments. It seems that WBGT acts as an appropriate index for assessing the heat stress in mining activities at outdoors.

Aknowledgment
We are thankful to our colleagues, the professors of occupational health departments, post-graduated students, and the personnel of environmental and occupational health center in the ministry of health, who provided expertise that greatly assisted the research. The authors declare that there is no conflict of interests.

Funding
All parts of this research has been supported by Institute for Environmental Research (GrantNo.94-01-46-28540).

Availability of data and materials
The dataset supporting the conclusions of this article will not be shared due to performing our next projects with this software.
Authors' contributions MA participated in the acquisition, analysis, and interpretation of data and helped to draft the manuscript. PN supervised the study in all steps (acquisition, analysis, and interpretation of data). MMM acquisition, analysis, and interpretation, FG participated in the acquisition, and interpretation, HAA provided Fuzzy AHP-TOPSIS Method and has been consulted by AA, and AS provided geographical information. All authors read and approved the final manuscript.