1 Background

Globally, wildfire accidents, particularly those caused by humans, have become more common in recent years (Plucinski et al. 2012; Page et al. 2014; Krawchuk and Moritz 2014). For example, from 2000 to 2013, there were 115,466 wildfires in China, the majority of which were due to human activities. Annually, there were approximately 8248 wildfires, with over 106,127 hectares of burned area and over 117 casualties on average (China Statistical Yearbook 2014). Wildfires can result in serious losses of forest resources and personal property. To effectively reduce and mitigate such losses, proactively developing suppression planning and emergency logistics responses within the context of Disaster Operations Management (DOM) is necessary. Measures that can drastically reduce the associated social, economic and environmental impacts caused by wildfires are especially crucial.

DOM, which was first introduced by Altay and Green III in 2006 (Altay and Green 2006), consists of the techniques for preparing a community and reducing the severity of damages caused by all disasters by developing communications systems, stockpiling inventory, building adequate structures, etc. (Hoyos et al. 2015). If executed properly, these techniques can make a community more resilient to natural disasters, (Guha-Sapir and Santos 2012). According to FEMA (2004), the DOM life cycle can be divided into four major stages: mitigation, preparedness, response and recovery. This four-phase life cycle provides a more focused view of emergency operations and management actions. Emergency responses involve the employment of resources and emergency procedures as guided by plans to preserve life, property, the environment, and the social, economic, and political structure of the community. The emphasis in emergency response operations is primarily placed on relief distribution, facility location and casualty transportation. Related emergency logistics planning includes the optimal pick-up and delivery schedules for vehicles within a considered time window and the optimal quantities and load types picked up and delivered on these routes. In terms of the emergency logistics planning for wildfires suppression, it includes the storage, transportation and delivery of rescue resources and the allocation and management of equipment and materials (e.g., fire-fighting forces and fire-fighting equipment). All these activities aim to reduce the damage caused by wildfires and assist with fire disaster relief operations.

Nevertheless, the decision-making process for emergency logistics planning has never been straightforward. This process differs greatly from its counterparts in normal business logistics because it involves a high level of urgency and uncertainty in terms of the number of people affected and in need of attention (Christie and Levary 1998; Pedraza-Martinez and Wassenhove 2012). Policy makers and technicians request scientific models to explain the damage caused by disasters and establish future scenarios of disaster risk evolution conditions (Rodrigues et al. 2014). Therefore, employing suitable tools and techniques to model this stochasticity in the decision-making process for effective preparedness and response to disasters is essential. Regarding wildfires, the need for such tools has led to the development of several prediction models (Martínez et al. 2009; Thompson and Calkin 2011; Ager et al. 2014) that focus on explaining spatial–temporal patterns with regard to certain variables (physiographic, infrastructural, socio-economic and weather-related) relating to the ignition of wildfires. However, in the literature, few studies incorporating wildfire propagation and logistics planning for disaster relief, revealing a research gap concerning appropriate solutions for these types of logistic problems.

In the above mentioned context, this paper aims to develop a two-layer emergency logistics system with a single depot and multiple demand sites for wildfire suppression based on the predicted trend in fire. The work described herein considers wildfires in the Daxingan Mountains in China as real investigation cases and applies the theories and methodologies of emergency logistics management, forest science and operations research. The novelty of this work arises from the improvement it makes on a well-established propagation model (i.e., the Wangzhengfei model) by rendering this model capable of accurately predicting the fire propagation behavior in the Daxingan Mountains to minimize the impact of the dynamic characteristics of fire behavior on the distribution of suppression resources. It ranks the severity of the fire sites in terms of indices representing the potential burned areas to assign the emergency priorities accordingly. Moreover, this work presents a new method for solving the allocation problem and vehicle routing problem (VRP) for the distribution of suppression recourses by distinguishing lower and higher propagation speeds.

To achieve this aim, the paper is organized as follows. Section 2 describes the relevant studies in the literature. Section 3 presents the fire propagation models, whereas Sect. 4 investigates the emergency logistics planning models (i.e., VRP) based on the outcomes from Sect. 3. In Sect. 5, real-world cases based on historical data concerning wildfires in the Daxingan Mountains are analyzed to demonstrate and verify the models. Finally, Sect. 6 concludes the paper.

2 Literature review

Emergency logistics is distinct from general business logistics in the following four respects (Caunhye et al. 2012): (a) Relief services are both urgent and highly diversified. (b) Accurate real-time information regarding urgent relief demands is nearly impossible to obtain. (c) The benefits of emergency logistics operations are always discounted because of non-emergency reactions. Often, emergency logistics decisions are focused more on efficiency and less on cost. (d) Governments intervene the market behavior to jointly participate in emergency logistics services. Emergency logistics has also raised numerous challenging issues (Sheu 2007): (a) The definition of emergency logistics remains ambiguous. (b) The timeliness of relief supply and distribution is hardly controllable in the context of emergencies. (c) Resource management for emergency logistics remains challenging. (d) Accurate, real-time relief demand information is required but almost inaccessible. Therefore, the existing studies are classified into two categories: (1) humanitarian or DOM during the response stage, especially for emergency logistics planning (Caunhye et al. 2016; Gutjahr and Nolz 2016; Wang et al. 2014) and (2) disaster dynamic prediction and evaluation of the emergency level (Sheu 2010).

2.1 Humanitarian or disaster operations management

Regarding studies on humanitarian or disaster operations, Caunhye et al. (2012) conducted a review of optimization models and classified the findings according to the type of operations undertaken (facility location, relief distribution and casualty transportation) and the model types, decisions, objectives and constraints. The urgent relief distribution problems are generally formulated as multi-commodity multi-modal flow problems with time windows (Rathi et al. 1992; Haghani and Oh 1996; Lei et al. 2015; Jin et al. 2015). Fiedrich et al. (2000) developed a dynamic combinatorial optimization model to determine the optimal resource rescuing schedule in order to minimize the total deaths during a search-and-rescue period. Barbarosoğlu et al. (2002) proposed a bi-level hierarchical decomposition approach for helicopter planning during a disaster relief service operation. In this approach, the top level was used to solve the tactical decision problems involving managing the helicopter fleet, crew assignment and the number of tours undertaken by each helicopter. Özdamar et al. (2004) noted that an emergency logistics plan can be obtained by solving a dynamic time-dependent transportation problem with the objective of minimizing the unsatisfied demands of all commodities through the planning horizon. Oloruntoba et al. (2016) reviewed and analyzed three social science and management theories (i.e., internationalization, behavioral and organizational theories) as worthy of consideration by scholars and practitioners in humanitarian operations research. Finally, Rodríguez-Espíndola et al. (2017) introduced a disaster preparedness system based on a combination of multi-objective optimization and geographical information systems to aid multi-organizational decision making.

To address the uncertain and dynamic nature of disasters, Chang et al. (2007) developed a scenario-based, two-stage stochastic programming model for the planning of flood emergency logistics with uncertainty regarding demand locations and quantities. The first stage consisted of grouping the disaster areas and classifying their emergency levels while minimizing the expected shipping distances. The second stage consisted of a location-allocation model, which was used to select local rescue based to open, the quantities of rescue equipment to send to storehouses and the routes that emergency transportation should follow to minimize costs. Balcik and Beamon (2008) combined the problems of resource allocation and relief distribution by proposing a multi-period Mixed Integer Planning (MIP) model to determine the schedules of a fixed set of vehicles and an equitable allocation of resources by minimizing transportation and unsatisfied or late-satisfied demand costs and maximizing the benefits to aid recipients for local distribution sites using a rolling horizon to address uncertainties in supply and demand. Subsequently, Rawls and Turnquist (2012) developed a multi-period version of the same approach using a two-stage, dynamic and scenario-based stochastic programming model for the pre-positioning and allocation of facilities and the distribution of commodities after a disaster event with the objective of minimizing total costs. Chu and Wang (2012) presented a new approach to address event probability uncertainties and analyze the probability distribution for constructing probable fire scenarios, and Jacobson et al. (2012) studied the probability distributions of aspects such as the survival time and service time of patients to develop a model for assigning different priority levels to disaster victims. Wohlgemuth et al. (2012) developed a multi-stage mixed integer model under variable demand and transport conditions in the context of disaster relief response planning and logistics. Huang and Song (2016) discussed an emergency logistics distribution routing problem based on uncertainty theory in which some parameters lack historical data and are given by experts’ estimations. Alem et al. (2016) developed a new two-stage stochastic network flow model to decide how to rapidly supply humanitarian aid to disaster victims considering the uncertainty about the exact nature and magnitude of the disaster.

Furthermore, Liberatore et al. (2013) established five major parameters based on which the uncertainties involved in humanitarian logistics are approached: (1) demand, including the size of the affected population and/or the quantities of required relief goods; (2) demand locations and (3) affected areas, indicating those directly related to the demographics of each location and the impact of the disaster; (4) supply, including considerations of the quality and availability of products in a post-disaster scenario; and finally, (5) the transportation network, including all possible damage to the infrastructure and the congestion effect. Iudin et al. (2015) applied this new approach to develop a cellular automaton forest-fire model related to the percolation methodology. Their model combined the dynamic and static percolation problems and exhibited critical fluctuations under certain conditions. Rezaei-Malek et al. (2016) developed a multi-objective, two-stage stochastic, non-linear, and mixed integer mathematical model for relief pre-positioning in disaster management based on a new utility level of the delivered relief commodities.

In DOM, many models are NP-hard (Sheu 2007; Caunhye et al. 2012); therefore, the algorithms for solving them are intensively studied (e.g., Hoyos et al. 2015; Özdamar and Ertem 2015; Yuan and Wang 2009). For example, Wei and Kumar (2007) proposed an ant-colony optimization method to solve the multi-commodity and vehicle dispatch problems that arise in disaster relief scenarios. Subsequently, Zhang et al. (2012) proposed a local search heuristic with MIP for solving the multiple-disaster, multiple-response emergency allocation and disaster relief problems, which were associated with the stochastic occurrence of a secondary disaster.

2.2 Disaster dynamic prediction

During rescue operations, fire disasters may continue to develop. Therefore, forecasting the evolution of the disasters (e.g., the fire or flood propagation) and incorporating it into the rescue decisions is necessary (Altay and Labonte 2014; Altay and Pal 2014; Oloruntoba 2010).

Regarding disaster evolution forecasting, García et al. (1995) introduced a logit model to predict the number of fire days in the Whitecourt Forest of Alberta. Mandallaz and Ye (1997) presented a general statistical methodology for predicting forest-fire occurrence using Poisson models. Wu et al. (2015) identified three fire environment zones among the boreal forest landscapes of northeastern China using analytical methods to examine the spatial clustering characteristics of the environmental variables of climate, vegetation, topography, and human activity. This work demonstrated how a developed fire environment zone map could be used to guide forest-fire management and fire regime prediction. As these forecasts were made at a macro level, the associated results might be used to formulate preparedness plans. Sheu et al. (2005) developed a model for forecasting relief demand and for clustering demand sites according to their individual uncertainty and urgency using a fuzzy clustering technique to allocate the daily consumption of relief commodities in disaster scenarios. The objective of this work was to minimize the costs of traveling from relief centers based on real-time analysis. The authors presented a real-time, micro-level forecast; hence, the results might be used as a basis for scheduling the operations of rescue vehicles and the delivery of emergency goods in real time. To properly support humanitarian decision makers, Charles et al. (2016) proposed a tooled methodology based on the definition of aggregate scenarios to reliably forecast demand using past disaster data and future trends.

Tzeng et al. (2007) developed a fuzzy multi-objective linear programming model for the design of relief delivery systems with the objectives of minimizing the total costs and total travel time and maximizing the minimal satisfaction during the planning period using a method of predicting the commodity demand at each established site and an uncertainty analysis of the achievement of each objective. This model combined emergency logistics planning and disaster spread prediction, and its dynamic planning decisions were based on the forecasted demands. Hu and Sheng (2015) conducted a similar study based on the principle of disaster spread in resource networks. They proposed a property model of resource nodes that incorporated the values of the resources, the disastrous energy of each node, the disaster spread path, and the disaster spread characteristics. To study the disaster spread dynamics, a state diagram was used to construct a microscopic behavior model for the resource nodes by treating the resource nodes as agents. Then, decision models were gradually constructed to determine the optimal time for disaster rescue and emergency resource preparation. Buscarino et al. (2015) disclosed a model that formalized an innovative methodology for simulating forest-fire propagation. This model was based on a multilayer network structure, which facilitated the establishment of flexible definitions of the spatial properties of the medium and the dynamical laws governing fire propagation. The dynamical core of each node in the network was represented by a hyperbolic reaction-diffusion equation that incorporated the intrinsic characteristics of the ignition behavior of trees.

Among the studies on disaster spread prediction (in the second category), statistical models are widely used in the risk assessment stage (Powell et al. 2016). For risk analysis of the impact of natural disasters on large-scale, critical infrastructure systems, electric power distribution systems, transportation systems, etc., Guikema (2009) proposed a statistical learning methodology that incorporated a diverse set of methods designed for performing inference analysis on large, complex data sets,. Matellini et al. (2013) developed a Bayesian network model to investigate the evolution of fires within dwellings and to assess the probabilities of the associated consequences. Holicky and Schleich (2000) studied fire protection systems and used a Bayesian network to model a building fire from its onset to the structural collapse of the building. This network included elements for the detection and automatic suppression of fires and a node representing fire brigade intervention. A fire scenario is characterized in terms of multiple interrelated variables, such as the ignition site, the nature of the fire, the fuel configuration, the number and locations of the occupants, the characteristics of the structure, and the ventilation, among others (Matellini et al. 2012). Several other studies have analyzed the emergency levels in affected areas by evaluating the degrees of severity associated with the affected areas (Sheu 2007; Ji and Zhu 2012).

Despite the efforts in previous works, the above review of the current literature reveals that there are few studies on incorporating dynamic disaster evaluation models into risk-based emergence logistics. To address this research gap, the kernel of this paper includes integrating models for predicting the fire propagation rates at fire sites, ranking the emergency levels of and determining the emergency priorities of the fire sites, and finally allocating resources and routing vehicles for suppression operations. The paper presents a new method developed for fighting wildfires in a typical forest area in China in terms of emergency logistics. By conducting careful calculations based on real data from wildfires that occurred in the Daxingan Mountains on March 19th, 2003, and on June 29th, 2010, we verify the model and demonstrate its advantages by comparing the findings with real-world observations.

3 Propagation model for wildfires in the daxingan mountains

A propagation model is used to calculate the impact of factors such as combustible materials, weather conditions and terrain on forest-fire proliferation. Using a suitable propagation model enables the spread pattern of wildfires in the mountain area of interest to be estimated, and the suppression or daily protection effort to then be effectively implemented based on the forecasted results. The most studied propagation models for this purpose are those for surface fires because they are the most common type of wildfires in mountain areas. There are two well-known classes of these models: real forest-fire models, in which adjacent trees catch fire in a step-by-step manner, and simplified versions, which involve instantaneous combustion (Iudin et al. 2015).

The systematic propagation models are the Rothermel model, the McArthur model, the Canadian propagation model and the Wangzhengfei propagation model (Sullivan 2009). The Rothermel model simulates flame front propagation according to the principle of the conservation of energy. This model has been widely applied in mountain areas, where combustible materials are evenly distributed and the effects of large clusters of combustible materials can be ignored (Rothermel 1983). However, because the Rothermel model is semi-experimental, certain parameters have to be obtained through experiments, and considerable input data are needed. In most areas of China, this model cannot be used because of data unavailability (Weise and Biging 1997; Perry 1998).

The McArthur model, a mathematical description of McArthur’s gauge for measuring fire risk, was developed by Noble et al. (1980). This model can be used to forecast fire weather situations or to predict certain parameters of fire behavior. Because it can be used only in areas where combustible materials exist in discrete units (e.g., grassplots or eucalyptus plantations), it has typically been applied in areas with a Mediterranean climate and hence is not suitable for areas with other climates (e.g., the Daxingan Mountains in China) (Chai et al. 1988; Yue et al. 2007; Miao et al. 2012).

In the Canadian forest-fire propagation model, combustible materials are divided into five categories based on Canadian vegetation types: conifers, broad-leaved trees, mingled forest, cutting bases and open ground. It is a statistical model, and the propagation rate equation was determined by inspecting the attributes of 290 forest fires. Although the Canadian propagation model can clearly reproduce the entire evolution of a forest fire and forecast the flame behavior and fire development pattern, it does not consider the heat transfer mechanisms of the combustible materials. Because it lacks a physical basis, the model’s accuracy becomes questionable when the actual conditions differ substantially from those corresponding to the experimental data used to construct the model.

The Wangzhengfei propagation model was developed based on real data from wildfires in the Daxingan Mountains and is suitable for other areas with slopes of less than \(60^{^{\circ }}\) (Yang et al. 2011). Because the Daxingan Mountains exhibit an undulating topography with slopes of less than \(15^{^{\circ }}\) accounting for 80% of the area, this model is selected to provide basic insights for the updated fire propagation model in this work. The Wangzhengfei propagation model is described as follows.

  • Model assumptions

  1. A1:

    The initial fire propagation rate varies linearly with temperature and wind power;

  2. A2:

    The target area exhibits an undulating topography, with no slope greater than \(60^{^{\circ }}\);

  3. A3:

    The target area contains regions of grassy marshland, secondary forest and coniferous forest;

  4. A4:

    The weather data are available in a timely manner.

  • Model structure

    $$\begin{aligned}&\hbox {Fire propagation model:}\,V_{F} = V_{0}K_{S}K_{W}K_{\varphi } = V_{0}K_{S} K_{\varphi }e^{0.1782V_w } \end{aligned}$$
    (1)
    $$\begin{aligned}&\hbox {Initial propagation speed model:}\,V_{0} =aT+bW+c \end{aligned}$$
    (2)

Here, \(V_{F}\) is the flame spread rate (m/s), \(V_{W}\) is the wind speed (m/s), \(V_{0}\) is the initial flame spread rate (m/min), \(K_{S}\) is a coefficient representing the different types of combustible materials, W is the wind force (grade), \(K_{W}\) is the coefficient of the wind force, \(K_{\varphi }\) is the coefficient of the geographic slope, \(\varphi \) is the slope, T is the temperature \(({^{\circ }}C)\); and a, b, and c are parameters. Based on the research by Wang (1992), parameters a, b, and c in the Daxingan Mountains are usually equal to 0.053, 0.048, and 0.275, respectively.

  • Variables and parameters

(1):

Initial spread speed

This quantity represents the physical and mechanical attributes of the combustible materials and their degrees of wetness/dryness; its value is typically estimated by burning samples of combustible materials and by the present temperature and wind force.

(2):

Coefficients for different types of combustible materials

These quantities represent the degree of inflammability of each material and its suitability for burning. The main types of combustible materials present in the Daxingan Mountains are marshy grassland, secondary forest and coniferous forest. The corresponding coefficients are therefore \(K_{s\,\textit{marshy}}=1.0\), \(K_{s\,\textit{secondary}} =0.7\), and \(K_{s\,\textit{coniferous}}=0.4\) (Wang 1992).

(3):

Coefficient of the wind force

A higher wind speed results in more rapid spreading of flames following the wind, i.e., a higher wind-induced yield value \(K_{w}\). For a grassland region with zero slope, the \(V_{F}\) values corresponding to wind forces of grades 1–12 are listed in Table 1 (Wen and Liu 1994).

Table 1 Wind speeds and wind forces for a grassland region with zero slope
(4):

Coefficient of the geographic slope

This parameter is one of the most important factors affecting fire propagation. The impact of the slope on the propagation speed \(K_{\phi }\) is summarized in Table 2 (Wen and Liu 1994).

Table 2 \(K_{\phi }\) values corresponding to different slopes

4 Relief allocation and vehicle routing model

Relief allocation and delivery models consider the demand at each fire site and the transportation of supplies to these sites by allocated vehicles. The vehicles begin their itineraries at depots and may serve several wildfire sites along their own routes. During this phase, the demand quantities and location uncertainties are the major factors influencing rational decisions.

Because allocation and delivery decisions determine the suppression efficiency, optimizing these decisions together based on different fire scenarios is necessary. During the suppression operations, various supplies are needed at different fire sites, and a rescue depot may store many different types of emergency resources. There are two main types of emergency resources that are needed to combat wildfires: fire-fighting forces and fire-fighting equipment. Emergency logistics involving wildfires include the storage, transportation and delivery of rescue resources and the use and management of the equipment and materials. Among the activities, the most fundamental are the transportation and scheduling of the fire-fighting forces.

4.1 Problem statement

Rescue commanders optimize the allocation of suppression resources and the routing of vehicles based on the rate of fire propagation at each site (10 m/min is the critical rate) (Yao and Wen 2002). In the case of rapid propagation \((V_{F}>10 \hbox { m/min})\), the rescue commanders may arrange several vehicles to transport suppression resources to each fire site to minimize total losses (Fig. 1). When the fire propagation rate is slow \((V_{F}<=10 \hbox { m/min})\), the commanders may send resources to several fire sites via one vehicle (Fig. 2) (Zhang 2007; Sun et al. 2007). This work investigates the resource allocation and vehicle dispatch problems for both cases.

Fig. 1
figure 1

Rescue operation for fast propagation

Fig. 2
figure 2

Rescue operation for slow propagation

4.2 Optimization model for rapid fire propagation

Suppose that a system for the transportation of suppression recourses for combating a wildfire consists of a single depot S (with capacity s) and I fire sites (demand points with demands \(d_i\>\,( i=1,2,{\ldots },I)\) for each site).

For the allocation and transport of the resources, both efficiency and equity should be considered; i.e., resources should be delivered to the demand sites as soon as possible, and each demand site should receive them on a just-in-need basis. During the suppression period, the decision makers must maximize the total benefit from all fire sites. Therefore, they should rank the severities of the fire sites based on the consequences of the fire accidents and damage attributes to allocate and deliver resources accordingly. Consequently, fire propagation rate becomes the main attribute of concern, as sites where the propagation is more rapid will suffer more serious damage and thus should be given a higher priority.

4.2.1 Model assumptions

  1. A1:

    Flames spread evenly from each fire site, with a spread radius that is proportional to the burn time;

  2. A2:

    The quantity of suppression resources is represented in combined units of fire-fighting resources in a way that the number of firemen and the amount of materials in each unit are fixed;

  3. A3:

    The suppression efficiency of all firemen is the same;

  4. A4:

    The rescue vehicles are identical in terms of their type and travel speed;

  5. A5:

    The flame propagation speed can be estimated and pre-defined using the model;

  6. A6:

    The supplies/resources available from the depot are sufficient to satisfy the demands at all fire sites.

4.2.2 Variables and parameters

\(C_1 \) is the loss of burning a unit of forest \((\hbox {hm}^{2})\) in monetary terms, \(C_2 \) is the time cost of a combined suppression resource unit, \(C_3 \) is the cost of transporting one combined resource unit, \(X_i \) is the number of combined resource units allocated to fire point i, \(t_{i{ 1}} \) is the time at which suppression begins at fire site i, \(t_{i{ 2}} \) is the time at which the fire at fire site i is extinguished, \(V_{Fi} \) is the fire propagation speed at fire site i, \(B(t_i )\) is the burned area at fire site i in period t, \(C(t_i )\) is the loss due to fire at site i in period t, N is the number of firemen in a combined resource unit, \(V_{i} \) is the speed of suppression, \(V_{T} \) is the travel speed of an emergency vehicle, \(d_i \) is the distance from the depot to fire site i, \(A_i \) is the number of firemen allocated to site i, and \(X_i^*\) is the optimal number of combined resource units to allocate to fire site i.

4.2.3 Model formulation

Let the time of ignition at fire site i be \(t_{i{ 0}}=0\); then, the time at which the firemen arrive at the fire site is \(t_{i{ 1}}= d_{i}{/}V_{T}\). From A1, we know that the change in the burned area \(B(t_{i})\) is continuous to the first order. In other words, the rate of increase is \(dB(t_{i})/dt_{i}\). Generally, the burning speed will decrease once the firemen begin the suppression operations. At fire site i, after fire-fighting begins at time \(t_{i{ 1}}\), the propagation rate will decrease until the fire is extinguished at time \(t_{i{ 2}}\). The burned area is represented by the triangle in Fig. 3 and can be calculated as follows:

Fig. 3
figure 3

The area of the burned forest at a fire site

$$\begin{aligned} \frac{dB(t_i )}{dt_i }=\left\{ \begin{array}{lll} V_{Fi} t_i,\quad 0\le t_i \le t_{i1}, &{} {i=1,2,...,I} \\ V_{Fi} t_i -V_i X_i (t_i -t_{i1} ),&{} t_{i1} \le t_i \le t_{i2}, &{} {i=1,2,...,I} \end{array}\right. \end{aligned}$$
(3)

when \(\frac{dB(t_i )}{dt_i }| {_{t_i =t_{i2} } =0},\)\(t_{i2} =t_{i1} +V_{Fi} t_{i1} /(V_i X_i -V_{Fi} ), {i=1,2,...,I} ,\) the corresponding area can be calculated as follows: \(B(t_{i2} )=V_{Fi} t_{i1}^2 /2+V_{Fi}^2 t_{i1}^2 /2(V_i X_i -V_{Fi} ),\quad {i=1,2,...,I} \).

In suppression, there are two kinds of costs: (1) the disposable costs of transporting suppression materials and (2) the costs of the suppression tools and the salaries of the workers. Thus, the total loss from a wildfire can be calculated as follows:

$$\begin{aligned} C(X_i )=\left( \frac{V_{Fi} t_{i{ 1}}^2 }{2}+\frac{V_{Fi}^2 t_{i{ 1}}^2 }{2(V_i X_i -V_{Fi} )}\right) C_{1}+\frac{C_{{ 2}}V_{Fi}t_{i{ 1}}X_{i}}{V_{i}X_{i} -V_{Fi} }+C_{3}X_{i}N_{i} ,\quad {i=1,2,...,I}\qquad \end{aligned}$$
(4)

The value \(X_{i}^{*}\)corresponding to the extremum of Eq. (4) is the necessary number of vehicles.

$$\begin{aligned} X_{i}^{*}= & {} {\frac{V_{Fi} }{V_{i}}+\frac{1}{V_i }\sqrt{(C_{{ 1}}V_{Fi}^{2} V_{i}t_{i{ 1}}^{2} +2C_{{ 2}}V_{Fi}^{2} t_{i{ 1}}^{2} )/2C_{{ 3}}N},}\quad {i=1,2,...,I} \end{aligned}$$
(5)
$$\begin{aligned} A_{i}= & {} N_{i}X_{i}^{*} = N_{i}\left( {\frac{V_{Fi}}{V_{i}}}+{\frac{1}{V_{i}}}\sqrt{(C_{{ 1}} V_{Fi}^{2}V_{i}t_{i{ 1}}^{2}+2C_{{ 2}}V_{iF}^{2}t_{i{ 1}}^{2})/2C_{{ 3}}N}\right) ,\quad {i=1,2,...,I}\nonumber \\ \end{aligned}$$
(6)

In this suppression scheme, the earliest time that the fire can be extinguished is calculated as follows:

$$\begin{aligned} T= & {} Max\left\{ {t_{12}, t_{22},...,{t_{{\mathrm{k}2}}}}\right\} \nonumber \\= & {} MAX\left\{ {t_{11}+\frac{V_{F1}t_{11}}{V_{i}X_{1} -V_{F1}},t_{21}+\frac{V_{F2}t_{21}}{V_{i}X_{2} -V_{F2}},...,{t_{I1}+\frac{V_{FI}t_{I1}}{V_{i}X_{I} -V_{FI}}}}\right\} ,\nonumber \\&\qquad \qquad {i=1,2,...,I} \end{aligned}$$
(7)

4.3 Optimization model for slow fire propagation

Because suppression resources must be transported within a certain time window and because the transportation cost is directly proportional to the distance that a vehicle travels, the allocation of suppression resources is a VRP concerning the routing of vehicles from the rescue depot to the fire sites under the constraints placed by the vehicles’ loading capacities (Wang et al. 2008; Zhu et al. 2007)

4.3.1 Model assumptions

  1. A1:

    The depot possesses sufficient vehicles, and each vehicle’s loading capacity is fixed;

  2. A2:

    The demands at the fire sites can be estimated based on the fire propagation rate;

  3. A3:

    Each vehicle’s travel speed is fixed, and the transport time is equal to the vehicle’s travel time (loading/unloading times are not taken into account);

  4. A4:

    A fire site with a faster propagation rate is considered to be suffering more severe fire damage and should be given a higher priority.

4.3.2 Variables and parameters

Let S represent the depot. \(D=\{i|i=1,2,...,I\}\) is the set of fire sites; \(C=\{m|m=1,2,...,M\}\) is the set of emergency vehicles; \(A=S\cup D\) is the set of all nodes, including the depot and the fire sites; \(C^{m}\) is the fixed cost of vehicle m; \(C_{ij}^m\) is the unit time cost of vehicle m traveling from i to j; \(d_{ij}\) is the distance between i and j; \(Q^{m}\) is the approved loading capacity of vehicle m; \(V_{Fi}\) is the fire propagation speed at fire site i; \(V_1\)is the fire-fighting speed of the firemen; \(u_i\) is the demand at fire site i (in terms of combined resource units); \(u_i =R_i /V_{i1} \); \(t_i^m \) is the time at which vehicle m arrives at fire site i; \(t_i^L \) is the latest time for the arrival of suppression materials at fire site i; \(t_{ij}^m \) is the time required for vehicle m to travel from i to j; and \(V_{ij}^m \) is the average travel speed of vehicle m from i to j.

$$\begin{aligned} x^{m}= \left\{ \begin{array}{ll} 1&{} {\hbox {vehicle }m \hbox {is used}} \\ 0&{} {\hbox {otherwise}} \end{array}\right. ; Y_{ij}^m = \left\{ \begin{array}{ll} 1\;\;\hbox {vehicle}\;m\;\hbox {from}\;\hbox {site}\;i\;\hbox {to}\;j\;\left( {i\ne j}\right) \\ 0\;\;\hbox {otherwise} \end{array} \right. . \end{aligned}$$

4.3.3 Model formulation

$$\begin{aligned}&\displaystyle Min{:}\,F_1 = \sum \nolimits _{m\in C} {\sum \nolimits _{i\in D} {t_i^m } } \end{aligned}$$
(8)
$$\begin{aligned}&\displaystyle Min{:}\,F_2 = \sum \nolimits _{m\in C} {x^{m}C^{m}} +\sum \nolimits _{m\in C} {\sum \nolimits _{i\in A} {\sum \nolimits _{j\in A} {C_j^m } } \frac{d_{ij} }{V_{ij}^m }y_{ij}^m } \end{aligned}$$
(9)
$$\begin{aligned}&\displaystyle \sum \nolimits _{i\in A} {y_{ij}^m } -\sum \nolimits _{i\in A} {y_{hj}^m } = 0,\forall m\in C,\forall h\in I \end{aligned}$$
(10)
$$\begin{aligned}&\displaystyle \sum \nolimits _{j\in D} {\sum \nolimits _{m\in C} {y_{Sj}^m } } > 0 \end{aligned}$$
(11)
$$\begin{aligned}&\displaystyle \sum \nolimits _{i\in A} {\sum \nolimits _{j\in D} {y_{ij}^{m} u_{j} } } \le Q^{m},\forall m\in C \end{aligned}$$
(12)
$$\begin{aligned}&\displaystyle \sum _{\hbox {m}\in C} {\sum _{i\in A} {\sum _{j\in D} {y_{ij}^m u_j \le Q_s } } } \end{aligned}$$
(13)
$$\begin{aligned}&\displaystyle t_i^m = t_j^m +t_{ji}^m y_{ji}^m ,\forall i,j\in A,\forall m\in C \end{aligned}$$
(14)
$$\begin{aligned}&\displaystyle t_i^m \le t_i^L ,\forall i\in D,\forall m\in C \end{aligned}$$
(15)
$$\begin{aligned}&\displaystyle y_{ij}^m = \{0,1\},\forall i,j\in A,\forall m\in C \end{aligned}$$
(16)
$$\begin{aligned}&\displaystyle x^{m}= \{0, 1\}, \forall m\in C \end{aligned}$$
(17)

In a suppression operation, both suppression speed and cost should be considered; thus, when the fire is propagating slowly, the optimization model may have two objectives. Equation (8), which is the primary objective, minimizes the total travel time of the emergency vehicles. Equation (9), which is the secondary objective, minimizes the total cost of the emergency vehicles. The constraints are as follows.

Equation (10) means that each vehicle will depart immediately after unloading at a fire site. Equation (11) means that the rescue depot S possesses sufficient vehicles to dispatch. Equation (12) means that the demand at any fire sites should be less than the vehicle capacity. Equation (13) means that the sum of the demands at all fire sites should be less than the inventory of depot S. Equation (14) means that the travel time from depot S to fire site j is equal to the combination of the time of the vehicle’s arrival at site iwith the travel time from i to j. Equation (15) means that the time of a vehicle’s arrival at a fire site should be earlier than the latest time at which the site requires the resources.

5 Case study

5.1 Case of rapid propagation

  1. 1.

    Fire attributes

According to the database of the Wildfires Prevention Headquarters of the Heilongjiang province and the Scientific Data Share Center of the China Meteorological Administration, a large-scale wildfire occurred in the Nanweng River area in the Daxingan Mountains at 10:35 on March 19th 2003. There were 4 fire sites in total: Wuoduhe (WDH), Xiaogushan (XGS), Naduli (NDL) and 597-Highland (597H). The areas are mainly covered by marshy grassland, brushwood and forest. The total affected area was \(133,008\hbox { hm}^{2}\), consisting of \(240\hbox { hm}^{2}\) of burned forestland and \(132,768\hbox { hm}^{2}\) of barren mountain grassland. The suppression operation was undertaken by 385 professional firemen and 1936 forest rangers. The fire was extinguished at 13:50 on May 23rd, 2003. Some of the detailed information that was obtained is shown in Tables 3 and 4. Other parameters were evaluated based on the interviews with the workers in the Forestry Bureau of the Daxingan Mountains and were as follows: \(c_{1}=15\) Yuan, \(c_{2}=1.3\) Yuan, \(c_{3}=0.82\) Yuan, \(v_{2}=100\hbox { km/h}\), \(v_{1}=12.5\hbox { m/min}\), and \(N=3\) persons.

Table 3 Meteorological data at the fire sites
Table 4 Geographical data at the fire sites
  1. 2.

    Solution results

According to the models presented above [i.e., Eqs. (3, 4)], the propagation rates at the four fire sites (WDH, XGS, NDL and 597H) were \(V_{F1}=179.38\hbox { m/min}\), \(V_{F2}=318.94\hbox { m/min}\), \(V_{F3}=105.84\hbox { m/min}\), and \(V_{F4}=188.64\hbox { m/min}\), respectively. The calculated allocation results for the combined suppression resource units and firemen teams are listed in Table 5. Comparing Table 5 with Table 6 reveals that with the optimized allocation of the combined suppression resource units and firemen teams, the wildfires could have been extinguished in a much shorter period of time. The shortest time required for putting out the fire at a single site is found to be 19.3 h. Compared with the actual period of 2 months, the optimized resource allocation and transportation scheme is much more efficient, albeit higher in cost. In the optimized analysis, 27,903 firemen are needed, which is 10,000 more than the number of firemen who participated in the real case, and the number of combined fire-fighting resource units is 9,301, which is 4 times greater than the actual number used. However, the higher cost may be well justified by the actual loss of forest resources suffered.

Table 5 Allocation scheme for emergency supplies
Table 6 Comparison of real and optimized results

5.2 Case of slow propagation

  1. 1.

    Description of the wildfire

At 10:40 on June 29th 2010, a wildfire was caused by lightning strikes in the Huzhong zone in the region administered by the Huzhong Forest Administrative Bureau. There were 7 fire sites: 1231-Highland (1231H), H311 Line (H311), X59 Line (X59), H59 Line (H59), LWM12 Line (LWM12), LWM3 Line (LWM3) and T6 Line (T6). The area involving the fires was approximately \(138.5\hbox { hm}^{2}\). The suppression activities were undertaken by 500 professional firemen and 140 forest rangers. The wildfire was extinguished on July 3rd. The distances between fire sites and the meteorological and geographical data, which were obtained from the Forestry Bureau of the Daxingan Mountains, are provided in Tables 7, 8 and 9.

Table 7 Distances between fire sites
Table 8 Meteorological data at the fire sites
Table 9 Geographical data at the fire sites

The Huzhong Forest Administrative Bureau fire station used 3 vehicles to send suppression materials to the 7 fire sites. The capacity of each vehicle was 9 tons, and the travel speed was 100 km/h. The fixed cost and unit time cost of each vehicle were 150 Yuan and 100 Yuan/h, respectively. One combined resource unit consisted of two instruments, one water cannon, one water pump and three firemen. The weight of one combined resource unit was 1000 kg, and the suppression speed was 2.5 m/min.

  1. 2.

    Results of the allocation model

The fire propagation rates at the 7 fire sites were calculated using the above propagation model (i.e., Eqs. (3, 4)), and the results are shown in Table 10.

Table 10 Fire propagation rates at the seven fire sites

Based on these speeds, the rescue priorities were ranked. A site with a higher priority ranking requires suppression materials more urgently. Therefore, the emergency vehicles should proceed from sites with higher scores to sites with lower ones. With the depot coded as 0, the demands calculated according to the suppression speed and the fire propagation rate at the fire sites are shown in Table 11 in terms of combined resource units.

Table 11 Demands for emergency materials at the seven fire sites
  1. 3.

    Vehicle routes

The VRP model could be effectively solved using an immune clonal algorithm (ICA) (Ma et al. 2009). The ICA, which is derived from the traditional evolutionary algorithm, introduces the mechanisms of avidity maturation, cloning and memorization. The performance of the corresponding operators is characterized by rapid convergence and a good global search capability. The property of rapid convergence to a global optimum of ICA is used to accelerate the searching for the most suitable subset among a number of vehicle routing plans (Zhang et al. 2005). The parameters of the algorithm were denoted as follows: initial population of an immune body, N=100; number of codes in each immune body, \(a_{i}(it)=21\); terminating generation, \(g_{max}=100\); clone ratio, \(q=5\); the expected remaining population of an immune body, \(N_{n}=20\).

Given that objective 1 (i.e., Eq. (8)) carries much more importance than objective 2 (i.e., Eq. (9)) in the decision-making process, we assigned them different weights, \(w_{1}=\upalpha =0.9\) and \(w_{2}=(1-\upalpha )=0.1\), to combine them into the integrated objective function \(F=\upalpha F_{1}+(1-\upalpha )F_{2}\). After 50 generations of calculations, the output results were obtained as follows.

1.:

The optimal routes of the 3 vehicles were found to be \(0\rightarrow 1\rightarrow 5\rightarrow 0\) (for Vehicle 1), \(0\rightarrow 2\rightarrow 6\rightarrow 0\) (for Vehicle 2), and \(0\rightarrow 3\rightarrow 4\rightarrow 7\rightarrow 0\) (for Vehicle 3). The total transportation time of all vehicles is 6.05 h, and the corresponding cost is 6720 Yuan.

2.:

The vehicle arrival time, transportation costs and travel distances on each path are listed in Table 12.

Table 12 Vehicle routes and the corresponding costs and distances

The results indicate that the vehicles can reach all the fire sites in less than 2 h without compromising the efficiency of suppression. Meantime, on each path, a site with more severe fire conditions will be given a higher priority. The vehicles always travel from high priority sites to the low priority sites.

6 Discussions and conclusions

This paper proposes a new two-layer emergency logistics system for wildfire suppression involving a single depot and multiple demand sites. The real case analysis substantiates that during a wildfire suppression operation, the plans for allocating and transporting suppression resources have a strong impact on the efficiency of the operations. Compared with the existing findings in Lee et al. (2012), the resource deployment is optimized and shifted from the planning sites with the highest fire burned areas to the planning sites with the highest standard response requirements (i.e., highest priority). As a result, the optimal resource allocation and vehicle routing plans based on our optimization model are a better solution compared to the existing findings in the literature. Thus, the proposed optimization models provide powerful insights into emergency logistics planning for wildfire suppression operations. Because of the high research demand on emergency logistics for dynamic disaster relief (taking into account risk evolution), the findings from this study can also be generalized to tackle other logistics resource optimization issues in dynamic operational environments, particularly those sharing the same feature of a single supply and multiple demands in logistics planning and operations (e.g., allocation of ambulances and police forces).

The major contributions of this work are the ability of the models to analyze the emergency logistics operations by considering different fire propagation patterns and the factors affecting those patterns. Fire propagation models are constructed to calculate the fire propagation rate and the damaged area at each fire site. The outcomes can serve as suitable criteria for determining the priority of demand for suppression resources. An emergency logistics delivery plan based on dynamic fire propagation behavior is applied to address wildfire suppression in real time. According to historical data from wildfires that occurred in the Daxingan Mountains, two case studies are analyzed to verify the proposed models and validate the findings that different logistics planning schemes should be adopted for wildfire accidents with different fire propagation rates. The optimal results based on minimizing the disaster losses and logistics costs in the case of slow propagation demonstrate that transporting suppression resources to several fire sites using a single vehicle is more effective. In this case, the VRP is the most important concern. By contrast, in the case of rapid fire propagation, each fire site must be served by one or more dedicated vehicles; therefore, matching the fire sites to the corresponding vehicles is the most important consideration.

The study left some problems unsolved. First, when combining objective 1 and objective 2 into the integrated objective function, we consider \(\upalpha \) to be 0.9 in the empirical study and analyzed the case in which objective 1 carries much more importance than objective 2 in the decision process. Assigning \(\upalpha \) different values in a sensitivity analysis to test its influence on the final decision could reveal the optimal solutions in dynamic and different scenarios and thus aid in better decision making. Second, in the fire propagation model, we use the flame spread rate as the dominant indicator to predict the fire development pattern; however, several other indicators (e.g., fire intensity and flaming area et al.) could possibly have some effects on the fire propagation trend. Therefore, investigating how their individual/combined effects influence the accuracy of the fire propagation model is a worthwhile future direction. Finally, in the relief allocation and vehicle routing model, the impact of the conditions of wildfires suppression paths on the optimization results can be further investigated. Given that the conditions of the suppression paths may also affect the total cost of the emergency vehicles, a comparative study of the relief allocation problem and VRP under different suppression path conditions based on the existing research would be interesting.

Our possible future research plan includes studies of the relief allocation problem and VRP for wildfire suppression with multiple depots and a single demand site or multiple depots and multiple demand sites. This research plan may also involve comparative analysis the suitability of different optimization algorithms (e.g., the ant-colony algorithm and genetic algorithm) for solving these kinds of relief allocation and vehicle routing issues.