Probabilistic assessment of annual repair rates in pipelines and of direct economic losses in water and sewage networks: application to Manizales,Colombia

A methodology for the development of fully probabilistic seismic risk assessments on water and sewage networks is presented in this paper together with a case study for the system of Manizales, Colombia. These kinds of assessments require the development of probabilistic seismic hazard analysis, the consideration of local site effects, when relevant, the assembly of databases to identify and characterize the exposed elements and the development and assignment of vulnerability models for each type of component. For the case of Manizales, a high-resolution exposure database has been developed (element by element, segment by segment) based on the information and data provided by the owner and operator of the network, Aguas de Manizales. Losses due to earthquakes are obtained after convoluting the hazard and vulnerability inputs in a fully probabilistic manner, using the state-of-the-art methodologies incorporated in the CAPRA risk assessment module. Several risk metrics such as the loss exceedance curve, the loss exceedance probabilities for different time frames and the average annual loss are obtained for the system as a whole as well as disaggregated by component. In addition, repair rates for the pipelines were also calculated. The risk results obtained in this study have been useful for the company in designing and implementing expansion and maintenance plans that explicitly account for seismic risk mitigation issues, as well as to explore and negotiate financial protection alternatives by means of risk transfer and retention schemes, thus becoming a valuable input in the continuous development of good disaster risk management practices in this city.

     

Probabilistic earthquake risk assessment using CAPRA: application to the city of Barcelona, Spain

The risk evaluation model CAPRA (Comprehensive Approach to Probabilistic Risk Assessment) is a techno-scientific methodology and information platform, composed of tools for evaluating and communicating risk at various territorial levels. The model allows evaluating losses on exposed elements using probabilistic metrics, such as the loss exceedance curve, the expected annual loss and the probable maximum loss, useful for multi-hazard risk analyses. In this article, the process of probabilistic seismic risk analysis is described, explaining the main features of the CAPRA modules of hazard, vulnerability and risk estimation applied to the city of Barcelona, Spain. In addition, according to the physical risk results and the information on the socioeconomic indicators of the city, this article presents the holistic evaluation of seismic risk, which is a valuable result to facilitate the integrated risk management by the different stakeholders involved in risk reduction decision making.     

Design and implementation of a voluntary collective earthquake insurance policy to cover low-income homeowners in a developing country

Understanding and evaluating disaster risk due to natural hazard events such as earthquakes creates powerful incentives for countries to develop planning options and tools to reduce potential damages. The use of models for earthquake risk evaluation allows obtaining outputs such as the loss exceedance curve, the expected annual loss and the probable maximum loss, which are probabilistic metrics useful for risk analyses, for designing strategies for risk reduction and mitigation, for emergency response strategies and for risk financing. This article presents, based on probabilistic risk models, the design and implementation of a risk transfer instrument to cover the private buildings of the city of Manizales, Colombia. This voluntary collective instrument provides financial protection to both, the estate tax payers and the low-income homeowners through a cross-subsidy strategy; besides, it promotes not only the insurance culture but also the solidarity of the community. The city administration and the insurance industry are promoting this program using the mechanism of the property tax payment. This collective insurance helps the government to access key resources for low-income householders recovery and improve disaster risk management at local level.     

Vulnerability index and capacity spectrum based methods for urban seismic risk evaluation. A comparison

This article contributes to the development and application of two latest-generation methods of seismic risk analysis in urban areas. The first method, namely vulnerability index method (VIM), considers five non-null damage states, defines the action in terms of macroseismic intensity and the seismic quality of the building by means of a vulnerability index. The estimated damage degree is measured by semi-empirical functions. The second method, namely capacity spectrum based method (CSBM), considers four no damage states, defines the seismic action in terms of response spectra and the building vulnerability by means of its capacity spectrum. In order to apply both methods to Barcelona (Spain) and compare the results, a deterministic and a probabilistic hazard scenario with soil effects are used. The deterministic one corresponds to a historic earthquake, while the probabilistic seismic ground motion has a probability of exceedence of 10% in 50 years. Detailed information on the building design has been obtained along years by collecting, arranging, improving, and completing the database of the dwellings of the city. A Geographic Information System (GIS) has been customized allowing storing, analysing, and displaying this large amount of spatial and tabular data of dwellings. The obtained results are highly consistent with the historical and modern evolution of the populated area and show the validity and strength of both methods. Although Barcelona has a low to moderate seismic hazard, its expected seismic risk is significant because of the high vulnerability of its buildings. Cities such as Barcelona, located in a low to moderate seismic hazard region, are usually not aware of the seismic risk. The detailed risk maps obtained offer a great opportunity to guide the decision making in the field of seismic risk prevention and mitigation in Barcelona, and for emergency planning in the city.     

PSHA of Van province for performance assessment using spectrally matched strong ground motion records

Within the framework of the performance based earthquake engineering, site specific earthquake spectra for Van province has been obtained. It is noteworthy that, in probabilistic seismic hazard assessment, as a first stage data from geological studies and records from the instrumental period were compiled to make a seismic source characterization for the study region. The probabilistic seismic hazard curves were developed based on selected appropriate attenuation relationships, at rock sites, with a probability of exceedance 2, 10 and 50% in 50 yrs period. The obtained results are compared with the spectral responses proposed for seismic evaluation and retrofit of building structure in Turkish Earthquake Code (2007), section 7. The acceleration response spectrums obtained from probabilistic seismic hazard analysis are matched to adjust earthquake accelerograms recorded during the 2011 Van earthquakes by using SeismoMatch v2.0 software. The aim of this procedure is to obtain a set of reasonable earthquake input motions for the seismic evaluation of existing buildings.     

Probabilistic seismic and tsunami hazard analysis for design criteria development of a coastal area in the city of Banda Aceh

Both seismic and tsunami hazards design criteria are essential input to the rehabilitation and long-term development of city of Banda Aceh Post Sumatra 2004 (M _w=9.3) disaster. A case study to develop design criteria for future disaster mitigation of the area is presented. The pilot study consists of probabilistic seismic and tsunami hazard analysis. Results of the probabilistic seismic hazard analysis indicates that peak ground acceleration at baserock for 10 and 2% probability of exceedance in 50 years is 0.3 and 0.55 g, respectively. The analysis also provides spectral values at short (T=0.2 s) and long period (T=1.0 s) motions. Some non-linear time-domain earthquake response analyses for soft, medium, and hard site-class were conducted to recommend design response spectra for each site-class. In addition, tsunami inundation maps generated from probabilistic tsunami hazard analysis were developed through tsunami wave propagation analysis and run-up numerical modeling associated with its probability of tsunamigenic earthquake source potential. Both the seismic and tsunami hazard curve and design criteria are recommended as contribution of this study for design criteria, as part of the disaster mitigation effort in the development process of the city. The methodology developed herein could be applied to other seismic and tsunami disaster potential areas.     

Non-Poisson probabilistic seismic hazard assessment

The development of the new seismic hazard map of metropolitan Tehran is based on probabilistic seismic hazard computation using the non-Poisson recurrence time model. For this model, two maps have been prepared to indicate the earthquake hazard of the region in the form of iso-acceleration contour lines. They display the non-Poisson probabilistic estimates of peak ground accelerations over bedrock for 10 and 63 % probability of exceedance in 50 years. To carry out the non-Poisson seismic hazard analysis, appropriate distributions of interoccurrence times of earthquakes were used for the seismotectonic provinces which the study region is located and then the renewal process was applied. In order to calculate the seismic hazard for different return periods in the probabilistic procedure, the study area encompassed by the 49.5–54.5°E longitudes and 34–37°N latitudes was divided into 0.1° intervals generating 1,350 grid points. PGA values for this region are estimated to be 0.30–0.32 and 0.16–0.17 g for 10 and 63 % probability of exceedance, respectively, in 50 years for bedrock condition.     

Comprehensive seismic hazard assessment of Tripura and Mizoram states

Northeast India is one of the most highly seismically active regions in the world with more than seven earthquakes on an average per year of magnitude 5.0 and above. Reliable seismic hazard assessment could provide the necessary design inputs for earthquake resistant design of structures in this region. In this study, deterministic as well as probabilistic methods have been attempted for seismic hazard assessment of Tripura and Mizoram states at bedrock level condition. An updated earthquake catalogue was collected from various national and international seismological agencies for the period from 1731 to 2011. The homogenization, declustering and data completeness analysis of events have been carried out before hazard evaluation. Seismicity parameters have been estimated using G–R relationship for each source zone. Based on the seismicity, tectonic features and fault rupture mechanism, this region was divided into six major subzones. Region specific correlations were used for magnitude conversion for homogenization of earthquake size. Ground motion equations (Atkinson and Boore 2003; Gupta 2010) were validated with the observed PGA (peak ground acceleration) values before use in the hazard evaluation. In this study, the hazard is estimated using linear sources, identified in and around the study area. Results are presented in the form of PGA using both DSHA (deterministic seismic hazard analysis) and PSHA (probabilistic seismic hazard analysis) with 2 and 10% probability of exceedance in 50 years, and spectral acceleration (T = 0. 2 s, 1.0 s) for both the states (2% probability of exceedance in 50 years). The results are important to provide inputs for planning risk reduction strategies, for developing risk acceptance criteria and financial analysis for possible damages in the study area with a comprehensive analysis and higher resolution hazard mapping.     

Seismic-event-based methodology to obtain earthquake-induced translational landslide regional hazard maps

A seismic-event-based methodology to generate earthquake-induced translational landslide maps using Newmark method is proposed. The steps are: (1) to construct a GIS-based geotechnical database; (2) to identify those areas that are susceptible to the occurrence of translational landslides based on available geological information; (3) to compute a static safety factor; (4) to compute the critical acceleration that defines the threshold acceleration required to cause a displacement; (5) to characterize the seismic hazard as a set of stochastic events, collectively exhaustive and mutually exclusive, that fully describes the hazard spatial distribution and annual frequency of occurrence (in accordance with the earthquake location, depth and magnitude) with the appropriate ground-motion prediction equations; (6) to compute the Newmark displacement; and finally, (7) to carry out a probabilistic translational landslide hazard analysis to estimate an exceedance rate of a given displacement. This methodology is applied to Mexico, and maps for return periods of 150 and 500 years are presented. Results shown in maps are estimations of where translational landslides may occur and should be useful to carry out local studies to elaborate recommendations of site specific hazard reduction plans as well as to calculate insurance rates. In addition, these results are useful to identify civil protection actions, risk management at regional and local level, and land use planning, as well as for promoting more detailed vulnerability and risk studies at different scales.     

Assessment of Global Seismic Loss Based on Macroeconomic Indicators

Most earthquake loss studies use a probabilistic approach in which predicted damages in various categories of structure and facilities in the region in concern are estimated and added together to obtain a total loss for particular intensity ranges. Such an approach requires a detailed inventory database of the structures and facilities in the region, which is not always readily available in many regions of the world. We have used an alternative means of estimating earthquake losses based on several macroeconomic indices such as the gross domestic product and population. Using published earthquake loss data for 1980–1995, the relations between GDP and earthquake loss have been formulated empirically for several intensity ranges. The world's land surface was divided into unit cells 0.5° lat. × 0.5° long. in size. The GDP of each cell was apportioned based on its population and GDP, and the population of the region to which it belongs. The predicted seismic loss of the cell was then estimated from the seismic hazard probability function, its GDP, and the empirical relation between GDP and seismic loss. A global seismic loss map is then compiled both for the intensity at 10% probability of exceedance and the probable maximum intensity. Employing readily available socioeconomic data as the basis for the vulnerability analysis, the method enables us to obtain seismic loss estimates for regions without the need for a detailed inventory of exposed structures or collateral geological information. Since such statistics are frequently compiled by the world's leading political and financial institutions, the seismic loss estimates can also be upgraded easily for the fast developing areas of the world.     

A Hazard Assessment Method for Potential Earthquake-Induced Landslides – A Case Study in Huaxian County, Shaanxi Province

The hazard assessment of potential earthquake-induced landslides is an important aspect of the study of earthquake-induced landslides. In this study, we assessed the hazard of potential earthquake-induced landslides in Huaxian County with a new hazard assessment method. This method is based on probabilistic seismic hazard analysis and the Newmark cumulative displacement assessment model. The model considers a comprehensive suite of information, including the seismic activities and engineering geological conditions in the study area, and simulates the uncertainty of the intensity parameters of the engineering geological rock groups using the Monte Carlo method. Unlike previous assessment studies on ground motions with a given exceedance probability level, the hazard of earthquake-induced landslides obtained by the method presented in this study allows for the possibility of earthquake-induced landslides in different parts of the study area in the future. The assessment of the hazard of earthquake-induced landslides in this study showed good agreement with the historical distribution of earthquake-induced landslides. This indicates that the assessment properly reflects the macroscopic rules for the development of earthquake-induced landslides in the study area, and can provide a reference framework for the management of the risk of earthquake-induced landslides and land planning.     

Comparative study of three probabilistic methods for seismic hazard analysis: case studies of Sochi and Kamchatka

Natural Hazards - This study examines the effect of the procedures used in three different probabilistic seismic hazard analysis (PSHA) methods for estimating the rates of exceedance of ground...     

Seismic risk in the city of Al Hoceima (north of Morocco) using the vulnerability index method,applied in Risk-UE project

Al Hoceima is one of the most seismic active regions in north of Morocco. It is demonstrated by the large seismic episodes reported in seismic catalogs and research studies. However, seismic risk is relatively high due to vulnerable buildings that are either old or don’t respect seismic standards. Our aim is to present a study about seismic risk and seismic scenarios for the city of Al Hoceima. The seismic vulnerability of the existing residential buildings was evaluated using the vulnerability index method (Risk-UE). It was chosen to be adapted and applied to the Moroccan constructions for its practicality and simple methodology. A visual inspection of 1102 buildings was carried out to assess the vulnerability factors. As for seismic hazard, it was evaluated in terms of macroseismic intensity for two scenarios (a deterministic and probabilistic scenario). The maps of seismic risk are represented by direct damage on buildings, damage to population and economic cost. According to the results, the main vulnerability index of the city is equal to 0.49 and the seismic risk is estimated as Slight (main damage grade equal to 0.9 for the deterministic scenario and 0.7 for the probabilistic scenario). However, Moderate to heavy damage is expected in areas located in the newer extensions, in both the east and west of the city. Important economic losses and damage to the population are expected in these areas as well. The maps elaborated can be a potential guide to the decision making in the field of seismic risk prevention and mitigation strategies in Al Hoceima.     

Local seismic hazard assessment in areas of weak to moderate seismicity—case study from Eastern Germany

Generally the seismic hazard of an area of interest is considered independent of time. However, its seismic risk or vulnerability, respectively, increases with the population and developing state of economy of the area. Therefore, many areas of moderate seismic hazard gain increasing importance with respect to seismic hazard and risk analysis. However, these areas mostly have a weak earthquake database, i.e., they are characterised by relative low seismicity and uncertain information concerning historical earthquakes. In a case study for Eastern Thuringia (Germany), acting as example for similar places in the world, seismic hazard is estimated using the probabilistic approach. Because of the lack of earthquakes occurring in the recent past, mainly historical earthquakes have to be used. But for these the actual earthquake sources or active faults, needed for the analysis, are imprecisely known. Therefore, the earthquake locations are represented by areal sources, a common practice. The definition of these sources is performed carefully, because their geometrical shape and size (apart from the earthquake occurrence model) influence the results significantly. Using analysis tools such as density maps of earthquake epicentres, seismic strain and energy release support this. Oversizing of areal sources leads to underestimation of seismic hazard and should therefore be avoided. Large location errors of historical earthquakes on the other hand are represented by several alternative areal sources with final superimposition of the different results. In a very similar way information known from macroseismic observations interpreted as source rather than as site effects are taken into account in order to achieve a seismic hazard assessment as realistic as possible. In very local cases the meaning of source effects exceeds those of site effects very likely. The influence of attenuation parameter variations on the result of estimated local seismic hazard is relatively low. Generally, the results obtained by the seismic hazard assessment coincide well with macroseismic observations from the thoroughly investigated largest earthquake in the region.     

Seismic hazard analysis of Laoag City,Northern Philippines for liquefaction potential assessment

A seismic hazard analysis was conducted in Laoag City, Northern Philippines to determine the design ground motion for liquefaction potential assessment of the area. Because the hazard analysis was done within the framework of liquefaction potential assessment, only those earthquakes with magnitude–distance combinations that are capable of generating liquefaction were considered in the study. Both probabilistic and deterministic approaches were used in the analysis. From the results of the probabilistic analysis, seismic hazard curves were generated from which the ground motion with a 10% probability of exceedance in 50years was obtained. This was then modified in consideration of the soft soil condition in the study area. Deaggregation was performed to determine the most likely earthquake to generate the said level of ground shaking.     

The comparison of the NDSHA, PSHA seismic hazard maps and real seismicity for the Italian territory

Rigorous and objective testing of seismic hazard assessments against the real seismic activity must become the necessary precondition for any responsible seismic risk estimation. Because seismic hazard maps seek to predict the shaking that would actually occur, the reference hazard maps for the Italian seismic code, obtained by probabilistic seismic hazard assessment (PSHA), and the alternative ground shaking maps based on the neo-deterministic approach (NDSHA), are cross-compared and tested against the real seismicity for the territory of Italy. The comparison between predicted intensities and those reported for past earthquakes shows that models generally provide rather conservative estimates, except for PGA with 10 % probability of being exceeded in 50 years, which underestimates the largest earthquakes. In terms of efficiency in predicting ground shaking, measured accounting for the rate of underestimated events and for the territorial extent of areas characterized by high seismic hazard, the NDSHA maps appear to outscore the PSHA ones.     

The importance of paleoseismology in seismic hazard studies for critical facilities

This paper is intended to provide a perspective on the use of paleoseismological studies in the seismic hazard assessment of critical facilities, such as dams, chemical/petrochemical facilities and nuclear power plants. In particular, the use of data obtained from paleoseismological studies for probabilistic seismic hazard analyses, when the required probabilities of exceedance are very low (e.g. 10^{− 6}–10^{− 7}) is considered. Recent revisions to the IAEA Safety Standards that provide guidance to Member States in their work related to the seismic safety of nuclear power plants are presented to illustrate the importance of this emerging discipline.     

Assessment of seismic hazard and liquefaction potential of Gujarat based on probabilistic approaches

Gujarat is one of the fastest-growing states of India with high industrial activities coming up in major cities of the state. It is indispensable to analyse seismic hazard as the region is considered to be most seismically active in stable continental region of India. The Bhuj earthquake of 2001 has caused extensive damage in terms of causality and economic loss. In the present study, the seismic hazard of Gujarat evaluated using a probabilistic approach with the use of logic tree framework that minimizes the uncertainties in hazard assessment. The peak horizontal acceleration (PHA) and spectral acceleration (Sa) values were evaluated for 10 and 2?% probability of exceedance in 50?years. Two important geotechnical effects of earthquakes, site amplification and liquefaction, are also evaluated, considering site characterization based on site classes. The liquefaction return period for the entire state of Gujarat is evaluated using a performance-based approach. The maps of PHA and PGA values prepared in this study are very useful for seismic hazard mitigation of the region in future.     

Probabilistic seismic hazard analysis for Bangalore

This article presents the results of probabilistic seismic hazard analysis (PSHA) for Bangalore, South India. Analyses have been carried out considering the seismotectonic parameters of the region covering a radius of 350 km keeping Bangalore as the center. Seismic hazard parameter ‘b’ has been evaluated considering the available earthquake data using (1) Gutenberg–Richter (G–R) relationship and (2) Kijko and Sellevoll (1989, 1992) method utilizing extreme and complete catalogs. The ‘b’ parameter was estimated to be 0.62 to 0.98 from G–R relation and 0.87 ± 0.03 from Kijko and Sellevoll method. The results obtained are a little higher than the ‘b’ values published earlier for southern India. Further, probabilistic seismic hazard analysis for Bangalore region has been carried out considering six seismogenic sources. From the analysis, mean annual rate of exceedance and cumulative probability hazard curve for peak ground acceleration (PGA) and spectral acceleration (Sa) have been generated. The quantified hazard values in terms of the rock level peak ground acceleration (PGA) are mapped for 10% probability of exceedance in 50 years on a grid size of 0.5 km × 0.5 km. In addition, Uniform Hazard Response Spectrum (UHRS) at rock level is also developed for the 5% damping corresponding to 10% probability of exceedance in 50 years. The peak ground acceleration (PGA) value of 0.121 g obtained from the present investigation is slightly lower (but comparable) than the PGA values obtained from the deterministic seismic hazard analysis (DSHA) for the same area. However, the PGA value obtained in the current investigation is higher than PGA values reported in the global seismic hazard assessment program (GSHAP) maps of Bhatia et al. (1999) for the shield area.