Multiperiod Equitable and Effi cient Allocation Strategy of Emergency Resources Under Uncertainty

Yanyan Wang 1 · Baiqing Sun 2

Abstract Equitable and effi cient allocation of emergency resources is critical to ensure the success of relief eff orts.The challenge comes largely from two aspects: the resources available for allocation are usually limited in quantity, especially in the early period of emergency response; and a large amount of uncertain information in the relief process aff ects the decision making of resource allocation. A multiperiod allocation model of emergency resources that takes into account both effi ciency and equity based on uncertain disaster information is proposed. Interval number and triangular fuzzy number are introduced to describe the diff erent sources of uncertainty (for example, demand, transportation time, and maximum transport amount), and the loss caused by unmet demand is used to quantify equity. Then, the deterministic transformation method of uncertain parameters is designed and the linear weighted sum method is applied to solve the proposed model. Finally, a computational case based on the 2017 Jiuzhaigou earthquake in Sichuan Province, China was conducted to validate the proposed model.The results show that the proposed model is feasible in the multiperiod allocation of emergency resources among multidisaster sites, and the f indings can help emergency managers to allocate emergency resources more scientif ically, equitably, and eff ectively under uncertainty.

Keywords China · Emergency resource allocation ·Equity and effi ciency · Uncertainty · Unmet demand

1 Introduction

In recent years, large-scale disasters such as earthquakes,tsunamis, volcanic eruptions, and f loods have occurred frequently in many regions of the world (Eshghi and Larson 2008; Adeagbo et al. 2016; Ogie and Pradhan 2019; Wang 2021), causing huge loss of life, property, and the economy(Barbarosoğlu and Arda 2004; Ahmadi et al. 2022). Especially in China, large-scale earthquake disasters such as the 1976 Tangshan earthquake and the 2008 Wenchuan earthquake caused catastrophic damage and impacts. According to the statistics from the China Earthquake Administration( 2022), from 2001 to 2021, 9353 earthquakes occurred in China, including 138 earthquakes with a magnitude of 6.0 or above (Fig. 1). In situations such as these, governments and researchers devote increasing attention to disaster relief to mitigate losses and impacts (Rising and Lurie 2013; Su et al. 2016).

In post-disaster emergency rescue work, material relief is the key to meet the survival needs of victims and the recovery of disaster-aff ected locations (Hu, Wang, et al. 2016;Wang et al. 2019). In addition to the effi cient allocation of resources, equitable allocation of available resources is also crucial in emergency relief (Erbeyoglu and Bilge 2020).Ignoring equity in emergency resource allocation leads to discontentment of victims and social unrest, increasing the burden of additional relief work and negative social impact(Ahmadi et al. 2022). After disasters, the key to eff ectively reducing the casualties and losses in the disaster area lies in the equitable and effi cient allocation of emergency resources(Kovacs and Spens 2009; Wang et al. 2019).

Fig. 1 Frequency of earthquakes in China from 2001 to 2021 (magnitude ≥ 6.0)

However, there are some challenging problems in the process of emergency resource allocation. First, the initial supply of emergency resources is limited, and there may be conf licts between diff erent rescue objectives; and second,information about resource demand, transportation time, and the maximum total amount of resources that can be transported from supply sites to aff ected locations is uncertain and inaccurate. Therefore, how to ensure the multiperiod equitable and effi cient allocation of emergency resources among multiple disaster-aff ected locations under uncertain information conditions is an important practical problem to be solved urgently in the f ield of emergency relief operations.

This study proposes a multiperiod emergency resource allocation model that balances equity and effi ciency under uncertainty. The objective of the proposed model is to minimize the total time of allocating emergency resources and the total losses caused by unmet demand. The purpose is to obtain an optimal multiperiod emergency resource allocation scheme by balancing the two objectives of equity and effi ciency, so as to provide support for emergency rescue decision making under uncertainty. The main contributions of this work are threefold: (1) the development of a multiperiod emergency resource allocation model under uncertainty that can maximize equity as well as effi ciency, and explore trade-off s between these two decision criteria; (2)the introduction of the interval number and triangular fuzzy number to describe diff erent sources of uncertainty, and the design of a deterministic transformation method to deal with the uncertainty; and (3) the illustration of the advantages of the proposed model in achieving equitable and effi cient resource allocation under uncertainty with a multi-angle case analysis.

Section 2 reviews the literature relevant to our work,Sect. 3 describes the research problem and proposes the notations and model formulation, Sect. 4 presents the solution method to solve the proposed model. A case study of the 2017 Jiuzhaigou earthquake in Sichuan Province, China is used to illustrate the applicability and potential advantages of the proposed model in Sect. 5. Finally, the conclusion is stated in Sect. 6.

2 Literature Review

Emergency resource allocation in humanitarian logistics has attracted increasingly more attention in recent research due to the increased frequency and destructive eff ects of disasters (Sheu and Pan 2014; Kamran et al. 2016; Moreno et al. 2016; Farahani et al. 2020). Constructing and solving models to provide decision support for emergency resource allocation has become the focus of existing research (Tzeng et al. 2007; Mete and Zabinsky 2010; Anaya-Arenas et al.2014; Özdamar and Ertem 2015; Wang 2021).

As the primary mission of humanitarian relief operations is saving lives, alleviating suff ering, and reducing property losses, the majority of the studies focus on effi ciency or eff ectiveness objectives/criteria (Hoyos et al. 2015; Wang 2021), that is, the shortest total allocation time (Wex et al.2014; Luan et al. 2020; Tang and Ye 2021), the lowest total transportation cost (Barbarosoğlu and Arda 2004; Minas et al. 2015; Zhou and Erdogan 2019; Zahedi et al. 2020;Han et al. 2021), the minimum total casualties (Fiedrich et al. 2000), and the shortest total distribution path (Sigal and Pritsker 1980; Yuan and Wang 2009; Hu et al. 2020).Tang and Ye ( 2021), for example, developed an emergency medical resource allocation model for multiple distribution centers and multiple demand points with the objective of minimizing the total allocation time of loading, unloading,and transportation. Han et al. ( 2021) constructed an optimization model of emergency supplies allocation for postdisaster relief operations with time window constraints to minimize total transportation costs. Fiedrich et al. ( 2000)proposed an optimization model for allocating emergency resources to minimize the total number of casualties during the search and rescue period. Hu et al. ( 2020) established a dynamic allocation model of emergency medical supplies considering the urgency of demand points with the goal of minimizing vehicle distribution paths.

Another criterion, equity, has been gradually paid attention to in recent years due to the need of realistic emergency relief decision making (Rennemo et al. 2014; Wang and Sun 2020). Equity and effi ciency are two important but conf licting goals in disaster relief operations (Feng et al.2017). Fast distribution may lead to inaccurate allocation of resources, while the allocation scheme with the lowest transportation cost will lead to unfair allocation of resources among multiple disaster sites (Wang et al. 2019). In disaster relief operations, the most effi cient allocation scheme may not be the most equitable, and inequities in the allocation of emergency resources may lead to increased social costs and negative impacts (Hu, Liu, et al. 2016). Especially in the case of resource shortage, equitable allocation is of great signif icance for improving the overall eff ects of emergency relief operations (Wang 2021). In the equity models, equity criteria are usually measured by factors such as demand satisfaction, deprivation cost, or deprivation time (Huang and Raf iei 2019). Holguín-Veras et al. ( 2013) introduced deprivation cost into an emergency resource allocation model for post-disaster relief, which mainly quantif ies equity through a series of negative impacts caused by a lack of resources. Huang and Raf iei ( 2019) used deprivation time–which ref lects the time of delay between sequential deliveries–as an equity criterion for allocation of emergency supplies. Chen and Fu ( 2018) developed a model for allocation of the f irst batch of emergency resources at the initial stage of large-scale disasters, taking the minimum total weighted envy value as the equity target and proportional equity as the constraint. He et al. ( 2021) took minimizing the comprehensive jealousy value as the equity objective, and then proposed an emergency resource allocation model with multiple rescue centers and multiple epidemic aff ected areas.

The allocation of emergency resources is usually a multiobjective optimization problem, and some researchers have gradually begun to consider both equity and effi ciency in recent research (Fu and Chen 2018; Wang et al. 2019; Wang 2021). Fu and Chen ( 2018) developed a bi-objective model for the allocation of the f irst batch of post-disaster emergency resources, in which effi ciency is measured by the latest delivery time and equity is measured by the demand satisfaction rate of the aff ected point. Wang ( 2021) designed an optimal allocation model of emergency resources for crossregional collaborative relief with the goal of minimizing total delivery time and total distribution costs (effi ciency)while maximizing coverage of resources (equity). However,these studies all centered on deterministic emergency rescue information conditions, ignoring the impact of the uncertainty of the information on emergency material distribution in disaster relief situations.

Uncertainty exists widely in post-disaster emergency rescue operations, which has an important impact on the equitable and effi cient allocation of emergency resources (Hu,Liu, et al. 2016). Existing studies mainly focus on the uncertainty of emergency demand (Beraldi et al. 2004; Zhang et al. 2010; Wang et al. 2014; Qin et al. 2017 ), and some studies also focus on the uncertainty of transportation time caused by road network damage (Lai et al. 2020; Wang and Sun 2020). Zhang et al. ( 2010) introduced interval numbers to represent the material demand of disaster aff ected locations and proposed an emergency resource allocation optimization model with the goal of minimizing the total costs.Wang et al. ( 2014) analyzed the problem of relief resource allocation under fuzzy demand conditions by constructing an optimization model with the goal of minimizing the total costs and transportation time. Qin et al. ( 2017) aimed at minimizing the total costs and proposed an optimization model of vehicle routing for emergency relief resource allocation,considering insuffi cient supplies and uncertain demand. Lai et al. ( 2020) applied the discrete scenario set to describe the uncertainty of transport time and established an optimization model of emergency resource allocation with effi ciency goal priority.

These studies have provided a sound theoretical basis for our work. The brief review of related research work indicates that it is crucial to consider effi ciency, equity, and uncertainty in emergency resource allocation. However, since the ultimate goal of disaster relief is to minimize casualties, the effi ciency criterion is often given the f irst attention by previous studies. Although some studies began to consider both effi ciency and equity, most of them lack detailed quantif ication and analysis of the sources of uncertainty in disaster information. In addition, previous emergency resource allocation problems in uncertain situations usually focused on the uncertainty of demand or time, which makes it easy to ignore other factors related to emergency rescue operations. Thus, we attempted in this study to develop a multiperiod and multi-objective optimization model for emergency resource allocation that simultaneously captures equity and effi ciency under uncertainty. In contrast to most existing works, three diff erent sources of uncertainty in emergency situations–demand, transportation time, and maximum transport amount–are simultaneously considered and quantif ied in this article. We also explored methods of managing uncertainty in an emergency context to help obtain highquality or optimal resource allocation schemes. We verif ied the validity and feasibility of the proposed model through a case study of the Jiuzhaigou earthquake, aiming to provide decision support for the multiperiod equitable and effi cient allocation of emergency resources in large-scale disaster relief operations.

3 Mathematical Model

This section describes the research problem, explains the notations for model formulation, and proposes the multiperiod equitable and effi cient allocation model of emergency resources.

3.1 Problem Description

Large-scale sudden disasters usually have the characteristics of large impact scope, strong destructiveness, and high rescue urgency. Emergency decision makers are unable to accurately obtain the required relief information because the relevant information–such as material demand, delivery time, and transportation volume–can only be determined in a certain range through limited information. To meet the reality of disaster relief, interval number is introduced to describe the uncertainty of emergency resource demand and transportation time, and triangular fuzzy number is used to represent the maximum total transport amount of resources per period under the condition of limited road network capacity. The equity objective is quantif ied by the loss caused by unmet demand and the effi ciency objective is measured by the allocation time. Therefore, the research problem of this study is the multiperiod equitable and effi -cient optimal allocation of large-scale disaster emergency relief resources under uncertainty (Fig. 2). Decision makers should consider both equity and effi ciency of resource allocation at each period to minimize the total losses caused by unmet demand for all disaster-aff ected locations as soon as possible within the shortest time.

3.2 Notations

The sets, parameters, and variables used in the model formulation are as follows.

3.2.1 Sets

Fig. 2 Research framework of multiperiod equitable and effi cient allocation of emergency resources under uncertainty

3.3 Model Formulation

The objective functions and constraints of the proposed model are as follows.

3.3.1 Objective Function

The objective function (1) minimizes the total losses caused by unmet demand at each aff ected location in all emergency periods to pursue the equity criterion. The objective function (2) minimizes the total time of resource allocation in all emergency periods to pursue the effi ciency criterion.

4.1 Deterministic Transformation Processing of the Uncertain Parameters

In this study, interval number and triangular fuzzy number in the corresponding objective functions and constraints are deterministically transformed according to the following def initions.

4.1.1 Deterministic Transformation of Interval Number

Constraint ( 3) guarantees that the amount of allocated resource to the aff ected location in a given period cannot exceed the sum of the new demand and the shortage remaining from previous periods. Constraint ( 4) ensures that the amount of allocated resource from the supply site in a given period cannot exceed the sum of the new supply and the inventory remaining from previous periods. Constraint ( 5)ensures that the resource allocation at each aff ected location is not higher than the preset maximum unmet rate. Constraint ( 6) gives an expression for measuring the loss caused by unmet demand at the aff ected point. Constraint ( 7) shows the transportation capacity limit under the condition of road network damage. Constraint ( 8) indicates the demand for each aff ected location should be met as much as possible in a given period. Constraint ( 9) def ines the logical relationship between the decision variables, indicating whether to allocate resources. Constraints ( 10) and ( 11), respectively,denote the dynamic changes in the supply and demand of emergency resources. Constraint ( 12) gives an expression for measuring the unmet demand. Constraint ( 13) requires that the binary variable must be equal to either 0 or 1. Constraint ( 14) gives the nonnegative constraints of the decision variables.

4 Model Solution Method

We f irst transform the uncertain parameters in the objective function and constraints of the proposed model into deterministic parameters, and then use the linear weighted sum method to solve the model.

4.2 Solution Step Design

The proposed model is a multi-objective model, which needs to be solved by using multi-objective optimization theory. The commonly used methods are: linear weighted sum method, ideal point method, minimax method, and hierarchical sorting method (Yang and Zhang 2008). In this study, the linear weighted sum method is used to solve the proposed model. The specif ic solution steps are as follows.

Step 1. Solve the extreme value of each objective function.importance ranking index. The importance binary contrast consistency scale matrix composed of the factor set is:

Then, the sum of each row in the consistency scale matrix is arranged in order from large to small, and the importance order of the factor set under the condition of ranking consistency is given.

(2) Quantitative determination of the weight of factor set.According to the consistency scale matrixZof importance ranking, mark the sequence numbersi=1,2,…,rin the order of importance from large to small, make a binary comparison on the importance of the qualitative ranking factor set, combined with the relationship between tone operator and quantitative scale, and then obtain the ordered binary comparison matrix.

5 Computational Case

In this section, a case study of the 2017 Jiuzhaigou earthquake is used to test the eff ectiveness and feasibility of the proposed model.

5.1 The Case of the 2017 Jiuzhaigou Earthquake in Sichuan Province, China

A magnitude 7.0 earthquake occurred at 21:19:46 p.m.,China Standard Time, on 8 August 2017 in Jiuzhaigou County, Aba Prefecture, Sichuan Province, China( 33.20◦N,103.82◦E). According to Chinese state offi cials,as of 20:00 on 13 August 2017, the earthquake caused 25 verif ied deaths, 525 people were injured, 6 people were missing, 176,492 people were aff ected (including tourists),and 73,671 houses were damaged to varying degrees (among which 76 collapsed). The basic infrastructure systems such as transportation, telecommunication, and power and water supply had also been damaged. In such emergency situations, decision makers face two main challenges: (1) the immediate need for a large number of various emergency resources at each disaster-aff ected location; and (2) some relief information from the aff ected locations is uncertain in terms of transportation time, demand for diff erent resources,and maximum total transport amount due to the damagedinfrastructure systems. Thus, in this context, equitable and effi cient allocation of emergency resources is crucial to reduce the destructive impact and improve the relief eff ect after the disaster.

Table 1 New supply of emergency resources at the supply sites at the beginning of each time period (7 days)

Table 2 New demand for emergency resources at the aff ected locations at the beginning of each time period (7 days)

Table 4 Disaster coeffi cient of aff ected locations in each time period (7 days)

In this case, we considered allocating two types of emergency resources–tents (m1) and drinking water (m2)–from two supply sites, Chengdu (CD) and Ziyang (ZY), to f ive severely affected locations–Jiuzhaigou County (JZG),Ruoergai County (REG), Hongyuan County (HY), Songpan County (SP), and Pingwu County (PW). We took seven days as the time period to analyze the resource allocation in the 28 days (a total of 4 emergency periods, nearly a month) after the disaster. The reason is that the four periods can usually ref lect the entire process of earthquake disaster relief, including the initial stage of emergency response(supply is less than demand), the middle stage (supply continues to increase), and the later stage (supply exceeds demand). The relevant data for the computational case were chosen using a combination of real and hypothetical data,since some disaster data were not published or could not be obtained through offi cial reports. For example, the new supply of emergency resources at the supply sites at the beginning of each period (Table 1) was estimated according to the resource storage before the disaster and the mobilization ability during the disaster. The new demand for emergency resources at the aff ected locations at the beginning of each period (Table 2) was estimated based on the number of victims and the severity of the disaster. Based on interviews conducted in September 2017 with f ive staffmembers of the volunteer service center who participated in the entire process of emergency resource allocation, the maximum total amount of emergency resources transported from the supply sites to the aff ected locations in each time period (Table 3) was estimated. The disaster coeffi cient ofaff ected locations in each time period (Table 4) was measured by disaster intensity and vulnerability of victims. The transportation time of allocating resources from the supply sites to the aff ected locations in each period (Table 5) was obtained via Baidu Maps. The loading and unloading times per unit of emergency resources (Table 6) were estimated based on the same f ive interviews. All types of emergency resources are in units of thousands. The weights of the two objective functions are calculated according to the fuzzy weighting method (see Sect. 4.2). In the qualitative ranking of importance, we set the two objective functions to have the same degree of importance. Then through calculation, we derived the weights of these two objective functions, which areω1=ω2=0.5(this also represents the consideration of both the effi ciency criterion and the equity criterion). The optimization levels of the deterministic transformation of uncertain parameters are set asα=0.9,β=0.9 , andγ=0.95 ,respectively, and the maximum unmet rate of resource allocation is set here asрdmn=40% . We solved the computational case in Matlab R2016a on a computer with an Intel(R)Core(TM) 1.90 GHz processor with 16.0 GB of RAM.

Table 3 Maximum total amount of emergency resources that can be transported from the supply sites to the aff ected locations in each time period(7 days)

Table 5 Transportation time(hour) of allocating resources from the supply sites to the aff ected locations in each time period (7 days)

Table 6 Loading time and unloading time (hour) per unit of emergency resources

Table 7 Total emergency resource allocation results, total losses, and total time

Fig. 3 Network structure diagram of overall emergency resource allocation. CD Chengdu; ZY Ziyang; JZG Jiuzhaigou County; REG Ruoergai County; HY Hongyuan County; SP Songpan County; PW Pingwu County. The data format in this f igure is ( m1 , m2 ), which respectively represents the overall allocation results of resource m1 and resource m2 from supply sites to aff ected locations in the entire emergency period. The value unit of the numbers is one thousand tents and one thousand boxes of drinking water.

Fig. 4 Relationship between disaster coeffi cient, loss, and resource allocation time. JZG Jiuzhaigou County; REG Ruoergai County; HY Hongyuan County; SP Songpan County; PW Pingwu County.

5.2 Computational Results

This section analyzes the following three computational results to illustrate the eff ectiveness and potential advantages of the proposed model.

5.2.1 Overall Emergency Resource Allocation SchemeTable 7 shows the overall allocation scheme of emergency resources from supply sites to aff ected locations, and Fig. 3 shows the allocation network structure for emergency resources, which is an optimal strategy for balancing effi -ciency and equity. It comprehensively considers multiple factors, including the disaster coeffi cient of aff ected locations, the maximum total amount of resources that can be transported under the condition of limited road network capacity, the transportation time, and the loading time and unloading time for allocating resources. Among them, we found that the loss from unmet demand of resources in Pingwu County (PW) is zero, which is mainly due to the fact that the resource needs of this aff ected location were fully met throughout the emergency rescue operations. In addition, Fig. 4 shows that, on the whole, the higher the disaster coeffi cient of the aff ected location, the greater the total loss caused by unmet demand, and the longer the required resource allocation time.

5.2.2 Emergency Resource Allocation Scheme in Each Period and Corresponding Loss and TimeThe allocation results, losses caused by unmet demand, and time associated with the allocation of all types of emergency resources in each time period are shown in Table 8. Figure 5 presents the emergency resource allocation network structure for each period. The trends of the losses and allocation time in each time period are shown in Fig. 6. Figure 7 shows the shortfalls and satisfaction rates of required resources at the end of each time period.

Figure 5 shows that the model allocates emergency resources to all disaster-aff ected locations in every emergency period as long as the locations have resource needs,so as to ensure that their needs are met to the greatest extent,thus preventing each disaster-aff ected location from having large losses due to unmet demand.

Table 8 and Fig. 6 show that in the entire emergency period (28 days from period 1 to period 4), the loss caused by unmet demand shows a trend of gradual decline, and the allocation time presents a trend that f irst increases and then decreases. This is consistent with the actual resource allocation in emergency situations. The required resourcesare often insuffi cient in the initial period of emergency relief operations, resulting in greater loss caused by unmet demand, and the allocation time is relatively short due to the small amount of resources that can be allocated initially.With the development of rescue activities, the supply of resources gradually increases. In order to make up for the shortage of resources in the early time period, the allocation of resources in the middle time period of the emergency response will be greatly increased, so that the loss is gradually reduced and the allocation time is gradually increased.In the later period of the emergency response operations, as the demands of the aff ected locations are gradually met, the loss is gradually reduced to zero, and the allocation time is also reduced due to the reduction of the amount of resources that need to be allocated.

Table 8 Emergency resource allocation results, loss, and time in each time period (7 days)

Fig. 5 Network structure diagram of emergency resource allocation in each time period (7 days). CD Chengdu; ZY Ziyang; JZG Jiuzhaigou County; REG Ruoergai County; HY Hongyuan County; SP Songpan County; PW Pingwu County. The data format in this f igure is ( m1 , m2 ), which respectively represents the allocation results of resource m1 and resource m2 from supply sites to aff ected locations in each time period (7 days). The value unit of the numbers is one thousand tents and one thousand boxes of drinking water.

Fig. 6 Allocation time and the loss caused by unmet demand in each time period (7 days)

Table 9 Emergency resource allocation results at each aff ected location in each time period (7 days)

Fig. 7 Shortfalls and satisfaction rates of resources at the end of each time period (7 days)

Fig. 8 Satisfaction rate of emergency resources at each aff ected location in each time period (7 days). Note: JZG Jiuzhaigou County; REG Ruoergai County; HY Hongyuan County; SP Songpan County; PW Pingwu County

Fig. 9 Trends in loss caused by unmet demand at each aff ected location in each time period (7 days). Note: JZG Jiuzhaigou County; REG Ruoergai County; HY Hongyuan County; SP Songpan County; PW Pingwu County

Fig. 10 Trends in the time (hour) required for allocating resources to each aff ected location in each time period (7 days). Note: JZG Jiuzhaigou County; REG Ruoergai County; HY Hongyuan County; SP Songpan County; PW Pingwu County

Figure 7 shows a large degree of resource shortfall at the disaster-aff ected locations in the f irst period, but through the resource allocation operations, the shortfall in the subsequent periods gradually decreases and the satisfaction rate gradually increases. When there is no shortfall of resources,the satisfaction rate reaches 100% (that is, the satisfaction rate of resourcesm2andm1reaches 100% in period 3 and period 4, respectively). This shows that the proposed model is eff ective for multiperiod scientif ic allocation of emergency resources under the condition of uncertain disaster information.

5.2.3 Loss, Time, and Satisfaction Rate of Resource Allocation at Each Aff ected LocationThe specific quantity of emergency resources allocated at each disaster-aff ected location in each time period and the corresponding loss and time are shown in Table 9. The resource satisfaction rate at each aff ected location in each time period is shown in Fig. 8. Figures 9 and 10 show the trends of the loss caused by unmet demand and the time for allocating resources to each aff ected location in each time period, respectively.

Figure 8 shows that the emergency resource satisfaction rate of all aff ected locations in each time period shows a gradual upward trend. Even in the initial period of emergency rescue (period 1), the resource satisfaction rate of all aff ected locations remains at 60% or above, and eventually the satisfaction rate of all aff ected locations reaches 100%.This shows that the proposed model is eff ective in improving the resource satisfaction rate of disaster-aff ected locations.

From Figs. 9 and 10, we can see that through continuous resource allocation, the loss caused by unmet demand at each aff ected location in each time period is gradually reduced, reaching zero at the end of the fourth period, which indicates that the resource shortfall of each disaster-aff ected location is gradually reduced, and the resource demand of all disaster-aff ected locations is f inally fully satisf ied. At the same time, the general trend in the time required for allocating resources is gradually increased. The reason is that in order to reduce the loss caused by unmet demand as far as possible, supply sites often need to allocate resources to disaster-aff ected locations that are far away, which means that,to a certain extent, the reduction in loss caused by unmet demand comes at the cost of increased allocation time.

Among the disaster-aff ected locations, the site with the least loss and the shortest allocation time is Pingwu County (PW). In each period of emergency relief operation, the resource demand of PW is minimal among the f ive aff ected locations, the transportation time between PW and the supply sites is the least, and the unloading time is also the shortest. So the supply sites–on the premise of guaranteeing equity in allocation in the shortest time–give priority to allocating all the required resources to PW. This is crucial for reducing the total loss and allocation time, and also indicates that the model can take into account both the effi ciency and equity of resource allocation in each time period, and minimize the loss caused by unmet demand of all aff ected locations in the shortest time.

6 Conclusion

From the computational results and analysis, four major conclusions and insights can be drawn:

(1) The reduction of loss usually comes at the expense of the increase of allocation time, and the optimal resource allocation strategy requires a scientif ic tradeoff between these two criteria.

(2) The disaster coeffi cient aff ects the allocation eff ect of emergency resources. The higher the disaster coeffi -cient of the aff ected location, the greater the total loss caused by unmet demand, and the longer the required allocation time.

(3) On the premise of balancing effi ciency and equity, the proposed model can maximize the satisfaction rate of resource allocation at each disaster-aff ected location,even at the initial period of emergency rescue.

(4) The model considers uncertain disaster information such as demand, transportation time, and maximum transport amount in resource allocation, which is more in line with the actual emergency rescue operation situation.

One recommendation for further research is to extend the proposed model to consider more realistic uncertain factors(for example, damage degree of road network and the impact of secondary disasters). With the increasing complexity of research problems, more effi cient model solving methods should be designed. Future research of these issues will be benef icial for optimal multiperiod allocation of large-scale emergency resources in real world situations.

AcknowledgmentsThis work was supported by the National Social Science Foundation of China (Grant No. 21CGL044).

Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source,provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.