Simulation-Based and Risk-Informed Assessment of the Effectiveness of Tsunami Evacuation Routes Using Agent-Based Modeling: A Case Study of Seaside, Oregon

Zhenqiang Wang · Gaofeng Jia

Abstract Typically,tsunami evacuation routes are marked using signs in the transportation network and the evacuation map is made to educate people on how to follow the evacuation route.However, tsunami evacuation routes are usually identified without the support of evacuation simulation, and the route effectiveness in the reduction of evacuation risk is typically unknown quantitatively.This study proposes a simulation-based and risk-informed framework for quantitative evaluation of the effectiveness of evacuation routes in reducing evacuation risk.An agentbased model is used to simulate the tsunami evacuation,which is then used in a simulation-based risk assessment framework to evaluate the evacuation risk.The route effectiveness in reducing the evacuation risk is evaluated by investigating how the evacuation risk varies with the proportion of the evacuees that use the evacuation route.The impacts of critical risk factors such as evacuation mode (for example, on foot or by car) and population size and distribution on the route effectiveness are also investigated.The evacuation risks under different cases are efficiently calculated using the augmented sample-based approach.The proposed approach is applied to the riskinformed evaluation of the route effectiveness for tsunami evacuation in Seaside,Oregon.The evaluation results show that the route usage is overall effective in reducing the evacuation risk in the study area.The results can be used for evacuation preparedness education and hence effective evacuation.

Keywords Agent-based model · Risk-informed evaluation · Route effectiveness · Simulation-based framework · Tsunami evacuation risk · Tsunami preparedness education

1 Introduction

The near-field earthquake-induced tsunami can result in substantial loss of life (Suppasri et al.2012, 2013).Evacuation to safety zones(for example,shelters)is regarded as one of the most effective ways to survive a tsunami (Mas et al.2012; Takabatake et al.2017).Proper evacuation planning is crucial to support effective evacuation and reduce the evacuation risk(UNESCO 2007).To produce an effective evacuation plan for evacuation risk management,many issues need to be considered such as a good understanding of the evacuation behavior, an accurate evacuation risk assessment, and effective risk mitigation (Wang 2021).For effective risk mitigation, many strategies could be adopted(Leo′n and March 2016)such as extensive civilengineered defenses (for example,seawalls), land use, and built environment measures (for example, building codes),and emergency readiness systems(for example,evacuation preparedness).Tsunami preparedness education is believed to effectively raise people’s awareness of tsunami risk,improve the evacuation,and ultimately reduce the casualty in tsunami-prone communities (L?vholt et al.2019).For instance, the huge death tolls in the 2004 Indian Ocean tsunami were not seen in the 2011 Tohoku event, which can be partially attributed to the good preparedness of the population in some coastal communities in Tohoku (L?vholt et al.2019).Preparedness education programs are held to increase people’s awareness of the future tsunami(Adiyoso and Kanegae 2012; City of Rockaway Beach,Oregon 2019) and familiarity with facilities that are designed for evacuation (Adiyoso and Kanegae 2012).To raise people’s tsunami awareness, community members(even the tourists) can be guided to visit tsunami museum,stranded vessel,and graveyard that are assigned as tsunami disaster educational facilities by the government (Adiyoso and Kanegae 2012).In addition,volunteer educators can go door-to-door to discuss tsunami awareness and safety with residents, hand out evacuation maps or brochures to residents(City of Rockaway Beach,Oregon 2019).To become familiar with evacuation facilities, people can visit the evacuation route and shelter building (Adiyoso and Kanegae 2012).People would easily find such facilities and evacuate to the shelter by following the evacuation route once a tsunami occurs.

Among various preparedness education programs,making tsunami hazard and evacuation maps is recommended to be the best way to start preparedness efforts(Dwelley Samant et al.2008).Evacuation routes are marked in the evacuation map for the public to make their evacuation plans (that is, route selection) (Kurowski et al.2011; DOGAMI 2013).In the evacuation, evacuees are encouraged to move to the nearest road and follow the route with evacuation signs to an assembly area, for example, horizontal or vertical shelter (Scheer et al.2011;DOGAMI 2013).However,the evacuation route usage can be highly uncertain due to the variability in the characteristics of evacuees.Different types of evacuees (for example, residents and visitors) may have different knowledge, perception, or even experience on the evacuation route, which would impact the evacuation decision and behavior,especially the path selection(Limanond et al.2011; Scheer et al.2011; Takabatake et al.2017).Typically, not all residents know their local evacuation routes(Dengler 2005),not to mention tourists.Learning about the evacuation map can make people (especially tourists)familiar with the evacuation route and effectively evacuate when the tsunami disaster strikes.Furthermore, some points need to be taken into account for the preparedness education about evacuation routes.Residents may have different views about the potential usefulness of evacuation routes (Dall’Osso and Dominey-Howes 2010).For example, some of the residents are concerned that evacuation routes will be blocked when many people follow the signs simultaneously(Bernard 2005).All of these can impact the evacuees’evacuation behavior(whether use the evacuation route or not,that is,there is uncertainty in the proportion of the evacuees that use the route), and ultimately impact the route effectiveness in reducing the evacuation risk in terms of casualty.

To gain high route effectiveness in reducing the evacuation risk through tsunami preparedness education(that is,educate people to effectively use the evacuation route),some important questions about the route effectiveness need to be answered.For example,does it always lead to a more effective evacuation by following the evacuation route compared to not following the route (for example,following the commonly used shortest path(Imamura et al.2012; Yamamoto and Li 2017))? In what cases, is following the route more effective? Do the evacuation mode(for example, on foot or by car), population size, and population distribution have a high impact on the route effectiveness? If so, how do they impact the route effectiveness?

The answers to the above questions associated with the effectiveness of route usage may be found through understanding how the evacuation map is made, and evacuation routes are identified.Evacuation routes are marked in the evacuation map or brochure (in print or online) for the public to make the evacuation planning (that is, which route to follow in evacuation) (Kurowski et al.2011;DOGAMI 2020a).Besides, web-based map interface by revising tsunami evacuation maps or interactive evacuation map viewer is available for the residents and visitors to custom design their own maps (DOGAMI 2020b).Evacuation maps for the at-risk communities along the coast are typically developed based on the state-developed inundation maps of a credible worst-case tsunami(Darienzo 2001;Gonza′lez et al.2001; Priest et al.2001; Gonza′lez et al.2005; Tanioka et al.2012).The inundation maps are produced through numerically efficient hydrodynamic model simulations with adequate, site-specific bathymetry and topography data (Gonza′lez et al.2001, 2005).Evacuation routes are typically identified by local and state officials,and reviewed by the department that develops the inundation maps (Gonza′lez et al.2001, 2005; Bernard 2005;Kurowski et al.2011),or determined in close collaboration with communities and with feedback from tsunami experts(Bernard 2005).Some efforts exist in the search of highpriority tsunami evacuation routes for more effective evacuation (NTHMP 2018) or in the identification of the optimal evacuation route using various advanced algorithms, for example, the modified shortest path algorithm(Shekhar et al.2012), the algorithm using GIS techniques(Schuster and Gomez 2013; Ai et al.2016; Trindade et al.2018), graph theory-based evacuation routing algorithm(Pe′roche et al.2014; Patel et al.2016), and discrete optimization algorithm (Forcael et al.2014), and so on.However,to the best of our knowledge,the development of the tsunami evacuation map and identification of the evacuation route lack the support from evacuation simulation,and few efforts are put in the quantitative evaluation of the existing tsunami evacuation route effectiveness(that is, investigating the impact of route usage on the evacuation risk based on tsunami evacuation simulation).In this context, it is difficult or even impossible to implement aneffective tsunami preparedness education on how to use the evacuation route effectively and hence improve the evacuation.

To address the above limitations, this study proposes a simulation-based and risk-informed approach to quantitatively evaluate the evacuation route effectiveness in reducing the evacuation risk for tsunami preparedness education.An agent-based model is used to simulate the tsunami evacuation, based on which the evacuation risk is quantified and assessed.The route effectiveness is evaluated by investigating how the evacuation risk varies with the proportion of the evacuees that use the evacuation route(that is,conditional evacuation risk assessment).The route effectiveness is also evaluated by investigating the impact of route usage on the evacuation risk under different values of critical risk factors (that is, the factors that have high impacts on the evacuation risk such as the proportion of the population that evacuate by car in the multi-modal evacuation, population size, and population distribution).The conditional evacuation risks are efficiently calculated using the augmented sample-based approach.The proposed approach is applied to the risk-informed evaluation of the route effectiveness for tsunami evacuation in Seaside,Oregon.However, the proposed approach is general and applicable to the evaluation of the route effectiveness in other regions with high tsunami risk around the world.Accordingly,the evaluation results can be used to improve the tsunami preparedness education(for example,guiding a more effective route usage) and hence the tsunami evacuation.

The remainder of this article is organized as follows.Section 2 describes the agent-based tsunami evacuation model for evacuation simulation.Section 3 presents the evaluation of route effectiveness through conditional evacuation risk assessment.Section 4 describes an illustrative example.Section 5 presents the results and discussions.Section 6 summarizes the research findings.

2 Agent-Based Tsunami Evacuation Model

For tsunami evacuation simulation,we use the agent-based evacuation model developed in Wang and Jia (2021a),which extends on the agent-based modeling (ABM)framework for a multi-modal near-field tsunami evacuation simulation proposed in Wang et al.(2016).The agentbased model used here is briefly reviewed.

The agent-based tsunami evacuation model (Fig.1)consists of the evacuation environment model (EEM),evacuation decision and behavior model (EBM), and evacuation performance model (EPM).The EEM includes the multi-hazard model (that is, seismic and tsunami hazards), transportation network, tsunami shelter, and population distribution.Any seismic event that can trigger a tsunami and the resulting tsunami hazard can be incorporated into the multi-hazard model.The traffic capacity of the damaged link(that is,road/bridge)is reduced based on the seismic damage state predicted from fragility curves.The time histories of the tsunami inundation are generated using ComMIT/MOST (Titov and Gonza′lez 1997).Different population sizes and distributions throughout the community are modeled to characterize the population mobility.The population size (including the residents and tourists) is estimated considering the actual residential population and the population mobility.The population distribution can be of high spatiotemporal variability(Mostafizi et al.2017).To model the uncertain population distribution throughout the community, different probability models such as the normal distribution and uniform distribution can be used (Mostafizi et al.2017).In evacuation simulation, the population is limited to the hazard area.In the EBM,the evacuation decision and behavior are modeled.The evacuation decision includes the departure time, the multi-modal evacuation (that is, evacuation on foot and by car), and the path selection (for example, following the evacuation route or using the shortest path)throughout the transportation network.For the evacuation behavior, individual behavior, interactions between evacuees, and interactions between evacuees and the damaged environment are modeled.For individual evacuation behavior,the variability in the preferred pedestrian speed is modeled considering different characteristics of pedestrians.The speed is adjusted according to the corresponding traffic density ahead for both the pedestrian and car.The interactions between evacuees include those occurring between pedestrians or cars (for example, collision avoidance through speed adjustment), and between the pedestrian and car (for example, the pedestrian-vehicle interaction is modeled by traffic stage transition where the traffic stage is determined by the volume ratio of pedestrians to vehicles).In addition, evacuees would interact with the damaged environment.For example, the traffic capacity reduction of the damaged link would cause traffic congestion and hence slow down the evacuees’ travel speed,which in turn impacts the traffic congestion level of the damaged link.The evacuees would need to reroute when some road is completely blocked or the bridge is completely damaged, which in turn impacts the traffic condition on the damaged transportation network.These dynamic interactions between evacuees and the damaged environment are also modeled.In the EPM,the quantity of interest such as the casualty rate (that is, the proportion of the number of casualties to the total population)is selected as the evacuation performance measure and the critical water depth is used to determine the casualty (that is, theevacuee is considered to become a casualty when the inundation depth exceeds the critical water depth).

Fig.1 Illustration of the agent-based tsunami evacuation model

The agent-based tsunami evacuation model is developed in NetLogo considering that the integrated ABM environment can use modern data with a high level of detail,capture the emergent phenomena, and tune the complexity of the agent’s behavior (Wilensky 2001).

3 Evaluation of Evacuation Route Effectiveness through Conditional Tsunami Evacuation Risk Assessment

To evaluate the evacuation route effectiveness(that is,how the evacuation route usage impacts the evacuation risk),the evacuation risk under different proportions of the evacuees that use the evacuation route is evaluated.In addition, the impacts of critical risk factors such as evacuation mode(for example, on foot or by car) and population size and distribution on the route effectiveness are also investigated.

3.1 Quantification of Tsunami Evacuation Risk

To quantify the tsunami evacuation risk, the simulationbased framework in Wang and Jia(2021b)is used.Figure 2 shows a flowchart for the risk quantification.In the evacuation risk quantification, various uncertainties associated with the evacuation process (including the aleatory and epistemic uncertainty) are considered, such as those associated with the seismic damage to the transportation network in the EEM, the evacuation decisions and behaviors of evacuees(for example,departure time)in the EBM,and the critical water depth in the EPM.Let x=[xe,xb,xp]= [x1, …, xi, …, xnx] ∈X represent all the input random variables in the space of X,where xe,xb,and xpdenote the input random variables in the EEM, EBM, and EPM,respectively, xiis the ithinput random variable, and nxis the total number of input random variables.The probability model p(x|θ)is used to quantify the uncertainty in x,which corresponds to probability density function (PDF) for continuous variables or probability mass function (PMF)for discrete variables.Here θ ∈Θ represents the distribution parameters that define the probability model p(·|θ),and Θ is the space of all potential values of θ.When distribution parameters θ are uncertain, the probability distribution p(θ) can be used to quantify the uncertainty in θ,which can be used to represent the epistemic uncertainty in the probability model for x.Based on the evacuation simulation, the evacuation performance can be evaluated for any given realization of input x.Let h(x) represent the performance measure of the evacuation system for given x.We can define h(x) as the risk consequence measure.Depending on the quantity of interest, metrics such as the evacuation time needed to get a certain proportion of the population in the community to the safety zone,the number of people evacuated within a certain time frame, and the proportion of the number of casualties to the total population within a certain time frame (that is, casualty rate)could be selected as the evacuation risk consequence measure.When h(x) corresponds to the risk consequence measure(for example,the casualty rate),propagation of the uncertainties in x under given probability model p(x|θ)leads to quantification of the evacuation risk (that is, conditional evacuation risk) H(θ)

Fig.2 Quantification of tsunami evacuation risk.EEM = evacuation environment model,EBM=evacuation decision and behavior model,EPM = evacuation performance model, x = all input random variables, xe, xb, xp = input random variables in EEM, EBM, and EPM, respectively, X = space of all potential values of x, h(x) = risk consequence measure, θ = distribution parameters that define the probability model p(·|θ),Θ=space of all potential values of θ,p(x|θ)=probability model that is used to quantify the uncertainty in x under given θ, p(θ) = probability distribution that is used to quantify the uncertainty in θ, H(θ) = tsunami evacuation risk under given θ, H = tsunami evacuation risk

This corresponds to the evacuation risk under the given selection of distribution parameter θ.When θ is uncertain,further propagating the uncertainties in θ gives the evacuation risk H

3.2 Conditional Tsunami Evacuation Risk

To evaluate the effectiveness of route usage, an input variable pr, which corresponds to the proportion of evacuees that use the evacuation route, is defined in x to characterize the amount of route usage.Probability model has been assigned to quantify the uncertainty in pr.The evacuation risk H already considered the uncertainty in pr.Therefore,to evaluate the effectiveness of route usage,here the evacuation risk conditional on different prvalues (that is, conditional on different proportion of evacuees that use the evacuation route) is used and investigated.This involves the evaluation of the conditional evacuation risk H(pr).Or more generally, the risk conditional on given value of scalar input random variable xi, that is, H(xi),needs to be evaluated, and for H(pr) it simply corresponds to H(xi) for the case of xibeing pr.To evaluate H(xi), the uncertainties in x～iand θ need to be propagated, which leads to

where x～irepresents the remaining of the input random variables excluding xi.

To gain more insights into the effectiveness of route usage, the impacts of critical risk factors such as evacuation mode (for example, on foot or by car) and population size on the route effectiveness are also investigated.This is realized by investigating the variation of the evacuation risk with prunder different values of the critical risk factors.Such investigation also involves the evaluation of the corresponding conditional evacuation risk, which can be written in general form as H(xij)with xij=[xi,xj]=[pr,xj]where xjrepresents one of the risk factors mentioned earlier.Similarly, to evaluate H(xij), the uncertainties in x～ijand θ need to be propagated, which leads to

where x～ijdenotes the remaining of the input random variables excluding xij.

The above conditional evacuation risks (that is, H(xi)and H(xij)) correspond to high dimensional integral, and they can be approximated using stochastic simulation, for example, Monte Carlo simulation (MCS) (Robert and Casella 2004).However, direct application of MCS to estimate H(xi) requires large number of model evaluations and calculation of H(xi) for different values of xiwould require repeating the above process and entails significant computational challenges, especially considering the fact that the evacuation model is expensive to run.Similar challenges hold for the evaluation of H(xij)for different xijvalues.

3.3 Efficient Calculation of Conditional Evacuation Risk Using the Augmented Sample-Based Approach and Kernel Density Estimation

To efficiently estimate the above conditional evacuation risks, the augmented sample-based approach proposed in Wang (2021) is used, which extends on the sample-based approach introduced in Jia and Taflanidis (2016) for efficient estimation of Sobol’ index.In the augmented approach, an augmented problem in terms of [x, θ] is first defined.Then, a joint auxiliary distribution π(x, θ) with respect to [x, θ] is defined, which is proportional to the integrand of the evacuation risk integral H

based on which, the marginal auxiliary distribution for xican be written as

Then H(xi) can be written as

Based on Eq.7, H(xi) for any given value of xican be efficiently estimated by: (1) generate a set of candidate samples {[xk, θk] , k = 1,…,N} from some joint proposal density q(x, θ), and evaluate the corresponding h(x); (2)estimate the evacuation risk H using

The coefficient of variation (CoV, denoted δCoV) of the estimate can be also established; (3) pass these candidate samples through some stochastic sampling algorithm such as rejection sampling (Jia et al.2017) to generate samples from the joint auxiliary PDF π(x, θ), denoted｛［，］，k ＝1，...，nπ｝; (4) project the generated samples to the space of xi,based on which an approximation to π(xi)(denoted as（xi）) is established using kernel density estimation (KDE) or the boundary corrected KDE for the input random variables that are bounded(Jia and Taflanidis 2016); (5) estimate the conditional evacuation risk under any value of xiefficiently using

Similarly, the conditional risk H(xij) can also be evaluated efficiently using the same set of nπsamples.More specifically, projection of these samples to the xijspace gives samples from the marginal density π(xij), and based on these marginal samples, an approximation to π(xij) (denoted as（xij）) is established using the multivariate KDE (Jia and Taflanidis 2014, 2016).Then（xij）can be calculated under any given value of xij=[xi,xj] using

4 Illustrative Example

We perform the simulation-based and risk-informed evaluation of evacuation route effectiveness in Seaside, Oregon.The evacuation simulation corresponds to the case study in Wang and Jia (2021a); however, here two improvements are made to facilitate the investigation of route effectiveness.First, besides searching the shortest path in terms of distance,following the evacuation route is also considered an option for path selection in the evacuation model, and then the proportion of the evacuees who use the evacuation route is used to characterize the amount of route usage and quantify the route effectiveness.Second,instead of assuming the same cross-section for all links as in Wang and Jia(2021a),here the actual road cross-section is used for both the evacuation route and all other links in the transportation network.According to the road crosssection, roads with different widths are classified into different functional classifications (that is, ‘‘Principal Arterial,’’ ‘‘Minor Arterial,’’ ‘‘Major Collector,’’ and ‘‘Local Road’’), and maximum car speed limits are defined based on the different classifications.This improvement is made for a more realistic traffic simulation and hence for more accurate insights into the route effectiveness evaluation in the study area.

Next, the study area, Seaside, Oregon is described briefly.Then characteristics of the transportation network(especially those related to the road cross-section),tsunami evacuation route, and evacuation zone (especially those related to the evacuation route) are introduced in detail.Then implementation details for the evacuation simulation and the evacuation risk assessment are presented,including the case studies that are defined to investigate the route effectiveness.

4.1 Study Area

The coastal community of Seaside, Oregon (Fig.3b),which has been selected in many tsunami evacuationrelated studies (Wang et al.2016; Mostafizi et al.2017, 2019), is selected as the study area.Seaside is regarded as one community with a high tsunami evacuation risk considering: (1) the high risk to the seismic and nearfield tsunami hazards due to the proximity to the Cascadia Subduction Zone (CSZ, shown in Fig.3a) (Priest et al.2016; Mostafizi et al.2017); (2) the high vulnerability to tsunami inundation due to the fairly flat topography with two rivers flowing through the city, the relatively far shelters (more than 1.5 km), and the bottlenecks in the transportation network caused by the 10 bridges over the two rivers (Mostafizi et al.2017).

Fig.3 The case study area.(a) Cascadia Subduction Zone (CSZ) (USGS 2020), and (b) Seaside, Oregon

4.2 Transportation Network

Figure 4 shows the transportation network, which consists of 700 nodes and 896 links.The road cross-section standards corresponding to the federal functional classification of roads throughout Seaside (Clatsop County 2019) can be found in the typical design of streets in Clatsop County(City of Seaside 2010).Table 1 summarizes the values of the parameters used to determine the road width occupied by car and pedestrian (City of Seaside 2010).For a given functional class,w represents the total road width(unit:m);wwminand wwrepresent the minimum and maximum road width (unit: m) that can be used by the pedestrian; wrdenotes the minimum width of each vehicle lane(unit:m);nlrepresents the minimum number of vehicle lanes used.To model the pedestrian-vehicle interaction, different traffic stages are defined.More specifically, according to the volume ratio between pedestrians and vehicles, different traffic stages(that is,vehicle-dominated,balanced,and pedestrian-dominated) are dynamically selected to determine the road widths occupied by pedestrians and cars on the same link (Wang and Jia 2021a).

Table 1 The parameters used to determine the road width occupied by the car or pedestrian

The well-known Greenshields’ model of speed and density is used to simulate the car speed adjustment.Here the maximum car speed limit in the Greenshields’model is taken as vcm= 64.4, 56.3, 48.3, and 40.2 (unit: kph) for Principal Arterial, Minor Arterial, Major Collector, and Local Road (ODOT 2020), respectively, and the traffic capacity of undamaged links is taken as ρcm=200 vehicle/km based on the study in Mostafizi et al.(2017, 2019).

4.3 Tsunami Evacuation Route and Evacuation Zone

Following the approach in Lonergan et al.(2015), first the evacuation routes are identified using the evacuation map(DOGAMI 2013) and shapefiles (Claptsop County 2008).Note that one additional evacuation route that goes through the bridge over the river in downtown is added considering that the evacuees nearby would intuitively choose this route to go across the river (Hajo Neis et al.2015).The identified evacuation routes are shown in Fig.5.

Inspired by the evacuation map (DOGAMI 2013) and the beat-the-wave evacuation map for tsunami hazards in Seaside (Priest et al.2016), the evacuation zone and evacuation direction are created throughout the community(Fig.5) to determine the selection of the evacuation route by evacuees in each zone as well as the evacuation flow direction in the evacuation model.Note that in DOGAMI(2013) no evacuation zone is provided and therefore the evacuation direction is not designed for any specific evacuation zone.In Priest et al.(2016), both evacuation zone and the corresponding evacuation direction are created where the evacuation zones have complex shapes(determined from the least-cost distance calculation).In the current study,both the evacuation route and corresponding evacuation direction are established, but they areestablished in a way that is different from DOGAMI(2013)and Priest et al.(2016).In particular, from the point of view of capturing the more possible selection and use of the evacuation route, the current study takes into account the following considerations about the creation of the evacuation zone and use of the evacuation route: (1) the region located between two main evacuation routes is divided into two evacuation zones such that the distance between each zone (centroid) and the adjacent route is close; (2) in each evacuation zone, the evacuees who decide to use evacuation routes will first move to the closest evacuation route and then follow the route to the shelter;(3) the evacuees who do not follow the evacuation route would search the shortest path; (4) when the bridge on some route collapses, the evacuees who plan to use the evacuation route would search the shortest path for evacuation.

Fig.5 Tsunami evacuation routes and zones in Seaside,Oregon.The brown, pink, and green polygons represent the evacuation zones on the beach,in the downtown area,and in the residential area, respectively

Fig.4 The transportation network in Seaside, Oregon

4.4 Implementation Details

The historical seismic event in 1700 is selected as the hypothetical input of the multi-hazard scenario in which Mw=9.0 and the focal depth of 40 km are considered.This scenario corresponds to the maximum inundation area that is used in the evacuation map in Seaside, Oregon (Group 2006).We select the ground-motion prediction equations(GMPEs) in Campbell and Bozorgnia (2008) to estimate the permanent ground deformation(PGD)at each road site and spectral acceleration (Sa) at each bridge site, respectively.Then the fragility functions of the road and bridge in HAZUS (DHS 2009) are used to estimate the probabilities for seismic damages to roads and bridges, respectively.According to the damage level (that is, none, slight,moderate, extensive, and complete), the residual traffic capacity of the road/bridge is calculated based on the residual rate of the traffic capacity defined in Wang and Jia(2021a).It is assumed that the evacuation occurs in summer and lasts for one hour after the occurrence of the earthquake.One hour is used because the tsunami would reach the run-up limit at around 55 min and the evacuation performance would become stable within one hour.

As for the input random variables, xe= [xdi] with i = 1,2,…, 896 and xdicorresponds to discrete random variable that represents the five seismic damage states of the links(that is,xdi=DSjwith j=1,2,…,5).xb=[t0,t,pc,T1P,T2P,T1M,T2M,T1C,T2C,T1L,T2L,vp,pr],and xp=[hc],where t0is the time when receiving the tsunami warning(unit:min);t is the departure time(unit:min);pcdenotes the proportion of the evacuees that use the car; T1P, T2P, T1M, T2M, T1C,T2C, T1L, and T2Ldenote the thresholds that determine the traffic stage corresponding to the road class Principal Arterial,Minor Arterial,Major Collector,and Local Road,respectively;vprepresents the pedestrian speed(unit:m/s);pris the proportion of the evacuees that follow the evacuation route(the proportion of the evacuees that search the shortest path is 1-pr); hcrepresents the critical depth (unit:m).The key continuous input random variables in x and corresponding distribution parameters in θ are presented in Table 2,where U and TN represent the uniform distribution and truncated normal distribution, respectively.The selection of the probability distributions and values can be found in Wang(2021).For the ith distribution parameter θi,the subscripts ‘‘l’’ and ‘‘u’’ represent the lower and upper bounds for the ith input random variable xi.For example,lt0and ut0are the lower and upper bounds for t0,respectively.For other θi,τtrepresents the delay time after receiving the tsunami warning and σtis the scale parameter; μcand σcdenote the mean and standard deviation of pc,respectively;μpand σpare the mean and standard deviation of vp,respectively.The values/distributions for the same θiare separated with semicolon (;), for example, for θicorresponding to τtand σt, τtfollows U[0,5] and σtfollows U[1,5], and they are separated with semicolon, written as U[0, 5]; U[1, 5].

In terms of the risk consequence measure h(x), we use the casualty rate in the community within an hour.Here,we consider the pedestrian casualty rate(PCR),that is,that number of casualties that evacuate on foot divided by the total population, car casualty rate (CCR), that is, the number of casualties that evacuate by car divided by the total population, and total casualty rate (TCR), that is, the number of all casualties divided by the total population.Then,the evacuation risk corresponds to the expected PCR,CCR, and TCR, respectively.

It has been shown that some risk factors could have relatively high impacts on the tsunami evacuation risk,including the car use (that is, pc) (Wang 2021; Wang and Jia 2021b), population size (Wang and Jia 2021a, 2021b)and distribution (Wang 2021).Here, the evacuation route effectiveness is evaluated under uncertain pc(that is,H(pr))and given values of pc(that is,H(pr,pc)),both of which are performed under different populations (denoted as ne) and population distributions throughout the three sub-areas(denoted as ρ1, ρ2, and ρ3corresponding to the percentage of the population on the beach, in the downtown area, and in the residential area, respectively).Six cases (denoted as C1—C6) are defined in Table 3 to consider different population sizes and distributions.The up-to-date residential population of Seaside is 6795 in 2018 according to U.S.Census (U.S.Census Bureau 2018); however, the net daytime residential population can decrease to around half of the population (Sleeter and Wood 2006) while the number of tourists can be up to 10,000 per day in the peaksummer(Mostafizi et al.2017),that is,population mobility exists.We select three population sizes (that is, ne= 5000,ne= 10,000, and ne= 15,000) to represent the population mobility.Here, ne= 5000 and ne= 15,000 are selected to represent the population size in the early summer and peak summer, respectively, while ne= 10,000 is used to model the population size during some time between the early summer and peak summer(Wang and Jia 2021a).Since the distribution varies with the time of the day, so two scenarios are defined to reflect this, corresponding to daytime and nighttime scenarios.Under the daytime scenario, the population is distributed across the three sub-areas with ρ1=0.4,ρ2=0.3,and ρ3=0.3 based on the study in Mostafizi et al.(2017).For the nighttime scenario,a larger proportion of evacuees would be in the downtown and the residential areas while fewer people would be on the beach.Here,ρ1=0.1, ρ2= 0.5, and ρ3= 0.4 are used.The population distributions and other related characteristics of the daytime and nighttime scenarios are summarized in Table 4.Note that all characteristics are applicable for the evacuations on foot and by car except ‘‘Preferred pedestrian speed.’’ Due to the automotive lighting, the car speed is considered not to be affected during the night.Besides,the three speeds in the pedestrian speed-density model, that is, 0.75, 1.5, and 3.83(unit:m/s)(Wang and Jia 2021a),are used here for the daytime scenario;to model the impact of nighttime(that is,worse lighting condition and shorter sight distance compared to daytime) on the pedestrian speed-density model,these three speeds are reduced by 20% for the nighttime scenario according to the study in Ouellette and Rea(1989), that is, reduced to 0.6, 1.2, and 3.06 (unit: m/s),respectively.

Table 2 Key continuous input random variables and associated distribution parameters

Table 3 Definition of evacuation simulation cases

Table 4 Characteristics of daytime and nighttime scenarios

5 Results and Discussions

The augmented sample-based approach presented in Sect.3.3 is used to efficiently estimate H(pr)and H(pr,pc)under the six cases.For each case, we select q(x, θ) = p(x|θ)p(θ)to generate N = 5,000 candidate samples, based on which the corresponding h(x)(that is,the PCR,CCR,and TCR)is evaluated.Under the selection of q(x,θ)and N=5,000,all the δCoVfor the estimate of ^H(that is, the expected PCR,CCR, and TCR) are below 2%.For each case, using the information from the N = 5,000 candidate samples, more than nπ=300,1,500,and 1,700 samples are generated from π(x,θ)when h(x)corresponds to the PCR,CCR,and TCR,respectively.Overall, for the cases with larger population size ne,there are more samples from π(x,θ)generated as a result of the relatively high evacuation risk.

5.1 Evaluation of Route Effectiveness under 5,000 Evacuees

This section investigates the route effectiveness under 5,000 evacuees.In particular, the variations of the evacuation risk (that is, expected PCR, CCR, and TCR) with route usage and with both route usage and evacuation mode under the daytime and nighttime scenarios(that is, C1 and C4) are investigated.

5.1.1 Variation of Evacuation Risk with Route Usage

The variations of the expected casualty rates with the proportion of the evacuees that use the evacuation route(that is, pr) for the PCR, CCR, and TCR under C1 and C4 are shown in Figs.6a—c, respectively.

Daytime scenario (C1).As shown in Fig.6a, the expected PCR almost does not change when pris relatively small(that is, smaller than 0.35) while it decreases slightly as prkeeps increasing(that is,only decreases from 2.8%to 2.5%when princreases from 0.35 to 1).This means if more pedestrians follow the evacuation route it will not reduce the PCR by much, hence for evacuation on foot following the evacuation route is not very effective in the risk reduction in this case.This low route effectiveness should be attributed to the close impacts of following the evacuation route and using the shortest path on the PCR under the relatively small population.Compared to following the evacuation route that is relatively wide and is distributed relatively evenly throughout the community, using the shortest path for pedestrians tends to induce a higher concentration of pedestrians on some links, especially when these links are relatively narrow (for example,Local Road).The higher concentration of pedestrians would cause more severe pedestrian congestion and hence more casualties due to the delay.This congestion is more obvious when the number of pedestrians is large.However,under ne= 5,000, the number of pedestrians is overall relatively small,and the pedestrian congestion level would not be much higher when using the shortest path.In other words, more pedestrians following the route (less pedestrians using the shortest path) would not help reduce the PCR by much under the relatively small population.

From Fig.6b, the expected CCR first decreases slightly and then decreases significantly as princreases, that is,from 21.2 to 20.7%,and to 16.3%when princreases from 0 to around 0.3, and to 1.Typically, evacuation by car tends to cause traffic congestion and hence casualties due to the delay caused by the congestion; this is true even under the relatively small population(that is,a small number of cars)(Mostafizi et al.2019; Wang and Jia 2021a).The traffic congestion would become more severe on some links with the high concentration of cars when using the shortest path,especially on the links with significant traffic capacity reduction due to the bridge failure or road blockage.In this context, evacuation by car following the route would alleviate the traffic congestion on some links, and the alleviation would become significant as more cars drive on the evacuation route (for example, when pris larger than 0.3).Then the overall evacuation by car would be improved a lot and the expected CCR would reduce significantly.

Due to the small variation of the PCR but the large variation of the CCR with pr, the variation of the TCR(Fig.6c) shows a similar trend to the CCR.The route effectiveness in reducing the evacuation risk (that is, the TCR) benefits from the route effectiveness for both evacuations on foot and by car, and further stems from the effectiveness in the congestion alleviation when using the relatively evenly distributed evacuation routes.

Nighttime scenario(C4)From Fig.6a,the PCRs under any given value of pras well as under the nighttime and daytime scenarios are very close, and hence the variations of the PCR with prunder the two scenarios are similar.Evacuation on foot is less likely to cause pedestrian congestion(Mostafizi et al.2019),especially under a relativelysmall number of pedestrians(Wang and Jia 2021a).In this case, the difference in the pedestrian speed or the population distribution across sub-areas between the daytime and nighttime scenarios would not change the pedestrian congestion level and the evacuation performance much(either when following the evacuation route or using the shortest path).Therefore, the PCRs are very close under any given value of prand the two scenarios.

Fig.6 Variation of expected casualty rate H(pr) as a function of proportion of evacuees that use the evacuation route pr for (a) pedestrian casualty rate (PCR), (b) car casualty rate (CCR), and (c) total casualty rate (TCR) under cases C1 and C4

The variations of the CCR with prunder the daytime and nighttime scenarios (Fig.6b) show a similar trend,although the CCR is lower for the nighttime scenario.Compared to the daytime scenario, the population distribution across sub-areas is overall closer to most of the shelters under the nighttime scenario.This would reduce the evacuation delay because of less traffic congestion (for example, only over a shorter evacuation distance) and hence reduce the CCR.However, as princreases, the difference in the expected CCR between the daytime and nighttime scenarios does not seem to change much.Ultimately, the evacuation route usage has a similar impact on the reduction of the CCR under the two scenarios.

The variations of the TCR with prunder the daytime and nighttime scenarios (Fig.6c) show a similar trend due to the similar variations of both the PCR and CCR with prunder the two scenarios.This indicates that following the evacuation route under different population distributions across sub-areas has similar effectiveness in reducing the evacuation risk (that is, the TCR) in this example.

In summary, as presented in Table 5, following the evacuation route can effectively reduce the evacuation risk under a relatively small population, especially for evacuation by car when a relatively large proportion of cars drive on the route in this example.This is true under different population distributions across sub-areas (that is, the daytime and nighttime scenarios).For tsunami evacuation planning,designation of evacuation routes has been applied in many areas of high tsunami risk such as the United States, Japan, and so on (L?vholt et al.2019).However,even many residents may not know their official evacuation routes (Dengler 2005).Therefore, encouraging the public to learn about and follow the evacuation route is highly recommended in evacuation preparedness education.The education can be implemented by reading the evacuationmap(City of Rockaway Beach,Oregon 2019),recognizing the evacuation route (Adiyoso and Kanegae 2012), and so forth.

Table 5 Major results under different simulation cases: Summary 1

5.1.2 Variation of Evacuation Risk with Route Usage and Evacuation Mode

Figure 7 shows the variations of the expected casualty rates H(pr, pc) as a function of [pr, pc] under C1 and C4.Note that the view angle for Fig.7a is different from those of Fig.7c, e to facilitate better illustration.To make detailed observations, the variations of the expected casualty rates with prunder several selected values of pc(that is, pc= 0,0.25, 0.5, 0.75, and 1) are also shown.

Daytime scenario (C1).From Fig.7a—e, the following observations can be made:(1)the PCR decreases as pcand princrease, and more reductions are observed as pcvaries;(2) the CCR increases as pcincreases and decreases as princreases, and much larger changes are found as pcvaries;(3)compared to the PCR, overall the variation of the CCR as a function of [pr, pc] is larger, especially the variation with pc; (4) due to the above variations of the PCR and CCR, the variation of the TCR as a function of [pr, pc]shows a similar trend to that of the CCR.

Under any selected value of pc(except pc= 1, which corresponds to no evacuation on foot),overall the variation of the PCR with pr(Fig.7b) is consistent with that under the uncertain pc(Fig.6a).On the other hand, as princreases,there is more reduction in PCR under the smaller value of pc.For example, the reductions in PCR with the increase of pr(from 0 to 1) are around 2.0% and 0.3%when pc= 0 and pc= 0.5, respectively.This means the evacuation route effectiveness in the reduction of the PCR is higher when fewer people evacuate by car(that is,more people evacuate on foot).As discussed for the PCR under C1 in Sect.5.1.1, using the shortest path for pedestrians tends to cause a higher PCR than following the route.This difference in the PCR would become larger under a larger number of pedestrians due to more severe pedestrian congestion on some links using the shortest path.In this case,the reduction of the PCR by following the route would become larger when more people evacuate on foot.

Overall, the variation of the CCR with prunder any selected non-zero value of pc(Fig.7d) is consistent with that under the uncertain pc(Fig.6b).However, the reduction of the CCR is larger under the larger value of pc,which indicates that the evacuation route usage has larger effectiveness in reducing the CCR when more people evacuate by car.This observation is similar to that for the PCR and can be explained similarly.

The variation of the TCR with prunder any selected value of pc(Fig.7f)is similar to that under the uncertain pc(Fig.6c).Furthermore, the route effectiveness in the reduction of the TCR is larger when more people evacuate by car while overall the TCR is also larger.More importantly, the reduction of the TCR by following the route is much smaller compared to the increase of the TCR as more people evacuate by car.For example, the largest reduction of the TCR is only 11.4% (that is, from 52.9% to 41.5%when princreases from 0 to 1)at pc=1 while the smallest increase of the TCR is 34.4%(that is,from 7.1%to 41.5%)at pr= 1.To reduce the evacuation risk effectively, evacuation by car should be discouraged besides encouraging to follow the evacuation route in this example.

Nighttime scenario (C4).As shown in Fig.7a, c, e, under the nighttime scenario, the variations of the expected casualty rates H(pr,pc)as a function of[pr,pc]demonstrate similarities to those under the daytime scenario, although the casualty rate under the same value of[pr,pc]is overall lower.

Fig.7 Variation of expected casualty rate H(pr, pc) as a function of proportion of evacuees that use the evacuation route and evacuate by car[pr,pc]for(a)pedestrian casualty rate(PCR)under cases C1 and C4, (c) car casualty rate (CCR) under C1 and C4, and (e) total

casualty rate (TCR) under C1 and C4, and variation of expected casualty rate H(pr)as a function of prunder several selected values of pcfor (b) PCR under C1 and C4, (d) CCR under C1 and C4, and(f) TCR under C1 and C4

As for the variations of the expected casualty rates with prunder selected values of pc(Fig.7b,d,f), similar trends can also be found between the nighttime and daytime scenarios, although the PCRs under the two scenarios are relatively close while the CCR and TCR under the nighttime scenario are lower.For example,the reductions of the PCR (Fig.7b) are, respectively, 2.0% and 1.3% at pc= 0 under the daytime scenario and nighttime scenario while the corresponding reductions of the CCR (Fig.7d) are 12.0% and 13.0% at pc= 1.This indicates that following the evacuation route has close effectiveness in the evacuation risk reduction under the daytime and nighttime scenarios.

Based on the above results under the daytime and nighttime scenarios (summarized in Table 5), the evacuation route usage is effective in reducing the evacuation risk.When more people evacuate by car, the reduction in the evacuation risk is higher as more people follow the evacuation routes.These hold under different population distributions across sub-areas (that is, the daytime and nighttime scenarios).However, the TCR increases significantly with car use and the increase is much larger than the risk reduction resulting from following the evacuation route.Existing research has shown that evacuation by car is more likely to cause a larger evacuation risk (Mostafizi et al.2019; Wang and Jia 2021a).Also, significant casualties due to the evacuation delay caused by a large amount of car use have occurred in historical tsunami evacuations(Murakami et al.2014).Therefore, for more effective tsunami evacuation risk reduction, following the evacuation route is encouraged and evacuation on foot instead of by car is recommended in this example.However,typically people would prefer to evacuate by car (Murakami et al.2014).Therefore,to address such issues,strategic planning can be carried out, for example, designing a tsunami evacuation drill to discourage people from evacuating by car through experiencing the evacuation delay caused by traffic congestion (Wang 2021).

5.2 Evaluation of Route Effectiveness under Different Population Sizes

This section investigates the route effectiveness under different population sizes.For each population size, the variations of the evacuation risk (that is, expected PCR,CCR,and TCR)with route usage and with both route usage and evacuation mode under the daytime and nighttime scenarios (C1—C6) are investigated.

5.2.1 Variation of Evacuation Risk with Route Usage

The variations of the expected casualty rates with the proportion of the evacuees that use the evacuation route(that is, pr) for PCR, CCR, and TCR under C1—C6 are shown in Fig.8a—c, respectively.

Fig.8 Variation of expected casualty rate H(pr) as a function of proportion of evacuees that use the evacuation route pr for (a) pedestrian casualty rate (PCR), (b) car casualty rate (CCR), and (c) total casualty rate (TCR) under cases C1—C6

Daytime scenario(C1-C3)The variations of the PCR with prunder ne= 10,000 and ne= 15,000 show some similarities to that under ne= 5,000 (Fig.8a), for example, the PCR would reduce as more pedestrians evacuate through following the evacuation route.However, the reduction ofthe PCR when princreases from 0 to 1 is larger under the relatively large population, that is, 0.3%, 2.8%, and 2.9%under ne= 5,000, 10,000, and 15,000, respectively.Unlike under the relatively small population(that is,ne=5,000 in Sect.5.1.1), the pedestrian congestion level on some links would be much higher when using the shortest path than following the route under the relatively large population(that is, ne= 10,000, ne= 15,000).In this case, the route usage would be more effective in the alleviation of pedestrian congestion and hence in the reduction of casualties.Note that the route effectiveness does not increase much (that is, from 2.8% to 2.9%) as the population increases from 10,000 to 15,000.This is because under the relatively large population, following the evacuation route for pedestrians tends to induce pedestrian congestion as severe as using the shortest path.In this context,following the route would not reduce the evacuation delay and casualty rate for pedestrians much more than using the shortest path.

By comparing the variations of the CCR with prunder different populations shown in Fig.8b, some similarities can be found such as the CCR decreases as princreases.However, overall the reduction of the CCR when princreases from 0 to 1 becomes smaller (though not by much) under the relatively large population, that is, 4.9%,4.5%, and 3.9% under ne= 5,000, ne= 10,000, and ne=15,000, respectively.This is because under a larger population, following the evacuation route for cars would also cause more severe traffic congestion as using the shortest path, which makes the alleviation of the traffic congestion by following the route smaller.

In comparison to the case for ne= 5,000, the variations of the TCR with prunder ne= 10,000, and ne= 15,000(Fig.8c) demonstrate some similar trends, for example,overall the TCR decreases with the increase of pr.As princreases from 0 to 1,however,overall the reduction of the TCR is large under the relatively large population, that is,5.2%,7.3%,and 6.8%under ne=5,000,ne=10,000,and ne= 15,000, respectively.This indicates the route effectiveness in the reduction of the TCR is overall higher under the relatively large population.In addition,the reduction in the TCR is smaller under the relatively large value of pr(for example, when pris larger than 0.9) and relatively large population.This indicates that the route effectiveness in the reduction of the TCR would decrease when most of the evacuees follow the route, especially under the relatively large population.

Nighttime scenario (C4-C6) From Fig.8a—c, under the same population, overall the variations of the PCR, CCR,and TCR with prunder the nighttime scenario are not very different from that under the daytime scenario, although the evacuation risk is overall lower under the nighttime scenario.In this context, the impact of the population size on the variations of the PCR, CCR, and TCR with prdemonstrates similar trends between the daytime and nighttime scenarios,for example,the route effectiveness in the reduction of the TCR is overall higher under the relatively large population (the reductions of the TCR under the nighttime scenario are 6.0%, 8.1%, and 7.3% under ne= 5,000, ne= 10,000, and ne= 15,000, respectively).

To summarize based on the above results (Tables 5 and 6), overall the evacuation route usage can also effectively reduce the evacuation risk under the relatively large population with different distributions across sub-areas, either for the pedestrian, car, or total in this example.For better evacuation preparedness education and hence for the improvement of evacuation, it is recommended to encourage the public to follow the evacuation route when the tsunami strikes.

Table 6 Major results under different simulation cases: Summary 2

5.2.2 Variation of Evacuation Risk with Route Usage and Evacuation Mode

Figure 9 shows the variations of the expected casualty rates H(pr, pc) as a function of [pr, pc] under C1—C6.Note thatthe view angles for Fig.9a, b are different from those of Fig.9c—f to facilitate better illustration.To make more detailed observations, the variations of the expected casualty rates with prunder several selected values of pc(that is, pc= 0, 0.25, 0.5, 0.75, and 1) are shown in Fig.10.

Fig.9 Variation of expected casualty rate H(pr, pc) as a function of proportion of evacuees that use the evacuation route and evacuate by car[pr,pc]for(a)pedestrian casualty rate(PCR)under cases C1—C3, (b) PCR under C4—C6, (c) car casualty rate (CCR) under C1—C3,(d) CCR under C4—C6, (e) total casualty rate (TCR) under C1—C3, and (f) TCR under C4—C6

Daytime scenario(C1—C3)As shown in Fig.9a,c,e,under the daytime scenario, the variations of the expected casualty rate as a function of[pr,pc]under ne=10,000 and ne=15,000 demonstrate similar trends to that under ne= 5000 for PCR,CCR,and TCR.However,under any given value of pr, as pcincreases, overall the reduction in the PCR or the increase in the CCR is larger while the increase in the TCR is smaller under the larger population.

Under the relatively large population(that is,ne=10,000 or ne= 15,000), the variations of the PCR with prunder different given values of pcshow similarities to those under ne=5000(Fig.10a)and can be explained similarly as discussed in Sect.5.1.2.Around medium values of pc(that is,not too small or too large values of pc),however,the variation of the PCR with prunder the relatively large population is different from that under ne=5000.For example,when pc=0.5,the reduction of the PCR is 0.3%,3.2%,and 4.2%under ne= 5000, 10,000, and 15,000, respectively.In this case,overall the route effectiveness in the reduction of the PCR is larger under the larger population.

Fig.10 Variation of expected casualty rate H(pr) as a function of proportion of evacuees that use the evacuation route pr under selected values of proportion of the evacuees that use the car pc for(a) pedestrian casualty rate (PCR) under cases C1—C3, (b) PCR under C4—C6, (c) car casualty rate (CCR) under C1—C3, (d) CCR under C4—C6,(e)total casualty rate(TCR)under C1—C3,and(f)TCR under C4—C6

As for the variations of the CCR with prunder different given values of pc, overall they show similar trends under different populations (Fig.10c) and the trends under the relatively large population (that is, ne= 10,000 or ne=15,000)can be explained similarly as in the case under ne=5000.However, under the relatively small value of pc(except pc=0),overall the amount of reduction in the CCR is larger under the relatively large population.For example,when pc= 0.25, the reduction of the CCR is 1.0%, 2.1%,and 2.5% under ne= 5000, ne= 10,000, and ne= 15,000,respectively.This means that the route effectiveness in the reduction of the CCR is larger under the larger population.

Under different populations,similarities can be observed in the variations of the TCR with prunder different given values of pc(Fig.10e).In particular, overall the route effectiveness in the reduction of the TCR is larger when more people evacuate by car(large value of pc)while overall the TCR is larger, which means the evacuation risk would increase when more people evacuate by car (fewer people evacuate on foot).As pcincreases,however,the reduction of the TCR (when princreases from 0 to 1) keeps increasing from 1.6%to 11.3%under ne=5,000 while fluctuates within[5.8%,8.4%],and[7.3%,9.2%]under ne=10,000,and ne=15,000,respectively.This means that the route effectiveness does not keep increasing with the car use and also does not change much under the relatively large population.

Nighttime scenario (C4—C6) Under the same population,overall the variations of the expected casualty rate H(pr,pc)as a function of [pr, pc] under the nighttime scenario(Fig.9b, d, f) demonstrate similarities to that under the daytime scenario (Fig.9a, c, e), although the casualty rate is overall lower under the nighttime scenario.In this case,the impact of the population size on the variations of the PCR, CCR, and TCR with [pr, pc] demonstrates similar trends between the daytime and nighttime scenarios.

Under the same population, overall the variations of the expected casualty rate with prunder any selected value of pcshow similar trends under the nighttime and daytime scenarios(Fig.10).This means that the route effectiveness in the reduction of the evacuation risk is close under the daytime and nighttime scenarios.Then, the impact of the population size on the variations of the PCR, CCR, and TCR with prunder the given value of pcshows similar trends between the daytime and nighttime scenarios.

Based on the above results under the daytime and nighttime scenarios (summarized in Table 5 and Table 6),overall the evacuation route usage is effective in reducing the evacuation risk (either for the pedestrian, car, or total)under different populations with different distributions across sub-areas in this example.In the evacuation preparedness education, it is recommended to encourage evacuees to follow the evacuation route in tsunami evacuation.In addition, evacuation by car tends to cause a larger evacuation risk.In this case, evacuation on foot is encouraged instead of by car.

6 Conclusion

This study proposed a simulation-based and risk-informed approach for the quantitative evaluation of evacuation route effectiveness in reducing evacuation risk for tsunami preparedness education.An agent-based model was used to simulate the tsunami evacuation.An augmented samplebased approach was used to efficiently estimate the evacuation risk under different proportions of the evacuees that use the evacuation route, as well as under different values of critical risk factors (that is, conditional evacuation risk assessment).Based on the results of the conditional evacuation risk, the route effectiveness in the evacuation risk reduction was evaluated.The results from the quantitative evaluation can provide important insights into the effective use of evacuation routes for an effective evacuation.Tsunami preparedness education programs incorporating the insights can be provided to make the public behave properly in evacuation drills and better prepare for future tsunami evacuation.

As an illustrative example, the evacuation route effectiveness in Seaside, Oregon was investigated using the proposed simulation-based and risk-informed approach.The route effectiveness was evaluated by examining how the proportion of the evacuees that use the evacuation route impacts the reduction in the evacuation risk.The route effectiveness was also evaluated under different proportions of the evacuees that evacuate by car.All the above evaluations of the route effectiveness were conducted under different population sizes and distributions.The results showed that overall the route usage is effective in reducing the evacuation risk.However,the evacuation risk would decrease less or even increase when most of the evacuees follow the route, especially under the relatively large population.Also, route usage is more effective in reducing the evacuation risk by car(that is,CCR)than the pedestrian evacuation risk (that is, PCR).However, it is important to note that as more people evacuate by car the overall evacuation risk (that is, TCR) also increases.The above results can be used to guide the improvement of the tsunami evacuation through preparedness education, for example,encouraging people to evacuate on foot instead of by car and follow the evacuation route.The car use could be discouraged by conducting the evacuation drill and following the evacuation route is encouraged by becoming familiar with the route (for example, reading the evacuation map).

Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing,adaptation,distribution and reproduction in any medium or format,as long as you give appropriate credit to the original author(s) and the source,provide a link to the Creative Commons licence,and indicate if changes were made.The images or other third party material in this article are included in the article’s Creative Commons licence,unless indicated otherwise in a credit line to the material.If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyrightholder.To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.