A GIS-Based System for Spatial-Temporal Availability Evaluation of the Open Spaces Used as Emergency Shelters: The Case of Victoria, British Columbia, Canada

: Canadian emergency management planners have historically ignored the self-motivated evacuation procedures of people who cannot initially choose the safest evacuation areas. In densely developed urban areas, open spaces can be seen as ideal evacuation areas and should thus be included in shelter planning. In this study, the public open spaces in Great Victoria were selected as the study area and evaluated using GIS technologies. A multi-criteria TOPSIS evaluation model was used to conduct comprehensive quantitative evaluations of the open spaces’ safety, accessibility, and availability. Through hybrid process, service area, and POI aggregation coupling analyses, a model is created that provides an overall evaluation at the district level. In addition to providing a model for evaluating open spaces as emergency shelters, applicable to most Canadian cities, this study emphasizes the importance and disadvantages of open space emergency shelters in Canada, which have heretofore been ignored by decision makers. In Great Victoria, we found that the distribution of open spaces does not match the dynamics of the population distribution, meaning that through inadequate preparation some districts lack a safe evacuation place—this in an area where people are at high risk of earthquake disasters and their subsequent effects.


Introduction
Severe earthquakes normally cause extensive but unpredictable loss of life and property, amid which evacuating people to shelter can be difficult. Damaged and collapsed buildings and interruptions in electricity, water, and natural gas supply can cause mass casualties after a severe quake. When an earthquake happens, the first evacuation step is that to stay away from buildings and move to open spaces. Since the local government may have no chance to give out evacuation directives, people need self-rescue, which is a self-motivated evacuation procedure. A proper well-planned self-motivated evacuation can reduce the casualties to a great extent [1]. Some earthquake videos reflect the problem of people not being fully aware of where to escape to, as no clear guiding to such shelters are available; in fact, around 28.9% of residents prioritize roads as evacuation sites [1][2][3]. These findings signify two potential issues: an uneven distribution or lack of safe open spaces in urban areas and unavailability of emergency shelter guidelines.
Earthquake emergency shelters provide people safe places away from active fault zones and other secondary disasters [4]. The earliest emergency shelter studies began in the 1920s in Japan and were developed in the 1970s by Hoshino [2]. However, most cities ISPRS Int. J. Geo-Inf. 2021, 10, 63 2 of 18 that have not suffered a severe earthquake in the last half century tend to minimize the importance of emergency shelter planning. Even so, many cities around the world have a high likelihood of suffering a severe earthquake. Especially in certain vast and sparsely populated countries, such as Canada, densely populated towns may lack emergency shelters. Urban planners should take the full advantage of open spaces as emergency shelters that can support post-earthquake evacuation and reduce casualties in a limited urban area [5].
Through the use of big data acquired from open resources, a quantitative emergency shelter evaluation model can be built to support future urban planning. The data-supported evaluation results adopt a scientific approach that takes different aspects into account, including safety, accessibility, availability, etc. Taking advantage of the analysis of the population distribution and emergency shelters in different periods, such as a weekday, weekends, and night time, the model is completed. Based on the results, planners can adjust and improve post-earthquake emergency evacuation planning to maximize the effective use of open spaces.

Emergency Shelter
Disaster relief (DR) shelters are different from other housing or buildings, with different requirements depending on the duration of the disaster [6]. The post-earthquake recovery process is divided into four stages: immediate relief (hour), immediate sheltering (day), temporary housing (week), and permanent housing (month) [7]. Service time decides the types of shelters and their supporting facilities, and some countries provide standards for building such shelters [8,9]. The emergency shelter service period is from hours to days and does not require permanent buildings. However, a large quantity, good distribution, and high accessibility are needed to support self-motivated evacuation [10,11]. In order to prevent using some private areas that are not open to public, public open spaces have been considered an ideal emergency open space amid natural disasters and especially the public parks are widely used to support evacuation after seismic events [12,13]. An eligible shelter must meet the standards of safety that would guard against disasters and other types of earthquake effects, such as high accessibility, so that residents and rescue workers can access the location, and the capability to at least accommodate 200 people at 2 square meters/person [14,15]. The government of the City of Victoria, Canada publishes reception centers on its website but notes that some places in the city are not covered [16]. They are some indoor long-term shelters, which are high maintenance and insufficient in quantity. Although decision makers and planners encourage people to prepare themselves for earthquakes, people still have no idea where the safest place is. The existing indoor emergency shelters are far from enough. Open spaces have historically been overlooked and not been included in post-disaster planning either. Open spaces as emergency shelters have not been evaluated, as they may also be at risk, and thus would not be able to support the post-earthquake emergency evacuation.

Conditions of Canadian Shelters
The high building code and sparse population has led Canadians to ignore emergency shelters in urban planning. Canada is vast and sparsely populated, with a national population density of only 3.41 per square meter, and its standard for seismic design is 7 [17,18], indicating that new buildings will not collapse under a 7-magnitude earthquake. However, research evaluating building damage in some parts of Victoria and Vancouver has shown that older neighborhoods may sustain up to 30% structural damage [19], indicating that open spaces as emergency shelters are needed because the indoor spaces are no longer safe. The design of Canadian shelters is mainly based on existing buildings, but an upgrade plan could provide long-term accommodation for political asylum seekers, the homeless young, and first nations (Canadian aborigines) [20,21]. However, these shelters have certain disadvantages that diminish the support they can offer during post-earthquake self-evacuation.
According to Schina and Onur et al. [19,22], Canada has little experience in establishing standards for emergency shelters. The Toronto and Vancouver Housing Authority has experience in settling the homeless, providing housing stability by providing daily services to deliver homeless people better living conditions and promote social equity [20]. However, these shelters are not suitable for earthquake response, because they do not consider the widespread effects of natural disasters.

Location Evaluation
Most emergency shelter evaluations give only the score of a location, ignoring the dynamic nature of the population distribution at different times [15]. A hybrid approach is needed to combine multi-criteria decision making methodology and coupling analysis with a view to producing a complete emergency shelter evaluation model.
The multi-criteria evaluation method is widely used to evaluate facility locations. Criteria selection may differ but generally include accessibility (road safety, distance to important facilities), safety (the earthquake disasters and effects), and effectiveness (area, capacity, and service area) [6,[23][24][25]. After the criteria are determined, they are applied through multi-criteria decision-making (MCDM) models [23]. Researchers have developed several methods, such as the analytical hierarchy process (AHP) and techniques for order preference by similarity to ideal solution (TOPSIS), to evaluate existing shelters. AHP has two main advantages: it is a clear systematic quantitative analytical method that is highly suitable for multi-criteria calculation, and it does not require large amounts of data [26]. However, criteria weight was strongly impacted by the subjective matrix, which is not the best method in this study. TOPSIS is popular in large-scale decision-making and location-based analysis since it is based on the distance to the alternative value [6,24].
In summary, this study is intended to provide safe places in Victoria, British Columbia, Canada, where people can wait on government-organized indoor relocation and sheltering arrangements, by building a model for evaluating emergency shelters and district-level shelter services based on Canadian national conditions. Based on the literature review, open spaces are resources for efficiency evacuation. GIS technology, new data, and a multiple criteria evaluation model are used to evaluate the availability of open spaces as voluntary emergency shelters. The distribution matching degree of populations and shelters is also being considered by the coupling analysis. The spatial aggregation of populations is used to identify possible population hotspots, which is represented by point of interest (POI) data. Multi-criteria evaluation and coupling analysis can provide a complete emergency evaluation model [27,28].

Study Area
Vancouver Island in Canada lies on the boundary of the Cascadia subduction zone, which is the most seismically active region in the world. The Juan de Fuca Plate and North American tectonic plate converge here, and the two plates are blocked together and build up the frictional stress, which is the primary cause of the earthquakes. The hundreds of years of enormous stress building up results in a "mega-thrust" earthquake when the stress finally reaches a threshold [27]. Some great historical earthquakes are shown in Figure 1a. According to the research of local geologists, in the next 50 years, 32% of the entire Vancouver Island on the coastal continental plate will suffer destructive earthquakes [16,28].
The six most densely populated regions were selected as the study area, which is known as the Great Victoria region (Figure 1b). The population density is 1275 people per square kilometers-over 300 times the Canadian average population density (4/km 2 ) [2]. As the most livable city in Canada, Great Victoria also attracts hundreds of thousands of tourists or elderly people who come to travel or settle down, and who are vulnerable groups of people after an earthquake happens. tourists or elderly people who come to travel or settle down, and who are vulnerable groups of people after an earthquake happens. Based on previous research by Schina and Onur et al. [19,22], rarely has Canada experienced establishing standards for emergency shelters, so the current shelters would provide limited support during a post-earthquake disaster. Toronto and Vancouver Housing Authority only provide daily service to deliver homeless people for better living conditions [19][20][21]. As Figure 2a shows, there are only seven shelters in Great Victoria, totaling 218 km 2 , and their distribution is only concentrated in Victoria downtown, with at least 197 km 2 having no shelter [29]. The design of shelters mainly relies on existing buildings, as the street images in Figure 2 show; however, their distribution and the amount are not nearly enough to support the post-earthquake evacuation and relocation, suggesting that open spaces need to be introduced in future planning. In summary, taking Great Victoria as a typical example, it would be valuable to improve Canadian emergency planning. Based on previous research by Schina and Onur et al. [19,22], rarely has Canada experienced establishing standards for emergency shelters, so the current shelters would provide limited support during a post-earthquake disaster. Toronto and Vancouver Housing Authority only provide daily service to deliver homeless people for better living conditions [19][20][21]. As Figure 2a shows, there are only seven shelters in Great Victoria, totaling 218 km 2 , and their distribution is only concentrated in Victoria downtown, with at least 197 km 2 having no shelter [29]. The design of shelters mainly relies on existing buildings, as the street images in Figure 2 show; however, their distribution and the amount are not nearly enough to support the post-earthquake evacuation and relocation, suggesting that open spaces need to be introduced in future planning. In summary, taking Great Victoria as a typical example, it would be valuable to improve Canadian emergency planning. tourists or elderly people who come to travel or settle down, and who are vulnerable groups of people after an earthquake happens. Based on previous research by Schina and Onur et al. [19,22], rarely has Canada experienced establishing standards for emergency shelters, so the current shelters would provide limited support during a post-earthquake disaster. Toronto and Vancouver Housing Authority only provide daily service to deliver homeless people for better living conditions [19][20][21]. As Figure 2a shows, there are only seven shelters in Great Victoria, totaling 218 km 2 , and their distribution is only concentrated in Victoria downtown, with at least 197 km 2 having no shelter [29]. The design of shelters mainly relies on existing buildings, as the street images in Figure 2 show; however, their distribution and the amount are not nearly enough to support the post-earthquake evacuation and relocation, suggesting that open spaces need to be introduced in future planning. In summary, taking Great Victoria as a typical example, it would be valuable to improve Canadian emergency planning.

Methodology
Taking advantage of advances in multiple new data, combined with GIS technology, this paper provides an evaluation model to evaluate the availability of open space as an emergency shelter in Canada. The hybrid approach consists of three major steps, and the methodology flow chart is shown in Figure 3.
Step 1: The emergency shelters were evaluated by using the weighted multi-criteria TOPSIS method (technique for order preference by similarity to ideal solution).
Step 2: The population distribution in the different periods was estimated (weekday, weekend, and night time).
Step 3: The shelters' service areas are calculated by the GIS network analysis tool, and the matching degree were measured by coupling analysis (Figure 3).

Methodology
Taking advantage of advances in multiple new data, combined with GIS technology, this paper provides an evaluation model to evaluate the availability of open space as an emergency shelter in Canada. The hybrid approach consists of three major steps, and the methodology flow chart is shown in Figure 3.
Step 1: The emergency shelters were evaluated by using the weighted multi-criteria TOPSIS method (technique for order preference by similarity to ideal solution).
Step 2: The population distribution in the different periods was estimated (weekday, weekend, and night time).
Step 3: The shelters' service areas are calculated by the GIS network analysis tool, and the matching degree were measured by coupling analysis (Figure 3).

Data Collection and Pre-Processing
This study contains several sets of open data. The POI data scraping is based on Python Script by using the Google API to extract the points data on Google Maps. The attributes used in this study are their location and function. Open spaces, roads, and building polygon data are downloaded from Open Street Map (OSM). The type and area of open spaces and building polygons are used in the next analysis. The NPP-VIIRS DNB

Data Collection and Pre-Processing
This study contains several sets of open data. The POI data scraping is based on Python Script by using the Google API to extract the points data on Google Maps. The attributes used in this study are their location and function. Open spaces, roads, and building polygon data are downloaded from Open Street Map (OSM). The type and area of open spaces and building polygons are used in the next analysis. The NPP-VIIRS DNB (Day/Night Band) is the artificial nighttime light data, which has a higher spatial and temporal resolution than traditional nighttime light data and can be used to estimate the active urban area [30]. In this study, the data are used to remove the forest area where there may be some POIs and buildings but there are no permanent residents. Victoria hazard data were obtained from the British Columbia Government Official Website geoscience map section; this geological survey is supported by the National Earthquake Hazard Reduction Program (NEHRP) [31]. Ground amplification motion, slides, and liquefaction hazards are the three types of natural hazard used to screen the open spaces that can be used as emergency shelters. Faulting and ground rupture hazards were not addressed in this program, because the above hazard area has covered the rupture area. The data in full detail is given in Appendix A.
The shape-file format data were obtained, projected, and adjusted into the same projection, the NAD83 UTM Zone 10N Transverse Mercator projection, which guarantees the accuracy of the area and direction. The software used in this research is ArcMap 10.6, provided by the Chinese University of Hong Kong computer lab.

Emergency Shelters Evaluation
Before evaluating, open spaces with a high risk of disaster must be removed first. According to the hazard evaluation data posted by NEHRP, over 20% of the open space areas are under a type of high-degree hazard, whether it is landslides, liquefaction, or ground amplification, and were thus removed. Since Great Victoria is not at risk from coastal communities facing the Pacific Ocean, the waves will only cause a low-risk slow rise in water levels of about 1.5 to 3 m [32]. Therefore, based on the DEM data and Victoria Tsunami Modelling and Mapping Report results, which is a tsunami flooding analysis, coastal areas that are lower than 12 mare defined as high tsunami risk [33]. The selected open space is defined as a potential available emergency shelter ( Figure 4).
(Day/Night Band) is the artificial nighttime light data, which has a higher spatial and temporal resolution than traditional nighttime light data and can be used to estimate the active urban area [30]. In this study, the data are used to remove the forest area where there may be some POIs and buildings but there are no permanent residents. Victoria hazard data were obtained from the British Columbia Government Official Website geoscience map section; this geological survey is supported by the National Earthquake Hazard Reduction Program (NEHRP) [31]. Ground amplification motion, slides, and liquefaction hazards are the three types of natural hazard used to screen the open spaces that can be used as emergency shelters. Faulting and ground rupture hazards were not addressed in this program, because the above hazard area has covered the rupture area. The data in full detail is given in Appendix A.
The shape-file format data were obtained, projected, and adjusted into the same projection, the NAD83 UTM Zone 10N Transverse Mercator projection, which guarantees the accuracy of the area and direction. The software used in this research is ArcMap 10.6, provided by the Chinese University of Hong Kong computer lab.

Emergency Shelters Evaluation
Before evaluating, open spaces with a high risk of disaster must be removed first. According to the hazard evaluation data posted by NEHRP, over 20% of the open space areas are under a type of high-degree hazard, whether it is landslides, liquefaction, or ground amplification, and were thus removed. Since Great Victoria is not at risk from coastal communities facing the Pacific Ocean, the waves will only cause a low-risk slow rise in water levels of about 1.5 to 3 m [32]. Therefore, based on the DEM data and Victoria Tsunami Modelling and Mapping Report results, which is a tsunami flooding analysis, coastal areas that are lower than 12 mare defined as high tsunami risk [33]. The selected open space is defined as a potential available emergency shelter (  The first step of the evaluation is to determine the criteria. In this study, emergency shelters serve for the period within the 24 h after an earthquake happens, and extra facilities or permanent buildings are not required. There are two principles of selecting the indicators that are independent of each other and quantifiable in estimation. The earthquake emergency shelter is evaluated on three aspects, including safety, accessibility, and The first step of the evaluation is to determine the criteria. In this study, emergency shelters serve for the period within the 24 h after an earthquake happens, and extra facilities or permanent buildings are not required. There are two principles of selecting the indicators that are independent of each other and quantifiable in estimation. The earthquake emergency shelter is evaluated on three aspects, including safety, accessibility, and availability, and each aspect includes different criteria. Based on previous research, there are three essential requirements (safety, accessibility, and capability), including six criteria (building collapsing, fire hazard, roads accessibility, rescue response of firefighting, medical service, and capacity) that are relevant to emergency shelter site evaluation after an earthquake [6,[23][24][25]. The following three paragraphs contain the criteria explanations.

1.
Criteria of Safety: People must avoid the high-risk places, such as areas at risk of building collapsing and fire. Buildings under high risks of liquefaction, ground amplification, and landslides have a high possibility of collapsing, because of the unstable geological structure [34,35]. Additionally, in Great Victoria, most individual houses are under three floors, and the wood frame structure buildings have less possibility of blocking the roads than the concrete buildings [36][37][38]. Therefore, those buildings higher than three floors and that have liquefication, high ground amplification, or a landslide risk have a high possibility of blocking the roads. The proportion of collapsing buildings within a 500-m buffer was used to evaluate the building collapsing hazard. Because the scope of gas station fire is unpredictable, the numbers of gas stations within the optimum service radius (500-m buffer) were counted to represent the fire hazard possible.

2.
Criteria of Accessibility: The accessibility condition includes criteria of road accessibility and rescue response speed [23,24]. The roads layer is overlaid on the natural hazard data layer, and the road conditions are divided into high, moderate, and low risks, which are given the weight of 1, 0.5, and 0, respectively. The sum of the road length timing their corresponding weight represents the roads accessibility within the 500-m buffer. TOPSIS approaches provide decision support to shelter location selection, and the calculation steps are orderly listed below, as a list from Process (1)- (7). Formula (1): A decision matrix is applied to build up the model and to evaluate the effectiveness of every open space. Because the calculated criteria values have different units, they need to be transformed to a dimensionless form through a mathematics model. Formula (2): The criteria are not equally important, thus relative weights are assigned to different criteria. Due to Canada's lack of a previous weight standard, for reference in this case the information entropy was used to determine the weight. Formulas (3) and (4): Finally, according to the distance to the worst and best alternative value, this study provides a score of the candidate sites, between 0 and 1; the higher score open space is more valuable as an emergency shelter during the post-earthquake evacuation process. where tj is the entropy value, and pij is the proportion of emergency shelter i of the total value of j.
To reduce the extreme values' impact, the largest and smallest value are removed in the entropy weight calculation where: J+ is the criteria having a positive impact, and J− is the criteria having a negative impact.
(6) Calculate the distance between the target worst or best alternative value: where D w i : the distance from the criteria value to the worst alternative, and D b i : the distance from the criteria value to the best alternative. (7) Open space evaluation score Si: The score is converted to a hundred-mark system, so the full mark is 100.

Service Area and Population Distribution Coupling Analysis
To calculate the service area of the emergency shelters, it must be understood that people's behavior influences the open space service efficiency, such as limited walking speed and choosing the familiar targets. Some researchers applied Simulex technology to simulate people's walking speed after an earthquake; the result shows that the walking speed is between 0.5 m/s and 2 m/s [2]. Other research, using videotape data, showed that the speed is between 2.3 m/s and 3 m/s, which does not consider vulnerable groups [39]. According to previous studies, the general speed interval is between 0.5 m/s and 3 m/s, the elderly people and young children move at slower speeds, and the average evacuation speed is assumed as 1.2 m/s (4.3 km/h). This speed value is lower than the average value in this study; the service area evaluation is thus adopting a conservative approach. During an evacuation, people avoid the high-risk areas and roads so that the high-risk roads are removed before processing the network analysis. A GIS-based network analyst tool was used to determine the theoretical 15-min service area.
(1) The coupling analysis is a quantity analysis that represents how the emergency shelter service area matches the population distribution [30,39]. Given the service area, the buildings in the service area compared with the total building area can represent a coupling degree. The population in the service area can be estimated according to the coupling degree. The opposite of coupling, the high building density functional zoning but without any emergency shelter, is defined as a gaping hole area. The calculation of the gaping hole method is similar to the coupling degree, which is the proportion of the building area in the gaping hole area of the total district building area. Some functional areas are hotspots that attract a large population in some period. Through the spatial join method, the building functions are defined by using POI data: educational buildings (including some art school, interest training organizations, etc.) and working buildings (e.g., office blocks, factories, etc.), business buildings (e.g., shopping mall, individual shop, etc.), recreation buildings (e.g., restaurants, board game bar, museums, antique stores, etc.), and residential buildings (e.g., house, apartment, condo, etc.). Combining the service area and population distribution results, the gaping hole also can be founded. The proportion of the building area in the service area of the total district building area, value c in Formula (5), represents the coupling degree or the proportion of the gaping hole. Next, Formula (6) In the day time, bx are all functional buildings, therefore all buildings except residential buildings are included to calculate the coupling degree. In the night time, light data are used to remove the vacant houses in the suburban area where the value digital number value is 0, and only bx are residential buildings. The buildings bx in the gaping hole estimation calculate the proportion of the gaping hole.
(2) The number of people in the service area or in the gaping hole area: P = c × district total population (6)

The Emergency Shelter Evaluation Score
There are a total 375 public open spaces in Great Victoria, and 60 of them are available as earthquake emergency shelters. Emergency shelters must meet the basic security requirements to be safer than roads, indoor spaces, and densely built areas. People should avoid evacuating to the other 315 open spaces, which may seem safe and open but are high-risk areas ( Figure 5). Each of the open spaces has a score; in order to simply represent the availability of emergency shelters, the evaluation scores were assigned into 5 categories, using the natural breaks classification method. A decision value from a low to a high level means that the open space is less suitable to very suitable as an emergency shelter, respectively (detailed information in Figure A2).
As Figure 5 shows, Saanich, the largest municipality, has the most emergency shelters, with 39 emergency shelters covering 1,111,865 m 2 and 15 of them are classified as low or less suitable conditions. The high fire risk and low rescue responding speed are two main reasons reducing the values of the emergency shelters. Victoria has 15 open spaces (114,398 m 2 ) that can be used as emergency shelters, generally uniformly distributed among communities, such as Fernwood, Oakland, and Downtown. All emergency shelters areas in Victoria are small, which is disadvantageous in a densely populated district. High Rock Park is the only emergency shelter available in Esquimalt (64,514 m 2 ), which is large in size and in the central area. As the highest point in Esquimalt, the landform reduces its accessibility. Fortunately, rescue response is high in this region, which may reduce casualties. The Langford and Colwood emergency shelters' surrounding areas are undeveloped, which offer outstanding accessibility and capability. They are high value emergency shelters, but they are far from the urban area so that they may have a low using efficiency. In Oak Bay, the shelter is in the northern part, but far from some of the populated coastal communities; it is a golf course (64,514 m 2 ), so although its safety and accessibility are high, the surrounding areas' functions mainly serve the purpose of golfing. Each of the open spaces has a score; in order to simply represent the availability of emergency shelters, the evaluation scores were assigned into 5 categories, using the natural breaks classification method. A decision value from a low to a high level means that the open space is less suitable to very suitable as an emergency shelter, respectively (detailed information in Figure A2). As Figure 5 shows, Saanich, the largest municipality, has the most emergency shelters, with 39 emergency shelters covering 1,111,865 m 2 and 15 of them are classified as low or less suitable conditions. The high fire risk and low rescue responding speed are two main reasons reducing the values of the emergency shelters. Victoria has 15 open spaces (114,398 m 2 ) that can be used as emergency shelters, generally uniformly distributed among communities, such as Fernwood, Oakland, and Downtown. All emergency shelters areas in Victoria are small, which is disadvantageous in a densely populated district. High Rock Park is the only emergency shelter available in Esquimalt (64,514 m 2 ), which is large in size and in the central area. As the highest point in Esquimalt, the landform reduces its accessibility. Fortunately, rescue response is high in this region, which may reduce casualties. The Langford and Colwood emergency shelters' surrounding areas are undeveloped, which offer outstanding accessibility and capability. They are high value emergency shelters, but they are far from the urban area so that they may have a low using efficiency. In Oak Bay, the shelter is in the northern part, but far from some of the populated coastal communities; it is a golf course (64,514 m 2 ), so although its safety and accessibility are high, the surrounding areas' functions mainly serve the purpose of golfing.

Service Area and Day Time Population Coupling
Due to the limitations associated with social and environmental conditions, the distribution of emergency shelters in different districts is uniform. Gaping holes are evident without open space available for use as emergency shelters, increasing the levels of evacuation uncertainty and risk.
As Figure 6 and Table 1 shows, Saanich has the largest service area and the highest population coupling degree (8.1%) during the daytime. Some 49.4% of the people are in the gaping hole; one such hole is the University of Victoria, where students are concentrated during weekdays. Another is Royal Oak, which has some shopping centers where people are concentrated on weekends. Emergency shelter service areas are highly overlapping in Victoria. There are five gaping holes in Victoria, covering 76.3% of the people, and the largest gaping hole is in the downtown area where is the center of business, recreation, and works. Therefore, whether it is weekdays or weekends, downtown Victoria always has a high population density. There is likewise a recreation hotspot associated with kayaking clubs. Because evacuation risk is unpredictable when people are on the water, it can be considered as a high-risk point. People living in or traveling the peninsula border of Esquimalt cannot reach a proper emergency shelter within 15 min. The coupling degree of Oak Bay is 0 because the golf course surrounding area has no functional buildings. The whole Oak Bay district has very limited functional buildings, and scattered working offices or recreation centers are in the southern part. There are some industrial parks in Esquimalt and Langford; therefore, on weekdays, the working zone is more likely to be densely populated. Langford downtown has no emergency shelter available, so it has the largest proportion of people in the gaping hole (83.9%). Colwood has the lowest population, and what is more, service areas are in sub-central areas of the districts, further reducing service availability. Due to the limitations associated with social and environmental conditions, the distribution of emergency shelters in different districts is uniform. Gaping holes are evident without open space available for use as emergency shelters, increasing the levels of evacuation uncertainty and risk.
As Figure 6 and Table 1 shows, Saanich has the largest service area and the highest population coupling degree (8.1%) during the daytime. Some 49.4% of the people are in the gaping hole; one such hole is the University of Victoria, where students are concentrated during weekdays. Another is Royal Oak, which has some shopping centers where people are concentrated on weekends. Emergency shelter service areas are highly overlapping in Victoria. There are five gaping holes in Victoria, covering 76.3% of the people, and the largest gaping hole is in the downtown area where is the center of business, recreation, and works. Therefore, whether it is weekdays or weekends, downtown Victoria always has a high population density. There is likewise a recreation hotspot associated with kayaking clubs. Because evacuation risk is unpredictable when people are on the water, it can be considered as a high-risk point. People living in or traveling the peninsula border of Esquimalt cannot reach a proper emergency shelter within 15 min. The coupling degree of Oak Bay is 0 because the golf course surrounding area has no functional buildings. The whole Oak Bay district has very limited functional buildings, and scattered working offices or recreation centers are in the southern part. There are some industrial parks in Esquimalt and Langford; therefore, on weekdays, the working zone is more likely to be densely populated. Langford downtown has no emergency shelter available, so it has the largest proportion of people in the gaping hole (83.9%). Colwood has the lowest population, and what is more, service areas are in sub-central areas of the districts, further reducing service availability. Images (a-f) they are, respectively, Esquimalt, Langford, Victoria, Saanich, Colwood, and Oak Bay, depicting the population density and service area coupling maps in the day time; the yellow and red polygons are, respectively, service areas and gaping holes. Images (a-f) they are, respectively, Esquimalt, Langford, Victoria, Saanich, Colwood, and Oak Bay, depicting the population density and service area coupling maps in the day time; the yellow and red polygons are, respectively, service areas and gaping holes.  Figure 7 and Table 2 show the nighttime population distribution and emergency shelter coupling degree. In Saanich, people in the gaping holes and service areas are similar, both around 30%. Around 40% of the people live in the semi-rural area, where the buildings are mainly individual houses and scattered around the farmland or national parks. Southeast Saanich has the largest gaping hole, which are mostly apartment buildings; therefore, it is the least vulnerable to any type of earthquake disaster and effect. Nearly 45% of the people in Victoria can get to emergency shelters on time. The riskiest region is from downtown to northwest Victoria, and the peninsula of Victoria, because of the high apartment building density. Esquimalt has the highest coupling degree, and there is no gaping hole because the emergency shelter is in the central downtown area, and the coastal residential are mostly individual houses. Oak Bay golf course covers 12.4% residential population at night; however, around 80% live in the southern part and are in the gaping hole. As an aging community with a long coastal line, Oak Bay has the highest pressure of emergency evacuation after an earthquake. Langford and Colwood have a high score of emergency shelters but the lowest coupling degree. Especially in Langford downtown, over 90% of people are in the gaping hole.   Figure 7 and Table 2 show the nighttime population distribution and emergency shelter coupling degree. In Saanich, people in the gaping holes and service areas are similar, both around 30%. Around 40% of the people live in the semi-rural area, where the buildings are mainly individual houses and scattered around the farmland or national parks. Southeast Saanich has the largest gaping hole, which are mostly apartment buildings; therefore, it is the least vulnerable to any type of earthquake disaster and effect. Nearly 45% of the people in Victoria can get to emergency shelters on time. The riskiest region is from downtown to northwest Victoria, and the peninsula of Victoria, because of the high apartment building density. Esquimalt has the highest coupling degree, and there is no gaping hole because the emergency shelter is in the central downtown area, and the coastal residential are mostly individual houses. Oak Bay golf course covers 12.4% residential population at night; however, around 80% live in the southern part and are in the gaping hole. As an aging community with a long coastal line, Oak Bay has the highest pressure of emergency evacuation after an earthquake. Langford and Colwood have a high score of emergency shelters but the lowest coupling degree. Especially in Langford downtown, over 90% of people are in the gaping hole.   6. Discussion

Spatial and Temporal Distribution of Emergency Shelters
In general, Langford, Colwood, and Oak Bay lack open spaces as emergency shelters although their scores are high. Esquimalt has better evacuation conditions than the other districts because of the shelter in the central urban area. The only disadvantage is the low accessibility, which leads to some people in the coastal area needing more time to get to the shelter. The emergency shelters in Victoria downtown are mostly of low suitability. The east part of the Victoria emergency shelter areas has the characteristics of being small, with a high score, and equally distributed. In Saanich, emergency shelters are along the Patricia Bay Highway. The fire hazard and low rescue responding speed are the main reasons that some shelters have low scores.
The population distribution and service area coupling degree in the night time are higher than in the day time, and the weekend coupling degree is better than the workdays degree. Because the working areas are normally with a high building density and lacking greening space or playgrounds, they become areas with the highest risk. Business and recreational areas attract more people on weekends. From Victoria downtown to Saanich uptown, the central belt of business and recreational facilities should be highlighted as a gaping hole. To be more specific towards a district-level evaluation, the north part of Saanich has only a few emergency shelters; it is a thinly populated, semi-rural area mainly occupied by agricultural land where people can safely evacuate [34]. The gaping holes in Saanich are in the southeast part-the greater Victoria sub-center-where no open space is available as an emergency shelter. There are three high-risk gaping holes in Victoria: Peninsula of Victoria, Victoria downtown, and northwest Victoria, due to the high population density in both the day time and at night. Esquimalt, being a peninsula, and Oak Bay, with its long coastal line, are both relatively isolated areas whose seaward facing may increase evacuation pressure. In Langford and Colwood, the developed urban area has no available open spaces suitable for use as emergency shelters, allowing only limited support of post-earthquake evacuation. The industrial parks of Langford correspond to gaping holes. Accordingly, if an earthquake happens on a weekday, people in these districts are at high evacuation risk.

Suggestions for Future Planning
Due to the limitations of historical city planning, a lack of urban disaster mitigation plans, and the social gap among different districts, urban emergency shelters are not included in the disaster prevention plans of greater Victoria. The open spaces used for early self-evacuation in greater Victoria have been evaluated based on the criteria scores ( Figure A2), which provide an overall picture for decision makers. Based on the results of this study, some insights, limitation, and suggestions can be provided to guide decision makers in making plans. Four points summarize the general evacuation conditions and suggestions: 1.
This study shows that in Great Victoria, public open spaces are an available resource for emergency sheltering during post-earthquake self-motivated evacuations. Emergency shelters can make up for shortages in existing shelters, releasing evacuation pressure, and reducing casualties. Furthermore, across the whole western coastal area of Canada it is also possible to apply open spaces as earthquake emergency shelters.

2.
Although the evaluation scores of the emergency shelters differ, they are all valuable for supporting evacuation. Highly scoring emergency shelters should be retained and developed carefully to take full advantage of them. By fully using these emergency shelters, the cost of evacuation planning can be reduced, consistent with the principle of integrating daily facilities with disaster prevention facilities. 3.
Residents can use the evacuation conditions map to find the nearest available safe place that can serve as an emergency shelter; they should stay in emergency shelters and wait for future rescue and relocation arrangements. Decision makers should also improve the whole first-time evacuation plan, focusing on high-risk gaping holes and densely populated areas.

4.
Evacuation plans should reflect population distribution dynamics, distinguishing and considering differences between day and night, weekday, and weekends. Therefore, the population concentrated in different areas, such as offices, industries, schools, etc.
(weekdays), recreation places, business, etc. (weekends), and residential (night), are different during different time periods. A population simulation can be introduced to estimate the population dynamics distribution [40]. Planners should pay more attention to this issue in the future.

Limitations and the Next Improvement
There are not only limitations in this contribution as mentioned in the methods, but also some limitations of the data resource and evaluation methodology. Due to the limited nature of the open resource data, some types of the earthquake effects and secondary disasters, such as debris and flooding, are underestimated. They are not high risks in this study but may highly influence other cities. Since the building data cannot support the debris simulation, the capability estimation is only acceptable under the quick approach and wide-scale conditions. In the aspect of methodology, since Canada lack research on emergency shelters, there is no reference in evaluating the open space as emergency shelters either. The entropy weighted method is based on the data itself, and it loses sight of how the social rules influence the evacuation behaviors. In sociology-related evaluation research, the professional grading by Canadian sociologists and policy makers is also necessary to include. Although the criteria selected in this study are applicable to most areas in Canada, still some criteria are ignored; for example, the gradient of the emergency shelter. In the hilly cities, some open spaces are built on rough ground, which can reduce their capability. Finally, because Canada lacks emergency shelter planning, the real usage efficiency of the emergency shelters cannot be estimated. We suggest that the signage needs to be built up near by the high-value emergency shelters. Sociologists can use this study results to further simulate human behaviors in the post-earthquake evacuation process. The evacuation simulation can be applied to further understand the possible shelter disadvantages of the surrounding environments. It also would be interesting to explore the Google thermal population thermal graphic, possibly applying it to the emergency evacuation plan making.

Conclusions
This study provides a first-time self-motivated post-earthquake emergency shelter evaluation model for Canada, which can be applied to most North American cities on the seismic belt. Public open spaces are used as emergency shelters, providing a clear safe evacuation destination for the public while taking full advantage of the open spaces. The advantages and disadvantages of each district are outlined so that the planners and decision makers can devise customized plans. Second, although Canadian geologists have conducted some research into earthquake prediction and hazard evaluation, meaningful scientific results have not been implemented for practical urban planning. Previous urban planning of greater Victoria has been limited by a lack of understanding of the open spaces' potential for use as emergency shelters, which may lead people to evacuate to high-risk highly appreciated.