Tra ﬃ c Accident Spatial Simulation Modeling for Planning of Road Emergency Services

: The appropriate locations of road emergency stations (RESs) can help to decrease the impact of tra ﬃ c accidents that cause around 50 million injuries per year worldwide. In this research, the appropriateness of existing RESs in the Khuzestan province, Iran, was assessed using an integrated fuzzy analytical hierarchy process (FAHP) and geographic information system (GIS) approach. The data used in this research were collected from di ﬀ erent sources, including the department of roads, the department of health, the statistics organization, forensics, police centers, the surveying and geological department, remotely-sensed and global positioning system (GPS) data of accident high crash zones. On the basis of previous studies and the requirements of the Ministry of Health and Medical Education, as well as the department of roads of Iran for the location of RESs, nine criteria and 19 sub-criteria were adopted, including population, safety, environmental indicators, compatible area in RES, incompatible area in RES, type of road, accident high crash zones, tra ﬃ c level and performance radius. The FAHP yielded the criteria weights and the ideal locations for establishing RESs using GIS analysis and aggregation functions. The resulting map matched the known road accident and high crash zones very well. The results indicated that the current RES stations are not distributed appropriately along the major roads of the Khuzestan province, and a re-arrangement is suggested. The ﬁnding of the present study can help decision-makers and authorities to achieve sustainable road safety in the case study area.


Introduction
Traffic accidents cause around 50 million injuries per year worldwide and are a significant concern for public health [1][2][3]. On the basis of findings by the World Health Organization (WHO), about 1.35 million people died worldwide because of road traffic accidents in 2016. Road traffic injuries are considered as one of the main reason for death the younger generation (5-29 years old) [4]. Although low-income and middle-income countries have only 50% of the world's vehicles, 90% of the world's fatalities on the roads occur in these countries [4]. Pre-hospital care has a significant role in handling and decreasing the traffic accident effects and saving the lives of people, and emergency medical services (EMS) play an indisputable role in this task [5,6]. EMS is a system that provides a variety

Case Study Area
Khuzestan Province, with an area of 63,633 square kilometers, is located in southwestern Iran ( Figure 2). This area is important in terms of the economy of the country, and for housing and industrial activities. According to the 2011 census, it consists of 24 counties, 61 cities and is home to 4,421,643 people. The sum total of the land routes of the province is 4727 kilometers, of which 2153 kilometers are main roads, and 2573 kilometers are secondary roads. The high number of recorded road accidents in this province is due to its topography, heavy traffic resulting from oil and gas industries and the fact that it is a center of trade and transportation in the southern part of Iran, as well as having access to pilgrimage sites and five roads connecting with neighboring countries.

Dataset
The present research is both descriptive and analytical in terms of its nature and methodology, and how it is applied in terms of the overall goal. Its results can be adapted in the programs and

Case Study Area
Khuzestan Province, with an area of 63,633 square kilometers, is located in southwestern Iran ( Figure 2). This area is important in terms of the economy of the country, and for housing and industrial activities. According to the 2011 census, it consists of 24 counties, 61 cities and is home to 4,421,643 people. The sum total of the land routes of the province is 4727 km, of which 2153 km are main roads, and 2573 km are secondary roads. The high number of recorded road accidents in this province is due to its topography, heavy traffic resulting from oil and gas industries and the fact that it is a center of trade and transportation in the southern part of Iran, as well as having access to pilgrimage sites and five roads connecting with neighboring countries.

Case Study Area
Khuzestan Province, with an area of 63,633 square kilometers, is located in southwestern Iran ( Figure 2). This area is important in terms of the economy of the country, and for housing and industrial activities. According to the 2011 census, it consists of 24 counties, 61 cities and is home to 4,421,643 people. The sum total of the land routes of the province is 4727 kilometers, of which 2153 kilometers are main roads, and 2573 kilometers are secondary roads. The high number of recorded road accidents in this province is due to its topography, heavy traffic resulting from oil and gas industries and the fact that it is a center of trade and transportation in the southern part of Iran, as well as having access to pilgrimage sites and five roads connecting with neighboring countries.

Dataset
The present research is both descriptive and analytical in terms of its nature and methodology, and how it is applied in terms of the overall goal. Its results can be adapted in the programs and

Dataset
The present research is both descriptive and analytical in terms of its nature and methodology, and how it is applied in terms of the overall goal. Its results can be adapted in the programs and policies of managers in road emergency departments. An extensive literature review was employed, along with interviews and questionnaires, to determine and extract the criteria needed for the research. The data used in this research were collected from different sources; traffic and road type data from the department of roads, the population information from the statistics organization, accident fatal data from forensics, performance radius of RESs and accident injuries data from the department of health and road accident data from police centers. Data for the incompatible area, compatible area and safety criteria were obtained from the surveying and geological department. Moreover, the topographic data were derived from a topographic 1:25,000 map and the shuttle radar topography mission (SRTM) digital elevation model (DEM) data with a spatial resolution of 30 m. We employed field observation data collected by global positioning system (GPS) to identify high crash zones and for validation purposes. After a preparation step, all datasets were imported into the GIS software and converted into raster format to form the basis of the considered criteria.

Methodology
The first step in the methodology was to recognize the criteria affecting the location of RESs. To this end, on the basis of previous studies [27,[34][35][36][37][38], several factors affecting the location of RESs were identified and submitted to 30 experts for initial ranking. The selected experts who represented the local and national university, the state RESs and EMS staff and operators filled the pairwise matrices of our study. They were asked to evaluate the importance of each criterion critically. Finally, nine casual criteria with one sub-criteria were specified to achieve the ultimate goal of the research. Then, the weights of each criterion, together with its sub-criteria, were computed using FAHP. In doing so, the final map of proposed sites for RESs was prepared by applying GIS-based aggregation functions. Figure 3 presents the main steps performed for the methodology. policies of managers in road emergency departments. An extensive literature review was employed, along with interviews and questionnaires, to determine and extract the criteria needed for the research. The data used in this research were collected from different sources; traffic and road type data from the department of roads, the population information from the statistics organization, accident fatal data from forensics, performance radius of RESs and accident injuries data from the department of health and road accident data from police centers. Data for the incompatible area, compatible area and safety criteria were obtained from the surveying and geological department. Moreover, the topographic data were derived from a topographic 1:25,000 map and the shuttle radar topography mission (SRTM) digital elevation model (DEM) data with a spatial resolution of 30 meters. We employed field observation data collected by global positioning system (GPS) to identify high crash zones and for validation purposes. After a preparation step, all datasets were imported into the GIS software and converted into raster format to form the basis of the considered criteria.

Methodology
The first step in the methodology was to recognize the criteria affecting the location of RESs. To this end, on the basis of previous studies [27,[34][35][36][37][38], several factors affecting the location of RESs were identified and submitted to 30 experts for initial ranking. The selected experts who represented the local and national university, the state RESs and EMS staff and operators filled the pairwise matrices of our study. They were asked to evaluate the importance of each criterion critically. Finally, nine casual criteria with one sub-criteria were specified to achieve the ultimate goal of the research. Then, the weights of each criterion, together with its sub-criteria, were computed using FAHP. In doing so, the final map of proposed sites for RESs was prepared by applying GIS-based aggregation functions. Figure 3 presents the main steps performed for the methodology.
. The fuzzy set theory (FST) mimics human judgment in that it deals with the uncertainty associated with measurement indicators [39]. In the Boolean logic, every element is considered as The fuzzy set theory (FST) mimics human judgment in that it deals with the uncertainty associated with measurement indicators [39]. In the Boolean logic, every element is considered as either true or false. However, in fuzzy logic, elements can be considered as both true and false within different degrees. Thus, each element can be in a set by a specific degree. The FST is practical when it is difficult to extract precise decisions by human judgment. Decision-makers and policymakers also desire natural language expressions instead of crisp numbers [40]. The FST has been a significant step towards representing vague and incomplete data because of its ability to deliver a methodology for computing directly with words. The FST makes it possible to use three key factors, namely FST, membership function and fuzzy numbers, to change vague data into practical data efficiently. Fuzzy numbers are natural generalizations of ordinary values. The triangular and trapezoidal fuzzy numbers are usually used to capture the vagueness of the parameters related to selecting the alternatives [26]. In the present study, we used triangular and trapezoidal fuzzy numbers to prioritize the suitability of land for RESs. These numbers are expressed with boundaries instead of crisp numbers to reflect the fuzziness, as decision-makers select the alternatives in the pairwise comparison matrices.

Integration Approach of Fuzzy AHP (Analytical Hierarchy Process)
The AHP method was proposed by Thomas L. Saaty (1980) [42]. This method is a widely used multiple criteria decision-making approach in GIS spatial analysis [43]. Although the AHP method aims to capture experts' knowledge, the conventional AHP still cannot take into account the human thinking style. The traditional AHP method has some limitations by using the exact values to express the decision-makers' judgments in a comparison of alternatives [25,44]. As the comparison of indicators in the traditional AHP method is based on expert judgments, some degree of inconsistency may transfer from uncertain comparisons to the results [39][40][41][42][43]45].
Moreover, experts are usually from different fields and may not be entirely sure about the

Integration Approach of Fuzzy AHP (Analytical Hierarchy Process)
The AHP method was proposed by Thomas L. Saaty (1980) [42]. This method is a widely used multiple criteria decision-making approach in GIS spatial analysis [43]. Although the AHP method aims to capture experts' knowledge, the conventional AHP still cannot take into account the human thinking style. The traditional AHP method has some limitations by using the exact values to express the decision-makers' judgments in a comparison of alternatives [25,44]. As the comparison of indicators in the traditional AHP method is based on expert judgments, some degree of inconsistency may transfer from uncertain comparisons to the results [39][40][41][42][43]45].
Moreover, experts are usually from different fields and may not be entirely sure about the importance of some of the indicators and their comparisons [25]. The inconsistency may be exacerbated when experts and decision-makers compare a wide range of factors related to different matters. Comparisons of decision-makers are the primary source to weigh the factors based on their importance regarding the specific problem. These comparisons are represented qualitatively by different linguistic variables. The relation among the quantitative values and linguistic variables are represented in Table 1. Table 1. Triangular fuzzy number of linguistic variables.

Linguistic Variables Triangular Fuzzy Numbers Reciprocal Triangular Fuzzy Numbers
Extremely strong (9,9,9) (1/9, 1/9, 1/9) Very strong ( The AHP method is often criticized because of its use of an unbalanced scale of comparisons and its limitation to adequately handle the associated uncertainty in the pairwise comparison process [45]. To overcome these shortcomings, the fuzzy analytical hierarchy process was developed for solving hierarchical problems. Most decision-makers frequently find that it is more reliable to give interval judgments than fixed value comparisons. This is because they are often unable to express explicit preferences [46]. Thus, fuzzy AHP is applied to solve hierarchical fuzzy problems. The integration of the AHP with the fuzzy set theory is helpful to reduce the subjectivity in the factor weightings. Therefore, fuzzy numbers were specified to make comparative pairwise and structuring of the resultant matrices. Let A be the pairwise comparison matrix that is prepared with the fuzzy numbers [41,47], then where each component of this matrix follows the condition that if a = l ij , m ij , u ij then a −1 ij = 1 l ij , 1/m ij , 1/u ij for 1,j, . . . , n and i j.
The phases of the extent analysis approach can be summarized as follows: where in this case d is the ordinate of the highest intersection point between µM1 and µM2. To perform this computation, the values of V (S j ≥ S i ) and V (S i ≥ S j ) with S i and S j have to be compared. As S j = (l j , m j , u j ) and S i = (l i , m i , u i ) are considered as two triangular fuzzy numbers ( Figure 5), the possibility degree of S j = (l j , m j , u j ) [ ] where in this case d is the ordinate of the highest intersection point between μM1 and μM2. To perform this computation, the values of V (Sj ≥ Si) and V (Si ≥ Sj) with Si and Sj have to be compared.
As Sj = (lj, mj, uj) and Si = (li, mi, ui) are considered as two triangular fuzzy numbers ( Figure 5), the possibility degree of Sj = (lj, mj, uj) ≥ S i= (li, mi, ui) is defined as For the possible degree of a convex fuzzy number to be larger than k, convex fuzzy Mi (I = 1,2, ...,k) numbers can be defined by Now we can have d(A i ) = min V (S i ≥ S k ) for k = 1, 2, . . . , n. At that point, the weight vector is given by Via normalization, the normalized weight vectors are where W is a non-fuzzy number.

Standardization of the Criteria Based on TFNs:
In most cases, after the criteria and the sub-criteria are determined they require standardization, since the criteria cannot be compared. Data standardization is a process in which the criteria are converted and re-scaled to be comparable. The common method for the standardization of the criteria is based on the FST. In the fuzzy operation, the inputs are captured and, using the corresponding membership functions such as a linear function, an appropriate grade is attributed to each [47]. However, there are other factors that increase the ability of the FST to be applied for the standardization of the criteria, such as uncertainty in decision-making and imprecise data [48]. Following the preparation of the desired layers, it was necessary to convert the criteria into raster format to apply the membership functions and standardize them. To this end, using the Euclidean distance function for continuous data as well as discrete data, each criterion was given a score in the range of zero to one based on the opinions of the experts. They were then converted into a raster layer.

Implementation of the Workflow
Determining the effective factors is considered one of the most important and fundamental steps in determining the locations of RESs. Therefore, the nine criteria and the 19 sub-criteria used in the present research were adopted based on previous studies and the requirements of the Ministry of Health and Medical Education, as well as the department of roads of Iran for the location of RESs. These nine criteria are explained in the following list: 1.
Population: The density of the population close to the road is a major factor affecting the location of RESs. Expansion of residential area may increase the road accident risk [13].

2.
Safety: The safety of the RESs is related to their proximity to the road outposts and their distance from faults or watercourses. This is because certain phenomena, such as faults, may cause instability, and proximity to police stations can be effective in creating order and security when necessary. 3.
Natural Factors: Natural factors, such as slope elevation and topography of the region, influence the location and ability to establish RESs. Moreover, road slope can have an effect on the road traffic level where drivers should be careful about their speed because of the high level of slope [49].

4.
The Incompatible Area in RESs: This term refers to the forms of land use that damages RESs when they are in close proximity with one-another, harming the mental health of people. They include gas stations and residential use, and the privacy of transmission lines of electricity, gas and phone. 5.
The Compatible Area in RESs: This term refers to forms of land use that deliver favorable results in the RESs. Proximity to Red Crescent centers and traffic police lead to synergies in accomplishing traffic missions where the presence of relief agencies beside one-another facilitates providing services and rescuing the injured. Furthermore, barren lands provide open space for relief helicopters to land, in addition to facilitating RES construction. 6.
Road Type: The road type is of prime importance in this kind of location analysis because roads with two-way crossings incur more accidents than expressways. It was reported that the number of road accidents on two-way crossings is higher than expressways [50]. 7.
Traffic Level: This is one of the important factors affecting the location of RESs. An increase in traffic leads to a growth in the number of accidents. 8.
High Crash Zones: According to the conducted studies, accidents are more likely to occur in certain parts of the roads, and these are known as high crash zones or road accident blackspots. The proximity of the RESs to blackspots is a crucial factor.
In the present research, the number of intensity-methods of accidents was employed to evaluate the high crash zones. In this method, an index, known as coefficient P, was determined for every position. This index is obtained from the following equation (from the Iran Road Maintenance Organization): The P index is obtained using the following formula: According to this equation: P represents the high crash zones, X shows the number of fatal accidents, Y indicates the number of accidents causing injury and Z is the number of accidents causing only damage to property. 9. Distance (or performance radius of existing RESs): According to the existing regulations (the organization regulation for the comprehensive coverage of pre-hospital emergency medical services), RESs should be established every 20 km or every 40 km at most [51][52][53]. In addition, the distance from existing stations should not be less than 15 km. The distance from the villages to the nearest RES should not exceed 40 km [53,54]. The Ministry of Health and Medical Education is required to provide emergency medical services to the villagers and nomads across the country through emergency stations in rural health centers that are more than 40 km from the pre-hospital emergency centers and are not located on the main roads of the country.

Results
From 2011 until 2017, on the basis of the data obtained from the Legal Medicine Organization and the emergency centers of the Khuzestan province, the number of the reported road accidents in this province was 6303, causing 877 deaths and 1210 injuries ( Table 2). The roads from Ahwaz to Andimeshk and Ahwaz to Ramhormoz had the highest death rate among all the roads, with 141 and 137, respectively. The lowest death rate of 14 was recorded on the road from Susangerd to Andimeshk. In this period, the number of EMS call-outs was 9755. The highest number of call-outs occurred in March, which coincides with the Iranian New Year holidays. The highest number of call-outs was at the Andimeshk and Ahwaz road emergency stations, with a total of 1539. The FAHP found that high crash zones were the primary factor among the nine casual criteria with a value of 0.279. Table 3 shows the weight of each criterion affecting the location of RESs in the Khuzestan province. The calculations and analysis of the achieved statistics identified only 40 high crash zones from among 400 sections (Table 2).  Table 3. Evaluation criteria of road emergency stations (RES) and fuzzy analytical hierarchy process (AHP) weights. After determining the weights, fuzzy maps ( Figure 6) were prepared for each of the criteria and sub-criteria, and the following results were obtained.  Finally, one the basis of the filled pairwise comparison matrices by the experts who were mentioned in the methodology section, all prepared fuzzy maps of the conditioning factors were overlaid with the respective FAHP weights in Arc GIS software. Thus, the final map was generated using the fuzzy maps for each factor and the corresponding weight for proposed points for establishing RESs (Figure 7). Finally, one the basis of the filled pairwise comparison matrices by the experts who were mentioned in the methodology section, all prepared fuzzy maps of the conditioning factors were overlaid with the respective FAHP weights in Arc GIS software. Thus, the final map was generated using the fuzzy maps for each factor and the corresponding weight for proposed points for establishing RESs (Figure 7).

Discussion
According to the accident statistics provided by the traffic police of the Khuzestan province during the 2011-2017 period, it was noted that 40 high crash zones exist in the connecting axes of the province's road network. In total, most of the accidents occurred in the axes of Ahwaz-Andimeshk, Ahwaz-Izeh and Abadan-Mahshahr during the 6 years. The highest number of fatal accidents occurred on the main roads leading to Ahwaz, Andimeshk and Abadan. At an increasing distance from Ahvaz city, which is the capital of the province, the density of fatal accidents decreases. This requires the serious attention of authorities to these critical areas. The obtained information yielded that the number of deaths in the categories of drivers/passengers of busses and drivers/passengers of heavy trucks in the Khuzestan province is higher than the average of the country. In contrast, in the categories of drivers/passengers of four-wheeled cars and pedestrians have lower shares (Figure 8). As illustrated in Figure 9, the number of reported road accidents witnessed a moderate increase; it reached to 1098 from 794. The number of accident injuries almost had the same trend; it increased from 150 to 187. In contrast, the number of accident deaths slightly reduced (from 135 to 125).

Whole country (Iran) (%)
Khuzestan Province (%) Figure 7. The resulting map of proposed points for establishing RES.

Discussion
According to the accident statistics provided by the traffic police of the Khuzestan province during the 2011-2017 period, it was noted that 40 high crash zones exist in the connecting axes of the province's road network. In total, most of the accidents occurred in the axes of Ahwaz-Andimeshk, Ahwaz-Izeh and Abadan-Mahshahr during the 6 years. The highest number of fatal accidents occurred on the main roads leading to Ahwaz, Andimeshk and Abadan. At an increasing distance from Ahvaz city, which is the capital of the province, the density of fatal accidents decreases. This requires the serious attention of authorities to these critical areas. The obtained information yielded that the number of deaths in the categories of drivers/passengers of busses and drivers/passengers of heavy trucks in the Khuzestan province is higher than the average of the country. In contrast, in the categories of drivers/passengers of four-wheeled cars and pedestrians have lower shares (Figure 8).

Discussion
According to the accident statistics provided by the traffic police of the Khuzestan province during the 2011-2017 period, it was noted that 40 high crash zones exist in the connecting axes of the province's road network. In total, most of the accidents occurred in the axes of Ahwaz-Andimeshk, Ahwaz-Izeh and Abadan-Mahshahr during the 6 years. The highest number of fatal accidents occurred on the main roads leading to Ahwaz, Andimeshk and Abadan. At an increasing distance from Ahvaz city, which is the capital of the province, the density of fatal accidents decreases. This requires the serious attention of authorities to these critical areas. The obtained information yielded that the number of deaths in the categories of drivers/passengers of busses and drivers/passengers of heavy trucks in the Khuzestan province is higher than the average of the country. In contrast, in the categories of drivers/passengers of four-wheeled cars and pedestrians have lower shares (Figure 8). As illustrated in Figure 9, the number of reported road accidents witnessed a moderate increase; it reached to 1098 from 794. The number of accident injuries almost had the same trend; it increased from 150 to 187. In contrast, the number of accident deaths slightly reduced (from 135 to 125).

Whole country (Iran) (%)
Khuzestan Province (%) Figure 8. Comparison between the Khuzestan province and the whole of Iran in the rate of deaths by road user category.
As illustrated in Figure 9, the number of reported road accidents witnessed a moderate increase; it reached to 1098 from 794. The number of accident injuries almost had the same trend; it increased from 150 to 187. In contrast, the number of accident deaths slightly reduced (from 135 to 125). Three elements, namely the driver, the road and the vehicle, have the most significant impact on the occurrence of road accidents in this region, based on the data obtained from the police center of the Khuzestan province. Because of their errors and interaction with other causes, for example, vehicles or roadway design, drivers are often the cause of crashes [55].
Furthermore, an inappropriate road network can increase the crash risk, as it determines how road users perceive their environment, meaning that the roadway offers instructions to the drivers on what they should be doing. Negative road engineering features include those where a road defect directly triggers a crash, or where some factor of the road environment misleads a driver and thereby creates human errors [56,57]. The results obtained in several recent studies [58][59][60] illustrate that human behaviors are the cause of accidents in almost 70% of cases (due to a lack of adherence to the laws and regulations, disregarding warning signs, a low level of education and a lack of familiarity with the minimum technical requirements of cars). In 15% of the cases, the cause of the accident is attributable to the structure of the road (for example, the shape of the bends, road surface slope and mountain passes). In a further 15% of cases, the cause is related to the vehicles (due to depreciation and malfunctioning). On the basis of objective evidence and the accident reports from the traffic police and the emergency department, the statistics indicate the cause of accidents in the Khuzestan province can be attributed to the driver, the road and the vehicle in 75%, 15% and 10% of cases, respectively.
After comparing the final map with the other existing maps of the practical factors, it was revealed that the ultimate map matched well with the maps of road traffic and high crash zones. It seems that areas with a high density of fatal accidents have a lack of RESs and, therefore, a distinct need for the establishment of new RESs ( Figure 10). This matter confirms the excellent performance of the FAHP method and the GIS in the identification of locations for the establishment of RESs. The results are consistent with those from Goli et al., 2015, andOhta et al., 2007 [36,61]. Three elements, namely the driver, the road and the vehicle, have the most significant impact on the occurrence of road accidents in this region, based on the data obtained from the police center of the Khuzestan province. Because of their errors and interaction with other causes, for example, vehicles or roadway design, drivers are often the cause of crashes [55].
Furthermore, an inappropriate road network can increase the crash risk, as it determines how road users perceive their environment, meaning that the roadway offers instructions to the drivers on what they should be doing. Negative road engineering features include those where a road defect directly triggers a crash, or where some factor of the road environment misleads a driver and thereby creates human errors [56,57]. The results obtained in several recent studies [58][59][60] illustrate that human behaviors are the cause of accidents in almost 70% of cases (due to a lack of adherence to the laws and regulations, disregarding warning signs, a low level of education and a lack of familiarity with the minimum technical requirements of cars). In 15% of the cases, the cause of the accident is attributable to the structure of the road (for example, the shape of the bends, road surface slope and mountain passes). In a further 15% of cases, the cause is related to the vehicles (due to depreciation and malfunctioning). On the basis of objective evidence and the accident reports from the traffic police and the emergency department, the statistics indicate the cause of accidents in the Khuzestan province can be attributed to the driver, the road and the vehicle in 75%, 15% and 10% of cases, respectively.
After comparing the final map with the other existing maps of the practical factors, it was revealed that the ultimate map matched well with the maps of road traffic and high crash zones. It seems that areas with a high density of fatal accidents have a lack of RESs and, therefore, a distinct need for the establishment of new RESs ( Figure 10). This matter confirms the excellent performance of the FAHP method and the GIS in the identification of locations for the establishment of RESs. The results are consistent with those from Goli et al., 2015, andOhta et al., 2007 [36,61].

Conclusions
In this study, the appropriateness of existing RESs in the Khuzestan province, Iran, was assessed using an integrated FAHP and GIS. One of the core purposes of adopting GIS in different applications is to make a support system for spatial decision-making because, in most cases, spatial issues are semi-structured. Optimum site selection is considered as one of these spatial problems. On the basis of the last reports of the Iranian road safety commission, as well as successful experiences from other countries, the following are some suggestions to decrease the number of road traffic deaths and injuries in the Khuzestan province: • According to the respondents, the highest impact on reducing the number of fatal injuries on the connecting roads of the province can be achieved by proper quality equipment of rescue posts and qualified medical staff, and it is thus proposed that these posts be equipped by the Ministry of Health and Medical Education immediately.
• The healthcare organization of the Khuzestan province should give priority to the standardization of the equipment and facilities of RESs.
• Improving the situation of media, which has an incredible impact on society, is considered as a powerful tool to provide people with awareness about road accident troubles.
• Using new technologies to record and analyze big data, as well as upgrading the existing roads to record violations, would effectively reduce the number of road violations by the drivers and, thus, lead to fewer road accidents.
• Conducting stricter laws on vehicle safety standard certificates, as the current situation of vehicles in Iran-of which most are made in Iran-is not acceptable, and they still require improvement in the safety topic.
• Improving the situation of roads in the Khuzestan province, as more than 80% of roads in this province are two-way crossing roads, with the province having heavy traffic.
The number of accidents causing death and injuries in the Khuzestan province is high enough to demand attention from decision-makers and authorities. Moreover, the appropriate locations of RESs not only can help to decrease the humanitarian aspects of road accident death and injuries, but also can reduce road crash costs on the economy. To conclude, the current RES stations are not

Conclusions
In this study, the appropriateness of existing RESs in the Khuzestan province, Iran, was assessed using an integrated FAHP and GIS. One of the core purposes of adopting GIS in different applications is to make a support system for spatial decision-making because, in most cases, spatial issues are semi-structured. Optimum site selection is considered as one of these spatial problems. On the basis of the last reports of the Iranian road safety commission, as well as successful experiences from other countries, the following are some suggestions to decrease the number of road traffic deaths and injuries in the Khuzestan province:

•
According to the respondents, the highest impact on reducing the number of fatal injuries on the connecting roads of the province can be achieved by proper quality equipment of rescue posts and qualified medical staff, and it is thus proposed that these posts be equipped by the Ministry of Health and Medical Education immediately.

•
The healthcare organization of the Khuzestan province should give priority to the standardization of the equipment and facilities of RESs. • Improving the situation of media, which has an incredible impact on society, is considered as a powerful tool to provide people with awareness about road accident troubles. • Using new technologies to record and analyze big data, as well as upgrading the existing roads to record violations, would effectively reduce the number of road violations by the drivers and, thus, lead to fewer road accidents.

•
Conducting stricter laws on vehicle safety standard certificates, as the current situation of vehicles in Iran-of which most are made in Iran-is not acceptable, and they still require improvement in the safety topic. • Improving the situation of roads in the Khuzestan province, as more than 80% of roads in this province are two-way crossing roads, with the province having heavy traffic.
The number of accidents causing death and injuries in the Khuzestan province is high enough to demand attention from decision-makers and authorities. Moreover, the appropriate locations of RESs not only can help to decrease the humanitarian aspects of road accident death and injuries, but also can reduce road crash costs on the economy. To conclude, the current RES stations are not distributed appropriately along the major roads of the Khuzestan province, and a re-arrangement is suggested.

Conflicts of Interest:
The authors declare no conflict of interest.