Study of the epicentral trends and depth sections for aftershocks of the 26th january 2001, Bhuj earthquake in Western India

The Geological Survey of India (GSI) established a twelve-station temporary microearthquake (MEQ) network to monitor the aftershocks in the epicenter area of the Bhuj earthquake (M _w7.5) of 26th January 2001. The main shock occurred in the Kutch rift basin with the epicenter to the north of Bhachao village, at an estimated depth of 25 km (IMD). About 3000 aftershocks (M _d ≥ 1.0), were recorded by the GSI network over a monitoring period of about two and half months from 29th January 2001 to 15th April 2001. About 800 aftershocks (M _d ≥ 2.0) are located in this study. The epicenters are clustered in an area 60 km × 30 km, between 23.3‡N and 23.6‡N and 70‡E and 70.6‡E. The main shock epicenter is also located within this zone. Two major aftershock trends are observed; one in the NE direction and other in the NW direction. Out of these two trends, the NE trend was more pronounced with depth. The major NE-SW trend is parallel to the Anjar-Rapar lineament. The other trend along NW-SE is parallel to the Bhachao lineament. The aftershocks at a shallower depth (<10km) are aligned only along the NW-SE direction. The depth slice at 10 km to 20 km shows both the NE-SW trend and the NW-SE trend. At greater depth (20 km–38 km) the NE-SW trend becomes more predominant. This observation suggests that the major rupture of the main shock took place at a depth level more than 20 km; it propagated along the NE-SW direction, and a conjugate rupture followed the NW-SE direction. A N-S depth section of the aftershocks shows that some aftershocks are clustered at shallower depth ≤ 10 km, but intense activity is observed at 15–38 km depth. There is almost an aseismic layer at 10–15 km depth. The activity is sparse below 38 km. The estimated depth of the main shock at 25 km is consistent with the cluster of maximum number of the aftershocks at 20–38 km. A NW-SE depth section of the aftershocks, perpendicular to the major NE-SW trend, indicates a SE dipping plane and a NE-SW depth section across the NW-SE trend shows a SW dipping plane. The epicentral map of the stronger aftershocksM ≥ 4.0 shows a prominent NE trend. Stronger aftershocks have followed the major rupture trend of the main shock. The depth section of these stronger aftershocks reveals that it occurred in the depth range of 20 to 38 km, and corroborates with a south dipping seismogenic plane.     

Estimation of static stress changes after the 2001 Bhuj earthquake: Implications towards the northward spatial migration of the seismic activity in Kachchh,Gujarat

Spatial-temporal patterns of aftershocks of the 2001 Mw7.7 Bhuj earthquake during 2001–2008 reveal a northward spatial migration of seismic activity in the Kachchh seismic zone, which could be related with the loading stresses caused by the continued occurrences of aftershocks on the north Wagad fault (NWF), the causative fault of the 2001-mainshock. Aiming at explaining the observed northward migration of activity, we modelled the Coulomb failure stress change (DCFS) produced by the 2001-mainshock, the 2006 Mw5.6 Gedi fault (GF) and the 2007 Mw4.5 Allah bund fault (ABF) events on optimally oriented plane. A strong correlation between occurrences of earthquakes and regions of increased DCFS is obtained on the associated three faults i.e. NWF, ABF and GF. Predicted DCFS on the GF increased by 0.9 MPa at 3 km depth, where the 7^th March 2006 Mw5.6 event occurred, whereas predicted DCFS on the ABF increased by 0.07 MPa at 30 km depth, where the 15^th December 2007 Mw4.5 event occurred. Focal mechanism solutions of three events on the ABF have been estimated using the iterative inversion of broadband data from 5–10 stations, which are also constrained by the first P-motion data from 8–12 stations. These focal mechanism solutions for the ABF events reveal a dominant reverse movement with a strike-slip component along a preferred northwest or northeast dipping plane (∼50–70°). Focal mechanisms of the events on all the three fault zones reveal an N-S oriented P- axis or maximum principal stress in the region, which agrees with the prevailing N-S compression over the Indian plate. It is apparent that the northward migration of the static stress changes from the NWF, resulting from the occurrence 2001 Bhuj mainshock, might have caused the occurrence of the events on the GF and ABF during 2006–08.     

以精定位背景地震活动与震源机制解研究郯庐断裂带中南段现今活动习性

郯庐断裂带是中国东部板内一条规模最大的强构造变形带与地震活动断裂带,其断裂结构与历史地震活动性具明显的分段活动性。文中通过沿郯庐断裂带中南段的历史地震活动性、精定位背景地震活动性与震源机制解分析,讨论了断裂带的深部几何结构与现今活动习性。现今地震活动在中段主要沿1668年郯城MS 8地震破裂带线性分布,线性条带在泗洪-诸城间延伸约340 km长,为1668年地震长期缓慢衰减的余震序列活动。大震地表破裂遗迹与精定位地震分布都揭示出郯庐断裂带中段的两条全新世活动断裂昌邑-大店断裂与安丘-莒县断裂以右阶斜列的形式共同参与了1668年郯城MS 8地震破裂。精定位震源剖面刻画出两条断裂结构面呈高角度相背而倾,其中昌邑-大店断裂倾向SE,安丘-莒县断裂倾向NW,两条断裂在深部没有合并汇聚。余震活动所代表的1668年地震震源破裂带是郯庐断裂带中现今尚未闭锁的安全段落,对应于高b值段。而未发生破裂的安丘以北段,小震活动不活跃,b值低,现今可能已成为应力积累的闭锁段。震源机制解揭示的断裂应力状态在中段以NE向主压应力为主,表现为右旋走滑活动性质,且存在少量正断分量,南段转为以NEE至近EW向为主,存在少量的逆冲分量。在中段与南段的转折处宿迁-嘉山段,主压应力方向垂直断裂带走向呈NWW向,反映出局部以挤压为主的应力特征,其中泗洪-嘉山段也是历史地震未破裂段,现今小震活动不活跃,因此该段可能更易于应力积累。精定位小震活动在郯庐断裂与北西向断裂相交汇处聚集,反映出北西向断裂的新活动性,以及郯庐断裂带现今的逆冲作用。在断裂带南端,精定位背景地震活动沿与其相交汇的襄樊-广济断裂带东段呈北西向线性分布,表明了该段的现今活动性。     

What triggers Koyna region earthquakes? Preliminary results from seismic tomography digital array

The cause for prolific seismicity in the Koyna region is a geological enigma. Attempts have been made to link occurrence of these earthquakes with tectonic strain as well as the nearby reservoirs. With a view to providing reliable seismological database for studying the earth structure and the earthquake process in the Koyna region, a state of the art digital seismic network was deployed for twenty months during 1996–97. We present preliminary results from this experiment covering an area of 60 × 80 km² with twenty seismic stations. Hypocentral locations of more than 400 earthquakes confined to 11×25 km² reveal fragmentation in the seismicity pattern — a NE — SW segment has a dip towards NW at approximately 45°, whilst the other two segments show a near vertical trend. These seismic segments have a close linkage with the Western Ghat escarpment and the Warna fault. Ninety per cent of the seismicity is confined within the depth range of 3–10 km. The depth distribution of earthquakes delimits the seismogenic zone with its base at 10 km indicating a transition from an unstable to stable frictional sliding regime. The lack of shallow seismicity between 0 and 3 km indicates a mature fault system with well-developed gouge zones, which inhibit shallow earthquake nucleation. Local earthquake travel time inversion for P- and S-waves show ≈ 2% higher velocity in the seismogenic crust (0–10 km) beneath the epicentral tract relative to a lower velocity (2–3%) in the adjoining region. The high P- and S-wave velocity in the seismogenic crust argues against the presence of high pressure fluid zones and suggests its possible linkage with denser lithology. The zone of high velocity has been traced to deeper depths (≈ 70 km) through teleseismic tomography. The results reveal segmented and matured seismogenic fault systems in the Koyna region where seismicity is possibly controlled by strain build up due to competent lithology in the seismic zone with a deep crustal root.     

Seismic slip rates,potential subsurface rupture areas and seismic potential of the Vienna Basin Transfer Fault

The Vienna Basin Transfer Fault (VBTF) is a slow active fault with moderate seismicity (I _max~8–9, M _max~5.7) passing through the most vulnerable regions of Austria and Slovakia. We use different data to constrain the seismic potential of the VBTF including slip values computed from the seismic energy release during the 20th century, geological data on fault segmentation and a depth-extrapolated 3-D model of a generalized fault surface, which is used to define potential rupture zones. The seismic slip of the VBTF as a whole is in the range of 0.22–0.31 mm/year for a seismogenic fault thickness of 8 km. Seismic slip rates for individual segments vary from 0.00 to 0.77 mm/year. Comparing these data to geologically and GPS-derived slip velocities (>1 mm/year) proofs that the fault yields a significant seismic slip deficit. Segments of the fault with high seismic slip contrast from segments with no slip representing locked segments. Fault surfaces of segments within the seismogenic zone (4–14 km depth) vary from 55 to 400 km². Empirical scaling relations show that these segments are sufficiently large to explain both, earthquakes observed in the last centuries, and the 4th century Carnuntum earthquake, for which archeo-seismological data suggest a magnitude of M ≥ 6. Based on the combination of all data (incomplete earthquake catalog, seismic slip deficits, locked segments, potential rupture areas, indications of strong pre-catalog earthquakes) we argue, that the maximum credible earthquake for the VBTF is in the range M _max = 6.0–6.8, significantly larger than the magnitude of the strongest recorded events (M = 5.7).     

Recent plumbing system of the Krakatau volcano revealed by teleseismic earthquake distribution

Spatial and temporal analysis of global seismological data 1964–2005 reveals a distinct teleseismic earthquake activity producing a columnar-like formation in the continental wedge between the Krakatau volcano at the surface and the subducting slab of the Indo-Australian plate. These earthquakes occur continuously in time, are in the body-wave (m _b) magnitude range 4.5–5.3 and in the depth range 1–100 km. The Krakatau earthquake cluster is vertical and elongated in the azimuth N30°E, suggesting existence of a deep-rooted fault zone cutting the Sunda Strait in the SSW-NNE direction. Possible continuation of the fault zone in the SW direction was activated by an intensive 2002/2003 aftershock sequence, elongated in the azimuth of N55°E. Beneath the Krakatau earthquake cluster, an aseismic gap exists in the Wadati-Benioff zone of the subducting plate at the depths 100–120 km. We interpret this aseismic gap as a consequence of partial melting inhibiting stress concentration necessary to generate stronger earthquakes, whereas the numerous earthquakes observed in the overlying lithospheric wedge beneath the volcano probably reflect magma ascent in the recent plumbing system of the Krakatau volcano. Focal depth of the deepest events (~100 km) of the Krakatau cluster constrains the location of the primary magma generation to greater depths. The ascending magmatic fluids stress fault segments within the Sunda Strait fault zone and change their friction parameters inducing the observed tectonic earthquakes beneath Krakatau.     

Source characteristics of the Yutian earthquake in 2008 from inversion of the co-seismic deformation field mapped by InSAR

On 21 March 2008, an Ms7.3 earthquake occurred at Yutian County, Xinjiang Uygur Autonomous Region, which is in the same year as 2008 Mw 7.9 Wenchuan earthquake. These two earthquakes both took place in the Bayar Har block, while Yutian earthquake is located in the west edge and Wenchuan earthquake is in the east. The research on source characteristics of Yutian earthquake can serve to better understand Wenchuan earthquake mechanism. We attempt to reveal the features of the causative fault of Yutian shock and its co-seismic deformation field by a sensitivity-based iterative fitting (SBIF) method. Our work is based on analysis and interpretation to high-resolution satellite (Quickbird) images as well as D-InSAR data from the satellite Envisat ASAR, in conjunction with the analysis of seismicity, focal mechanism solutions and active tectonics in this region. The result shows that the 22 km long, nearly NS trending surface rupture zone by this event lies on a range-front alluvial platform in the Qira County. It is characterized by distinct linear traces and a simple structure with 1–3 m-wide individual seams and maximum 6.5 m width of a collapse fracture. Along the rupture zone are seen many secondary fractures and fault-bounded blocks by collapse, exhibiting remarkable extension. The co-seismic deformation affected a big range 100 km × 40 km. D-InSAR analysis indicates that the interferometric deformation field is dominated by extensional faulting with a small strike-slip component. Along the causative fault, the western wall fell down and the eastern wall, that is the active unit, rose up, both with westerly vergence. The maximum subsidence displacement is ～2.6 m in the LOS, and the maximum uplift is 1.2 m. The maximum relative vertical dislocation reaches 4.1 m, which is 10 km distant from the starting rupture point to south. The 42 km-long seismogenic fault in the subsurface extends in NS direction as an arc, and it dipping angle changes from 70° near the surface to 52° at depth ～10 km. The slip on the fault plane is concentrated in the depth range 0–8 km, forming a belt of length 30 km along strike on the fault plane. There are three areas of concentrating slip, in which the largest slip is 10.5 m located at the area 10 km distant from the initial point of the rupture.     

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

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

On the aftershock sequence of a 4.6 mb earthquake of the Garhwal Himalaya

Locally recorded data for eighteen aftershocks of a magnitude(m_b) 4.6 earthquake occurring near Ukhimath in the Garhwal Himalaya were analysed. A master event technique was adopted to locate seventeen individual aftershock hypocentres relative to the hypocentre of the eighteenth aftershock chosen as the master event. The aftershock epicentres define an approximately 30 km² rupture zone commensurate with the magnitude of the earthquake. The distribution of epicentres within this zone and the limited amount of first motion data support the view that a group of parallel, sub-vertical, sinistral strike-slip faults oriented N46°, transverse to the regional NW-SE trend of the Garhwal Himalaya, was involved in this seismic episode. Since the estimated focal depth range for aftershocks of this sequence is 3–14 km, we infer that this transverse fault zone extends through the upper crustal layer to a depth of 14 km at least.     

The Mw = 6.0, 7 September 1999 Athens Earthquake

On 7 September 1999 at 11:56 GMT a destructive earthquake (M_w = 6.0) occurred close to Athens (Greece). The rupture process is examined using data from the Cornet local permanent network, as well as teleseismic recordings. Data recorded by a temporary seismological network were analyzed to study the aftershock sequence. The mainshock was relocated at 38.105°N, 23.565°E, about 20 km northwest of Athens. Four foreshocks were also relocated close to the mainshock. The modeling of teleseismic P and SH waves provides a well-constrained focal mechanism of the mainshock (strike = 105°, dip = 55° and rake = -80°) at a depth of 8 km and a seismic moment M₀ = 1.010²⁵ dyn·cm. The obtained fault plane solution represents normal faulting indicating an almost north-south extension. More than 3500 aftershocks were located, 1813 of which present RMS < 0.1 s and ERH, ERZ < 1.0 km. Two main clusters were distinguished, while the depth distribution is concentrated between 2 and 11 km. Over 1000 fault plane solutions of aftershocks were constrained, the majority of which also correspond to N–S extension. No surface breaks were observed but the fault plane solution of the mainshock is in agreement with the tectonics of the area and with the focal mechanisms obtained by aftershocks. The hypocenter of the mainshock is located on the deep western edge of the fault plane. The relocated epicenter coincides with the fringe that represents the highest deformation observed on the differential interferometric image. The calculated source duration is 5 sec, while the estimated dimensions of the fault are 15 km length and 10 km width. The source process is characterized by unilateral eastward rupture propagation, towards the city of Athens. An evident stop phase observed in the recordings of the Cornet local stations is interpreted as a barrier caused by the Aegaleo Mountain.     

Source parameters of 1<Superscript>st</Superscript> April 2015 Chamoli earthquake (Mw 4.8) vis-à-vis seismotectonics of the region

The paper presents a detailed analysis of 1^st April 2015 earthquake, whose epicenter (30.16° N, 79.28° E) was located near Simtoli village of Chamoli district, Uttarakhand. The focal depth is refined to 7 km by the grid search technique using moment tensor inversion. The source parameters of the earthquake as estimated by spectral analysis method suggested the source radius of ~1.0 km, seismic moment as 1.99E+23 dyne-cm with moment magnitude (Mw) of 4.8 and stress drop of 69 bar. The fault plane solution inferred using full waveform inversion indicated two nodal planes, the northeast dipping plane having strike 334° and dip 5° and the southwest dipping plane with dip 86° and strike 118°. The parallelism of the nodal plane striking 334° with dip 5° as indicated in depth cross sections of the tectonic elements suggested the north dipping Main Boundary Thrust (MBT) to be the causative fault for this earthquake. Spatio-temporal distribution of earthquakes during the period 1960-2015 showed seismic quiescence during 2006-2010 and migration of seismicity towards south.     

The seismic anisotropy beneath Southern Sakhalin: Evidence from the S-wave parameters of deep-focus earthquakes

The measurements of the parameters of split shear (S) waves from local deep-focus earthquakes recorded in 2005–2007 by a network of 12 seismic stations in Southern Sakhalin are presented. The results revealed the heterogeneous distribution of the anisotropic properties beneath Southern Sakhalin. The azimuths of the fast S-wave polarization beneath the stations in the central part of the peninsula are oriented along the NNW and NNE-NE directions normal to and along the Kuril Trench. Beneath the stations located along the western and eastern coasts, the azimuths of the fast S-wave polarization change their direction from NNW in the northern area to E-SE in the southern area. The highest anisotropy degree (up to 0.9–1.5%) is recorded beneath the central part of Southern Sakhalin. The maximum values of the discrepancy in the arrival time of the split S-waves are observed when the azimuth of the fast S-wave is oriented along the NNE beneath the active fault zones. The analysis of the variations of the S-wave lag time shows their weak depth dependence. The highest anisotropy is assumed in the upper layers of the medium (down to a depth of about 250 km). The results obtained for the dominating wave frequency of 1–5 Hz represent mainly the medium-scale anisotropy of the top of the studied region.     

3D modeling of earthquake cycles of the Xianshuihe fault,southwestern China

We perform 3D modeling of earthquake generation of the Xianshuihe fault, southwestern China, which is a highly active strike-slip fault with a length of about 350 km, in order to understand earthquake cycles and segmentations for a long-term forecasting and earthquake nucleation process for a short-term forecasting. Historical earthquake data over the last 300 years indicates repeated periods of seismic activity, and migration of large earthquake along the fault during active seismic periods. To develop the 3D model of earthquake cycles along the Xianshuihe fault, we use a rate- and state-dependent friction law. After analyzing the result, we find that the earthquakes occur in the reoccurrence intervals of 400–500 years. Simulation result of slip velocity distribution along the fault at the depth of 10 km during 2694 years along the Xianshuihe fault indicates that since the third earthquake cycle, the fault has been divided into 3 parts. Some earthquake ruptures terminate at the bending part of the fault line, which may means the shape of the fault line controls how earthquake ruptures. The change of slip velocity and displacement at 10 km depth is more tremendous than the change of the shallow and deep part of the fault and the largest slip velocity occurs at the depth of 10 km which is the exact depth of the seismic zone where fast rupture occurs.     

Modeling of source parameters of the 15 December 2015 Deogarh earthquake of M<Subscript>w</Subscript> 4.0

We herein present source parameters and focal mechanism of a rare cratonic upper crustal earthquake of M_w4.0, which occurred at 8 km depth (centroid depth) below a region near Deogarh, Jharkhand. For our study, we used broadband waveform data from a seismic network of 15 three-component seismographs in the eastern Indian craton. The average seismic moment, moment magnitude and source radius are estimated to be 1.1 × 10¹⁵ N-m, 4.0 and 180.6 m, respectively. The high average stress drop of 14.27 MPa could be attributed to its lower-crustal origin. The mean corner frequency is calculated to be 4.1 Hz. To study the source mechanism, we perform a deviatoric constrained full waveform moment tensor inversion of multiple point sources on the band-passed (0.06 – 0.14 Hz) broadband displacement data of the Deogarh event, using ISOLA software. The best fit is obtained for the source at 8 km centroid depth, with a moment magnitude 3.7, and a right-lateral strike-slip mechanism with strike 162°, dip 72° and rake 169°. The P-axis orients N24°E, which is parallel to the direction of the absolute plate motion direction of Indian plate, while T-axis orients E-W, which is parallel to the strike of the pre-existing Damodar Graben (DG) of Gondwana age. The occurrence of this earthquake is attributed to the neotectonic reactivation of a fault associated with the E-W trending DG shear zone.     

A paleo-seismological study of the Dauki fault at Jaflong,Sylhet, Bangladesh: Historical seismic events and an attempted rupture segmentation model

A paleo-seismological study was conducted at Jaflong, Sylhet, Bangladesh, which is on the eastern part of the Dauki fault. The geomorphology around Jaflong is divided into the Shillong Plateau, the foothills, the lower terraces, and the alluvial plain from north to south. Because the foothills and lower terraces are considered to be uplifted tectonically, an active fault is inferred to the south of the lower terraces. This fault, which branches from the Dauki fault as a foreland migration, is known as the Jaflong fault in this paper. The trench investigation was conducted at the southern edge of the lower terrace. The angular unconformity accompanied by folding, which is thought to be the top of the growth strata, was identified in the trench. An asymmetric anticline with a steep southern limb and gentle northern limb is inferred from the back-tilted lower terrace and the folding of the gravel layer parallel to the lower terrace surface. The timing of the seismic event which formed the folding and unconformity is dated to between AD 840 and 920.The trench investigation at Gabrakhari, on the western part of the Dauki fault, revealed that the Dauki fault ruptured in AD 1548 (Morino et al., 2011). Because the 1897 great Indian earthquake (M ⩾ 8.0; Yeats et al., 1997) was caused by the rupture of the Dauki fault (Oldham, 1899), it is clear that the Dauki fault has ruptured three times in the past one thousand years. The timing of these seismic events coincides with that of the paleo-liquefactions confirmed on the Shillong Plateau. It is essential for the paleo-seismological study of the Dauki fault to determine the surface ruptures of the 1897 earthquake. The Dauki fault might be divided into four rupture segments, the western, central, eastern, and easternmost segments. The eastern and western segments ruptured in AD 840–920 and in 1548, respectively. The 1897 earthquake might have been caused by the rupture of the central segment.     

Characteristics of ring seismicity in different depth ranges before large and great earthquakes in the Sumatra region

Characteristics of the seismicity in depth ranges 0–33 and 34–70 km before ten large and great (M _w = 7.0−9.0) earthquakes of 2000–2008 in the Sumatra region are studied, as are those in the seismic gap zones where no large earthquakes have occurred since at least 1935. Ring seismicity structures are revealed in both depth ranges. It is shown that the epicenters of the main seismic events lie, as a rule, close to regions of overlap or in close proximity to “shallow” and “deep” rings. Correlation dependences of ring sizes and threshold earthquakes magnitudes on energy of the main seismic event in the ring seismicity regions are obtained. Identification of ring structures in the seismic gap zones (in the regions of Central and South Sumatra) suggests active processes of large earthquake preparation proceed in the region. The probable magnitudes of imminent seismic events are estimated from the data on the seismicity ring sizes.     

青藏高原东北缘推覆体构造稳定性定量化评价——以龙门山推覆体为例

2008年5月12日汶川特大地震震害调查及分析表明,目前以活动断裂和历史地震调查为重点的工程区域构造稳定性评价方法存在漏判与误判特大地震问题,从而为工程安全埋下重大安全隐患。以龙门山活动推覆体为例,在已有研究成果基础上,利用岩体结构控制论、拜尔利定律等普适性原理对龙门山地壳岩体结构力学特征、控震结构面的抗剪强度与地震震级的线性相关性、地震震级与抗震设防烈度的关系进行了定量研究,对评价区域构造稳定性的关键问题进行了探讨。结果表明,推覆体型活动地块边界带中的滑脱层是对推覆体区域构造稳定性起主要控制作用的构造结构面--控震结构面,地震震级与滑脱层的埋深、抗剪切强度存在显著相关性：8级地震的震源深度接近20 km、7级地震的震源深度接近14 km、6级地震的震源深度接近10 km,据此对研究区及邻近的古地震进行了深度核定,圈定了龙门山活动推覆体-岷山地块的6级以上强震可能发生的范围、对应Ⅶ－Ⅺ度的抗震设防烈度范围。此研究成果弥补了以往根据活动断裂-发震断裂-历史最大震级与对应地震烈度评价工程区域构造稳定性,因历史地震资料疏漏不全、活动断裂带研究平面与深度范围局限以及忽视区域构造稳定性的岩体力学实质而导致评价结果常常出现误判与漏判的诸多缺陷。     

Parameters of Coseismic Reverse‐ and Oblique‐Slip Surface Ruptures of the 2008 Wenchuan Earthquake,Eastern Tibetan Plateau

Abstract: On May 12th, 2008, the Mw7.9 Wenchuan earthquake ruptured the Beichuan, Pengguan and Xiaoyudong faults simultaneously along the middle segment of the Longmenshan thrust belt at the eastern margin of the Tibetan plateau. Field investigations constrain the surface rupture pattern, length and offsets related to the Wenchuan earthquake. The Beichuan fault has a NE-trending right-lateral reverse rupture with a total length of 240 km. Reassessment yields a maximum vertical offset of 6.5±0.5 m and a maximum right-lateral offset of 4.9±0.5 m for its northern segment, which are the largest offsets found; the maximum vertical offset is 6.2±0.5 m for its southern segment. The Pengguan fault has a NE-trending pure reverse rupture about 72 km long with a maximum vertical offset of about 3.5 m. The Xiaoyudong fault has a NW-striking left-lateral reverse rupture about 7 km long between the Beichuan and Pengguan faults, with a maximum vertical offset of 3.4 m and left-lateral offset of 3.5 m. This pattern of multiple co-seismic surface ruptures is among the most complicated of recent great earthquakes and presents a much larger danger than if they ruptured individually. The rupture length is the longest for reverse faulting events ever reported.     

The Lassee segment of the Vienna Basin fault system as a potential source of the earthquake of <Emphasis Type="Italic">Carnuntum</Emphasis> in the fourth century <Emphasis Type="SmallCaps">a.d.</Emphasis>

The Vienna Basin fault system is a slow moving (1–2 mm/y) active sinistral fault extending from the Alps through the Vienna Basin into the Carpathians. It comprises an array of NE-striking sinistral strike-slip segments, which differ both by their kinematic and seismologic properties. Among these, the Lassee segment 30 km east of Vienna is of particular interest for seismic hazard assessment as it shows a significant seismic slip deficit. The segment is located about 8 km from the Roman city of Carnuntum, for which archaeological data indicate a destructive earthquake in the fourth century a.d. (local intensity about 9 EMS-98). Mapping of the Lassee segment using 2D seismic, GPR, tectonic geomorphology and Pleistocene basin analysis shows a negative flower structure at a releasing bend of the Vienna Basin fault. The hanging wall of the flower structure includes a Quaternary basin filled with up to 100-m thick Pleistocene growth strata. Faults root in the basal detachment of the Alpine-Carpathian floor thrust at about 8 km depth. The active faults east of the flower structure offset a Middle Pleistocene terrace of the Danube River forming an up to 20-m high composite fault scarp. High-resolution GPR (40, 500 MHz) mapped at least four distinct surface-breaking faults along this scarp including three faults, which are covered by about 2 m of post-tectonic strata. The youngest fault offsets these strata and coincides with a 0.5-m high scarp. This scarp may be interpreted as the product of a single surface-breaking earthquake, provided that the mapped fault offset formed during coseismic surface rupture. Data indicate that the Lassee segment may well be regarded the source of the fourth century earthquake. The interpretation is in line with local attenuation relations indicating a source close to the damaged site, observed fault dimensions and the fault offsets recorded by GPR and morphology.     

Slip distribution and source parameters of the 20 July 2017 Bodrum-Kos earthquake (Mw6.6) from GPS observations

Greek-Turkish boundary near the cities Kos and Bodrum has been shaken on July 20, 2017 by a Mw6.6 earthquake. The mainshock is located offshore and did not generate an on-land surface rupture. Analyzing pre- and post-earthquake continuous/survey-type static GPS observations, we investigated co-seismic surface displacements at 20 sites to characterize source parameters and slip-distribution of the mainshock. Fault plane solutions as well as co-seismic slip distribution have been acquired through the inversion of co-seismic GPS displacements modeling the event as elastic dislocations in a half space. Fault plane solution shows a southward dipping normal-type fault segment extending a depth down to ~12 km, which remains within the brittle upper crust. Results from the distributed slip inversion show that the mainshock activated a ~65 km fault section, which has three high slip patches, namely western, central and eastern patches, where the coseismic slips reach up to 13, 26, and 5 cm, respectively. This slip pattern indicates that the pre-earthquake coupling, which is storing the slip deficit, occurred on these three patches.