Modeling Multi-Objective Optimization with Updating Information on Humanitarian Response to Flood Disasters

: Unpredictable natural disasters brought by extreme climate change compound difﬁculties and cause a variety of systemic risks. It is thus critical to provide possibilistic scheduling schemes that simultaneously involve emergency evacuation and relief allocation. But the existing literature seldom takes emergency evacuation and relief supplies as a joint consideration, nor do they explore the impact of an unpredictable ﬂood disaster on the scheduling scheme. A multi-stage stochastic programming model with updating information is constructed in this study, which considers the uncertainty of supply and demand, road network, and multiple types of emergency reliefs and vehicles. In addition, a fuzzy algorithm based on the objective weighting of two-dimensional Euclidean distance is introduced, through moderating an effect analysis of the fuzzy number. Computational results show that humanitarian equity for allocating medical supplies in the fourth period under the medium and heavy ﬂood is about 100%, which has the same as the value of daily and medical supplies within the ﬁrst and third period in the heavy scenarios. Based on verifying the applicability and rationality of the model and method, the result also presents that the severity of the ﬂood and the fairness of resources is not a simple cause-and-effect relationship, and the consideration of survivor is not the only factor for humanitarian rescue with multi-period. Speciﬁcally, paying more attention to a trade-off analysis between the survival probability, the timeliness, and the fairness of humanitarian service is essential. The work provides a reasonable scheme for updating information and responding to sudden natural disasters ﬂexibly and efﬁciently.


Introduction
Various disasters (e.g., natural disasters, unconventional events, public health emergencies, etc.) have frequently occurred in recent years. According to data from Emergency Events Database [1], there were 999 natural disasters in Europe, focusing on the last 20 years (2001-2020). The most prominent ones are floods, storms, and extreme temperatures. It had killed over 150,000 people, affected over 11 million others, and cost over 217 billion dollars. Flood disaster characterized by the suddenness and interruption of disaster, poses a significant threat to the lives and welfare of entire populations [2][3][4]. Therefore, designing the multi-objective optimization by updating flood information, which consists of emergency evacuation [5][6][7] and supplies distribution [8] in dynamic flood scenarios, is a big challenge for decision-makers.
When a disaster occurs, the immediate evacuation of a person is pivotal in minimizing the possible risk of losses. Shahparvari et al. [9] aimed to maximize the number of evacuees with time windows and disruption risks under uncertainties. The succession of humanitarian evacuation depends on the capacity of available teams and fleets. Besides, relief items are only one kind of resource required for coping with the flood [3], which conduces to support evacuation activities and has increasingly drawn attention. Elçi and Noyan [10] discussed the state of the road through emergency transportation time in different scenarios.
Meanwhile, the research has turned its attention to the inter-connectivity and interdependency in various types of networks [11], network recovery and emergency response [12], as well as facilities location-routing problems with uncertainties [13]. In 2020, for instance, a large excess stock of supplies did not normally deliver as infrastructure had been damaged or destroyed [14,15]. But most of the literature assumes that no limitation exists in the state of the transportation network, the availability of vehicles and work crews, which in turn heavily obstructs the accessibility of schemes [16,17]. Moreover, a secondary disaster of the flood may occur, road network and evacuee population could be very hard to analyze, especially when considering the flood information. It is thus vital to research the problem of emergency evacuation and relief allocation simultaneously.
The existing work seldom takes emergency evacuation and relief supplies as a joint consideration, nor do they explore the impact of an unpredictable flood disaster on the scheduling scheme. Questions in research are thus proposed: How to develop a multi-objective collaborative model with updated flood information that involves humanitarian evacuation and supplies distribution, according to the actual requirement of the disaster?
How to determine the optimal decisions to ensure that the trapped people and emergency reliefs are quickly and timely allocated to the nodes in the flood event?
Motivated by the literature and the actual situation of humanitarian rescue, the contexts of this study are included as follows: Firstly, a multi-stage stochastic programming model with the changed flood information is constructed. It not only applies Bayes's theorem to discuss the variability of disasters, but also introduces the interval number and the triangular fuzzy number to describe the fuzzy uncertainty of variables. Secondly, in an attempt to systematically discuss the fuzzy uncertainty of functions and constraints in the model, the deterministic transformation of the fuzzy number, and a fuzzy algorithm that based on the objective weighting of two-dimensional Euclidean distance can be mainly involved in the method. Thirdly, the effectiveness and feasibility of the model and method can be verified by a case study, which is significantly impacted by heavy rains.
In brief, the contributions are as follows: In the first place, unlike the traditional study about the flood disaster, this work not only simultaneously considers emergency evacuation and relief allocation, but only integrates it into a multi-objective collaborative model with updating information in emergency response to the disasters. In addition, a methodology seamlessly and simultaneously maximizes problem-solving ability, which could effectively and efficiently dispatch emergency crews and supplies under the uncertainties of humanitarian demands and supplies.
This research is organized as follows: Section 2 briefly introduces an overview of the literature. This is followed by Section 3, where a mathematical model is given based on problem description and assumption. Section 4 demonstrates the processes of an efficient method. Results and discussion of case study can be drawn in Section 5, whilst Section 6 introduces the conclusion and future work.

Literature Review
Two streams of emergency evacuation and allocation and the main differences of the research have been presented.
After disaster strikes, emergency evacuation, a complex and systematic project, is essential in influencing humanitarian response. Scholars have independently studied or combined emergency evacuation behavior (e.g., the evacuee demographics, route preference, route diversion, warning signals, etc.) and modeling [5][6][7]. For example, Bayram [7] analyzed the traffic assignment models, research objectives in evacuation modeling (e.g., the total evacuation or clearance time minimization, total shelter cost, shelter coverage, etc.), and the evacuee behavior issues (e.g., gender, ethnicity, perceived risk, etc.). Meanwhile, few authors have applied technologies to obtain real-time traffic network, which is easily damaged by emergencies. Liebig et al. [18] studied dynamic route planning based on real-time traffic information. Rodríguez-Espíndola and Gaytán [19] used a GIS analysis to discard the facilities and identify road failures. Research is also largely related to work teams that could quickly relocate evacuees to protected or less vulnerable shelters. Vkbari and Salman [20] provided a method to generate a synchronized work schedule for the road-clearing teams. Li et al. [21] integrated the logistic support scheduling with repair crew scheduling and routing problem.
To make emergency rescue more practical, more and more attention is being paid to allocating humanitarian supplies, which ensures the evacuated individuals' daily lives. It not only uses the aggregated coverage model, maximum coverage model, and P-center model, but also systematically discusses the capacity constraints of shelters or vehicles. However, previous research often applies a single type of tool, and ignores the availability of the network in the hardest-hit areas after severe destructive flooding. It is argued that multitype vehicles could be dispatched according to the characteristics of humanitarian resources. Alem et al. [22] developed a two-stage stochastic network flow model for decisions on the fleet size of multiple types of vehicles over a dynamic multi-period horizon.
Apart from focusing on one of the above issues, researchers have gradually paid attention to the impact of facilities networks on humanitarian allocation in the post-disaster. For example, Tofighi et al. [8] aimed to minimize the cost of stock pre-positioning, and the facility located in the first stage, whereas they sought to minimize total distribution time, the total cost of unused inventories and the unmet demands in the second stage. Renkli and Duran [23] used a scenario-based formulation to minimize the weighted probability of road blockage, whereas Xu et al. [24] studied a balance between the demands and restoration vehicle reliefs based on real-time disruption scenarios. Several authors have also constructed multi-stage stochastic models to coordinate the procession of emergency allocation. Rennemo et al. [25] proposed three-stage stochastic facility routing model to maximize the combination of demand fulfillment and unused budget.
In summary, the possibilistic scheduling schemes for emergency evacuation and relief distribution are rarely discussed in the existing studies. The main differences in this research can be seen as follows: Most of the literature takes an earthquake, typhoon, snow disaster, and other emergencies as the background, but there are few reviews on flood disasters. Emergency rescue, which involves the evacuation and relief allocation in the changed flood scenarios, is studied. In addition, some reviews often regard cost or time as the common objective. The performance of victims and the fairness of reliefs allocation are rarely considered simultaneously. A multi-objective collaborative model of emergency rescue with updated information is established, and the triangular fuzzy number and the interval number can be used to discuss the uncertainties of the demands and supplies. Furthermore, the realistic case is applied to verify the reliability of the approach, which focuses on minimizing emergency service time, the unfulfilled rate of relief and maximizing the saving of lives.

Problem Description
As an essential aspect of disaster response, emergency evacuation and allocation have been challenged by many factors [26,27], including the intensity of the disaster, the uncertainty of demand, and the availability of networks and vehicles. Based on the realistic circumstance of a disaster, and the multi-type of reliefs and vehicles, a multi-objective optimization problem with multiple periods in the changed flood scenarios is presented. To simplify the transmission dynamics of the model, the affected nodes, alternative roads, humanitarian service teams and vehicles are considered in this study. In the first place, victims in the trapped nodes are serviced by work crews, which are dispatched by the rescue nodes. And residents in the more serious areas are moved to the resettlement nodes.
In an attempt to maximize the saving of lives and property, based on the feasibility of the service path, the ready availability of emergency evacuation is improved through the work crews of the rescue nodes. Besides, the multi-type vehicles, which take more responsibility for delivering daily and medical resources to the trapped and resettlement nodes, are assigned by emergency distribution center. It is how headquarters perform these activities in the stipulated time, space, and resource constraints, to minimize the timeliness and fairness of emergency service.

Notations and Assumptions
The graph G = (N, A) of which N (A j ∪ B e ∪ D i ∪ D l ∈ N) is the set of nodes, and A is the set of links. There are various nodes of emergency work teams A j ∈ N, where each crew j ∈ A j ; A set of the trapped points in the flood B e ∈ N, where each node e ∈ B e ; Demand nodes (D i ∈ N) of which they are either the trapped nodes (i ∈ D i ) or the resettlement nodes (w ∈ D w ∈ D i ); A set of emergency distribution centers D l ∈ N, where each reliefs center l ∈ D l . In addition, the work denotes the path of acr (l, i) as c li , marks the road of acr (j, e) as c je . At the same time, the observed time variables t ∈ D t can be presented. The indices, sets, parameters and decision variables are given below (Table 1): Table 1. Relevant parameters are presented in this model.

Sets and Indices Definitions
j ∈ A j A set of the rescue nodes of work teams; e ∈ B e A set of the trapped sites in the disaster; The interrupted or damaged road; c ∈ C Categories of reliefs; v ∈ L Sets of emergency service vehicles; ϕ ∈ δ The possible set of flood scenarios;

H
The fixed time between two adjacent observation moments; d li The distance from the node l to node i; ε The critical value of the waterlogged road which meets the traffic condition, 0<ε<1; γ The threshold value of whether the road network is passable, δ <γ<1; δ The critical value of the waterlogged road, which is restorative, ε< δ <1; η 1 The corrected parameter of deterioration coefficient; The correction factor of repairability; ∂, τ, ϑ Coefficients affected by the connectivity of road network; S j t je The survival probability of the trapped node e that was impacted by emergency rescue node j; ∂ t Results of the observed information ∂ at time t; p ϕ The prior probability of flood disaster; The probability of the flood scenario ϕ, which can be understood as the posterior probability of disaster by using the Bayesian theorem; The occurred probability of observable information ∂ under the flood scenario ϕ at time t; a ke The waterlogged degree of road c je where emergency team k that dispatched from the rescue point j to node e; d c The fuzzy demand for reliefs c per capita, The fuzzy capacity of emergency crew k, Number of victims in the trapped node e at time t; N i t The actual number of people in node i at time t; N w t The upper limit that node w can accommodate victims at time t;

R lic
Quantity of relief supplies c distributed by emergency distribution center l to node i, R lic = (r − lic , r + lic ). The upper and lower limits of the interval numbers are denoted by r − lic and r + lic , respectively; Q v The loaded capacity of each type v of emergency service vehicle; M lv Number of the type v of vehicles that dispatched by emergency distribution center l;

R l
The capacity of emergency distribution center l (l ∈ D l ) for dealing with reliefs. R l = (r − l , r + l ), r − l and r + l parameters explain the upper and lower limits of the interval numbers; M j Number of emergency rescue teams deployed from the node j;

Sets and Indices Definitions a li
The waterlogged degree of c li between node l and node i; Rescuing time of emergency rescue team k from node j to the node e; t livc Traveling time of the type v of vehicle τ to transport reliefs c from node l and node i; The actual time of emergency rescue team k that aims to rescue the node e; t li The restorative time from node l and node i; The rehabilitation time of road network from node j and node e; t e The time point when the survival probability of victims in the node e is 0.5; v lic v The average speed of the type v of vehicle τ, which contributes to transporting the supplies c from node l to node i; The remained quality of reliefs supplies c after arriving the node i by the type v of vehicle τ; x ke x ke = 1, if emergency work team k is responsible for the trapped node e; 0, otherwise; x keg x keg = 1, if emergency crew k completes humanitarian work for the trapped node e, and still moves to the next node g; 0, otherwise; y je y je = 1, if the trapped node e is processed by emergency rescue crew, which is assigned by the node j; 0 otherwise; x livc x livc = 1, if a node i ∈ D i is visited by the type v of vehicle τ that loading reliefs c, which is assigned by reliefs center l; 0, otherwise; Considering the complexity and uncertainty of the flood disaster, a clear set of assumptions underpinning the model ought to be illustrated to simplify its scope and a set of variables. The main assumptions are as follows: 1.
The probability (p ∂|ϕ ) of the certain observed information in the flood.

2.
Mixing of the same kind of resources is allowed, but the mixing of multiple supplies is unworkable.

3.
Information regarding the road network in the hardest-hit areas can be done by satellite images, helicopters or drones. 4.
The number of people in the affected area is known, and the demand at each node has a linear function of its capacity.

5.
Each crew departing from the rescue nodes must visit at least one trapped node through its route, and each trapped node should be visited by at least a team at the end of the task.

Model Formulation
To integrate emergency rescue efficiently and precisely into a multi-objective collaborative model with updated information. Considering the fuzzy uncertainty of emergency demand and supply as well as the road network, the triangular fuzzy number, Bayes's theorem and the interval number are introduced to construct the upper-level and lower-level modeling with multiple periods.

The Upper-Level Modeling
As an important aspect of the design of the upper-level model is to rescue the trapped people, its thus necessary to maximize the saving of lives in the flood scenario ϕ, which can be expressed as follow: It is indicated that the prior probability of flood disaster and the occurred probability of the observed information that denoted by the p ϕ and p ∂ t |ϕ , respectively. The occurred probability p ϕ t which can be understood as the posterior probability in flood is presented: In view of the index function of the survival probability of victims, which is proposed by Fiedrich et al. [28]. A piecewise function is illustrated as Equation (3), which shows the probability of evacuating people in a flood disaster.
The above function is inflected by t e and t ke variables. If one of the condition t ke ≤ t e , 1 − 0.5(t ke /t e ) 2 is feasible for the survival probability of victims in flood disasters. S e (t ke ) is calculated by 0.5e 2(1−t ke /t e ) , otherwise.
Disaster-formative environment refers to the conditions and surroundings where flood disasters occur, determined by the factors that mainly result from the combination of climate variables and underlying surfaces. The connectivity of road networks which impacted by disaster-formative environments may be struck or damaged. And the transferred time of the rescue team is easily affected by a part of the waterlogged road. Ajam et al. [29] analyzed the issue of minimizing the total latency of the nodes with a single work team. Based on analyzing various parameters such as t jke , t je , η 1 , η 2 and a ke , the most important aspects of the time of work crews that assigned to the trapped nodes can be captured: ∑ e∈B e ∑ k∈M x ke ≤ M j , ∀j ∈ A j , ∀e ∈ B e , ∀k ∈ M Constraints (5) and (6) demonstrate that each rescue team is sent to a maximum of one node at a time, and emergency crews dispatched do not exceed the maximum number.
∑ k∈M ∑ e∈B e f k x ke ≥ ∑ t∈D t ∑ e∈B e N e t , ∀e ∈ B e , ∀t ∈ D t , ∀k ∈ M (7) Constraint (7) expresses that the capability of emergency crew shall not be lower than the population of victims. The number of individuals that resettled to each node is not greater than its maximum, as illustrated in Constraint (8).

The Lower-Level Modeling
Satisfying emergency demand of nodes is crucial to the success of humanitarian relief allocation, as a lack of relief supplies may cause suffering and life loss for victims [30,31]. Fairness is included in the lower-level model that uses the performance measure related to the filling rate. This max-min objective function could be classified as an equality-based method [32]. Equation (13) is to maximize the minimum level of dissatisfaction on the condition of the flood scenario ϕ, thereby balancing the fill rate across them.
With the consideration of various variables such as the observation moments and waterlogged scenarios, the objective function of minimizing the service time of relief supplies is illustrated in Equation (14): Many methodologies have been applied to discuss the damaged network [33,34] and vehicle travel time [35]. For instance, Baskaya et al. [34] explored the different route distances between relief centers and affected locations to reveal the disruption levels of road infrastructure. Humanitarian allocation greatly depends on the available fleet of emergency service vehicles and their capacity, the arrival time of vehicles is mainly affected by the distance of the flood-stricken areas. When the t livc is identified in Equation (15), which is associated with d li , v lic v and the coefficients that impacted by the connectivity of road network in the flood disaster. S.t.: Constraints (16) and (17) are mainly analyzed for the capacity of emergency distribution vehicles and service centers. ∑ l∈D l ∑ i∈D i x livc ≤ M lv , ∀c ∈ C, ∀i ∈ D i , ∀l ∈ D l , ∀v ∈ L The distribution of relief supplies to the affected nodes does not exceed its demand, as depicted in Constraint (20). 0 < ∑ l∈D l ∑ c∈C R lic ≤ ∑ t∈D t ∑ i∈D i ∑ c∈C N i t d c , ∀l ∈ D l , ∀c ∈ C,∀i ∈ D i , ∀t ∈ D t (20) Constraint (21) defines the limited time of distribution service.
Constraint (22) refers to the relationship between resources within a vehicle when it goes through two nodes before and after.
Constraint (23) demonstrates the domain of decision variables.

Methodology
An efficient task does depend on the coordination between emergency evacuation and reliefs allocation, but addressing the multi-stage stochastic programming model with updating information in parallel gives a challenge in comparison to the single-objective case. Several studies in recent years have focused on the model and algorithm to improve the capability of humanitarian service after disruptive events, such as the heuristic algorithm [36,37], the random search algorithm [38] or other approaches [39]. Nevertheless, these methodologies are less used in a multi-objective collaborative model with the updated flood information. Instead, higher stability and reliability of the fuzzy algorithm based on the objective weighting of two-dimensional Euclidean distance is applied, considering the blurred boundaries, the urgency of reliefs, and the mutual and interconnected influences among objectives.
To deterministically analyze the fuzzy uncertainty of functions and constraints, which are transformed by the interval number and the triangular fuzzy number, the deterministic transformation method of fuzzy number is first utilized. A triangular fuzzy number is seen as a fuzzy set in the studied domain, and the interval number is included in the set, which consists of all real numbers in the range of a closed interval. The objective function (13) and some constraints can be updated as follows: Moreover, based on the discussion of deterministically transforming the fuzzy number within the model, the specific ideas and processes of an applied algorithm is seen in Figure 1. Main steps in the approach are: (1) Based on the discussion and definition of simulating level, the function of three objectives in the mathematical model is rephrased as: (2) Finding out the maximum and minimum of objectives, at the given constraints of modeling, the upper and lower bounds of them noted as sup{z i (x)} and inf{z i (x)}, respectively.
(4) According to various factors such as the severity of a disaster, emergency supply and demand in each period, the decision preference coefficients are jointly determined by humanitarian experts and deciders. Denoting the coefficients of decision preference as w 1 and w 2 . It is noteworthy that the effective solution of modeling could be as close as possible to the ideal value.
where b and c mean the preferred genera of the negative and positive solutions of the fuzzy, respectively. (5) Combining the solution method and the preference coefficients of deciders, the multi-objective collaborative model can be transferred into a single one.
(6) By constructing the Lagrange function, the above model is analyzed by Equation (37). To obtain the multiple objectives, it is necessary to find the partial derivative for w i , x i , σ, namely.

Data for Numerical Studies
In this subsection, the emergency rescue of rare sudden heavy rainfall in A District, Henan Province, China, in July 2021 is an analysis example. The rainfall volume of the case has broken the historical record in mainland China [40,41], which has seriously threatened people's lives and the safety of the road. As shown in Figure 2   Flood process segmentation rules identify the beginning and ending times of flood phases considering both natural and social aspects and divide the flooding process into four phases: latency, onset, development, and recovery [42]. Increases in rainfall return periods and rainfall durations in the last three stages could lead to more severe flooding. It is of great significance to analyze the impact of the road network of the affected nodes on the runoff process in combination with the characteristics of the flood. The distance among the flood-stricken areas is seen in Figure 4, whereas the damaged level of road network in an extreme rainstorm is defined in Figure 5. Two commodity types, which require multi-type of tools (emergency daily and medical vehicles), are supplied by the relief centers (DL, SD, DG). The fuzzy demand per capita for daily relief and the medical item is [1,4,5] and [0, 2, 3], and the average weight of each pack is 5 kg and 0.36 kg.  In urban flood disasters, the disaster drivers primarily refer to heavy rain expressed by rainfall duration and accumulated rainfall [41]. At the same time, to accurately describe the difference in rainfall throughout the flooding process, the scenarios of the flood are divided into mild, moderate, and severe groups. By comparing with the official data released by the government [41,43,44] and the result of expert consultation, the relevant parameters settings for p ∂|ϕ are represented in Table 2. It assumes that humanitarian managers not only may make a direct decision without observing the flood information at the first time, but also can carry out the rescue efforts. In the case of data for some variables that cannot be published, according to the review of relevant literature [40,41], the investigation and consultation carried out in the scope of the research, the relevant parameters are set in a way that combines the simulated data are seen in Figure 6 and the Table 3, respectively.

Ideal Weight of Sub-Objectives for Humanitarian Rescue
To obtain the scheduling schemes for the flood disaster, the first step is to discuss the ideal weight of sub-objectives for humanitarian rescue. The objectives of the model with multiple periods in the flood scenario ϕ are calculated, with the constraints of decision preference (W1 = 0.10, W2 = 0.50, W3 = 0.40; W1 = 0.40, W2 = 0.30, W3 = 0.30; W1 = 0.70, W2 = 0.20, W3 = 0.10) which follow ∑ i∈n w i = 1. In the light, medium and heavy scenarios, the trend in the total mortality and unfulfilled rate in the first, second and fifth periods can be randomly selected (Please see Figure 7). In the process of humanitarian rescue, the coefficient of decision preference has an important influence on the choice of scheme. As the weight of the survival probability gradually increases, the determination to minimize the mortality rate rises progressively. It increases the probability of survival within periods under the low, medium, and heavy disasters. Additionally, there is a tendency to fall and rise, when the unsatisfied logistics service gradually decreases. When the weight of the survival coefficient is greater than the unfulfilled coefficient, it implies that relief agencies focus on minimizing the mortality rate. Emergency performance is related to the capability of work teams and the updated information. On the contrary, emergency deciders are inclined to the fairness of relief allocation, with the consideration of the flood, victims, and tools information.
Emergency distribution time for each period in the low, medium, and heavy flood disaster is depicted in Figure 8, panels (a), (b) and (c), respectively. Following the decrease in service time preference, the requirement of planners for minimizing time is not very strong. It indicates that the target of emergency distribution time in multiple periods progressively decreases and then increases. What's more, there is a different tendency from the first period to the fifth one in the same scenario. The reason why humanitarian time in the severe flood is significantly lower than others is that the value of p ϕ t is generally smaller. In terms of giving the preference of subobjectives, deciders should be scientific in choosing the "degree", and pay more attention to the dynamic information about the flood disaster, actual demand, and transportation network. And with the consideration of the uncertainty and different durations of the whole flood process [45], the planners could reduce barriers to exchange between regions and promote the circulation efficiency of elements in various areas [42].

Humanitarian Scheme Impacted by the Flood Disaster
The second step is to explore humanitarian schemes affected by flood information. Figure 9 presents that when W1 = 0.70, the sum of the probability of death in the flood disaster is generally larger than the remaining scenarios. It depicts that the worse the disaster, the greater threat to people's lives. Regarding scenario W1 = 0.70, the minimum mortality is located in the third period of the light disaster, the second period of the medium and heavy flood.  Figure 10 shows that when W1 = 0.10, W2 = 0.50 and W3 = 0.40, the unmet rate in the low and medium flood is bigger. Expecting that W1 = 0.70 in the heavy flood scenario, the objective trend with the constraint of W1 = 0.10 and W1 = 0.40 is more similar. The cause lies in the fact that, the unfulfilled rate is influenced by the damaged network, decision preference, and relief supplies in each cycle. As seen in Figure 11, increasing the time preference (W3) rises the objective. Although the variability of change for emergency distribution time in the low flood is greater, a similar tendency of it in the medium and heavy flood scenario can be found. The reason is that emergency distribution time in the periods under the low scenario is impacted by time preference, the length and accessibility of the network.  According to the model's results with multiple periods in the same flood scenario, when the preference coefficient of morality is 0.4, the value of it is smaller than the result with the constraint of W1 = 0.70, and larger than the rate of death when W1 = 0.10. Compared to the preference of W1 = 0.70 and W1 = 0.10, the unmet rate and emergency distribution time under the condition that W1 = 0.40 is the smallest, and its curve is more similar to the structure of the result with the constraints of W1 = 0.70. It's not only relevant to the deciders' preference, but also indirectly related to the severity of the disaster, the rescue tools, and road network. The objective of the model with the constraint of W1 = 0.40 is relatively greater than that in the preferences of W1 = 0.70 and W1 = 0.10. To further explore the optimal strategy of the multi-objective collaborative model with the condition of W1 = 0.40, W2 = 0.30 and W3 = 0.30, we start to analyze the rescue performance in the flood scenario ϕ.
As can be seen from Figure 12, there is a noticeable difference in the survival probability under the updated flood scenarios, whereas the survival probability of nodes rescued by the crews in the first period is indicated in Figure 13. The more severe the flood disaster, the higher the mortality rate of the affected nodes. The finding is not only easy to understand, but also fits with a realistic scenario. Affected by the rain island effect, the frequency of heavy rain increases in the concentrated areas [45], where the wetland area and surface water infiltration as well as stormwater absorption capacity may be reduced. It is thus important to make more efficient flood control, such as constructing flood control facilities and technologies in accordance with the principles of adapting differentiated behaviors to a local location. Besides, since emergency demand is impacted by the number of victims in different scenarios, the next step is to analyze the node's demand with multiple periods under the flood scenario ϕ. As shown in Figure 14, vague demand for medical and daily reliefs in the first period varies across the demand nodes.   On the basis of the information about the changed flood, traffic network, supplies of reliefs (shown in Table 4) and emergency tools and demand, humanitarian allocation scheme with multi-period in the updated flood scenarios can be illustrated in Table 5. Humanitarian equity for allocating medical supplies in the fourth period under the medium and heavy flood is about 100%, which has the same as the value of daily and medical supplies within the first and third period in the heavy scenarios. It represents that the severity of the flood and the fairness of resources is not a simple cause-and-effect relationship. In other words, emergency equity is impacted by the combined effect of variables (e.g., the supplies, the flood scenarios and road network, etc.) in models. Comparing the results in the low flood disaster, the time in the second period is the lowest, while the maximum time is located in the fifth period under the medium scenario. It enables the determination of an equilibrium strategy for multi-commodity, and a timely plan within various periods. And attention could be paid to the time evolution characteristics and variation laws, the alternative concepts and technologies on flood disasters.

Conclusions and Future Work
As the suddenness and destruction of the flood disasters, the limited capacity of the organization, the possibilitic and adaptive scheme is crucial to humanitarian response. Emergency evacuation and resource allocation are formulated by a multi-stage stochastic programming problem. From the view of systematic optimization, the work considers the abruptness of flood disasters, the vulnerability of roads, and the uncertainty and urgency of emergency demand, a multi-objective collaborative model with dynamic flood scenarios and multiple periods is constructed. It not only implies the interval number, the triangle fuzzy number and Bayes's theorem, but also involves the quantity and capacity limitation of tools, the constraints of network accessibility and service continuity.
Based on the fuzzy number clarification method and model solution, multiple objectives of the model, such as the survival probability, timeliness, and fairness of humanitarian rescue, have been resolved and balanced throughout computational results. Humanitarian equity for allocating medical supplies in the fourth period under the medium and heavy flood is about 100%, which has the same as the value of daily medical supplies within the first and third period in the heavy scenarios. But the severity of the flood and the fairness of resources is not a simple cause-and-effect relationship. Moreover, our study reveals that the p ϕ t parameter, and the decision preference usually affects the model's performance of the model. Although the shortage of humanitarian service time is helped to improve the survival rate of victims, it is pointed out that time is not the only factor to be considered. The reason is that the single consideration of sub-objectives is one-sided. The "degree" of humanitarian choice, which shall be grasped, plays a mutually reinforcing role. It provides a reference for planners to prioritize the high-risk nodes and improve emergency-response capacity of flood-stricken areas.
The next step is to quantify the risk of flood disaster, explore the emergency location and reservation problems in the first stage, and further analyze humanitarian service for the victims and supplies in the next stage. Secondly, in an effect to build a network model with dynamic changes for hurricanes, earthquakes, snow disasters and other emergencies. It can explore the economic measurement of the psychological trauma or shadow of victims by the rescue performance. And a more comprehensive analysis will be conducted to test the objectives in more instances.

Data Availability Statement:
The data that support the findings of this study is available upon request to the corresponding author.

Conflicts of Interest:
The author declares that there is no conflict of interest regarding the publication of this paper.