Quick-Response Model for Pre-and Post-Disaster Evacuation and Aid Distribution: The Case of the Tula River Flood Event

: In humanitarian logistics, efficiently evacuating people during disasters poses a significant challenge. Comprehensive plans are essential for moving individuals from affected areas, encompassing both pre-and post-disaster phases. These plans must optimize resource usage, including public transportation, and consider those with mobility challenges. This paper focuses on evacuating and assisting individuals in temporary shelters, specifically in the Tula River region, Hidalgo State, Mexico. Despite the recurring nature of this challenge, there is a noticeable gap in the scientific literature addressing quick-response methods for decision-makers to adapt existing spaces as temporary shelters and efficiently evacuate people from risk areas. To bridge this gap, we introduce a methodology aiming to minimize evacuation and aid distribution costs. Leveraging established algorithms like integer linear programming, the model determines shelter activation, while the vehicle routing problem assesses aid delivery strategies. Our research identified optimal evacuation routes from 13 affected areas to 34 shelters and aid distribution from collection centers to affected zones. We analyzed the average evacuation times for different demand scenarios: original, increased by 10%, and decreased by 10%, considering the transport units allocated and the distances on Google Maps. This study also evaluates the costs associated with humanitarian aid distribution under varying collection strategies involving state and municipal governments.


Introduction 1.Problem Statement
Currently, hazard events pose a significant threat to humanity.According to the United Nations, 90% of extreme events are linked to climate change, indicating a potential increase in their frequency year after year [1].These events, stemming from destructive origins, result in substantial economic impacts and human losses.In 2021, 401 disasters occurred worldwide, with the Pacific region bearing the brunt of these incidents [2].Various types of hazard events exist, with earthquakes, droughts, storms, floods, volcanic activity, extreme temperatures, landslides, and forest fires among the most common.Among these, storms and floods were the most frequent, accounting for 223 out of the 401 disasters in 2021 [3].As shown in Figure 1, floods ranked as the most common hazard events in Latin America and the Caribbean between 2000 and 2019 [4].
Due to the increased frequency of hazard events, the humanitarian supply chain (HSC) has gained significant attention in social, academic, and governmental circles due to its relevance in managing HSC operations, logistics coordination, traditional logistics, performance metrics, and modeling [5].The primary aim of the HSC is to study various problems, such as facility location, aid distribution, inventory management, and mass evacuation of people, with the overarching goal of safeguarding affected individuals and providing essential resources for their survival.Due to the increased frequency of hazard events, the humanitarian supply chain (HSC) has gained significant attention in social, academic, and governmental circles due to its relevance in managing HSC operations, logistics coordination, traditional logistics, performance metrics, and modeling [5].The primary aim of the HSC is to study various problems, such as facility location, aid distribution, inventory management, and mass evacuation of people, with the overarching goal of safeguarding affected individuals and providing essential resources for their survival.
Facility location pertains to the positioning of warehouses and distribution centers for emergency resources [6,7].Distribution studies delineate the formation of the distribution network, the choice of transportation modes, and the required capacity to assist people [8,9].Inventory-related issues revolve around determining resource levels, availability, cost minimization, and demand forecasting [10].Additionally, addressing the mass evacuation of people encompasses shelter location, evacuation point concentration, transport selection, traffic management, and the time required for evacuation [11,12].
Despite the growing interest in the HSC, the pre-disaster phase has received relatively less attention.However, this phase is of great significance in operations and disaster management, particularly concerning the mass evacuation of people [6,13].In the context of post-disaster decision-making, numerous proposals exist in the literature for optimizing the distribution of humanitarian aid using various modeling approaches.Nevertheless, they are often designed for diverse scenarios with numerous facilities, leading to complex decision-making processes and unresolved issues.Consequently, future research should explore the characterization of victims' needs and the systematic examination of multilevel service networks [12].
Effective evacuation planning in disaster scenarios must account for both the predisaster and post-disaster phases, ensuring the efficient utilization of resources, including public transportation, and addressing the needs of individuals with mobility difficulties and disabilities [14][15][16][17][18]. Notably, various HSC initiatives have been undertaken in Mexico across different domains, such as food banks, vaccine distribution, and food and donation collection, using diverse solution approaches.This study focuses on a case study methodology for evacuating people in a flood scenario.It employs a mathematical optimization model utilizing integer linear programming (ILP).Additionally, this study utilizes the vehicle routing problem (VRP) to assess various strategies for delivering humanitarian aid.
In the State of Hidalgo, Mexico, the Tula River poses a continuous flood threat to the municipalities that it traverses.Although annual flood events are not uncommon in this Facility location pertains to the positioning of warehouses and distribution centers for emergency resources [6,7].Distribution studies delineate the formation of the distribution network, the choice of transportation modes, and the required capacity to assist people [8,9].Inventory-related issues revolve around determining resource levels, availability, cost minimization, and demand forecasting [10].Additionally, addressing the mass evacuation of people encompasses shelter location, evacuation point concentration, transport selection, traffic management, and the time required for evacuation [11,12].
Despite the growing interest in the HSC, the pre-disaster phase has received relatively less attention.However, this phase is of great significance in operations and disaster management, particularly concerning the mass evacuation of people [6,13].In the context of post-disaster decision-making, numerous proposals exist in the literature for optimizing the distribution of humanitarian aid using various modeling approaches.Nevertheless, they are often designed for diverse scenarios with numerous facilities, leading to complex decision-making processes and unresolved issues.Consequently, future research should explore the characterization of victims' needs and the systematic examination of multilevel service networks [12].
Effective evacuation planning in disaster scenarios must account for both the predisaster and post-disaster phases, ensuring the efficient utilization of resources, including public transportation, and addressing the needs of individuals with mobility difficulties and disabilities [14][15][16][17][18]. Notably, various HSC initiatives have been undertaken in Mexico across different domains, such as food banks, vaccine distribution, and food and donation collection, using diverse solution approaches.This study focuses on a case study methodology for evacuating people in a flood scenario.It employs a mathematical optimization model utilizing integer linear programming (ILP).Additionally, this study utilizes the vehicle routing problem (VRP) to assess various strategies for delivering humanitarian aid.
In the State of Hidalgo, Mexico, the Tula River poses a continuous flood threat to the municipalities that it traverses.Although annual flood events are not uncommon in this region, one of the most severe floods in recent memory occurred on 6 September 2021, which was classified as a natural disaster.Unfortunately, this disaster claimed the lives of 15 people and left over 31,000 homes severely affected [19,20].According to government authorities, the disaster was precipitated by heavy rains in Mexico City and the State of Mexico, leading to the saturation of the drainage system.As a result, the National Water Commission activated the Joint Operation Protocol for Deep Drainage of the Valley of Mexico, causing the Tula River to swell, overflow its banks, and significantly increase its flow.This overflow affected several municipalities, including Tula de Allende, Ixmiquilpan, Tlaxcoapan, Tezontepec de Aldama, Chilcuautla, Tasquillo, Tlahuelilpan, Tepeji del Río, and Mixquiahuala [21].
The drainage system of the metropolitan area in the Valley of Mexico is interconnected through a network of tunnels and canals, facilitating the discharge of water into the Mezquital Valley via the Tula River.However, this setup has adverse consequences for managing urban stormwater.The problem extends beyond just engineering considerations and encompasses a broader perspective that involves public policy [22].Effectively managing the flood risk stemming from the drainage system of the metropolitan area of the Valley of Mexico requires a comprehensive approach.This approach should encompass geographical and interdisciplinary dimensions, including the construction of defense structures and diversion gates, river realignment, and dredging.Moreover, it should integrate additional measures aimed at reducing the severity of the damage, such as improved damage control and mitigation efforts.Public awareness and preparedness are crucial aspects, ensuring that the population is well-informed about flood risks and is aware of evacuation and safety zones to protect human lives [23].

Problem Description
This paper presents a comprehensive strategy to evacuate residents from the municipalities impacted by the Tula River prior to a disaster, and to offer humanitarian assistance in shelters after the disaster.This approach aims to minimize costs and collaboratively plan the essential resources, including determining the number and placement of shelters based on calculations of the affected population and establishing routes for delivering aid to these shelters.Importantly, this strategy is grounded in real-world data, incorporating distances and costs, which can serve as valuable decision-making tools for both state and federal governments.
As will be observed in Section 2, there are few works that provide real-world case studies combining both pre-and post-disaster stages.In fact, only Mollah et al. [24] work presents such a case, but their solution method is complex.Although the solution methods utilized in this work have been used before, the motivation here is to explore a straightforward approach that allows for the planning of which shelters to open, how to evacuate the most vulnerable individuals, and how to distribute humanitarian aid.Following a review of the state of the art, the most significant contributions can be summarized as follows: • A simple decision-making model is proposed.

•
Real demographic and social data are employed.

•
Consideration is given to the most vulnerable individuals, specifically those over 60 years of age and people with physical disabilities.

•
Both the pre-and post-disaster stages are taken into account.
The methodology involved data collection, building the evacuation model, solving the model, and generating options for aid distribution; the research methodology is detailed in Section 3. Section 4 describes the results and discussion, while Section 5 presents the conclusions.In the next section, the literature review is presented.

Literature Review
In recent years, there has been a substantial surge in publications pertaining to humanitarian supply chains (HSCs).These publications have addressed various key aspects, including facility location, distribution models, inventory management, and mass evacuation procedures.Furthermore, different mathematical models have been put forth for the evacuation of people, employing a range of algorithmic solutions tailored to the unique characteristics of the specific problem at hand.Within the context of HSCs, there exist conditions of uncertainty, particularly in relation to the demand for assistance.Addressing such uncertainty is crucial when devising evacuation models.
In the field of humanitarian logistics, various models have been developed to address diverse objectives.These objectives include the allocation of ambulances under uncertain conditions, the minimization of vehicle numbers, travel cost reduction, limiting maximum latency, optimizing shelter locations, and efficient distribution of aid.Additionally, challenges related to facility location, distribution, inventory management, and mass evacuation have been tackled using multi-objective optimization and multi-period distribution planning, often employing metaheuristic techniques like the Non-Dominated Sorting Genetic Algorithm III, simulated annealing (SA), and variable neighborhood search (VNS) due to the inherent complexity [24][25][26][27][28][29].
The humanitarian supply chain (HSC) is characterized by four distinct phases: prevention, preparedness, response, and reconstruction and recovery.Decisions made in the pre-disaster phase significantly influence those in the post-disaster phase.The location of distribution centers, the availability of inventory, and the positioning of shelters all play pivotal roles in the distribution of relief supplies and the evacuation of people.Concerning both the pre-disaster and post-disaster phases, problems related to stock prepositioning, facility placement, evacuation planning, and hospital and distribution center locations have been addressed.These complex challenges have been approached using various metaheuristics, including classical approaches (CAs), pattern search algorithms (PSAs), genetic algorithms (GAs), and Non-Dominated Sorting Genetic Algorithm III [15,[28][29][30].
One of the most critical issues within the HSC domain is the efficient mass evacuation of people for providing timely care to victims.The literature offers proposals to minimize the costs of rescue operations while giving priority to injured individuals.These proposals also cover scheduling and assigning rescue teams, taking into account factors like fatigue and deprivation time.Metaheuristic methods, such as genetic algorithms (GAs) and particle swarm optimization (PSO), have been employed to solve these models [16][17][18].
Among the recent contributions, Shehadeh and Tucker [31] focused on determining the location of humanitarian aid warehouses in both the pre-and post-disaster phases, considering uncertainty, and positioning inventory.Seraji et al. [32] optimized the placement of emergency shelters and coordinated the movement of aid vehicles using two-stage multi-objective metaheuristics.Table 1 provides a summary of the articles found in the literature that explore approaches in both pre-and post-disaster scenarios, along with the solution methods, some of which have been applied in real-world cases.Other works involving individuals with disabilities include the one by Yazdani and Haghani [34], where they presented a case study in Australia for evacuating the elderly using optimization algorithms combined with behavioral training.The same authors also introduced an emergency plan for relocating vulnerable individuals to shelters in response to heatwaves, presenting a method that they tested in a set of randomly generated locations [35].
While more methods with broader scope exist, involving larger networks and traffic modeling, and considering human behavior [36], in the case of the Tula River, the network only encompasses 50 possible shelters, 9 assembly points, and 2 distribution centers for aid, also without different route options.Thus, integer programming models are sufficient for making quick decisions.The concept of nearest and safest allocation, as outlined by Southworth [37] and Barrett et al. [38], is employed, and this planning is conducted in advance of flooding.

Research Methodology
The methodology to address the problem proposed in this work is composed of four stages: (1) data collection, (2) evacuation model, (3) model solution, and (4) distribution of humanitarian aid (Figure 2).
Logistics 2024, 8, x FOR PEER REVIEW 6 of 21 Step (4)-Distribution of Humanitarian Aid: This phase considers the use of the VRP technique to determine the collection routes for the distribution of humanitarian aid from the distribution centers to the shelters.In addition, the routes for the evacuation of people are determined and the time required for said evacuation is calculated.

Data Collection
In the initial phase, the municipalities situated closest to the Tula River were identified.Subsequently, the 13 municipalities that were deemed to be the most severely affected were chosen.The municipalities considered in the State of Hidalgo were Alfajayucan, Atitalaquia, Atotonilco de Tula, Ixmiquilpan, Mixquiahuala, Tasquillo, Tepeji del Río, Tezontepec de Aldama, Tlahuelilpan, Tlaxcoapan, Tula de Allende and Zimapán, and for the State of Mexico we considered the municipality Apaxco.
The method for estimating the number of people requiring evacuation relied on data sourced from the 2020 Population and Housing Census, which was conducted by the National Institute of Statistics and Geography (INEGI 2021) [39].Notably, this analysis considered two specific groups: individuals aged 60 years or older, and those with physical disabilities, as these populations are the most vulnerable during a disaster.The calculation of the number of people to be evacuated for each municipality involved summing the percentages of people aged over 60 and people with physical disabilities, and then multiplying the total by the municipality's overall population.Additional details can be found in Table 2.  Step (1)-Data Collection: The initial phase involves complete data collection to establish a solid foundation of understanding of the problem for further analysis.The data collected cover the calculation of the population that requires evacuation, the identification of possible shelters, the location of distribution centers for humanitarian aid, and the calculation of distances between affected areas and shelters.
Step (2)-Evacuation Model: This step involves the formulation of the optimization model using integer linear programming.It begins with the calculation of transportation costs for the evaluation of people and the distribution of aid, as well as the calculation of the budget for the qualification of the shelters.These costs are the basis for the optimization of the model; subsequently, the variables involved in the mathematical model are established, and the objective function and the set of restrictions are formulated.
Step (3)-Model Solution: The third phase includes the search for tools to provide a solution to the proposed model and establish the procedure to perform the sensitivity analysis.
Step (4)-Distribution of Humanitarian Aid: This phase considers the use of the VRP technique to determine the collection routes for the distribution of humanitarian aid from the distribution centers to the shelters.In addition, the routes for the evacuation of people are determined and the time required for said evacuation is calculated.

Data Collection
In the initial phase, the municipalities situated closest to the Tula River were identified.Subsequently, the 13 municipalities that were deemed to be the most severely affected were chosen.The municipalities considered in the State of Hidalgo were Alfajayucan, Atitalaquia, Atotonilco de Tula, Ixmiquilpan, Mixquiahuala, Tasquillo, Tepeji del Río, Tezontepec de Aldama, Tlahuelilpan, Tlaxcoapan, Tula de Allende and Zimapán, and for the State of Mexico we considered the municipality Apaxco.
The method for estimating the number of people requiring evacuation relied on data sourced from the 2020 Population and Housing Census, which was conducted by the National Institute of Statistics and Geography (INEGI 2021) [39].Notably, this analysis considered two specific groups: individuals aged 60 years or older, and those with physical disabilities, as these populations are the most vulnerable during a disaster.The calculation of the number of people to be evacuated for each municipality involved summing the percentages of people aged over 60 and people with physical disabilities, and then multiplying the total by the municipality's overall population.Additional details can be found in Table 2. Regarding the selection of temporary shelters, public spaces like auditoriums, schools, and churches were considered, with a preference for locations situated at a significant distance from the Tula River.Initially, 50 such sites were identified.However, after a thorough analysis of infrastructure and distances, the list was refined to 37.These selected facilities were subsequently integrated into the mathematical optimization model.For a comprehensive list of these shelters and their respective capacities, please refer to Table A1 in Appendix A.
To determine the locations of distribution centers, two distinct strategies were evaluated: The first strategy involved utilizing municipal presidencies as collection points, from which aid would be dispatched to various shelters.This approach was made feasible by strategically stockpiling a certain amount of food during the disaster preparedness phase.In the second strategy, the team considered the Government Palace, the National System for the Integral Development of the Family (DIF), and the Tula Presidency as distribution centers.This approach assumes that there exists a disaster relief fund that could be employed to directly dispatch aid kits to the shelters without necessitating intermediary distribution centers.The shelters, affected areas, and distribution centers are shown on the map in Figure 3.
To calculate the distances between the affected areas and the temporary shelters, we relied on real-time data obtained through Google Maps coordinates.This methodology was also applied to determine the distances between the distribution centers and the shelters.More detailed information is accessible through the provided repository [40].

Evacuation Model
In the second phase, the transportation costs were calculated.These costs were derived from selecting a particular type of transportation (NISSAN/URVAN/2015) with a capacity of 15 passengers.To compute the cost per trip, the fuel cost was used, which depended on the distance traveled, and a fixed cost was added for renting the vehicle.The fuel consumption per kilometer was determined based on data provided by the Ministry of Environment and Human Resources, while the fixed cost was obtained through price consultations with transportation providers [41].
The budget calculation for equipping a shelter involved selecting a closed facility and determining the necessary items for its preparation.The essential requirements included portable toilets, sinks, a semi-industrial kitchen, aluminum cookware, basic utensils, plastic sinks, mattresses, bedspreads, blankets, waste containers, tents, tables, and chairs.The calculation was made to accommodate approximately 150 to 200 people, resulting in a budget of MXN 159,147.Detailed information can be found in Table A2 in Appendix A.
The formulation of the optimization model employed integer linear programming (ILP) for evacuating people.In this model, the decision variables included the assignment of individuals to be evacuated from each affected area to each temporary shelter and the quantity of aid to be dispatched from each distribution center to each temporary shelter.The primary objective was to minimize transportation costs.
points, from which aid would be dispatched to various shelters.This approach was made feasible by strategically stockpiling a certain amount of food during the disaster preparedness phase.In the second strategy, the team considered the Government Palace, the National System for the Integral Development of the Family (DIF), and the Tula Presidency as distribution centers.This approach assumes that there exists a disaster relief fund that could be employed to directly dispatch aid kits to the shelters without necessitating intermediary distribution centers.The shelters, affected areas, and distribution centers are shown on the map in Figure 3.The model is as follows: Equation ( 1) represents the objective function that minimizes the evacuation and the humanitarian aid distribution costs, considering the distances from the affected municipalities to the shelters and the distances from the distribution centers to the shelters.Additionally, the rental cost of the transport units is considered for evacuating people to the shelters.Equation (2) ensures that the number of people evacuated from municipality i is less than the capacity of shelter j.Equation (3) establishes that the total number of people who need to be evacuated must be sent to a shelter.Equation ( 4) ensures that the number of humanitarian aid kits that must be sent from distribution center k to shelter j is equal to the number of people sent to this shelter.Equation (5) establishes that the number of humanitarian aid kits that must be sent from distribution center k to shelter j is less than the capacity of the distribution center.Equation (6) allows that the costs of evacuating people and the distribution of humanitarian aid are less than the budget allocated to manage a flood disaster.Finally, Equation (7) establishes that the variable y j is 1 if the shelter opens and 0 otherwise.
Table 3 provides the notation used in the ILP model proposed for the evacuation problem.

Model Solution
To address the ILP model for the evacuation described in Equations ( 1)-( 7), we employed Lingo 17.0 software, a tool specifically designed for constructing and solving mathematical optimization models.
To conduct sensitivity analysis, the number of people to be evacuated was used as an input parameter in the Lingo model.After obtaining the optimal response with the original data, the model was solved with a 10% variation in demand (both increasing and decreasing).This was performed to examine how the total cost of evacuation and humanitarian aid delivery behaves under different demand scenarios.

Distribution of Humanitarian Aid
In this research, a vehicle routing problem (VRP) was used for determining the routes of the collection centers.The VRP can be represented as the following graph theory problem: let G = (V, A) be a complete graph, where V = {0, 1, ..., n} is the set of vertices and A is the set of arches.The vertices j = 1, ..., n correspond to customers, each with a known non-negative demand, while vertex 0 corresponds to the deposit.A non-negative cost is associated with each arc and represents the cost of traveling from vertex i to vertex j.The VRP consists of finding a set k of simple circuits, each one corresponding to a vehicular route with a minimum cost defined as the sum of the costs of the arcs of the circuits, in such a way that (a) each circuit visits vertex 0, (b) each vertex j is visited by exactly one circuit, and (c) the sum of the demands of the vertices visited by a circuit does not exceed the capacity of the vehicle [42].
The VRP that we used considered the allocation of a vehicle at each distribution center, specifically at the State of Hidalgo Government, DIF Hidalgo, and the Presidency of Tula.The investigation encompassed 34 shelters characterized by normal demand.To obtain real coordinates and estimated travel times, while accounting for traffic conditions, Bing Maps data were utilized.The analysis assumed that each humanitarian aid kit weighed approximately 1 kg and that articulated vehicles with a capacity of 25,000 kg were employed, allowing each vehicle to transport 25,000 kits.These vehicles were assumed to possess a fuel efficiency of 100 km per 35 L of diesel [43], resulting in a cost of MXN 8.28 per kilometer traveled.Additionally, a delivery cost of MXN 500 was incorporated for each humanitarian aid drop-off.The VRP model was implemented with a single vehicle at each distribution center.Subsequently, the number of vans was determined through sensitivity analysis.To estimate the hours required for evacuation, the research considered the deployment of 5, 10, 15, and 20 vans.Furthermore, a 10% variation in demand was considered, encompassing both decreases and increases.
This methodology provides a well-structured guide to address the identified problem in a systematic and effective way.It should be noted that this methodology was designed for quick responses in pre-and post-disaster evacuation and aid distribution, with few temporary shelters and a small population.

Evacuation Model
To solve the ILP for the evacuation model proposed in Equations ( 1)-( 7), Lingo 17.0 software was used.Solving the model through the branch-and-bound solver included in Lingo, the optimal allocation was obtained from each of the affected municipalities (concentration centers for victims) to the different shelters, and at the same time the optimal distribution of aid was determined.It should be noted that a sensitivity analysis was carried out by modifying the demand parameters, with a variation of 10% above and below, and these were used as input parameters in Lingo.The objective function with the initial demand was MXN 2,969,813.00,with decreasing demand by 10% it was MXN 2,635,202.00,and with increasing demand by 10% it was MXN 3,250,149.00.It can be seen that with a small modification of the demand, the cost of evacuation and distribution of aid increases by 9% (MXN 280,336.00)with respect to the initial demand because, to meet the demand, three more shelters must be opened (see Figure 4).
The solutions for the three different scenarios of demand are shown in the repository [40]; only 34 shelters were opened, and the shelter information can be found in Table 4. Considering the solution in the repository for the original demand, the people who need to be evacuated from a central location can be determined; for example, the Tula Presidency moves 4500 people (p) to the Tula Sports Unit, 4000 p to the Tula Municipal Auditorium, 4000 p to the Infonavit San Marcos Auditorium, 2000 p to the Cathedral of Tula, 1342 p to the Benito Juárez Elementary School, and 2000 p to the Miguel Hidalgo y Costilla Elementary School, giving a total of 17,842 p.

Evacuation Time and Vans Required
Based on the best results in the provided repository [40], the required evacuation time was calculated, considering variations in the number of vans for the transfer of people; it was considered that each van has a transfer capacity of 15 people.In order to know the number of hours required for evacuation, 5, 10, 15, and 20 vans were considered, as well as considering a 10% variation in the demand (decrease and increase).Since the municipalities of Tula, Tepeji, and Ixmiquilpan are the ones with the largest numbers of people to evacuate, they are also the ones that require the longest time.For example, if 5 vans are assigned to the Tepeji Presidency, considering the average demand, 62.04 h is required to evacuate, while with 10 vans 31.0244h is required, and with 20 vans 15.51 h is required, changing proportionally.The total time varies between 8.65% and 11.38% less when the demand goes down and increases by between 25.48% and 29.34% when the demand increases (see Table 5).The solutions for the three different scenarios of demand are shown in the repository [40]; only 34 shelters were opened, and the shelter information can be found in Table 4. Considering the solution in the repository for the original demand, the people who need to be evacuated from a central location can be determined; for example, the Tula Presidency moves 4500 people (p) to the Tula Sports Unit, 4000 p to the Tula Municipal Auditorium, 4000 p to the Infonavit San Marcos Auditorium, 2000 p to the Cathedral of Tula, 1342 p to the Benito Juárez Elementary School, and 2000 p to the Miguel Hidalgo y Costilla Elementary School, giving a total of 17,842 p.To provide an adequate response to the evacuation problem in the face of a natural disaster, time is a key factor.In this sense, the numbers of vans required to evacuate in smaller timeframes (3, 5, and 8 h) were calculated.Table 5 shows that when demand decreases, the number of trucks needed decreases by between 10.65% and 11.11%, while if demand increases, the number of trucks needed increases by between 24.71% and 25.12%.The number of vans required to evacuate the population using the average demand would be 418 for a time of 5 h, while 3 h would require 693 vans, and for 8 h 263 vans would be needed; for more information, consider Table A5.Concerning the strategy evaluation for using different distribution centers using the VRP, the VRP Spreadsheet Solver v3.8 was used.Two different strategies were evaluated.The first strategy made use of three collection centers: two in Pachuca (DIF Hidalgo and State of Hidalgo Government) and the third center in the Tula Presidency, each with one articulated vehicle with a capacity of 25,000 kits; the solution is shown in Tables A6-A8 in Appendix A, where the routes to be followed, approximate time, and total cost each time that the aid is delivered (MXN 19,646) are shown.The second strategy made use of only two collection centers: DIF Hidalgo, with one vehicle, and a second center in the Tula Presidency with two vehicles.The routes are shown in Tables A7-A9.The second strategy was the best, with a cost of MXN 18,611 each time that the kits are delivered.

Conclusions
Currently, hazard events represent a great threat to humanity.A further increase in their incidence is expected due to climate change.This is why the HSC is gaining more and more interest in academia, government, and society.Some relevant problems are the location of shelters, the distribution of aid, inventories, and the mass evacuation of people.
In the literature reviewed here, each author adapted the model and solution strategies to the specific case that they were trying to solve.In this article, shelters are selected, people are evacuated to the shelters and, finally, strategies to distribute aid are explored considering demographic and real geographic data for the evacuation of people with mobility difficulties and disabilities.
The flood that occurred on 6 September 2021 was considered to be an extreme event, and year after year there are flooding problems.The main idea of this study was to provide strategies for municipal and state governments to make decisions about how many shelters to open and how to distribute the aid.The methodology consisted of calculating the population to be evacuated, identifying the possible shelters that could be set up, deciding where it would be convenient to locate the distribution centers, calculating the distance between locations and the expenses incurred to adapt a temporary shelter, and formulating the optimization model to decide which shelters would be opened and how people would be evacuated.Finally, a VRP with different strategies to distribute humanitarian aid was presented.
An analysis was conducted with three types of demand-the original, increased by 10%, and decreased by 10%-to evaluate the adequacy of the planned shelters.Furthermore, we calculated the number of vans needed to execute the evacuation.Our hope, as natives of this region, is that this research serves as a foundation for decision-making and can be employed to conduct similar analyses in other communities during extreme events.Some of the limitations of this work include the reliance on Google Maps distances for the integer linear programming (ILP), which do not account for real-time traffic conditions.Future research may involve considering the evacuation of other vulnerable individuals in the region, such as those in hospitals, prisons, etc., in the event of other extreme events.Additionally, it could explore how aerial evacuations would be conducted in situations where all routes are obstructed and strategize on the organization of brigades within shelters and/or evacuation points.

Data Availability Statement:
The data are available at the repository https://doi.org/10.6084/m9.figshare.23680881.v1(accessed on 12 December 2023) [40] and can be requested from the corresponding author when necessary.

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

Logistics 2024, 8 , 21 Figure 1 .
Figure 1.Most common hazard events in Latin America and the Caribbean from 2000 to 2019, by number of occurrences.Source: [4].

Figure 1 .
Figure 1.Most common hazard events in Latin America and the Caribbean from 2000 to 2019, by number of occurrences.Source: [4].

Figure 3 .
Figure 3. Shelters (blue dots), municipal presidencies (red dots) and distribution centers (green dots).Source: own elaboration in Gurobi-PythonTo calculate the distances between the affected areas and the temporary shelters, we relied on real-time data obtained through Google Maps coordinates.This methodology

Figure 4 .
Figure 4. Solution of the evacuation model in Lingo.Source: own elaboration.

Figure 4 .
Figure 4. Solution of the evacuation model in Lingo.Source: own elaboration.

Table 1 .
Models proposed for the attention of natural disasters.

Table 2 .
Calculation of the evacuee population.

Table 2 .
Calculation of the evacuee population.

Table 3 .
Notation used in the linear programming model for the evacuation problem.

Table 5 .
Time for evacuation, original demand.
a Five vans assigned to every presidency; b hours needed for evacuation; c time in hours; d vans required.

Table A1 .
Selection of properties for the location of shelters.Source: own work.

Table A2 .
Items required to outfit a shelter.

Table A3 .
Numbers of trucks required for evacuation within 3, 5, and 8 h.

Table A4 .
Routes and approximate times for delivery of kits by vehicle 1-first strategy.

Table A5 .
Routes and approximate times for delivery of kits by vehicle 2-first strategy.

Table A6 .
Routes and approximate times for delivery of kits by vehicle 3-first strategy.

Table A7 .
Routes and approximate times for delivery of kits by vehicle 1-second strategy.

Table A8 .
Routes and approximate times for delivery of kits by vehicle 2-second strategy.

Table A9 .
Routes and approximate times for delivery of kits by vehicle 3-second strategy.