Overview of Taiwan Earthquake Loss Estimation System

The National Science Council (NSC) of Taiwan started the HAZ-Taiwan project in 1998 to promote researches on seismic hazard analysis, structural damage assessment, and socio-economic loss estimation. The associated application software, “Taiwan Earthquake Loss Estimation System (TELES)”, integrates various inventory data and analysis modules to fulfill three objectives. First, it helps to obtain reliable estimates of seismic hazards and losses soon after occurrence of large earthquakes. Second, it helps to simulate earthquake scenarios and to provide useful estimates for local governments or public services to propose their seismic disaster mitigation plans. Third, it helps to provide catastrophic risk management tools, such as proposing the seismic insurance policy for residential buildings. This paper focuses on the development and application of analysis modules used in early loss estimation system. These modules include assessments of ground motion intensity, soil liquefaction potential, building damage and casualty.     

Analysis of input data and seismic hazard assessment for the low seismicity area ‘Belgium,The Netherlands and NW Germany’

From the earthquake catalogue of the region of interest 2–8°E, 49–52°N, only events with intensity > 3 are considered. The fore/aftershocks and the induced seismicity are removed. The completeness and other properties of the remaining data set are investigated. The seismic data, as well as the geological, tectonic, and other geophysical data of the area are used for the determination of the seismotectonic zones. For each zone, the cumulative intensity-frequency relation (taking into account the completeness of the catalogue), the attenuation depending on the direction, the upper bound of intensity and the average depth are calculated. When the seismic activity is changing within a zone, this zone is divided into subzones with an almost constant activity and a b value equal to the value of the whole zone. When necessary for each subzone, a different attenuation coefficient can be used, to take into account more regional effects. Then the seismic hazard is assessed with the modified McGuire program. Annual probability versus intensity for several sites and hazard intensity and probability maps are calculated. Finally, an error discussion of the whole procedure for one test site is given.     

A method for producing digital probabilistic seismic landslide hazard maps

The 1994 Northridge, California, earthquake is the first earthquake for which we have all of the data sets needed to conduct a rigorous regional analysis of seismic slope instability. These data sets include: (1) a comprehensive inventory of triggered landslides, (2) about 200 strong-motion records of the mainshock, (3) 1:24 000-scale geologic mapping of the region, (4) extensive data on engineering properties of geologic units, and (5) high-resolution digital elevation models of the topography. All of these data sets have been digitized and rasterized at 10 m grid spacing using ARC/INFO GIS software on a UNIX computer. Combining these data sets in a dynamic model based on Newmark's permanent-deformation (sliding-block) analysis yields estimates of coseismic landslide displacement in each grid cell from the Northridge earthquake. The modeled displacements are then compared with the digital inventory of landslides triggered by the Northridge earthquake to construct a probability curve relating predicted displacement to probability of failure. This probability function can be applied to predict and map the spatial variability in failure probability in any ground-shaking conditions of interest. We anticipate that this mapping procedure will be used to construct seismic landslide hazard maps that will assist in emergency preparedness planning and in making rational decisions regarding development and construction in areas susceptible to seismic slope failure.     

South Australian earthquakes, 1980–92

The seismicity of South Australia over the period 1980–92 is presented as a follow‐up to earlier studies. The South Australian seismic network has undergone a significant expansion in the last decade, and with it an increase in the number and precision of located earthquakes. The distribution of recent seismic activity is similar to the historical pattern of earthquakes and the previous instrumental seismicity maps, all of which show the three main areas as being the Flinders‐Mt Lofty Ranges, Eyre Peninsula, and the southeast. The one notable exception in the recent study is the presence of earthquake activity in the Musgrave Block, a previously aseismic region. Intensity characteristics are reported for earthquakes that were sufficiently widely felt. Fault plane solutions for three Flinders Ranges earthquakes (previously unpublished) are also presented; the focal mechanisms are consistent with predominant northeast‐southwest compression. The seismic moment method was used to estimate the seismic risk for the major population centres in terms of probability of exceedance of seismic intensity within a given period. These estimates are based on the recurrence parameters and intensity attenuation function for the region. The results place Adelaide close to the AS2121 ‐ 1979 Earthquake Code Zone I/Zone 2 boundary.     

Probabilistic Assessment of Earthquake Insurance Rates for Turkey

A probabilistic model is presented to obtain a realistic estimate of earthquake insurance rates for reinforced concrete buildings in Turkey. The model integrates information on seismic hazard and information on expected earthquake damage on engineering facilities in a systematic way, yielding to estimates of earthquake insurance premiums. In order to demonstrate the application of the proposed probabilistic method, earthquake insurance rates are computed for reinforced concrete buildings constructed in five cities located in different seismic zones of Turkey. The resulting rates are compared with the rates currently charged by the insurance companies. The earthquake insurance rates are observed to be sensitive to the assumptions on seismic hazard and damage probability matrices and to increase significantly with increasing violation of the code requirements.     

Seismic vulnerability assessment for Montreal

In Canada, Montreal is the second city with the highest seismic risk. This is due to its relatively high seismic hazard, old infrastructures and high population density. The region is characterised by moderate seismic activity with no recent record of a major earthquake. The lack of historical strong ground motion records for the region contributes to large uncertainties in the estimation of hazards. Among the sources of uncertainty, the attenuation function is the main contributor and its effect on estimates of risks is investigated. Epistemic uncertainty was considered by obtaining damage estimates for three attenuation functions that were developed for Eastern North America. The results indicate that loss estimates are highly sensitive to the choice of the attenuation function and suggest that epistemic uncertainty should be considered both for the definition of the hazard function and in loss estimation methodologies. Seismic loss estimates are performed for a 2% in 50 years seismic threat, which corresponds to the design level earthquake in the national building code of Canada, using HAZUS-MH4 for the Montreal region over 522 census tracts. The study estimated that for the average scenario roughly 5% of the building stock would be damaged with direct economic losses evaluated at 1.4 billion dollars for such a scenario. The maximum number of casualties would result in approximately 500 people being injured or dead at a calculated time of occurrence of 2?pm.     

Effect of Aftershocks on Earthquake Hazard Estimation: An Example from the North Anatolian Fault Zone

In order to investigate the effect of aftershocks on earthquake hazard estimation, earthquake hazard parameters (_m, b and M_max) have been estimated by the maximum likelihood method from the main shocks catalogue and the raw earthquakes catalogue for the North Anatolian Fault Zone (NAFZ). The main shocks catalogue has been compiled from the raw earthquake catalogue by eliminating the aftershocks using the window method. The raw earthquake catalogue consisted of instrumentally detected earthquakes between 1900 and 1992, and historical earthquakes that occurred between 1000–1900. For the events of the mainshock catalogue the Poisson process is valid and for the raw earthquake catalogue it does not fit. The paper demonstrates differences in the hazard outputs if on one hand the main catalogues and on the other hand the raw catalogue is used. The maximum likelihood method which allows the use of the mixed earthquake catalogue containing incomplete (historical) and complete (instrumental) earthquake data is used to determine the earthquake hazard parameters. The maximum regional magnitude (M_max, the seismic activity rate (_m), the mean return period (R) and the b value of the magnitude-frequency relation have been estimated for the 24°–31° E, 31°–41° E, 41°–45° E sections of the North Anatolian Fault Zone from the raw earthquake catalogue and the main shocks catalogue. Our results indicate that inclusion of aftershocks changes the b value and the seismic activity rate _m depending on the proportion of aftershocks in a region while it does not significantly effect the value of the maximum regional magnitude since it is related to the maximum observed magnitude. These changes in the earthquake hazard parameters caused the return periods to be over- and underestimated for smaller and larger events, respectively.     

Probabilistic assessment of seismic hazard of dam sites in Jordan

A Probabilistic method is used to evaluate the seismic hazard of nineteen embankment dam sites in Jordan. A line source model developed by McGuire (1978) is used in this study. An updated earthquake catalogue covering the period from 1 A.D. to 1991 A.D. is used for this purpose. This catalogue includes all earthquakes that occurred in Jordan and adjacent areas, more specifically between latitudes 27.0°–35.5° N and longitudes 32.0°–39.0° E.Nine distinct seismic sources of potential seismic activities are identified. The seismic hazard parameters are determined using the method suggested by Kijko and Sellevoll (1989).The Peak Ground Acceleration (PGA) is selected as a measure of ground motion severity. Esteva (1974) attenuation relationship is used in evaluating PGA values at each dam site. Analysis is carried out for 50%, 90%, and 95% probability that is not being exceeded in a life time of 50, 100, and 200 years.Results of analysis indicate that PGA values are higher for dam sites closer to the Dead Sea Fault. This fault is believed to be responsible for most earthquake activities in Jordan and vicinity. The highest PGA value is found to be for Al-Karama dam site.     

Seismic risk assessment of buildings in urban areas: a case study for Denizli,Turkey

This study aims to carry out a seismic risk assessment for a typical mid-size city based on building inventory from a field study. Contributions were made to existing loss estimation methods for buildings. In particular, a procedure was introduced to estimate the seismic quality of buildings using a scoring scheme for the effective parameters in seismic behavior. Denizli, a typical mid-size city in Turkey, was used as a case study. The building inventory was conducted by trained observers in a selected region of Denizli that had the potential to be damaged from expected future earthquakes according to geological and geotechnical studies. Parameters that are known to have some effect on the seismic performance of the buildings during past earthquakes were collected during the inventory studies. The inventory includes data of about 3,466 buildings on 4,226 parcels. The evaluation of inventory data provided information about the distribution of building stock according to structural system, construction year, and vertical and plan irregularities. The inventory data and the proposed procedure were used to assess the building damage, and to determine casualty and shelter needs during the M6.3 and 7.0 scenario earthquakes, representing the most probable and maximum earthquakes in Denizli, respectively. The damage assessment and loss studies showed that significant casualties and economic losses can be expected in future earthquakes. Seismic risk assessment of reinforced concrete buildings also revealed the priorities among building groups. The vulnerability in decreasing order is: (1) buildings with 6 or more stories, (2) pre-1975 constructed buildings, and (3) buildings with 3–5 stories. The future studies for evaluating and reducing seismic risk for buildings should follow this priority order. All data of inventory, damage, and loss estimates were assembled in a Geographical Information System (GIS) database.     

Near-field probabilistic seismic hazard analysis with characteristic earthquake effects

In conventional seismic hazard analysis, uniform distribution over area and magnitude range is assumed for the evaluation of source seismicity which is not able to capture peculiar characteristic of near-fault ground motion well. For near-field hazard analysis, two important factors need to be considered: (1) rupture directivity effects and (2) occurrence of scenario characteristic ruptures in the nearby sources. This study proposed a simple framework to consider these two effects by modifying the predictions from the conventional ground motion model based on pulse occurrence probability and adjustment of the magnitude frequency distribution to account for the rupture characteristic of the fault. The results of proposed approach are compared with those of deterministic and probabilistic seismic hazard analyses. The results indicate that characteristic earthquake and directivity consideration both have significant effects on seismic hazard analysis estimates. The implemented approach leads to results close to deterministic seismic hazard analysis in the short period ranges (T < 1.0 s) and follows probabilistic seismic hazard analysis results in the long period ranges (T > 1.0 s). Finally, seismic hazard maps based on the proposed method could be developed and compared with other methods.     

Seismic hazard assessment in Jordan and its vicinity

The earthquake hazard in Jordan and its vicinity is assessed on the basis of probabilistic methods. For this purpose, an updated earthquake catalog is compiled which covers the period between AD 1–1989. The earthquakes lie between latitudes 27.0°-35.5° N and longitudes 32.0°-39.0° E. Thirteen seismic zones are defined on a regional seismic and tectonic map presented for the area. Point-source and line-source models are used. The seismic hazard parameters, namely, theb-parameter (of the Gutenberg-Richter relation),m ₁ (the upper bound magnitude), and ₄ (the annual rate of occurrence of earthquakes with local magnitudeM _L 4.0) are calculated for each zone. The results of the seismic hazard assessment are displayed as iso-acceleration contours expected to be exceeded during typical economic life times of structures, i.e. 50 and 100 years. For each model, two seismic hazard maps are derived. In order to determine the importance of the South-eastern Mediterranean zone and the north part of the Red Sea zone from a seismic hazard point of view for Jordan, one seismic hazard map which corresponds to 50 years' economic life for every model, excluding the seismicity of these zones, is derived.     

Spatial variation of seismicity parameters across India and adjoining areas

An attempt has been made to quantify the variability in the seismic activity rate across the whole of India and adjoining areas (0–45°N and 60–105°E) using earthquake database compiled from various sources. Both historical and instrumental data were compiled and the complete catalog of Indian earthquakes till 2010 has been prepared. Region-specific earthquake magnitude scaling relations correlating different magnitude scales were achieved to develop a homogenous earthquake catalog for the region in unified moment magnitude scale. The dependent events (75.3%) in the raw catalog have been removed and the effect of aftershocks on the variation of b value has been quantified. The study area was divided into 2,025 grid points (1°×1°) and the spatial variation of the seismicity across the region have been analyzed considering all the events within 300 km radius from each grid point. A significant decrease in seismic b value was seen when declustered catalog was used which illustrates that a larger proportion of dependent events in the earthquake catalog are related to lower magnitude events. A list of 203,448 earthquakes (including aftershocks and foreshocks) occurred in the region covering the period from 250 B.C. to 2010 A.D. with all available details is uploaded in the website .     

Potential of future seismogenesis in Hebei Province (NE China) due to stress interactions between strong earthquakes

Northeast China, a densely populated area, is affected by intense seismic activity, which includes large events that caused extensive disaster and tremendous loss of life. For contributing to the continuous efforts for seismic hazard assessment, the earthquake potential from the active faults near the cities of Zhangjiakou and Langfang in Hebei Province is examined. We estimate the effect of the coseismic stress changes of strong (M ⩾ 5.0) earthquakes on the major regional active faults, and mapped Coulomb stress change onto these target faults. More importantly our calculations reveal that positive stress changes caused by the largest events of the 1976 Tangshan sequence make the Xiadian and part of Daxing fault, thus considered the most likely sites of the next strong earthquake in the study area. The accumulated static stress changes that reached a value of up to 0.4 bar onto these faults, were subsequently incorporated in earthquake probability estimates for the next 30 years.     

A procedure to estimate the seismic hazard in an urban area: an application to Acireale (Eastern Sicily)

A “standard procedure” to characterize the seismic hazard of a given area was proposed. It is based on a multidisciplinary approach implying: (1) the knowledge of the seismic history of the area; (2) detailed geological surveys; (3)seismic noise measurements; (4) simulations of earthquake scenarios. The downtown of Acireale, a typical baroque town located on Eastern Sicily, was chosen as the “test area”. A catalog of the local seismogenic faults (able to generate earthquakes in historical times) has been compiled, as well as a seismic catalog for the effects of both local and regional earthquakes. The analysis of both catalogs allowed us to make the following conclusions: (1) the most important seismogenic faults affecting the Acireale municipality do not affect the downtown, while the related local earthquakes attenuate their energy (and intensity) in short (few km) distances; (2) the highest seismic intensity (degree X) experienced in Acireale downtown was caused by the 1693 regional earthquake; (3) over the last 140 years, the downtown has experienced the highest intensity value of VII only once, while six times the intensity was VI. On the whole, this implies a moderate seismic hazard. The estimation of the seismic hazard has been also approached by the experimental method of recording seismic noise. Measurements have been performed at seven different sites, where drills gave detailed information on the shallow subsurface geology to obtain HV (horizontal/vertical) spectral ratios. On the whole, the highest site amplification factor was moderate (about 7). A further investigation based on synthetic seismograms (and spectra) produced by simulating two given earthquake scenarios was also performed. The two scenarios are, respectively, representative of the largest expected earthquake in the area (the 1693 shock) and of a moderate (magnitude ca. 5.5) local earthquake (as the 1818 one). Moderate to strong locally expected accelerations were evidenced.     

Perceptible earthquakes in the broad Aegean area

A probabilistic estimate of seismic hazard can be obtained from the spatial distribution, of earthquake sources, their frequency–magnitude distribution and the rate of attenuation of strong ground motion with distance. We calculate the earthquake perceptibility, i.e. the annual probability that a particular level of ground shaking will be generated by earthquakes of particular magnitude, by weighting frequency–magnitude data with the predicted felt area for a given level of ground shaking at a particular magnitude. This provides an earthquake selection criterion that can be used in the anti-seismic design of non-critical structures. We calculate the perceptibility, at a particular value of isoseismal intensity, peak ground acceleration and velocity, as a function of source magnitude and frequency for the broad Aegean area using local attenuation laws. We use frequency–magnitude distributions that were previously obtained by combining short-term catalogue data with tectonic moment rate data for 14 tectonic zones in Greece with sufficient earthquake data, and where contemporary strain rates are available from satellite data. Many of the zones show a ‘characteristic earthquake’ distribution with the most perceptible earthquake equal to the maximum magnitude earthquake, but a relatively flat perceptibility between magnitudes 6 and 7. The maximum perceptible magnitude is in the fastest-deforming region in the middle of the Aegean sea, and tends to be systematically low on the west in comparison to the east of the Aegean sea. The tectonic data strongly constrain the long-term recurrence rates and lead to low error estimates (±0.2) in the most perceptible magnitudes.     

大连开发区场地地震设计参数

根据工作区潜在震源区划分、潜在震源区地震活动性参数和地震动衰减关系, 进行场地地震危险性分析, 得到不同概率水平下场区相应地震烈度和基岩水平加速度峰值及其反应谱。根据场地工程地震条件划分不同地质单元及相应的场地类别, 进行不同概率水平的地震反应分析计算, 确定地震动设计参数。     

Seismic Hazard and Loss Estimation for Central America

A new methodology of seismic hazard and loss estimation has been proposed by Chen et al. (Chen et al., 1998; Chan et al., 1998) for the study of global seismic risk. Due to its high adaptability for regions of different features and scales, the methodology was applied to Central America. Seismic hazard maps in terms of both macro-seismic intensity and peak ground acceleration (PGA) at 10% probability of exceedance in 50 years are provided. The maps are all based on the global instrumental as well as historical seismic catalogs and available attenuation relations. Employing the population-weighted gross domestic product (GDP) data, the expected earthquake loss in 50 years for Central America is also estimated at a 5' latitude × 5' longitude resolution. Besides the seismic risk index, a measure of the relative loss or risk degree is calculated for each individual country within the study area. The risk index may provide a useful tool to help allocations of limited mitigation resources and efforts for the purpose of reduction of seismic disasters. For expected heavy loss locations, such as the Central American capital cities, earthquake scenario analysis is helpful in providing a quick overview of loss distribution assuming a major event occurs there. Examples of scenario analysis are given for San Jose, capital of Costa Rica, and Panama City, capital of Panama, respectively.     

Evaluating the seismic hazard to the National Capital (Delhi) Region,India, from moderate earthquakes using simulated accelerograms

The National Capital Region (NCR) of India is exposed to high seismic hazard and risk due to a great earthquake in the central seismic gap of Himalaya and/or due to moderate-size earthquake within NCR. The high population density, rapid growth of infrastructure, and old engineering structures in the region make it more vulnerable to the human as well as economic loss due to moderate-size earthquakes also. The evaluation of seismic hazard is the first step to prepare a proper mitigation plan for the region. The aim of this paper is to evaluate the seismic hazard and risk due to moderate-size earthquakes in the vicinity of NCR. For this purpose, a suit of accelerograms have been generated from hypothetical moderate-size earthquakes (M 5.5 and 6.0) in the region. A basic fault-plane solution is assumed for this purpose. The ranges of the different parameters like depth of focus and stress drop values have been used in order to examine the effect of these parameters on hazard. The accelerograms have been synthesized using two basic velocity models, one representing a hard site and the other a site with a significant low-velocity cover. These two velocity models represent the ridge area and trans-Yamuna river area in the NCR. The decay of peak ground acceleration (pga) values with distance is fast in trans-Yamuna region (with low-velocity surface layer of 100 m) as compared to that of ridge area (with low-velocity surface layer of 1 m). Also, the decay of pga becomes slower if we increase the depth of focus from 10 to 20 km. The response spectra (5% damping) of the synthetic accelerograms for the three periods T = 0.4s, 0.75s, and 1.25s have been estimated and presented in the form of decay curves. The amplifications as a function of epicentral distance and stress drop have also been estimated. We note that the amplifications in 100-m layer case do not occur uniformly at all the distances, rather it is dependent on the angle of incidence of energy into the layers. The pga values are generally amplified by more than twice with increase in stress drop from 100 to 400 bars. The seismic exposure of the population in Delhi city has been presented. The results presented in this study may serve as an important input in the planning of mitigation and disaster management programs in the National Capital Region.     

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.