The focal mechanism of two Peruvian earthquakes

Summary The focal mechanisms of two Peruvian earthquakes, January 15, 1958, and January 13, 1960, are determined using data from both the first motion ofP and the polarization ofS. The fault motions correspond to strike slip motion on neighboring faults corresponding in strike to geophysical evidences of regional faulting. The motion is left lateral in one of the shocks, right lateral in the other.     

Classical focal mechanism techniques for body waves

The initial motion of primary body waves and polarization directions of secondary body waves have been applied successfully to the study of the mechanism operating at the earthquake focus. Equal area plots of these body wave characteristics resulted in radiation patterns that were compared to patterns due to theoretical focal mechanism sources. Such an approach indicated that a double couple force is the source mechanism operating at earthquake foci. This can be physically represented by faulting at the focus. This seems likely because of the relation between earthquakes and fault motion in shallow earthquakes. The possibility of other mechanisms operating in deep focus earthquakes has not been ruled out.The technique of solving for fault plane orientation/motion in equal area plots ofP-wave first motion andS-wave polarization, has been a powerful tool in areas where fault motions cannot be directly observed in studying tectonics. Such an approach has been used to test the theory of global tectonics. This approach has resulted in the confirmation of the suggested mechanics of plate motion, and the results of plate motion, such as spreading of the sea floor from mid-oceanic ridges and underthrusting of lithospheric plates at the sites of oceanic trenches.     

Free and forced shallow-water oscillations in a rotating channel of parabolic cross-section

Free and forced oscillations of shallow water in an infinitely long rotating channel of parabolic cross-section are analyzed. The pure cross-channel oscillations ofChrystal (1905) and solutions for zero rotation first discussed byProudman (1925) andHidaka (1932) are special asymptotic solutions for the free modes of this model. However, for increasingly large, along-shore wave number, our solutions donot uniformly approach those ofReid (1958) andBall (1967) for a single shore-line and semi-infinite ocean. A method of computing eigen frequencies and eigen functions for the general problem is described, and a sufficient number of these are exhibited graphically to permit visualization of the transitions between the asymptotic regions.The forced problem consists of an incoming wave-train or surge generated at the center of the channel. Amplitude and transports near the shore are computed for a wide range of dimensionless incoming-wave frequencies and rotational frequencies.     

Subduction zone geometry and stress patterns in the tyrrhenian sea

Joint hypocenter determination is performed for intermediate and deep earthquakes of the Tyrrhenian Sea region.This analysis allowed us to obtain a catalogue of 70 well-located events in this peculiar Benioff zone, which is characterized by quite low seismic activity, compared to the Pacific deep earthquake regions. The method used for the analysis is that ofFrohlich (1979), a variant of the successive approximation technique, which allows use of a great number of events and stations but saves computer memory. The results show a spoon-shaped Benioff zone, dipping NW in the Tyrrhenian Sea to 500km depth. 32 reliable fault-plane solutions have been determined using these new earthquake locations, confirming the predominance of down-dip compression in the central part of the slab and more complex motion along the borders of the zone, as previously suggested byGasparini et al. (1982).     

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.     

On waves generated by a moving pressure field,II: Cnoidal wave effects

Summary In a recent paperHunt andTanner [3]²) investigated the waves generated by a steadily moving two-dimensional pressure distribution, which was zero ahead of the disturbance and a constantp ₀, tehind it, these regions being joined smoothly by a cubic function. Only those solutions with supercritical flow in both regions were considered, these were found to lead to an asymmetric solitary wave.This result is now extended to take account of the possibility of subcritical flow in either or both the regions, that is when there is a cnoidal wave train either behind and or ahead of the main solitary wave crest.The wave profiles are determined by the iterative method employed in the previous paper. This together with the wave drag associated with each system is computed for various values ofp ₀/U ², where is the fluid density andU a typical velocity.     

地震的震源机制解的计算和分析

本文根据大柴旦周围40余个地震监测台站记录到的P波初动及部分S波初动等信息计算了2008年以来大柴旦50余次M3.5级以上的地震,其中合理的得出了46个震源机制解,通过与哈佛给出的结果分析了两次6级地震震源机制解的可靠性,同时讨论了中小震的震源机制解特性,并对区域的应力场进行了估算,得到了较为合理的结果,为该区域的地震研究提供一定参考.     

P in the shadow zone of the earth's core—Part I

Summary Regional variations have been indicated in the slope of theP travel-time curve in the shadow zone of the earth's core. Further study is needed since the uncertainties of the slope are large, especially for the observations from North American stations. There is no significant difference between themean slope of theP travel-time curve in the 95°102.9 range and those obtained byJeffreys, andJeffreys andBullen. However, there is a significant difference between themean slope in the 103° to 135° range as obtained in this study, and those obtained byJeffreys andBullen, and in a later revision byJeffreys. Themean travel-time curve ofP in the shadow zone of the earth's core should be lowered. A trial travel-time table is given with amean slope of 4.41 sec/deg. This table is in close agreement with the times obtained byGutenberg andRichter, and with the trial travel-times ofLehmann. Under the assumption of diffraction the longitudinal wave velocity has been determined to be 13.7 km/sec at the core-mantle boundary.This paper was presented at the Annual Meeting of the Seismological Society of America Reno, Nevada, 1966.     

Über die Bindung der Aerosolteilchen auf der Oberfläche der Wolkenelemente

Summary The importance ofStefan's flow by collection of insoluble aerosol particles on the surface of the evaporating droplets by condensation has recently been emphasized mainly in connection with the possible wash-out mechanism in the atmosphere. A simple theory is deduced for the collection efficiency of insoluble particles on droplets under theStefan's force and a comparision is made with the binding of aerosol particles byBrown's motion and microturbulence of air flow. In general we must take in account the influence ofStefan's flow by calculating the wash-out efficiency.     

Inversion of source process and related studies of the Taiwan Strait earthquake using genetic algorithm

Applying genetic algorithm to inversion of seismic moment tensor solution and using the data of P waveform from digital network and initial motion directions of P waves of Taiwan network stations, we studied the moment tensor solutions and focal parameters of the earthquake of M=7.3 on 16 September of 1994 in Taiwan Strait and other four quakes of M _L≥5.8 in the near region (21°–26°N, 115°–120°E). Among the five earthquakes, the quake of M=7.3 on September 16, 1994 in Taiwan Strait is the strongest one in the southeastern coast area since Nan’ao earthquake of M=7.3 in 1918. The results show that moment tensor solution of M=7.3 earthquake is mainly double-couple component, and is normal fault whose fault plane is near NW. The strike of the fault plane resembles that of the distributive bands of earthquakes before the main event and fracture pattern shown by aftershocks. The tension stress axis of focal mechanism is about horizontal, near in NE strike, the compressive stress axis is approximately vertical, near in NWW strike. It seems that this quake is controlled by the force of Philippine plate’s pressing Eurasian plate in NW direction. But from the viewpoint of P axis of near vertical and T axis of near horizontal, it is a normal fault of strong tensibility. There are relatively big difference between focal mechanism solution of this quake and those of the four other strong quakes. The complexity of source mechanism solution of these quakes represents the complexity of the process of the strait earthquake sequences. Contribution No. 98A01001, Institute of Geophysics, State Seismological Bureau, China. The subject is supported and helped by Academician Yun-Tai CHEN, Profs. Qing-Yao HONG, Zhen-Xing YAO, Tian-Yu ZHENG, Yao-Lin SHI, Ji-An XU, Bo-Shou HUANG and colleague Mei-Jian AN, Xue-Reng DING, Rui-Feng LIU. De-Chong ZHANG and Ming Li provided the digital data warm-heartedly. Lin-Ying WANG offered us the catalogue of earthquakes in southeastern coastal area in China. Xi-Li WANG and Tong-Xia BAI provided us the issued annual reports data. The authors would like to express their gratitude to all of these people. This paper is sponsored by the National Natural Science Foundation of China and Scientific and Technological Commission of Shantou, Guangdong Province.     

新疆伽师强震群的震源破裂特征

利用国内外短周期P波初动符号及全球数字地震台网 (GDSN)的宽频带数字记录资料 ,研究了 1997年 1月 2 3日至 1997年 11月 4日伽师强震群的震源机制解和伽师强震群中的 5个强震、后续地震及其周边地震的震源破裂过程。从震源机制解来看 ,伽师强震群主要有走滑和正倾两种破裂类型 ,主压应力轴方向主要为NNE向和近垂直 ,而主张应力轴为NW向并近水平 ,与区域构造应力场方向存在差异 ,具有明显的局部特征。从破裂过程看 ,伽师强震群的破裂过程相对简单 ,破裂面积不大 ,上升时间较短 ,是由一点向四周快速扩散的脆性破裂 ,无明显的伸展方向 ,与阿图什地震完全不同。研究结果表明 ,伽师强震群与震源附近的地壳结构在垂向和横向上的非均匀变化有着密切关系 ,而阿图什地震与塔里木盆地的现今构造运动关系密切。伽师强震群是在震源区附近地壳上部垂向和水平向力共同作用下发生的多次沿NNE向的快速脆性破裂 ,从而形成了以张性破裂和左旋走滑为主的震源特征     

利用CAP方法和瑞利面波振幅谱联合反演宁强5.3级地震震源深度

2018年9月12日19时6分,陕西省汉中市宁强县发生5.3级地震,不同机构给出的震源深度结果相差较大。为进一步确定宁强5.3级地震的震源深度,基于区域速度模型,首先利用CAP方法反演得到该地震的震源机制解,然后采用瑞利面波振幅谱和CAP深度误差函数联合反演,进一步测定了此次地震的矩心深度。结果显示:CAP方法得到的陕西宁强5.3级地震矩心深度约为12km,瑞利面波振幅谱测定的矩心深度为13km,结合引入的误差函数联合反演,最终确定陕西宁强5.3级地震的矩心深度为13km左右,表明此次地震仍属于发生于上地壳的地震。     

Wave propagation in two layered medium under initial stresses

Summary In this paper, the frequency equation for phase velocity of waves propagated in a laminated medium consisting of two eleastic layers of finite thickness under initial stresses, has been obtained. It has been shown that when wave length becomes very small compared to the thickness of each layer, the wave approaches two Rayleigh waves at the two outer surfaces with the possibility of Stoneley waves at the interface. The propagation ofSH-waves in the composite medium under initial stresses has also been discussed. A particular case has been taken to find the velocity of Love wave in the homogeneous half space under initial compressive stresses.Biot's incremental deformation theory has been used.     

The reservoir-associated earthquakes of April 1983 in western Thailand: Source modeling and implications for induced seismicity

On 22 April 1983, a very large area of Thailand and part of Burma were strongly shaken by a rare earthquake (m _b=5.8,M _s=5.9). The epicenter was located at the Srinagarind reservoir about 190 km northwest of Bangkok, a relatively stable continental region that experienced little previous seismicity. The main shock was preceded by some foreshocks and followed by numerous aftershocks. The largest foreshock ofm _b=5.2 occurred 1 week before the main shock, and the largest aftershock ofm _b=5.3 took place about 3 hours after the main shock. Focal mechanisms of the three largest events in this earthquake sequence have been studied by other seismologists using firts-motion data. However, the solutions for the main shock and the largest aftershock showed significant inconsistency with known surface geology and regional tectonics. We reexamined the mechanisms of these three events by using teleseismicP-andS-waveforms and through careful readings ofP-wave first motions. The directions of theP axes in our study range from NNW-SSE to NNE-SSW, and nodal planes strike in the NW-SE to about E-W in agreement with regional tectonics and surface geology. The main shock mechanism strikes 255°, dips 48°, and slips 63.5°. The fault motions during the main shock and the foreshock are mainly thrust, while the largest aftershock has a large strike-slip component. The seismic moment and the stress drop of the mainshock are determined to be 3.86×10²⁴ dyne-cm and 180 bars, respectively. The occurrence of these thrust events appears to correlate with the unloading of the Srinagarind reservoir. The focal depths of the largest foreshock, the main shock, and the largest aftershock are determined to be 5.4 km, 8 km, and 22.7 km, respectively, from waveform modeling and relative location showing a downward migration of hypocenters of the three largest events during the earthquake sequence. Other characteristics of this reservoir-induced earthquake sequence are also discussed.     

CodaQ as a combination of scattering and intrinsic attenuation: Numerical simulations with the boundary integral method

Numerical modelling ofSH wave seismograms in media whose material properties are prescribed by a random distribution of many perfectly elastic cavities and by intrinsic absorption of seismic energy (anelasticity) demonstrates that the main characteristics of the coda waves, namely amplitude decay and duration, are well described by singly scattered waves in anelastic media rather than by multiply scattered waves in either elastic or anelastic media. We use the Boundary Integral scheme developed byBenites et al. (1992) to compute the complete wave field and measure the values of the direct waveQ and coda wavesQ in a wide range of frequencies, determining the spatial decay of the direct wave log-amplitude relation and the temporal decay of the coda envelope, respectively. The effects of both intrinsic absorption and pure scattering on the overall attenuation can be quantified separately by computing theQ values for corresponding models with (anelastic) and without (elastic) absorption. For the models considered in this study, the values of codaQ ^–1 in anelastic media are in good agreement with the sum of the corresponding scatteringQ ^–1 and intrinsicQ ^–1 values, as established by the single-scattering model ofAki andChouet (1975). Also, for the same random model with intrinsic absorption it appears that the singly scattered waves propagate without significant loss of energy as compared with the multiply scattered waves, which are strongly affected by absorption, suggesting its dominant role in the attenuation of coda waves.     

Propagation ofSH-waves in a heterogeneous visco-elastic layer overlying a homogeneous voigt medium

Summary The propagation ofSH-waves in a vosco-elastic layer having variable material properties has been studied to explain the damping phenomenon of the wave motion. The damping coefficient may be calculated and the effect of heterogeneity on such type of wave motion may be observed from the obtained frequency equation.     

Diffraction ofSH-waves originating from a moving point source by a rigid quarter space

Summary ApplyingWiener-Hopf techniques the diffraction ofSH-waves originating from a moving point source due to the presence of a step change in the elevation of the surface of a homogeneous half-space have been obtained. An interesting feature of resonance associated with the motion of a load at a velocityv ₀ ^–1 approachingv ₀ ^–1 has been observed.     

Estimation of crustal structure in central and west Japan from waveform modeling of regionalPL waves

An inverse method of modeling the regionalPL waveform with the predominant period of about 20 s was developed to estimate the averageS-velocity structure of the upper crust. Applicability of the waveform modeling was confirmed by the results of the numerical experiments: thePL waveform is most sensitive to theS velocity in the upper crust, whereas it is not affected significantly by errors involved in the focal mechanism solution and focal depth determination when thePL wave is well developed. The method was applied to the observed seismograms recorded in central Japan from the earthquakes with epicentral distances 300–500 km. As a result, distinct regional differences were found in the upper crustalS velocity; in particular, between the southern Shikoku district, west Japan, and the southern Chubu district, central Japan, and between the mountainous and the coastal areas in the southern Chubu district. These differences are in agreement with the general features ofP-velocity structures obtained by explosion experiments and by analyses of natural earthquakes. Our method is effective to the extent that the crustal structure along the propagation path can be assumed a horizontally layered structure; it is not applicable when the sensitivity of thePL waveform to the error in the focal mechanism solution is exceptionally high.     

用P波、S波初动和振幅比计算新疆伽师两次强地震震源机制解

针对哈佛大学与USGS给出的1998年8月2日新疆伽师6.1级地震的震源机制解的不同,利用Snoke开发的利用P、S波初动和P、S波的振幅比联合计算震源机制解的程序对这次强震的震源机制解进行了重新计算,同时也计算了8月27日6.4级地震的震源机制解。通过与哈佛大学和USGS的结果对比以及余震分布的分析,表明结果是合理的,同时讨论了这种方法的实用性。     

2021年阿克塞M5.5地震重定位及构造意义

2021年8月26日甘肃省酒泉市阿克塞县发生M5.5地震,这次地震是发生在祁连山地震构造带西段的一次显著地震。利用区域台网记录的宽频带地震波形数据,通过CAP方法反演阿克塞M5.5主震及其M3.7余震的震源机制解,进一步利用双差定位方法对研究区2021年1月1日—11月29日间的地震事件进行重定位,以此分析此次地震的发震构造及其意义。结果表明阿克塞M5.5地震是一次以逆冲性质为主的地震事件;重定位后,地震事件呈明显的集中分布特征,阿克塞M5.5主震震源深度为14.1 km,余震序列震源深度大多分布于15～25 km。综合分析震源机制解、重定位结果以及区域构造背景,认为阿克塞地震的发震断层为党河南山南缘断裂,未来应重点关注祁连山西段发生中强地震的危险。