A Typhoon Shelter Selection and Evacuee Allocation Model: A Case Study of Macao (SAR), China

: Typhoon disaster represent one of the most prominent threats to public safety in the Macao Special Administrative Region (SAR) of China and can cause severe economic losses and casualties. Prior to the landing of typhoons, a ﬀ ected people should be evacuated to shelters as soon as possible; this is crucial to prevent injuries and deaths. Various models aim to solve this problem, but the characteristics of disasters and evacuees are often overlooked. This study proposes a model based on the inﬂuence of a typhoon and its impact on evacuees. The model’s objective is to minimize the total evacuation distance, taking into account the distance constraint. The model is solved using the spatial analysis tools of Geographic Information Systems (GIS). It is then applied in Macao to solve the evacuation process for Typhoon Mangkhut 2018. The result is an evacuee allocation plan that can help the government organize evacuation e ﬃ ciently. Furthermore, the number of evacuees allocated to shelters is compared with shelter capacities, which can inform government shelter construction in the future.


Introduction
One major consequence of global warming is the increased intensity and frequency of extreme weather events such as typhoons [1,2]. Constructing shelters important for disaster management, as shelters can house affected people before and during disasters. In this way, casualties can be reduced. The southeastern coastal area of China is highly vulnerable to typhoons, and the threat to people and property in this area is increasing [3]. As one of the wealthiest and most booming regions in China, the Macao Special Administrative Region (SAR) has been seriously affected by typhoon disasters (especially Typhoon Hato, in 2017, and Typhoon Mangkhut, in 2018). Therefore, constructing shelters, planning efficient evacuation and planning evacuee allocation are of great importance to the Macao SAR government.
Engineering techniques can mitigate the impacts of disasters such as earthquakes and typhoons [4]. However, in server cases, these techniques can no longer protect people. Therefore, the establishment of shelters and the allocation of evacuees to them is effective [5,6]. To determine shelter locations, A mathematical model, using the spatial analysis technology of GIS, enables the selection of the shelters to which evacuees from all buildings are allocated. The 2D version of Flood Map [30,31] was used to determine the affected range and the number of evacuees. Sea level and digital elevation model (DEM) data were used in 2D Flood Map to determine the inundation depth and range every fifteen minutes. The sea level data was the Macao government's observation data, from Barra Station. This provides sea level every 5 min. The digital elevation data used for inundation modeling was provided by the Macao government; the resolution was 5 m and covered the whole city area. The inundation depth and range generated per 15 min using 2D Flood Map; buildings located in highly risky areas (maximum inundation depth > 1 m) were then selected from all buildings in Macao. To select shelters and allocate evacuees, the first step is to establish a road network for the city. This can be done from the path data, using a network analysis tool in ArcGIS. The second dataset is the shelter layer, which lists the capacities and locations of shelters. The third dataset is the 'demand' layer, which lists both building locations and the number of evacuees per building. The mathematical model determines how evacuees are allocated to the shelters.
In the following, a mathematical model of evacuee allocation to typhoon shelter is developed. A solver model is also described, built with the GIS ModelBuilder tool to solve the mathematical model. There are three assumptions: (1) Residents of a given building will be allocated to the same shelter; (2) Residents will go to their assigned shelters along the shortest possible route; (3) Residents will be allocated to shelters if the maximum inundation depths of their buildings exceed 1 m.

Mathematical Model for Evacuee Allocation
The objective of the mathematical model is to minimize the total evacuation distance from all buildings, to ensure affected people reach their shelters as soon as possible.
Our notation is as follows: j: Index of building, equal to 1, 2, …, M k: Index of shelter, equal to 1, 2, …, N dkj: Shortest evacuation route distance from building j to shelter k Bkj: Whether building j is assigned to shelter k (1 = yes, 0 = no) Dj: The maximum evacuation distance that residents in building j can travel Xk: Whether shelter k is selected (1 = yes, 0 = no) To select shelters and allocate evacuees, the first step is to establish a road network for the city. This can be done from the path data, using a network analysis tool in ArcGIS. The second dataset is the shelter layer, which lists the capacities and locations of shelters. The third dataset is the 'demand' layer, which lists both building locations and the number of evacuees per building. The mathematical model determines how evacuees are allocated to the shelters.
In the following, a mathematical model of evacuee allocation to typhoon shelter is developed. A solver model is also described, built with the GIS ModelBuilder tool to solve the mathematical model. There are three assumptions: (1) Residents of a given building will be allocated to the same shelter; (2) Residents will go to their assigned shelters along the shortest possible route; (3) Residents will be allocated to shelters if the maximum inundation depths of their buildings exceed 1 m.

Mathematical Model for Evacuee Allocation
The objective of the mathematical model is to minimize the total evacuation distance from all buildings, to ensure affected people reach their shelters as soon as possible.
Our notation is as follows: j: Index of building, equal to 1, 2, . . . , M k: Index of shelter, equal to 1, 2, . . . , N d kj : Shortest evacuation route distance from building j to shelter k B kj : Whether building j is assigned to shelter k (1 = yes, 0 = no) D j : The maximum evacuation distance that residents in building j can travel X k : Whether shelter k is selected (1 = yes, 0 = no) Equation (1) is the objective to minimize the total evacuation distance of all affected buildings; d kj is the shortest distance from affected building j to shelter k, and is calculated based on the route data. Affected buildings are defined as buildings that are inundated in the flood with a depth > 1 m. The buildings' distribution and the number for each building were obtained from the Macao government. The inundation depth was calculated by modeling the coastal flooding of Typhoon Mangkhut 2018 using the 2D version of Flood Map [30,31]. As previously explained, to calculate the inundation range and depth, DEM data with resolution of 5 m, basic geographic data, and the sea level data per five minutes were used. All data were provided by the government, with sea level data from the Barra observation station in Macao. The resolution of DEM was high and the sea level data was detailed, allowing a high-quality inundation result.
Decision variable B kj indicates whether the building j will be allocated to shelter k (1 allocated, 0 otherwise). Equation (2) is the distance constraint that a shelter will not be selected for a building if the distance between them is greater than D j . This represents the furthest residents in the building can travel. Equation (3) ensures that evacuees from a given building can only be allocated to one shelter; X k is a decision variable that indicates whether the shelter k will be selected by any building under the typhoon scenario (1 = selected, 0 = not selected).

Solver Model Made by GIS Model Builder
To solve the mathematical model, many approaches can be used, including the spatial analysis technology of GIS and optimum algorithm. The mathematical model proposed in this paper involves just one objective and the constraints are not very complex. Therefore, the ModelBuilder tool of GIS was used to solve the mathematical model. This is built based on the spatial analysis technology of GIS. The model named "shortest path" was built using the ModelBuilder tool as shown in Figure 2. Firstly, the path network was generated using GIS based on the Macao path data. Then, a network analysis tool named "make location-allocation layer" was employed, inputting the path network data. Meanwhile, shelter and building location data were provided and the model was solved. The results show which shelter will be selected by evacuees of each building for refuge. Equation (1) is the objective to minimize the total evacuation distance of all affected buildings; dkj is the shortest distance from affected building j to shelter k, and is calculated based on the route data. Affected buildings are defined as buildings that are inundated in the flood with a depth > 1 m. The buildings' distribution and the number for each building were obtained from the Macao government. The inundation depth was calculated by modeling the coastal flooding of Typhoon Mangkhut 2018 using the 2D version of Flood Map [30,31]. As previously explained, to calculate the inundation range and depth, DEM data with resolution of 5 m, basic geographic data, and the sea level data per five minutes were used. All data were provided by the government, with sea level data from the Barra observation station in Macao. The resolution of DEM was high and the sea level data was detailed, allowing a high-quality inundation result.
Decision variable Bkj indicates whether the building j will be allocated to shelter k (1 allocated, 0 otherwise). Equation (2) is the distance constraint that a shelter will not be selected for a building if the distance between them is greater than Dj. This represents the furthest residents in the building can travel. Equation (3) ensures that evacuees from a given building can only be allocated to one shelter; Xk is a decision variable that indicates whether the shelter k will be selected by any building under the typhoon scenario (1 = selected, 0 = not selected).

Solver Model Made by GIS Model Builder
To solve the mathematical model, many approaches can be used, including the spatial analysis technology of GIS and optimum algorithm. The mathematical model proposed in this paper involves just one objective and the constraints are not very complex. Therefore, the ModelBuilder tool of GIS was used to solve the mathematical model. This is built based on the spatial analysis technology of GIS. The model named "shortest path" was built using the ModelBuilder tool as shown in Figure 2. Firstly, the path network was generated using GIS based on the Macao path data. Then, a network analysis tool named "make location-allocation layer" was employed, inputting the path network data. Meanwhile, shelter and building location data were provided and the model was solved. The results show which shelter will be selected by evacuees of each building for refuge.

Case Study
Macao, which is located in South China, is frequently affected by tropical cyclones. The typhoon observation record begins in 1953, but the effect of tropical cyclones on Macao was minimal until 2017, when Typhoon Hato caused the worst damage in 64 years-killing 10 injuring dozens and sweeping away 11,470 million MOP [32]. Although the government then put effort into typhoon disaster mitigation, Typhoon Mangkhut subsequently occurred in 2018, causing 40 injuries and making nearly 6000 people homeless. During typhoon disasters, most injuries and deaths are caused

Case Study
Macao, which is located in South China, is frequently affected by tropical cyclones. The typhoon observation record begins in 1953, but the effect of tropical cyclones on Macao was minimal until 2017, when Typhoon Hato caused the worst damage in 64 years-killing 10 injuring dozens and sweeping away 11,470 million MOP [32]. Although the government then put effort into typhoon disaster mitigation, Typhoon Mangkhut subsequently occurred in 2018, causing 40 injuries and making nearly 6000 people homeless. During typhoon disasters, most injuries and deaths are caused by inundation of houses. Therefore, providing reasonable allocation schemes for evacuees in Macao is important to avoid losses caused by typhoon disasters.
Typhoon Mangkhut affected the entire territory of Macao. Its wind speed reached 42-48 m/s (a magnitude of [14][15] and it caused nearly 108 mm of precipitation. It caused a serious storm surge, with a maximum sea level rise of 5.21 m. The storm surge caused inundation, and the maximum inundation depth in Macao was up to 2.01 m. During the typhoon, 598 accidents happened and were dealt with by the government. Additionally, 5650 residents was evacuated, of which 1346 lived in shelters. In this paper, Typhoon Mangkhut was selected to simulate an inundation situation and determine how to allocate evacuees to shelters.
Macao has a land area of 32.9 square kilometers; this consists of the Macao Peninsula, Taipa, and Coloane. Besides, a 1 km 2 parcel of land in neighboring Hengqin island also falls under the regional government's jurisdiction. Figure 3a indicates the location of Macao in China and Figure 3b illustrates its administrative divisions.
Sustainability 2020, 12, x FOR PEER REVIEW 5 of 13 by inundation of houses. Therefore, providing reasonable allocation schemes for evacuees in Macao is important to avoid losses caused by typhoon disasters. Typhoon Mangkhut affected the entire territory of Macao. Its wind speed reached 42-48 m/s (a magnitude of [14][15] and it caused nearly 108 mm of precipitation. It caused a serious storm surge, with a maximum sea level rise of 5.21 m. The storm surge caused inundation, and the maximum inundation depth in Macao was up to 2.01 m. During the typhoon, 598 accidents happened and were dealt with by the government. Additionally, 5650 residents was evacuated, of which 1346 lived in shelters. In this paper, Typhoon Mangkhut was selected to simulate an inundation situation and determine how to allocate evacuees to shelters.
Macao has a land area of 32.9 square kilometers; this consists of the Macao Peninsula, Taipa, and Coloane. Besides, a 1 km 2 parcel of land in neighboring Hengqin island also falls under the regional government's jurisdiction. Figure 3a indicates the location of Macao in China and Figure 3b illustrates its administrative divisions.   Figure 4b shows the locations of 6 shelters and 390 buildings in Taipa and Coloane. The sizes of shelter symbols correspond to their capacities, which are listed in Table 1. These 17 shelters were constructed by Macao's government and are operated during typhoon disasters. The residents of the buildings need to be allocated to their designated shelters when the inundation depth is more than 1 m. The figure shows that most of the buildings are in the Macao Peninsula, while the buildings in Taipa and Coloane are much less dense.
There are 11 shelters in the Macao Peninsula and 6 shelters in Taipa and Coloane. Although there are more shelters in the Macao Peninsula than in Taipa and Coloane combined, the shelters in the Macao Peninsula can house only 10,728 people in total. Meanwhile, the shelters in Taipa and Coloane can house 13,117 people in total.   Figure 4b shows the locations of 6 shelters and 390 buildings in Taipa and Coloane. The sizes of shelter symbols correspond to their capacities, which are listed in Table 1. These 17 shelters were constructed by Macao's government and are operated during typhoon disasters. The residents of the buildings need to be allocated to their designated shelters when the inundation depth is more than 1 m. The figure shows that most of the buildings are in the Macao Peninsula, while the buildings in Taipa and Coloane are much less dense. Taipa Refugee center 28 13 Taipa and Coloane Social Service Center 28 14 Olympic Sports Center (Taipa) 1883   Taipa Refugee center 28 13 Taipa and Coloane Social Service Center 28 14 Olympic Sports Center (Taipa) 1883 15 Macao East Asian Games Dome (Coloane) 10,878 16 Macao Academy of Public Security Forces 50 17 Training Base of Macao Flying Eagle Association 250

Results and Discussion
In the case study area of Macao, this section shows the inundation result for Typhoon Mangkhut 2018, using a flooding estimation model. The affected buildings are determined and inputted into the mathematical model. In this way, the allocation scheme for affected buildings is generated.

Affected Building
The resolution of the DEM data was 5 m, with sea level data observed by the Barra station of Macao every five minutes. This data was applied in the 2D version of Flood Map [30,31]

Results and Discussion
In the case study area of Macao, this section shows the inundation result for Typhoon Mangkhut 2018, using a flooding estimation model. The affected buildings are determined and inputted into the mathematical model. In this way, the allocation scheme for affected buildings is generated.

Affected Building
The resolution of the DEM data was 5 m, with sea level data observed by the Barra station of Macao every five minutes. This data was applied in the 2D version of Flood Map [30,31] Figure 5; Figure 6 show the changes in inundation range and depth per four hours for the Macao Peninsula and for Taipa and Coloane. It can be seen that the inundation depth of the Macao Peninsula is deeper than that of Taipa and Coloane.
In Figure 5, the inundation depth is more than one meter in most of the Macao Peninsula, except the central area. The flooding is most serious after 4 h. Similarly, in Figure 6, flooding in the northern part of Taipa is more visible than in the south. It can be seen that the west coast is the most affected area in Taipa and Coloane. Comparing Figures 5 and 6, it is clear that until 4 h after the start of the flood, the inundation depth in the west coastal area of the Macao Peninsula was more than 1 m. Meanwhile, in the west coastal area of Taipa, it was less than one meter.
For these shelters, the government will make additional preparation to ensure evacuee safety. Moreover, as evacuees cannot walk on paths where the inundation depth is more than one meter, all of the evacuees should be evacuated as soon as possible upon receiving warnings, before the typhoon arrives. Compared with other typhoon location-allocation studies, the results in this paper could help the government to obtain knowledge about high risk areas and the people who need to be evacuated. Based on this information, the evacuees could be estimated.    In this paper, it is assumed that the people in buildings where the inundation depth reaches more than one meter during the inundation process should be evacuated. Figure 7 shows the maximum inundation range where the depth reaches more than one meter. The buildings in this range were Sustainability 2020, 12, 3308 8 of 13 selected from among all of the buildings and are shown in the figure. There were a total of 2,258 buildings from which evacuation to shelters was required. Notably, Macao Peninsula was much more affected than Taipa and Coloane in terms of the number of evacuated buildings. It can also be seen that some shelters are in the area where the inundation depth is more than one meter. For these shelters, the government will make additional preparation to ensure evacuee safety. Moreover, as evacuees cannot walk on paths where the inundation depth is more than one meter, all of the evacuees should be evacuated as soon as possible upon receiving warnings, before the typhoon arrives. Compared with other typhoon location-allocation studies, the results in this paper could help the government to obtain knowledge about high risk areas and the people who need to be evacuated. Based on this information, the evacuees could be estimated.

Evacuee Allocation
With the solution of the mathematical model obtained using GIS ModelBuilder, the scheme for typhoon shelter selection and evacuees allocation was obtained. This is shown in Figure 8. The buildings and their assigned shelters are connected with different colored lines, according to the shelters allocated. Evacuees in these buildings should be allocated to shelters, as the inundation depth is more than one meter. It can be seen that all of the shelters are selected in the Macao Peninsula, while two shelters are not selected on the offshore island (Taipa and Coloane). The two unselected shelters are Taipa and Coloane Social Service Center, and Macao East Asian Games Dome (Coloane).

Evacuee Allocation
With the solution of the mathematical model obtained using GIS ModelBuilder, the scheme for typhoon shelter selection and evacuees allocation was obtained. This is shown in Figure 8. The buildings and their assigned shelters are connected with different colored lines, according to the shelters allocated. Evacuees in these buildings should be allocated to shelters, as the inundation depth is more than one meter. It can be seen that all of the shelters are selected in the Macao Peninsula, while two shelters are not selected on the offshore island (Taipa and Coloane). The two unselected shelters are Taipa and Coloane Social Service Center, and Macao East Asian Games Dome (Coloane). This does not imply that these shelters are suboptimal in terms of construction and capacity. In fact, there are many buildings (where the inundation depth is not > 1 m) near shelter 13 (Taipa and Coloane Social Service Center). If these buildings are affected by inundation, the evacuees can quickly reach shelter 13. Shelter 15 is the Macao East Asian Games Dome (Coloane), whose capacity is the greatest of all 17 shelters. Thus, these two shelters can potentially be of service when the other shelters' capacities are insufficient.
In Figure 8, it can also be observed that there are some intersecting lines. This is because of the paths' condition and paths' connectivity. For example, assume that the straight distance between building A and shelter A is greater than between building A and shelter B. However, the route, which comprises different paths between building A and shelter A, could be shorter than that between building A and shelter B. For the Macao Peninsula, the average evacuation distance is 573 m, while the maximum evacuation distance is 2081 m and the minimum evacuation distance is 53 m. For Taipa and Coloane, the average evacuation distance is 1139 m, while the maximum evacuation distance is 2724 m and the minimum evacuation distance is 28 m. building A and shelter A is greater than between building A and shelter B. However, the route, which comprises different paths between building A and shelter A, could be shorter than that between building A and shelter B. For the Macao Peninsula, the average evacuation distance is 573 m, while the maximum evacuation distance is 2081 m and the minimum evacuation distance is 53 m. For Taipa and Coloane, the average evacuation distance is 1139 m, while the maximum evacuation distance is 2724 m and the minimum evacuation distance is 28 m.  Table 2 presents the numbers of buildings, people, and evacuees allocated to each shelter. Notably, the number of people assumes that all of the people in the buildings go to shelters. However, typhoons are unlike other disasters, such as earthquakes. Well-built and tall buildings tend not be affected much by inundation resulting from typhoon disasters. Thus, not all of the people in buildings will go to shelters. According to the simulation results for the inundation, and experience, only residents living on the ground floor should be allocated to shelters. According to data provided by Statistics and Census Bureau, the average height of the buildings is 28 m, which equals around 9 floors. It can thus be estimated that the population living on a building's ground floor is only around one-ninth of the total population. Table 2 indicates that only shelter 5, 7, 9, 13, 14 and 15 can accommodate the evacuees allocated. However, there are many other ways to influence evacuee's behavior, so that some people may go to other cities for a short period prior to the typhoon. Therefore, there are usually not so many people to house in the shelters. Even so, the government still needs to enhance the shelters' capacities.   Table 2 presents the numbers of buildings, people, and evacuees allocated to each shelter. Notably, the number of people assumes that all of the people in the buildings go to shelters. However, typhoons are unlike other disasters, such as earthquakes. Well-built and tall buildings tend not be affected much by inundation resulting from typhoon disasters. Thus, not all of the people in buildings will go to shelters. According to the simulation results for the inundation, and experience, only residents living on the ground floor should be allocated to shelters. According to data provided by Statistics and Census Bureau, the average height of the buildings is 28 m, which equals around 9 floors. It can thus be estimated that the population living on a building's ground floor is only around one-ninth of the total population. Table 2 indicates that only shelter 5, 7, 9, 13, 14 and 15 can accommodate the evacuees allocated. However, there are many other ways to influence evacuee's behavior, so that some people may go to other cities for a short period prior to the typhoon. Therefore, there are usually not so many people to house in the shelters. Even so, the government still needs to enhance the shelters' capacities.

Discussion
As mentioned, this paper considers the inundation caused by a typhoon; this has usually been overlooked in other literature. Inundation was modeled using DEM data with high resolution (5 m). In general, the resolution of DEM data is higher and the modeling is better, but the modeling time is longer. Overall, the 5 m resolution of DEM data is better for modeling inundation. Using the estimated inundation, the number of evacuees in each building was estimated. The number of evacuees allocated to each of the 17 Macao government shelters was then compared with the shelter capacities. In the case of Typhoon Mangkhut 2018, residents in buildings where inundation was less than 1 m did not evacuate. Additionally, residents above the ground floor were safe in their homes even if their buildings were affected, as the maximum inundation depth was 2.01 m. Thus, the government could carry out rescue and evacuation more accurately, and emergency resources could be used in a more effective way. For Typhoon Mangkhut 2018, 11 shelters in the Macao Peninsula and 4 shelters in Taipa and Coloane were selected. The other two shelters were not selected by evacuees in this paper, as they were relatively far away from the affected buildings. However, these shelters remain important; indeed, they will be essential for accommodating surge evacuees if other shelters are full. The selected shelters would be crowded if all evacuees chose to come, especially shelter 1. Most of the shelters' capacities are smaller than the number of evacuees allocated to them. Even though the government could help some evacuees to go to shelter 13 and 15, combined capacity of all 17 shelters is still inadequate. In this paper, the capacity constraint was not considered in the mathematical model as it has been done in other studies. These studies assumed that if a shelter could not accommodate all members of a community, it would not be selected by the community, even if it were nearer than other larger shelters. To achieve system-optimal behavior, some evacuees must go further to shelters; however, people might exhibit self-centered behavior and to go to nearer shelters overlooking the shelter's capacity. Therefore, in this paper, the capacity constraint was not considered in the mathematical model; instead, the shelters' capacities were considered after the result of the mathematical model was obtained. The number of evacuees allocated to each shelter and the capacities of each shelter are shown in Table 3. This indicates that most shelters cannot satisfy the needs of the evacuees near them, especially shelter 1. On the other hand, the number of evacuees per building was estimated based on average building height because of the lack of other data. In reality, the ground floors of some buildings are used as stores, which was not considered in this paper. Furthermore, all of the evacuees were assumed to go to shelters. In reality, they could opt for other safe places, such as their friends' houses. Thus, the number of evacuees allocated to shelters obtained from the model may be greater than the true number of evacuees. However, as shown in Table 3, some shelters have a much lower capacity than the number of evacuees allocated to them. Therefore, it is suggested that the government should consider expanding existing shelters or build more shelters on suitable sites. For example, the capacities of shelter 1, 6 and 12 are much too low for the number of evacuees allocated to them, so they need to be expanded. Although the capacities of shelters 16 and 17 are also insufficient according to our results, the gaps are much smaller; therefore, they can maintain their status, as some evacuees will not choose to go to shelters. By comparing the capacity of each shelter with the number of evacuees allocated to it, the government can achieve better planning of shelter construction and improvement.
In future work, the number of evacuees can be estimated more realistically, by considering the usage of some buildings' ground floors. Moreover, surveys could help determine the number of evacuees willing to go to shelters. Another assumption was that the objective of the model was to minimize total evacuation distance rather than evacuation time. Walking speeds and the proportions of people of different ages were not considered in this paper. This factor also should be considered in future work.

Conclusions
The aim of this paper was to present work on allocating affected people to shelters constructed by the government. The number of affected people was determined using the flooding simulation model. Our findings can inform future work on affected people allocation. They can also help the government to plan the construction of new shelters and improvement of existing shelters.
In this paper, a mathematical model designed to minimize total evacuation distance was proposed, with a distance constraint. A flooding simulation model was used to simulate the inundation process caused by Typhoon Mangkhut 2018; the affected buildings from which residents should be evacuated were then identified. In the model, all of the existing shelters in the Macao peninsula were occupied, while two shelters in Taipa and Coloane were not. Moreover, considering the average building height, the number of evacuees was estimated. The results show if all affected people chose to go to shelters, the government would need to improve the capacities of the 17 existing shelters.
In further work, a number of improvements could be made. Firstly, the number of people that live in the ground floor of each building should be determined, and the number of affect people can be estimated in more detail. In this way, the capacities of shelters can be assessed to help the government improve current shelters. Additionally, the willingness of affected people to evacuate should be considered, since some people will not go to shelters even if they are affected by a typhoon. Furthermore, the influence of older people and children on the evacuation process can be analyzed. Finally, dynamic changes in road networks and buildings due to damage caused by a typhoon can be considered.