Problems in the application of the SSHAC probability method for assessing earthquake hazards at Swiss nuclear power plants

From 2000 to 2004 a large scale probabilistic seismic hazard analysis (PEGASOS) was created and performed as a research project, sponsored by the Swiss NPP utilities to improve the assessment methodology for seismic risk in support of the plant-specific seismic PRAs. The project followed the methodology of the SSHAC [Senior Seismic Hazard Analysis Committee (SSHAC), 1997. Recommendations for Probabilistic Seismic Hazard Analysis: Guidance on Uncertainty and Use of Experts. NU-REG/CR-6372] procedures at its most elaborate way—level 4. Before practical implementation was to occur, a detailed review was performed including validation tests and analysis of uncertainty propagation. This paper presents the main results of the review. The review revealed that current PSHA (Probabilistic Seismic Hazard Analysis) methodology as based on logic trees, in conjunction with the SSHAC procedures, potentially leads to a significant overestimation of the seismic hazard in areas with low seismic activity due to the inherent possibilities of unconstrained accumulation of uncertainties. The preliminary results of the project did not pass any of our logical geological–scientific benchmark tests used in our attempts to perform a validation process of the PEGASOS analysis method. Some of the problems encountered are of generic nature and shall be studied carefully before making the decision of whether or not the Swiss nuclear power industry will adopt the recommended use of SSHAC-procedures as a basis for the evaluation of the seismic hazard for individual nuclear power plant seismic PRA without correction.     

白龙江引水工程代古寺水库及其周围地区地震危险性评估

白龙江引水工程是我国拟建的一项重大战略工程,而代古寺水库是该工程的水源枢纽。代古寺水库及其周围地区(本文研究区)活动断层发育、大地震频发,故亟需开展可靠的地震危险性评估,为该研究区内的工程建设和运营保驾护航。由于传统评估方法物理依据不足,难以正确评估研究区的地震危险性,故本文采用了基于地震物理预测的地震危险性评估新方法。研究结果表明,该研究区位于海原地震区,未来100年内该研究区的地震危险性主要源于海原地震区的下一次M_S8.5标志性地震。根据断层地震活动、发震潜力与展布特征,我们预判了该标志性地震的可能发震断层和震中位置;应用地震烈度衰减关系,考虑不同震中位置,分别计算了其产生的地震烈度。为确保“百年大计”的白龙江引水工程代古寺水库水资源枢纽安全,我们建议该研究区的抗震设防烈度不宜低于8度。     

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 procedure to estimate the seismic hazard in an urban area: an application to Acireale (Eastern Sicily)

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

A database for worldwide seismicity quantification

The first step in a seismicity analysis usually consists of defining the seismogenic units, seismic zones or individual faults. The worldwide delimitation of these zones involves an enormous effort and is often rather subjective. Also, a complete recording of faults will not be available for a long time yet. The seismicity model presented in this paper therefore is not based on individually defined seismic zones but rather on the assumption that each point in a global 1/2° grid of coordinates represents a potential earthquake source. The corresponding seismogenic parameters are allocated to each of these points. The earthquake occurrence frequency, one of the most important parameters, is determined purely statistically by appropriately spreading out the positions of past occurrences. All the other significant seismicity characteristics, such as magnitude-frequency relations, maximum possible magnitude and attenuation laws including the dependence on focal depth are determined in a global 1/2° grid of co-ordinates. This method of interpreting seismicity data allows us to establish a transparent, sufficiently precise representation of seismic hazard which is ideally suited for computer-aided risk analyses.     

孕震构造块体与相应地震区划分方法

可靠地划分地震区可奠定地震预测与地震危险性评价的地质基础,具有十分重要的意义。笔者等通过研究分析指出板内孕震构造块体侧向边界可由区域性大断层或由区域性大断层与板块边界界定,底边界为康拉德面或低速高导层;板间孕震构造块体为俯冲板块,可由区域性大断层和(或)板块边界约束;在同一个孕震构造块体和同一轮地震周期的地震具有内在联系。因此,地震区可定义为代表相应孕震构造块体地震活动的区域,其可表征该块体内源自锁固段破裂的地震活动。基于笔者等提出的孕震构造块体和相应地震区边界确定原则,把全球两大地震带(环太平洋地震带和欧亚地震带)划分为62个地震区;每个地震区的分区方案均通过了多锁固段脆性破裂理论的检验,这说明方案可靠。进而,笔者等归纳总结了地震区划分方法。     

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.     

Assessment of the industrial utilization of scoria materials in central Harrat Rahat, Saudi Arabia

Although basaltic volcanic scoria deposits are widespread in and around the scoria cones present in the Harrats, very few occurrences have been evaluated and still fewer have been exploited. The scoria from a quarry in central Harrat Rahat was investigated and assessed for its industrial utilization. Cubes were prepared from several concrete mixes using scoria as lightweight coarse and fine aggregates in different percentages. The compressive strength values of the cubes were found to be acceptable and satisfy the ASTM (1995) requirement for structural concrete. The scoria was also assessed for its utilization as a cement additive. Pozzolanic activity was tested according to the Italian standards and found to be acceptable. The strength activity index with Portland cement and the effectiveness of scoria admixture in controlling alkali-silica reactions were tested according to ASTM (1995) standards. Mortar cubes were specially prepared for these studies using different mixes and different storage procedures. The results satisfied the ASTM (1995) requirements as cement additive. The utilization of scoria as a heat-insulating material was tested and the results were found to satisfy the ASTM (1995) requirements. This fact suggests it could be utilized in the manufacture of the building blocks. It is recommended to evaluate the other scoria deposits, exploit the economically feasible ones and utilize them for different industrial applications. The manufacturing of heat-insulating concrete or building blocks using scoria is of prime importance as an energy saver.     

Cautious neotectonic hypotheses for assessing the seismic hazard in northeastern Italy

A simplified tectonic scheme for hazard purposes was recently adopted for northeastern Italy, introducing large generalized seismogenic areas containing systems of complex geometry faults. This scheme considers only major faults with documented seismic activity. In the present analysis, a different tectonic scheme, with linear elements as seismogenic sources, is presented. The assessment of the regional seismic hazard is done with the fault rupture model, its most important advantage being the recognition that the length of fault rupture during an earthquake is an important consideration in probabilistic calculations of seismic hazard. Moreover, some structures with no associated seismicity but with notable neotectonic activity are considered, and their contribution to the results investigated. Important uncertainties such as those in the maximum possible magnitude of future earthquakes, in the location of the fault, in the focal depth, and in the attenuation law are accounted for in the calculations and their influence studied. The results identify a seismic belt running from Lake Garda to Friuli and along the Yugoslav coast and are very similar to those already known for Friuli, with the largest values corresponding to the zone around Gemona. Some slight differences in the shape of the areas of equal acceleration are probably due to the delineation of the seismic sources of the proposed model. For a cautious elaboration, some neotectonic lines without present seismicity were added into the fault model. Their contribution is negligible in the areas of highest acceleration, but increases remarkably in the areas where acceleration is not expected to exceed the medium values.     

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.     

The role of small‐scale fold and fault development in seismogenic zones: example of the western Huércal‐Overa Basin (eastern Betic Cordillera,Spain)

The NW–SE shortening between the African and the Eurasian plates is accommodated in the eastern Betic Cordillera along a broad area that includes large N‐vergent folds and kilometric NE–SW sinistral faults with related seismicity. We have selected the best exposed small‐scale tectonic structures located in the western Huércal‐Overa Basin (Betic Cordillera) to discuss the seismotectonic implications of such structures usually developed in seismogenic zones. Subvertical ESE–WNW pure dextral faults and E–W to ENE–ESW dextral‐reverse faults and folds deform the Quaternary sediments. The La Molata structure is the most impressive example, including dextral ESE–WNW Neogene faults, active southward‐dipping reverse faults and associated ENE–WSW folds. A molar M1 assigned to Mimomys savini allows for precise dating of the folded sediments (0.95–0.83 Ma). Strain rates calculated across this structure give ～0.006 mm a^?1 horizontal shortening from the Middle Pleistocene up until now. The widespread active deformations on small‐scale structures contribute to elastic energy dissipation around the large seismogenic zones of the eastern Betics, decreasing the seismic hazard of major fault zones. Copyright © 2009 John Wiley & Sons, Ltd.     

基于震害统计的城市地下空间地震安全性评价

建筑物的地震安全性是城市规划和建设过程首先要回答的问题。我国城市地震安全性评价的方法理论多针对地上建筑物,而对地下空间的地震安全性研究较为薄弱,严重滞后于城市发展对地下空间的需求。活断层是诱发地震、导致建筑物破坏的的直接因素。考虑到空间关系上,地下空间与断层之间的交互关系为相交或相离。因此,本文将地下空间分为两类:与断层相交的地下空间称为跨断层地下空间,远离断层的地下空间称为远离断层地下空间。本文尝试将断裂带同震地表破裂、地震峰值加速度、地震烈度等地表地震安全性评价考量的要素与地下空间埋藏深度建立联系,并在此基础上总结基于震害统计的地下空间地震安全性评价方法。最后,本文选取地下空间利用需求较高的深圳和北京地区为实例进行介绍。     

青海玉树地区活动断裂与地震

青海玉树是巴颜喀拉地块西南边界上的典型历史强震区。最新的活动断裂遥感解译与地表调查结果表明,该区新构造期间主要发育清水河断裂带、玉树断裂带、阿布多断裂带和杂多断裂带4条NW向左旋走滑活动断裂带。其中,构成玉树—鲜水河—小江断裂系尾端构造的玉树活动断裂带是该区活动性最显著的岩石圈断裂。该断裂是由当江断裂、结古—结隆断裂和巴塘断裂3条斜接的主干断层和夹杂其间的多条次级断裂所共同构成的Z型左旋剪切张扭性变形带。它在上新世以来和晚第四纪期间的左旋走滑速率为4.0~5.4mm/a,调节了该区大部分的块体挤出与旋转变形,并构成该区大震活动的主要控震构造。历史强震梳理和古地震研究揭示,玉树主干走滑断裂带自约14530a BP以来至少发生了包括2010年地震在内的共11次大地震,原地重复间隔平均在千年以上,最长达近3000a。1738年玉树西北地震之后,玉树—甘孜断裂带的主干断层表现为平均间隔为50~100a的低频、串联式分段破裂过程,并且大震活动存在从东南向西北迁移的趋势。通过对玉树断裂未来大地震危险性进行综合地质判定认为,该区至少仍存在6段未来百年内大地震危险程度不同的地震空区,潜在的大地震震级为Mw6.6~7.3,其中危险性相对较高的段落主要是当江断裂带的当江—拉则段和结古—结隆断裂带上的结隆—叶卡诺段与桑卡—相古段。     

Paleoseismicity and seismic hazard in southern Patagonia (Argentina-Chile; 50°–55°S) and the role of the Magallanes-Fagnano transform fault

Patagonia, including the island of Tierra del Fuego, lies in southernmost South America at the junction of the South American, Antarctic, and Scotia tectonic plates. Historical and instrumental records have documented several local earthquakes of damaging magnitude, posing a threat to the rapidly growing population of 300,000 and the expanding industrial and service infrastructure. Short and inaccurate instrumental records of local seismic events and a diffuse epicenter distribution not clearly related to the recognized seismogenic structures have hindered an adequate evaluation of the seismic hazard for this region. To improve this situation, a paleoseismological study was carried out on two gravelly strandplains on the Atlantic coast of Patagonia. Surveying combined ground-probing radar, vertical electric sounding, and seismic refraction. Coseismic normal faults buried beneath the strandplain bodies were revealed and related to the morphology of the strandplains. The faults have probable ages between 0.9 and 6.4 kyr BP and a recurrence rate of about 1 kyr. The more likely source for these structures is the Magallanes-Fagnano fault, a continental transform fault that crosses Tierra del Fuego. The distance of more than 300 km from the buried coseismic structures to the trace of the Magallanes-Fagnano fault argues for high-magnitude earthquake activity on this fault throughout the Holocene. Urban development on soft glacial and alluvial substrates increases the hazard.     

Seismicity, seismotectonics and seismic hazard of Italy

A first generation of probabilistic seismic hazard maps of the Italian country are presented. They are based on seismogenic zoning deriving from a kinematic model of the structural tectonic units and on an earthquake catalogue with the foreshock and aftershock events filtered out. The following ground motion parameters have been investigated and mapped using attenuation equations based on strong-motion recordings of Italian earthquakes: peak ground acceleration and velocity; Arias intensity; strong motion duration; and the pseudovelocity and pseudoacceleration spectral values at 14 fixed frequencies both for the vertical and the largest horizontal component. A Poissonian model of earthquake occurrence is assumed as a default and the hazard maps are presented in terms of ground motion values expected to be exceeded at a 10% probability level in 50 years (return period 475 years) according to the requirement of Eurocode 8 for the seismic classification of national territories, as well as in terms of exceedance probabilities of selected ground motion values. Finally, as a tentative study, the use of hybrid methods (implementing both seismogenic zones and structures), renewal processes (including earthquake forecasting) and the influence of site effects (as the basis for the planning of earthquake scenarios) were explored.     

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.     

On the Bayesian analysis of the earthquake hazard in the North-East Indian peninsula

The Bayesian extreme-value distribution of earthquake occurrences has been used to estimate the seismic hazard in 12 seismogenic zones of the North-East Indian peninsula. The Bayesian approach has been used very efficiently to combine the prior information on seismicity obtained from geological data with historical observations in many seismogenic zones of the world. The basic parameters to obtain the prior estimate of seismicity are the seismic moment, slip rate, earthquake recurrence rate and magnitude. These estimates are then updated in terms of Bayes’ theorem and historical evaluations of seismicity associated with each zone. From the Bayesian analysis of extreme earthquake occurrences for North-East Indian peninsula, it is found that for T = 5 years, the probability of occurrences of magnitude (M _w = 5.0–5.5) is greater than 0.9 for all zones. For M _w = 6.0, four zones namely Z1 (Central Himalayas), Z5 (Indo-Burma border), Z7 (Burmese arc) and Z8 (Burma region) exhibit high probabilities. Lower probability is shown by some zones namely␣Z4, Z12, and rest of the zones Z2, Z3, Z6, Z9, Z10 and Z11 show moderate probabilities.     

3-D seismic structure of the Kachchh,Gujarat, and its implications for the earthquake hazard mitigation

Several pieces of studies on the January 26, 2001, Bhuj earthquake (Mw 7.6) revealed that the mainshock was triggered on the hidden unmapped fault in the western part of Indian stable continental region that caused a huge loss in the entire Kachchh rift basin of Gujarat, India. Occurrences of infrequent earthquakes of Mw 7.6 due to existence of hidden and unmapped faults on the surface have become one of the key issues for geoscientific research, which need to be addressed for evolving plausible earthquake hazard mitigation model. In this study, we have carried out a detailed autopsy of the 2001 Bhuj earthquake source zone by applying three-dimensional (3-D) local earthquake tomography (LET) method to a completely new data set consisting of 576 local earthquakes recorded between November 2006 and April 2009 by a seismic network consisting of 22 numbers of three-component broadband digital seismograph stations. In the present study, a total of 7560 arrival times of P-wave (3820) and S-wave (3740) recorded at least 4 seismograph stations were inverted to assimilate 3-D P-wave velocity (Vp), S-wave velocity (Vs), and Poisson’s ratio (σ) structures beneath the 2001 Bhuj earthquake source zone for reliable interpretation of the imaged anomalies and its bearing on earthquake hazard of the region. The source zone is located near the triple junction formed by juxtapositions of three Indian, Arabian, and Iranian tectonic plates that might have facilitated the process of brittle failure at a depth of 25 km beneath the KRB, Gujarat, which caused a gigantic loss to both property and persons of the region. There may be several hidden seismogenic faults around the epicentral zone of the 2001 Bhuj earthquake in the area, which are detectable using 3-D tomography to minimize earthquake hazard for a region. We infer that the use of detailed 3-D seismic tomography may offer potential information on hidden and unmapped faults beneath the plate interior to unravel the genesis of such big damaging earthquakes. This study may help in evolving a comprehensive earthquake risk mitigation model for regions of analogous geotectonic settings, elsewhere in the world.     

Seismogenic structures along continental convergent zones: from oblique subduction to mature collision

We summarize seismogenic structures in four regions of active convergence, each at a different stage of the collision process, with particular emphases on unusual, deep-seated seismogenic zones that were recently discovered. Along the eastern Hellenic arc near Crete, an additional seismogenic zone seems to occur below the seismogenic portion of the interplate thrust zone—a configuration found in several other oblique subduction zones that terminate laterally against collision belts. The unusual earthquakes show lateral compression, probably reflecting convergence between the subducting lithosphere's flank and the collision zone nearby. Along oblique zones of recent collision, the equivalence between space and time reveals the transition from subduction to full collision. In particular, intense seismicity beneath western Taiwan indicates that along the incipient zone of arc–continent collision, major earthquakes occur along high-angle reverse faults that reach deep into the crust or even the uppermost mantle. The seismogenic structures are likely to be reactivated normal faults on the passive continental margin of southeastern China. Since high-angle faults are ineffective in accommodating horizontal motion, it is not surprising that in the developed portion of the central Taiwan orogen (<5 Ma), seismogenic faulting occurs mainly along moderate-dipping (20–30°) thrusts. This is probably the only well-documented case of concurrent earthquake faulting on two major thrust faults, with the second seismogenic zone reaching down to depths of 30 km. Furthermore, the dual thrusts are out-of-sequence, being active in the hinterland of the deformation front. Along the mature Himalayan collision zone, where collision initiated about 50 Ma ago, current data are insufficient to distinguish whether most earthquakes occurred along multiple, out-of-sequence thrusts or along a major ramp thrust. Intriguingly, a very active seismic zone, including a large (M_w=6.7) earthquake in 1988, occurs at depths near 50 km beneath the foreland. Such a configuration may indicate the onset of a crustal nappe, involving the entire cratonic crust. In all cases of collision discussed here, the basal decollement, a key feature in the critical taper model of mountain building, appears to be aseismic. It seems that right at the onset of collision, earthquakes reflect reactivation of high-angle faults. For mature collision belts, earthquake faulting on moderate-dipping thrust accommodates a significant portion of convergence—a process involving the bulk of crust and possibly the uppermost mantle.     

滇西南2014年景谷中-强震群的地质构造成因——茶房-普文断裂带贯通过程的构造响应

2014年10—12月期间,云南景谷接连发生了Ms6.6、Ms5.8、Ms5.9三次中-强地震。为确定地震的地质构造成因,在地表调查的基础上,综合该区的地质构造情况、烈度与余震分布、震源机制解等资料,确定此次震群活动的宏观震中位于永平盆地东南侧山地,发震断层为地质与地貌表现不显著的NW向右旋走滑断层。此次震群活动及余震迁移过程指示,由于断层斜接部位岩桥的临时阻碍,Ms6.6地震破裂在向南东扩展过程中发生短暂停滞,突破障碍后进一步引发了Ms5.8和Ms5.9地震,这符合震源破裂沿NW向发震断裂分段破裂的行为。区域活动断裂的遥感解译结果发现,发震断层位置恰好处于NW向右旋走滑的茶房断裂与普文断裂之间,区域上属于该断裂带的不连贯部位,指示此次中-强震群活动应该是茶房-普文断裂带贯通过程的构造活动表现。结合思茅地块的历史地震资料发现,思茅地块地震活动多以小于等于6.8级为主,发震构造多为NW向断裂。指示在现今构造应力场作用下,该区NW向断裂的活动性相对NE向断裂更加显著,属于该区主要控震构造,应在今后的地震地质工作中给予更多关注。