Near-regional CMT and multiple-point source solution of the September 5, 2012, Nicoya,Costa Rica Mw 7.6 (GCMT) earthquake

We use acceleration data from the Observatorio Vulcanologico y Sismologico, Universidad Nacional de Costa Rica (OVSICORI-UNA) and Laboratorio de Ingenieria Sismica, Universidad de Costa Rica (LIS-UCR) seismic network for the relocation and moment-tensor solution of the September 5, 2012, 14:42:03.35 UTC, Nicoya, Costa Rica earthquake (Mw 7.6 GCMT). Using different relocation methods we found a stable earthquake hypocenter, near the original OVSICORI-UNA location in the Nicoya Peninsula, NW Costa Rica at Lat 9.6943°N, Lon 85.5689°W, depth 15.3 km, associated with the subduction of the Cocos plate under Caribbean plate. Acceleration records at OVSICORI-UNA and LIS-UCR stations (94–171 km), at 0.03 < f < 0.06 Hz were used in the waveform inversion for a single-point centroid moment tensor (CMT). Using spatial grid search the centroid position was found at the depth of 30 km, situated at Lat 10.0559°N, Lon 85.4778°W, i.e. of about 41 km NNE from the epicenter. The centroid time is 14:42:18.89 UTC, i.e. 15.54 s later relative to the location-based origin time. The nodal plane (strike 318°, dip 27° and rake 115°) is the fault plane that agrees with the geometry of the subducted slab at Nicoya, NNW Costa Rica. Increasing the maximum studied frequency from 0.06 to 0.15 Hz, the multiple point source inversion model leads to two subevents. The first one was located near the centroid and the second subevent was situated 20 km along strike and 10 km down dip from the first subevent and 6 s later. The uncertainty of the source model was carefully examined using complementary inversion methods, viz the iterative deconvolution and non-negative least squares.     

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.     

2019年5月18日松原M5.1地震构造机制分析

基于区域地震台网的数字化波形资料,使用ISOLA方法对2019年5月18日吉林松原M5.1地震进行矩张量反演,研究地震的震源机制,并且收集了地震序列中M_L2.5以上地震的震源机制解,采用FMSI（focal mechanism stress inversion）方法反演震中区构造应力场。结果显示：松原M5.1地震的矩震级为4.9,矩心深度为6 km,双力偶分量为91.5%,主压应力P轴方位角、倾角分别为76°和3°,主张应力T轴方位角、倾角分别为166°和16°,震源机制解显示典型的构造地震特征;震中区构造应力场理论应力轴σ₁方位角、倾伏角分别为88.0°和0.9°,σ₂方位角、倾伏角分别为178.2°和9.6°,σ₃方位角、倾伏角分别为352.5°和80.4°,这一结果与区域构造应力场一致。推断认为区域构造应力场触发了2019年松原M5.1地震活动,地震震源机制解的北西向节面与震中区附近的第二松花江断裂现今活动性质完全一致,认为第二松花断裂可能是松原M5.1地震的发震断层。     

青海玉树Mw6.9级地震震源破裂过程

针对2010年青海玉树藏族自治州发生的Mw6.9(Ms7.1)级地震,利用地震波形资料和InSAR获取的同震位移资料,根据同震形成的地表位移干涉图,构建三段式断层模型,反演重建地震的破裂过程。研究显示本次地震断层面走向为119°,倾角79°,滑动角-2.2°,最大滑动量达到200cm,震源深度12.5km,地震标量地震矩为2.18×1026dyn·cm。震源破裂特征表明,玉树地震主要是沿甘孜—玉树断裂发生的左旋走滑破裂事件,反映了印度板块向北的推挤作用下,青藏高原东部不同次级块体东向不均匀挤出的运动学特征。     

Teleseismic finite-fault inversion of two <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> = 6.4 earthquakes along the East Anatolian Fault Zone in Turkey: the 1998 Adana and 2003 Bingöl earthquakes

The East Anatolian Fault Zone is a continental transform fault accommodating westward motion of the Anatolian fault. This study aims to investigate the source properties of two moderately large and damaging earthquakes which occurred along the transform fault in the last two decades using the teleseismic broadband P and SH body waveforms. The first earthquake, the 27 June 1998 Adana earthquake, occurred beneath the Adana basin, located close to the eastern extreme of Turkey’s Mediterranean coast. The faulting associated with the 1998 Adana earthquake is unilateral to the NE and confined to depths below 15 km with a length of 30 km along the strike (53°) and a dipping of 81° SE. The fixed-rake models fit the data less well than the variable-rake model. The main slip area centered at depth of about 27 km and to the NE of the hypocenter, covering a circular area of 10 km in diameter with a peak slip of about 60 cm. The slip model yields a seismic moment of 3.5?×?10¹⁸ N-m (M_w???6.4). The second earthquake, the 1 May 2003 Bingöl earthquake, occurred along a dextral conjugate fault of the East Anatolian Fault Zone. The preferred slip model with a seismic moment of 4.1?×?10¹⁸ N-m (M_w???6.4) suggests that the rupture was unilateral toward SE and was controlled by a failure of large asperity roughly circular in shape and centered at a depth of 5 km with peak displacement of about 55 cm. Our results suggest that the 1998 Adana earthquake did not occur on the mapped Göksun Yakap?nar Fault Zone but rather on a SE dipping unmapped fault that may be a split fault of it and buried under the thick (about 6 km) deposits of the Adana basin. For the 2003 Bingöl earthquake, the final slip model requires a rupture plane having 15° different strike than the most possible mapped fault.     

The 20 May 2016 Petermann Ranges earthquake: centroid location,magnitude and focal mechanism from full waveform modelling

Ground velocity records of the 20 May 2016 Petermann Ranges earthquake are used to calculate its centroid-moment-tensor in the 3?D heterogeneous Earth model AuSREM. The global-centroid-moment-tensor reported a depth of 12?km, which is the shallowest allowed depth in the algorithm. Solutions from other global and local agencies indicate that the event occurred within the top 12?km of the crust, but the locations vary laterally by up to 100?km. We perform a centroid-moment-tensor inversion through a spatiotemporal grid search in 3?D allowing for time shifts around the origin time. Our 3?D grid encompasses the locations of all proposed global solutions. The inversion produces an ensemble of solutions that constrain the depth, lateral location of the centroid, and strike, dip and rake of the fault. The centroid location stands out with a clear peak in the correlation between real and synthetic data for a depth of 1?km at longitude 129.8° and latitude –25.6°. A collection of acceptable solutions at this centroid location, produced by different time shifts, constrain the fault strike to be 304?±?4° or 138?±?1°. The two nodal planes have dip angles of 64?±?5° and 26?±?4° and rake angles of 96?±?2° and 77?±?5°, respectively. The southwest-dipping nodal plane with the dip angle of 64° could be seen as part of a near vertical splay fault system at the end of the Woodroffe Thrust. The other nodal plane could be interpreted as a conjugate fault rupturing perpendicular to the splay structure. We speculate that the latter is more likely, since the hypocentres reported by several agencies, including the Geoscience Australia, as well as the majority of aftershocks are all located to the northeast of our preferred centroid location. Our best estimate for the moment magnitude of this event is 5.9. The optimum centroid is located on the 20?km surface rupture caused by the earthquake. Given the estimated magnitude, the long surface rupture requires only ～4?km of rupture down dip, which is in agreement with the shallow centroid depth we obtained.     

Estimates of source parameters of <Emphasis Type="Italic">M</Emphasis>4.9 Kharsali earthquake using waveform modelling

This paper presents the computation of time series of the 22 July 2007 M 4.9 Kharsali earthquake. It occurred close to the Main Central Thrust (MCT) where seismic gap exists. The main shock and 17 aftershocks were located by closely spaced eleven seismograph stations in a network that involved VSAT based real-time seismic monitoring. The largest aftershock of M 3.5 and other aftershocks occurred within a small volume of 4 × 4 km horizontal extent and between depths of 10 and 14 km. The values of seismic moment (M _∘) determined using P-wave spectra and Brune’s model based on f ² spectral shape ranges from 10¹⁸ to 10²³ dyne-cm. The initial aftershocks occurred at greater depth compared to the later aftershocks. The time series of ground motion have been computed for recording sites using geometric ray theory and Green’s function approach. The method for computing time series consists in integrating the far-field contributions of Green’s function for a number of distributed point source. The generated waveforms have been compared with the observed ones. It has been inferred that the Kharsali earthquake occurred due to a northerly dipping low angle thrust fault at a depth of 14 km taking strike N279°E, dip 14° and rake 117°. There are two regions on the fault surface which have larger slip amplitudes (asperities) and the rupture which has been considered as circular in nature initiated from the asperity at a greater depth shifting gradually upwards. The two asperities cover only 10% of the total area of the causative fault plane. However, detailed seismic imaging of these two asperities can be corroborated with structural heterogeneities associated with causative fault to understand how seismogenesis is influenced by strong or weak structural barriers in the region.     

Source parameters and focal mechanisms of local earthquakes: Single broadband observatory at ISM Dhanbad

A three-component broadband seismograph is in operation since January 2007 at the Indian School of Mines (ISM) campus. We have used the broadband seismograms of two local earthquakes M <3 recorded by this single station to illustrate its efficacy in understanding the source processes and tectonics in Dhanbad area. Source parameters and fault plane solutions are obtained through waveform inversion. It is observed that these two earthquakes occurred in the lower crust at a depth of 26 km by strike slip faulting. North-south compressional and east-west tensional stresses are dominant in the area, and the lower crust is the source area for the local earthquakes.     

PRELIMINARY RESULTS OF EARTHQUAKE MOMENT TENSOR INVERSION IN WESTERN CHINA

PRELIMINARY RESULTS OF EARTHQUAKE MOMENT TENSOR INVERSION IN WESTERN CHINA     

基于D—InSAR的九寨沟地震同震形变场反演分析

利用覆盖九寨沟地区的RadarSat—2数据与Sentinel—1A数据,采用精轨数据进行定轨,消除轨道误差,并结合合成孔径差分(D-InSAR)方法中的双轨差分技术,获取2017年8月8日Mw7. 0级地震的同震形变场。结果表明,视线方向(LOS)最大沉降量约为20 cm,隆起量达9 cm。基于弹性半空间形变模型反演该地震的断层滑动分布,得出该地震断层滑动以左旋走滑为主,走向为330°,倾角为32°,滑动角为-170°,同震滑动分布主要集中在4～12 km深度处,最大滑动量位于9 km处,约为6. 14 m,平均滑动量为0. 57 m。反演获得的地震标量矩为4. 06E+18N·m,震级Mw约为6. 4,深度为19. 5 km。     

Tectonic structures across the East African Rift based on the source parameters of the 20 May 1990 M7.2 Sudan earthquake

Earthquakes in Kenya are common along the Kenya Rift Valley because of the slow divergent movement of the rift and hydrothermal processes in the geothermal fields. This implies slow but continuous radiation of seismic energy, which relieves stress in the subsurface rocks. On the contrary, the NW-SE trending rift/fault zones such as the Aswa-Nyangia fault zone and the Muglad-Anza-Lamu rift zone are the likely sites of major earthquakes in Kenya and the East African region. These rift/fault zones have been the sites of a number of strong earthquakes in the past such as the M _w = 7.2 southern Sudan earthquake of 20 May 1990 and aftershocks of M _w = 6.5 and 7.1 on 24 May 1990, the 1937 M _s = 6.1 earthquake north of Lake Turkana close to the Kenya-Ethiopian border, and the 1913 M _s = 6.0 Turkana earthquake, among others. Source parameters of the 20 May 1990 southern Sudan earthquake show that this earthquake consists of only one event on a fault having strike, dip, and rake of 315°, 84°, and ?3°. The fault plane is characterized by a left-lateral strike slip fault mechanism. The focal depth for this earthquake is 12.1 km, seismic moment M _o = 7.65 × 10¹⁹ Nm, and moment magnitude, M _w = 7.19 (?7.2). The fault rupture started 15 s earlier and lasted for 17 s along a fault plane having dimensions of ?60 km × 40 km. The average fault dislocation is 1.1 m, and the stress drop, ?σ, is 1.63 MPa. The distribution of historical earthquakes (M _w ≥ 5) from southern Sudan through central Kenya generally shows a NW-SE alignment of epicenters. On a local scale in Kenya, the NW–SE alignment of epicenters is characterized by earthquakes of local magnitude M _l ≤ 4.0, except the 1928 Subukia earthquake (M _s = 6.9) in central Kenya. This NW–SE alignment of epicenters is consistent with the trend of the Aswa-Nyangia Fault Zone, from southern Sudan through central Kenya and further southwards into the Indian Ocean. We therefore conclude that the NW–SE trending rift/fault zones are sites of strong earthquakes likely to pose the greatest earthquake hazard in Kenya and the East African region in general.     

THE HIGH RESOLUTION SEISMIC TOMOGRAPHIC IMAGE IN QINGHAI—TIBET PLATEAU AND ITS DYNAMIC IMPLICATIONS

THE HIGH RESOLUTION SEISMIC TOMOGRAPHIC IMAGE IN QINGHAI—TIBET PLATEAU AND ITS DYNAMIC IMPLICATIONSeasthenospherehadbe     

青海玉树Mw6.9级地震震源破裂过程

针对2010年青海玉树藏族自治州发生的Mw6.9（Ms7.1）级地震,利用地震波形资料和InSAR获取的同震位移资料,根据同震形成的地表位移干涉图,构建三段式断层模型,反演重建地震的破裂过程。研究显示本次地震断层面走向为119°,倾角79°,滑动角-2.2°,最大滑动量达到200cm,震源深度12.5km,地震标量地震矩为2.18×1026dyn·cm。震源破裂特征表明,玉树地震主要是沿甘孜—玉树断裂发生的左旋走滑破裂事件,反映了印度板块向北的推挤作用下,青藏高原东部不同次级块体东向不均匀挤出的运动学特征。     

Co‐seismic Faults and Geological Hazards and Incidence of Active Fault of Wenchuan Ms 8.0 Earthquake,Sichuan, China

Abstract: There are two co-seismic faults which developed when the Wenchuan earthquake happened. One occurred along the active fault zone in the central Longmen Mts. and the other in the front of Longmen Mts. The length of which is more than 270 km and about 80 km respectively. The co-seismic fault shows a reverse flexure belt with strike of N45°–60°E in the ground, which caused uplift at its northwest side and subsidence at the southeast. The fault face dips to the northwest with a dip angle ranging from 50° to 60°. The vertical offset of the co-seismic fault ranges 2.5–3.0 m along the Yingxiu-Beichuan co-seismic fault, and 1.5–1.1 m along the Doujiangyan-Hanwang fault. Movement of the co-seismic fault presents obvious segmented features along the active fault zone in central Longmen Mts. For instance, in the section from Yingxiu to Leigu town, thrust without evident slip occurred; while from Beichuan to Qingchuan, thrust and dextral strike-slip take place. Main movement along the front Longmen Mts. shows thrust without slip and segmented features. The area of earthquake intensity more than IX degree and the distribution of secondary geological hazards occurred along the hanging wall of co-seismic faults, and were consistent with the area of aftershock, and its width is less than 40km from co-seismic faults in the hanging wall. The secondary geological hazards, collapses, landslides, debris flows et al., concentrated in the hanging wall of co-seismic fault within 0–20 km from co-seismic fault.     

Multi-scale imaging of a slow active fault zone: contribution for improved seismic hazard assessment in the Swiss Alpine foreland

Seismic hazard assessment of slow active fault zones is challenging as usually only a few decades of sparse instrumental seismic monitoring is available to characterize seismic activity. Tectonic features linked to the observed seismicity can be mapped by seismic imaging techniques and/or geomorphological and structural evidences. In this study, we investigate a seismic lineament located in the Swiss Alpine foreland, which was discussed in previous work as being related to crustal structures carrying in size the potential of a magnitude M 6 earthquake. New, low-magnitude (?2.0 ≤ M_L ≤ 2.5) earthquake data are used to image the spatial and temporal distribution of seismogenic features in the target area. Quantitative and qualitative analyses are applied to the waveform dataset to better constrain earthquakes distribution and source processes. Potential tectonic features responsible for the observed seismicity are modelled based on new reinterpretations of oil industry seismic profiles and recent field data in the study area. The earthquake and tectonic datasets are then integrated in a 3D model. Spatially, the seismicity correlates over 10–15 km with a N–S oriented sub-vertical fault zone imaged in seismic profiles in the Mesozoic cover units above a major decollement on top of the mechanically more rigid basement and seen in outcrops of Tertiary series east of the city of Fribourg. Observed earthquakes cluster at shallow depth (<4 km) in the sedimentary cover. Given the spatial extend of the observed seismicity, we infer the potential of a moderate size earthquake to be generated on the lineament. However, since the existence of along strike structures in the basement cannot be excluded, a maximum M 6 earthquake cannot be ruled out. Thus, the Fribourg Lineament constitutes a non-negligible source of seismic hazard in the Swiss Alpine foreland.     

Earthquake source parameters at the sumatran fault zone: Identification of the activated fault plane

Fifteen earthquakes (M_w 4.1–6.4) occurring at ten major segments of the Sumatran Fault Zone (SFZ) were analyzed to identify their respective fault planes. The events were relocated in order to assess hypocenter uncertainty. Earthquake source parameters were determined from three-component local waveforms recorded by IRIS-DMC and GEOFON broadband lA networks. Epicentral distances of all stations were less than 10°. Moment tensor solutions of the events were calculated, along with simultaneous determination of centroid position. Joint analysis of hypocenter position, centroid position, and nodal planes produced clear outlines of the Sumatran fault planes. The preferable seismotectonic interpretation is that the events activated the SFZ at a depth of approximately 14–210 km, corresponding to the interplate Sumatran fault boundary. The identification of this seismic fault zone is significant to the investigation of seismic hazards in the region.     

On the initiation and movement mechanisms of a catastrophic landslide triggered by the 2008 Wenchuan (M<Subscript>s</Subscript> 8.0) earthquake in the epicenter area

The Niumiangou landslide (~7.5 × 10⁶ m³) was the largest that occurred in the town of Yingxiu (the epicentral area) during the 2008 Wenchuan earthquake. This landslide originated on a steep slope (~30°) that was located directly above the rupture surface of the responsible fault and then traveled ~2 km after flowing down the axes of two gently sloping (<12°) valleys. Evidence at the site indicates that the landslide materials were highly fluidized and underwent rapid movement. To examine the initiation and movement mechanisms of this landslide, we performed a detailed field survey, conducted laboratory tests on samples taken from the field, and analyzed the seismic motion. We conclude that the landside materials were displaced due to seismic loading during the earthquake and that liquefaction may have been triggered in saturated layers above the sliding surface with progressive downslope sliding, which resulted in the high mobility of the displaced materials. The liquefaction of colluvial deposits along the travel path due to loading by the sliding mass enhanced the mobility of the displaced mass originating in the source area. Using an energy-based approach, we estimated the dissipated energy in our cyclic loading test and the possible energy dissipated to the soil layer on the slope by the earthquake. We infer that the seismic energy available for the initiation of the slope failure in the source area may have greatly exceeded the amount required for the initiation of the liquefaction failure. The slope instability might have been triggered several seconds after the arrival of seismic motion.     

Source parameters of March 10 and September 13, 2007, United Arab Emirates earthquakes

On March 10 and September 13, 2007 two earthquakes with moment magnitudes 3.66 and 3.94, respectively, occurred in the eastern part of the United Arab Emirates (UAE). The two events were widely felt in the northern Emirates and Oman and were accompanied by a few aftershocks. Ground motions from these events were well recorded by the broadband stations of Dubai (UAE) and Oman seismological networks and provide an excellent opportunity to study the tectonic process and present day stress field acting in this area. In this study, we report the focal mechanisms of the two main shocks by two methods: first motion polarities and regional waveform moment tensor inversion. Our results indicate nearly pure normal faulting mechanisms with a slight strike slip component. We associated the fault plane trending NNE–SSW with a suggested fault along the extension of the faults bounded Bani Hamid area. The seismicity distribution between two earthquake sequences reveals a noticeable gap that may be a site of a future event. The source parameters (seismic moment, moment magnitude, fault radius, stress drop and displacement across the fault) were also estimated from displacement spectra. The moment magnitudes were very consistent with waveform inversion. The recent deployment of seismic networks in Dubai and Oman reveals tectonic activity in the northern Oman Mountains that was previously unknown. Continued observation and analysis will allow for characterization of seismicity and assessment of seismic hazard in the region.     

The Woods Reef (New South Wales) earthquake of 14 November 1990: Focal mechanism derived from amplitude ratios and synthetic seismograms

On the morning of 15 November 1990 local time, Armidale and the area to the west of Armidale was shaken by a magnitude 3.2 earthquake. The epicentre was located at 30.39° S, 150.88° E and the depth of focus at 12 ± 7 km. As the epicentre was close to the Peel Fault an attempt was made to constrain the focal mechanism of this earthquake. The conventional method, which is based on the analysis of P wave polarities, was not applicable because the event was not strong enough. In an alternative method, the amplitudes of various seismic phases recorded at a number of stations well distributed in azimuth were compared with theoretical amplitudes calculated with the reflectivity method for a point shear dislocation in a layered medium. The differences between observed and calculated amplitudes were minimized as a function of fault strike, fault dip and direction of the slip vector. The analysis indicates that none of the possible fault planes had the strike of the Peel Fault. The solution suggests predominantly strike slip motion along two possible, steeply dipping fault planes. The inferred direction of the maximum compressional stress. is east‐west which is in good agreement with other estimates of the stress field for eastern Australia.     

Waveform inversion and focal mechanisms of two weak earthquakes in Cordillera Principal (Argentina) between 35° and 35.5° S

Only few (six) focal mechanism, in CMT Catalog, have been so far known for intraplate shallow events in the Andean chain close to Chile–Argentina state border at latitudes ∼35° S. We add two more mechanisms, depths and moment magnitudes by carefully analyzing full waveforms of weak events recorded by broad-band stations of the Chile Argentina Geophysical Experiment (southern profile). The moment magnitudes of both events (Mw = 3.6 and 3.7) are lower than the duration magnitudes (Md = 4.0 and 4.29) reported by NEIC. The source depth, constrained by waveforms for one of the studied events (5.5–8.5 km) seems to be considerably shallower than the hypocenter depth located by means of arrival times (∼20 km). The waveform analysis was complemented by first-motion polarities which resulted in an uncertainty assessment of the focal mechanism. Event 1 (2001-11-03) has a strike-slip mechanism with a small normal component and almost vertical nodal planes in the north-south and east-west directions. The north-south nodal plane could be related to the Calabozos faults system. Event 2 (2002-02-16) has a strike-slip mechanism with a small thrust component. The latter event (its subhorizontal nodal plane) could be associated with the El Diablo-El Fierro fault system. Dextral strike-slip solutions are consistent with recent studies in the area.