Methodology and Application of Spatial Vulnerability Assessment for Evacuation Shelters in Disaster Planning

Evacuation shelters are the most important means for safeguarding people in hazardous areas and situations, and thus minimizing losses, particularly those due to a disaster. Therefore, evacuation shelter assignment and evacuation planning are some of the critical factors for reducing vulnerability and increasing resilience in disaster risk reduction. However, an imbalance of shelter distribution and spatial heterogeneity of a population are the critical issues limiting the accessibility of evacuation shelters in real situations. In this study, we propose a methodology for spatial assessment to reduce vulnerability and evaluate the spatial distribution of both shelter demand and resources, considering spatial accessibility. The method was applied to the case study of Mabi, in the context of a disaster caused by the 2018 flooding. We applied this approach to evaluate the area and identified the vulnerability of the evacuation shelters and the residents. The proposed method revealed that 54.55% of the designated evacuation shelters and 59% of the total population were physically vulnerable to the flood. The results highlight, using GIS maps, that the total shelter capacity was significantly decreased to 43.86%. The outcome assessment addressed specific vulnerable shelters and the imbalance between the demand for and resources of each shelter. Accordingly, this study provides practical information and a valuable reference for supporting local governments and stakeholders to improve future disaster planning, prevention, and preparedness.


Introduction
Globally, disasters are forecasted to occur more frequently and to be more severe in the near future, due to climate change, increasingly unplanned urbanization, and the rising concentration of people and assets in hazardous areas [1,2]. Over the past decade, extreme climate events have caused an enormous number of fatalities and a significant amount of economic damage. At least 15,490 natural disasters have been reported, since 1990, in the Emergency Management Database (EMDAT-CRED). These disasters, taken together, have caused more than 32.6 million deaths, affected over eight billion people, and caused a total damage of at least USD 3.5 trillion, worldwide [3]. To address these severe issues, in 2015, three global policy frameworks, i.e., the Sendai Framework for Disaster Risk Reduction (DRR), the Sustainable Development Goals (SDGs), and the Paris Agreement on Climate Change, were established by the United Nations. Among these three frameworks, vulnerability reduction is shelters have commonly been included in the aforementioned studies. However, so far, discussions on demand for shelter capacity have not fit the actual situation due to limited data analysis that capture this in reality [19]. Numerous studies have assumed that 25% of residents would evacuate to emergency shelters; therefore, they have treated it as an input datum using a simple formula (e.g., 25% × the residential population) and estimated the number of people and evaluated the deficiency in terms of shelter capacity [20][21][22]. However, in reality, the evacuation of residents to an emergency shelter is impacted by several factors.
Spatial distribution is one of the significant issues that requires further discussion in relation to the estimation of demand for and the limited supply of emergency shelters. Considering a population in different scenarios, particularly when a disaster occurs, can aid in the estimation of the affected population in terms of evacuation planning and response. Previous studies and proposed models have considered spatial population distribution through urban and community analyses, using GIS, to analyze emergency shelter demand [23][24][25]. Although most research has placed more emphasis on the spatial distribution of a population, only limited research has taken the spatial distribution of existing emergency shelters into consideration [26]. However, both factors of distribution are critical issues. They influence the availability of shelters, and also the accessibility of shelters to a population. An imbalance between the spatial distribution of people and shelters could lead to and become a vulnerability that must be assessed in disaster and evacuation planning.
Furthermore, additional aspects that notably influence the accessibility of emergency shelters are the evacuation route conditions and distance. These are important factors influencing a decision to evacuate and the selection of an evacuation destination for people [27]. Routes inundated by floods or covered by collapsed buildings due to an earthquake decrease the efficiency of an evacuation and impact the travel cost for evacuees. Moreover, traveling along an evacuation route to reach an emergency shelter is considered to be a challenging task, especially for individuals who belong to a vulnerable population group, for example, the elderly [28][29][30]. Previous studies have also investigated the relationship between gender-specific variations in terms of socioeconomic status and evacuation, revealing that women were more likely to evacuate than men [31,32].
To estimate the accessibility of facilities, according to previous literature, there are two common methods for measuring approachability, i.e., a distance-based function such as Euclidean distance or network distance, and a time-based function. According to these methods, several models and optimizations have been developed to deliver the best all-purpose evacuation route by minimizing the total evacuation distance and potential devastating damages caused by a disaster [7,24,33,34]. Ye et al. proposed a mathematical formulation to minimize distance by assuming that shelters become available at different periods [24]. Another mathematic approach was proposed by Zhao et al. which considered the different degrees of house damage to identify the optimal locations for evacuation shelters [7]. However, most simulation models have usually been developed on the basis of different demand circumstances and criteria among the various available candidates. Such methods limit administrators with respect to applying valuable unknown criteria in reality. In addition, most methods have concentrated more on developing a new algorithm to simulate evacuation procedures, rather than considering the feasibility of evacuating residents efficiently in practice. Consequently, several studies have revealed that it is not yet clear how the accessibility of disaster sheltering should be evaluated in practice [10,35].
In this study, we propose a novel methodology for the evaluation of accessibility, integrating the spatial distribution of both a population and the evacuation shelters for assessing and reducing the vulnerability of disaster shelter planning. The objectives of the study are as follows: First, to investigate and evaluate the spatial distributions in the case study area and estimate the deficiency of an evacuation shelter based on the spatial evacuation demand; second, to analyze the spatial accessibility of the evacuation shelters based on the evacuation route conditions and distance, by comparing the real road network of the normal situation with the real disaster situation; and last, to identify the vulnerability of existing evacuation shelters and populations in terms of disaster risk reduction. As a systematic Sustainability 2020, 12, 7355 4 of 22 approach for analysis of spatial availability, integrating all aspects of spatial distribution and spatial accessibility, including distance, and evaluating the physical and social vulnerabilities, the proposed methodology is described in Section 2. In Section 3, we demonstrate the evidentiary foundation through a case study using GIS technology. A discussion based on the algorithm, suggestions for future studies, and the summary are presented in Section 4. The overall target of this study is to support decision makers by providing a useful reference for developing policies and strategies to reduce disaster vulnerability.

Methodology
According to the current literature, effective shelter assignment and proficient evacuation planning require an understanding of the suitability of the shelter location and the distribution between the shelter and the population. GIS technology is a significant tool that supports researchers for understanding critical information and visualizing through mapping. However, to analyze and assess the vulnerability of disaster shelter planning, the accessibility of a shelter based on an imbalance in the distribution of the shelter demand and resources is a critical issue that requires further discussion.
Our assumption is that vulnerability during a disaster is caused by the spatial heterogeneity of the imbalance between demand and capacity resources, which could be concealed in summarized statistics and due to limited accessibility. To examine this hypothesis, we structured the method as follows ( Figure 1): (1) analysis of the spatial distribution of shelter planning resources and population to assess the spatial availability of shelters; (2) estimation of the accessibility based on the spatial availability and distribution; (3) estimation of the shelter demand to assess the shelter deficiency based on the spatial accessibility of the normal situation as compared with the disaster situation; and (4) analysis of the vulnerability of the evacuation shelters in terms of vulnerability reduction in disaster preparedness and management. We describe the details of the methodology in the following sections.
Sustainability 2020, 12, x FOR PEER REVIEW 4 of 23 reduction. As a systematic approach for analysis of spatial availability, integrating all aspects of spatial distribution and spatial accessibility, including distance, and evaluating the physical and social vulnerabilities, the proposed methodology is described in Section 2. In Section 3, we demonstrate the evidentiary foundation through a case study using GIS technology. A discussion based on the algorithm, suggestions for future studies, and the summary are presented in Section 4. The overall target of this study is to support decision makers by providing a useful reference for developing policies and strategies to reduce disaster vulnerability.

Methodology
According to the current literature, effective shelter assignment and proficient evacuation planning require an understanding of the suitability of the shelter location and the distribution between the shelter and the population. GIS technology is a significant tool that supports researchers for understanding critical information and visualizing through mapping. However, to analyze and assess the vulnerability of disaster shelter planning, the accessibility of a shelter based on an imbalance in the distribution of the shelter demand and resources is a critical issue that requires further discussion.
Our assumption is that vulnerability during a disaster is caused by the spatial heterogeneity of the imbalance between demand and capacity resources, which could be concealed in summarized statistics and due to limited accessibility. To examine this hypothesis, we structured the method as follows ( Figure 1): (1) analysis of the spatial distribution of shelter planning resources and population to assess the spatial availability of shelters; (2) estimation of the accessibility based on the spatial availability and distribution; (3) estimation of the shelter demand to assess the shelter deficiency based on the spatial accessibility of the normal situation as compared with the disaster situation; and (4) analysis of the vulnerability of the evacuation shelters in terms of vulnerability reduction in disaster preparedness and management. We describe the details of the methodology in the following sections.

Spatial Distribution of the Shelter Demand and Resources
Understanding the size, density, and spatial distribution of both a population and designated or existing evacuation shelters is essential to disaster preparation planning and resource allocation. Understanding the size, density, and spatial distribution of both a population and designated or existing evacuation shelters is essential to disaster preparation planning and resource allocation. Census data are generally used with specific geographic coverage areas based on units of the census tract scale, such as country-, region-, or municipal-scale units. According to the administrative units, they have a low spatial resolution and do not reflect the dynamic distribution in a practical situation. In addition, the actual population density can fluctuate significantly within an area. The same value for population density could be a feature of an urban area and also a mountain area, where the resident is normally a part of the administrative area [36].
Therefore, in this study, an area-based mapping method, in which the original administrative unit is divided into smaller spatial units, was implemented to obtain the information and estimate the spatial distribution of the case study area. Within the boundaries, the case study area of the city polygon, which represented the shape and location, contained various data features, and the rectangular cell girds were employed to capture the data at the municipal polygon scale into a smaller spatial grid unit, as illustrated in Figure 2. Census data are generally used with specific geographic coverage areas based on units of the census tract scale, such as country-, region-, or municipal-scale units. According to the administrative units, they have a low spatial resolution and do not reflect the dynamic distribution in a practical situation. In addition, the actual population density can fluctuate significantly within an area. The same value for population density could be a feature of an urban area and also a mountain area, where the resident is normally a part of the administrative area [36]. Therefore, in this study, an area-based mapping method, in which the original administrative unit is divided into smaller spatial units, was implemented to obtain the information and estimate the spatial distribution of the case study area. Within the boundaries, the case study area of the city polygon, which represented the shape and location, contained various data features, and the rectangular cell girds were employed to capture the data at the municipal polygon scale into a smaller spatial grid unit, as illustrated in Figure 2. Because the development of GIS technology allows us to use digital maps that recombine and process data, this technique was generated and implemented to be combined with and integrated into a spatial GIS environmental system.

Spatial Accessibility Estimation
Regarding an efficient estimation of the evacuation shelter deficiency, the accessibility based on the spatial distribution of both the population and shelters, including the evacuation distance and route, were taken into account in this research. However, several studies have estimated that 25 percent of the total number of census data could evacuate to emergency shelters, and therefore have applied this information as the input data to estimate shelter demand. However, in a real situation of a disaster, where and how people can reach shelters are also critical issues.
In the literature, the concept of spatial accessibility is described as a cumulative opportunity that defines how many potential opportunities there are or how many destinations can be reached from a particular origin within a threshold travel time or distance [37,38]. Thus, previous studies have often applied the Euclidean distance, which focuses on the simple radius as a linear distance from shelters, to measure the accessibility. However, this method is still, impractical for community disaster evacuation in practice. In addition, based on the post-disaster reports and surveys, evacuation route and distance are the significant factors impacting decisions regarding evacuation and the selection of an evacuation route. Additionally, previous research and guidelines have suggested that the desired distance of an evacuation point should be appropriate for residents, particularly vulnerable populations.
Therefore, this study applied the distance-based approach of an actual road network to analyze the actual accessibility of shelters for a population. To present an origin and perform an accessibility Because the development of GIS technology allows us to use digital maps that recombine and process data, this technique was generated and implemented to be combined with and integrated into a spatial GIS environmental system.

Spatial Accessibility Estimation
Regarding an efficient estimation of the evacuation shelter deficiency, the accessibility based on the spatial distribution of both the population and shelters, including the evacuation distance and route, were taken into account in this research. However, several studies have estimated that 25 percent of the total number of census data could evacuate to emergency shelters, and therefore have applied this information as the input data to estimate shelter demand. However, in a real situation of a disaster, where and how people can reach shelters are also critical issues.
In the literature, the concept of spatial accessibility is described as a cumulative opportunity that defines how many potential opportunities there are or how many destinations can be reached from a particular origin within a threshold travel time or distance [37,38]. Thus, previous studies have often applied the Euclidean distance, which focuses on the simple radius as a linear distance from shelters, to measure the accessibility. However, this method is still, impractical for community disaster evacuation in practice. In addition, based on the post-disaster reports and surveys, evacuation route and distance are the significant factors impacting decisions regarding evacuation and the selection of an evacuation route. Additionally, previous research and guidelines have suggested that the desired distance of an evacuation point should be appropriate for residents, particularly vulnerable populations.
Therefore, this study applied the distance-based approach of an actual road network to analyze the actual accessibility of shelters for a population. To present an origin and perform an accessibility Sustainability 2020, 12, 7355 6 of 22 calculation, in this study, the centroid point feature or the center of an area from the grid units was generated and represented the population points within each area. At this stage, a road network analysis was conducted, using the GIS platform to generate evacuation routes and the estimated evacuation cost of navigation for each population point. Accordingly, several evacuation guidelines and studies have suggested the creation of suitable shelters at the nearest location as the best solution for all evacuees [39,40]. Hence, the approach in this study employed a "closest facility" analysis to create the travel route between the population points and designated evacuation sites. This application was based on Dijkstra's algorithm for finding the shortest paths by implementing two types of path-finding algorithms [41]. The first type is the exact shortest path. The algorithm searches for the distance from the starting point to every other vertex, until it reaches the destination point. The second is a hierarchical path solver, which has a faster performance and gives the shortest path possible.
As a starting point, by implementing the "closest-facility" analysis, the population points were applied as the incident features in the GIS application. The location points of the designated evacuation shelters were factored in as facility features using the closest facility analysis technique. Furthermore, according to previous literature, the use of cars or other vehicles was not recommended for evacuation purposes, because they represented a threat to the physical safety of evacuees [42,43]. In addition, during a disaster, in heavy rain or when it is dark, older people generally do not have the confidence to drive; thus, shelters with a convenient location and within walking distance could give the elderly greater encouragement and support for evacuating [44]. Thus, it was assumed, in this investigation, that walking would be the primary method for accessing the closest shelters. The travel mode of walking was applied, with a fixed walking speed of 1.4 m/s, which was considered to be the average rate of a normal adult. In the analysis of the case study area, a two-kilometer network distance was considered to be appropriate for the accessibility measurement, according to the guidelines of each country [45], which suggested that the desirable distance of evacuation on foot should be less than 2 km in all cases.

Shelter Demand Estimation and Deficiency Assessment
While shelters and populations can be mapped within spatial contexts, accessibility, which is impacted by evacuation decisions, is strongly constrained by the road network and other conditions. Therefore, in this step, the travel cost along the evacuation route, which is the distance and duration of the route, was investigated to measure the accessibility of the evacuation shelters. The road network distance between each point, representing residential units and shelters, was computed, including the estimated travel time. Then, each evacuation shelter, represented by a shelter ID, was distinguished from all the population access points, denoted by population IDs, within the threshold distance appropriate for the case study area. As the output of this process, a table of each evacuation shelter, with the population IDs, the size of the population within each grid area, and the travel distance and time was provided (Table 1). Then, the total access area for each shelter was summarized, including the total estimated population demand, which was also calculated by merging the total size of the accessed population points that can reach the evacuation shelter using the following formula: where D i is the estimated population demand for shelter, i, and P j,i is the size of the population of each area, j, that can access shelter, i.
To evaluate the capacity of each shelter, the ratio of the accommodation capacity was also examined by dividing the designated capacity of each shelter by the estimated population demand using the following equation: where R i is the ratio of the capacity of shelter i and C i is the designated capacity of shelter i.

Vulnerability Analysis
Since the concept of vulnerability varies among people due to their different circumstances, taking into consideration previous research, the main concepts were classified into two dimensions, i.e., physical and social perspectives. In engineering science research, the physical aspect is usually considered to be a function of a given process intensity towards physical structures. Therefore, physical vulnerability is often associated with the susceptibility of the element at risk, which comprises its geography, location, place, settlement patterns, and physical structures [46][47][48]. Conversely, the social dimension is complex and multifaceted and typically describes the socioeconomic and demographic factors that affect the resilience of communities. It describes a population's ability to prepare for, respond to, cope with, recover from, and ultimately adapt to environmental hazards [49].
Therefore, for this stage, to analyze the vulnerability of the evacuation shelters in the case study area, the location of previous hazardous areas was taken into account and applied to prevent location from being a barrier feature in the evaluation of the accessibility of the shelter during an actual disaster event in the case study area. Then, the new evacuation route to the closest evacuation shelter was once more generated, and the travel cost was calculated for each evacuation shelter. To evaluate the vulnerability of each shelter, with consideration of a previous actual disaster in the case study area, this study investigated the physical vulnerability in terms of the location, capacity, and travel cost impacted by the disaster. Furthermore, the total size of the access area for each shelter and the estimated population demand after the disaster were also calculated using the following equation: where D f i is the estimated population demand for shelter i after the disaster, and P f j,i is the size of the population for each area j that accesses shelter i after the disaster.

Case Study Area
The case study area is the town of Mabi, which is located in Kurashiki City, Okayama prefecture, in the western region of Japan, as shown in Figure 3. Japan is one of several countries that is vulnerable to natural disasters due to its topography. Japan has experienced numerous types of disaster, and recently, in 2018, it suffered the worst flooding in decades in the western part of the country. The damage affected many people, killing 225 residents and destroying their houses, properties, and livelihoods. In particular, more than 8 million people were required to evacuate across 23 prefectures [50][51][52]. This disaster caused enormous damage and imposed massive operational difficulties on the local government.
The town of Mabi, one of the most severely affected areas, is geographically enclosed by two large rivers, namely, Oda river and Takahashi river, which overflowed and covered almost one-third of the area of the town, causing 51 deaths [53,54]. The area of Mabi is 44.12 km 2 , with a total population of 22,594 and a population density of 512.1 people/km 2 . Historically, Mabi has repeatedly experienced devastating floods due to the overflow of the Oda river; however, the water level has never exceeded 50 cm since the 1970s [55]. Nonetheless, under Japanese law, hazardous maps for disasters are required by the city authority, which considers the area to be high risk. However, by the time the official hazardous map was distributed in 2016, most of the high-risk flooding areas were already occupied Sustainability 2020, 12, 7355 8 of 22 due to suburbanization [56]. Due to the heavy rain from the end of June to the beginning of July 2018, the dikes along the Oda and Takahashi rivers were breached on 6 July 2018. The intensity of the 2018 heavy rain was extreme, because it did not last for a short period of one to six hours but continued for an extensive period of more than 24 h. For that reason, a flood was triggered, and the town of Mabi, therefore, was extensively inundated. However, the city office had been preparing for an emergency since the morning of 5 July, and the first evacuation was ordered for the high-risk areas around the Oda river. According to the report [55], it was revealed that the residents did not start evacuating to a safe place in time; therefore, 51 residents died in the water, and 40 of these victims drowned in their own houses. Flood depths exceeding 5 m were found all over areas of Mabi [54]. The flood lasted for 5 days, and the flood depths were very high (exceeding 5 m all over areas of Mabi). More than 4600 families were affected and required immediate temporary accommodations, and numerous residents had to continue to stay at shelters until the end of August. In Japan, the mayors of the municipalities were also required to designate both public and private facilities (e.g., schools, junior high schools, and community centers) as "designated shelters" by the Japanese law. Furthermore, to ensure a smooth and prompt evacuation, these facilities were also designated as "emergency evacuation sites" and "shelters" for temporary stays. Due to the catastrophic loss due to the flood in 2018, an evaluation of vulnerability reduction in an effort to "Build Back Better", i.e., enhance disaster preparedness for an effective response within the Sendai Framework, is urgently needed.

GIS-Related Data
The data used in this research and analysis of both spatial and non-spatial data were obtained from several sources. Detailed information about the data, sources, and application in this study is presented in Table 2.
Geographic data on census statistics were obtained from the government statistic portal site of Japan [56]. GIS boundary data for the case study area at the scale of the city area and the 5th-level grids, as the smallest polygon mesh of approximately 250 m grids (7.5 s latitude and 11.25 s longitude), were collected in a shape-file format. In addition, census data operating with this 250 m grid resolution were also provided for every 5 years. To investigate the base statistics, the latest census data, for 2015, from the e-stat portal of Okayama prefecture were obtained as a CSV file for processing.
People older than 65 years are usually included in the definition of vulnerable groups who are often frail and require regular and continuous health care [57,58]. This case study area, Mabi, has essentially been driven to reduce disaster risk by a considerable increase in the elderly population; therefore, the population of this vulnerable group of residents aged over 65 years old was also investigated in this study.  [59]. The list of designated evacuation sites was obtained as a CSV file from the latest update published in March 2018, by the Kurashiki city authority's department.
Flooding data on the 2018 disaster event area were acquired from the Geospatial Information Authority of Japan (GSI), who published estimated inundation maps from 7 July 2018, for investigation and comparison in this study. These collected data were based on various pieces of information, such as the image on July 7, which confirmed the endpoint of the submerge range and estimated the flooded surface from the height of that point using elevation data [60].

Results
In this study, ArcGIS 10.7.1 was used to process, map, and analyze the data. The point location of all emergency centers, i.e., the latitude and longitude (X and Y) coordinates of these designated evacuation sites, were obtained from Google Maps and added to the ArcGIS software.
The census data were merged with areas with the same key code into a 250 m polygon mesh grid layer. Using the Mabi town boundary layer, this grid layer was clipped and joined with the Mabi town feature layer. Within these total grid areas, 337 grids contained a population of at least one resident. The rest (413 grid areas or 45% of the study area) comprised zero population data, which were areas in which mountains, roads, lakes, or rivers were located, as shown in Figure 3.

Spatial Distribution
In Mabi, more than 30% of the total population is considered to be vulnerable, and 31.68% and 13.06% of the elderly population are over 65 years old and 75 years old, respectively (Table 3). Additionally, within the vulnerable group of elderly people over 65 years old, 54.39% (3893 residents) are female seniors. Previous studies have revealed that the elderly tend to have a higher mortality rate during a disaster [61]. Focusing on age is significant, because the elderly problem is more critical in Japan than in other countries. For the overall population, 22 facilities were identified as the designated evacuation shelters in Mabi. Table 3 illustrates that all of these designated evacuation shelters have the capacity to accommodate 3110 people, which are available for 13.76% of the total population.
Technical planning according to previous disaster research on emergency shelter planning [18,19] has indicated that 25% of both total populations and vulnerable residents would use public shelters in the case of an evacuation. In the case study area, a shelter capacity of 13.76% for the total population is insufficient. However, for the vulnerable group of residents aged over 65, the shelter capacity of 3110 is accessible for 43.45% of them, which is more than 25%, thus meeting the technical planning standard. However, according to a report on the 2018 flood disaster, many residents could not access the evacuation shelters due to the condition of the evacuation route and a shortage and poor distribution of shelters [55]. Figure 4 illustrates the spatial distribution of both the residential population and the designated evacuation shelters in the study area. The spatial distribution of the population ranges from 1 to 409, with a mean value of the population of 30.12 per grid area. The population density of each grid area is classified by a natural break and shown in a different color, from light yellow to dark red, as shown in Figure 4. The high population density (red) areas are mostly located in the northeastern part of the case study area, while the western part of the town mostly has a low population density (yellow) and areas with no population, such as the areas with mountains, lakes, or rivers. However, the shelter capacity distribution ranges from 20 to 340, with a mean and median value of 141.36 and 150 accommodations, respectively. The size of each shelter capacity is illustrated and distinguished from the others by using blue symbols with different scales. The figure also shows the location of the designated evacuation shelters, where the shelters with greater capacities (larger blue symbols) are mostly located in the center and northeast of the town.

Accessibility Assessment
In the case study area, during the flood situation in 2018, an evacuation was ordered by the local authorities for the town of Mabi, including areas with no flooding. Accordingly, the study concentrated on the accessibility and shortage of evacuation shelters in terms of the evacuation demand, resources, and evacuation route distance. Therefore, the demand of the population in the areas with no flooding was also investigated to evaluate the availability and accessibility of the evacuation shelter. Figure 5 shows the evacuation routes, along with the actual road network, to the closest designated evacuation shelters, from every population point in each grid area. To determine the accessibility for each population point and conduct a road network analysis, 303 of 337 population points, which constitutes 89.91% of them, are located within a walking distance of 2 km from the designated evacuation shelters, as shown in the brown evacuation route in the figure. However, the blue evacuation routes in the figure illustrate the remaining 34 population points, which are located within a walking distance of 7 km from the evacuation shelters, and 23, 3, 3, and 5 population points are estimated to be within 3, 4, 6, and 7 km, respectively. These long-distance evacuation routes are located in the western part of the town, and most areas are covered by mountains and forests. Particularly at the edge of the boundary in the southwestern part, the evacuation route exceeds the boundary of the case study area due to the unavailability of a road network through the mountain zone within the border.

Accessibility Assessment
In the case study area, during the flood situation in 2018, an evacuation was ordered by the local authorities for the town of Mabi, including areas with no flooding. Accordingly, the study concentrated on the accessibility and shortage of evacuation shelters in terms of the evacuation demand, resources, and evacuation route distance. Therefore, the demand of the population in the areas with no flooding was also investigated to evaluate the availability and accessibility of the evacuation shelter. Figure 5 shows the evacuation routes, along with the actual road network, to the closest designated evacuation shelters, from every population point in each grid area. To determine the accessibility for each population point and conduct a road network analysis, 303 of 337 population points, which constitutes 89.91% of them, are located within a walking distance of 2 km from the designated evacuation shelters, as shown in the brown evacuation route in the figure. However, the blue evacuation routes in the figure illustrate the remaining 34 population points, which are located within a walking distance of 7 km from the evacuation shelters, and 23, 3, 3, and 5 population points are estimated to be within 3, 4, 6, and 7 km, respectively. These long-distance evacuation routes are located in the western part of the town, and most areas are covered by mountains and forests. Particularly at the edge of the boundary in the southwestern part, the evacuation route exceeds the boundary of the case study area due to the unavailability of a road network through the mountain zone within the border. Sustainability 2020, 12, x FOR PEER REVIEW 12 of 23  Table 4 illustrates the accessibility of each designated evacuation shelter in Mabi, showing the total area within a 2 km walking distance from the closest shelter, the points at the longest distance from a shelter, as well as the fastest and longest time durations. Shelter Number 14 is the most accessible for the population points of 29 areas, whereas shelter Number 4 has the lowest accessibility and is accessible for the population points in only two areas. The closest emergency shelters are capable of being reached from the population points in less than 8 min. However, the longest walking duration to an evacuation shelter is almost 24 min, and it is to shelter Number 14, which is also situated at the longest distance from the population points (nearly 2 km). Each designated evacuation shelter is capable of being reached from the population points in less than 8 min.   Table 4 illustrates the accessibility of each designated evacuation shelter in Mabi, showing the total area within a 2 km walking distance from the closest shelter, the points at the longest distance from a shelter, as well as the fastest and longest time durations. Shelter Number 14 is the most accessible for the population points of 29 areas, whereas shelter Number 4 has the lowest accessibility and is accessible for the population points in only two areas. The closest emergency shelters are capable of being reached from the population points in less than 8 min. However, the longest walking duration to an evacuation shelter is almost 24 min, and it is to shelter Number 14, which is also situated at the longest distance from the population points (nearly 2 km). Each designated evacuation shelter is capable of being reached from the population points in less than 8 min.   Table 5 outlines the capacity of each designated evacuation shelter in Mabi. The results for the estimated shelter demand (D i ) for the closest facilities within a 2 km walking distance using network analysis show that the estimated demand at almost every shelter massively outweighs their capacity; only the lowest estimated demand, at shelter Number 4, with an estimated access area of two population points, displays a high ratio in terms of the shelter capacity (R i ), i.e., 450%. However, the capacity ratios of the other shelters are considerably low, particularly the capacity of shelter Number 13, which was designated for 30 accommodations, while the estimated demand through the closest analysis of this shelter was considered to be 1561 residents. Table 5. The accessibility of each designated evacuation shelter, as distinguished by the total estimated population demand (D i ), the capacity ratio of each shelter (R i ), and the total estimated demand of the population over 65 years old (D i,O65 ).

Shelter Number
Shelter Capacity Number of Accessed Area D i R i D i, O65   1  160  13  2116  8%  509  2  180  12  1302  14%  388  3  300  16  1918  16%  510  4  180  2  40  450%  18  5  180  18  1633  11%  447  6  280  15  1232  23%  352  7  210  3  555  38%  238  8  270  14  296  91%  104  9  160  6  969  17%  424  10  160  28  711  23%  275  11  70  13  1463  5%  421  12  50  7  912  5%  353  13  30  17  1561  2%  560  14  40  29  1368  3%  423  15  140  10  665  21%  233  16  40  9  1283  3%  333  17  80  27  896  9%  318  18  50  10  456  11%  147  19  30  25  659  5%  240  20  340  11  1074  32%  341  21  140  6  385  36%  128  22  20  12  353  6%  122 A comparison of the designated capacity and the estimated shelter demand is demonstrated in Figure 6. The result shows that even in the western part, with a low population density, an evacuation shelter space deficit still exists in the area. 353 6% 122 A comparison of the designated capacity and the estimated shelter demand is demonstrated in Figure 6. The result shows that even in the western part, with a low population density, an evacuation shelter space deficit still exists in the area.  Figure 7 shows the evacuation route for the 2018 flooding disaster in the case study area, including the spatial distribution of the population points and designated evacuation shelters. The evacuation distances of less than 2 km and more than 2 km are illustrated in brown and dark blue, respectively. The figure also demonstrates that, of a total of 22 designated evacuation shelters in the study area, after the flood disaster, only 10 shelters remained available and could be accessed along the road network. Furthermore, in the southern and southwestern parts of Mabi, two designated shelters (shelter Numbers 9 and 17) remained available for the residents after the flood. However, Figure 6. Comparison of the shelter planning capacity and the estimated shelter demand in Mabi. Figure 7 shows the evacuation route for the 2018 flooding disaster in the case study area, including the spatial distribution of the population points and designated evacuation shelters. The evacuation distances of less than 2 km and more than 2 km are illustrated in brown and dark blue, respectively. The figure also demonstrates that, of a total of 22 designated evacuation shelters in the study area, after the flood disaster, only 10 shelters remained available and could be accessed along the road network. Furthermore, in the southern and southwestern parts of Mabi, two designated shelters (shelter Numbers 9 and 17) remained available for the residents after the flood. However, due to the long distance of the evacuation routes (more than 2 km), the residents from this area would have had difficulty accessing these shelters.

Vulnerability of Shelters
The results also show that there are 269 population points located outside the flooding area, which can access a shelter through the road network. However, 88 of these population points have a long evacuation route of more than 2 km in order to access a shelter, due to the mountainous conditions of their geographical location and the unavailability of a road connection. From all 337 population points within the study area, 68 population points are identified as vulnerable to flood, because they are located in the high hazard zone, and there is no road connection for the flood evacuation.
Because the aim of this study was to evaluate the physical vulnerability of the study area in terms of the location, distance, and accessibility of the shelters, a simulation using GIS was employed, which demonstrated that 12 of 22 (54.55%) designated evacuation shelters were physically vulnerable to flood disaster, because they were located in the flooding area. These emergency shelters were not accessible due to the unavailability of a road network and were isolated by the remaining network (Table 6).
Furthermore, Table 6 presents the remaining capacity of 10 evacuation shelters after flooding, and shows that capacity was reduced to 1340 (43.86% of the total capacity in Mabi). The demand ratio for accessibility to six of these shelters was significantly increased (shelter Numbers 2, 10, 15, 16, 20, and 22), and they were located near the flooding area ( Figure 8). The contrast between the estimated demand and the estimated demand after the flooding, including the planning shelter capacity, is shown in Figure 8.
Sustainability 2020, 12, x FOR PEER REVIEW 15 of 23 due to the long distance of the evacuation routes (more than 2 km), the residents from this area would have had difficulty accessing these shelters. The results also show that there are 269 population points located outside the flooding area, which can access a shelter through the road network. However, 88 of these population points have a long evacuation route of more than 2 km in order to access a shelter, due to the mountainous conditions of their geographical location and the unavailability of a road connection. From all 337 population points within the study area, 68 population points are identified as vulnerable to flood, because they are located in the high hazard zone, and there is no road connection for the flood evacuation.
Because the aim of this study was to evaluate the physical vulnerability of the study area in terms of the location, distance, and accessibility of the shelters, a simulation using GIS was employed, which demonstrated that 12 of 22 (54.55%) designated evacuation shelters were physically vulnerable to flood disaster, because they were located in the flooding area. These emergency shelters were not accessible due to the unavailability of a road network and were isolated by the remaining network (Table 6).
Furthermore, Table 6 presents the remaining capacity of 10 evacuation shelters after flooding, and shows that capacity was reduced to 1340 (43.86% of the total capacity in Mabi). The demand ratio for accessibility to six of these shelters was significantly increased (shelter Numbers 2,10,15,16,20,and 22), and they were located near the flooding area ( Figure 8). The contrast between the estimated demand and the estimated demand after the flooding, including the planning shelter capacity, is shown in Figure 8. Table 6. The vulnerability of the shelters.

Shelter
Before Flood After Flood   1  160  13  2116  0  0  2  180  12  1302  30  3057  3  300  16  1918  0  0  4  180  2  40  2  40  5  180  18  1633  0  0  6  280  15  1232  0  0  7  210  3  555  0  0  8  270  14  296  0  0  9  160  6  969  6  969  10  160  28  711  30  775  11  70  13  1463  0  0  12  50  7  912  0  0  13  30  17  1561  0  0  14  40  29  1368  29  1368  15  140  10  665  12  708  16  40  9  1283  23  2860  17  80  27  896  19  632  18  50  10  456  0  0  19  30  25  659  0  0  20  340  11  1074  17  2205  21  140  6  385  0  0  22  20  12  353  13 406 Sustainability 2020, 12, x FOR PEER REVIEW 17 of 23  Table 7 illustrates the estimated number of vulnerable populations that are directly impacted by the 2018 flooding disaster. The results show that 59% of the total population is vulnerable, because they are located in the flooded area. More than 50% of the elderly population (3976 elderly residents) is located in the inundated area, whereas 2148 and 1828 are female and male seniors, respectively, who are expected to be evacuated to a safe location when flooding occurs. In addition, the results indicate that among the flooded senior population, more than 1500 people are elderly, aged over 75 years old. These residents are the most vulnerable because of their physical incapacity, and also due to their vulnerable location to the flooding and the fact that there is no road network to the shelter. This indicates that this senior group would have the most difficulty evacuating in the case of disaster situations.   Table 7 illustrates the estimated number of vulnerable populations that are directly impacted by the 2018 flooding disaster. The results show that 59% of the total population is vulnerable, because they are located in the flooded area. More than 50% of the elderly population (3976 elderly residents) is located in the inundated area, whereas 2148 and 1828 are female and male seniors, respectively, who are expected to be evacuated to a safe location when flooding occurs. In addition, the results indicate that among the flooded senior population, more than 1500 people are elderly, aged over 75 years old. These residents are the most vulnerable because of their physical incapacity, and also due to their vulnerable location to the flooding and the fact that there is no road network to the shelter. This indicates that this senior group would have the most difficulty evacuating in the case of disaster situations.

Discussion
This study proposed a methodology for evaluating the accessibility of evacuation shelters and assessing the vulnerability of existing evacuation shelters and populations regarding disaster risk reduction. Because vulnerability during a disaster is increased by the spatial heterogeneity of an imbalance between population demand and shelter capacity resources, our approach integrated both spatial distribution of shelter demand and resource into the spatial accessibility estimation and shelter deficiency assessment. The above analysis and results indicate that the proposed method, in this study, can be used to map the spatial distribution, and also provide a statistical estimation of shelter demand and the vulnerability of shelters and populations. In most of the previous researches, a specific population ratio was applied to calculate and evaluate shelter demand [20][21][22]. However, this is not realistic, given the diversity of disaster situations. In actual catastrophe circumstances, several factors can influence the real demand and accessibility of shelters, including where the population and shelters are located. Therefore, the proposed algorithm integrated the aspects associated with the location, particularly the spatial distribution of both shelters and population to evaluate the imbalance between shelter demand and existing resources of shelters in a real disaster situation. Applying this algorithm in the case study area in Japan, the method revealed that the shelter capacity was reduced to 43.86% by the flood. In addition, GIS-based mapping was also an essential tool for the assessment of the spatial distribution of evacuation shelters and was used to further support shelter assignments and evacuation planning. The specific locations of the closest shelters to each population area were identified and visualized, using the appropriate distance, as suggested by the guidelines for the case study area, since it has been reported that the residents in Mabi did not know where to go during the disaster [55]. Therefore, our results indicate the nearest designated evacuation and travel route on the GIS map, which can support and enhance the awareness of location and accessibility of shelters for disaster preparation.
According to previous studies based on statistical planning [18,19], in the case of an evacuation, there should be sufficient public shelters available for 25% of both the total population and vulnerable residents. However, our results, in the case study area, revealed that the total shelter capacity planning was available to 13.76% of the total population and 43.45% of the senior group in the town of Mabi. Without considering spatial distribution, it is possible that planning assessment does not reflect the real situation of where residents and evacuation shelters are located when a disaster occurs, including how people can access the shelters. After applying the method, in the case study area, our findings show the spatial heterogeneity of the diverse distribution of population and evacuation shelters. The numerous areas with high population density are located in the northeastern, middle, and southeastern parts of the town, whereas, the areas with high capacity of evacuation shelters are mostly located only in the northeastern and middle part of the case study area. The visualized results confirm that spatial distribution is the critical issue that impacts the availability of shelter and accessibility of residents.
Moreover, the condition and distance of the evacuation routes were also considered in conjunction with this approach. The distance-based function of the actual road network was taken into account in the study, rather than using the Euclidean distance, to estimate the accessibility of facilities, because, in the real world, travel routes could be covered, collapsed, or disconnected due to a disaster, which could decrease the efficiency and impact the evacuation and accessibility of shelters [25]. By comparing the results with the 2018 flooding situation, we also observed that due to the disconnected road network, the shelter demand around the flooding area, based on our study, was significantly increased and enormously exceeded the designated shelter capacity. However, it is not an easy task to build a new shelter to provide greater accommodation capacity in a city where there is limited space. Moreover, it is challenging for the local authorities to determine the designated capacity of evacuation shelters for all the evacuees in a large-scale disaster. On the basis of the statistical results for the case study area, our findings indicated the shelter demand and the location of the remaining sites due to the flood, especially the locations with a demand that exceeds its capacity. Our findings suggest and offer information that could be used by city authorities and disaster administration for planning and managing essential supplies to meet the demand, such as food and sanitation facilities, for example, toiletries and bathroom supplies. In particular, emergency healthcare or medical support should be provided for those considered to be vulnerable groups in this research, who were directly impacted by the disaster situation. The estimated 3976 seniors who experienced the disaster and survived may have lasting mental health impacts [62,63], which, in this case, could require both urgent care and extended plan support. Therefore, this study provides vital evidence for policymakers and practitioners for long term disaster management.
Furthermore, the method also showed that people located around the area of the border of the case study area had difficulties evacuating due to the mountainous condition of their geographic location and the long evacuation distance to the shelters. The results for these specific areas provide beneficial evidence that the local government needs to manage and prepare for evacuations of these areas, specifically the elderly population in these areas, who tend to hesitate to evacuate if the evacuation distances are long [29,30]. On the basis of our findings, it is recommended that some other options for shelters be arranged in those areas, or cooperation with the neighboring towns should be coordinated. Providing more convenient facilities with easier access, particularly facilities within walking distance, would reduce the risk and stress of evacuation [25,26]. In addition, informal shelters, such as parks, temples, local residences, and community halls, could be used to provide accommodation for evacuees during a disaster.
In the future, this study could be further improved. In this study, only people aged 65 and older, as vulnerable groups, were taken into account, since Japan, as a "super-aging" society, has an extremely high proportion of elderly citizens, both in the rural and urban areas. Therefore, additional vulnerable populations, such as children, adolescents, people living with disabilities, and disabled persons, could be included in future research to make the findings more generalizable. Moreover, this research mainly calculated the evacuation routes based on the evacuation distance, which were suggested as guidelines for evacuation, because the objective of the study was to investigate the primary accessibility of evacuation shelters. Differences in terms of altitude were not taken into the account in this study, and this also constitutes a limitation of this research and an area for future improvement. The evacuation time, considering more real-world aspects, should be calculated in the future. Additionally, the hazardous area data applied in this work were from the previous 2018 flooding disaster and the population data were from 2015. Hence, the population distribution could have been altered by the time of this analysis. Additionally, the assessment of the distribution and shelter accessibility could also be investigated in other disaster scenarios, such as earthquakes or tsunamis. Furthermore, this methodology could be implemented in a wide variety of studies due to the advantage of it being a concise technique and the requirement for a minimum number of data, which are usually available for countries like Japan and other countries. A future study could implement the method to assess the impact of further disaster conditions and identify the vulnerability of another area or at a broader scale.

Conclusions
This paper explored a novel approach to evaluate accessibility and shelter deficiency, by integrating both spatial distribution of a population and designated evacuation shelters in a GIS environment. Because the population was unevenly distributed and located in different areas and shelter capacity planning was limited and also unequally allocated, to address these issues and evaluate the vulnerability of shelter and population, evacuation shelter demand based on accessibility was proposed and calculated for the case study area.
By using small spatial units to obtain the spatial distribution of the population and shelters in Mabi, the approach helped to estimate and identify impacted residents, particularly the vulnerable groups when a disaster occurred. In addition, the estimation of shelter demand evidently enhanced the estimation accuracy as compared with using common census data, and based on the accessibility analysis led to identifying the vulnerable shelters regarding insufficient capacity. Additionally, this approach can be conducted using GIS platform technology which provides a practical tool and offers flexibility in terms of data availability for investigating an area to support decision makers and serve as a useful reference for developing and improving emergency management strategies that could reduce disaster vulnerability.
In summary, we believe the assessment in this research could provide a reference for enhancing disaster preparedness, and therefore achieve an effective disaster response in the future. Designated evacuation resources could be reasonably allocated and organized efficiently on the basis of the results from this study. The provision of information regarding vulnerable areas and vulnerable shelters could support local authorities and policymakers in making decisions regarding sustainable reduction of vulnerability.