Earthquake hazard and risk assessment based on unified scaling law for earthquakes: Altai–Sayan Region

We apply the general concept of seismic risk analysis based on morphostructural analysis of the territory, pattern recognition of earthquake-prone nodes, and the Unified Scaling Law for Earthquakes, USLE, in another seismic region of Russia to the west from Lake Baikal, i.e., Altai–Sayan Region. The USLE generalizes the empirical Gutenberg–Richter relationship making use of apparently fractal distribution of earthquake sources of different size: \( \log_{10} N\left( {M,L} \right)\, = \,A\, + \,B \cdot \left( {5\, - \,M} \right)\, + \,C \cdot \log_{10} L, \) where N (M, L) is the expected annual number of earthquakes of a certain magnitude M within an seismically prone area of linear dimension L. The local estimates of A, B, and C allow determination of the expected maximum credible magnitude in a given time interval and the associated spread around ground shaking parameters (e.g., peak ground acceleration, PGA, or macroseismic intensity, I₀). Compilation of the corresponding seismic hazard map of Altai–Sayan Region and its rigorous testing against the available seismic evidences in the past is used to model regional maps of specific earthquake risks for population, cities, and infrastructures.     

Hybrid attenuation model for estimation of peak ground accelerations in the Kutch region, India

The Kutch region of Gujarat in India is the locale of one of the most devastating earthquake of magnitude (M _w) 7.7, which occurred on January 26, 2001. Though, the region is considered as seismically active region, very few strong motion records are available in this region. First part of this paper uses available data of strong motion earthquakes recorded in this region between 2006 and 2008 years to prepare attenuation relation. The developed attenuation relation is further used to prepare synthetic strong motion records of large magnitude earthquakes using semiempirical simulation technique. Semiempirical simulation technique uses attenuation relation to simulate strong ground motion records of any target earthquake. The database of peak ground acceleration obtained from simulated records is used together with database of peak ground acceleration obtained from observed record to develop following hybrid attenuation model of wide applicability in the Kutch region: $$ \begin{aligned} \ln \left( {\text{PGA}} \right) & = - 2.56 + 1.17 \, M_{\text{w}} - \, 0.015R - 0.0001\ln \left( {E + 15} \right) \\ &\quad 3.0 \le M_{\text{w}} \le 8.2;\quad 12 \le R \le 120;\quad {\text{std}} . {\text{ dev}}.(\sigma ): \pm 0.5 \\ \end{aligned} $$ ln ( PGA ) = ? 2.56 + 1.17 M w ? 0.015 R ? 0.0001 ln ( E + 15 ) 3.0 ≤ M w ≤ 8.2 ; 12 ≤ R ≤ 120 ; std . dev . ( σ ) : ± 0.5 In the above equation, PGA is maximum horizontal ground acceleration in gal, M _w is moment magnitude of earthquake, R is hypocentral distance, and E is epicentral distance in km. The standard deviation of residual of error in this relation is 0.5. This relation is compared with other available relations in this region, and it is seen that developed relation gives minimum root mean square error in comparison with observed and calculated peak ground acceleration from same data set. The applicability of developed relation is further checked by testing it with the observed peak ground acceleration from earthquakes of magnitude (M _w), 3.6, 4.0, 4.4, and 7.7, respectively, which are not included in the database used for regression analysis. The comparison demonstrates the efficacy of developed hybrid attenuation model for calculating peak ground acceleration values in the Kutch region.     

Probabilistic seismic hazard at the archaeological site of Gol Gumbaz in Vijayapura,south India

Probabilistic seismic hazard analysis (PSHA) is carried out for the archaeological site of Vijayapura in south India in order to obtain hazard consistent seismic input ground-motions for seismic risk assessment and design of seismic protection measures for monuments, where warranted. For this purpose the standard Cornell-McGuire approach, based on seismogenic zones with uniformly distributed seismicity is employed. The main features of this study are the usage of an updated and unified seismic catalogue based on moment magnitude, new seismogenic source models and recent ground motion prediction equations (GMPEs) in logic tree framework. Seismic hazard at the site is evaluated for level and rock site condition with 10% and 2% probabilities of exceedance in 50 years, and the corresponding peak ground accelerations (PGAs) are 0.074 and 0.142 g, respectively. In addition, the uniform hazard spectra (UHS) of the site are compared to the Indian code-defined spectrum. Comparisons are also made with results from National Disaster Management Authority (NDMA 2010), in terms of PGA and pseudo spectral accelerations (PSAs) at T = 0.2, 0.5, 1.0 and 1.25 s for 475- and 2475-yr return periods. Results of the present study are in good agreement with the PGA calculated from isoseismal map of the Killari earthquake, \({\hbox {M}}_{\mathrm{w}} = 6.4\) (1993). Disaggregation of PSHA results for the PGA and spectral acceleration (\({\hbox {S}}_{\mathrm{a}}\)) at 0.5 s, displays the controlling scenario earthquake for the study region as low to moderate magnitude with the source being at a short distance from the study site. Deterministic seismic hazard (DSHA) is also carried out by taking into account three scenario earthquakes. The UHS corresponding to 475-yr return period (RP) is used to define the target spectrum and accordingly, the spectrum-compatible natural accelerograms are selected from the suite of recorded accelerograms.     

Modeling of strong ground motions for 1991 Uttarkashi, 1999 Chamoli earthquakes, and a hypothetical great earthquake in Garhwal–Kumaun Himalaya

In this study, stochastic finite fault modeling is used to simulate Uttarkashi (1991) and Chamoli (1999) earthquakes using all available source, path, and site parameters available for the region. These two moderate earthquakes are recorded at number of stations of a strong motion network. The predicted peak ground accelerations at these stations are compared with the observed data and the ground motion parameters are constrained. The stress drop of Uttarkashi and Chamoli earthquakes is constrained at 77 and 65?bars, respectively, whereas the quality factor Q _C is 112 $ f^{0.97} $ and 149 $ f^{0.95} $ for these two regions. The high-frequency attenuation parameter Kappa is in the range 0.04?C0.05. The constrained ground motion parameters are then used to simulate Mw 8.5 earthquake in central seismic gap region of Himalaya. Two scenarios are considered with epicenter of future great earthquake at locations of Uttarkashi and Chamoli earthquakes using above constrained parameters. The most vulnerable towns are the towns of Dehradun and Almora where expected PGA is in excess of 600?cm/s² at V_S30 520?m/s when the epicenter of the great earthquake is at the location of Uttarkashi (1991) earthquake. The towns of Shimla and Chandigarh can expect PGA close to 200?cm/s². Whereas when the epicenter of the great earthquake is at the location of Chamoli (1999) earthquake, the towns of Dehradun and Almora can expect PGA of around 500 and 400?cm/s², respectively, at V_S30 620?m/s. The National Capital Region, Delhi can expect accelerations of around 80?cm/s² in both the cases. The PGA contour maps obtained in this study can be used to assess the seismic hazard of the region and identify vulnerable areas in and around central Himalaya from a future great earthquake.     

Assessment of the statistical earthquake hazard parameters for NW Turkey

We performed a probabilistic analysis of earthquake hazard input parameters, NW Turkey covers Gelibolu and Biga Peninsulas, and its vicinity based on four seismic sub-zones. The number of earthquakes with magnitude M ≥ 3.0 occurred in this region for the period between 1912 and 2007 is around 5130. Four seismic source sub-zones were defined with respect to seismotectonic framework, seismicity and fault geometry. The hazard perceptibility characterization was examined for each seismic source zone and for the whole region. The probabilities of earthquake recurrences were obtained by using Poisson statistical distribution models. In order to determine the source zones where strong and destructive earthquakes may occur, distribution maps for a, b and a/b values were calculated. The hazard scaling parameters (generally known as a and b values) in the computed magnitude–frequency relations vary in the intervals 4.28–6.58 and 0.59–1.13, respectively, with a RMS error percentage below 10 %. The lowest b value is computed for sub-zone three indicating the predominance of large earthquakes mostly at Gelibolu (Gallipoli) and north of Biga Peninsula (southern Marmara region), and the highest b value is computed for sub-zone two Edremit Bay (SW Marmara region). According to the analysis of each seismic sub-zone, the greatest risk of earthquake occurrence is determined for the triangle of Gelibolu–Tekirda? western part of Marmara Sea. Earthquake occurrence of the largest magnitude with 7.3 within a 100-year period was determined to be 46 % according to the Poisson distribution, and the estimated recurrence period of years for this region is 50 ± 12. The seismic hazard is pronounced high in the region extending in a NW–SE direction, north of Edremit Bay, west of Saros Bay and Yenice Gönen (southern Marmara region) in the south. High b values are generally calculated at depths of 5–20 km that can be expressed as low seismic energy release and evaluated as the seismogenic zone.     

Probabilistic seismic hazard maps for the sultanate of Oman

This study presents the results of the first probabilistic seismic hazard assessment (PSHA) in the framework of logic tree for Oman. The earthquake catalogue was homogenized, declustered, and used to define seismotectonic source model that characterizes the seismicity of Oman. Two seismic source models were used in the current study; the first consists of 26 seismic source zones, while the second is expressing the alternative view that seismicity is uniform along the entire Makran and Zagros zones. The recurrence parameters for all the seismogenic zones were determined using the doubly bounded exponential distribution except the zones of Makran, which were modelled using the characteristic distribution. Maximum earthquakes were determined and the horizontal ground accelerations in terms of geometric mean were calculated using ground-motion prediction relationships developed based upon seismic data obtained from active tectonic environments similar to those surrounding Oman. The alternative seismotectonic source models, maximum magnitude, and ground-motion prediction relationships were weighted and used to account for the epistemic uncertainty. Hazard maps at rock sites were produced for 5?% damped spectral acceleration (SA) values at 0.1, 0.2, 0.3, 1.0 and 2.0?s spectral periods as well as peak ground acceleration (PGA) for return periods of 475 and 2,475?years. The highest hazard is found in Khasab City with maximum SA at 0.2?s spectral period reaching 243 and 397?cm/s² for return periods 475 and 2,475 years, respectively. The sensitivity analysis reveals that the choice of seismic source model and the ground-motion prediction equation influences the results most.     

Seismic hazard map of Coimbatore using subsurface fault rupture

This study presents the future seismic hazard map of Coimbatore city, India, by considering rupture phenomenon. Seismotectonic map for Coimbatore has been generated using past earthquakes and seismic sources within 300 km radius around the city. The region experienced a largest earthquake of moment magnitude 6.3 in 1900. Available earthquakes are divided into two categories: one includes events having moment magnitude of 5.0 and above, i.e., damaging earthquakes in the region and the other includes the remaining, i.e., minor earthquakes. Subsurface rupture character of the region has been established by considering the damaging earthquakes and total length of seismic source. Magnitudes of each source are estimated by assuming the subsurface rupture length in terms of percentage of total length of sources and matched with reported earthquake. Estimated magnitudes match well with the reported earthquakes for a RLD of 5.2% of the total length of source. Zone of influence circles is also marked in the seismotectonic map by considering subsurface rupture length of fault associated with these earthquakes. As earthquakes relive strain energy that builds up on faults, it is assumed that all the earthquakes close to damaging earthquake have released the entire strain energy and it would take some time for the rebuilding of strain energy to cause a similar earthquake in the same location/fault. Area free from influence circles has potential for future earthquake, if there is seismogenic source and minor earthquake in the last 20 years. Based on this rupture phenomenon, eight probable locations have been identified and these locations might have the potential for the future earthquakes. Characteristic earthquake moment magnitude (M _w ) of 6.4 is estimated for the seismic study area considering seismic sources close to probable zones and 15% increased regional rupture character. The city is divided into several grid points at spacing of 0.01° and the peak ground acceleration (PGA) due to each probable earthquake is calculated at every grid point in city by using the regional attenuation model. The maximum of all these eight PGAs is taken for each grid point and the final PGA map is arrived. This map is compared to the PGA map developed based on the conventional deterministic seismic hazard analysis (DSHA) approach. The probable future rupture earthquakes gave less PGA than that of DSHA approach. The occurrence of any earthquake may be expected in near future in these eight zones, as these eight places have been experiencing minor earthquakes and are located in well-defined seismogenic sources.     

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 hazard assessment of the city of Hamedan and its vicinity, west of Iran

This article presents the results of deterministic and probabilistic seismic hazard analyses (DSHA and PSHA) of the city of Hamedan and its neighboring regions. This historical city is one of the developing cities located in the west of Iran. For this reason, the DSHA and PSHA approaches have been used for the assessment of seismic hazards and earthquake risk evaluation. To this purpose, analyses have been carried out considering the historic and instrumented earthquakes, geologic and seismotectonic parameters of the region covering a radius of 100?km, keeping Hamedan as the center. Therefore, in this research, we studied the main faults and fault zones in the study area and calculated the length and distance of faults from the center of Hamedan. In the next step, we measured the maximum credible earthquake (MCE) and peak ground acceleration (PGA) using both DSHA and PSHA approaches and utilized the various equations introduced by different researchers for this purpose. The results of DSHA approach show that the MCE-evaluated value is 7.2 Richter, which might be created by Nahavand fault activities in this region. The PGA value of 0.56?g will be obtained from Keshin fault. The results of PSHA approach show that the MCE-evaluated value is 7.6 Richter for a 0.64 probability in a 50-year period. The PGA value of 0.45?g will be obtained from Keshin fault. Seismic hazard parameters have been evaluated considering the available earthquake data using Gutenberg?CRichter relationship method. The ??a?? and ??b?? parameters were estimated 5.53 and 0.68, respectively.     

Seismicity of Gujarat

Paper describes tectonics, earthquake monitoring, past and present seismicity, catalogue of earthquakes and estimated return periods of large earthquakes in Gujarat state, western India. The Gujarat region has three failed Mesozoic rifts of Kachchh, Cambay, and Narmada, with several active faults. Kachchh district of Gujarat is the only region outside Himalaya-Andaman belt that has high seismic hazard of magnitude 8 corresponding to zone V in the seismic zoning map of India. The other parts of Gujarat have seismic hazard of magnitude 6 or less. Kachchh region is considered seismically one of the most active intraplate regions of the World. It is known to have low seismicity but high hazard in view of occurrence of fewer smaller earthquakes of M????6 in a region having three devastating earthquakes that occurred during 1819 (M _w7.8), 1956 (M _w6.0) and 2001 (M _w7.7). The second in order of seismic status is Narmada rift zone that experienced a severely damaging 1970 Bharuch earthquake of M5.4 at its western end and M????6 earthquakes further east in 1927 (Son earthquake), 1938 (Satpura earthquake) and 1997 (Jabalpur earthquake). The Saurashtra Peninsula south of Kachchh has experienced seismicity of magnitude less than 6.     

Region-specific deterministic and probabilistic seismic hazard analysis of Kanpur city

A seismic hazard map of Kanpur city has been developed considering the region-specific seismotectonic parameters within a 500-km radius by deterministic and probabilistic approaches. The maximum probable earthquake magnitude (M _max) for each seismic source has been estimated by considering the regional rupture characteristics method and has been compared with the maximum magnitude observed \(\left ({M_{\max }^{\text {obs}}}\right )\), \(M_{\max }^{\text {obs}} +0.5\) and Kijko method. The best suitable ground motion prediction equations (GMPE) were selected from 27 applicable GMPEs based on the ‘efficacy test’. Furthermore, different weight factors were assigned to different M _max values and the selected GMPE to calculate the final hazard value. Peak ground acceleration and spectral acceleration at 0.2 and 1 s were estimated and mapped for worst-case scenario and 2 and 10% probability of exceedance for 50 years. Peak ground acceleration (PGA) showed a variation from 0.04 to 0.36 g for DSHA, from 0.02 to 0.32 g and 0.092 to 0.1525 g for 2 and 10% probability in 50 years, respectively. A normalised site-specific design spectrum has been developed considering three vulnerable sources based on deaggregation at the city center and the results are compared with the recent 2011 Sikkim and 2015 Nepal earthquakes, and the Indian seismic code IS 1893.     

Seismicity and source parameters of moderate earthquakes in Sikkim Himalaya

In this study, we accurately relocate 360 earthquakes in the Sikkim Himalaya through the application of the double-difference algorithm to 4?years of data accrued from a eleven-station broadband seismic network. The analysis brings out two major clusters of seismicity??one located in between the main central thrust (MCT) and the main boundary thrust (MBT) and the other in the northwest region of Sikkim that is site to the devastating Mw6.9 earthquake of September 18, 2011. Keeping in view the limitations imposed by the Nyquist frequency of our data (10?Hz), we select 9 moderate size earthquakes (5.3????Ml????4) for the estimation of source parameters. Analysis of shear wave spectra of these earthquakes yields seismic moments in the range of 7.95?×?10²¹ dyne-cm to 6.31?×?10²³ dyne-cm and corner frequencies in the range of 1.8?C6.25?Hz. Smaller seismic moments obtained in Sikkim when compared with the rest of the Himalaya vindicates the lower seismicity levels in the region. Interestingly, it is observed that most of the events having larger seismic moment occur between MBT and MCT lending credence to our observation that this is the most active portion of Sikkim Himalaya. The estimates of stress drop and source radius range from 48 to 389?bar and 0.225 to 0.781?km, respectively. Stress drops do not seem to correlate with the scalar seismic moments affirming the view that stress drop is independent over a wide moment range. While the continental collision scenario can be invoked as a reason to explain a predominance of low stress drops in the Himalayan region, those with relatively higher stress drops in Sikkim Himalaya could be attributed to their affinity with strike-slip source mechanisms. Least square regression of the scalar seismic moment (M ₀) and local magnitude (Ml) results in a relation LogM ₀?=?(1.56?±?0.05)Ml?+?(8.55?±?0.12) while that between moment magnitude (M _w ) and local magnitude as M _w ?=?(0.92?±?0.04)Ml?+?(0.14?±?0.06). These relations could serve as useful inputs for the assessment of earthquake hazard in this seismically active region of Himalaya.     

Strong ground motion attenuation in Aswan area, Egypt

Recent and paleo seismicity indicate that moderate seismic activity is relatively large for Aswan area. This is a warning on the possibility of occurrence of earthquakes in the future too. No strong motion records are available in Aswan area for engineers to rely upon. Consequently, the seismological modeling is an alternative approach till sufficient instrumental records around Aswan become available. In the present study, we have developed new ground motion attenuation relationship for events spanning 4.0?≤? M _w?≤?7.0 and distance to the surface projection of the fault up to 100 km for Aswan based on a statistically simulated seismological model. We generated suites of ground motion time histories using stochastic technique. The ground motion attenuation relation describes the dependence of the strength of the ground motions on the earthquake magnitude and distance from the earthquake. The proposed equation for peak ground acceleration (PGA) for the bed rock is in the form of: $ {\mathbf{log}}{\text{ }}\left( {{\mathbf{PGA}}/{\mathbf{gal}}} \right){\text{ }} = {\mathbf{1}}.{\mathbf{24}} + {\mathbf{0}}.{\mathbf{358}}{M_{\mathbf{w}}} - {\text{ }}{\mathbf{log}}\left( {\mathbf{R}} \right){\text{ }}-{\text{ }}{\mathbf{0}}.{\mathbf{008}}{\text{ }}{\mathbf{R}}{\text{ }} + {\text{ }}{\mathbf{0}}.{\mathbf{22}}{\text{ }}{\mathbf{P}} $ . Where PGA is the peak ground acceleration in gal (cm/s²); M_w, its moment magnitude; R is the closest distance between the rupture projection and the site of interest; and the factor P is a dummy variable. It is observed that attenuation of strong motion in Aswan is correlated with those used before in Egypt.     

Numerical evaluation of seismic response of shallow foundation on loose silt and silty sand

This study includes the results of a set of numerical simulations carried out for sands containing plastic/non-plastic fines, and silts with relative densities of approximately 30?40% under different surcharges on the shallow foundation using FLAC 2D. Each model was subjected to three ground motion events, obtained by scaling the amplitude of the El Centro (1940), Kobe (1995) and Kocaeli (1999) Q12earthquakes. Dynamic behaviour of loose deposits underlying shallow foundations is evaluated through fully coupled nonlinear effective stress dynamic analyses. Effects of nonlinear soil structure interaction (SSI) were also considered by using interface elements. This parametric study evaluates the effects of soil type, structure weight, liquefiable soil layer thickness, event parameters (e.g., moment magnitude of earthquake (M _w ), peak ground acceleration PGA, PGV/PGA ratio and the duration of strong motion (D _5?95) and their interactions on the seismic responses. Investigation on the effects of these parameters and their complex interactions can be a valuable tool to gain new insights for improved seismic design and construction.     

Sensitivity analysis of seismic hazard for the northwestern portion of the state of Gujarat, India

We test the sensitivity of seismic hazard to three fault source models for the northwestern portion of Gujarat, India. The models incorporate different characteristic earthquake magnitudes on three faults with individual recurrence intervals of either 800 or 1600 years. These recurrence intervals imply that large earthquakes occur on one of these faults every 266–533 years, similar to the rate of historic large earthquakes in this region during the past two centuries and for earthquakes in intraplate environments like the New Madrid region in the central United States. If one assumes a recurrence interval of 800 years for large earthquakes on each of three local faults, the peak ground accelerations (PGA; horizontal) and 1-Hz spectral acceleration ground motions (5% damping) are greater than 1 g over a broad region for a 2% probability of exceedance in 50 years' hazard level. These probabilistic PGAs at this hazard level are similar to median deterministic ground motions. The PGAs for 10% in 50 years' hazard level are considerably lower, generally ranging between 0.2 g and 0.7 g across northwestern Gujarat. Ground motions calculated from our models that consider fault interevent times of 800 years are considerably higher than other published models even though they imply similar recurrence intervals. These higher ground motions are mainly caused by the application of intraplate attenuation relations, which account for less severe attenuation of seismic waves when compared to the crustal interplate relations used in these previous studies. For sites in Bhuj and Ahmedabad, magnitude (M) 7 3/4 earthquakes contribute most to the PGA and the 0.2- and 1-s spectral acceleration ground motion maps at the two considered hazard levels.     

Source parameters of some recent earthquakes in the Gulf of Aqaba, Egypt

Gulf of Aqaba is recognized as an active seismic zone where many destructive earthquakes have occurred. The estimation of source parameters and coda Q attenuation are the main target of this work. Fifty digital seismic events in eight short-period seismic stations with magnitude 2.5–5.2 are used. Most of these events occurred at hypocentral depths in the range of 7–20 km, indicating that the activity was restricted in the upper crust. Seismic moment, M _o, source radius, r, and stress drop, Δσ, are estimated from P- and S-wave spectra using the Brune’s seismic source model. The average seismic moment generated by the whole sequence of events was estimated to be 4.6E?+?22 dyne/cm. The earthquakes with higher stress drop occur at 10-km depth. The scaling relation between the seismic moment and the stress drop indicates a tendency of increasing seismic moment with stress drop. The seismic moment increases with increasing the source radius. Coda waves are sensitive to changes in the subsurface due to the wide scattering effects generating these waves. Single scattering model of local earthquakes is used to the coda Q calculation. The coda with lapse times 10, 20, and 30 s at six central frequencies 1.5, 3, 6, 12, 18, 24 Hz are calculated. The Q _c values are frequency dependent in the range 1–25 Hz, and are approximated by a least squares fit to the power law [ $ {Q_c}(f) = {Q_o}{(f/{f_o})^\eta } $ ]. The average of Q _c values increases from 53?±?10 at 1.5 Hz to 700?±?120 at 24 Hz. The average of Q _o values ranges from 13?±?1 at 1.5 Hz to 39?±?4 at 24 Hz. The frequency exponent parameter η ranges between 1.3?±?0.008 and 0.9?±?0.001.     

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).     

Magnitude conversion problem for the Turkish earthquake data

Earthquake catalogues which form the main input in seismic hazard analysis generally report earthquake magnitudes in different scales. Magnitudes reported in different scales have to be converted to a common scale while compiling a seismic data base to be utilized in seismic hazard analysis. This study aims at developing empirical relationships to convert earthquake magnitudes reported in different scales, namely, surface wave magnitude, M _S, local magnitude, M _L, body wave magnitude, m _b and duration magnitude, M _d, to the moment magnitude (M _w). For this purpose, an earthquake data catalogue is compiled from domestic and international data bases for the earthquakes occurred in Turkey. The earthquake reporting differences of various data sources are assessed. Conversion relationships are established between the same earthquake magnitude scale of different data sources and different earthquake magnitude scales. Appropriate statistical methods are employed iteratively, considering the random errors both in the independent and dependent variables. The results are found to be sensitive to the choice of the analysis methods.     

Seismic input at the archaeological site of Kancheepuram in Southern India

A probabilistic seismic hazard assessment at Kancheepuram in Southern India was carried out with the scope of defining the seismic input for the vulnerability assessment of historical and monumental structures at the site, in terms of horizontal Uniform Hazard Spectra and a suite of spectrum-compatible natural accelerograms to perform time-history analysis. The standard Cornell?CMcGuire and a zone-free approach have been used for hazard computations after the compilation of a composite earthquake catalogue for Kancheepuram. Epistemic uncertainty in the seismic hazard was addressed within a logic-tree framework. Deaggregation of the seismic hazard for the peak ground acceleration shows low seismicity at Kancheepuram controlled by weak-to-moderate earthquakes with sources located at short distances from the archaeological site. Suites of natural accelerograms recorded on rock have been selected by imposing a custom-defined compatibility criterion with the probabilistic spectra. The site of Kancheepuram is characterized by a seismicity controlled by weak-to-moderate earthquakes with sources at short distances from the site, the PGA expected for 475- and 2,475-year return period are, respectively, 0.075 and 0.132?g. The Indian code-defined spectra (DBE and MCE) tend to underestimate spectral ordinates at low periods. On the other hand, the PGA are comparable and the spectral ordinates for longer periods from the probabilistic study are significantly lower.     

Seismic hazard across Bulgaria and neighbouring areas: regional and site-specific maximum credible magnitudes and earthquake perceptibility

A probabilistic seismic hazard assessment is developed here using maximum credible earthquake magnitude statistics and earthquake perceptibility hazard. Earthquake perceptibility hazard is defined as the probability a site perceives ground shaking equal to or greater than a selected ground motion level X, resulting from an earthquake of magnitude M, and develops estimates for the most perceptible earthquake magnitude, M _P(max). Realistic and usable maximum magnitude statistics are obtained from both whole process and part process statistical recurrence models. These approaches are extended to develop relationships between perceptible earthquake magnitude hazard and maximum magnitude recurrence models that are governed by asymptotic and finite return period properties, respectively. Integrated perceptibility curves illustrating the probability of a specific level of perceptible ground motion due to all earthquakes over the magnitude range extending from ?∞ to a magnitude M _i are then developed from reviewing site-specific magnitude perceptibility. These lead on to achieving site-specific annual probability of exceedance hazard curves for the example cities of Sofia and Thessaloniki for both horizontal ground acceleration and ground velocity. Both the maximum credible earthquake magnitude M ₃ and the most perceptible earthquake magnitude M _P(max) are of importance to the earthquake engineer when approaching anti-seismic building design. Both forms of hazard are illustrated using contoured hazard maps for the region bounded by 39°–45°N, 19°–29°E. Patterns are observed for these magnitude hazard estimates—especially M _P(max) specific to horizontal ground acceleration and horizontal ground velocity—and compared to inferred patterns of crustal deformation across the region. The full geographic region considered is estimated to be subject to a maximum credible earthquake magnitude M ₃—estimated using cumulative seismic moment release statistics—of 7.53 M _w, calculated from the full content of the adopted earthquake catalogue, while Bulgaria’s capital, Sofia, is estimated a comparable value of 7.36 M _w. Sofia is also forecast most perceptible earthquake magnitudes for the lowest levels considered for horizontal ground acceleration of M _PA(50) = 7.20 M _w and horizontal ground velocity of M _PV(5) = 7.23 M _w for a specimen focal depth of 15 km.