Earthquake Evacuation Choice and Management in a Developing Archipelagic Country—A Case Study of Surigao City, Philippines

: The Philippines is a developing archipelagic country in Southeast Asia. The country is susceptible to multiple natural disasters, speciﬁcally earthquakes. This implies the signiﬁcance of understanding earthquake evacuation choice in order to design effective planning and management of evacuation to minimize chaos, damage cost, and the loss of lives. This study investigated the determinants of earthquake evacuation and proposed earthquake evacuation planning and management in the Philippines, featuring the case study of Surigao City. The study used the primary dataset of 1055 observations gathered in 2019 and applied the nested logit model (NLM) to investigate the potential factors of earthquake evacuation decisions. We considered three output variables: evacuation choice, evacuation duration, and travel mode choice. We found that residents were more likely to evacuate their homes upon receiving an earthquake warning and move to a public shelter or open space. Additionally, respondents were more inclined to leave their homes when their houses suffered from moderate to severe/complete damage or when electricity and water supply were cut-off. Respondents were most likely to walk to evacuation centers as the majority of residents initially moved to the nearest open space immediately after an earthquake and stayed in an open space for less than 6 hours. No correlation was found between personal and household income factors with evacuation choice and travel mode choice. Furthermore, the study used the Analytical Hierarchy Process (AHP) to determine areas suitable for earthquake evacuation using insights from local government ofﬁcials and planners. The areas identiﬁed for earthquake evacuation were developed to support evacuation planning and management.


Introduction
The Philippines is a developing archipelagic country with 7,641 islands, located in Southeast Asia [1]. The country experienced an average economic growth of 6.6% annually from 2017 to 2019 [2] with a total population of 104.9 million in 2018 [1]. According to the World Risk Report 2018, the Philippines is the world's third-highest disaster-prone country, with an index value of 25.14% [1]. At least 60% of the country's total land area is prone to multiple disasters, and 74% of the total population is exposed to these disaster impacts [3]. According to official data from the Philippine Institute of Volcanology and Seismology (PHIVOLCS), the Philippines is geographically located near 54 active volcanoes, and earthquakes frequently occur up to eighteen times per day, on average [4]. Its proximity to the Eurasian and Pacific tectonic plates makes it more vulnerable to seismic activities that can generate high-magnitude earthquakes in the most seismic-prone areas (e.g., Manila and Davao) [5]. An earthquake with a 5.0-6.0 M is likely to occur over most of the Philippines in the next 50 years with a near-100% probability, and an earthquake with a 7.0 M will almost certainly occur anywhere across the entire country, also with a near-100% probability [5]. Historically, three destructive earthquakes happened, on average, per year in the Philippines [4]. These earthquakes caused various degrees of damage to various structures and utilities. The epicenter of an earthquake mainly triggers tsunamis and coastal flooding, thereby causing damage to properties and infrastructure and causing human and animal deaths.
Some of these earthquakes measured greater than a magnitude (M) 6, as listed in Table 1. The earthquake in Surigao City on 10 February 2017 at 10:03 PM (Pacific Standard Time) led to the third-highest damage cost of PhP 720 million ( 51.82 PhP = 1 USD) among other earthquakes in the Philippines since 1990. This was a 6.7 M earthquake, which led to ground shaking in northeastern Mindanao, specifically in Surigao del Norte. The epicenter of the earthquake was located 16 km offshore northwest of Surigao City, and the Surigao Strait (a strait between the Bohol Sea and the Leyte Gulf of the Philippine Sea) has a depth of around 10 km in the southern Philippines [6]. Several municipalities were hit by the earthquake, though Surigao City suffered from the strongest ground shaking, with a PHIVOLCS Earthquake Intensity Scale (PEIS) of VII (Destructive). Some buildings, roads, and bridges incurred considerable damage. Official data of the Department of Public Works and Highways (DPWH) showed that four roads and five bridges suffered different levels of damage with an estimated cost of more than PhP 100M for rehabilitation and repair [7]. Figure 1 shows the locations of the damaged roads and bridges caused by the 2017 Surigao earthquake, and the types of evacuation facilities are highlighted with various colors. The number of each facility type is provided in Table 2. Schools accounted for the largest share of evacuation destination choice, which may be due to their capacity and accessibility (shorter distance). Gymnasiums/covered courts accounted for the secondlargest proportion of evacuation sites due to their capacity and accessibility to barangay halls and other essential facilities.
Negative impacts of earthquakes can be mitigated through effective earthquake disaster management planning based on lessons learned from global experience for local best practices. The Republic Act or "The Philippine Disaster Risk Reduction and Management Act of 2010" was issued to develop policies and planning for risk assessment, preearthquake warning, capacity building, awareness-raising, risk reduction, and preparedness for effective response and prompt recovery [8]. Evacuation planning for earthquake disasters is an important component of disaster preparedness in the disaster management cycle. This component deals with necessary preparations, such as what to do, where to go, or whom to call during an emergency. Facilities and resources are also considered in the three stages of the evacuation process: (1) before evacuation (transmission of evacuation warnings and instructions), (2) during evacuation (evacuation guidance and routes), and (3) after the evacuation process (evacuation centers) [9].
A thorough understanding of the determinants of natural disaster evacuation (e.g., destination, travel mode, and evacuation duration) is needed to design disaster management and planning programs. Natural disaster evacuation (e.g., earthquakes, tsunamis, flooding, and hurricanes) is significantly affected by various factors. A received warning, a distance closer to the threat, and structural damage are statistically associated with earthquake evacuation choice [10][11][12]. A study of earthquake evacuation using the Immersive Virtual Reality (IVR) tool in Auckland showed that building occupants were influenced by other people during earthquake evacuation, specifically the authorities [13]. Only 13% of evacuees did not pay attention to the instructions of the authorities and other people around them and tended to move out of the building [13]. A similar study in rural areas showed that built-up environments and disaster risk perception were positively associated with evacuation choice behavior [14]. A behavioral study of tsunami evacuation facility choice revealed that evacuees tended to move to designated facilities close to their locations with high altitudes and floors and away from the sea [15]. A recent study of tsunami evacuation in Indonesia indicated that a strong relationship among people in local communities is a key factor to decrease the number of fatalities; 83% of respondents decided to evacuate because they saw other people evacuating [16]. A study of tsunami evacuation behavior and travel mode choice in the Cascadia subduction zone showed that two-thirds of respondents evacuated, and older citizens and females were more likely to travel by foot during evacuation [17]. A study on flooding evacuation in Quezon City, Philippines, showed that evacuation choice was mainly influenced by household characteristics (i.e., gender, educational level, presence of children, years in residence, house ownership, number of dwelling stories, and type of house material) and hazard-related factors (i.e., distance from flood source, level of flood damage, and source of warning) [18]. Another study of hurricane evacuation decisions revealed that the destination choice was most significantly caused by hazard severity, income, the type of emergency, age, ethnicity, education level, and pet ownership [19][20][21]. Distance to the evacuation destination, the population at the destination, and destination location type were the potential variables of hurricane evacuation decisions in South Carolina [22]. Another study of hurricane evacuation choice in northern New Jersey, using stated preference (SP) household survey data, revealed that not having a motorized vehicle was the strongest predictor of evacuation choice [23]. Also, low-incomes, low education levels, and high population densities were positively associated with travel modes for evacuation, other than private vehicles [23]. These findings can help policymakers and planners prepare planning and management for natural disaster evacuation. The evacuation planning should be specific to different areas and to each natural disaster, e.g., earthquakes, tsunamis, flooding, and hurricanes.
Given the need for effective earthquake evacuation planning, this study investigated the potential determinants of evacuation choice using revealed preference (RP) data and proposed specific earthquake evacuation planning in Surigao City, Philippines. This city was featured as the case study as it had one of the strongest earthquakes in recent years. Residents, therefore, would have a more intact recollection of their experiences during evacuation during the night. The Surigao earthquake also caused the third-highest damage cost in the Philippines since 1990 (see Table 1). Sociodemographic characteristics, earthquake-related damage variables, and housing utilities and assets were taken into account as explanatory variables. Evacuation choice, evacuation duration, and travel mode choice for evacuation were considered as the output variables. This study used the primary RP data of 1055 observations collected in 2019 through a face-to-face post-event survey. The nested logit model was applied for statistical data analysis. The novelty of this study is in that it took an exhaustive set of factors of earthquake evacuation choice in a developing archipelagic country into account. The findings of this study provide significant contributions to earthquake evacuation planning both in the local and national levels, thereby contributing to the mitigation of damage cost and loss of lives and the increase of community resilience and sustainability.
The rest of this paper is structured as follows. Section 2 is the methodology of the study, including the study area, research design and hypotheses, descriptive statistics of the data sample, and the conceptual frameworks for statistical data analysis. Section 3 consists of the model estimation results and discussion. Section 4 proposes earthquake evacuation planning based on the Analytical Hierarchy Process (AHP) for Surigao City. The last section summarizes the results and suggests future research directions to advance human knowledge in earthquake evacuation planning and management.

Study Area and Data Source
The Philippines is an archipelagic country, with a land area of roughly 300,000 km 2 [1]. The country consists of three regions: Luzon (north), Visasyas (center), and Mindanao (south). Surigao City is located in the Surigao del Norte province, the northeastern tip of Mindanao. The city consists of 54 barangays, with a total area of 245.34 km 2 [27]. The population of the city was 154,137 in 2015 [28], with an average annual growth rate of 1.83% from 2010 to 2015 [27]. The study area covered the mainland part of the city, which hosts the majority of the city's population. The mainland consists of 33 barangays (Barangay is the smallest administrative division in the Philippines.) that are clustered into rural-coastal (i.e., Cabongbongan, Capalayan, Day-Asan, Nabago, Orok, San Isidro), rural-inland (i.e., Anomar, Balibayon, Bonifacio, Danao, Mabini, Mapawa, Mat-I, Poctoy, Quezon, San Roque, Serna, Silop, Sukailang, Trinidad), sub-urban (i.e., Cagniog, Ipil, Lipata, Mabua, Punta Bilar, Rizal, Sabang, Togbongon), and urban areas (i.e., Canlanipa, Luna, San Juan, Taft, Washington) [29]. A map of the study area is illustrated in Figure 2. The study was also limited to the mainland area as it is situated near the active Surigao Fault, thereby making this area prone to earthquake events, as illustrated in Figure 2. This study covered an area of 139.746 km 2 or 57% of the city's total land area [29].
This study hypothesized that earthquake-related damage variables, dwelling utilities and assets, and socioeconomic characteristics affected evacuation choice, travel mode choice, and evacuation duration. For the earthquake-related damage variables, severe house damage and water and electricity cutoffs might force residents to evacuate from their homes to an open space or a public shelter to get various supplies and help. For the dwelling utilities and assets, people living in concrete-made houses might not evacuate and instead decide to stay home. For the socioeconomic characteristics, people with small children, elderly people, and PWDs might decide to stay home or evacuate to an open space near their homes.
The sampling technique used a combination of stratified and convenience sampling. After determining the sample size, the samples were clustered based on the clustering of barangays within the mainland of Surigao City: urban, sub-urban, rural-inland, and ruralcoastal. Convenience sampling was used during the actual data collection as explained in the succeeding section. These sampling techniques were found to be economical, fast, and simple to conduct. The sample size was computed using Equation (1) [30]: where n is the sample size; z is the z-score corresponding to the desired level of confidence (probability of error); p is the estimated proportion of the population which might be based on prior research, pilot study, estimates from an experienced researcher(s); q = 1 − p; and e is the tolerable margin of error or precision of the estimate. Sustainability 2021, 13, x FOR PEER REVIEW 6 of 26 For a level of confidence of 95%, the corresponding z score is 1.96. The margin of error was selected to be e = 0.03 in our study. Given no prior studies, the most conservative estimate of 0.50 for p was used [30]. Therefore, the sample size was at least n = 1067 based on Equation (1). The survey took place for 10 days from 1 to 11 February 2019. The actual sample size collected was 1079. The proportion of respondents by barangay was determined on the basis of the population distribution. Of the sample data, 74 observations were collected in rural-coastal areas, 226 in rural-inland areas, 210 in suburban areas, and 569 in urban areas. The distribution of the sample data is listed in Table 3. Some questionnaire forms with incomplete information were deliberately removed to avoid data inconsistency. After data polishing, there were 1055 observations used for data modeling. For a level of confidence of 95%, the corresponding z score is 1.96. The margin of error was selected to be e = 0.03 in our study. Given no prior studies, the most conservative estimate of 0.50 for p was used [30]. Therefore, the sample size was at least n = 1067 based on Equation (1). The survey took place for 10 days from 1 to 11 February 2019. The actual sample size collected was 1079. The proportion of respondents by barangay was determined on the basis of the population distribution. Of the sample data, 74 observations were collected in rural-coastal areas, 226 in rural-inland areas, 210 in suburban areas, and 569 in urban areas. The distribution of the sample data is listed in Table 3. Some questionnaire forms with incomplete information were deliberately removed to avoid data inconsistency. After data polishing, there were 1055 observations used for data modeling.  Table 4 illustrates the earthquake evacuation choice caused by the 6.7M 2017 Surigao earthquake. Of all the respondents, 62% decided to leave home, while 38% stayed home. Of those in the leave-home choice group, 52% of the respondents evacuated to open spaces, and 10% of the respondents moved to public shelters for evacuation. Of the evacuees in open spaces, 37% of respondents stayed for less than 6 h, while 15% stayed for at least 6 h. Of the evacuees in public shelters, 7% decided to stay there for less than one day, and 4% stayed for at least one day. Table 5 shows the travel mode choice during the evacuation. As mentioned, 62% of the respondents left home during the evacuation. The majority of evacuees traveled by walking (53%), followed by two-/three-wheeler modes (6%) and the other modes (2%). These survey results are consistent with a study of submarine landslide earthquake evacuation in Indonesia [16].  Table 6 presents the descriptive statistics of variables related to damage caused by the 2017 Surigao earthquake and dwelling utilities and assets. The majority of respondents (75%) did not receive any earthquake warning to evacuate. About 16% of residents received a warning through megaphones made by the local authorities or through word-of-mouth from other residents, and 9% received a warning through telecommunication. According to the data sample, 47% of houses did not suffer from the earthquake, 46% of houses suffered moderate damage, and only 7% of houses sustained severe/complete damage. The majority (67%) could not access electricity after the earthquake hit, while 36% of houses incurred the water supply cut-off. For dwelling utilities and assets, about nine-tenths of respondents lived in their own houses. The majority of houses are one-story houses (84%) and are made of concrete (66%). Most respondents have access to safe water (95%), electricity (97%), and toilets (98%). Based on survey data, 97% of residents have no generator set (genset). The majority of households have no motorcycles (59%), cars (95%), multicabs (98%), trucks (100%), or other vehicles (93%). Table 7 presents the sociodemographic characteristics of the respondents. The majority of them were female (62%), married (65%), and 18-29 years old (26%). At least 80% of the respondents had a high school diploma or higher. In terms of occupation category, 37% of the respondents were retired/unemployed/student, 31% were employees, and 32% were self-employed. Just over half (56%) of respondents were not household heads. Most of the respondents belonged to families with four members (24%), followed by five-member families (20%), and three-member families (16%). Based on the data sample, the majority of residents belonged to households with no children aged ≤ 5 years (65%), no elderly citizens (81%), no PWDs (96%), and no pets (46%). One-fourth of the respondents lived in Surigao City for 20-29 years, and three-fifths of them had a personal income of less than 5000 PhP per month. About three-fifths of the households had a combined income of less than 15,000 PhP per month.

Nested Logit Model
The multinomial logit model (MNL) is widely applied for discrete choice modeling because of its simplicity, ease of estimation, availability of estimation software, and potential to add new alternatives [31]. However, it implies proportional substitution across alternatives due to its performance of independence from irrelevant alternatives (IIA) [31]. The nested logit model can handle this limitation. The nested logit model is generalized from the standard logit model, which was initially proposed by McFadden in 1978 [32]. 'Nests' indicate that some alternatives may be joined in several groups or clusters. Error terms in the same nest present some correlation, but error terms of different nests are not correlated. The nested logit model is a member of the generalized extreme value (GEV) family that provides a variety of substitution patterns. The nested logit model has been applied in various fields, such as energy, transportation, housing, and telecommunication, among others [31]. We applied the nested logit model for statistical data analysis in this study because it is preferable to partition the finite set of alternatives into subsets, called 'nests'. The probability of a chosen alternative by an individual for a two-level nested logit model can be written as [31]: where i (i ∈ B k ) is the index of a chosen alternative; B k is a chosen nest; λ k is a measure of the degree of independence in observed utility among the alternatives in nest k, and K is the finite number of nests. j is a set of alternatives in a nest B l . The measure of independence is used for a measure of correlation (1−λ k ). V ni (V ni = β i x ni ) is the observed term of a utility function of a chosen alternative i made by an individual n. x ni is a column vector of variables (including a constant, earthquake-related damage, dwelling utilities and assets, and socioeconomic characteristics), and β i is a column vector of the corresponding coefficients. The nested logit probability of a chosen alternative can be written as a product of two standard logit probabilities: a marginal probability (the probability of a chosen nest B k ) and a conditional probability (the probability of a chosen alternative i within a nest B k ), as expressed below [31]: where W nk W nk = β k x nk is a variable that varies over a nest k, but not over alternatives within the nest. β k is a column vector of parameter estimates of a nest k. I nk is the inclusive value or inclusive utility of nest B k . The inclusive value (I nk = ln ∑ j∈ B k e Y nj +λ k ) is used to link the upper and lower models by bringing information from the lower model to the upper model. A measure of the degree of independence λ k can be excluded if all the coefficients are not generic (common) over nests because λ k in each nest is used to differentiate coefficients over nests [33]. λ k was excluded in our study because all variables are generic for all alternatives. Therefore, all parameter estimates are alternative-specific.   The explanatory variables of household car number, household motorcycle number, household multicab number, household truck number, the number of other vehicles in the household, household size, the number of children aged less than 5 years, the number of elderly people aged at least 65 years, the number of people with disabilities, the number of pets, individual income, and household income were arranged as the continuous variables.
The other variables were arranged as categorical variables. The dummy coding technique was used to investigate the impacts of the category variables on the alternatives.

Earthquake Evacuation and Duration Choice
The three-level nested logit model of evacuation and duration choice was developed and then used to estimate the percentage shares of alternatives. The estimated percentage shares of the alternatives are listed in Table 8. The estimated and actual values were the same for the alternatives of stay home, stay in an open space for less than 6 h, and stay at a public shelter for less than one day. The estimated and actual values were comparable for the other alternatives. The model estimation results of the three-level nested logit model are presented in Table 9. The McFadden R 2 of the model was 0.009, which is much lower than a decent value. The developed model did not fit the data well. This could be due to the inclusion of the sub-nest models at the bottom level (evacuation durations at a public shelter and R programming scripts were written to estimate the parameters of the two nested logit models. The package "maxLik" was used to estimate the parameters based on the standard maximum likelihood technique [34]. Two estimation approaches were used to estimate the parameters: sequential and simultaneous estimation approaches. The former is consistent, though not as efficient as the latter [31]. The former is used when problems arise in the simultaneous estimation approach, and parameter estimates in the sequential estimation approach are used as the starting values in the simultaneous estimation approach [31]. The forward selection approach was applied to deliberately remove parameters that were not statistically significant at the 0.10 significance level.

Earthquake Evacuation and Duration Choice
The three-level nested logit model of evacuation and duration choice was developed and then used to estimate the percentage shares of alternatives. The estimated percentage shares of the alternatives are listed in Table 8. The estimated and actual values were the same for the alternatives of stay home, stay in an open space for less than 6 h, and stay at a public shelter for less than one day. The estimated and actual values were comparable for the other alternatives. The model estimation results of the three-level nested logit model are presented in Table 9. The McFadden R 2 of the model was 0.009, which is much lower than a decent value. The developed model did not fit the data well. This could be due to the inclusion of the sub-nest models at the bottom level (evacuation durations at a public shelter and an open space). The first and second columns contain the explanatory variables and their corresponding categories/description. The third column contains the parameter estimates of the top model for the leave home alternative, and the stay home alternative was used as the reference. The fourth column contains the parameter estimates of the middle model for the public shelter alternative, and the open space alternative was used as the reference. The last two columns contain the parameter estimates of the bottom models. Stay at an open space for less than 6 h and a public shelter for less than 1 day were used as the references for the open space and public shelter models, respectively. The parameter estimates can be interpreted as follows. The intercept coefficients were included in the model to capture the average unobserved effect; there was no interpretable meaning [31].    As hypothesized, earthquake-related damage variables affected the decision on evacuation choice and duration. It is worth noting that earthquake warnings via telecommunications led residents to leave home and evacuate to a public shelter rather than an open space.
However, evacuees to an open space were more likely to stay for at least 6 h. However, this factor had no impact on the decision to stay for a certain duration at a public shelter. It is interesting to note that residents warned via the megaphone were more prone to evacuate to a public shelter rather than an open space. Evacuees to the open space, however, were more likely to stay for at least 6 h. Similar to residents informed via telecommunications, the megaphone method did not affect the decision of staying for a certain duration at the public shelter.
Coefficients of the home damage levels caused by the earthquake indicated that residents were more prone to leave home when it was severely damaged. Under a severe home damage case, people had a higher baseline likelihood to stay at the open space for at least 6 h and the public shelter for at least one day. A similar finding showed that the housing damage factor had a considerable impact on evacuation in rural areas [14]. Severe home damage was not statistically significant for the decision on evacuation location type. Under a moderate home damage case, residents were more likely to evacuate to the public shelter than the open space, and those evacuating to the open space were more likely to stay for longer than 6 h. However, this damage level did not affect the duration decision of staying at the public shelter. A similar study using the stated preference (SP) data showed that home damage intensity was positively associated with evacuation choice, and that public shelters were more preferable to open spaces [12].
Coefficients of the electricity and water supply cutoff implied that people preferred leaving home when electricity and water access were cut off due to the earthquake. However, the electricity cutoff did not affect the decision on evacuation location type and evacuation duration. For the water supply cutoff case, residents were likely to evacuate to the open space rather than the public shelter and stay in the open space for at least 6 h. Similar to the electricity supply cutoff factor, the water supply cutoff did not influence the decision on the duration of staying at the public shelter.
The availability of home utilities and assets, home ownership, number of floors of the house, genset availability, number of motorcycles, number of cars, number of multicabs, and number of trucks were not statistically significant at the 0.1 significance level for all the alternatives, nests, and sub-nests. Residents living in a timber-made house evacuated to the public shelter and preferred to stay for at least one day. However, the timber house did not affect the decisions on leaving or not leaving home, the evacuation location type, and the duration of staying in the open space. Residents living in houses with safe water and electricity access were likely to leave home when there was an earthquake. People residing in houses with safe water access were likely to evacuate to open spaces rather than to public shelters, but this factor was not statistically significant for the decision on evacuation duration. Electricity access did not affect the decision on evacuation location type and duration. It is worth noting that the presence of a toilet at home influenced the residents' decision of staying home rather than leaving home. However, evacuees were likely to evacuate to a public shelter rather than an open space. The availability of a toilet at home did not affect the decision on evacuation duration for the two evacuation location types. The presence of other vehicles in the household discouraged residents to leave home when there was an earthquake, but it did not influence the decision on evacuation duration for residents leaving home. Contrary to the tsunami evacuation behavior observed in a Cascadia Subduction Zone City, vehicle disability and lack of transportation modes discouraged residents to evacuate in an event of a tsunami [17].
In terms of sociodemographic characteristics, it is worth noting that gender, marital status, educational background, the presence of elderly people and pets, residential period, and individual income were not statistically significant for all the alternatives, nests, and sub-nests at the 0.10 significance level. The influences of the socioeconomic characteristics on evacuation choice and duration were less than the earthquake-related damage factors. A similar study empirically found that sociodemographic factors had lower impacts than subjective perception factors (i.e., built-up environments and disaster risk perception) on earthquake evacuation behavior, and gender and marital status were not statistically significant at the 0.05 significance level [14]. Contrary to the finding of Ao et al., residents with higher-education backgrounds tended to escape the building to seek refuge [14]. The 18-29 age group was less likely to evacuate to a public shelter than an open space, and those moving to the public shelter were likely to stay at least one day. Other age groups were not statistically significant at the 0.10 significance level. In terms of occupation type, self-employed residents were less likely to move to a public shelter than an open space and were likely to stay at a public shelter longer than one day. Household heads were more likely to leave home, and residents with the presence of more household members preferred an open space. Residents evacuating to a public shelter were likely to stay for at least one day. The presence of children aged less than 5 years in the household discouraged residents to evacuate, probably due to difficulties in moving small children. Households with the presence of PWDs pushed family members to evacuate to an open space and stay there for at least 6 h. People with high household income were more likely to stay at the public shelter for at least one day.

Earthquake Evacuation and Travel Mode Choice
The streamlined model of earthquake evacuation and travel mode choice based on the two-level nested logistic regression was developed. The developed model was then used to estimate the percentage shares of earthquake evacuees and travel mode types used and compare them with the actual ones, as can be seen in Table 10. The estimated and actual values were the same for the alternatives of stay home and leave home by walking and were comparable for the alternatives of leave home by two-and three-wheeler and leave home by the other modes. The model estimation results of the two-level nested logit model are presented in Table 11. The McFadden R 2 of the model was 0.481, which implied that the developed model provided a good fit for the data. The first and second columns contain the explanatory variables and their corresponding categories/descriptions. The third column contains the parameter estimates of the upper model for the leave home alternative, and the choice of stay home was used as the reference. The last two columns contain the parameter estimates of the lower model for the walking and two-and three-wheeler alternatives, respectively. The choice of the other modes was used as the reference. The parameter estimates can be interpreted as follows. Intercept coefficients have no interpretable meaning.
For earthquake-related damage variables, residents were more likely to evacuate upon receiving an earthquake warning, and warning via telecommunication was more influential than through megaphone/word-of-mouth. The telecommunication warning method did not influence the travel mode choice for evacuation, but residents who were warned via the megaphone/word-of-mouth were least likely to travel by walking. Residents were more likely to evacuate in the case of severe/complete home damage, followed by moderate damage. When homes were severely damaged, residents were least likely to evacuate by walking. Residents were more prone to evacuate when water and electricity supply were cut off. The impact of the electricity supply cutoff was statistically greater than the water supply cutoff on the decision on evacuation, i.e., stay home or leave home. In the case of water supply cutoff, the residents were most likely to leave home by two-and three-wheeler.
Home utilities and assets, home ownership, home-building material, electricity access, availability of genset and toilet, and the number of household cars, multicabs, and trucks were not statistically significant at the 0.10 significance level in affecting evacuation and travel mode choice. Residents living in one-story homes preferred not to evacuate. The number of floors did not affect travel mode choice for an earthquake evacuation. Residents with safe water access were likely to evacuate. It is surprising to note that the presence of household motorcycles and other vehicles discouraged people to evacuate but did not affect the travel mode choice decision for the earthquake evacuation. Table 11. Earthquake evacuation and travel mode choice coefficients (standard error).   Sociodemographic characteristics of gender, household size, number of children aged less than 5 years, presence of elderly people and pets, residential period, and individual and household incomes were not statistically significant at the 0.10 significance level in affecting the evacuation and travel mode choice. Contrary to tsunami evacuation behavior observed in Cascadia Subduction Zone City, evacuees with high income were more likely to travel by vehicle than walking, but females had a higher baseline likelihood than males to evacuate by foot [17]. Married people were least likely to evacuate by walking. The 30-39 and 50-59 age groups were more likely to evacuate than the other age groups. The 40-49 age group was most likely to evacuate by two-and three-wheeler, followed by walking. Ao et al. showed that the impact of the age factor on earthquake evacuation choice was not statistically significant [14]. A study of earthquake and tsunami evacuation behavior revealed that older people were more prone to evacuate by foot [17]. Residents with a high school degree were more likely to stay home. Employees and household heads were most likely to evacuate by walking. The presence of PWDs in households discouraged the residents to evacuate, probably due to difficulties in moving. The same finding was also confirmed: that households with the presence of PWDs may decide against evacuating [17]. Household evacuees with PWDs were most likely to travel by two-and threewheeler, followed by walking.

Earthquake Disaster Management and Planning
To mitigate chaos and damage cost and to enhance safety, it is important to understand the determinants of evacuation choice, which contribute to making policy and planning of areas and locations that are best suited for evacuation. The modeling results of the two nested logit modes informed us of the potential factors of evacuation choice, evacuation duration, and travel mode choice. The results showed that house damage levels and household utility cutoff (e.g., safe water supply and electricity) noticeably influenced earthquake evacuation decisions. However, these modeling results cannot rank geographic criteria for the selection of earthquake evacuation locations. The third stage of the evacuation process is knowing where to evacuate. To further support disaster management and planning, AHP was applied to understand the ideas of local government officials and planners. We considered five factors (criteria), i.e., tsunami inundation, terrain elevation, proximity to the road network, distance to the nearest fault line, and terrain slope. The factors and their corresponding sub-criteria, score, and weight are listed in Table 12. The actual values of the factors for every evacuation area were extracted from GIS and the multiple relevant government agencies. The importance weights of the factors were derived from the assessment of earthquake preparation and resource. Because Surigao City is a coastal city, it is not surprising that the important weight of the tsunami inundation factor appeared to be the highest (35%). It is followed by terrain elevation (19%), proximity to the road network (17%), distance to the nearest fault line (16%), and terrain slope (13%). These important weights were used to support designing the areas best suited for the management and planning of earthquake evacuation. Figure 5 illustrates the identified areas for evacuation based on the criteria of various factors. Areas are highlighted according to suitability: green (most suitable), yellow (suitable), and red (not suitable). The categorized areas for earthquake evacuation based on the terrain elevation are shown in Figure 5a. About 75% of locations are most suitable for earthquake evacuation. According to the existing fault line factor as illustrated in Figure 5b, 99% of locations belong to the suitable and most suitable categories for evacuation. In terms of proximity to the road network, 98% of locations fall within the suitable and most suitable categories for earthquake evacuation, as shown in Figure 5c. In terms of the terrain slope factor, which can be seen in Figure 5d, the majority of locations are most suitable for earthquake evacuation. Figure 5e illustrates areas categorized for earthquake evacuation based on the tsunami inundation. The majority (79%) of locations are most suitable for evacuation, and the rest are not suitable. Figure 5f illustrates the location categories for evacuation based on the combined factors. The majority (84%) of locations are suitable or most suitable for evacuation. The illustrated location categories for earthquake evacuation can assist policymakers, planners, and the local authorities to manage evacuation sites and develop guidelines and practices to minimize chaos, damage cost, and risks of injuries and fatalities in Surigao city, Philippines. (e) (f)

Summary of Findings and Recommendations
The Philippines is a developing archipelagic country susceptible to earthquakes. Surigao City was hit by a 6.7 M earthquake in 2017, which caused the third-highest damage cost in the country since 1990. Lessons and experiences learned from the 6.7 M Surigao earthquake are helpful to manage and plan strategies to minimize chaos, damage cost, and loss of lives. This study investigated the determinants of earthquake evacuation choice and duration and travel mode choice during an evacuation. This study took into account earthquake-related damage, dwelling utilities and assets, and socioeconomic characteristics as the input variables. This study surveyed 1055 residents who experienced the 6.7 M Surigao earthquake based on the revealed preference approach. The survey results indicated that 62% of residents decided to evacuate, while 52% moved to open spaces and 10.14% moved to public shelters. About half of the respondents traveled by walking because the majority of residents moved to open spaces immediately after the earthquake and stayed there for less than 6 h (37%).
Nested logit modeling was applied to explore potential determinants of earthquake evacuation choice, duration, and travel mode. The results are concluded as follows. As hypothesized, the earthquake-related damage variables significantly affected the evacuation choice and duration and had greater impacts than housing utilities and assets and socioeconomic characteristics. It is interesting to note that the residents were likely to evacuate upon warning through telecommunication and preferred public shelters to open spaces. Locals were more prone to evacuate when their houses incurred moderate and severe/complete damage and when the electricity and water supply were cut off. Residents were most likely to walk, as it is conceivable that the majority of residents moved immediately after the earthquake struck to open spaces near their houses/locations (with shorter distances). Employees and household heads appeared to be most likely to evacuate by walking. On the other hand, people with a high school degree were found to be less likely to evacuate their homes. Similarly, the presence of PWDs in households discouraged the residents to evacuate; this is likely due to difficulties in moving. The residential period was not statistically associated with evacuation choice and duration and travel

Summary of Findings and Recommendations
The Philippines is a developing archipelagic country susceptible to earthquakes. Surigao City was hit by a 6.7 M earthquake in 2017, which caused the third-highest damage cost in the country since 1990. Lessons and experiences learned from the 6.7 M Surigao earthquake are helpful to manage and plan strategies to minimize chaos, damage cost, and loss of lives. This study investigated the determinants of earthquake evacuation choice and duration and travel mode choice during an evacuation. This study took into account earthquake-related damage, dwelling utilities and assets, and socioeconomic characteristics as the input variables. This study surveyed 1055 residents who experienced the 6.7 M Surigao earthquake based on the revealed preference approach. The survey results indicated that 62% of residents decided to evacuate, while 52% moved to open spaces and 10.14% moved to public shelters. About half of the respondents traveled by walking because the majority of residents moved to open spaces immediately after the earthquake and stayed there for less than 6 h (37%).
Nested logit modeling was applied to explore potential determinants of earthquake evacuation choice, duration, and travel mode. The results are concluded as follows. As hypothesized, the earthquake-related damage variables significantly affected the evacuation choice and duration and had greater impacts than housing utilities and assets and socioeconomic characteristics. It is interesting to note that the residents were likely to evacuate upon warning through telecommunication and preferred public shelters to open spaces. Locals were more prone to evacuate when their houses incurred moderate and severe/complete damage and when the electricity and water supply were cut off. Residents were most likely to walk, as it is conceivable that the majority of residents moved immediately after the earthquake struck to open spaces near their houses/locations (with shorter distances). Employees and household heads appeared to be most likely to evacuate by walking. On the other hand, people with a high school degree were found to be less likely to evacuate their homes. Similarly, the presence of PWDs in households discouraged the residents to evacuate; this is likely due to difficulties in moving. The residential period was not statistically associated with evacuation choice and duration and travel mode choice. It is surprising to observe that the availability of household vehicles for all vehicle types did not show a significant impact on the travel mode choice. It is also worth noting that the personal income factor did not affect earthquake evacuation choiceand travel mode choice. After interpreting the modeling results and discussion, AHP was applied to develop decision-making criteria for ranking locations and facilities to support earthquake management and planning. The importance weight of tsunami inundation became the highest (35%), followed by terrain elevation (19%), proximity to the road network (17%), distance to the nearest fault line (16%), and terrain slope (13%). According to the combined factors, about 84% of locations are suitable and most suitable for earthquake evacuation in Surigao city, Philippines. The location categories for earthquake evacuation will help policymakers and planners as well as the local authorities to manage evacuation sites and develop plans and guidelines for future earthquakes in the city.
Even though this study provided some novel findings to deal with sustainability and resiliency issues, there were some inherent limitations. First, the study used the selfreporting approach that cannot avoid a lack of bias-free, honest, and accurate reporting. Second, the study did not include psychological factors such as disaster risk perception and subjective norms. Third, the authors did not consider what should constitute significant supplies for evacuees as important factors to develop evacuation planning. Fourth, the study used the revealed preference data based on the respondents' experience with the 2017 Surigao earthquake that struck during the nighttime at 10:03 p.m.. The impacts of determinants of earthquake evacuation during the nighttime might be different from the daytime. Future studies should focus on the exploration of the difficulties that evacuees encountered and the primary supplies they need during evacuation. Furthermore, the impacts of earthquake evacuation guidance and signs on evacuation choice and duration and travel mode choice should be investigated to make evacuation planning more exhaustive and compelling. Route choice and evacuation facility recommendations should not be ignored. The impacts of psychological, social, and risk perception factors on human evacuation behavior should be taken into account for future research. Analysis of satisfaction with the authorities' services should be considered as a key factor to improve the effectiveness of management and planning for future earthquake evacuation to reduce damage cost and the number of injuries and deaths. Funding: This research received funding from the Engineering Research for Development and Technology (ERDT) under the Department of Science and Technology (DOST), Philippines. Also, special thanks go to DLSU Science Foundation for the publication support.