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.     

Probabilistic seismic hazard analysis at a strategic site in the Bay of Bengal

Probabilistic seismic hazard analysis (PSHA) along the route of an offshore pipeline for the transport of oil in the Bay of Bengal has been performed, in order to set up design parameters and identify possible geohazards. The complexity of geological and seismotectonic setting of the region where the pipeline is planned to be installed is the result of the interaction of the Indian, Eurasian and Burmese tectonic plates. In order to properly account for the intricate way by which these plates interact, a large area extending 450 km from the pipeline route has been considered for the compilation of a comprehensive earthquake catalogue, spanning the period 1663–2012 AD. Differently from earlier PSHA analyses conducted in the region based on assuming two-dimensional polygons as seismogenic provinces, this study adopted a seismotectonic source model which also includes for the first time a linear tectonic lineament representing the northward extension of the Sunda mega thrust, responsible for the large Sumatra–Andaman earthquake of 26 December 2004. Hazard computations have been performed over a grid of sites spaced 0.045° covering a rectangular area which contains the pipeline. Epistemic uncertainty in the hazard computations has been taken into account by a logic tree framework, incorporating different seismotectonic source models, maximum cut-off magnitude and ground-motion prediction equations. Horizontal median uniform hazard spectra and median uniform hazard spectra plus and minus one sigma on stiff ground have been calculated at the selected sites for different return periods. Peak ground acceleration with 10 % probability of exceedance in 50 years has been compared with values from previous hazard studies available for Bangladesh.     

Seismotectonic regionalization of the Red Sea area and its application to seismic risk analysis

A seismological evaluation of the Red Sea margin is presented in this contribution based on the concept of seismotectonic regionalization. The geology and the tectonic structure are critically reviewed to define regions of homogeneous seismicity in the study area, and available seismicity data are implemented to estimate the seismic parameters of the region. The results of the study are applied to evaluate the seismic hazard of an offshore platform site.     

First-order regional seismotectonic model for South Africa

A first-order seismotectonic model was created for South Africa. This was done using four logical steps: geoscientific data collection, characterisation, assimilation and zonation. Through the definition of subunits of concentrations of earthquake foci and large neotectonic and structural domains, seismotectonic structures, systems and domains were created. Relatively larger controls of seismicity exist between the Great Escarpment and the coast. In the south, this region is characterised by large aeromagnetic anomalies and large EW trending faults. In the west, it is characterised by the NW–SE trending Wegener stress anomaly, radial-trending dykes and earthquake clusters. In the east, it is characterised by a large neotectonic domain where several large historical earthquakes occurred. In the centre of South Africa, several clusters of earthquake activity are found, often related to mining activity. Further north, seismicity is related to both mining activity and neotectonic deformation. This work contributes to the development of a seismotectonic model for South Africa by (1) bringing together, digitally, several data sets in a common GIS platform (geology, geophysics, stress, seismicity, neotectonics, topography, crustal and mantle structure and anisotropy), (2) understanding the significance of data sets for seismotectonic zonation and limitations thereof and (3) obtaining a reasonable regional model for use in seismic hazard assessments.     

Probabilistic seismic hazard assessment in the northeastern part of Algeria

This work involves updating the evaluation of seismic hazard in Northeast Algeria by a probabilistic approach. This reassessment attempts to resolve inconsistencies between seismic zoning in regional building codes and is further motivated by the need to refine the input data that are used to evaluate seismic hazard scenarios. We adopted a seismotectonic model that accounts for differences in interpretations of regional seismicity. We then performed a probabilistic assessment of regional seismic hazard in Northeast Algeria. Based on a homogeneous earthquake catalog and geological and seismotectonic data gathered in the first part of the study, a seismotectonic zoning map was created and seven risk areas were identified. For each area, peak ground acceleration hazard maps were produced. Details of the calculations are provided, including hazard curves at periods of 0.1, 0.2, 0.33, 0.5, 1.0, and 2.0 s and uniform hazard spectra at urban locations in the area, including Sétif, Constantine, Kherrata, Bejaia, and Jijel.     

A seismic hazard map of India and adjacent areas

We have produced a probabilistic seismic hazard map showing peak ground accelerations in rock for India and neighboring areas having a 10% probability of being exceeded in 50 years. Seismogenic zones were identified on the basis of historical seismicity, seismotectonics and geology of the region. Procedures for reducing the incompleteness of earthquake catalogs were followed before estimating recurrence parameters. An eastern United States acceleration attenuation relationship was employed after it was found that intensity attenuation for the Indian region and the eastern United States was similar. The largest probabilistic accelerations are obtained in the seismotectonic belts of Kirthar, Hindukush, Himalaya, Arakan-Yoma, and the Shillong massif where values of over 70% g have been calculated.     

Earthquake-induced liquefaction hazard mapping at national-scale in Australia using deep learning techniques

Australia is a relatively stable continental region but not tectonically inert, having geological conditions that are susceptible to liquefaction when subjected to earthquake ground motion. Liquefaction hazard assessment for Australia was conducted because no Australian liquefaction maps that are based on modern AI techniques are currently available. In this study, several conditioning factors including Shear wave velocity (Vs30), clay content, soil water content, soil bulk density, soil thickness, soil pH, distance from river, slope and elevation were considered to estimate the liquefaction potential index (LPI). By considering the Probabilistic Seismic Hazard Assessment (PSHA) technique, peak ground acceleration (PGA) was derived for 50 yrs period (500 and 2500 yrs return period) in Australia. Firstly, liquefaction hazard index (LHI) (effects based on the size and depth of the liquefiable areas) was estimated by considering the LPI along with the 2% and 10% exceedance probability of earthquake hazard. Secondly, ground acceleration data from the Geoscience Australia projecting 2% and 10% exceedance rate of PGA for 50 yrs were used in this study to produce earthquake induced soil liquefaction hazard maps. Thirdly, deep neural networks (DNNs) were also exerted to estimate liquefaction hazard that can be reported as liquefaction hazard base maps for Australia with an accuracy of 94% and 93%, respectively. As per the results, very-high liquefaction hazard can be observed in Western and Southern Australia including some parts of Victoria. This research is the first ever country-scale study to be considered for soil liquefaction hazard in Australia using geospatial information in association with PSHA and deep learning techniques. This study used an earthquake design magnitude threshold of M_w 6 using the source model characterization. The resulting maps present the earthquake-triggered liquefaction hazard and are intending to establish a conceptual structure to guide more detailed investigations as may be required in the future. The limitations of deep learning models are complex and require huge data, knowledge on topology, parameters, and training method whereas PSHA follows few assumptions. The advantages deal with the reusability of model codes and its transferability to other similar study areas. This research aims to support stakeholders’ on decision making for infrastructure investment, emergency planning and prioritisation of post-earthquake reconstruction projects.     

Seismological hazard assessment for a site in Northern Germany, an area of low seismicity

Potential sites of nuclear waste deposits in the Federal Republic of Germany are situated in areas of low seismicity. Nevertheless, seismic hazard assessment has to be performed for a very long time period in order to prove the facilities of the repositories able to withstand seismic induced loads; even though there is a considerable debate whether or not it is possible to quantify the seismic risk in such an area.

A combination of deterministic and probabilistic methods is used to assess the seismic hazard for a site in Northern Germany, fulfilling the standards of the German building code for nuclear power plants. As an example, the site of the former iron ore mine Konrad is investigated. The deterministic method is based on the assumption that the strongest earthquakes inside a tectonic region can happen everywhere there, also near the selected site. For the probabilistic method, several models describing the seismicity in an area of 200 km around the site are used to show the influence of the variability in input parameters, like the maximum intensity of each source region on the exceeding probability of the site intensity. It can be shown that the seismic hazard of a site in an area of low seismicity is mainly caused by the effects of distant but strong source regions and the background seismicity for very low probabilities.

Probabilistic evaluation has the advantage of quantifying the seismic risk. But deterministic and probabilistic methods together seem a practical tool for mutual control of the results and to overcome the weakness of each approach alone. The historical German earthquake catalog with an observation period of about 1200 years is the basis for the input data for a probabilistic model. From a deep knowledge of geological development and structural geology, the time history of the surrounding faults is developed. Indications were found that the nearest and most important fault was active at least 5 Ma ago. The combination of both seismicity and tectonics provides the basis for a long term prognostic with probabilities of exceedance in the order of 10⁻⁵ per year.

For the investigated site the following parameters were derived: site intensity as a function of exceeding probability; site acceleration; strong motion duration; site dependent response spectra for the surface and the underground inside the mine.     

Earthquake hazard assessment from the recent space and surface measurements of the Earth's gravity field

The distinctive nature of the structure of the Earth's gravity field over high seismicity areas, as observed in some limited scale studies, indicates that it is possible to use these associative patterns to outline areas of high seismicity or high earthquake hazard potential from a knowledge of the finer structure of the Earth's gravity field. The global scale investigation of such relationships as a function of the geophysical characteristics of various tectonic provinces, and the parameterisation of such relationships by hypocentral depth and earthquake intensity/energy data, become even more attractive because of the availability of (1) satellite-determined gravity models which provide global information on the long wavelength components of the gravity field, (2) satellite altimetry data which provide oceanwide information on the detailed geoidal structure, (3) surface gravity data which provide information on the short wavelength components in the areas of surface gravity coverage, and (4) bathymetric and topographic data which, though still somewhat limited in spite of their recent extensions of coverage, are available in more and more areas and provide information on the tectonic and morphological environments of an area to enable its gravity data conversion to some standard environments for direct comparisons of underlying structures. Additionally, the rapidly mounting geological and geophysical evidence of considerable intraplate tectonic activity, not so fashionable until recently, makes the study of these correlative patterns even more attractive and productive from a scientific viewpoint. This paper presents the major elements of the theoretical formulation for conducting such investigations.     

Tectonic stress,seismicity, and seismic hazard in the southeastern Carpathians

Intermediate-depth earthquakes in the Vrancea region occur in response to stress generation due to descending lithosphere beneath the southeastern Carpathians. In this article, tectonic stress and seismicity are analyzed in the region on the basis of a vast body of observations. We show a correlation between the location of intermediate-depth earthquakes and the predicted localization of maximum shear stress in the lithosphere. A probabilistic seismic hazard assessment (PSHA) for the region is presented in terms of various ground motion parameters on the utilization of Fourier amplitude spectra used in engineering practice and risk assessment (peak ground acceleration, response spectra amplitude, and seismic intensity). We review the PSHA carried out in the region, and present new PSHA results for the eastern and southern parts of Romania. Our seismic hazard assessment is based on the information about the features of earthquake ground motion excitation, seismic wave propagation (attenuation), and site effect in the region. Spectral models and characteristics of site-response on earthquake ground motions are obtained from the regional ground motion data including several hundred records of small and large earthquakes. Results of the probabilistic seismic hazard assessment are consistent with the features of observed earthquake effects in the southeastern Carpathians and show that geological factors play an important part in the distribution of the earthquake ground motion parameters.     

A seismic source zone model for the seismic hazard assessment of the Italian territory

We designed a new seismic source model for Italy to be used as an input for country-wide probabilistic seismic hazard assessment (PSHA) in the frame of the compilation of a new national reference map.

We started off by reviewing existing models available for Italy and for other European countries, then discussed the main open issues in the current practice of seismogenic zoning.

The new model, termed ZS9, is largely based on data collected in the past 10 years, including historical earthquakes and instrumental seismicity, active faults and their seismogenic potential, and seismotectonic evidence from recent earthquakes. This information allowed us to propose new interpretations for poorly understood areas where the new data are in conflict with assumptions made in designing the previous and widely used model ZS4.

ZS9 is made out of 36 zones where earthquakes with M_w > = 5 are expected. It also assumes that earthquakes with M_w up to 5 may occur anywhere outside the seismogenic zones, although the associated probability is rather low. Special care was taken to ensure that each zone sampled a large enough number of earthquakes so that we could compute reliable earthquake production rates.

Although it was drawn following criteria that are standard practice in PSHA, ZS9 is also innovative in that every zone is characterised also by its mean seismogenic depth (the depth of the crustal volume that will presumably release future earthquakes) and predominant focal mechanism (their most likely rupture mechanism). These properties were determined using instrumental data, and only in a limited number of cases we resorted to geologic constraints and expert judgment to cope with lack of data or conflicting indications. These attributes allow ZS9 to be used with more accurate regionalized depth-dependent attenuation relations, and are ultimately expected to increase significantly the reliability of seismic hazard estimates.     

A new probabilistic seismic hazard study of Spain

In this paper, we present a probabilistic seismic hazard analysis (PSHA) for mainland Spain that takes into account recent new results in seismicity, seismic zoning, and strong ground attenuation not considered in the latest PSHA of the Spanish Building Code. Those new input data have been obtained as a three-step project carried out in order to improve the existing hazard map for mainland Spain. We have produced a new earthquake catalogue for the area, in which the earthquakes are given in moment magnitude through specific deduced relationships for our territory based on intensity data (Mezcua et al. in Seismol Res Lett 75:75–81, 2004). In addition, we included a new seismogenetic zoning based on the recent partial zoning studies performed by different authors. Finally, as we have developed a new strong ground motion model for the area García Blanco (2009), it was considered in the hazard calculation together with other attenuations gathered from different authors using data compatible with our region. With this new data, a logic tree process is defined to quantify the epistemic uncertainty related to those parts of the process. A sensitivity test has been included in order to analyze the different models of ground motion and seismotectonic zonation used in this work. Finally, after applying a weighting scheme, a mean hazard map for PGA, based on rock type condition for 10% exceedance probability in 50 years, is presented, including 15th and 85th percentile hazard maps. The main differences with the present official building code hazard map are analyzed.     

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.     

Source-to-site distance distributions for area type of seismic sources used in PSHA applications

The use of recent ground motion prediction equations in probabilistic seismic hazard analysis (PSHA) with area type of seismic sources requires defining the probability distributions of various source-to-site distance metrics with finite fault rupture taken into account. This task is rendered very difficult due to large epistemic uncertainties involved in specifying the details of the causative faults for area sources of diffused seismicity. However, it may generally be possible to constrain the strike and dip angles for fault ruptures in area sources from regional seismotectonic and geological information. This paper proposes to estimate the various finite fault distance measures from a site to a location in an area source by averaging the distances for several fault rupture scenarios with randomly distributed strike and dip over specified ranges. To consider the spatial distribution of the seismicity, the paper then provides the guidelines for defining the distance distributions by assigning suitable weight factors to the distance estimates for a grid of locations in the source area. The PSHA computation based on the distance distributions thus defined is shown to provide quite realistic and objective estimate of the hazard.     

近断层强地震动预测中的有限断层震源模型

提出了近断层强地震动预测中建立活断层上设定地震有限断层震源模型的方法和步骤.首先,根据地震地质和地震活动性调查以及地球物理勘探等资料,确定活断层的空间方位和滑动类型; 然后,根据地震定标律确定活断层的宏观震源参数; 第三,将高强体模型与k平方滑动模型相结合,产生断层破裂面上的混合滑动分布.在此基础上,预测了与1994年Northridge地震断层类型、矩震级(M_w6.7)基本一致的设定地震的有限断层震源模型.最后,将预测的有限断层震源模型与基于地震学的、使用动力学拐角频率的地震动随机合成方法相结合,预测了1994年Northridge地震近断层12个基岩台站的加速度时程,并和实际记录进行了对比.结果表明,用上述方法和步骤建立的有限断层震源模型是可行、实用的.      

Systematization of active faults for the assessment of the seismic hazard

A systematization of active faults has been developed based on the progress of scientists from the leading countries in the world in the study of seismotectonics and seismic hazard problems. It is underlain by the concept of the fault-block structure of the geological-geophysical environment governed by the interaction of differently oriented active faults, which are divided into two groups—seismogenic and nonseismogenic faults. In seismogenic fault zones, the tectonic stress accumulated is relieved by means of strong earthquakes. Nonseismogenic fault zones are characterized by creep displacement or short-term, oscillatory, and reciprocal movements, which are referred to local superintense deformations of the Earth’s crust (according to the terminology used by Yu.O. Kuz’min). For a situation when a strong earthquake happens, a subgroup of seismodistributing faults has been identified that surround the seismic source and affect the distribution of the seismic waves and, as a consequence, the pattern of the propagation of the coseismic deformations in the fault-block environment. Seismodistributing faults are divided into transit and sealing faults. Along transit faults, secondary coseismic effects (landfalls, landslides, ground fractures, liquefaction, etc) are intensified during earthquakes. In the case of sealing faults, enhancement of the coseismic effects can be observed on the disjunctive limb nearest to the epicenter, whereas, on the opposite limb, the intensity of such effects appreciably decreases. Seismogenic faults or their systems are associated with zones of earthquake source origination (ESO), which include concentrated seismicity regions. In such zones, each earthquake source is related to the evolution of a fault system. ESO zones also contain individual seismogenic sources being focuses of strong earthquakes with M of ≥5.5 in the form of ruptures, which can be graphically represented in 2D or 3D as a surface projection of the source. Depending on the type of data based on which they are identified, individual seismogenic sources are divided into geological-geophysical and macroseismic sources. The systematization presented is the theoretical basis for and the concept of the relational database that is being developed by the authors as an information system for the generation of seismotectonic GIS projects required for the subsequent analysis of the seismic hazard and the assessment of the probability of the origination of macroseismic earthquake effects in a predetermined location.     

Seismic hazard assessment of Chennai city considering local site effects

Chennai city suffered moderate tremors during the 2001 Bhuj and Pondicherry earthquakes and the 2004 Sumatra earthquake. After the Bhuj earthquake, Indian Standard IS: 1893 was revised and Chennai city was upgraded from zone II to zone III which leads to a substantial increase of the design ground motion parameters. Therefore, a comprehensive study is carried out to assess the seismic hazard of Chennai city based on a deterministic approach. The seismicity and seismotectonic details within a 100 km radius of the study area have been considered. The one-dimensional ground response analysis was carried out for 38 representative sites by the equivalent linear method using the SHAKE91 program to estimate the ground motion parameters considering the local site effects. The shear wave velocity profile was inferred from the corrected blow counts and it was verified with the Multichannel Analysis of Surface Wave (MASW) test performed for a representative site. The seismic hazard is represented in terms of characteristic site period and Spectral Acceleration Ratio (SAR) contours for the entire city. It is found that structures with low natural period undergo significant amplification mostly in the central and southern parts of Chennai city due to the presence of deep soil sites with clayey or sandy deposits and the remaining parts undergo marginal amplification.     

Improving the level of seismic hazard parameters in Saudi Arabia using earthquake location

Saudi Arabia is characterized as largely aseismic; however, the tectonic plate boundaries that surround it are very active. To improve characterization of seismicity and ground motion hazard, the Saudi Arabian Digital Seismic Network (SANDSN) was installed in 1998 and continues to be operated by the Saudi Geological Survey (SGS) and King Abdulaziz City for Science and Technology (KACST). This article describes research performed to improve seismic hazard parameters using earthquake location and magnitude calibration of the high-quality SANDSN data. The SANDSN consists of 38 seismic stations, 27 broadband, and 11 short period. All data are telemetered in real time to a central facility at KACST in Riyadh. The SANDSN stations show low background noise levels and have good signal detection capabilities; however, some stations show cultural noise at frequencies above 1.0 Hz. We assessed the SANDSN event location capabilities by comparing KACST locations with well-determined locations derived from ground truth or global observations. While a clear location bias exists when using the global average iasp91 earth model, the locations can be improved by using regional models optimized for different tectonic source regions. The article presents detailed analysis of some events and Dead Sea explosions where we found gross errors in estimated locations. New velocity models we calculated that should improve estimated locations of regional events in three specific regions include (1) Gulf of Aqabah—Dead Sea region, (2) Arabian Shield, and (3) Arabian Platform. Recently, these models were applied to the SANDSN to improve local and teleseismic event locations and to develop an accurate magnitude scale for Saudi Arabia. The Zagros Thrust presents the most seismic hazard to eastern Saudi Arabia because of the frequent occurrence of earthquakes. Although these events are 200 km or further from the Arabian coast, wave propagation through sedimentary structure of the Gulf causes long-duration ground motions for periods between 3 and 10 s. Such ground motions could excite response in large engineered structures (e.g., tall buildings and long bridges) such as was experienced after the November 22, 2005 Qeshm Island earthquake off the southern coast of Iran.