Earthquake rupture complexity due to dynamic nucleation and interaction of subsidiary faults

We numerically study the dynamic interaction of propagating cracks. It is assumed that propagating cracks can nucleate and drive subsidiary cracks because of shear strain enhancement near the propagating crack tips. The critical strain fracture criterion is assumed in the analysis. Intense interaction is expected to occur among the cracks. All the cracks are assumed to be parallel and antiplane strain deformation is assumed in the computation.In the interaction of two non-coplanar cracks, a strain shadow is formed in the neighborhood of each crack because of the strain release by the introduction of the crack. The growth of each crack is accelerated when the propagating tips of each crack are outside of the strain shadow of the other crack. In general, the crack tips enter the strain shadow, and the crack tips decelerate. The calculation shows that only one of the two cracks can continue to grow, and the other's growth is decelerated and arrested. If we can assume that the suite of cracks interact in a pairwise manner only, then this may suggest that only a limited number of cracks can continue to grow during the final stage of the rupture process. Hence the crack interaction causes complexity in dynamic earthquake faulting. The concepts of barrier and asperity have been employed by many researchers for the interpretation of complex seismic wave data. However, the physical realities of such concepts are obscure. Our calculations show that dynamic crack interactions can produce barriers and asperities in some cases; the crack tip deceleration or arrest due to the interactions among non-coplanar cracks can be interpreted as being due to a barrier. The dynamic coalescence among the coplanar cracks can be regarded as an asperity.Umeda found a localized area that strongly radiates high-frequency seismic waves in the epicentral areas of some large shallow earthquakes. He defined this as an earthquake bright spot. Our analysis implies that only a limited number of cracks continue to grow when many interactive cracks nucleate, and that all other cracks stop extending soon after nucleation. Hence, if the nucleation and termination of several cracks occur in a localized area, it will be observed seismologically as an earthquake bright spot. This is because it is theoretically known that the sudden termination of crack growth and dynamic crack coalescence efficiently emits high-frequency elastic waves.     

Blind prediction of the site effects at Ashigara Valley,Japan, and its comparison with reality

Weak and strong ground motions were numerically predicted for three stations of the Ashigara Valley test site. The prediction was based on the records from a rock-outcrop station, one weak-motion record from a surface-sediments station, and the standard geotechnical model. The data were provided by the Japanese Working Group on the Effects of Surface Geology as a part of an international experiment. The finite-difference method for SH waves in a 2-D linear viscoelastic medium (a causalQ model) was employed.Comparison with the real records shows that at two stations the predictions fit better than at the third one. Strangely, the two better predictions were for stations situated at larger distances from the reference rock station (one station was on the surface, the other in a borehole). The strong ground motion (the peak acceleration of about 200 cm s^–2) was not predicted qualitatively worse than the weak motion (8 cm s^–2). A less sophisticated second prediction (not submitted during the experiment), in which we did not attempt to fit the available weak-motion record at the sedimentary station, agrees with the reality significantly better.     

An integrable case of a rotational motion analogous to that of Lagrange and Poisson for a gyrostat in a Newtonian force field

The scope of the present paper is to provide analytic solutions to the problem of the attitude evolution of a symmetric gyrostat about a fixed point in a central Newtonian force field when the potential function isV ⁽²⁾.We assume that the center of mass and the gyrostatic moment are on the axis of symmetry and that the initial conditions are the following: (t ₀)=₀, (t ₀)=₀, (t ₀)=(t ₀)=₀, ₁(t ₀)=0, ₂(t ₀)=0 and ₃(t ₀)= ₃ ⁰ .The problem is integrated when the third component of the total angular momentum is different from zero (B ₁ 0). There now appear equilibrium solutions that did not exist in the caseB ₁=0, which can be determined in function of the value ofl ₃ ^r (the third component of the gyrostatic momentum).The possible types of solutions (elliptic, trigonometric, stationary) depend upon the nature of the roots of the functiong(u). The solutions for Euler angles are given in terms of functions of the timet. If we cancel the third component of the gyrostatic momentum (l ₃ ^r =0), the obtained solutions are valid for rigid bodies.     

Automatic equipment to monitor groundwater radon

Prototype instrumentation, able to automatically measure groundwater radon content variations, is presented. The equipment is made of stainless steel and has spherical valves with automatic and pneumatic control. The deemanation of the gases from the water is obtained by evacuating a suitable expansion chamber. The instrumentation can make discrete sampling ranging from 1 per hour to 1 per 99 hours. The equipment was tested in the laboratory: the efficiency was measured by means of a²⁶⁶Ra solution. A mean value of (0.65±0.07) count/s/Bq was obtained. A calibration test was carried out by comparing countings from the automatic equipment with those obtained by the standard laboratory cell. Results of an operational check over a period of approximately one year indicate that variations in radon at the calibration site are attributable more to meteorological than to tectonic causes.     

1983年菏泽5.9级地震前后沂水泉氡的变化

沂水泉位于山东省沂水一汤头断裂带上。多年来的观测分析表明,影响沂水泉氡值变化的最大干扰是降雨。由于这种干扰存在着“记忆”滞后影响,所以用一般的相关分析难以将这种影响排除干净。本文对最小二乘法为建立模型的判据,用动态灰箱分析法拟合并推估观测系统的变化。在分析中用了四个状态变量：长趋势成分Ｌ（ｔ）、年周期变化成分Ｐ（ｔ）、降雨影响部分Ｒ（ｔ）和随机变化部分Ｓ（ｔ）。经过分析表明,沂水泉氡的长趋势部分１     

华北典型箕状断陷盆地层序地层学模式及其与油气赋存关系

本文根据层序地层学的基本原理,立足东部箕状断陷盆地的特点,以冀中拗陷和二连盆地为重点研究对象,提炼出陆相断陷盆地层序地层学模式,在缓坡带,以风暴浪基面作为“陆架边角”的参照物,层序地层特征相似于海相被动大陆边缘背景;在陡坡带则截然不同,并以湖面与边界断裂的交点取代“陆架边角”,各体系域所含沉积体系与海相模式不尽相同,在盆地形成初期的层序,低水位体系域发育,湖泊消亡以后,主要发育高水位体系域。此外,本文还探索了层序地层与油气赋存的关系,指出富砂的“盆底扇”是今后勘探开发的重要方向。     

A Step Towards Evaluation of the Seismic Response Reduction Factor in Multistorey Reinforced Concrete Frames

A seismic nonlinear time-history analysis was made for four-, six-, and eight-storey reinforced concrete buildings. These buildings were made as three-dimensional space frame structures with shear walls in both orthogonal directions. They have five bays with 4.8 m spacing each in the horizontal direction, and three bays with 4.2 m spacing each in the transversal direction. The frames were designed according to the Jordanian Seismic Code of practice for Seismic Zones 4, 3, 2, and 1 as proposed for Jordan by several authors. Time-history analysis was made using the El Centro (N-S) earthquake record of May 1940 as an actual earthquake excitation. The response reduction factor (R) that primarily consists of two factors that are the ductility reduction (Rµ) and the overstrength (), is obtained. It has been seen that the seismic zoning has a slight effect on the ductility reduction factor for different buildings, since it ranges from Zone 4 to Zone 1 as 2.37 to 2.52, 1.72 to 1.78, and 1.14 to 1.18 for four-, six-, and eight-storey buildings, respectively. Moreover, it is observed that, for different buildings and different seismic zones, the ductility reduction factor (Rµ) is slightly different from the system ductility factor (µ) especially for higher values of µ (i.e., Rµ µ). The response reduction factor, called overstrength (), was evaluated. The overstrength factor was found to vary with seismic zones (Z) , number of stories, and design gravity loads. However, the dependency on seismic zones was the strongest. The average overstrength of these buildings in Zones 4 and 1 was 2.61 and 6.94, respectively. The overstrength increased as the number of storeys decreased: overstrength of a four-storey building was higher than an eight-storey building by 36% in Zone 4, and 39% in Zone 1. Furthermore, buildings of the three heights had an average overstrength 165.9% higher in Zone 1 than in Zone 4. These observations have a significant implications for the seismic design codes which currently do not take into account the variation of the response reduction factor, R (i.e., ductility reduction factor times overstrength).     

地震短临综合预报方法的研究

本文在《中国震例》资料基础上,详细分析了华北地区１７次中强以上地震前的短临异常,得到华北地区中强以上地震短临异常的三个综合特征,在归纳出短临异常综合特征的基础上,提出了两个判定孕震过程、前兆异常由中期向短期过渡的定量的综合标志。据中期异常的追踪分析和短临异常综合分析相结合的原则,定义并计算了综合预报指标Ｓ值,以此值作为是否发生中强以上地震的判据。本文的重点是在分析短临异常综合特征的基础上,结合实用     

合肥遥测台网地震速报系统

论述了适用于遥测地震台网大震速报的一套计算机处理系统,系统设计以《全国遥测台网观测技术规范》和《近震分析》理论为依据,实现了数据输入,震相选择,数据处理,地震定位,震级计算,结果输出及贮存等多项功能,达到了快速确定地震三要素,提高速报速度和精度之目的,是一个较完善的,实用的速报系统。     

Possible creep-related tilt precursors obtained in the Central Apennines (Italy) and in the Southern Caucasus (Georgia)

Daily averaged tilt component data from two sites of the Central Apennines (Italy) and of the Southern Caucasus (Georgia), respectively, revealed intermediate-term tilts as possible precursors to earthquakes (M=3.0÷4.7) which occurred in the above-mentioned seismic areas within a distance of 50 km from the sites. A good temporal correlation as well as a fair spatial correspondence between these residual tilts (with amplitude and duration of some microradians and months, respectively) and main shocks were pointed out, by removing both secular trends and seasonal thermoelastic effects from the raw tilts. An attempt was made to justify the above-mentioned results, based on the assumption that the observed intermediate-term preseismic tilts are the manifestation of aseismic creep episodes of comparable duration in the fault materials of thrust faults close to the tilt sites. The mechanism refers to a strain field slowly propagating from the preparation focal area to the tilt site, through crustal blocks separated by weak transition zones. This propagation is thought to be the cause of the local aseismic fault slip recorded by the tiltmeters. Previously, both discrete structures and strain propagation effects were revealed in the Central Apennines and are thought also to exist in the Southern Caucasus. As in the past, the rheological properties of fault materials are revealed as viscoelastic ones. In fact, creep equations obtained by applying several viscoelastic models on our data, proved to fit quite well some of the observed tilt precursors, producing viscosity and rigidity values very similar to those reported in literature.Professor Petr Viktorovich Manjgaladze died during the writing of this paper