The stress field in Australia from composite fault plane solutions of the strongest earthquakes in the continent

The purpose of this study is to make an estimate of the stress field in the Australian continent from composite fault plane solutions of the strongest events in the area. 21 seismic events with magnitude greater than 5.0 have been chosen from a database of 400 Australian earthquakes since 1902 for performing composite fault plane solutions. The choice is based on the availability of data for first P-wave arrivals. Using data from ISC bulletins we have calculated composite fault plane solutions for six major earthquake zones in the country. An individual fault plane solution for at least one event in each zone has been computed, too. The space distribution of the focal mechanisms confirms a dominant east-west compression stress in the continent, with different tilt in different areas. The main compression axes in the Tennant Creek area (Northern Australia) are with clear north-south direction. The general direction of the plate motion is north-northeast.     

2015年尼泊尔M_S8.1强震对中国大陆静态应力触发影响的初探

本研究基于地震危险性分析中地震重复原则和构造外推原则,以及一段时期内活动断裂上发生地震机制是相似或相近的假设,利用中国已发生的震级为MS5.0以上的地震震源机制解,将其中的一个节面视为实际地质断层参数和静态应力触发计算中的接收断层面,计算了2015年4月25日尼泊尔MS8.1地震对中国大陆地区的静态应力触发情况.研究结果表明:2015年4月25日尼泊尔MS8.1地震,触发了其后续的2次强余震的发生;本次地震对中国大陆产生的应力变化量值很小,其产生的应力加载主要集中在其邻近的西藏和新疆地区的部分断层上,而其余地区的断层则主要受到了应力卸载作用.     

Relative location and source mechanism of inland earthquakes in Northern Egypt

We have relocated 259 inland earthquakes in northern Egypt using the double-difference hypocenter technique. Among this dataset we are able to determine source mechanisms of 200 events using P-wave polarities and amplitude ratios as well. The studied earthquakes have been recorded by the Egyptian National Seismological Network from October 1997 to December 2006 with local magnitude (M_L) varies between 1.5 and 5.0. Three earthquake dislocations have been defined namely: Dahshour, southeast Beni-Suef, and Cairo-Suez district. Earthquake activities tend to occur in clusters along the first dislocation (Dahshour) however, relatively scattered along the second (southeast Beni-Suef) and the third (Cairo-Suez district) dislocations. At Dahshour dislocation three distinct clusters have been distinguished. Source mechanism solutions of Dahshour earthquakes displayed normal faulting with a strike-slip component to strike-slip faulting with a minor normal dip-slip component. Most of earthquake focal mechanism orientations are varying from NE-SW to NW-SE. The fault plane solutions of Beni-Suef earthquakes represented normal faulting with a strike-slip component. If the NNW-SSE striking plane has been chosen to be the actual fault plane, some solutions would indicated normal faulting with a sinistral strike-slip motion and other reflect normal faulting with a dextral strike-slip component. The fault plane solutions of Cairo-Suez district earthquakes are compatible with E-W to ENE-WSW striking normal fault with a dextral strike-slip motion.     

华北强震断层面解和震源深度特征

本文利用远震P波及SH波和区域地震P_(nl)波波形反演方法,测定了华北9次强震的断层面解、震源深度及地震矩.结合其他作者测定的7次地震的结果,分析了华北16次强震断层面解及华北地震震源深度特征.结果表明,华北多数强震为走滑型地震,也有个别正断层型及逆断层型地震;16次地震多数发生在地壳中部10—25km深度范围内.     

华北强震断层面解和震源深度特征

本文利用远震P波及SH波和区域地震P_(nl)波波形反演方法,测定了华北9次强震的断层面解、震源深度及地震矩.结合其他作者测定的7次地震的结果,分析了华北16次强震断层面解及华北地震震源深度特征.结果表明,华北多数强震为走滑型地震,也有个别正断层型及逆断层型地震;16次地震多数发生在地壳中部10-25km深度范围内.     

2010年以来四川地区中强地震震源机制反演及深度确定

2008年5月12日四川龙门山断裂带发生了汶川8.0级地震,之后四川境内发生了两次7.0级地震（其中一个是芦山地震）,为了研究汶川地震之后龙门山断裂带及周边区域的地震活动性,本研究收集了国家地震台网和四川区域地震台网2010年1月1日-2017年12月31日四川地区发生的17次M ≥ 5.0地震以及120多次5.0 > M ≥ 4.0地震的波形资料,利用波形拟合法反演了震源机制解及区域应力场.反演结果显示,位于龙门山断裂带上的地震,震源机制以逆冲型为主,鲜水河断裂带地震震源机制以走滑型为主,而川滇块体西南部的理塘断裂、金沙江断裂附近,震源机制解以正断层为主.根据震源机制解反演得到的龙门山地区、鲜水河地区的主压应力场方向为WNW、近EW向.川滇块体的巴塘、理塘等地区,其主压应力轴方向为12°左右,接近SN向,且仰角接近40°左右.本研究利用面波振幅谱特征对震源深度进行了精确定位,定位结果与中国地震台网中心（CENC）,美国地震调查局（USGS）,国际地震中心（ISC）等机构地震目录进行了对比.结果显示,四川地区强震震源深度主要分布在20 km以上的中上地壳.龙门山地区震源优势分布在10~20 km,鲜水河断裂地震震源深度在10 km左右,川滇块体西南部的理塘断裂,巴塘断裂,金沙江断裂等地区,震源深度一般在5~10 km范围.

     

2008年盈江地震序列的震源参数和震源机制相关系数研究

基于云南数字地震台网记录的2008年盈江地震序列的数字波形资料, 采用波谱分析方法和Brune震源模型, 得到盈江地震序列的震源参数。 利用两次地震事件的相同台站的震源谱参数(零频振幅)计算谱振幅相关系数, 据此对地震的震源机制进行聚类分组, 并收集和对比地震序列已知的震源机制解结果, 发现每组内震源机制解P轴的相关性较好, 且相关系数越大, P轴的方位角就越接近。 将盈江地区划分为三个研究区, 联合震源参数的应力降和聚类分组中每组的平均震源机制解结果研究地震序列发生过程中应力的释放水平和应力场方向变化特征。 结果表明: 不同阶段震源机制解类型的变化和转化特征一定程度上反映了孕震过程中区域应力场随时间的变化特征, 并且地震震源机制解类型在时间段上的集中并向区域构造应力场方向转换的现象可能是发生强震的标志。 震源机制解分组类型和对应类型的地震的应力降有一定依赖关系, 震源机制解类型反映的应力场与区域应力场接近的地震应力降高, 震源机制解类型反映的应力场与区域应力场差距较大的地震应力降普遍较低。     

Strong motion records and synthetic isoseismals of the Pyrgos,Peloponnisos, Southern Greece,Earthquake Sequence of March 26, 1993

In this work, several seismological observations are presented in order to explain characteristic features of the earthquake sequence which occurred in March 1993 in southwestern Greece, very close to the city of Pyrgos.Fault plane solutions of the largest fore- and aftershocks and the main shock, as well as the directions at which the maximum ground accelerations were recorded suggest that this earthquake sequence has been developed by rupturing three distinct focal planes with different focal mechanisms. The first focal plane, located in the off-shore area, strikes NW-SE, dips SE and includes most of the foreshock activity. The foreshock activity migrated to the northeastern part of the city of Pyrgos and took place on planes with a predominant direction NE-SW. The main shock ofM _x =5.5 occurred in a focal plane located between the two above-mentioned areas. Strong motion records of significant shocks of the sequence show peak acceleration values on components consistent with the relevant fault plane solutions.Furthermore, the observed macroseismic field has been compared with synthetic isoseismals computed by using a certain velocity model and the focal mechanism parameters of the main shock.     

2011年腾冲中强地震序列震源机制和发震构造研究

2011年腾冲地区连续发生了3次中强地震,本文运用全波形模拟的方法研究了这3次中强地震序列的震源机制解,并对地震序列进行了重定位。笔者结合由震源机制解及重定位结果对其发震构造进行了分析,认为3次中强地震及序列均位于同一发震构造,呈北西走向,破裂面倾角陡直,在近水平的北北东向压应力作用下作右旋水平走滑错动,断层较浅仅限于上地壳,活动过程有向深部发展的趋势;震源区附近无明显已知构造与发震构造相对应,发震构造可能为未知隐伏构造。     

The aftershock sequence of the 5 August 2014 Orkney earthquake (M<Subscript>L</Subscript> 5.5), South Africa

More than 1000 aftershocks were recorded within a month after the occurrence of the M_L 5.5, 5 August 2014 Orkney earthquake. The events were relocated using the double difference method as part of an effort to identify the fault which might be the source of the events. A north–south trend of seismicity was revealed by the relocated events, with a diffuse cluster to the north of the main event. A depth profile shows these two clusters: one at a depth of about 2 km to the north of the main event and the other at depth between 3 and 6 km south of the main event. Focal mechanism solutions of 18 aftershocks were determined using first motion polarities from seismic stations of the Council for Geoscience cluster networks. Stress inversion analysis results from the focal mechanism solutions show a dominant extensional stress field in the region; the main event had a strike-slip fault plane solution. This is consistent with the regional stress field which is predominantly related to the East African rift system. It is possible that the occurrence of the main event triggered seismicity on shallower faults within the mining horizons oriented in a different direction to the fault on which the main event occurred. The area has a complex heterogeneous faulting structure as indicated by the observed low p values and complex focal mechanism solutions.     

Earthquake focal mechanisms and stress inversion in the Irpinia Region (southern Italy)

The goal of this study was to estimate the stress field acting in the Irpinia Region, an area of southern Italy that has been struck in the past by destructive earthquakes and that is now characterized by low to moderate seismicity. The dataset are records of 2,352 aftershocks following the last strong event: the 23 November 1980 earthquake (M 6.9). The earthquakes were recorded at seven seismic stations, on average, and have been located using a three-dimensional (3D) P-wave velocity model and a probabilistic, non-linear, global search technique. The use of a 3D velocity model yielded a more stable estimation of take-off angles, a crucial parameter for focal mechanism computation. The earthquake focal mechanisms were computed from the P-wave first-motion polarity data using the FPFIT algorithm. Fault plane solutions show mostly normal component faulting (pure normal fault and normal fault with a strike-slip component). Only some fault plane solutions show strike-slip and reverse faulting. The stress field is estimated using the method proposed by Michael (J Geophys Res 92:357–368, 1987a) by inverting selected focal mechanisms, and the results show that the Irpinia Region is subjected to a NE–SW extension with horizontal σ ₃ (plunge 0°, trend 230°) and subvertical σ ₁ (plunge 80°, trend 320°), in agreement with the results derived from other stress indicators.     

2015年7月3日皮山6.5级地震发震构造初步研究

基于新疆区域数字地震台网记录,采用CAP（Cut and Paste）方法反演了2015年7月3日皮山6.5级主震和部分M_S3.6以上余震的震源机制解和震源深度;采用HypoDD方法重新定位了序列中M_L2.5以上地震序列的震源位置,并利用小震分布和区域应力场拟合了可能存在的发震断层面参数.基于上述研究,综合分析了皮山6.5级地震序列的震源深度、震源机制和震源破裂面特征,探讨可能的发震构造.结果显示,利用CAP方法得到的最佳双力偶机制解节面I：走向280°/倾角60°/滑动角90°;节面Ⅱ：走向100°/倾角30°/滑动角90°,矩心深度19 km,表明该地震为一次逆冲型地震事件.大部分M_S3.6以上余震震源机制与主震具有一定的相似性.双差定位结果显示,M_L2.5以上的余震序列主要分布在主震的西南方向,深度主要分布在0～15 km范围内,余震分布显示出与发震构造泽普隐伏断裂一致的倾向南西的特征.利用小震分布和区域应力场拟合得到发震断层参数为走向104°/倾角34°/滑动角94°,该结果与主震震源机制解中节面Ⅱ的滑动角较为接近,绝大多数余震发生在断层面附近10 km左右的区域.根据本研究得到的震源机制、精定位结果以及利用小震分布和区域应力场拟合得到的断层面的参数,结合震源区地质构造情况,初步给出了此次皮山6.5级地震的发震模式.     

2013年7月22日岷县漳县6.6级地震序列震源机制解及其特征分析

针对2013年7月22日甘肃岷县漳县发生的M_s6.6地震序列,采用CAP(Cut and Paste)方法反演了岷县漳县6.6级地震及部分强余震的震源机制解。结果显示,6.6级主震最佳双力偶解节面I走向189°,倾角51°,滑动角142°;节面II参数走向305°,倾角61°和滑动角46°,主震为逆冲兼走滑型,矩心震源深度均为7km;8次M_s≥3.0余震震源机制解向NE倾的节面II的优势倾角约为52°,表现出逆冲分量大的特性。结合震区的活动构造、余震及烈度分布,判定节面II代表了相应地震的发震断层面,地震序列震源机制的特性反映了与临潭一宕昌断裂带相似的活动特征,分析认为,岷县漳县6.6级地震的发生与临潭-宕昌断裂的活动可能存在一定的相关性。     

Evidence of multiple strain fields beneath the eastern flank of Mt. Etna volcano (Sicily,Italy) deduced from seismic and geodetic data during 2003–2004

We carried out a study of the seismicity and ground deformation occurring on Mt. Etna volcano after the end of the 2002–2003 eruption and before the onset of the 2004–2005 eruption. Data were recorded by the permanent local seismic network run by Istituto Nazionale di Geofisica e Vulcanologia – Sezione di Catania and by geodetic surveys carried out in July 2003 and July 2004 on the GPS network. Most earthquakes were grouped in two main clusters located in the northeastern and southeastern sectors of the volcano. The areal distribution of seismic energy associated with the recorded earthquakes allowed us to highlight the main seismogenic areas of Mt. Etna. In order to better understand the kinematic processes of the volcano, 3D seismic locations were used to compute fault plane solutions, and a selected dataset was inverted to determine stress and strain tensors. The focal mechanisms in the northeastern sector show clear left-lateral kinematics along an E-W fault plane, consistent with events occurring along the Pernicana Fault system. The fault plane solutions in the southeastern sector show mainly right-lateral kinematics along a NNE and ENE fault plane and left lateral-kinematics along NW fault planes that together suggest roughly E-W oriented compression. Surface ground deformation affecting Mt. Etna measured by GPS surveys highlighted a marked inflation during the same period and exceptionally strong seawards motion of its eastern flank. The 2D geodetic strain tensor distribution was calculated and the results show mainly ENE-WSW extension coupled with WNW-ESE contraction, indicating right-lateral shear along a NW-SE oriented fault plane. The different deformation of the eastern sector of the volcano, as measured by seismicity and ground deformation, must be interpreted by considering the different depths of the two signals. Seismic activity in the southeastern sector of volcano is located between 3 and 8 km b.s.l. and can be associated with a very strong additional E-W compression induced by a pressurizing source just westwards and at the same depth, located by inverting GPS data. Ground deformation, in contrast, is mainly affected by the shallower dynamics of the fast moving eastern flank which produces a shallower opposing E-W extension. The entire dataset shows that two different processes affect the eastern flank at the same time but at different depths; the boundary is clearly located at a depth of 3 km b.s.l. and could represent the décollement surface for the mobile flank.     

The 2014 Mihoub earthquake (Mw4.3), northern Algeria: empirical Green’s function analysis of the mainshock and the largest aftershock

On November 15, 2014, an Mw4.3 earthquake occurred 2 km west of Mihoub village, 60 km SE of Algiers. In this study, we retrieve the relative source-time functions of the mainshock and largest aftershock (Mw3.9) for rupture analysis using the empirical Green’s function method. The two events are nearly colocated with a smaller aftershock (Mw3.5), which is treated as the empirical Green’s function. Moreover, these three events have similar focal mechanisms, suggesting that deconvolution is well posed in this case. The three events were recorded by nine stations of the Algerian permanent network. We use mainly P-wave data. The focal mechanism solution shows dominant reverse faulting with a strong strike-slip component. The two nodal planes align almost E-W, dipping to the south, and NNE-SSW, dipping to the NW, respectively; the fault and auxiliary planes cannot be resolved from hypocenter locations alone because too few aftershocks were recorded by the permanent network. The results show unilateral rupture propagation to the ENE and complex rupture with multiple episodes for the mainshock. The largest aftershock shows similar behavior with slightly less pronounced directivity at some sites. The rupture directivity for the mainshock is estimated at about N66° E, and the rupture velocity is Vr = 0.66β. The E-W nodal plane of the best-fit focal mechanism is the preferred fault plane because it best agrees with the directivity direction and is consistent with the E-W faulting that dominates in the region.     

阿拉善左旗M_S5.8地震前后地震活动性分析和重定位

2015年4月15日15时39分,内蒙古自治区阿拉善左旗(39.78°N,106.34°E)发生M_S5.8级地震,震源深度为10km.本文通过收集中国地震科学台阵探测——喜马拉雅Ⅱ期布设在南北地震带北段的34个流动台站连续观测数据,首先以中国地震局台网中心目录给出的138个地震事件作为模板,采用匹配定位方法(MatchLocate)对阿拉善左旗M_S5.8地震前后(2015年1月1日至2015年7月31日)小震进行了检测和识别,检测后给出的地震事件高达1688个,相比于台网地震目录的地震事件数目增加了10倍多.然后利用双差层析成像方法(tomoDD)对研究区域内阿拉善M_S5.8地震前后(2015-01-01—2015-07-31)1级以上的地震序列进行重定位,重定位后的结果显示:主震震中为39.8023°N、106.4388°E,震源深度为20.597km,余震序列展布优势方向为NEE方向,且在深度方向呈近垂直分布,结合野外地质调查结果和震源机制解分析,认为该次地震事件为一次走滑型事件,其发震断裂可能是一条近E-W向的隐伏断层.     

Application of microseismic multiplet analysis to the Basel geothermal reservoir stimulation events

The stimulation of a geothermal well in Basel, Switzerland produced a distribution of microseismic event locations with an overall alignment in the direction of the maximum horizontal stress. Fault plane solutions of individual larger events indicated movements on fracture planes at an angle to the maximum horizontal stress that could not be reliably interpreted from the event locations. To obtain higher resolution images of the microseismic event locations, events with similar waveforms have been identified by multiplet analysis. A number of receivers were used in the multiplet processing to ensure each multiplet is represented by a unique group of waveforms. The location accuracy within each multiplet has been significantly improved using cross‐correlation to refine the shear‐wave traveltime picks. The distribution of events within each multiplet can be interpreted as being due to movements on a single fracture or a number of near parallel fractures. It is shown that whilst the overall distribution of events is around the direction of the maximum horizontal stress, the individual multiplets representing fracture planes have a variety of azimuths and dips.     

华北地区1900-1970年5级以上地震发震位置与震源机制研究

本文通过整理世界地震台网1900-1970年世界台站震相报告,对该时期发生在我国华北地区震级5级以上强震的震源位置及震源机制（断层面解）进行了测量.具体开展了以下三点工作：（1）使用1900-1970年之间Shide Circulars（BAASSC,1900-1912）,《国际地震资料汇编》（ISS）和EHB Bulletin的震相到时,对该时间范围内的22个地震事件重新定位,得到了21个地震事件的可靠结果;（2）使用1933年至1970年之间ISS的P波初动,对该时间范围内的15个地震事件求解震源机制解,得到了其中12个地震事件的震源机制解的可靠结果,从而丰富了华北地区基于全球地震台网观测的1900-1970年历史强震目录震源机制资料;（3）基于华北地区地震地质及区域构造应力场研究结果,推断了华北地区M≥5.0的部分缺乏地震资料观测的历史地震断层面参数,并对这种方法推断的断层面参数的可靠性进行了讨论.     

Analysis of induced seismicity for stress field determination and pore pressure mapping

The focal mechanisms of some one hundred microseismic events induced by various water injections have been determined. Within the same depth interval, numerous stress measurements have been conducted with the HTPF method. When inverted simultaneously, the HTPF data and the focal plane solutions help determine the complete stress field in a fairly large volume of rock (about 15×10⁶ m³). These results demonstrate that hydraulically conductive fault zones are associated with local stress heterogeneities. Some of these stress heterogeneities correspond to local stress concentrations with principal stress magnitudes much larger than those of the regional stress field. They preclude the determination of the regional stress field from the sole inversion of focal mechanisms. In addition to determining the regional stress field, the integrated inversion of focal mechanisms and HTPF data help identify the fault plane for each for each of the focal mechanisms. These slip motions have been demonstrated to be consistent with Terzaghi's effective stress principle and a Coulomb friction law with a friction coefficient ranging from 0.65 to 0.9. This has been used for mapping the pore pressure in the rock mass. This mapping shows that induced seismicity does not outline zones of high flow rate but only zones of high pore pressure. For one fault zone where no significant flow has been observed, the local pore pressure has been found to be larger than the regional minimum principal stress but no hydraulic fracturing has been detected there.     

用滑动方向拟合法反演唐山余震区的平均应力场

将每个震源机制解的两个节面交替取作尝试的断层面,使计算的断层面上的剪切力方向与观测的断层面的滑动方向夹角最小,从而根据29个震源机制解结果反演了唐山余震区三个分区的平均主应力方向和中等主应力相对大小。本文结果说明,唐山余震区的最大主应力轴水平,取近东西方向,较之1976年唐山大震前的方向可能水平地顺时针转动了大约30。对同样地区用平均P,B,T轴推断主应力方向的结果与用上述方法所得结果基本一致。所作的数值试验说明,在特定情况下,用平均P,B,T轴推断构造应力主轴方向时有可能出现系统偏差。