Local site effects in Latur, India: an empirical study based on the aftershocks of the Killari earthquake of September 29, 1993

The site amplification is estimated at five seismic stations of the Latur region using the horizontal to vertical spectral ratios of 33 aftershocks of the main Killari earthquake of September 29, 1993 (UTC). Spectral amplifications, ranging from a factor of 2–6 are found to vary with frequency at different places. Significant amplification is found at four sites within the Latur region, at Basavakalyan, Kasgi, Killari, and Mudgad Eakoji villages. Our results show a positive correlation between the site amplification and the damage pattern in area. The pattern and the nature of the site amplification estimated in the present study corroborates also with the analytical models and the borehole data indicating alternating layers of unconsolidated sediments and basaltic rocks.     

Rupture mechanism of the 1993 Killari earthquake, India: constraints from aftershocks and static stress change

The Killari earthquake of September 29, 1993 (Mw=6.2) in peninsular India triggered several aftershocks that were recorded by a network of 21 stations. We computed the change in regional static stress caused by coseismic slip on the earthquake rupture and correlated it with the aftershocks with a view to constrain some of the rupture parameters of this earthquake. We evaluated the six available estimates of fault plane solutions for this earthquake and concluded that reverse slip on a 42° dipping, N112° trending fault, which extends up to the surface from a depth of 7 km, produces maximum correlation between the increased static stress and aftershock distribution. Our analysis suggests that the majority of coseismic slip occurred on the part of the rupture that lies in the depth range of 3–6.5 km.     

3-D seismic structure of the source area of the 1993 Latur, India, earthquake and its implications for rupture nucleations

The Latur earthquake (M_w 6.1) of 29 September 1993 is a rare stable continental region (SCR) earthquake that occurred on a previously unknown blind fault. In this study, we determined detailed three-dimensional (3-D) P- and S-wave velocity (V_p, V_s) and Poisson's ratio (σ) structures by inverting the first P- and S-wave high-quality arrival time data from 142 aftershocks that were recorded by a network of temporary seismic stations. The source zone of the Latur earthquake shows strong lateral heterogeneities in V_p, V_s and σ structures, extending in a volume of about 90 × 90 × 15 km³. The mainshock occurred within, but near the boundary, of a low-V_p, high-V_s and low-σ zone. This suggests that the structural asperities at the mainshock hypocenter are associated with a partially fluid-saturated fractured rock in a previously unknown source zone with intersecting fault surfaces. This might have triggered the 1993 Latur mainshock and its aftershock sequence. Our results are in good agreement with other geophysical studies that suggest high conductivity and high concentration of radiogenic helium gas beneath the source zone of the Latur earthquake. Our study provides an additional evidence for the presence of fluid related anomaly at the hidden source zone of the Latur earthquake in the SCR and helps us understand the genesis of damaging earthquakes in the SCR of the world.     

First observation of microspherule from the infratrappean Gondwana sediments below Killari region of Deccan LIP,Maharashtra (India) and possible implications

A rare occurrence of a microspherule has been found in the infratrappean sediments, encountered below 338 m thick Deccan volcanic cover in KLR-1 scientific borehole, drilled in the epicentral zone of the 1993 Killari earthquake (Maharashtra, India). Palynological studies of the sediments indicate their age as Early Permian (Asselian, 298–295 Ma) for deposition. Transmission electron microscope studies reveal that the spherule from the infratrappeans, is having a similar composition to that of the Neoarchean amphibolite to granulite facies mid crustal basement. The spherule is non-spherical in nature, containing mostly FeO (10.70 ± 0.20 wt.%), CaO (13.8 ± 0.5 wt.%), Al₂O₃ (7.78 ± 0.30 wt.%), MgO (6.47 ± 0.3 wt.%), SiO₂ (47.46 ± 0.50 wt.%), TiO₂ (2.47 ± 0.3 wt.%), K₂O (1.89 ± 0.20 wt.%), and Cl (0.33 ± 0.05 wt.%). Since the Fe composition of the spherule is almost same as the basement rock (10.5 wt.%), and the chlorine content is also in the same range as the basement (0.04–0.24 wt.%), it would suggest possibility of an extraterrestrial impact over the Indian terrain during the erstwhile Gondwana sedimentation period that may be associated with the Permian–Triassic mass extinction, the most severe one in the Earth's history.     

井间地震时频域联合速度和衰减系数同步层析反演

井间地震资料中蕴含丰富的地下储层岩性和物性信息,利用井间地震直达波初至信息,通过层析反演可以得到两井之间的速度剖面。通过对井间地震直达波振幅衰减信息的处理,可以实现井间的吸收衰减层析成像。基于井间地震直达波射线路径的一致性,提出了井间地震时域直达波走时层析反演和频域质心频率衰减层析同步反演方法,增加了约束条件,增强了抗干扰能力。对模型和实际资料的处理结果表明,速度和吸收衰减同步反演方法,提高了速度层析反演精度,得到了对油气更为敏感的衰减属性,增强了井间地震在储层描述和油气预测中的应用能力。     

The 2nd century AD earthquake in central Italy: archaeoseismological data and seismotectonic implications

The 2nd century AD earthquake in central Italy is only known by an epigraph that mentions restorations to a damaged weighing-house at the ancient locality of Pagus Interpromium. The available seismic catalogues report this event with the conventional date of 101 AD, a magnitude M _aw of 6.3, and an epicentral location at the village of San Valentino in Abruzzo Citeriore, in the province of Pescara. In order to improve the knowledge of the damage pattern, we gathered all the archaeological data collected during modern excavations at sites located in the area, which were presumably struck by the earthquake. This information is mainly represented by (1) stratigraphic units due to the sudden collapse of buildings over still frequented floors; (2) stratigraphic units demonstrating restoration or re-building of edifices; (3) stratigraphic units formed as the result of the abandonment of sites or of their lack of frequentation for decades or centuries. Only stratigraphic evidence consistent with an earthquake occurrence during the 2nd century AD has been considered. The most recent archaeological material found in a collapsed unit is a coin of Antoninus Pius, dated at 147–148 AD. This may represent a post quem date very close to the occurrence of the earthquake. The gathered information, plus the stratigraphic data that excluded the earthquake occurrence at some sites, has allowed us to roughly delineate an area of possible damage, including the Sulmona Plain and surrounding areas. Comparisons between the possible 2nd century damage distribution and (i) the damage patterns of more recent historical events that have struck the investigated area, (ii) the distribution of virtual intensities obtained by simulating an earthquake having an epicenter in the Sulmona Plain and applying an intensity attenuation relationship and (iii) a shaking scenario obtained by modelling the activation of the major active fault of the Sulmona Plain area (the Mt. Morrone fault) have revealed consistency between the ancient earthquake and the activation of this fault. Since no other historical events can be attributed to this active fault, we conclude that the time that has elapsed since the last fault activation should be in the order of 1,850 years, i.e. a time span that is very close to the recurrence interval of Apennine seismogenic sources. Moreover, considering the fault length, the causative source may be responsible for earthquakes with M up to 6.6–6.7. The comparison between the presumed 2nd century damage and the shaking scenario suggests that the magnitude mentioned is consistent with the presumed effects of the ancient earthquake. Finally, considering that Sulmona (the most important town in the region investigated) is located in the middle of the Mt. Morrone fault hanging wall, we consider it as the probable epicentral area. Therefore, to summarise the information on the 2nd century AD earthquake, we can conclude that (i) it occurred shortly after 147–148 AD; (ii) a magnitude M _w 6.6–6.7 can be attributed to it and (iii) the probable macroseismic epicentral area was Sulmona.     

Relocation of aftershocks, focal mechanisms and stress inversion: Implications toward the seismo-tectonics of the causative fault zone of Mw7.6 2001 Bhuj earthquake (India)

The HYPODD relocation of 1172 aftershocks, recorded on 8–17 three-component digital seismographs, delineate a distinct south dipping E–W trending aftershock zone extending up to 35 km depth, which involves a crustal volume of 40 km × 60 km × 35 km. The relocated focal depths delineate the presence of three fault segments and variation in the brittle–ductile transition depths amongst the individual faults as the earthquake foci in the both western and eastern ends are confined up to 28 km depth whilst in the central aftershock zone they are limited up to 35 km depth. The FPFIT focal mechanism solutions of 444 aftershocks (using 8–12 first motions) suggest that the focal mechanisms ranged between pure reverse and pure strike slip except some pure dip slip solutions. Stress inversion performed using the P and T axes of the selected focal mechanisms reveals an N181°E oriented maximum principal stress with a very shallow dip (= 14°). The stress inversions of different depth bins of the P and T axes of selected aftershocks suggest a heterogeneous stress regime at 0–30 km depth range with a dominant consistent N–S orientation of the P-axes over the aftershock zone, which could be attributed to the existence of varied nature and orientation of fractures and faults as revealed by the relocated aftershocks.     

The rupture process during the 1999 Düzce, Turkey, earthquake from joint inversion of teleseismic and strong-motion data

Teleseismic and strong-motion data are inverted to determine the rupture process during the November 1999 Düzce earthquake in NW Turkey. The fault geometry, rise time and rupture velocity are determined from the aftershock distribution and preliminary inversions of the teleseismic data. Joint inversion of the teleseismic and strong-motion data is then carried out for the slip distribution. We obtain the strike 264°, dip 64°, rake −172°, seismic moment 5.0×10¹⁹ N m (M_w 7.1), and average stress drop 7 MPa. This earthquake was characterized by bilateral fault rupture and asymmetric slip distribution. Two asperities (areas of large slip) are identified, the eastern one being 1.5 times larger than the western one. The derived slip distribution is consistent with the aftershock distribution, surface rupture and damage. The point of rupture initiation in this Düzce earthquake coincided with the eastern tip of the aftershock distribution of the August 1999 Izmit earthquake.     

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.     

The Vrancea earthquake of March 4, 1977 — a multiple seismic event and its seismotectonic implications

A study concerning the multiplicity character of Vrancea earthquake of March 4, 1977 based on records of the Romanian seismic network is presented. Four separate shocks within a 9 s interval were recognized of which three with M = 6.0–6.9. Some explanations regarding the macroseismic intensity are made. Also analysed are the shaking duration, the influence of the Carpathian Mountains on the shape of the isoseismal curves and the focal mechanism for the first two main shocks.     

测井与地震资料联合反演预测K2灰岩及其含水性

阳煤五矿主要可采煤层15号煤层顶板发育的K₂灰岩不是良好的地震波反射界面,常规地震剖面很难连续追踪。测井曲线上的K₂灰岩表现为高密度和高视电阻率异常,采用密度与视电阻率两种测井曲线融合生成拟密度曲线,基于模型反演得到地层岩性数据体,从而识别灰岩的赋存形态与厚度;采用概率神经网络反演的方法,优选出9种地震属性,构成神经网络训练样本,对灰岩的孔隙度和视电阻率进行神经网络反演,预测灰岩的富水性。     

Observations of the performance of buildings during the 1985 Mexico earthquake,and structural design implications

Summary The high amplification of the 1985 Mexico City earthquake is explained by the large number of strong and nearly continuous cycles of 2 s period motion lasting for more than 30 s. The type of damage caused by the earthquake—particularly to engineered multi-storey buildings of high flexibility is described, and methods of adding damping and stiffness elements to reduce and resist earthquake demand forces are discussed.     

基于叠前深度偏移的AVO反演在神狐海域试采区的应用

为了更好地研究神狐海域试采区天然气水合物的展布形态和赋存状况,基于高精度层析速度反演后的叠前深度偏移地震资料进行AVO反演技术研究。为在天然气水合物富集区有效地进行AVO研究,AVO处理前运用叠前深度偏移技术,通过高精度层析反演不断迭代速度模型;然后通过偏移来改善资料的质量,结合实测钻井资料进行AVO正演分析、反演可行性分析、属性交会分析;最后将具有高横向分辨能力的地震资料和高垂向分辨能力的测井资料与地质相结合进行叠前同时反演。该方法可反演出纵、横波阻抗等属性,直接进行岩性和流体识别,结果表明,在AVO反演中引入叠前深度偏移技术,为反演提供了较高质量的叠前共反射点道集和精细的层速度模型,改善了AVO技术在试采区的应用效果,提高了反演的可靠性。反演出的各种属性为AVO分析提供了丰富的信息,能够比较准确地反映地下水合物矿体的展布情况,可用于后续水合物储层的分析及综合研究,对未来的水合物勘探开采具有重要的现实意义。     

On the validity of time-predictable model for earthquake generation in north-east India

North-east India is seismically very active and has experienced many widelydistributed shallow, large earthquakes. Earthquake generation model for the region was studied using seismicity data [(1906–1984) prepared by National Geophysical Data Centre (NGDC), Boulder Colorado, USA]. For establishing statistical relations surface wave magnitudes (M _s≥5·5) have been considered. In the region four seismogenic sources have been identified which show the occurrences of atleast three earthquakes of magnitude 5·5≤M _s≤7·5 giving two repeat times. It is observed that the time interval between the two consecutive main shock depends on the preceding main shock magnitude (M _p) and not on the following main shock magnitude (M _f) revealing the validity of time predictable model for the region. Linear relation between logarithm of repeat time (T) and preceding main shock magnitude (M _p) is established in the form of logT=cM _p+a. The values ofc anda are estimated to be 0–36 and 1–23, respectively. The relation may be used for seismic hazard evaluation in the region.     

Possibility of a major earthquake in the tokai district, Japan and its pre-estimated seismotectonic effects

Crustal deformations, tsunamis and seismic intensity are pre-estimated for a large hypothetical earthquake, which it is feared may occur in the Tokai district along the Nankai trough. The long-term seismic quiescence since 1854, as well as the high rate of the present crustal movements in the district, form the principal evidence for the risk of the approaching catastrophe. The location and the mode of faulting in this earthquake are hypothesized in reference to the source mechanisms of the recent and historical earthquakes there. The fault parameters thus assumed are as follows: dip direction: N30°W; dip angle: 25°; fault dimension : 100 km × 70 km; dislocation: 4 m (reverse dip-slip: 3.8 m; right-lateral strike-slip: 1.3 m). The following are the principal conclusions: (a) the eastern part of the epicentral region including the Point Omaezaki will rise up about 100 cm, whereas the western part covering Ise and Mikawa bays will subside about 10–30 cm; (b) the coast extending from Omaezaki to the Shima peninsula will receive tsunami waves as high as 3 m in maximum, which may be locally amplified by the factor 2 or more on the rias coast along the Shima peninsula; (c) the Tokai coastal region with thick alluvium layers may suffer seismic damages as severe as those experienced in the 1854 Ansei I earthquake.     

Normalized impedance function and the straightforward inversion scheme for magnetotelluric data

This paper investigates the performance of normalized response function obtained by normalizing the Cagniard impedance function by a suitable factor and then rotating the phase by 45‡ to make it purely real for homogeneous half-space and equal to the square root of the half-space resistivity. Two apparent resistivity functions based on respectively the real and imaginary parts of this response function are proposed. The apparent resistivity function using the real part contains almost the same information as that yielded by the Cagniard expression while the one using the imaginary part qualitatively works as an indicator of the number of interfaces in the earth model. The linear straightforward inversion scheme (SIS), developed by the authors employing the concept of equal penetration layers, has been used to validate the proposed apparent resistivity functions. For this purpose, several synthetic and field models have been examined. Five synthetic models are studied to establish the veracity of the new functions and two well-studied published field data sets are inverted through SIS for comparison. We noticed that the new function and SIS compliment each other and lead to better understanding of the data information and model resolution.     

Source study and tectonic implications of the 1995 Ventimiglia (border of Italy and France) earthquake (ML=4.7)

We analyse the source process and the aftershock distribution of the April 21, 1995, Ventimiglia, M_L=4.7 earthquake using the records of permanent high dynamic broad-band seismic stations and a temporary network deployed on land and at sea few hours after the earthquake. This event occurred on the western Mediterranean coast, near the border between Italy and France, at a depth of 9 km, at a point where Alpine tectonic units and Late Oligocene extensional structure overlap and are currently undergoing compressional stress. The focal solutions of the mainshock and three aftershocks depict a dominant reverse faulting with an important strike-slip component, which underlines two nodal planes: a NW–SE-dipping north fault and a NE–SW-dipping south fault. We operate a careful re-location of the aftershocks using a master-event technique and data from the temporal network and obtain a predominant NW–SE alignment. Then, we analyse the rupture process using an empirical Green function approach. We find that the mainshock broke a 0.5 to 1 km fault length and that the rupture propagated during 0.1–0.2 s probably in a SE direction. Those two arguments, together with the recent fault trace that exists close to the epicentre, leads us to propose that this event expresses the reactivation of an old transverse NW–SE structure with a dextral movement. This study thus emphasizes the role of inherited, deep-rooted, transcurrent features in the tectonic reactivation of this passive margin. It also underlines the importance of combining short-period and broad-band seismology to better resolve and understand regional tectonic processes in areas of moderate seismic activity and complex geology.     

Analysis of soil effects and distribution of damage from the Pyrgos 1993 (Greece) earthquake

Summary On 26 March, 1993, a moderate magnitude earthquake (M _s=5.5) occurred at 3 km epicentral distance from the town of Pyrgos, in Southern Greece, causing extensive damage to masonry houses. To explain the variability of seismic intensity over the town and to propose measures against future seismic activity, a microzonation study was undertaken which combined geological, geophysical and geotechnical investigations, site specific analyses of seismic ground response and detailed recording of structural damage. The analytical predictions of ground response are correlated to soil conditions and then used to identify (micro-)zones of sites with similar seismic response. Furthermore, they are compared to quantitative estimates of damage distribution over the town. It is concluded that the peak ground acceleration, normalized against the input peak seismic acceleration, is a function of the local soil conditions as well as the seismic excitation characteristics. Hence, it cannot be defined uniquely at a site, without reference to the seismic excitation. However, the normalized peak ground velocity and the acceleration response spectra are mainly functions of the soil conditions and can be used as criteria for the practical definition of (micro-)zones. The distribution of damage in various parts of the town is at least partially attributable to local soil effects. The small epicentral distance of the earthquake, connected with the direction of the fault rupture, as well as the quality and techniques of construction, are additional factors that may have influenced the extent and distribution of damage.     

福建水口水库地区微震双差法重新定位

本文采用“双差地震定位法”重新定位了福建数字地震台网正式运行以来发生的水口库区附近的地震,结果表明,记录的地震分二个时段,分别发生于不同的小区域内1999年1月~12底,发生于N26.35°~N26.40°,E118.67°~E118.70°,约5×3km2范围内,震源深度集中在9~10km范围内,整个震源的分布为一近水平的南北展布,具有一定的构造地震特征;2000年初至2005年末,地震主要发生于N26.39°~N26.41°,E118.62°~E118.64°,约2×2km2范围内,震源深度约2~4km,由于记录的地震震级小、个数少、时间跨度大,认为应是浅表应力局部调整的结果。