Seismic hazard in Northern Algeria using spatially smoothed seismicity. Results for peak ground acceleration

Seismic hazard in terms of peak ground acceleration (PGA) has been evaluated in northern Algeria using spatially smoothed seismicity data. We present here a preliminary seismic zoning in northern Algeria as derived from the obtained results.Initially, we have compiled an earthquake catalog of the region taking data from several agencies. Afterwards, we have delimited seismic areas where the b and m_max parameters are different. Finally, by applying the methodology proposed by Frankel [Seismol. Res. Lett. 66 (1995) 8], and using four complete and Poissonian seismicity models, we are able to compute the seismic hazard maps in terms of PGA with 39.3% and 10% probability of exceedance in 50 years.A significant result of this work is the observation of mean PGA values of the order of 0.20 and 0.45 g, for return periods of 100 and 475 years, respectively, in the central area of the Tell Atlas.     

Earthquake Hazard Assessment in the Oran Region (Northwest Algeria)

This paper deals with the probabilistic seismic hazard analysis carried out in the Oran region, situated in the Northwest of Algeria. This part of Algeriawas historically struck by strong earthquakes. It was particularly affected during theOctober 9, 1790 Oran earthquake of intensity X. The main purpose of this work is to assessseismic hazard on rocks in order to provide engineers and planners with a basic tool for seismicrisk mitigation. The probabilistic approach is used in order to take into account uncertaintiesin seismic hazard assessment. Seismic sources are defined in the light of the most recentresults obtained from seismotectonics analyses carried out in North Algeria.Source parameters such as b-values, slip rate and maximum magnitude are assessed for eachseismic source. The attenuation of ground shaking motion with distance is estimated byusing attenuation relationships developed elsewhere throughout the world (Sadigh et al., 1993; Ambraseys and Bommer, 1991). The two relationships agree well with the local data. Differentchoices of source parameter values and attenuation relationships are assigned weights in alogic tree model. Results are presented as relationships between values of peak groundacceleration (PGA) and annual frequency of exceedance, and maps of hazard for returnperiods of 200 years and 500 years. A maximum peak ground acceleration of 0.42 g is obtainedfor the Oran site for a return period of 500 years.     

Seismic hazard evaluation at Algiers using Benouar's earthquake catalogue

This paper presents the evaluation of seismic hazard at the site of Algiers (capital of Algeria). In order to implement earthquake-resistant design codes, it is usually necessary to know the maximum dynamic load which a particular structure might experience during its economic life, or alternatively, the most probable return period of a specified design load. The evaluation of the seismic hazard at the site, based on peak ground motion acceleration and using Cornell's method and Benouar's earthquake Maghreb catalogue, in terms of return period, probability of exceedance of PGA, design ground motion and a response spectrum, is carried out for the City of Algiers and its surroundings. The response spectrum for Algiers presented in this paper is the first one realized in Algeria using revised Algerian data.     

Some aspects of the seismic hazard in Panama City

The definition of earthquake sources in the Panama region on the basis of both tectonics and average seismicity rates, have recently led to the concept of a microplate surrounded by seismically active areas. The effects of these earthquakes on the place where the most important concentration of investments and population is located, the capital city of Panama, are analyzed in this paper using statistical approaches.The parameters of Gumbel's Type-I distribution of extreme values for a continuous interval of 60 yr annual maximum magnitudes were used to make probabilistic estimations of the seismic hazard in Panama City. An earthquake with magnitude 7.5 is capable of producing a modified Mercalli intensity VII in Panama City, provided the source distance is of the order of 100 km. This earthquake has a probability of occurrence of 69% in 50 yr.     

Seismic Hazard Analysis for the Bangalore Region

Indian peninsular shield, which was once considered to be seismically stable, is experiencing many earthquakes recently. As part of the national level microzonation programme, Department of Science and Technology, Govt. of India has initiated microzonation of greater Bangalore region. The seismic hazard analysis of Bangalore region is carried out as part of this project. The paper presents the determination of maximum credible earthquake (MCE) and generation of synthetic acceleration time history plot for the Bangalore region. MCE has been determined by considering the regional seismotectonic activity in about 350 km radius around Bangalore city. The seismotectonic map has been prepared by considering the faults, lineaments, shear zones in the area and historic earthquake events of more than 150 events. Shortest distance from the Bangalore to the different sources is measured and then peak ground acceleration (PGA) is calculated for the different source and moment magnitude. Maximum credible earthquake found in terms of moment magnitude is 5.1 with PGA value of 0.146 g at city centre with assuming the hypo central distance of 15.88 km from the focal point. Also, correlations for the fault length with historic earthquake in terms of moment magnitude, yields (taking the rupture fault length as 5% of the total fault length) a PGA value of 0.159 g. Acceleration time history (ground motion) and a response acceleration spectrum for the corresponding magnitude has been generated using synthetic earthquake model considering the regional seismotectonic parameters. The maximum spectral acceleration obtained is 0.332 g for predominant period of 0.06 s. The PGA value and synthetic earthquake ground motion data from the identified vulnerable source using seismotectonic map will be useful for the PGA mapping and microzonation of the area.     

Seismic Hazard and Risk Assessment for Tulbagh,South Africa: Part I – Assessment of Seismic Hazard

In this part of our study the probabilistic seismic hazard analysis (PSHA) for Tulbagh was performed. The applied procedure is parametric and consists essentially of two steps. The first step is applicable to the area in the vicinity of Tulbagh and requires an estimation of the area-specific parameters, which, in this case, is the mean seismic activity rate, , the Gutenberg-Richter parameter, b, and the maximum regional magnitude, m_max. The second step is applicable to the Tulbagh site, and consists of parameters of distribution of amplitude of the selected ground motion parameter. The current application of the procedure provides an assessment of the PSHA in terms of peak ground acceleration (PGA) and spectral acceleration (SA). The procedure permits the combination of both historical and instrumental data. The historical part of the catalogue only contains the strongest events, whereas the complete part can be divided into several subcatalogues, each assumed complete above a specified threshold of magnitude. In the analysis, the uncertainty in the determination of the earthquake was taken into account by incorporation of the concept of `apparent magnitude'. The PSHA technique has been developed specifically for the estimation of seismic hazard at individual sites without the subjective judgement involved in the definition of seismic source zones, when the specific active faults have not been mapped or identified, and where the causes of seismicity are not well understood. The results of the hazard assessment are expressed as probabilities that specified values of PGA will be exceeded during the chosen time intervals, and similarly for the spectral accelerations. A worst case scenario sketches the possibility of a maximum PGA of 0.30g. The results of the hazard assessment can be used as input to a seismic risk assessment.     

Multiple source and attenuation relationships for evaluation of deterministic seismic hazard: logic tree approach considering local site effects

The most important seismic hazard parameters required to demarcate seismic zones are the peak horizontal acceleration (PHA) and spectral acceleration (SA). The two approaches for evaluation of seismic hazard are the probabilistic seismic hazard analysis and the deterministic seismic hazard analysis (DSHA). The present study evaluates the seismic hazard of the South Indian Peninsular region based on the DSHA methodology. In order to consider the epistemic uncertainties in a better manner, a logic tree approach was adopted in the evaluation of seismic hazard. Two types of seismic sources and three different attenuation relations were used in the analysis. The spatial variation of PHA (mean and 84th percentile values) and SA values for 1 Hz and 10 Hz at bedrock level (84th percentile values) for the entire study area were evaluated and the results are presented here. The surface level peak ground acceleration (PGA) values will be different from that of the bedrock level values due to the local site conditions. The PGA values at ground surface level were evaluated for four different National Earthquake Hazard Reduction Program site classes by considering the non-linear site response of different soil types. The response spectra for important cities in South India were also prepared using the deterministic approach and the results are presented in this paper.     

Global Attenuation Relationship for Estimating Peak Ground Acceleration

Peak Ground Acceleration (PGA) is a very important ground motion parameter which is used to define the degree of ground shaking during an earthquake. It is also very helpful for designing earthquake resistant structure. The PGA can be estimated by attenuation relationships using magnitude, distance, source type etc of a ground motion. In the past, several researchers have developed over 450 attenuation relationships for predicting PGA for a specific region. In the present study an attempt has been made to develop an attenuation relationship on the basis of these available previous relationships in rock site which will be applicable for any region of the world. In the present study, PGA has been expressed as a function of moment magnitude and hypo-central distance in rock site. Chi-square test have also been performed with available earthquake data in American and Indian region for verifying the accuracy of the generated attenuation relationship. Using multiple regression and Genetic Algorithm (GA) the attenuation relationship equations have also been generated. These equations will be very helpful for performing seismic hazard analysis and predicting earthquake force in any region of the world.     

First Order Seismic Microzonation of Delhi,India Using Geographic Information System (GIS)

A first order seismic microzonation map of Delhi is prepared using five thematic layers viz., Peak Ground Acceleration (PGA) contour, different soil types at 6 m depth, geology, groundwater fluctuation and bedrock depth, integrated on GIS platform. The integration is performed following a pair-wise comparison of Analytical Hierarchy Process (AHP), wherein each thematic map is assigned weight in the 5-1 scale: depending on its contribution towards the seismic hazard. Following the AHP, the weightage assigned to each theme are: PGA (0.333), soil (0.266), geology (0.20), groundwater (0.133) and bedrock depth (0.066). The thematic vector layers are overlaid and integrated using GIS. On the microzonation theme, the Delhi region has been classified into four broad zones of vulnerability to the seismic hazard. They are very high (> 52%), high (38–52%), moderate (23–38%) and less ( < 23%) zones of seismic hazard. The “very high” seismic hazard zone is observed where the maximum PGA varies from 140 to 210 gal for a finite source model of M_w 8.5 in the central seismic gap. A site amplification study from local and regional earthquakes for Delhi region using Delhi Telemetry Network data shows a steeper site response gradient in the eastern side of the Yamuna fluvial deposits at 1.5 Hz. The ‘high’ seismic hazard zone occupies most of the study area where the PGA value ranges from 90 to 140 gal. The ‘moderate’ seismic hazard zone occurs on either side of the Delhi ridge with PGA value varying from 60 to 90 gal. The ‘less’ seismic hazard zone occurs in small patches distributed along the study area with the PGA value less than 60 gal. Site response studies, PGA distribution and destruction pattern of the Chamoli earthquake greatly corroborate the seismic hazard zones estimated through microzonation on GIS platform and also establishes the methodology incorporated in this study.     

Probabilistic seismic hazard assessment for the Horn of Africa based on seismotectonic regionalisation

A probabilistic seismic hazard analysis (PSHA) for the Horn of Africa is presented. Our seismicity database consists of a revised and up-to-date regional catalogue compiled from different agencies, checked for completeness with respect to time and homogenized with respect to magnitude (M_s). The seismic source zones are based on our present day knowledge of the regional seismotectonics. Among the results we present regional hazard maps for 0.01 annual probability for intensity and Peak Ground Acceleration (PGA) and hazard curves and response spectra for six economical significant sites within the region. The model uncertainties with respect to seismicity are analysed in a novel approach and form part of a sensitivity analysis that quantifies our PSHA modelling uncertainties.

For 0.01 annual probability we find randomly oriented horizontal PGA that exceed just 0.2 g and MM-scale intensity VIII in the Afar depression and southern Sudan. Uncertainties amount to 20% g PGA in some cases, mainly due to attenuation uncertainties. Intensity uncertainties seldom exceed 0.5 intensity units. Relatively large seismic hazard is found for Djibouti (VIII for 0.01 annual probability), slightly lower for the port of Massawa (between VII and VIII for 0.01 annual probability) and low for the port of Assab (between VI and VII for 0.01 annual probability).     

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.     

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.     

Application of the Spatially Smoothed Seismicity and Monte Carlo Methods to Estimate the Seismic Hazard of Eritrea and the Surrounding Region

The region of interest is characterized by incomplete data sets and little information about the tectonic features. Therefore, two methodologies for estimating seismic hazard were used in order to elucidate the robustness of the results: the method of spatially smoothed seismicity introduced by Frankel (1995) and later extended by Lapajne et al. (1997) and a Monte Carlo approach presented by Ebel and Kafka (1999). In the first method, fault-rupture oriented elliptical Gaussian smoothing was performed to estimate future activity rates along the causative structures. Peak ground accelerations were computed for a grid size of 15 km × 415 km assuming the centre of the grids as epicentres, from which the seismic hazard map was produced. The attenuation relationship by Ambraseys et al. (1996) was found suitable for the region under study. PGA values for 10% probability of exceedence in 50 years (return period of 475 years) were computed for each model and a combined seismic hazard map was produced by subjectively assigning weights to each of these models. A worst-case map is also obtained by picking the highest value at each grid point from values of the four hazard maps. The Monte Carlo method is used to estimate seismic hazard, for comparison to the results from our previous approach. Results obtained from both methods are comparable except values in the first set of maps estimate greater hazard in areas of low seismicity. Both maps indicate a higher hazard along the main tectonic features of the east African and Red Sea rift systems. Within Eritrea, the highest PGA exceeded a value 25% of g, located north of Red Sea port of Massawa. In areas around the capital, Asmara, PGA values exceed 10% of g.     

An estimate of probabilistic seismic hazard for five cities in Greece by using the parametric-historic procedure

A probabilistic procedure was applied to assess seismic hazard for the sites of five Greek cities (Athens, Heraklion, Patras, Thessaloniki and Volos) using peak ground acceleration as the hazard parameter. The methodology allows the use of either historical or instrumental data, or a combination of both. It has been developed specifically for the estimation of seismic hazard at a given site and does not require any specification of seismic sources or/and seismic zones. A new relation for the attenuation of peak ground acceleration was employed for the shallow seismicity in Greece. The computations involved the area- and site-specific parts. When assessing magnitude recurrence for the areas surrounding the five cities, the maximum magnitude, m_max, was estimated using a recently derived equation. The site-specific results were expressed as probabilities that a given peak ground acceleration value will be exceeded at least once during a time interval of 1, 50 and 100 years at the sites of the cities. They were based on the maximum peak ground acceleration values computed by assuming the occurrence of the strongest possible earthquake (of magnitude m_max) at a very short distance from the site and using the mean value obtained with the help of the attenuation law. This gave 0.24 g for Athens, 0.53 g for Heraklion (shallow) and 0.39 g Heraklion (intermediate-depth seismicity), 0.30 g for Patras, 0.35 g for Thessaloniki and 0.30 g for Volos. In addition, the probabilities of exceedance of the estimated maximum peak ground acceleration values were calculated for the sites. The standard deviation of the new Greek attenuation law demonstrates the uncertainty and large variation of predicted peak ground acceleration values.     

Seismic hazard in Greece. II. Ground acceleration

In a previous paper (Makropoulos and Burton, 1985) the seismic hazard in Greece was examined in terms of magnitude recurrence using Gumbel's third asymptotic distribution of extreme values and concepts of the physical process of strain energy release. The present study extends the seismic hazard methods beyond magnitude to the estimation of expectations of levels of peak ground acceleration exceedance thus allowing for a direct comparison between these two methodologies as well as establishing information relevant to design and planning criteria.The limited number of strong motion records do not permit regional study of attenuation of ground vibration in Greece. An average formula is derived from eight well known formulae which resulted from worldwide studies, this is: a = 2164 e^0.70m (r+20)^−1.80 cm s⁻² where a is peak ground acceleration, m is earthquake magnitude and r is hypocentral distance in kilometres. This formula agrees with the observed values of peak ground acceleration values recorded in Greece.Acceleration seismic hazard is calculated at each of six chosen cities. Values of maximum acceleration with probability 70% of not been exceeded in the next 25, 50, 100, and 200 years are obtained along with corresponding values of velocity and displacement. The same detailed acceleration evaluation is then applied to the whole area of Greece by dividing it into cells of 0.5° lat × 0.5° long, and the results are illustrated through isoacceleration maps.Differences in magnitude and acceleration hazard maps reflect the fact that in acceleration hazard assessment the focal distance from a particular place in an important factor. The cities of Heraklion and Rodhos have the lowest acceleration hazard although the expected earthquakes may have large magnitude. Intermediate depth earthquakes characterise these two cities. Acceleration estimates, unlike magnitude hazard parameters, refer to a particular place and not to an area around it. Hence, even if two places have similar earthquake depth distributions, the hazards may differ significantly because of the different spatial distribution of the foci. This is observed in the case of Athens and Corinth. These cities have almost the same magnitude hazard, but the acceleration hazard is much lower for Athens where the hazard is mainly due to more distant earthquakes.The isoacceleration maps for Greece as a whole also define areas of high seismic hazard. These are the areas around Cephalonia and Leukas Islands in the Ionian Sea and the eastern Sporadhes, Lesbos Islands and Chalkidiki in the Northern Aegean Sea. At the 70% probability level the maximum acceleration is expected to be around 0.2g within the next 50 years. The areas where the maximum acceleration at the 70% probability level is expected to reach a value of 0.3g in the next 200 years are around Cephalonia and Leukas Islands and near the Dardanelles.     

Seismic Hazard Assessment for the Taiwan Region on the Basis of Recent Strong-Motion Data and Prognostic Zonation of Future Earthquakes

The paper describes an integrated approach to seismic hazard assessment, which was applied for the Taiwan region. First, empirical modelsfor ground motion estimation in the region were obtained on the basisof records from recent (1993-1999) earthquakes. The databaseincludes strong-motion data collected during the recent Chi-Chiearthquake (M=7.6, 21 September 1999) and large (M=6.8)aftershocks. The ground-motion database was also used for evaluationof generalised site amplification functions for typical soil classes(B, C and D). Second, the theoretical seismic catalogue (2001–2050)for the Taiwan region had been calculated using the 4D-model(location, depth, time) for dynamic deformation of the Earth' crustand 5D-model (location, depth, time, magnitude) for seismic process.The models were developed on the basis of available geophysical andgeodynamic data that include regional seismic catalogue. Third, theregion & site & time-dependent seismic analysis, which is basedon schemes of probable earthquake zones evaluated from the theoreticalcatalogue, regional ground motion models, and local site responsecharacteristics, has been performed. The seismic hazard maps arecompiled in terms of Peak Ground Acceleration (PGA) and ResponseSpectra (RS) amplitudes. The maps show distribution of amplitudesthat will not be exceeded with certain probability in condition oftypical soil classes during all possible earthquakes that may occur inthe region during time period of 2003–2025. The approach allowsintroducing new parameter that describes dependency of seismichazard on time, so-called 'period of maximum hazard'. Theparameter shows the period, during which every considered sitewill be subjected by the maximum value of ground motioncharacteristic (PGA or RS).     

An overview on the seismic zonation and microzonation studies in India

The present study presents a review on the progressive development of the seismic zonation map of India both from official agencies and also from independent individual studies. The zonation map have been modified and updated regularly with the occurrence of major destructive earthquakes over the years in the Indian subcontinent with the addition of new data. This study discusses the criteria chosen for the progressive zonation and the major earthquakes that were responsible for retrospection of the earlier published maps. The seismic zonation maps of India have also been prepared by various independent workers by adopting different approaches to achieve the purpose of the zonation. Despite the endeavors from various sources to provide a solution for the problem of earthquake hazards in India, there were many limitations on the zonation map as it gives the picture at a regional scale mostly on the bedrock level without addressing the local site conditions. But nevertheless, the seismic zonation map gives basic guidelines for any region to know the hazard scenario and if any city or urban population is under threat from seismic point of view, further site specific seismic microzonation may be carried out. In the International scenario, the Global Seismic Hazard Assessment Program (GSHAP) in 1999 prepared a hazard map for world in terms of peak ground acceleration (PGA) with a 10% probability of exceedance in 50 years, but it turned out to be an underestimation of the hazard parameter when compared with the observed PGA. To tackle the problem of seismic hazards, there was a need to have a detail study on the local site conditions in terms of its geological, geophysical and geotechnical properties. With the advent of better instrumentation and knowledge on the mechanics of earthquakes, it was possible to identify zones of hazards at a local level and this gives rise to the study of seismic microzonation. Seismic microzonation work has been carried out in India in some of the strategic important mega cities and industrial build up that has the potential of being damaged from future earthquakes, as has been shown in the past. Though the microzonation map is not the final output map, as it can still be updated at later stage with more input data, it does provide a more realistic picture on the site specific seismic hazard.     

Different approaches to the seismic hazard of Sannio-Matese (Southern Italy)

The seismic hazard of Sannio-Matese is calculated using a new seismogenetic zoning of southern Italy of seven areas and different, azimuth-dependent, attenuation laws. Various approaches (Gumbel's first and third asymptotic distribution, Cornell) lead to similar results for the different exposure times considered (100 and 200 years) and probability levels (37 and 68%). The present seismic regulation proposed in 1980 by the Italian National Council for Researches, and based on a different approach, is, in general, confirmed by the results.     

Probabilistic seismic hazard of Pakistan, Azad-Jammu and Kashmir

The seismic hazard study for Pakistan and Azad Jammu and Kashmir has been conducted by using probabilistic approach in terms of peak ground acceleration (PGA) in m/s² and also seismic hazard response spectra for different cities. A new version of Ambraseys et al. (Bull Earthq Eng 3:1–53, 2005) ground acceleration model is used, and parameterization is based on most recent updated earthquake catalogs that consisted of 14,000 events. The threshold magnitude was fixed at M _w 4.8, but seismic zones like northern Pakistan–Tajikistan, Hindukush and northern Afghanistan–Tajikistan border had M _w 5.2. The average normalized ‘a’ and ‘b’ values for all zones are 6.15 and 0.95, respectively. Seismicity of study area was modeled, and ground motion was computed for eight frequencies (0.025, 0.1, 0.2, 0.5, 1.0, 1.5, 2.0, 2.5 s) for different annual exceedance rates of 0.02, 0.01, 0.005, 0.002 and 0.001 (return periods 50, 100, 200, 500 and 1,000 years) for stiff rocks at the gridding of 0.1° × 0.1°. Seismic hazard maps based on computed PGA for 0.02, 0.01 and 0.002 annual exceedance are prepared. These maps indicate the earthquake hazard of Pakistan and surrounding areas in the form of acceleration contour lines, which are in agreement with geological and seismotectonic characteristics of the study area. The maximum seismic hazard values are found at Muzaffarabad, Gilgit and Quetta areas.     

Sensitivity analysis of seismic hazard for Western Liguria (North Western Italy): A first attempt towards the understanding and quantification of hazard uncertainty

The use of logic trees in probabilistic seismic hazard analyses often involves a large number of branches that reflect the uncertainty in the selection of different models and in the selection of the parameter values of each model. The sensitivity analysis, as proposed by Rabinowitz and Steinberg [Rabinowitz, N., Steinberg, D.M., 1991. Seismic hazard sensitivity analysis: a multi-parameter approach. Bull. Seismol. Soc. Am. 81, 796–817], is an efficient tool that allows the construction of logic trees focusing attention on the parameters that have greater impact on the hazard.In this paper the sensitivity analysis is performed in order to identify the parameters that have the largest influence on the Western Liguria (North Western Italy) seismic hazard. The analysis is conducted for six strategic sites following the multi-parameter approach developed by Rabinowitz and Steinberg [Rabinowitz, N., Steinberg, D.M., 1991. Seismic hazard sensitivity analysis: a multi-parameter approach. Bull. Seismol. Soc. Am. 81, 796–817] and accounts for both mean hazard values and hazard values corresponding to different percentiles (e.g., 16%-ile and 84%-ile). The results are assessed in terms of the expected PGA with a 10% probability of exceedance in 50 years for rock conditions and account for both the contribution from specific source zones using the Cornell approach [Cornell, C.A., 1968. Engineering seismic risk analysis. Bull. Seismol. Soc. Am. 58, 1583–1606] and the spatially smoothed seismicity [Frankel, A., 1995. Mapping seismic hazard in the Central and Eastern United States. Seismol. Res. Lett. 66, 8–21]. The influence of different procedures for calculating seismic hazard, seismic catalogues (epicentral parameters), source zone models, frequency–magnitude parameters, maximum earthquake magnitude values and attenuation relationships is considered. As a result, the sensitivity analysis allows us to identify the parameters with higher influence on the hazard. Only these parameters should be subjected to careful discussion or further research in order to reduce the uncertainty in the hazard while those with little or no effect can be excluded from subsequent logic-tree-based seismic hazard analyses.