A novel agent-based model for tsunami evacuation simulation and risk assessment

Tsunami evacuation is an effective way to save lives from the near-field tsunami. Realistic evacuation simulation can provide valuable information for accurate evacuation risk assessment and effective evacuation planning. Agent-based modeling is ideal for tsunami evacuation simulation due to its capability of capturing the emergent phenomena and modeling the individual-level interactions among agents and the agents’ interactions with the environment. However, existing models usually neglect or simplify some important factors and/or mechanisms in tsunami evacuation. For example, uncertainties in seismic damages to the transportation network are not probabilistically considered (e.g., by simply removing the damaged links (roads/bridges) from the network). Typically a relatively small population (i.e., evacuees) is considered (due to computational challenges) while neglecting population mobility. These simplifications may lead to inaccurate estimation of evacuation risk. Usually, only single traffic mode (e.g., on foot or by car) is considered, while pedestrian speed adjustment and multi-modal evacuation (e.g., on foot and by car) are not considered concurrently. Also, pedestrian–vehicle interaction is usually neglected in the multi-modal evacuation. To address the above limitations, this study proposes a novel and more realistic agent-based tsunami evacuation model for tsunami evacuation simulation and risk assessment. Uncertainties in seismic damages to all links in the transportation network as well as uncertainties in other evacuation parameters are explicitly modeled and considered. A novel and more realistic multi-modal evacuation model is proposed that explicitly considers the pedestrian–vehicle interaction, walking speed variability, and speed adjustment for both the pedestrian and car according to traffic density. In addition, several different population sizes are used to model population mobility and its impact on tsunami evacuation risk. The proposed model is applied within a simulation-based framework to assess the tsunami evacuation risk assessment for Seaside, Oregon.

     

A GIS-based analysis of constraints on pedestrian tsunami evacuation routes: Cascais case study (Portugal)

Tsunami hazard in coastal areas susceptible to flooding, although reduced (in terms of probability of occurrence), may pose a high risk. Therefore, in these areas, a detailed evacuation planning of the affected population is required as a risk mitigation measure. The knowledge and enforcement of evacuation routes may reduce the population vulnerability, making it more resilient and reducing risk. This paper presents a GIS approach for modelling evacuation routes based on the optimal path search problem, of the graph theory, which is implemented on ArcCasper tool. The methodology proposed considers the elements involved in the evacuation process, the worst credible tsunami inundation scenario (hazard extent and travel time), the number of people that needs to be evacuated in different time scenarios, the safe areas or destination points of the evacuation routes, the roads network characteristics and finally the time available to evacuate. The knowledge of those elements allows predicting some possible outcomes of the evacuation, such as the arrival time of the evacuees to a shelter and the identification of congestion hot spots resulting from the application of a flocking model which simulates the path to be used by evacuees avoiding obstacles. The municipality of Cascais was used to test the methodology proposed in this study. Cascais is one of the largest urban centres located about 25 km west of Lisbon, Portugal, with a high density of infrastructure along the coastline whereby most of the population and economic activities are exposed to a tsunami. The results, presented in the form of maps, allow identifying the optimal evacuation routes as well as the unfeasible routes. This crucial information could be used to the evacuation optimization regarding the location of meeting points and vertical shelters as well as to improve the accessibility of the areas to be evacuated.     

Method to determine the locations of tsunami vertical evacuation shelters

The 2004 Indian Ocean tsunami and the 2011 Great Tohoku Japan earthquake and tsunami focused a great deal of the world??s attention on the effect of tsunamis on buildings and infrastructure. When a tsunami impacts structures in a coastal community, the structures are often not strong enough to withstand the forces and may collapse. Therefore, to maximize the survival probability, people evacuate to higher ground or move outside the inundation zone. However, this is not always possible because of short warning times for near-field tsunamis. Thus, sheltering-in-place or ??sheltering-near-place?? using vertical evacuation should be considered as an alternative approach to lateral evacuation from a tsunami inundation zone. This paper presents the method and results of a study to develop and demonstrate a methodology that applied genetic optimization to determine optimal tsunami shelter locations with the goal of reducing evacuation time, thereby maximizing the probability of survival for the population in a coastal community. The City of Cannon Beach, Oregon, USA, was used as an illustrative example. Several cases were investigated ranging from a single shelter to multiple shelters with locations of high elevation already in place near the city. The method can provide decision-support for the determination of locations for tsunami vertical evacuation shelters. The optimum location of the shelter(s), which was found to vary depending on the number of shelters considered, can reduce the evacuation time significantly, thereby reducing the number of fatalities and increasing the safety of a community.     

Hazard assessment of underwater landslide-generated tsunamis: a case study in the Padang region, Indonesia

Submarine landslides can generate local tsunamis with high run-ups, posing a hazard to human lives and coastal facilities. Both ancient (giant Storegga slide off Norwegian coast, 8200 B. P.) and recent (Papua New Guinea, 1998) events show high potential danger of tsunamigenic landslides and the importance of mitigation efforts. This contribution presents newly discovered landslides 70 km off Padang (Western Sumatra, Indonesia) based on recent bathymetry measurements. This highly populated city with over 750,000 inhabitants exhibits high tsunami vulnerability due to its very low elevation. We model tsunamis that might have been induced by the detected landslide events. Estimations of run-up heights extrapolated from offshore tsunami amplitudes for Padang and other locations in the northern Mentawai fore-arc basin yield maximum values of about 3 m. We also provide a systematic parametric study of landslide-induced tsunamis, which allows us to distinguish potentially dangerous scenarios for Padang. Inside the fore-arc basin, scenarios involving volumes of 0.5–25 km³ could endanger Padang. Apart from slide volume, the hazard distribution mainly depends on three landslide parameters: distance to Padang, water depth in the generation region, and slide direction.     

Agent-Based Modeling and Analysis of Hurricane Evacuation Procedures for the Florida Keys

The unique geography of the Florida Keys presents both high risk of hurricane landfall and exceptional vulnerability to the effects of a hurricane strike. Inadequate hurricane shelters in the Keys make evacuation the only option for most residents, but the sole access road can become impassable well in advance of a major storm. These extraordinary conditions create challenges for emergency managers who must ensure that appropriate emergency plans are in place and to ensure that an orderly exodus can occur without stranding large numbers of people along an evacuation route with inadequate shelter capacity. This study attempts to answer two questions: (1) What is the minimum clearance time needed to evacuate all residents participating in an evacuation of the Florida Keys in advance of a major hurricane for 92,596 people – a population size calculated based on the 2000 US Census population data, census undercounts, and the number of tourists estimated to be in the area? (2) If a hurricane makes landfall in the Keys while the evacuation is in progress, how many residents will need to be accommodated if the evacuation route becomes impassable? The authors conducted agent-based microsimulations to answer the questions. Simulation results suggest that it takes 20 h and 11 min to 20 h and 14 min to evacuate the 92,596 people. This clearance time is less than the Florida state mandated 24-h clearance time limit. If one assumes that people evacuate in a 48-h period and the traffic flow from the Keys would follow that observed in the evacuation from Hurricane Georges, then a total of 460 people may be stranded if the evacuation route becomes impassable 48 h after an evacuation order is issued. If the evacuation route becomes impassable 40 h after an evacuation order is issued, then 14,000 people may be stranded.     

Changes in population evacuation potential for tsunami hazards in Seward, Alaska, since the 1964 Good Friday earthquake

Pedestrian evacuation modeling for tsunami hazards typically focuses on current land-cover conditions and population distributions. To examine how post-disaster redevelopment may influence the evacuation potential of at-risk populations to future threats, we modeled pedestrian travel times to safety in Seward, Alaska, based on conditions before the 1964 Good Friday earthquake and tsunami disaster and on modern conditions. Anisotropic, path distance modeling is conducted to estimate travel times to safety during the 1964 event and in modern Seward, and results are merged with various population data, including the location and number of residents, employees, public venues, and dependent care facilities. Results suggest that modeled travel time estimates conform well to the fatality patterns of the 1964 event and that evacuation travel times have increased in modern Seward due to the relocation and expansion of port and harbor facilities after the disaster. The majority of individuals threatened by tsunamis today in Seward are employee, customer, and tourist populations, rather than residents in their homes. Modern evacuation travel times to safety for the majority of the region are less than wave arrival times for future tectonic tsunamis but greater than arrival times for landslide-related tsunamis. Evacuation travel times will likely be higher in the winter time, when the presence of snow may constrain evacuations to roads.     

Tsunami inundation in Napier, New Zealand, due to local earthquake sources

Deterministic analysis of local tsunami generated by subduction zone earthquakes demonstrates the potential for extensive inundation and building damage in Napier, New Zealand. We present the first high-resolution assessments of tsunami inundation in Napier based on full simulation from tsunami generation to inundation and demonstrate the potential variability of onshore impacts due to local earthquakes. In the most extreme scenario, rupture of the whole Hikurangi subduction margin, maximum onshore flow depth exceeds 8.0 m within 200 m of the shore and exceeds 5.0 m in the city centre, with high potential for major damage to buildings. Inundation due to single-segment or splay fault rupture is relatively limited despite the magnitudes of M_W 7.8 and greater. There is approximately 30 min available for evacuation of the inundation zone following a local rupture, and inundation could reach a maximum extent of 4 km. The central city is inundated by up to three waves, and Napier Port could be inundated repeatedly for 12 h. These new data on potential flow depth, arrival time and flow kinematics provide valuable information for tsunami education, exposure analysis and evacuation planning.     

Maritime tsunami evacuation guidelines for the Pacific Northwest coast of Oregon

Recent tsunamis affecting the West Coast of the USA have resulted in significant damage to ports and harbors, as well as to recreational and commercial vessels attempting to escape the tsunami. With the completion of tsunami inundation simulations for a distant tsunami originating from the Aleutian Islands and a locally generated tsunami on the Cascadia subduction zone (CSZ), the State of Oregon is now able to provide guidance on the magnitudes and directions of the simulated currents for the Oregon coast and shelf region. Our analyses indicate that first wave arrivals for an Aleutian Island event would take place on the north coast,?~?3 h 40 min after the start of the earthquake,?~?20 min later on the southern Oregon coast. The simulations demonstrated significant along-coast variability in both the tsunamis water levels and currents, caused by localized bathymetric effects (e.g., submarine banks and reefs). A locally generated CSZ event would reach the open coast within 7–13 min; maximum inundation occurs at?~?30–40 min. As the tsunami current velocities increase, the potential for damage in ports and harbors correspondingly increases, while also affecting a vessels ability to maintain control out on the ocean. Scientific consensus suggests that tsunami currents?<?1.54 m/s are unlikely to impact maritime safety in ports and harbors. No such guidance is available for boats operating on the ocean, though studies undertaken in Japan suggest that velocities in the region of 1–2 m/s may be damaging to boats. In addition to the effects of currents, there is the added potential for wave amplification of locally generated wind waves interacting with opposing tsunami currents in the offshore. Our analyses explore potential wave amplification effects for a range of generic sea states, ultimately producing a nomogram of wave amplification for a range of wave and opposing current conditions. These data will be useful for US Coast Guard and Port authorities as they evaluate maritime tsunami evacuation options for the Oregon coast. Finally, we identify three regions of hazard (high, moderate, and low) across the Oregon shelf, which can be used to help guide final designation of tsunami maritime evacuation zones for the coast.     

Tsunami evacuation behavior of coastal residents in Kochi Prefecture during the 2014 Iyonada Earthquake

The 2014 Iyonada Earthquake, which occurred at 02:06 JST on 14 March, measured 6.2 on the Richter scale and originated in the Seto Inland Sea of Japan. To elucidate tsunami evacuation behavior, we examined two coastal communities in Kochi Prefecture, Okitsu and Mangyo, where residents evacuated to high ground in anticipation of a tsunami. In the event of a Nankai megathrust earthquake and tsunami, it is expected that a huge tsunami will be generated and these communities will be severely damaged. Before the Iyonada Earthquake, we had previously collected data about tsunami preparedness and evacuation plans from the residents of these communities, and after the earthquake, we conducted in-depth interviews and questionnaire surveys with the residents regarding the actual evacuation behaviors that they took. This enabled us to compare evacuation plans with evacuation behaviors. Results indicate that many residents responded quickly to the earthquake, either by immediately evacuating to emergency shelters on high ground or by preparing themselves for evacuation. Additionally, the earthquake revealed great differences between the prior evacuation plans and the actual situation of residents’ evacuation, such as specific triggers that significantly led residents to evacuate and the use of vehicles in evacuation.     

Spatial and temporal variations in the incidence of dust events over Iran

Large near-field tsunamis pose a significant threat to the Canadian West Coast due to its proximity to the circum-Pacific belt where a significant tsunami-inducing earthquake event from the Cascadia subduction zone is expected. This study investigated the risks associated with such an event in terms of pedestrian evacuation needs and plans for the Town of Tofino, a small community located on the West Coast of Vancouver Island. The population-at-risk within the hazard zone and its ability to evacuate to safety is evaluated using anisotropic path-distance modelling. Mitigation measures, such as vertical evacuation buildings, are quantitatively evaluated. Site-specific inundation modelling was not performed as part of this study; tsunami hazard and safe zones were computed using a range of run-ups varying between 3 and 25 m. It was established that up to 80% of the population is within the maximum hazard zone considered. This evacuation modelling exercise indicates that a maximum of 13% of the population would have insufficient time to reach safety when using a mobility-impaired ambulatory speed. The use of three vertical evacuation buildings can reduce the risk of losing population in this category by 99%. Although some conservative assumptions were used (vertical datum at higher high water, reductions in safe zones by generalization process and mobility-impaired evacuation speeds), the evacuation potential is likely overestimated due to the coarseness of the topographic data used in the evacuation modelling and from an overestimated first wave arrival time. This is the first Canadian study which used anisotropic evacuation modelling to evaluate the vulnerability of a Canadian community to tsunami inundation.

     

Simulation-based optimization of emergency evacuation strategy in ultra-high-rise buildings

Emergency evacuation in high-rise buildings is a crucial problem. The evacuation strategy of using stairs and evacuation elevators should be optimized. In this paper, simulation-based optimization method is used to optimize the evacuation strategy of using stairs and elevators in high-rise buildings. The stair simulation is based on a cellular automata model, and several typical pedestrians’ walk preferences are considered in this model. In the simulation, evacuation elevators can arrive at the refuge floors, and the scheduling of the elevators is optimized based on the GA algorithm. The simulation-based optimization is designed as a two-level problem: The upper level is a strategy level; the lower level is an operation level. In the study case, the evacuation strategy of a 100-floor ultra-high-rise office building is optimized. We find that if evacuees follow the traditional stair evacuation strategy, the evacuation time is 42.6 min. After optimization, the evacuation time of optimal strategy by using both stairs and elevators is 25.1 min. Compared with the traditional stair evacuation strategy, the efficiency of evacuation is improved by 41.1%. It is also found that the merging behavior in stairwells will decrease the velocity of the pedestrian flow. Stairs are still the main egress, and evacuation elevators are an assistant egress during high-rise building evacuation.     

A method for estimating casualties due to the tsunami inundation flow

This study develops a method for estimating the number of casualties that may occur while people evacuate from an inundation zone when a tsunami has inundated an area. The method is based on a simple model of hydrodynamic forces as they affect the human body. The method uses a Tsunami casualty index (TCI) computed at each grid point of a numerical tsunami model to determine locations and times within the tsunami inundation zone where evacuation during the tsunami inundation is not possible and therefore where casualties are likely to occur. The locations and times can be combined with information about population density to compute the potential number of casualties. This information is useful in developing tsunami evacuation routes that avoid such locations. To illustrate the method, it is applied to the Seattle waterfront in Washington State, USA, that is under the threat of possible tsunami disasters due to Seattle Fault earthquakes. Preliminary results suggest that the tsunami casualties may occur within the Seattle waterfront for 15 min, during the time interval from 3 to 18 min after a large Seattle Fault tsunami is generated when the background tide level is mean high water.     

Tsunami impact to Washington and northern Oregon from segment ruptures on the southern Cascadia subduction zone

This paper explores the size and arrival of tsunamis in Oregon and Washington from the most likely partial ruptures of the Cascadia subduction zone (CSZ) in order to determine (1) how quickly tsunami height declines away from sources, (2) evacuation time before significant inundation, and (3) extent of felt shaking that would trigger evacuation. According to interpretations of offshore turbidite deposits, the most frequent partial ruptures are of the southern CSZ. Combined recurrence of ruptures extending ~490 km from Cape Mendocino, California, to Waldport, Oregon (segment C) and ~320 km from Cape Mendocino to Cape Blanco, Oregon (segment D), is ~530 years. This recurrence is similar to frequency of full-margin ruptures on the CSZ inferred from paleoseismic data and to frequency of the largest distant tsunami sources threatening Washington and Oregon, ~M _w 9.2 earthquakes from the Gulf of Alaska. Simulated segment C and D ruptures produce relatively low-amplitude tsunamis north of source areas, even for extreme (20 m) peak slip on segment C. More than ~70 km north of segments C and D, the first tsunami arrival at the 10-m water depth has an amplitude of <1.9 m. The largest waves are trapped edge waves with amplitude ≤4.2 m that arrive ≥2 h after the earthquake. MM V–VI shaking could trigger evacuation of educated populaces as far north as Newport, Oregon for segment D events and Grays Harbor, Washington for segment C events. The NOAA and local warning systems will be the only warning at greater distances from sources.     

A 3-tier tsunami vulnerability assessment technique for the north-west coast of Peninsular Malaysia

The 2004 tsunami that struck the Sumatra coast gave a warning sign to Malaysia that it is no longer regarded as safe from a future tsunami attack. Since the event, the Malaysian Government has formulated its plan of action by developing an integrated tsunami vulnerability assessment technique to determine the vulnerability levels of each sector along the 520-km-long coastline of the north-west coast of Peninsular Malaysia. The scope of assessment is focused on the vulnerability of the physical characteristics of the coastal area, and the vulnerability of the built environment in the area that includes building structures and infrastructures. The assessment was conducted in three distinct stages which stretched across from a macro-scale assessment to several local-scale and finally a micro-scale assessment. On a macro-scale assessment, Tsunami Impact Classification Maps were constructed based on the results of the tsunami propagation modelling of the various tsunami source scenarios. At this stage, highly impacted areas were selected for an assessment of the local hazards in the form of local flood maps based on the inundation modelling output. Tsunami heights and flood depths obtained from these maps were then used to produce the Tsunami Physical Vulnerability Index (PVI) maps. These maps recognize sectors within the selected areas that are highly vulnerable to a maximum tsunami run-up and flood event. The final stage is the development of the Structural Vulnerability Index (SVI) maps, which may qualitatively and quantitatively capture the physical and economic resources that are in the tsunami inundation zone during the worst-case scenario event. The results of the assessment in the form of GIS-based Tsunami-prone Vulnerability Index (PVI and SVI) maps are able to differentiate between the various levels of vulnerability, based on the tsunami height and inundation, the various levels of impact severity towards existing building structures, property and land use, and also indicate the resources and human settlements within the study area. Most importantly, the maps could help planners to establish a zoning scheme for potential coastline development based on its sensitivity to tsunami. As a result, some recommendations on evacuation routes and tsunami shelters in the potentially affected areas were also proposed to the Government as a tool for relief agencies to plan for safe evacuation.     

Hurricane preparedness and planning in coastal public school districts in the United States

One hundred school districts were surveyed along the Atlantic and Gulf of Mexico coasts from North Carolina to Texas. Nearly all had recent experience with a tropical storm or hurricane and had hurricane plans in place. About half teach hurricane preparedness to students and 85 % train staff members in hurricane preparedness. Sources of information about cyclone threats were the National Hurricane Center (91 %), local television news (74 %), The Weather Channel (67 %), and the internet (67 %). Only 36 % would cancel classes for a hurricane warning but 89 % would cancel classes for a mandatory evacuation. Most districts (75 %) would use schools as storm shelters, and 92 % would use school busses to assist in community evacuations. Districts with a higher percentage of Hispanic population provided hurricane information in Spanish. Larger school districts were less likely to cancel classes in the middle of the day for a storm threat. Districts with higher home values were less likely to use school busses for evacuations, and smaller school districts were less likely to provide schools as storm shelters. There were no other significant associations between hurricane preparedness of the districts and district demographic variables of poverty, percent black, percent Hispanic, population, district size, or median home values.     

Assessing building vulnerability to earthquake and tsunami hazard using remotely sensed data

Quantification of building vulnerability to earthquake and tsunami hazards is a key component for the implementation of structural mitigation strategies fostering the essential shift from post-disaster crisis reaction to preventive measures. Facing accelerating urban sprawl and rapid structural change in modern urban agglomerations in areas of high seismic and tsunami risk, the synergetic use of remote sensing and civil engineering methods offers a great potential to assess building structures up-to-date and area-wide. This paper provides a new methodology contextualizing key components in quantifying building vulnerability with regard to sequenced effects of seismic and tsunami impact. The study was carried out in Cilacap, a coastal City in Central Java, Indonesia. Central is the identification of significant correlations between building characteristics, easily detectable by remote sensing techniques, and detailed in situ measurements stating precise building vulnerability information. As a result, potential vertical evacuation shelters in the study area are detected and a realistic vulnerability assessment of the exposed building stock is given. These findings obtained allow for prioritization of intervention measures such as awareness and preparedness strategies and can be implemented in local disaster management.     

Community variations in population exposure to near-field tsunami hazards as a function of pedestrian travel time to safety

Efforts to characterize population exposure to near-field tsunami threats typically focus on quantifying the number and type of people in tsunami-hazard zones. To develop and prioritize effective risk-reduction strategies, emergency managers also need information on the potential for successful evacuations and how this evacuation potential varies among communities. To improve efforts to properly characterize and differentiate near-field tsunami threats among multiple communities, we assess community variations in population exposure to tsunamis as a function of pedestrian travel time to safety. We focus our efforts on the multiple coastal communities in Grays Harbor and Pacific Counties (State of Washington, USA), where a substantial resident and visitor population is threatened by near-field tsunamis related to a potential Cascadia subduction zone earthquake. Anisotropic, path distance modeling is conducted to estimate travel times to safety, and results are merged with various population data, including residents, employees, public venues, and dependent-care facilities. Results suggest that there is substantial variability among communities in the number of people that may have insufficient time to evacuate. Successful evacuations may be possible in some communities assuming slow walking speeds, are plausible in others if travel speeds are increased, and are unlikely in another set of communities given the large distances and short time horizon. Emergency managers can use these results to prioritize the location and determine the most appropriate type of tsunami risk-reduction strategies, such as education and training in areas where evacuations are plausible and vertical-evacuation structures in areas where they are not.     

Recurrence intervals of major paleotsunamis as calibrated by historic tsunami deposits in three localities: Port Alberni, Cannon Beach, and Crescent City, along the Cascadia margin, Canada and USA

We review geologic records of both historic and prehistoric tsunami inundations at three widely separated localities that experienced significant damage from the 1964 Alaskan tsunami along the Cascadia margin. The three localities are Port Alberni, Cannon Beach, and Crescent City, representing, respectively, the north, central, and south portions of the study area (1,000 km in length). The geologic records include anomalous sand sheets from marine surges that are hosted in supratidal peaty mud deposits. Paleotsunami sand sheets that exceed the thickness, continuity and/or extent of the 1964 historic tsunami are counted as major paleotsunami inundations. Major paleotsunamis (6–7 in number) at each locality during the last 3,000 years demonstrate mean recurrence intervals of 450–540 years, and within-cluster intervals (three events each) of 270–460 years. It has been 313 years since the last major paleotsunami from a great Cascadia earthquake in AD 1700. We compare the dated sequences of major paleotsunami inundations to the nearest regional records of coastal coseismic subsidence in Willapa Bay in the central margin, Waatch/Neah Bay in the northern margin, and Coquille in the southern margin. Similar numbers of events from both types of records suggest that the major paleotsunamis are locally derived (near-field) from ruptures of the Cascadia margin megathrust fault zone, rather than from transoceanic tsunamis (far-field) originating at other subduction zones around the Pacific Rim. Given the catastrophic hazard of the near-field Cascadia margin tsunamis, we propose a basic rule for reminding the general public of the need for self-initiated evacuation following a great earthquake at the Cascadia margin.     

Tsunamis from the Arica-Tocopilla source region and their effects on ports of Central Chile

The last great earthquake in northern Chile took place in 1877, and the ensuing tsunami affected not only that region but also Central Chile. For example, the Bay of Concepción, which is located 1,500 km south of the tsunami source, experienced an inundation height of around 3 m. Ports are important in the Chilean economy, due to the fact that a large percentage of Chilean exports (excluding copper) use ports located in Central Chile. With this in mind, the authors investigated the potential effect of an 1877-like tsunami on the main ports of Central Chile. To do this, the dispersive wave model Non-hydrostatic Evolution of Ocean WAVEs was used. In addition, the first tsunami forecast model for Talcahuano, inside the Bay of Concepción, was developed by means of numerical simulation of several events of different moment magnitudes. The results showed that most of the important ports (Valparaiso, San Antonio, San Vicente and Coronel) had inundation heights on the order of just 1 m, while inundation levels in Talcahuano reached up to 3.5 m. The forecast model for Talcahuano uses only earthquake magnitude, focal depth and tide level to determine tsunami inundation heights. In addition, the tsunami arrival time was computed to be 3 h, and the maximum tsunami amplitude takes place at 4 h and 45 min after the earthquake.     

Advancing tsunami risk assessment by improving spatio-temporal population exposure and evacuation modeling

Tsunamis are among the most destructive and lethal of coastal hazards. These are time-specific events, and despite directly affecting a narrow strip of coastline, a single occurrence can have devastating effects and cause massive loss of life, especially in urbanized coastal areas. In this work, in order to consider the time dependence of population exposure to tsunami threat, the variation of spatio-temporal population distribution in the daily cycle is mapped and analyzed in the Lisbon Metropolitan Area. High-resolution daytime and nighttime population distribution maps are developed using ‘intelligent dasymetric mapping,’ that is, applying areal interpolation to combine best-available census data and statistics with land use and land cover data. Workplace information and mobility statistics are considered for mapping daytime distribution. In combination with a tsunami hazard map, information on infrastructure, land use and terrain slope, the modeled population distribution is used to assess people’s evacuation speed, applying a geospatial evacuation modeling approach to the city of Lisbon. The detailed dynamic population exposure assessment allows producing both daytime and nighttime evacuation time maps, which provide valuable input for evacuation planning and management. Results show that a significant amount of population is at risk, and its numbers increase dramatically from nighttime to daytime, especially in the zones of high tsunami flooding susceptibility. Also, full evacuation can be problematic in the daytime period, even if initiated immediately after a major tsunami-triggering earthquake. The presented approach greatly improves tsunami risk assessment and can benefit all phases of the disaster management process.