# Mathematical Models in Humanitarian Supply Chain Management: A Systematic Literature Review

Authors: Muhammad Salman Habib; Young Hae Lee; Muhammad Saad Memon
Journal: Mathematical Problems in Engineering
 (2016)

---

## Abstract

In the past decade the humanitarian supply chain (HSC) has attracted the attention of researchers due to the increasing frequency of disasters. The uncertainty in time, location, and severity of disaster during predisaster phase and poor conditions of available infrastructure during postdisaster phase make HSC operations difficult to handle. In order to overcome the difficulties during these phases, we need to assure that HSC operations are designed in an efficient manner to minimize human and economic losses. In the recent times, several mathematical optimization techniques and algorithms have been developed to increase the efficiency of HSC operations. These techniques and algorithms developed for the field of HSC motivate the need of a systematic literature review. Owing to the importance of mathematical modelling techniques, this paper presents the review of the mathematical contributions made in the last decade in the field of HSC. A systematic literature review methodology is used for this paper due to its transparent procedure. There are two objectives of this study: the first one is to conduct an up-to-date survey of mathematical models developed in HSC area and the second one is to highlight the potential research areas which require attention of the researchers.


---

## Body

## 1. Introduction

The humanitarian supply chain (HSC) has become an important issue for academia and professionals since the Asian tsunami which occurred in 2004. After the disaster of the 2004 tsunami in Indian Ocean, excessive relief goods blocked the airports and warehouses in the affected regions. Humanitarian relief providing agencies had to struggle a lot to sort out required goods and distribute them in timely and economic way. Before the disaster of the 2004 tsunami, the HSC was not considered as important as it is today and most of the work in the field of disaster relief before 2004 was in the context of the commercial supply chain. Usually it is thought that natural disasters are low frequency and high consequence incidents, yet there are some parts of the world that are hit by disaster several times a year. For example, Central America and the Caribbean are frequent victims of hurricanes. Similarly the central plains of USA, named as “Tornado Alley,” are repeatedly hit by a number of tornados. Japan has a long history of earthquakes because it lies in the region of tectonic plates called “the Ring of Fire.”Disasters can also be caused by human beings. Technological disasters that include chemical spill, radioactive radiation, road, air, maritime accidents, and groundwater contamination have also caused significant property damage and loss of life. Japan’s Fukushima nuclear power station failure after earthquake in 2011 and Ukraine’s Chernobyl nuclear disaster on April 26, 1986, are the examples of horrible accidents of modern human times. Even after four years the Fukushima nuclear power plant is still producing highly toxic water.Figure1 shows that the number of disasters is increasing over the last 55 years. According to Thomas and Kopczak [1] the market for disaster relief will keep on growing due to increasing number of disasters. It is estimated that the frequency of occurrence of these disasters will increase five times in the next 50 years [2]. Keeping in view this forecast, delivering humanitarian aid can be perceived as a significant future global industry.Figure 1
Reported natural disasters between 1960 and 2015. Source: EM-DAT (the International Disaster Database).The humanitarian supply chain is defined as follows:… the process of evacuating people from disaster stricken areas to safe places and planning, implementing and controlling the efficient, cost-effective flow of goods, meanwhile collecting related information from the point of supply to the point of consumption for the purpose of alleviating the sufferings of vulnerable people.
The conceptual framework of HSC is depicted in Figure2.Figure 2
Conceptual framework of humanitarian supply chain.In the last decade, numerous studies have been published in the field of the HSC. They have addressed this topic by different titles. However, their ultimate objective is to help affected people of disaster stricken areas. One can find the research of HSC by different names likedisaster relief operation, disaster relief supply chain, emergency relief operation, andemergency management. A humanitarian operation for a disaster can be divided into four different phases, namely, mitigation, preparedness, response, and recovery [2]. The predisaster phase covers the mitigation and preparedness phase. Mitigation includes the steps to reduce vulnerability to disaster impact such as injuries and loss of life and property, while preparedness includes educating communities on how a disaster can affect them so that they can adopt a proactive approach. The postdisaster phase covers the response and recovery phases. The response phase addresses immediate threats to minimize economic and human losses, while the recovery phase supports the restoration of all the damage caused by the disaster.Research related to the literature reviews in the field of HSC is done under different names with different aims. The research of Altay and Green III [2] focuses on disaster operation life cycle. They explained all the activities required to execute in each phase of the disaster. Simpson and Hancock [3] reviewed the implementation of operation research techniques in the field of disaster response for the past 50 years by providing a detailed network of citations. Natarajarathinam et al. [4] provided insights to manage the supply chain in the times of crisis and proposed a five-dimensional framework to classify the literature. Caunhye et al. [5] reviewed the optimization models in the area of facility location, relief distribution, and casualty transportation. Dasaklis et al. [6] focused on the role of logistics in HSC to control epidemic outbreaks. John et al. [7] divided the HSC in predisaster and postdisaster phases and explained how postdisaster HSC operational issues can be treated analogous to the operational and tactical issues in commercial supply chain. Galindo and Batta [8] extended the article of Altay and Green III [2] with the new advancements of OR/MS in disaster operation management. Abidi et al. [9] provided the systematic literature review for the performance measurements in the HSC by proposing performance measurement directions with input and output criteria. They also explained how to design, deploy, and disseminate performance measurement and humanitarian supply chains. Özdamar and Ertem [10] provided a review about the models for response and recovery phase of disaster only. Anaya-Arenas et al. [11] provided a review of research about relief distribution networks in the HSC by categorizing them according to objective function, model constraints, and solution methodology. Hoyos et al. [12] provided a review of the work in the field of HSC that used the OR models with stochastic part. Zheng et al. [13] provided a survey of evolutionary optimization techniques being used for disaster relief operations. The above reviews focus on a certain phase of disaster or a certain type of modelling technique for HSC. Currently, the HSC is lacking a literature review that is holistic in approach, covering all the phases of disaster, and provides a summary of modelling techniques and solution methodologies.The rest of the paper is organized as follows: Section2 explains the methodology of the systematic literature review. Section 3 provides a detailed review of mathematical models published in the last decade in the field of the HSC. Section 4 targets our research questions and provides the details of mathematical modelling techniques and solution methodologies used in HSC and points out the unexplored areas in the HSC to provide an agenda for the future research. Finally, the conclusion is presented in Section 5.


## 2. Methodology of Systematic Literature Review

A review is called systematic if it is based on clearly formulated questions, relevant studies, evaluated quality, and the synthesized results. It is the categorical and organized approach that differentiates a systematic review from a traditional review. The purpose of a systematic review is to summarize the best available research on a particular topic. It follows a transparent procedure to collect, analyze, and synthesize the results of relevant research. The whole process is explicitly defined in order to maintain the transparency [14]. A systematic review methodology consists of four steps: planning, searching, screening, and extraction. In the following section we discussed how these four steps have been implemented in this research.Planning. In the planning phase research questions are framed. The questions should be designed in such a way that the problem to be addressed is specified in the form of clear and unambiguous questions [14]. In this study, the following research questions are developed:RQ 1. What is the current status of research on the mathematical models in the field of HSC?
RQ 2. What are the unexplored areas in the field of HSC?Searching. The key terms to collect the research papers related to HSC are developed based on the above research questions. In early 2015, the papers are collected which are focusing on the following keywords: “disaster response model,” “disaster relief operation,” “emergency response,” “emergency relief operations,” “humanitarian supply chain,” and “humanitarian relief operations.” In order to find more precise research papers, a search using a Boolean connector (AND, OR, AND NOT) was also conducted. This search was done in four research databases: the Web of Science, Scopus, Elsevier, and Google Scholar. The period of publication we determined was from 2005 to 2015, because the majority of work in the perspective of humanitarian supply chain was published after 2004 when tsunami in the Indian Ocean hit many countries. Our search process was able to collect 1,487 articles, which, after removing duplicates and research papers from the commercial supply chain, left 452 articles.Screening. To make sure of the objectivity of this research, it was necessary to define inclusion and exclusion criteria. We determined the following inclusion and exclusion criteria.Inclusion criteria: As the objective of this paper focuses on the mathematical modelling techniques in the field of HSC, only those papers are included which proposed any type of mathematical technique. We included only peer reviewed journal articles, except a few conference papers as they were important to answer the above developed research questions.
Exclusion criteria: articles written in the perspective of healthcare management field and mathematical models developed in the perspective of commercial supply chains are out of the scope of this paper.From 452 articles, total 140 articles are selected based on the following two points: whether the article’s theme is HSC and whether the article uses any mathematical technique. Among these 140 articles, 94 articles were finally selected to provide significant insights into developed research questions. Figure3 illustrates the systematic screening process.Figure 3
Paper screening methodology (modified from Gimenez and Tachizawa [99]).Extraction. In the extraction phase the selected papers were divided according to the important operations of HSC. These selected papers are discussed in the Discussion and Implications section of this research paper.


## 3. Research in the Humanitarian Supply Chain

Before going into the details of HSC research, it is necessary to have a brief look at all the processes of the HSC. Disaster relief operation starts with an alert warning announced by disaster management authorities. An alert warning is not always a part of the HSC because some disasters like earthquakes and bomb blasts occur so instantly that it is impossible to issue an alert warning. After the alert, the mitigation and preparedness phases start. In the mitigation phase, measures are taken to reduce the severity of the disaster while in the preparedness phase relief distribution centers are located and people are shifted to safe places. After occurrence of the disaster, the next task is to assess damage caused by the disaster. The damaged transportation links are repaired and network is designed to distribute relief goods to the affected people following the completion of assessment process. After relief goods are dispatched from a central distribution center, the next task is to track the deliveries until they are delivered to the disaster stricken areas. Finally, feedback is collected from all regional distribution centers and relief goods are dispatched according to the updated demand information. The phases of the HSC described here are not as simple as they appear, because poor conditions of communication systems and infrastructure, lack of resources, and uncertainty in terms of time and severity of the disaster make the disaster relief operation much more complex. Research in the field of HSC is divided into three major categories:(i)
facility location,
(ii)
network design and relief distribution,
(iii)
mass evacuation.HSC operations are classified in these three areas because most of the research about HSC operations comes under the umbrella of these three categories. All other processes of the HSC have been discussed under these categories. The category of facility location covers the mathematical models made for resource allocation, predisaster facility location, postdisaster facility location, location-allocation, and maximal covering models. The category of network design and relief distribution covers the relief distribution planning, vehicle routing, assessment routing, causality transportation, international relief distribution, and location-routing models. The category of mass evacuation covers traffic control planning and mass evacuation models in pre- and postdisaster scenarios. In the next section, research on each category is discussed in detail.


### 3.1. Facility Location in Humanitarian Supply Chain

After the occurrence of a disaster when assessment of damage is complete, next step is to provide relief to the affected people. For this purpose, a relief distribution channel is designed which includes central distribution centers, warehouses, and regional relief distribution centers. The challenge in this process is to locate these facilities in such a way that demand of entire disaster stricken area is fulfilled with minimal delivering cost and maximum service level. Factors that increase the complexity of facility location task are as follows:(i)
Very short lead time and sudden surge in demand just after the disaster.
(ii)
Uncertainty in timing and scale of the disaster.
(iii)
Very high stakes associated with the timely delivery of relief goods.
(iv)
Damaged infrastructure.Research in the section of facility location for a disaster can be classified into two categories. The classification is with respect to the timeline of the disaster, whether the model is developed for a predisaster scenario or postdisaster scenario.In the predisaster scenario, a strategic approach is adopted in which researchers plan for the best locations of the facilities in order to distribute relief goods at minimal cost and maximum service level. Some authors along with the facility location planning also consider allocation of resources from central distribution centers to warehouses and regional distribution centers. Such problems are called location-allocation problems. For instance, Dekle et al. [15] constructed a model to locate the disaster recovery centers in the predisaster context and maximized the covering of relief goods in the disaster target area. Chang et al. [16] developed a model for resource distribution in urban flood disaster scenario. The authors considered the location of the warehouse, prioritization of facility allocation, shortage and penalties for surplus with the objective function of minimizing transportation cost, facility setup cost, and transportation cost of rescue equipment. McCall [17] presented a model with the prepositioning of assistance pack-up kits during disaster while considering the constraints of model facility capacity, number of kits for prepositioning before disaster, and unsatisfied demands with the objective function of minimizing victim nautical miles and shortages. Akgün et al. [18] adopted a different approach for facility location model and suggested that prepositioning supplies in the preparation phase is desirable to locate supplies near the expected disaster area so that aid can be distributed within a minimal time. On the other hand it is also important that if distribution centers are close to expected disaster area, then these distribution centers may get destroyed. The goal of the study is to choose locations so that reliability can be maximized. Just after the occurrence of a disaster, there is an urgent need of basic emergency supplies and this surge in demand continues for the first 72 hours. Keeping this point in view, Rawls and Turnquist [19] and Lin et al. [20] developed a model for the short term disaster response.On the other side of the disaster timeline there are some researchers who consider the postdisaster scenario for the facility location problem. In this postdisaster scenario, it is assumed that, after the occurrence of disaster, the assessment of the damage has been completed and all the major information required for a relief operation is available. For instance, Balcik and Beamon [21] suggested a facility location problem in the postdisaster scenario. The problem is a variant of maximal covering location model with budget and capacity constraints. They determined a number of locations of distribution centers and the amount of relief goods at each center to meet the demand. Bozorgi-Amiri et al. [22] presented a stochastic model that first determines the location of aid distribution centers and then allocates relief goods to the affected area with an objective function of minimizing the costs for predisaster setup that includes costs of procurement, transportation, holding, and shortage. Horner and Downs [23] proposed a model that is a variant of warehouse location problem in which they used warehouses to allocate the relief goods to the affected areas. In fact, it would be better to consider them as an intermediate distribution facility model with the objective function of minimizing the costs of distributing of relief goods. Zhang et al. [24] developed a resource allocation model while considering the constraints of multiresources and multidepots, with the objective function of minimizing the cost of the total time for dispatching emergency resources. Hong et al. [25] made a comparison of two robust mathematical models. The authors compared the robust integer facility location with the robust continuous facility location model and concluded that robust integer facility location works best in normal conditions, while the robust continuous facility location works best in the hours of disaster. Abounacer et al. [26] adopted an exact solution approach. In a location-transportation problem, the authors first determined the location and number of distribution centers and then designed a path to distribute the aid from distribution centers to the demand points. Barzinpour and Esmaeili [27] suggested a location-allocation model with the objective of maximizing cumulative coverage of population and of minimizing total cost using goal programming.Some authors deal with the location-allocation problem for medical services required during a disaster. There can be different objective functions for medical services, such as deciding the number of dispensing sites required, maximal covering objectives, minimum number of staff required, or maximizing demand satisfaction. Murali et al. [28], Jia et al. [29], and Lee et al. [30] considered location-allocation problem in terms of medical services. Details of all the models for facility location in HSC are given in Table 1.Table 1
Facility location models in HSC.
Authors
Objective function
Constraints/decision
Problem type
Balcik and Beamon [21]
Maximize (demand coverage by distribution centers)
Budget constraint,  inventory level at  distribution centers
Maximal covering location model
Bozorgi-Amiri et al. [22]
Minimize (costs for predisaster setup, procurement, transportation, holding, shortage)
Capacity for relief distribution canter, commodity flow, supply and demand
Location-allocation model
Horner and Downs [23]
Minimize (costs of distributing relief goods)
Demand fulfilment constraint, number of distribution centers
Intermediate distribution facility model
Dekle et al. [15]
Minimize (facilities for each area with a given distance)
Identify the location of the facility for each area
Covering location model
Hong et al. [25]
Minimize (total logistics cost)
Distance between warehouse and facility, number of facilities, demand
Facility location model
Chang et al. [16]
Minimize (transportation cost, facility setup cost, distance of rescue equipment cost)
Number of facilities and their capacity, prioritization of facility allocation, storage, shortage, penalties for surplus
Location allocation model
McCall [17]
Minimize (victim nautical miles, shortage)
Facility capacity, number of kits for prepositioning before disaster, unsatisfied demands
Facility location model
Rawls and Turnquist [91]
Minimize (costs of facility opening, unsatisfied demand, transportation)
Location and inventory level decision at each facility
Location-allocation model
Zhang et al. [24]
Minimize (cost of the total time of dispatching emergency resources)
Equilibrium of supply and demand for primary disaster, equilibrium of supply and demand for potential secondary disaster, resources available for secondary disaster
Location-allocation model
Akgün et al. [18]
Minimize (risk for unsatisfied demand)
Response time, distance between facility and disaster point
Facility location model
Barzinpour and Esmaeili [27]
Maximize (cumulative coverage of population)Minimize (total cost)
Demand and supply, transportation capacity, facility storage capacity
Location-allocation model
Abounacer et al. [26]
Minimize (distance from distribution center to demand point, number of facilities, unsatisfied demand)
Daily working hours, supply and demand, vehicle capacity
Location-transportation model
Rawls and Turnquist [19]
Minimize  (costs of commodity acquisition, stocking decision, transportation, shortage, holding)
Demand, number of facilities, inventory level
Dynamic allocation model
Murali et al. [28]
Maximize (number of people who receive medication)
Supply and demand, distance, facility capacity
Maximal covering location model
Lin et al. [20]
Minimize (shortage penalty cost, delayed delivery penalty cost, shipping cost, unfairness of service cost)
Number of depots, vehicles, travel time, delivery items quantity
Facility location model


### 3.2. Network Design and Relief Distribution in Humanitarian Supply Chain

Various mathematical models with little variation in constraints and objective functions have been developed in the research area of network design and relief distribution for the HSC. Some of the constraints are explained here. Transportation capacity constraint defines the capacity of the vehicle being used in operation. The number of vehicles defines how many vehicles are available for relief operation. Another important constraint is fleet composition which defines whether the company has homogenous fleet or heterogeneous fleet. The number of depots constraints defines whether the problem considers a single depot or multiple depots, because network design of single depot problem will be totally different from multiple depot problem. In some models a penalty function is assumed for any unsatisfied demand. There are various types of objective functions considered in network design and relief distribution but the general theme of objective functions is to increase responsiveness and cost efficiency. In the relief goods distribution process considering only cost minimization is not an appropriate objective function, because in extreme conditions when demand increases exponentially and shortages occur, the primary goal is to save human lives.HSC network design and relief distribution models are very complex in nature. As we continue adding details, models become more complex. Although adding details makes the model more practical, the difficulty level to solve the model increases rapidly. The HSC network design and relief distribution models mostly consider multiproduct and multiperiod with a deterministic approach. However, stochastic approach seems to be more practical due to high level of uncertainty associated with HSC operations. The problem with the consideration of multiproduct, multiperiod, and multiobjective with stochastic approach is that it becomes hard to obtain global optimal solution. Very few authors have suggested models that are multiperiod, multiobjective, and stochastic in nature.In the last decade work in the field of humanitarian network design and relief distribution has been done in various spectrums. Some researchers considered routing problem along with location of the distribution center in their problem called location-routing model. Ahmadi et al. [31] suggested a multidepot location-routing model. They determined the location of local depots and developed a routing model in order to deliver the relief goods from local depots to the affected people. They added a factor of network failure considerations, like road destruction and penalty cost for unsatisfied demands, and developed a neighborhood search algorithm. Campbell et al. [32] developed a model in which they considered two objective functions, the first one for a travelling salesman problem to minimize the maximum arrival time and the second one for the vehicle routing problem to minimize the average arrival time. Jabbarzadeh et al. [33] presented a robust network design model for blood supplies with an objective function of minimizing cost for holding, facility location, and transportation. Their model can be used for blood facility locations and allocation of blood for the injured people in a disaster stricken area. Wohlgemuth et al. [34] addressed pickup and delivery problem in which locations can perform both functions of receiving and sending goods. They considered the routing and scheduling for forwarding agencies in disaster relief operation with the objective of avoiding delays in delivery and maximizing equipment utilization. Shen et al. [35] analyzed the problem of vehicle routing from the perspective of a bioterrorism emergency. They proposed 2-stage model for disaster. In the first stage, which is called planning stage, routes are proposed. In the second stage, when some information is available, plans are modified according to the updated information. To formulate the problem they used mixed-integer linear programming and proposed an approximation heuristic. Bozorgi-Amiri et al. [36] suggested a supplies prepositioning and aid distribution model using robust stochastic programming approach. They considered uncertainty factor in supplies and cost of procurement with the objective function of minimizing the total cost of relief chain and maximizing the satisfaction level.Some researchers have developed models solely for relief distribution operation. For instance, Tzeng et al. [37] constructed a relief distribution model utilizing a fuzzy multiobjective programming method. They pursued the efficiency by making an objective of minimizing cost and travel time and fairness by using the objective of maximizing minimal satisfaction level. Sheu [38] suggested a postdisaster demand oriented emergency logistics operational model. In the model for centralized distribution, a new approach is adopted that considers the factors of postdisaster psychology, resilience of survivor and emotional cognition, which affects the attitude of survivor. They considered the objective function of maximizing the collective resilience of survivors during emergency logistics operation. Afshar and Haghani [39] recommended an integrated mathematical model that controls the flow of different type of relief goods in large amounts during disaster. As relief goods should reach the disaster stricken areas as soon as possible to save the maximum number of people, objective function of maximizing the survival rate and minimizing the cost is considered. Lin et al. [40] proposed a logistics model in which they prioritized deliveries during disaster and made the response time minimum. They introduced the factor of a soft time window with multiperiod routing and explained the results for not prioritizing delivery during disaster relief operation. Vitoriano et al. [41] suggested a model for logistics operation with the objectives of minimizing cost, minimizing the maximum ransack probability of the road, and maximizing the minimum reliability of a link using goal programming. Chen et al. [42] developed an equipment distribution model just after the disaster, which was based on geographic information system (GIS) with the objective of minimizing the relief and resource allocation time. Liberatore et al. [43] took a new angle for distribution of emergency goods. They suggested that after disaster infrastructures like roads and bridges are unsafe. They made a model for the recovery of damaged elements of distribution network with the objective function of maximizing the demand satisfaction.Whenever donations are given to a disaster stricken country they are received at seaports or airports and transferred to large warehouses. From large warehouses they are sorted out and transported to the central warehouses. At central warehouses, demands from regional distribution centers are collected and relief goods are dispatched accordingly. From regional distribution centers, relief goods are delivered to the affected people. As conditions in a disaster stricken area are most of the time unfavorable, the management of distribution from regional distribution center to the affected people plays a critical role in demand satisfaction. Balcik et al. [44] research is among one of the pioneers in the field of last mile distribution. They studied making decisions for relief supply allocation to the distribution centers and determining the delivery schedule with the objective of minimizing the transportation cost and penalty cost for unsatisfied demand. Wohlgemuth et al. [34] worked on the planning and logistics operation for the last mile distribution in an emergency relief chain with the objective function of minimizing the total travel time and number of vehicles. Ahmadi et al. [31] constructed a model for location of depot and routing for last mile distribution in an earthquake stricken area with the objective of minimizing the distribution cost and cost for opening a local depot.After occurrence of a disaster, transportation of injured people and casualties to the hospitals is another important task of the HSC. In this area, Dean and Nair [45] developed a model with the purpose of effectively evacuating victims and transporting them to different hospitals, not overwhelming any single hospital with the objective function of maximizing the number of expected survivals from incident. Wang et al. [46] built an agent based simulation model for an urban area and simulated the emergency response for a mass casualty disaster using geographic information system (GIS). Wilson et al. [47] developed a combinatorial model for casualty processing with the objective function of minimizing fatalities and sufferings and maximizing efficiency. They employed a variable neighborhood descent metaheuristic to solve the model. Salman and Gül [48] provided a multiperiod mixed-integer programming model that optimizes capacity allocation and casualty transportation decisions with the objectives of minimizing travel and waiting time and cost of establishing a new facility. Apte et al. [49] presented a study for Columbia for which they developed a tool to assist planners in selecting location of casualty collection point in a disaster stricken area with the objective function of maximizing weighted throughput of casualties and minimizing travelling time to shelter for casualties.Some other research works in the HSC include Hu and Sheu [50] who proposed a novel reverse logistics system in the postdisaster scenario. In their model, they included the reverse logistics cost and psychological cost with the objective of minimizing logistical cost, environmental and operational cost, and psychological cost. Similarly, Chiou and Lai [51] presented an integrated model that contains rescue path model, postdisaster traffic assignment model, and traffic controlled arc selection model. In their model objectives of minimizing travel time for rescue path, total detour time and number of police men are considered. When a disaster occurs, aid from donor countries and international organizations is sent to the disaster stricken country. For the effective management of these relief goods Adıvar and Mert [52] and Camacho-Vallejo et al. [53] suggested models with the objectives of minimizing total response time for delivering and cost for procurement. Details of all the models discussed above for network design and disaster relief distribution in HSC are given in Table 2.Table 2
Relief distribution models in HSC.
Authors
Objective function
Constraints/decision
Problem type
Ahmadi et al. [31]
Minimize (distribution time, penalty cost of unsatisfied demand, fixed costs of opening local depot)
Arrival and destination, number of vehicles, demand, working time, depot capacity
Multidepot location-routing model
Yi and Kumar [78]
Minimize (weighted sum of unsatisfied demand)
Flow of wounded people, number of unserved wounded people, vehicle load and capacity, number of vehicles
Multicommodity network flow model
Vitoriano et al. [96]
Minimize (time, cost)Maximize (equity, reliability)
Supply and demand balance at each node, vehicle type, subcycle elimination, vehicle capacity
Relief distribution model
Tzeng et al. [37]
Minimize (transportation cost, travel distance)Maximize (minimum satisfaction)
Shipment period, selection of depot, uncertain demand
Relief distribution model
Chen et al. [42]
Minimize (decision making and equipment transportation time)
Balance of inflow and outflow at each node, vehicle routing
Relief equipment distribution model
Wang et al. [72]
Minimize (travelling time, relief  distribution cost)Maximize (route reliability)
Vehicle arrival and destination, quantity of relief, demand and supply of relief, vehicle capacity
Multiobjective open location-routing model
Jabbarzadeh et al. [33]
Minimize (costs of locating blood facilities, transportation, and holding)
Location and number of facilities,  quantity of blood required at each facility, blood inventory level at the end of each period
Robust network design model
Balcik et al. [44]
Minimize (logistic costs, penalty cost, and shortage cost)
Demand fulfilment, vehicle capacity
Last mile relief  distribution model
Tirado et al. [98]
Minimize (deviation of delivered aid with respect to the planned amount)
Dynamic flow balance at each node, flow balance for vehicle, vehicle availability, vehicle capacity, amount of load
Lexicographical dynamic flow model
Liberatore et al. [43]
Maximize (demand satisfaction)
Arrival time, total served demand, maximum ransack probability, arc reliability
Humanitarian aid distribution model
Campbell et al. [32]
Minimize (maximum and minimum average arrival time)
Subtour elimination, vehicle route destination, arrival time
Travelling salesman problem (TSP) and vehicle routing problem (VRP)
Sheu [38]
Maximize (collective resilience of survivors during emergency logistics operations)
Population size, number of affected areas, setup cost, transportation cost, relief demand and supply
Relief distribution model
Afshar and Haghani [39]
Minimize (total amount of weighted unsatisfied demand)
Commodity flow, vehicular flow, facility location, capacities for temporary facilities
Relief distribution model
Huang et al. [79]
Minimize (sum of arrival times to beneficiaries)
Number of vehicles, flow balance, subtour elimination, arrival time
Assessment routing model
Bozorgi-Amiri et al. [36]
Minimize (total cost of the relief chain, sum of the maximum shortages)Maximize (satisfaction level)
Commodity flow, capacity limits of distribution centers,  number of distribution centers
Relief distribution model
Özdamar and Demir [64]
Minimize (estimated total travel time)
Commodity flow balance, unmet demands, inventory level at warehouse, vehicle capacity, number of vehicles, number of routes
Vehicle routing model
Hu and Sheu [50]
Minimize (logistical costs, environmental and operational risk costs, and psychological costs)
Recycled amounts for use, stocks of the debris amounts stocked, recycled, transported, and disposed, debris transportation
Postdisaster debris reverse logistics model
Lin et al. [40]
Minimize (penalty function, unsatisfied demands, and total travel time)
Maximum service level, fairness, vehicle capacity, working hours
Relief distribution model
Wohlgemuth et al. [34]
Minimize (delays in delivery time) Maximize (equipment utilization)
Vehicle capacity, subtour elimination, time window, time consistency
Last mile relief  distribution model
Vitoriano et al. [41]
Minimize (operation cost, maximum ransack probability)Maximize (reliability in a link)
Availability of goods, vehicles flow, vehicle load,  budget
Humanitarian aid distribution system
Shen et al. [35]
Minimize (unsatisfied demand)
Route feasibility, time, service, demand flow
Vehicle routing model
Chiou and Lai [51]
Minimize (travel time of rescue path, total detour travel time, number of unconnected trips of nonvictims, and number of police officers)
Access reliability, traffic capacity, degree of damage of transportation facility
Optimal rescue path and traffic control model
Berkoune et al. [71]
Minimize (total duration of all trips)
Relief goods availability, supply and demand, daily work time, vehicle type, vehicle capacity
Relief distribution model
Adıvar and Mert [52]
Maximize (minimum credibility with respect to every item)Minimizing (total cost of procurement plus transportation)
Flow conservation, capacity limitations of the available transportation assets, time period in which the relief item arrives, available number of vehicles
International relief planning model
Camacho-Vallejo et al. [53]
Minimize (total response time for delivering aid, cost of transportation)
Available space in each storage center, relief goods quantity,  demand and supply
International aid distribution model


### 3.3. Mass Evacuation

The type of evacuation method in the time of disaster depends on the factors of the location and size of the area to be evacuated. Population density and infrastructure of a location can vary from place to place; because in urban areas the population density is high, it needs different type of planning as compared to a place where population density is low. Meanwhile, evacuation teams have to consider people with special needs during evacuation planning. Mass evacuation models can be divided mainly into the following three categories:(i)
Public transport evacuation model and private transport evacuation model: In the public transport evacuation model most of the time pickup stops, shelters, and bus routes to reach depot are optimized, while in the private transport evacuation model the main issue is to make a smooth flow of traffic.
(ii)
Urban area evacuation model: Urban areas have high population density, so in urban evacuation model flow conservation, inflow capacity, street capacity, and lane consistencies are optimized.
(iii)
No-notice evacuation and short-notice evacuation: No-notice and short-notice evacuations have a difference in terms of the predictability of occurrence for the disaster. A no-notice evacuation needs to start immediately after the warning, while a short-notice evacuation may have a lead time of 24–72 hours before the disasters strike [54]. This lead time is the essential difference between no-notice and short-notice evacuations. In these evacuation category objective functions of maximizing number of evacuees and flow rate, minimizing the cost for transportation and time for evacuation is considered.Most of the models in this mass evacuation area are developed for public transport evacuation. Sheu and Pan [55] considered public evacuation as a part of their integrated emergency supply network model. In their model they considered the objective of minimizing travel distance, operational cost, and psychological cost. They integrated shelter networks, medical networks, and distribution networks to make a centralized emergency supply network. Whenever a disaster like flood is going to occur, governments evacuate the people from risky areas ahead of time. In such situations the most common problem is the availability of a sufficient number of bus drivers. Viewing this problem Morgul et al. [56] suggested two stochastic models to determine extra drivers needed during an emergency evacuation operation with the objective of minimizing cost for the unsatisfied demand and costs for the hired extra board staff. Naghawi and Wolshon [57] also made a bus-based evacuation model in which they evaluated the impact of a transit bus-based evacuation on the operation of a regional road network.With regard to private evacuation research, the study of Chiu and Zheng [58] used a cell transmission model in which certain groups, for example, doctors, have priority during the evacuation process with the objective of minimizing total prioritized travel time of all mobilization superiority groups. Hsu and Peeta [59] proposed an information based model in order to make private evacuation operations successful. By using behavior-consistent information strategies, they addressed the demand and supply interactions with the objective of minimizing absolute difference between the desired proportions and the predicted proportions of evacuees taking evacuation routes. In another paper, Hsu and Peeta [60] suggested a stage- based evacuation operation on the basis of the evacuation risk zone. They made the evacuation risk zone on the basis of disaster characteristics, traffic demand pattern, and network supply conditions.There are some models which are made with the purpose to deliver relief to the affected persons and to evacuate those people who are still stuck in the disaster. Among them, Özdamar and Yi [61] suggested a model in which vehicles exploit any foreseeable opportunity by greedy neighborhood search method for disaster relief or evacuation operation purpose. Kongsomsaksakul et al. [62] suggested a bilevel location-allocation model for route selection and demand allocation for the shelters. In the upper level problem is a location problem in which authority’s decision is modelled, while in the lower level evacuees’ decision is modelled using a game theory concept and this problem is solved by genetic algorithm. Another such model is developed by Najafi et al. [63] with the objective of minimizing the total unserved injured people, unsatisfied demands, and vehicles utilized during the distribution of relief and evacuation operation in an earthquake hit area. Özdamar and Demir [64] suggested a work with the purpose of relief distribution and evacuation activity. They proposed a clustering algorithm called a hierarchical cluster and route procedure in order to coordinate vehicle routing. In another paper, Yi and Özdamar [65] suggested a mixed-integer multicommodity network flow model with the objective function of minimizing the weighted sum of unsatisfied demand and weighted sum of wounded people waiting for evacuation.Evacuation for disaster in an urban area is much more difficult than in a nonurban area, because high density of population and limited flow capacity of streets restrict evacuation process. The situation becomes worse, if it is a short-notice evacuation. Chen et al. [66] for evacuation in an urban area proposed that traffic signals can facilitate the evacuation process. In a simulation model, they used the signal timing for smoothening the traffic flow and facilitating evacuation process for a short-notice disaster in an urban area. Sayyady and Eksioglu [67] developed a mixed-integer linear programming model with the objective of minimizing evacuation time and number of causalities for an urban area. Bretschneider and Kimms [68] made a model that reorganizes a traffic routing in particular area for evacuation purpose with the objective of minimizing weighted sum of flow for evacuation process. Yuan and Puchalsky [69] presented a simulation model for evacuation in the city of Philadelphia using a dynamic sequential assignment method and simulated the interaction of private and public transport. Kirby et al. [70] suggested the concept of a regional hub reception center (RHRC) for evacuation in urban areas in case of a nuclear blast disaster. Basically, RHRC is a temporary shelter which provides basic needs to evacuees and registers them for allocation to permanent shelters.Details of all the models discussed above for mass evacuation process in HSC are given in Table3.Table 3
Mass evacuation models in HSC.
Authors
Objective function
Constraints/decision
Problem type
Sheu and Pan [55]
Minimize (travel distance, operational cost, psychological cost)
Number of evacuees, number of affected people, evacuee flow in shelter, distance travelled by  the evacuee
Centralized supply network model
Bish [97]
Minimize (total duration of the evacuation)
Flow balance at demand nodes, vehicle capacity, shelter capacity, evacuees delivery to shelter, multidepot, multitrip
Bus-based evacuation planning model
Bretschneider and Kimms [68]
Minimize (weighted sum of flow during evacuation)
Flow conservation, inflow capacity, street capacity, lane consistency
Urban areas evacuation model
Sayyady and Eksioglu [67]
Minimize (total evacuation time, number of casualties)
Vehicle capacity, flow conservation for citizens, number of vehicles, citizens pickup location
Urban areas evacuation model
Chiu and Zheng [58]
Minimize (total prioritized travel time of all mobilization priority groups)
Flow conservation at the source cells, maximum service flow rate, shelter capacity
Cell transmission model for evacuation
Naghawi and Wolshon [57]
Maximize (total number of evacuees)
Evacuation routes, flows
Bus-based evacuation model
Hsu and Peeta [60]
Maximize (total number of evacuees)
Number of vehicles, resource availability
Stage-based evacuation model
Özdamar and Yi [61]
Minimize (total service delay)
Demand and supply balance, flow of evacuees, number of injured people served, number of vehicles, vehicle capacity
Disaster relief distribution and evacuation model
Yi and Özdamar [65]
Minimize (weighted sum of unsatisfied demand, weighted sum of wounded people waiting)
Vehicle capacity, type of vehicle, number of vehicles, number of wounded people waiting, number of injured people served
Location-distribution model
Najafi et al. [63]
Minimize (total unserved injured people, unsatisfied demands, vehicles utilized)
Transportation of commodities, vehicle type, hospital capacity
Postdisaster distribution and evacuation model
Morgul et al. [56]
Minimize (cost for the unsatisfied demand, costs for the hired extra board staff)
Quality of service
Bus-based evacuation planning model
Kongsomsaksakul et al. [62]
Minimize (travel time)
Flow conservation, production constraint, flow and demand
Distribution and assignment model
Hsu and Peeta [59]
Minimize (absolute difference between the desired proportions and the predicted proportions of evacuees taking evacuation routes)
Linguistic message, demand conservation, evacuee route choice
Behavior-consistent information based network traffic control model
Özdamar and Demir [64]
Minimize (total travel time, efficient vehicle utilization)
Material flow balance, unsatisfied demand, inventory level at warehouse, vehicle capacity
Postdisaster distribution and evacuation model


## 3.1. Facility Location in Humanitarian Supply Chain

After the occurrence of a disaster when assessment of damage is complete, next step is to provide relief to the affected people. For this purpose, a relief distribution channel is designed which includes central distribution centers, warehouses, and regional relief distribution centers. The challenge in this process is to locate these facilities in such a way that demand of entire disaster stricken area is fulfilled with minimal delivering cost and maximum service level. Factors that increase the complexity of facility location task are as follows:(i)
Very short lead time and sudden surge in demand just after the disaster.
(ii)
Uncertainty in timing and scale of the disaster.
(iii)
Very high stakes associated with the timely delivery of relief goods.
(iv)
Damaged infrastructure.Research in the section of facility location for a disaster can be classified into two categories. The classification is with respect to the timeline of the disaster, whether the model is developed for a predisaster scenario or postdisaster scenario.In the predisaster scenario, a strategic approach is adopted in which researchers plan for the best locations of the facilities in order to distribute relief goods at minimal cost and maximum service level. Some authors along with the facility location planning also consider allocation of resources from central distribution centers to warehouses and regional distribution centers. Such problems are called location-allocation problems. For instance, Dekle et al. [15] constructed a model to locate the disaster recovery centers in the predisaster context and maximized the covering of relief goods in the disaster target area. Chang et al. [16] developed a model for resource distribution in urban flood disaster scenario. The authors considered the location of the warehouse, prioritization of facility allocation, shortage and penalties for surplus with the objective function of minimizing transportation cost, facility setup cost, and transportation cost of rescue equipment. McCall [17] presented a model with the prepositioning of assistance pack-up kits during disaster while considering the constraints of model facility capacity, number of kits for prepositioning before disaster, and unsatisfied demands with the objective function of minimizing victim nautical miles and shortages. Akgün et al. [18] adopted a different approach for facility location model and suggested that prepositioning supplies in the preparation phase is desirable to locate supplies near the expected disaster area so that aid can be distributed within a minimal time. On the other hand it is also important that if distribution centers are close to expected disaster area, then these distribution centers may get destroyed. The goal of the study is to choose locations so that reliability can be maximized. Just after the occurrence of a disaster, there is an urgent need of basic emergency supplies and this surge in demand continues for the first 72 hours. Keeping this point in view, Rawls and Turnquist [19] and Lin et al. [20] developed a model for the short term disaster response.On the other side of the disaster timeline there are some researchers who consider the postdisaster scenario for the facility location problem. In this postdisaster scenario, it is assumed that, after the occurrence of disaster, the assessment of the damage has been completed and all the major information required for a relief operation is available. For instance, Balcik and Beamon [21] suggested a facility location problem in the postdisaster scenario. The problem is a variant of maximal covering location model with budget and capacity constraints. They determined a number of locations of distribution centers and the amount of relief goods at each center to meet the demand. Bozorgi-Amiri et al. [22] presented a stochastic model that first determines the location of aid distribution centers and then allocates relief goods to the affected area with an objective function of minimizing the costs for predisaster setup that includes costs of procurement, transportation, holding, and shortage. Horner and Downs [23] proposed a model that is a variant of warehouse location problem in which they used warehouses to allocate the relief goods to the affected areas. In fact, it would be better to consider them as an intermediate distribution facility model with the objective function of minimizing the costs of distributing of relief goods. Zhang et al. [24] developed a resource allocation model while considering the constraints of multiresources and multidepots, with the objective function of minimizing the cost of the total time for dispatching emergency resources. Hong et al. [25] made a comparison of two robust mathematical models. The authors compared the robust integer facility location with the robust continuous facility location model and concluded that robust integer facility location works best in normal conditions, while the robust continuous facility location works best in the hours of disaster. Abounacer et al. [26] adopted an exact solution approach. In a location-transportation problem, the authors first determined the location and number of distribution centers and then designed a path to distribute the aid from distribution centers to the demand points. Barzinpour and Esmaeili [27] suggested a location-allocation model with the objective of maximizing cumulative coverage of population and of minimizing total cost using goal programming.Some authors deal with the location-allocation problem for medical services required during a disaster. There can be different objective functions for medical services, such as deciding the number of dispensing sites required, maximal covering objectives, minimum number of staff required, or maximizing demand satisfaction. Murali et al. [28], Jia et al. [29], and Lee et al. [30] considered location-allocation problem in terms of medical services. Details of all the models for facility location in HSC are given in Table 1.Table 1
Facility location models in HSC.
Authors
Objective function
Constraints/decision
Problem type
Balcik and Beamon [21]
Maximize (demand coverage by distribution centers)
Budget constraint,  inventory level at  distribution centers
Maximal covering location model
Bozorgi-Amiri et al. [22]
Minimize (costs for predisaster setup, procurement, transportation, holding, shortage)
Capacity for relief distribution canter, commodity flow, supply and demand
Location-allocation model
Horner and Downs [23]
Minimize (costs of distributing relief goods)
Demand fulfilment constraint, number of distribution centers
Intermediate distribution facility model
Dekle et al. [15]
Minimize (facilities for each area with a given distance)
Identify the location of the facility for each area
Covering location model
Hong et al. [25]
Minimize (total logistics cost)
Distance between warehouse and facility, number of facilities, demand
Facility location model
Chang et al. [16]
Minimize (transportation cost, facility setup cost, distance of rescue equipment cost)
Number of facilities and their capacity, prioritization of facility allocation, storage, shortage, penalties for surplus
Location allocation model
McCall [17]
Minimize (victim nautical miles, shortage)
Facility capacity, number of kits for prepositioning before disaster, unsatisfied demands
Facility location model
Rawls and Turnquist [91]
Minimize (costs of facility opening, unsatisfied demand, transportation)
Location and inventory level decision at each facility
Location-allocation model
Zhang et al. [24]
Minimize (cost of the total time of dispatching emergency resources)
Equilibrium of supply and demand for primary disaster, equilibrium of supply and demand for potential secondary disaster, resources available for secondary disaster
Location-allocation model
Akgün et al. [18]
Minimize (risk for unsatisfied demand)
Response time, distance between facility and disaster point
Facility location model
Barzinpour and Esmaeili [27]
Maximize (cumulative coverage of population)Minimize (total cost)
Demand and supply, transportation capacity, facility storage capacity
Location-allocation model
Abounacer et al. [26]
Minimize (distance from distribution center to demand point, number of facilities, unsatisfied demand)
Daily working hours, supply and demand, vehicle capacity
Location-transportation model
Rawls and Turnquist [19]
Minimize  (costs of commodity acquisition, stocking decision, transportation, shortage, holding)
Demand, number of facilities, inventory level
Dynamic allocation model
Murali et al. [28]
Maximize (number of people who receive medication)
Supply and demand, distance, facility capacity
Maximal covering location model
Lin et al. [20]
Minimize (shortage penalty cost, delayed delivery penalty cost, shipping cost, unfairness of service cost)
Number of depots, vehicles, travel time, delivery items quantity
Facility location model


## 3.2. Network Design and Relief Distribution in Humanitarian Supply Chain

Various mathematical models with little variation in constraints and objective functions have been developed in the research area of network design and relief distribution for the HSC. Some of the constraints are explained here. Transportation capacity constraint defines the capacity of the vehicle being used in operation. The number of vehicles defines how many vehicles are available for relief operation. Another important constraint is fleet composition which defines whether the company has homogenous fleet or heterogeneous fleet. The number of depots constraints defines whether the problem considers a single depot or multiple depots, because network design of single depot problem will be totally different from multiple depot problem. In some models a penalty function is assumed for any unsatisfied demand. There are various types of objective functions considered in network design and relief distribution but the general theme of objective functions is to increase responsiveness and cost efficiency. In the relief goods distribution process considering only cost minimization is not an appropriate objective function, because in extreme conditions when demand increases exponentially and shortages occur, the primary goal is to save human lives.HSC network design and relief distribution models are very complex in nature. As we continue adding details, models become more complex. Although adding details makes the model more practical, the difficulty level to solve the model increases rapidly. The HSC network design and relief distribution models mostly consider multiproduct and multiperiod with a deterministic approach. However, stochastic approach seems to be more practical due to high level of uncertainty associated with HSC operations. The problem with the consideration of multiproduct, multiperiod, and multiobjective with stochastic approach is that it becomes hard to obtain global optimal solution. Very few authors have suggested models that are multiperiod, multiobjective, and stochastic in nature.In the last decade work in the field of humanitarian network design and relief distribution has been done in various spectrums. Some researchers considered routing problem along with location of the distribution center in their problem called location-routing model. Ahmadi et al. [31] suggested a multidepot location-routing model. They determined the location of local depots and developed a routing model in order to deliver the relief goods from local depots to the affected people. They added a factor of network failure considerations, like road destruction and penalty cost for unsatisfied demands, and developed a neighborhood search algorithm. Campbell et al. [32] developed a model in which they considered two objective functions, the first one for a travelling salesman problem to minimize the maximum arrival time and the second one for the vehicle routing problem to minimize the average arrival time. Jabbarzadeh et al. [33] presented a robust network design model for blood supplies with an objective function of minimizing cost for holding, facility location, and transportation. Their model can be used for blood facility locations and allocation of blood for the injured people in a disaster stricken area. Wohlgemuth et al. [34] addressed pickup and delivery problem in which locations can perform both functions of receiving and sending goods. They considered the routing and scheduling for forwarding agencies in disaster relief operation with the objective of avoiding delays in delivery and maximizing equipment utilization. Shen et al. [35] analyzed the problem of vehicle routing from the perspective of a bioterrorism emergency. They proposed 2-stage model for disaster. In the first stage, which is called planning stage, routes are proposed. In the second stage, when some information is available, plans are modified according to the updated information. To formulate the problem they used mixed-integer linear programming and proposed an approximation heuristic. Bozorgi-Amiri et al. [36] suggested a supplies prepositioning and aid distribution model using robust stochastic programming approach. They considered uncertainty factor in supplies and cost of procurement with the objective function of minimizing the total cost of relief chain and maximizing the satisfaction level.Some researchers have developed models solely for relief distribution operation. For instance, Tzeng et al. [37] constructed a relief distribution model utilizing a fuzzy multiobjective programming method. They pursued the efficiency by making an objective of minimizing cost and travel time and fairness by using the objective of maximizing minimal satisfaction level. Sheu [38] suggested a postdisaster demand oriented emergency logistics operational model. In the model for centralized distribution, a new approach is adopted that considers the factors of postdisaster psychology, resilience of survivor and emotional cognition, which affects the attitude of survivor. They considered the objective function of maximizing the collective resilience of survivors during emergency logistics operation. Afshar and Haghani [39] recommended an integrated mathematical model that controls the flow of different type of relief goods in large amounts during disaster. As relief goods should reach the disaster stricken areas as soon as possible to save the maximum number of people, objective function of maximizing the survival rate and minimizing the cost is considered. Lin et al. [40] proposed a logistics model in which they prioritized deliveries during disaster and made the response time minimum. They introduced the factor of a soft time window with multiperiod routing and explained the results for not prioritizing delivery during disaster relief operation. Vitoriano et al. [41] suggested a model for logistics operation with the objectives of minimizing cost, minimizing the maximum ransack probability of the road, and maximizing the minimum reliability of a link using goal programming. Chen et al. [42] developed an equipment distribution model just after the disaster, which was based on geographic information system (GIS) with the objective of minimizing the relief and resource allocation time. Liberatore et al. [43] took a new angle for distribution of emergency goods. They suggested that after disaster infrastructures like roads and bridges are unsafe. They made a model for the recovery of damaged elements of distribution network with the objective function of maximizing the demand satisfaction.Whenever donations are given to a disaster stricken country they are received at seaports or airports and transferred to large warehouses. From large warehouses they are sorted out and transported to the central warehouses. At central warehouses, demands from regional distribution centers are collected and relief goods are dispatched accordingly. From regional distribution centers, relief goods are delivered to the affected people. As conditions in a disaster stricken area are most of the time unfavorable, the management of distribution from regional distribution center to the affected people plays a critical role in demand satisfaction. Balcik et al. [44] research is among one of the pioneers in the field of last mile distribution. They studied making decisions for relief supply allocation to the distribution centers and determining the delivery schedule with the objective of minimizing the transportation cost and penalty cost for unsatisfied demand. Wohlgemuth et al. [34] worked on the planning and logistics operation for the last mile distribution in an emergency relief chain with the objective function of minimizing the total travel time and number of vehicles. Ahmadi et al. [31] constructed a model for location of depot and routing for last mile distribution in an earthquake stricken area with the objective of minimizing the distribution cost and cost for opening a local depot.After occurrence of a disaster, transportation of injured people and casualties to the hospitals is another important task of the HSC. In this area, Dean and Nair [45] developed a model with the purpose of effectively evacuating victims and transporting them to different hospitals, not overwhelming any single hospital with the objective function of maximizing the number of expected survivals from incident. Wang et al. [46] built an agent based simulation model for an urban area and simulated the emergency response for a mass casualty disaster using geographic information system (GIS). Wilson et al. [47] developed a combinatorial model for casualty processing with the objective function of minimizing fatalities and sufferings and maximizing efficiency. They employed a variable neighborhood descent metaheuristic to solve the model. Salman and Gül [48] provided a multiperiod mixed-integer programming model that optimizes capacity allocation and casualty transportation decisions with the objectives of minimizing travel and waiting time and cost of establishing a new facility. Apte et al. [49] presented a study for Columbia for which they developed a tool to assist planners in selecting location of casualty collection point in a disaster stricken area with the objective function of maximizing weighted throughput of casualties and minimizing travelling time to shelter for casualties.Some other research works in the HSC include Hu and Sheu [50] who proposed a novel reverse logistics system in the postdisaster scenario. In their model, they included the reverse logistics cost and psychological cost with the objective of minimizing logistical cost, environmental and operational cost, and psychological cost. Similarly, Chiou and Lai [51] presented an integrated model that contains rescue path model, postdisaster traffic assignment model, and traffic controlled arc selection model. In their model objectives of minimizing travel time for rescue path, total detour time and number of police men are considered. When a disaster occurs, aid from donor countries and international organizations is sent to the disaster stricken country. For the effective management of these relief goods Adıvar and Mert [52] and Camacho-Vallejo et al. [53] suggested models with the objectives of minimizing total response time for delivering and cost for procurement. Details of all the models discussed above for network design and disaster relief distribution in HSC are given in Table 2.Table 2
Relief distribution models in HSC.
Authors
Objective function
Constraints/decision
Problem type
Ahmadi et al. [31]
Minimize (distribution time, penalty cost of unsatisfied demand, fixed costs of opening local depot)
Arrival and destination, number of vehicles, demand, working time, depot capacity
Multidepot location-routing model
Yi and Kumar [78]
Minimize (weighted sum of unsatisfied demand)
Flow of wounded people, number of unserved wounded people, vehicle load and capacity, number of vehicles
Multicommodity network flow model
Vitoriano et al. [96]
Minimize (time, cost)Maximize (equity, reliability)
Supply and demand balance at each node, vehicle type, subcycle elimination, vehicle capacity
Relief distribution model
Tzeng et al. [37]
Minimize (transportation cost, travel distance)Maximize (minimum satisfaction)
Shipment period, selection of depot, uncertain demand
Relief distribution model
Chen et al. [42]
Minimize (decision making and equipment transportation time)
Balance of inflow and outflow at each node, vehicle routing
Relief equipment distribution model
Wang et al. [72]
Minimize (travelling time, relief  distribution cost)Maximize (route reliability)
Vehicle arrival and destination, quantity of relief, demand and supply of relief, vehicle capacity
Multiobjective open location-routing model
Jabbarzadeh et al. [33]
Minimize (costs of locating blood facilities, transportation, and holding)
Location and number of facilities,  quantity of blood required at each facility, blood inventory level at the end of each period
Robust network design model
Balcik et al. [44]
Minimize (logistic costs, penalty cost, and shortage cost)
Demand fulfilment, vehicle capacity
Last mile relief  distribution model
Tirado et al. [98]
Minimize (deviation of delivered aid with respect to the planned amount)
Dynamic flow balance at each node, flow balance for vehicle, vehicle availability, vehicle capacity, amount of load
Lexicographical dynamic flow model
Liberatore et al. [43]
Maximize (demand satisfaction)
Arrival time, total served demand, maximum ransack probability, arc reliability
Humanitarian aid distribution model
Campbell et al. [32]
Minimize (maximum and minimum average arrival time)
Subtour elimination, vehicle route destination, arrival time
Travelling salesman problem (TSP) and vehicle routing problem (VRP)
Sheu [38]
Maximize (collective resilience of survivors during emergency logistics operations)
Population size, number of affected areas, setup cost, transportation cost, relief demand and supply
Relief distribution model
Afshar and Haghani [39]
Minimize (total amount of weighted unsatisfied demand)
Commodity flow, vehicular flow, facility location, capacities for temporary facilities
Relief distribution model
Huang et al. [79]
Minimize (sum of arrival times to beneficiaries)
Number of vehicles, flow balance, subtour elimination, arrival time
Assessment routing model
Bozorgi-Amiri et al. [36]
Minimize (total cost of the relief chain, sum of the maximum shortages)Maximize (satisfaction level)
Commodity flow, capacity limits of distribution centers,  number of distribution centers
Relief distribution model
Özdamar and Demir [64]
Minimize (estimated total travel time)
Commodity flow balance, unmet demands, inventory level at warehouse, vehicle capacity, number of vehicles, number of routes
Vehicle routing model
Hu and Sheu [50]
Minimize (logistical costs, environmental and operational risk costs, and psychological costs)
Recycled amounts for use, stocks of the debris amounts stocked, recycled, transported, and disposed, debris transportation
Postdisaster debris reverse logistics model
Lin et al. [40]
Minimize (penalty function, unsatisfied demands, and total travel time)
Maximum service level, fairness, vehicle capacity, working hours
Relief distribution model
Wohlgemuth et al. [34]
Minimize (delays in delivery time) Maximize (equipment utilization)
Vehicle capacity, subtour elimination, time window, time consistency
Last mile relief  distribution model
Vitoriano et al. [41]
Minimize (operation cost, maximum ransack probability)Maximize (reliability in a link)
Availability of goods, vehicles flow, vehicle load,  budget
Humanitarian aid distribution system
Shen et al. [35]
Minimize (unsatisfied demand)
Route feasibility, time, service, demand flow
Vehicle routing model
Chiou and Lai [51]
Minimize (travel time of rescue path, total detour travel time, number of unconnected trips of nonvictims, and number of police officers)
Access reliability, traffic capacity, degree of damage of transportation facility
Optimal rescue path and traffic control model
Berkoune et al. [71]
Minimize (total duration of all trips)
Relief goods availability, supply and demand, daily work time, vehicle type, vehicle capacity
Relief distribution model
Adıvar and Mert [52]
Maximize (minimum credibility with respect to every item)Minimizing (total cost of procurement plus transportation)
Flow conservation, capacity limitations of the available transportation assets, time period in which the relief item arrives, available number of vehicles
International relief planning model
Camacho-Vallejo et al. [53]
Minimize (total response time for delivering aid, cost of transportation)
Available space in each storage center, relief goods quantity,  demand and supply
International aid distribution model


## 3.3. Mass Evacuation

The type of evacuation method in the time of disaster depends on the factors of the location and size of the area to be evacuated. Population density and infrastructure of a location can vary from place to place; because in urban areas the population density is high, it needs different type of planning as compared to a place where population density is low. Meanwhile, evacuation teams have to consider people with special needs during evacuation planning. Mass evacuation models can be divided mainly into the following three categories:(i)
Public transport evacuation model and private transport evacuation model: In the public transport evacuation model most of the time pickup stops, shelters, and bus routes to reach depot are optimized, while in the private transport evacuation model the main issue is to make a smooth flow of traffic.
(ii)
Urban area evacuation model: Urban areas have high population density, so in urban evacuation model flow conservation, inflow capacity, street capacity, and lane consistencies are optimized.
(iii)
No-notice evacuation and short-notice evacuation: No-notice and short-notice evacuations have a difference in terms of the predictability of occurrence for the disaster. A no-notice evacuation needs to start immediately after the warning, while a short-notice evacuation may have a lead time of 24–72 hours before the disasters strike [54]. This lead time is the essential difference between no-notice and short-notice evacuations. In these evacuation category objective functions of maximizing number of evacuees and flow rate, minimizing the cost for transportation and time for evacuation is considered.Most of the models in this mass evacuation area are developed for public transport evacuation. Sheu and Pan [55] considered public evacuation as a part of their integrated emergency supply network model. In their model they considered the objective of minimizing travel distance, operational cost, and psychological cost. They integrated shelter networks, medical networks, and distribution networks to make a centralized emergency supply network. Whenever a disaster like flood is going to occur, governments evacuate the people from risky areas ahead of time. In such situations the most common problem is the availability of a sufficient number of bus drivers. Viewing this problem Morgul et al. [56] suggested two stochastic models to determine extra drivers needed during an emergency evacuation operation with the objective of minimizing cost for the unsatisfied demand and costs for the hired extra board staff. Naghawi and Wolshon [57] also made a bus-based evacuation model in which they evaluated the impact of a transit bus-based evacuation on the operation of a regional road network.With regard to private evacuation research, the study of Chiu and Zheng [58] used a cell transmission model in which certain groups, for example, doctors, have priority during the evacuation process with the objective of minimizing total prioritized travel time of all mobilization superiority groups. Hsu and Peeta [59] proposed an information based model in order to make private evacuation operations successful. By using behavior-consistent information strategies, they addressed the demand and supply interactions with the objective of minimizing absolute difference between the desired proportions and the predicted proportions of evacuees taking evacuation routes. In another paper, Hsu and Peeta [60] suggested a stage- based evacuation operation on the basis of the evacuation risk zone. They made the evacuation risk zone on the basis of disaster characteristics, traffic demand pattern, and network supply conditions.There are some models which are made with the purpose to deliver relief to the affected persons and to evacuate those people who are still stuck in the disaster. Among them, Özdamar and Yi [61] suggested a model in which vehicles exploit any foreseeable opportunity by greedy neighborhood search method for disaster relief or evacuation operation purpose. Kongsomsaksakul et al. [62] suggested a bilevel location-allocation model for route selection and demand allocation for the shelters. In the upper level problem is a location problem in which authority’s decision is modelled, while in the lower level evacuees’ decision is modelled using a game theory concept and this problem is solved by genetic algorithm. Another such model is developed by Najafi et al. [63] with the objective of minimizing the total unserved injured people, unsatisfied demands, and vehicles utilized during the distribution of relief and evacuation operation in an earthquake hit area. Özdamar and Demir [64] suggested a work with the purpose of relief distribution and evacuation activity. They proposed a clustering algorithm called a hierarchical cluster and route procedure in order to coordinate vehicle routing. In another paper, Yi and Özdamar [65] suggested a mixed-integer multicommodity network flow model with the objective function of minimizing the weighted sum of unsatisfied demand and weighted sum of wounded people waiting for evacuation.Evacuation for disaster in an urban area is much more difficult than in a nonurban area, because high density of population and limited flow capacity of streets restrict evacuation process. The situation becomes worse, if it is a short-notice evacuation. Chen et al. [66] for evacuation in an urban area proposed that traffic signals can facilitate the evacuation process. In a simulation model, they used the signal timing for smoothening the traffic flow and facilitating evacuation process for a short-notice disaster in an urban area. Sayyady and Eksioglu [67] developed a mixed-integer linear programming model with the objective of minimizing evacuation time and number of causalities for an urban area. Bretschneider and Kimms [68] made a model that reorganizes a traffic routing in particular area for evacuation purpose with the objective of minimizing weighted sum of flow for evacuation process. Yuan and Puchalsky [69] presented a simulation model for evacuation in the city of Philadelphia using a dynamic sequential assignment method and simulated the interaction of private and public transport. Kirby et al. [70] suggested the concept of a regional hub reception center (RHRC) for evacuation in urban areas in case of a nuclear blast disaster. Basically, RHRC is a temporary shelter which provides basic needs to evacuees and registers them for allocation to permanent shelters.Details of all the models discussed above for mass evacuation process in HSC are given in Table3.Table 3
Mass evacuation models in HSC.
Authors
Objective function
Constraints/decision
Problem type
Sheu and Pan [55]
Minimize (travel distance, operational cost, psychological cost)
Number of evacuees, number of affected people, evacuee flow in shelter, distance travelled by  the evacuee
Centralized supply network model
Bish [97]
Minimize (total duration of the evacuation)
Flow balance at demand nodes, vehicle capacity, shelter capacity, evacuees delivery to shelter, multidepot, multitrip
Bus-based evacuation planning model
Bretschneider and Kimms [68]
Minimize (weighted sum of flow during evacuation)
Flow conservation, inflow capacity, street capacity, lane consistency
Urban areas evacuation model
Sayyady and Eksioglu [67]
Minimize (total evacuation time, number of casualties)
Vehicle capacity, flow conservation for citizens, number of vehicles, citizens pickup location
Urban areas evacuation model
Chiu and Zheng [58]
Minimize (total prioritized travel time of all mobilization priority groups)
Flow conservation at the source cells, maximum service flow rate, shelter capacity
Cell transmission model for evacuation
Naghawi and Wolshon [57]
Maximize (total number of evacuees)
Evacuation routes, flows
Bus-based evacuation model
Hsu and Peeta [60]
Maximize (total number of evacuees)
Number of vehicles, resource availability
Stage-based evacuation model
Özdamar and Yi [61]
Minimize (total service delay)
Demand and supply balance, flow of evacuees, number of injured people served, number of vehicles, vehicle capacity
Disaster relief distribution and evacuation model
Yi and Özdamar [65]
Minimize (weighted sum of unsatisfied demand, weighted sum of wounded people waiting)
Vehicle capacity, type of vehicle, number of vehicles, number of wounded people waiting, number of injured people served
Location-distribution model
Najafi et al. [63]
Minimize (total unserved injured people, unsatisfied demands, vehicles utilized)
Transportation of commodities, vehicle type, hospital capacity
Postdisaster distribution and evacuation model
Morgul et al. [56]
Minimize (cost for the unsatisfied demand, costs for the hired extra board staff)
Quality of service
Bus-based evacuation planning model
Kongsomsaksakul et al. [62]
Minimize (travel time)
Flow conservation, production constraint, flow and demand
Distribution and assignment model
Hsu and Peeta [59]
Minimize (absolute difference between the desired proportions and the predicted proportions of evacuees taking evacuation routes)
Linguistic message, demand conservation, evacuee route choice
Behavior-consistent information based network traffic control model
Özdamar and Demir [64]
Minimize (total travel time, efficient vehicle utilization)
Material flow balance, unsatisfied demand, inventory level at warehouse, vehicle capacity
Postdisaster distribution and evacuation model


## 4. Discussion and Implications

In this section, the questions we have developed in the planning phase of this systematic literature review study are discussed with their implications.RQ 1. What Is the Current Status of Research on the Mathematical Models in the Field of HSC? After reviewing the selected articles, our first observation is that most of the focus of researchers in the HSC is on preparedness and the response phase, considering facility location problem, relief distribution problem, and evacuation problems. Research in the mitigation and recovery phase of the HSC, which includes debris management and repair of effected infrastructure, has a very small portion, particularly the long term recovery after disaster is the most neglected area. These areas need the attention of researchers.The second observation is that the proportion of deterministic models is more as compared to stochastic or fuzzy models. However, stochastic or fuzzy models are more practical due to an uncertainty factor involved in HSC operations. In the preparedness and response phase, demand fulfilment and location of demand are major factors that influence the efficiency of the operation. That is why most of the researchers considered both of these factors as an uncertain component in their research. In the response phase, roads are damaged, so relief distribution operation route reliability is a major concern. Viewing the importance of this factor, many researchers considered the route reliability as uncertain in their research. A few researchers also considered the meeting needs on time, demand, and supply as uncertain in the research for HSC.Another observation in this study is that some authors have used algorithms and heuristics in their research. Heuristics enable us to find quality solutions for many problems. Although their main disadvantage is that they provide near optimal solutions (not global optimal), they are very useful. In the recent times genetic algorithm, particle swarm optimization, and ant colony optimization are the different methods that have facilitated the disaster relief operation. Berkoune et al. [71], Kongsomsaksakul et al. [62], Jia et al. [29], Wang et al. [72], Yang et al. [73], and Hamedi et al. [74] used genetic algorithms in HSC research. Some other authors used the particle swarm optimization (PSO) in different scenarios of disaster. PSO is metaheuristic, which is a population based optimization method. It uses a number of candidate solutions which are called particles. PSO can search for many particles, but the disadvantage is that it does not guarantee the optimal solution. Bozorgi-Amiri et al. [22], Cheng et al. [75], and Gan et al. [76] utilized PSO for their research in HSC. Ant colony optimization (ACO) is another useful optimization technique which imitates the behavior of ants which live in colonies and communicate with each other for the shortest path in search of food. ACO has been successfully implemented in vehicle routing, travelling salesman problem, assignment problem, and scheduling in disasters. Yan and Shih [77] and Yi and Kumar [78] considered ACO approach. Other metaheuristic algorithms include Tabu search, variable neighborhood search, and simulated annealing. Tabu search is a metaheuristic search method. It takes a search for a problem in its immediate neighborhood to find improved solutions. Tabu search enhances the performance of the local search but its main drawback is that it sticks in suboptimal regions. Huang et al. [79], Wohlgemuth et al. [34], Shen et al. [35], and Sayyady and Eksioglu [67] used the Tabu search algorithm in their research. Simulated annealing (SA) is metaheuristic for global optimization problem. For some problems, SA can be a good option given that the ultimate goal is to find an acceptable solution in a given amount of time. Murali et al. [28] used a simulated annealing algorithm for facility location problem in HSC. Variable neighborhood search (VNS) is another metaheuristic method to get solution in a global optimization problem. VNS looks for distant neighborhood and moves for a new one if there is any improvement in solution. Ahmadi et al. [31] and Özdamar and Yi [61] used the VNS approach in their research.Some research papers used simulation methodology. The authors using simulation methodology mostly utilize GIS related applications and spatial data analysis technique. The main goal in a simulation model is to forecast the possible scenarios. For example, Chang et al. [16] and Horner and Downs [23] developed simulation models. Simulation models are also developed for the evacuation process. Among them, most of the simulation models are developed to smoothen the traffic flow during evacuation process. Naghawi and Wolshon [57], Yuan and Puchalsky [69], and Noreña et al. [80] developed simulation models in different evacuation scenarios.In order to know about the current status and trends of research in HSC we evaluated each paper on the basis of four points of disaster phase, uncertain components considered in the model, model formulation technique, and solution methodology used by the author. These four points will define the current research status of HSC. According to these points the research papers considered for review in the areas of facility location problem, network design and relief distribution problem, and mass evacuation problem are given in Tables4, 5, and 6, respectively.Table 4
Facility location models in HSC.
Author
Disaster phase
Uncertain component
Model formulation
Solution technique
Balcik and Beamon [21]
Preparation and response phase(maximal covering location model)
Demand and location of demand
Mixed-integer linear programming
GAMS/CPLEX
Bozorgi-Amiri et al. [22]
Preparation and response phase(location-allocation model)
Demand and location,  cost of procurement and transportation
Mixed-integer nonlinear programming
(Robust optimization)Particle swarm optimization metaheuristic
Chang et al. [16]
Preparation and response phase(location-allocation model)
Demand and location of demand
Mixed-integer programming
Sample average approximation schemeESRI ArcGIS 9.x to perform spatial data analysis
Rawls and Turnquist [91]
Preparation and response phase(location-allocation model)
Demand and route reliability
Two models(1) Mixed-integer linear program(2) Mixed-integer nonlinear program
Lagrangian L-shaped method heuristic algorithm
Rawls and Turnquist [19]
Preparation and response phase(dynamic allocation model)
Short term demand and demand location
Mixed-integer programming
Computation of the solution with CPLEX
Horner and Downs [23]
Mitigation phase(warehouse location model)
—
Mixed-integer programming
C++ programming, spatial data managed by GIS and TransCAD simulation
Dekle et al. [15]
Mitigation phase(maximal covering location model)
__
Integer programming
Excel IP solver, ArcVIEW, GIS software simulation
Zhang et al. [24]
Response phase(resource allocation model)
—
Integer programming
Local search heuristic
Barzinpour and Esmaeili [27]
Preparation and response phase(location-allocation model)
—
Mixed-integer linear programming
Goal programming
Rawls and Turnquist [92]
Preparation phase(facility location model)
Demand
Mixed-integer linear program
Lagrangian L-shaped method heuristic algorithm
Jia et al. [29]
Response phase(maximal covering model)
—
Integer programming
Genetic algorithm and Lagrange relaxation heuristic
Murali et al. [28]
Response phase(maximal covering location model)
Demand
Mixed-integer nonlinear programming
Coded simulated annealing heuristic in C++ using Microsoft Visual Studio 5.0
Lin et al. [20]
Response phase (location-allocation model)
—
Mixed-integer programming
Two-phase heuristic approach
Tricoire et al. [93]
Recovery phase(maximal covering location model)
Demand
Mixed-integer nonlinear programming
Epsilon-constraint algorithm
Yushimito et al. [94]
Preparation phase(facility location model)
—
Mixed-integer nonlinear programming
Voronoi-based heuristic algorithm
Abounacer et al. [26]
Response phase(location-transportation model)
—
Mixed-integer programming
Exact algorithmTable 5
Network design and relief distribution models in HSC.
Author
Disaster phase
Uncertain component
Model formulation
Solution technique
Sheu et al. [95]
Response phase(resource allocation model)
Fuzziness of demand and supply
Integer programming and fuzzy linear programming
Hybrid fuzzy-optimization
Ahmadi et al. [31]
Preparation and response phase(location-routing model)
Road destruction location due to earthquake
Mixed-integer nonlinear programming
Metaheuristic(variable neighborhood search) coded in MATLAB
Yi and Kumar [78]
Response phase(multicommodity network flow model)
Route reliability
Mixed-integer multicommodity network flow model
Metaheuristic(ant colony optimization)
Vitoriano et al. [96]
Response phase(relief distribution model)
Ransack probability of the link
Mixed-integer linear programming
Goal programming
Tzeng et al. [37]
Response phase(relief distribution model)
Demand
Fuzzy multiobjective programming
Fuzzy programming
Chen et al. [42]
Response phase(vehicle routing model)
—
Mixed-integer nonlinear programming
GIS based simulation
Wang et al. [72]
Response phase(location-routing model)
Route reliability after earthquake
Mixed-integer nonlinear programming
Comparison of NSGA and NGDA algorithms using MATLAB
Jabbarzadeh et al. [33]
Response phase(aid distribution model)
Blood demand
Mixed-integer nonlinear programming
Robust optimization
Liberatore et al. [43]
Recovery phase(relief distribution model)
Route reliability and complete distribution without any loss
Mixed-integer nonlinear programming
RecHADS algorithm
Campbell et al. [32]
Response phase(vehicle routing model)
—
Mixed-integer programming
Heuristic(insertion and local search)
Sheu [38]
Response phase(relief distribution model)
Effect of emotional contagion on the attitude of survivor
Mixed-integer linear programming
Structural equation modelling
Afshar and Haghani [39]
Response phase(integrated supply chain model)
—
Mixed-integer programming
CPLEX
Huang et al. [79]
Response phase(assessment routing model)
—
Integer programming
Tabu search heuristic coded in C++
Bozorgi-Amiri et al. [36]
Preparation and response phase(relief distribution model)
Demand and location of demand
Mixed-integer nonlinear programming
Robust stochastic optimization approach, compromise programming
Özdamar and Demir [64]
Preparation phase(vehicle routing model)
—
Mixed-integer nonlinear programming
Hierarchical cluster and route heuristic
Lin et al. [40]
Response phase(vehicle routing model)
—
Mixed-integer programming
Decomposition and assignment heuristic
Wohlgemuth et al. [34]
Response phase(pickup and delivery model)
Demand
Mixed-integer linear programming
Tabu search heuristic
Vitoriano et al. [41]
Response phase(relief distribution model)
Route reliability
Mixed-integer nonlinear programming
Goal programming
Shen et al. [35]
Response phase(vehicle routing model)
—
Mixed-integer programming
Tabu search heuristic
Chiou and Lai [51]
Response phase(shortest path, traffic assignment, and traffic control model)
Route reliability
Mixed-integer nonlinear programming
Fuzzy set theory and genetic algorithm
Berkoune et al. [71]
Response phase(relief distribution model)
—
Mixed-integer programming
Genetic algorithm(greedy heuristic, set enumeration heuristic)
Adıvar and Mert [52]
Response phase(international disaster relief distribution model)
Meeting the needs on time
Mixed-integer linear programming
Fuzzy programmingTable 6
Mass evacuation models in HSC.
Author
Disaster phase
Uncertain component
Model formulation
Solution technique
Sheu and Pan [55]
Response phase(centralized emergency supply network model)
—
Mixed-integer linear programming
Numerical case study on LINGO
Bish [97]
Preparation phase(bus-based evacuation model)
—
Mixed-integer linear programming
Heuristic algorithm
Bretschneider and Kimms [68]
Response phase(evacuation model)
—
Mixed-integer programming
LP-based heuristic
Sayyady and Eksioglu [67]
Response phase(no-notice urban evacuation model)
—
Mixed-integer linear programming
Simulation package, Dynasmart-P, Tabu search algorithm
Chiu and Zheng [58]
Response phase(dynamic traffic evacuation model)
—
Linear programming
Algorithm
Naghawi and Wolshon [57]
Response phase(citizen-assisted evacuation model)
—
TRANSIMS agent based simulation
GIS and TRANSIMS based simulation
Hsu and Peeta [60]
Preparation and response phase(stage-basedevacuation model)
—
Mixed-integer nonlinear programming
Branch-and-bound algorithm
Özdamar and Yi [61]
Response phase(evacuation and logistics support model)
—
Mixed-integer programming
(Greedy neighborhood search, Path-Builder  heuristic)Implemented in C++
Yi and Özdamar [65]
Response phase(location-distribution model for logistics support and evacuation operations)
Occurrence of disaster
Mixed-integer linear programming
Routing algorithm
Najafi et al. [63]
Response phase(integrated HSC model)
Number of injured people, demand
Mixed-integer linear programming
Exact methodology
Kongsomsaksakul et al. [62]
Response phase(location-allocation model for flood evacuation planning)
—
Mixed-integer nonlinear programming, bilevel programming
Genetic algorithm
Özdamar and Demir [64]
Preparation phase(aid distribution and evacuation model)
—
Mixed-integer nonlinear programming
Multilevel clustering algorithm
Hsu and Peeta [59]
Response phase(information based control evacuation model)
Demand
Mixed-integer nonlinear programming
Fuzzy programmingRQ 2. What Are the Unexplored Areas in the Field of HSC? After reviewing the selected research papers, the following gaps have been suggested which will provide a future research agenda for researchers in the field of HSC.In the current era, global warming is changing the temperature of the Earth. Consequently, weather is becoming severe and causing more natural disasters with the passage of time. Adoption of green supply chain in HSC is inevitable in order to minimize the effect of global warming. Transportation of relief goods is the major part of HSC that requires use of fuel. Just after the disaster for the first 72 hours in response phase of HSC, primary goal of minimizing delivery times is considered. During these first 72 hours of response phase the factor, which type of fuel is being used and how it would affect the environment, is ignored. Most of the time, the urgency of the disaster requires the use of fuels that result in an increased carbon percentage in environment. A lot of research has been done in commercial supply chain with the purpose of sustainability but in this perspective the HSC is still an unexplored research area.After the occurrence of disaster two main tasks are initiated: the first one is to evacuate the injured people from disaster location and the second one is to assess damage caused by disaster. Special designed equipment is required for both of these tasks. For evacuation of injured people, one needs specially designed equipment because the environment after a disaster is not always suitable for human beings. One example of such incident is the Fukushima nuclear power plant accident where it was impossible to enter into the plant due to harmful radiations. In this situation, robots can be used to control any activity which can be harmful for human beings. The DARPA Robotics Challenge (DRC) is a contest in which teams from all over the world participate. This is a competition of robot systems and software teams seeking to develop robots capable of assisting humans in responding to natural and man-made disasters. Although it requires very high expertise in the field of robotics, it needs time to make interdisciplinary collaboration research work to serve mankind. Some authors like Hirose and Fukushima [81], Tanaka et al. [82], and Rosen [83] have made contributions in this area. The second important task in the initial stages of the disaster is to assess the damage. For this purpose, currently satellite generated pictures are used but they take a lot of time and they provide general information. In the last few years, unmanned aerial vehicles (UAV) have helped in postdisaster damage assessment operation. Though few contributions in this area include Bendea et al. [84], Thiels et al. [85], and Tuna et al. [86], still this area has a lot of potential to serve HSC.Consideration of gender related needs is a factor not yet considered in HSC models. Gender determines responsibilities and powers associated with being male and female. As basic health services are part of emergency relief, both women and men should be catered in a gender sensitive manner. Emergency health services personnel should be trained to consider gender related needs. Along with this, the factor of special-risk populations that include children, the elderly, and people with disabilities should also be considered. By considering these factors, researchers can adopt more realistic approach for HSC research.Debris management is another unexplored area in the field of HSC. In the short term, debris hampers transportation of relief goods to disaster stricken areas while, in the long run, it poses serious threats to the environment and human health. Until debris is cleared, relief operation cannot be performed in an effective manner. According to Van Wassenhove and Pedraza Martinez [87], after the Hattian earthquake, in spite of the availability of excessive aid from international agencies, relief operations were not as effective as they should have been. This was due to the unavailability of debris removing equipment that caused huge bottlenecks in terms of access. In some models of the medical supply chain, the assumption that roads are not blocked after the disaster and therefore all hospitals are accessible seems unrealistic. In the area of debris management, Çelik et al. [88] addressed debris management from the HSC point of view. Debris management consists of three phases. The first one is debris clearance, in which debris is cleared from roads and placed on the road sides to start the relief operation. The second phase is debris collection, in which contractors are given contracts to pick up the debris from road sides. In the third phase, debris is assigned to debris processing units. Plant locations for debris processing and debris collection contracts are the areas that require special attention of the researchers.Contracts in the commercial supply chain are very effective, because by using different types of contracts, all stakeholders can optimize their profit. The first step to design a contract is to minimize an overall cost in collaboration and the second step is to allocate the reduced cost among stakeholders. In the HSC, primary goal is saving lives and providing a better service to the affected people. Suitable contracting in HSC relief may help in both minimizing cost and increasing level of service. Contract management in HSC is another dimension that can provide research agenda for researchers. Only few research papers that consider contracts for HSC include Liang et al. [89] and Sheu and Pan [90].In the commercial supply chain for a product, either push strategy or pull strategy is selected. In some cases a combination of both push and pull strategy, with a defined boundary, is used. In the HSC context, a very little contribution has been made regarding selection of strategy for relief goods. Until now, no dedicated model has been introduced regarding push-pull boundary for relief goods. For disasters that are noninstantaneous like hurricane, flood, famine, and volcanic eruptions the concept of push-pull boundary can be used on short term basis. We can estimate tentative disaster time of occurrence by using technology. For example, in case of hurricane meteorologist can forecast the hurricane route and time when hurricane will hit a particular area. During this estimated time, which is named as “lead time for disaster,” the concept of push-pull boundary can be implemented. Utilization of this concept in HSC will decrease shortages and inventory holding cost, which would be helpful to perform disaster relief operations in a cost effective manner. Viewing this property of lead time, general guidelines for defining push-pull boundary for each type of disaster need to be defined.In a disaster relief operation, procedures never go according to the way they are planned. Customs clearance at airport or seaports, fuel shortages, vehicle shortages, and blocked roadways are few among many factors that can cause disturbance in the planning of the disaster relief operation. The situation can be worse if a country is landlocked because it will increase lead time and the uncertainty factor. As a result, one distribution center can have excessive inventory of relief goods which may contain food or medicine. It is obvious that with the passage of time this food or medicine will decay. Some vaccines require specific temperatures and have a particular expiry date. In this aspect the research on the issue, how to handle such goods that have short life and involve quick decay, is still absent.Contrary to labor of commercial supply chain, labor in HSC is not similar throughout the relief operation. HSC requires skilled logisticians for planning, paid workers for normal routine responsibilities, and volunteers in the time of disaster. Each type of labor has its own parameters. As disasters are unpredictable, it is really tough to retain qualified personnel for the long term. Similarly, lack of trained staff and workers is another important issue. Volunteers may join or leave the disaster relief operation in dynamic fashion. If this versatile nature of labor is not managed properly relief operation cannot be performed in an effective manner. So, human resource management in the perspective of HSC is also an area that requires attention of researchers.


## 5. Conclusion

The ultimate goal of this research is to address the following questions: (1) What is the current status of research on the mathematical models in the field of HSC? (2) What are the unexplored areas in the field of HSC? To answer the first question, research in HSC was broadly classified into three classes of facility location, relief distribution, and mass evacuation. All other topics of HSC with little variation have been discussed under these three classes. A summary of modelling techniques, solution methodologies, with their pros and cons, has been presented. It is observed that HSC research focuses more on preparation and response phase as compared to mitigation and recovery phase. Very few models have adopted an integrated approach in the HSC because it increases computational difficulties. However, recent advances in optimization algorithms have enabled researchers to solve such complex models. It was noticed that most of the models developed in the HSC are deterministic in nature but some authors adopted stochastic or fuzzy approach. A summary of the uncertain factors considered in stochastic and fuzzy models is provided. Finally, in the second question the unexplored research areas in the HSC are provided for future research agenda. It is noticed that many researchers assumed unrealistic situations in their model formulation. Unrealistic assumptions cause a limitation in implementation of these models in real life; that is why it is necessary to develop practical and feasible models so that research can be filtered down into policy, practice, and procedures.A limitation of this review study is that it only considers the research that uses any type of mathematical technique. Other papers that have proposed any type of framework for HSC and theoretical in nature are out of scope of this study. Secondly, this study only considers the research papers written after 2004 because most of the research in HSC was done after Asian tsunami 2004. Finally, this study can be fruitful for the researchers and practitioners to utilize the proposed modelling techniques and solution methodologies and suggest suitable solution methodologies for unexplored research areas in the field of humanitarian supply chain.



---

*Source: 3212095-2016-02-17.xml*