Earthquake damage assessment for deterministic scenarios in Iquique,Chile

Risk evaluation and loss analysis is key in foreseeing the impact of disasters caused by natural hazards and may contribute effectively in improving resilience in a community through the pre-evaluation of preparedness and mitigation actions. The pilot study presented herein is for the Chilean city of Iquique, which is located at the core of a seismic gap that extends from south Perú to north Chile, and has strategic geopolitical and economic importance for the country. The region was hit April 1, 2014, by an \(M_\mathrm{w}\) 8.2 earthquake that caused only moderate damage, but seismological evidence suggests that there is still a potential for a much larger event in the region. Therefore, a careful damage assessment study is fundamental to anticipate the possible physical, social, and economic consequences that Iquique may face in the future. In this work, the HAZUS-MH platform was adapted and used to simulate a set of ten plausible physics-based future seismic scenarios with magnitudes ranging from \(M_\mathrm{w}\) 8.40 to \(M_\mathrm{w}\) 8.98, which were proposed based on an analysis of interplate locking and the residual slip potential remaining after the April 1, 2014, earthquake. Successful application of this damage assessment methodology relies on the construction of a comprehensive exposure model that takes into account regional features and a good characterization of the physical vulnerabilities. For Iquique, a large body of public and local data was used to develop a detailed inventory of physical and social assets including an aggregated building count, demographics, and essential facilities. To characterize the response of the built environment to seismic demand, appropriate HAZUS fragility curves were applied, and outcomes were validated against the damage observed in the 2014 earthquake. After satisfactory testing, a deterministic earthquake damage assessment study was carried out for the collection of predictive scenarios aimed to estimate their expected impacts. This analysis provides data for future evaluations of different physical and social mitigation measures for the city.     

A deterministic approach for preparation of seismic hazard maps in North East India

A method of seismic zonation based on deterministic modeling of rupture plane is presented in this work. This method is based on the modeling of finite rupture plane along identified lineaments in the region using the semi-empirical technique, of Midorikawa [(1993) Tectonophysics 218:287–295]. The modeling procedure follows ω² scaling law, directivity effects, and other strong motion parameters. The technique of zonation is applied for technoeconomically important NE part of Brahmaputra valley that falls in the seismic gap region of Himalaya. Zonation map prepared for Brahmaputra valley for earthquakes of magnitude M > 6.0 show that approximately 90,000 km² area fall in the highly hazardous zone IV, which covers region that can have peak ground accelerations of order more than 250 cm/s². The zone IV covers the Tezu, Tinsukia, Dibrugarh, Ziro, North Lakhimpur, Itanagar, Sibsagar, Jorhat, Golaghat, Wokha, Senapati, Imphal, and Kohima regions. The Pasighat, Daring, Basar, and Seppa region belong to zone III with peak ground accelerations of the order 200–250 cm/s². The seismic zonation map obtained from deterministic modeling of the rupture is consistent with the historical seismicity map and it has been found that the epicenter of many moderate and major earthquakes fall in the identified zones.     

Probabilistic tsunami hazard assessment for Messina Strait Area (Sicily, Italy)

The general modular Bayesian procedure is applied to provide a probabilistic tsunami hazard assessment (PTHA) for the Messina Strait Area (MSA), Italy. This is the first study in an Italian area where the potential tsunamigenic events caused by both submarine seismic sources (SSSs) and submarine mass failures (SMFs) are examined in a probabilistic assessment. The SSSs are localized on active faults in MSA as indicated by the instrumental data of the catalogue of the Italian seismicity; the SMFs are spatially identified using their propensity to failure in the Ionian and Tyrrhenian Seas on the basis of mean slope and mean depth, and using marine geology background knowledge. In both cases the associated probability of occurrence is provided. The run-ups were calculated at key sites that are main cities and/or important sites along the Eastern Sicily and the Southern Calabria coasts where tsunami events were recorded in the past. The posterior probability distribution combines the prior probability and the likelihood calculated in the MSA. The prior probability is based on the physical model of the tsunami process, and the likelihood is based on the historical data collected by the historical catalogues, background knowledge, and marine geological information. The posterior SSSs and SMFs tsunami probabilities are comparable and are combined to produce a final probability for a full PTHA in MSA.     

A probabilistic tsunami hazard assessment for the Makran subduction zone at the northwestern Indian Ocean

A probabilistic tsunami hazard assessment is performed for the Makran subduction zone (MSZ) at the northwestern Indian Ocean employing a combination of probability evaluation of offshore earthquake occurrence and numerical modeling of resulting tsunamis. In our method, we extend the Kijko and Sellevoll’s (1992) probabilistic analysis from earthquakes to tsunamis. The results suggest that the southern coasts of Iran and Pakistan, as well as Muscat, Oman are the most vulnerable areas among those studied. The probability of having tsunami waves exceeding 5 m over a 50-year period in these coasts is estimated as 17.5%. For moderate tsunamis, this probability is estimated as high as 45%. We recommend the application of this method as a fresh approach for doing probabilistic hazard assessment for tsunamis. Finally, we emphasize that given the lack of sufficient information on the mechanism of large earthquake generation in the MSZ, and inadequate data on Makran’s paleo and historical earthquakes, this study can be regarded as the first generation of PTHA for this region and more studies should be done in the future.     

Subaqueous landslide-triggered tsunami hazard for Lake Zurich,Switzerland

Subaqueous landslides can induce potentially damaging tsunamis. Tsunamis are not restricted to the marine environment, but have also been documented on lakes in Switzerland and worldwide. For Lake Zurich (central Switzerland), previous work documented multiple, assumedly earthquake-triggered landslides. However, no information about past tsunamis is available for Lake Zurich. In a back-analysis, we model tsunami scenarios as a consequence of the earthquake-triggered landslides in the past. Furthermore, on the basis of a recent map of the earthquake-triggered subaqueous landslide hazard, we present results of a tsunami hazard assessment. The subaqueous landslide progression, wave propagation and inundation are calculated with a combination of open source codes. Although no historic evidence of past tsunamis has been documented for Lake Zurich, a tsunami hazard exists. However, only earthquakes with long return periods are assumed to cause considerable tsunamis. An earthquake with an exceedance probability of 0.5% in 50 years (corresponding to an earthquake with a return period of 9975 years) is assumed to cause tsunamigenic landslides on most lateral slopes of Lake Zurich. A hypothetical tsunami for such an event would create damage especially along the shores of the central basin of Lake Zurich with estimated peak flow depths of up to ~?4.6 m. Our results suggest that for an earthquake with an exceedance probability of 10% in 50 years (i.e., mean return period of 475 years), no considerable tsunami hazard is estimated. Even for a worst-case scenario, the cities of Zurich and Rapperswil, located at the northern and southern ends of the lake, respectively, are assumed to experience very little damage. The presented first-order results of estimated wave heights and inundated zones provide valuable information on tsunami-prone areas that can be used for further investigations and mitigation measures.     

A deterministic seismic hazard analysis for shallow earthquakes in Greece

The maximum expected ground motion in Greece is estimated for shallow earthquakes using a deterministic seismic hazard analysis (DSHA). In order to accomplish this analysis the input data include an homogeneous catalogue of earthquakes for the period 426 BC–2003, a seismogenic source model with representative focal mechanisms and a set of velocity models. Because of the discrete character of the earthquake catalogue and of errors in location of single seismic events, a smoothing algorithm is applied to the catalogue of the main shocks to get a spatially smoothed distribution of magnitude. Based on the selected input parameters synthetic seismograms for an upper frequency content of 1 Hz are computed on a grid of 0.2° × 0.2°. The resultant horizontal components for displacement, velocity, acceleration and DGA (Design Ground Acceleration) are mapped. The maps which depict these results cannot be compared with previously published maps based on probabilistic methodologies as the latter were compiled for a mean return period of 476 years. Therefore, in order to validate our deterministic analysis, the final results are compared with PGA estimated from the maximum observed macroseismic intensity in Greece during the period 426 BC–2003.Since the results are obtained for point sources, with the frequency content scaled with moment magnitude, some sensitivity tests are performed to assess the influence of the finite extent of fault related to large events. Sensitivity tests are also performed to investigate the changes in the peak ground motion quantities when varying the crustal velocity models in some seismogenic areas. The ratios and the relative differences between the results obtained using different models are mapped and their mean value computed. The results highlight the importance in the deterministic approach of using good and reliable velocity models.     

A local solution for deterministic and stochastic transport equations

A local method is developed for finding the hydraulic head at an arbitrary point of a soil deposit with conductivity varying randomly in space. The method is said to be local since it delivers the hydraulic head at an arbitrary point of a soil deposit directly rather than extracting it from a field solution. The local method always converges to the exact solution, is ideal for parallel computation, and is simple to implement. The method is applied to solve locally one-dimensional transport equations with mixed boundary conditions, calculate corresponding effective conductivity, and examine size effect in specimens with hydraulic conductivity varying randomly in space.     

Misattributed tsunami: Chile, Sumatra and the subduction model

Tsunami intensity is poorly correlated with earthquake magnitude. The distribution of aftershocks that immediately followed the 2010 Maule (Chile), the 2004 Sumatra–Andaman and the 2005 Nias (Indonesia) events supports the view that faulting within an accretionary wedge or an outer rise can sometimes disrupt the seafloor more effectively than a megathrust even if the associated seismicity is minor. Monitoring offshore faults would thus seem an effective way to supplement modes of tsunami early warning which hinge on instrumental earthquake detection or wave height and period.     

Shallow landslide hazard assessment using a three-dimensional deterministic model in a mountainous area

Shallow landslides are common in mountainous areas after intense rainfall. Of all landslide hazard assessment methods, deterministic methods provide the best quantitative information on landslide hazard. However, they require a large amount of detailed in situ data, derived from laboratory tests and field measurements, and therefore it is difficult to apply them over large areas. One of the most important input parameters is soil depth. For large areas, it is impossible to obtain soil depth through field measurements. To overcome this difficulty, a statistics-based regression analysis is used to evaluate soil depths. All the terrain attributes that control soil depths are selected as influential factors. By using multi-linear regression, the soil depths at each location can be predicted. Slope stability analysis can then be performed using deterministic methods with the evaluated soil depths. The study area is divided into slope units. For each slope unit, Monte-Carlo simulation and a GIS-based 3D limit equilibrium model are used to locate the critical slip surface and calculate the corresponding safety factor. The effectiveness of the proposed method has been tested by applying it to a mountainous area in Japan.     

A structured approach to enhance flood hazard assessment in mountain streams

An evidence-based flood hazard analysis in mountain streams requires the identification and the quantitative characterisation of multiple possible processes. These processes result from specific triggering mechanisms on the hillslopes (i.e. landslides, debris flows), in-channel morphodynamic processes associated with sudden bed changes and stochastic processes taking place at critical stream configurations (e.g. occlusion of bridges, failure of levees). From a hazard assessment perspective, such possible processes are related to considerable uncertainties underlying the hydrological cause-effect chains. Overcoming these uncertainties still remains a major challenge in hazard and risk assessment and represents a necessary condition for a reliable spatial representation of process intensities and the associated probabilities. As a result of an accurate analysis of the conceptual flaws present in the procedures currently employed for hazard mapping in South Tyrol (Italy) and Carinthia (Austria), we propose a structured approach as a means to enhance the integration of hillslope, morphodynamic and stochastic processes into conventional flood hazard prediction for mountain basins. To this aim, a functional distinction is introduced between prevailing one-dimensional and two-dimensional process propagation domains, i.e., between confined and semi- to unconfined stream segments. The former domains are mostly responsible for the generation of water, sediment and wood fluxes, and the latter are where flooding of inactive channel areas (i.e. alluvial fans, floodplains) can occur. For the 1D process propagation domain, we discuss how to carry out a process routing along the stream system and how to integrate numerical models output with expert judgement in order to derive consistent event scenarios, thus providing a consistent quantification of the input variables needed for the associated 2D domains. Within these latter domains, two main types of spatial sub-domains can be identified based on the predictability of their dynamics, i.e., stochastic and quasi-deterministic. Advantages and limitations offered by this methodology are finally discussed with respect to hazard and risk assessment in mountain basins.     

Probabilistic tsunami hazard assessment along Oman coast from submarine earthquakes in the Makran subduction zone

The Sultanate of Oman is among the Indian Ocean countries that were subjected to at least two confirmed tsunamis during the twentieth and twenty-first centuries: the 1945 tsunami due to an earthquake in the Makran subduction zone in the Sea of Oman (near-regional field tsunami) and the Indian Ocean tsunami in 2004, caused by an earthquake from the Andaman Sumatra subduction zone (far - field tsunami). In this paper, we present a probabilistic tsunami hazard assessment for the entire coast of Oman from tectonic sources generated along the Makran subduction zone. The tsunami hazard is assessed taking into account the contribution of small- and large-event magnitudes. Results of the earthquake recurrence rate studies and the tsunami numerical modeling for different magnitudes were used through a logic-tree to estimate the tsunami hazard probabilities. We derive probability hazard exceedance maps for the Omani coast considering the exposure times of 100, 250, 500, and 1000 years. The hazard maps consist of computing the likelihood that tsunami waves exceed a specific amplitude. We find that the probability that a maximum wave amplitude exceeds 1 m somewhere along the coast of Oman reaches, respectively, 0.7 and 0.85 for 100 and 250 exposure times, and it is up to 1 for 500 and 1000 years of exposure times. These probability values decrease significantly toward the southern coast of Oman where the tsunami impact, from the earthquakes generated at Makran subduction zone, is low.     

A visibility-based assessment of tsunami evacuation signs in Seaside,Oregon

Tsunami risk mitigation programs often include iconic evacuation signage to direct locals and visitors to safety during a tsunami event. This paper examines sign placement in Seaside, Oregon, from a visibility perspective. It leverages existing visibility analysis methodologies characterize the visibility of the community’s evacuation signage and reveals patterns in the viewable landscape. Additionally, we develop a topologically 3D approach to visibility analysis using raw LiDAR datasets. This applied work situates a discussion on existing patterns of visibility, how to improve existing signage placement, 2D and 3D representation of landscape, and the importance of visibility analysis. This work aims to stimulate discussion and development of hazard research that incorporates a visibility perspective.     

A method for tsunami risk assessment: a case study for Kamakura,Japan

This paper presents a methodology for tsunami risk assessment, which was applied to a case study in Kamakura, Japan. This methodology was developed in order to evaluate the effectiveness of a risk-reducing system against such hazards, also aiming to demonstrate that a risk assessment is possible for these episodic events. The tsunami risk assessment follows these general steps: (1) determination of the probability of flooding, (2) calculation of flood scenarios, (3) assessment of the consequences and (4) integration into a risk number or graph. The probability of flooding was approximated based on the data provided by local institutes, and the flood scenarios were modeled in 1D using the Simulating WAves till SHore model. Results showed that a tsunami in Kamakura can result in thousands of casualties. Interventions such as improvements in evacuation systems, which would directly reduce the number of casualties, would have a large influence in risk reduction. Although this method has its limits and constraints, it illustrates the value it can add to existing tsunami risk management in Japan.     

Seismic hazard assessment for Mymensingh,Bangladesh

Mymensingh municipality lies in one of the most earthquake-prone areas of Bangladesh. The town was completely destroyed during the Great Indian Earthquake of 12 June 1897, for which the surface-wave magnitude was 8.1. In this study the 1897 Great Indian Earthquake was used as a scenario event for developing seismic microzonation maps for Mymensingh. For microzonation purposes SPT data from 87 boreholes were collected from different relevant organizations. To verify those data ten boreholes of depth up to 30 m were drilled. Intensity values obtained for different events were calibrated against attenuation laws to check applicability to the study area. Vibration characteristics at diverse points of the study area were estimated by employing the one-dimensional wave-propagation software SHAKE. SHAKE discretizes the soil profile into several layers and uses an iterative technique to represent the non-linear behavior of the soil by adjusting the material properties at each iteration step. The required input information includes depth, shear wave velocity, damping factor, and unit weight of each soil layer. The liquefaction resistance factor and the resulting liquefaction potential were estimated to quantify the severity of liquefaction. Quantification of secondary site effects and the weighting scheme for combining the various seismic hazards were heuristic, based on judgment and expert opinion.     

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 systematic approach for the assessment of flooding hazard and risk associated with a landslide dam

Inundation caused by landslide dams may occur in the upstream and downstream of the dams. A proper flooding hazard assessment is required for reaction planning and decision-making to mitigate possible flooding hazards caused by landslide dams. Both quick and detailed procedures can be used to evaluate inundation hazards, depending on the available time and information. This paper presents a systematic approach for the assessment of inundation hazards and risks caused by landslide dam formation and breaches. The approach includes the evaluation of dam-breach probability, assessment of upstream inundation hazard, assessment of downstream inundation hazard, and the classification of flooding risk. The proposed assessment of upstream inundation estimates the potential region of inundation and predicts the overtopping time. The risk level of downstream flooding is evaluated using a joint consideration of the breach probability of a landslide dam and the level of flooding hazard, which is classified using a flooding hazard index that indicates the risk of potential inundation. This paper proposes both quick and detailed procedures for the assessments of inundation in both the upstream and downstream of a landslide dam. An example of a landslide dam case study in southern Taiwan was used to demonstrate the applicability of the systematic approach.     

Shallow landslides in pyroclastic soils: a distributed modelling approach for hazard assessment

An effective assessment of shallow landslide hazard requires spatially distributed modelling of triggering processes. This is possible by using physically based models that allow us to simulate the transient hydrological and geotechnical processes responsible for slope instability. Some simplifications are needed to address the lack of data and the difficulty of calibration over complex terrain at the catchment's scale. We applied two simple hydrological models, coupled with the infinite slope stability analysis, to the May 1998 landslide event in Sarno, Southern Italy. A quasi-dynamic model (Barling et al., 1994) was used to model the contribution to instability of lateral flow by simulating the time-dependent formation of a groundwater table in response to rainfall. A diffusion model [Water Resour. Res. 36 (2000) 1897] was used to model the role of vertical flux by simulating groundwater pressures that develop in response to heavy rainstorms. The quasi-dynamic model overestimated the slope instability over the whole area (more than 16%) but was able to predict correctly slope instability within zero order basins where landslides occurred and developed into large debris flows. The diffusion model simulated correctly the triggering time of more than 70% of landslides within an unstable area amounting to 7.3% of the study area. These results support the hypothesis that both vertical and lateral fluxes were responsible for landslide triggering during the Sarno event, and confirm the utility of such models as tools for hazard planning and land management.