近断层速度脉冲地震动引起的边坡永久位移预测模型

基于小波分析方法,从NGA数据库的3 551条地震记录中选取189条速度脉冲地震动,地震动均转换成发生最强脉冲的方向。基于Newmark方法,分析了近断层速度脉冲地震动作用引起的边坡永久位移值。结果表明：近断层速度脉冲地震动对边坡产生特殊的破坏作用,表现在滑动位移值大、滑动体破坏力强等方面;边坡永久位移值与速度脉冲地震动的峰值速度具有高度相关性,位移值较大时尤为明显。建立了基于单变量形式的峰值速度及双变量形式的峰值速度、峰值加速度两种边坡永久位移预测模型,模型简单实用,与回归数据具有很好的相关性,前者更适用于预测对实际工程影响较大的永久位移值,且离散性较小。提出的预测模型为考虑近断层地震动速度脉冲特性影响的边坡永久位移值的概率地震灾害分析提供了基础。     

Near-fault strong motion complexity of the 2000 Tottori earthquake (Japan) from a broadband source asperity model

The October 6/2000 Tottori earthquake that occurred in central Japan was an intermediate size strike-slip event that produced a very large number of near field strong motion recordings. The large amount of recorded data provides a unique opportunity for investigating a source asperity model of the Tottori earthquake that, combined with a hybrid strong motion simulation technique, is able to reproduce the observed broadband frequency near-fault ground motion.

We investigated the optimum source asperity parameters of the Tottori earthquake, by applying a Genetic Algorithm (GA) inversion scheme to optimise the fitting between simulated and observed response spectra and Peak Ground Acceleration (PGA) values. We constrained the initial model of our inversion by using the heterogeneous slip distribution obtained from a kinematic inversion of the source of previous studies. We used all the observed near-fault ground motions (−100 m) from the borehole strong motion network of Japan (KiK-Net), which are little affected by surficial geology (site effects).

The calculation of broadband frequency strong ground motion (0.1–10 Hz) is achieved by applying a hybrid technique that combines a deterministic simulation of the wave propagation for the low frequencies and a semi-stochastic modelling approach for the high frequencies. For the simulation of the high frequencies, we introduce a frequency-dependent radiation pattern model that efficiently removes the dependence of the pattern coefficient on the azimuth and take-off angle as the frequency increases. The good agreement between the observed and simulated broadband ground motions shows that our inversion procedure is successful in estimating the optimum asperity parameters of the Tottori earthquake and provides a good test for the strong ground motion simulation technique.

The ratio of background stress drop to average asperity stress drop from our inversion is nearly 50%, in agreement with the theoretical asperity model of Das and Kostrov [Das, S., Kostrov, B.V., 1986. Fracture of a single asperity on a finite fault: a model for weak earthquakes? Earthquake Source Mechanics, AGU, pp. 91–96.], and an empirical ratio of asperities to rupture area [Seismol. Res. Lett. 70 (1999) 59–80.].

The simulated radiation pattern is very complex for epicentral distances within half the fault length, but it approaches the radiation of a double-couple point source for larger distances.

The rupture velocity and rise time have a significant influence on the Peak Ground Velocity (PGV) distribution around the fault. An increase in rupture velocity produces a similar effect on the ground motion as a reduction in rise time.     

Seismic performance evaluation of dam-reservoir-foundation systems to near-fault ground motions

Ground motion records obtained in recent major strong earthquakes have provided evidence that ground motions recorded near the near-fault regions differ in many cases from those observed further away from the seismic source. As the forward directivity and fling effect characteristics of the near-fault ground motions, they have the potential to cause more considerable damage to structures during an earthquake. Therefore, understanding the influence of near-fault ground motions on the performance of structures is critical to mitigate damage and perform effective response. This paper presents results of a study aimed at evaluating the effects of near-fault and far-fault ground motions on seismic performance of concrete gravity dams including dam-reservoir-foundation interaction. Koyna gravity dam is selected as a numerical application. Four different near-fault ground motion records with an apparent velocity pulse are used in the analyses. The earthquake ground motions recorded at the same site from other events that the epicenter far away from the site are employed as the far-fault ground motions. The seismic performance evaluation method based on the demand-capacity ratio, the cumulative overstress duration and the spatial extent of overstressed regions is presented. The concrete damaged plasticity model including the strain hardening or softening behavior is employed in nonlinear analyses. Nonlinear seismic damage analyses of the selected concrete dam subjected to both near-fault and far-fault ground motions are performed. The results obtained from the analyses show the effects of near-fault ground motions on seismic performance of concrete gravity dams and demonstrate the importance of considering the near-fault ground excitations.     

Destructive near-fault strong ground motion from the 2016 Kumamoto prefecture,Japan, M7.3 earthquake

The sequence of the 2016 Kumamoto, Japan, earthquake, which included an initial M6.5 foreshock on April 14, followed by a larger M7.3 mainshock on April 16, and subsequently occurred high aftershock activity, caused significant damage in Kumamoto and neighboring regions. The near-field strong motion record by strong motion network (K-NET and KiK-net) and the intensity meter network demonstrated clearly the characteristics of the strong ground motion developed by the shallow (H = 12 km), inland earthquake comprising short-time duration (<15–20 s) but large (>1G) ground accelerations. The velocity response spectra of the near-fault motion at Mashiki and Nishihara showed large levels (>300–550 cm/s) in the short-period range (T = 1–2 s), several times larger than that of the near-field record of the destructive 1995 Kobe earthquake (M7.3) and that of the 2004 Mid-Niigata earthquake (M6.8). This period corresponds to the collapse vulnerability of Japanese wooden-frame houses, and is the major cause of severe damage during the Kumamoto earthquake. The response spectra also showed extremely large levels (>240–340 cm/s) in the long-period (T > 3 s) band, which is potentially disastrous for high-rise buildings, large oil storage tanks, etc. to have longer resonant period. Such long-period motion was, for the most parts, developed by the static displacement of the fault movement rather than by the seismic waves radiating from the source fault. Thus, the extreme near-fault long-period motion was hazardous only close to the fault but it attenuated very rapidly away from the fault.     

Near-field probabilistic seismic hazard analysis with characteristic earthquake effects

In conventional seismic hazard analysis, uniform distribution over area and magnitude range is assumed for the evaluation of source seismicity which is not able to capture peculiar characteristic of near-fault ground motion well. For near-field hazard analysis, two important factors need to be considered: (1) rupture directivity effects and (2) occurrence of scenario characteristic ruptures in the nearby sources. This study proposed a simple framework to consider these two effects by modifying the predictions from the conventional ground motion model based on pulse occurrence probability and adjustment of the magnitude frequency distribution to account for the rupture characteristic of the fault. The results of proposed approach are compared with those of deterministic and probabilistic seismic hazard analyses. The results indicate that characteristic earthquake and directivity consideration both have significant effects on seismic hazard analysis estimates. The implemented approach leads to results close to deterministic seismic hazard analysis in the short period ranges (T < 1.0 s) and follows probabilistic seismic hazard analysis results in the long period ranges (T > 1.0 s). Finally, seismic hazard maps based on the proposed method could be developed and compared with other methods.     

近断层强地震动预测中的有限断层震源模型

提出了近断层强地震动预测中建立活断层上设定地震有限断层震源模型的方法和步骤.首先,根据地震地质和地震活动性调查以及地球物理勘探等资料,确定活断层的空间方位和滑动类型; 然后,根据地震定标律确定活断层的宏观震源参数; 第三,将高强体模型与k平方滑动模型相结合,产生断层破裂面上的混合滑动分布.在此基础上,预测了与1994年Northridge地震断层类型、矩震级(M_w6.7)基本一致的设定地震的有限断层震源模型.最后,将预测的有限断层震源模型与基于地震学的、使用动力学拐角频率的地震动随机合成方法相结合,预测了1994年Northridge地震近断层12个基岩台站的加速度时程,并和实际记录进行了对比.结果表明,用上述方法和步骤建立的有限断层震源模型是可行、实用的.      

Source complexity of the 4 March 2010 Jiashian,Taiwan, Earthquake determined by joint inversion of teleseismic and near field data

The Jiashian earthquake (M_L 6.4) occurred on 4 March 2010. It was the largest inland event in southern Taiwan of 2010. The mainshock location was unexpected since it occurred in an area with relatively low background seismicity. In addition, reports of earthquake focal mechanisms do not fit with any known active fault geometry. In order to understand the origin of this earthquake, especially its rupture process, we perform a joint source inversion by using teleseismic body wave, GPS coseismic displacements and near field ground motion data. In this study, we considered a northwest–southeast trending fault with a northeast dip retrieved from GPS coseismic data and aftershocks distribution. To analyze the detailed slip distribution in space and time, we used near field 3D Green’s functions provided by spectral-element method and a full time–space inversion technique. We find a complex rupture process with several slip patches distributed inside two main asperities. The slip map reveals a mean slip of 12.9 cm for a maximum slip of 27.3 cm leading to a M_w 6.47 for this event. The rupture initiates in the deepest portion of the fault at 20 km depth, and propagated upward up to 2 km depth to form the two asperities. The source time function of this event revealed two pulses corresponding to the two asperities, for a total duration time of about 16 s. Most aftershocks occurred near the upper boundary of the deepest asperity while no aftershocks are located close to the shallowest one. We infer that the locations of these slip patches are related to the surrounding fault systems that may have restricted the rupture propagation during the earthquake.     

基于实验结果讨论断层破裂与强震物理过程的若干问题

基于断层摩擦滑动实验、含凹凸体断层的变形破坏实验、断层撕裂扩展的实验、交叉断层的变形实验等多种实验结果并结合前人的工作 ,讨论了与断层破裂与强震物理过程相关的若干问题。研究表明 ,断层的整体滑动引起其两侧块体弹性应变的释放 ,是强震发生的原因 ,因此构造活动区具有较大尺度、结构连续且简单、介质均匀的断层 (或断层段 )是产生强震必备的构造条件 ,深部新生断层(盲断层 )向上撕裂扩展产生强震 ,尚需“弱层”提供“解耦”条件以便断层发生整体滑动。强震孕育过程中包含着凹凸体的破裂 ,断层面上凹凸体的尺度、强度及数量决定着前震活动的特征、强震动态破裂过程以及前兆现象。由断层分割的块体通过边界断层的交替滑动、以“框动”的方式运动 ,因此块体周边的断层上强震活动具有交替性。     

Fault architecture and deformation mechanisms in exhumed analogues of seismogenic carbonate-bearing thrusts

Faults in carbonates are well known sources of upper crustal seismicity throughout the world. In the outer sector of the Northern Apennines, ancient carbonate-bearing thrusts are exposed at the surface and represent analogues of structures generating seismicity at depth. We describe the geometry, internal structure and deformation mechanisms of three large-displacement thrusts from the km scale to the microscale. Fault architecture and deformation mechanisms are all influenced by the lithology of faulted rocks. Where thrusts cut across bedded or marly limestones, fault zones are thick (tens of metres) and display foliated rocks (S-CC′ tectonites and/or YPR cataclasites) characterized by intense pressure-solution deformation. In massive limestones, faulting occurs in localized, narrow zones that exhibit abundant brittle deformation. A general model for a heterogeneous, carbonate-bearing thrust is proposed and discussed. Fault structure, affected by stratigraphic heterogeneity and inherited structures, influences the location of geometrical asperities and fault strain rates. The presence of clay minerals and the strain rate experienced by fault rocks modulate the shifting from cataclasis-dominated towards pressure-solution-dominated deformation. Resulting structural heterogeneity of these faults may mirror their mechanical and seismic behaviour: we suggest that seismic asperities are located at the boundaries of massive limestones in narrow zones of localized slip whereas weak shear zones constitute slowly slipping portions of the fault, reflecting other types of “aseismic” behaviour.     

A dislocation model of aseismic fault slip under nonuniform friction

A model is proposed for studying the mechanical behaviour of faults during their interseismic periods. The model considers a plane fault surface in an elastic medium, subject to a uniform shear stress which increases slowly with time. A1-D friction distribution is assumed on the fault, characterized by asperities and a weaker zone. The traction vector on the fault plane has an arbitrary orientation: in particular, it can be nonperpendicular to the asperity borders. Aseismic fault slip takes place when the applied stress exceeds the frictional resistance: slip starts in weak zones and is confined by asperities, where it propagates at increasing velocity. Propagation into asperities is characterized by a dislocation front, advancing perpendicularly to the asperity border. Fault slip does not take prate in the direction of traction, except when traction is perpendicular or parallel to the asperity border. The propagation of such aseismic dislocations produces a stress redistribution along the fault and can play a key role in determining the conditions which give rise to earthquakes.     

The 27 August 2010 Mw 5.7 Kuh-Zar earthquake (Iran): field investigation and strong-motion evidence

In this study, seismological aspects and field observation of the 2010 Kuh-Zar earthquake has been investigated. The Kuh-Zar earthquake, of magnitude 5.7 (Mw), occurred in northeastern Iran on August 27, 2010. The area is surrounded by branches of the active faults which are coated by the quaternary alluvium. During the past several decades, this area has been struck by a number of earthquakes. This earthquake with a moderate magnitude caused a higher rate of damage contrasted with previous earthquakes of the same magnitude range in Iran. Fortunately, the source of the Kuh-Zar earthquake was in a sparsely populated area, and therefore, it caused a few loss of life with the highest observed intensity of shaking VII (modified Mercalli intensity) in the Kuh-Zar village. The shock killed 4 people, injured 40, and destroyed more than 12 villages. According to the field observation, the mechanism of this shock is defined as a left-lateral strike slip. We also checked out the properties of strong ground motions in this earthquake using the records availed by Iranian strong motion network. At KUZ station, about 7 km east of the epicenter, the recorded PGA and PGV in both horizontal and vertical components were remarkably large for an event of this size, and visual inspection of the velocity time history reveals a pulse-like shape. Unfortunately no other recording stations were located close enough to the fault to capture a directivity pulse. Finally, according to the strong-motion properties and observed information, ShakeMaps of the earthquake have been generated by the native intensity observations and the recorded strong motions.     

Near-field and far-field simulation of accelerograms of Sikkim earthquake of September 18, 2011 using modified semi-empirical approach

The semi-empirical approach for modeling of strong ground motion given by Midorikawa (Tectonophysics 218:287?C295, 1993) has been modified in the present paper for component wise simulation of strong ground motion. The modified approach uses seismic moment in place of attenuation relation for scaling of acceleration envelope. Various strong motion properties like directivity effect and dependence of peak ground acceleration with respect to surface projection of source model have been studied in detail in the present work. Recently, Sikkim earthquake of magnitude 6.9 (M _w ) that occurred on September 18, 2011 has been recorded at various near-field and far-field strong motion stations. The modified semi-empirical technique has been used to confirm the location and parameters of rupture responsible for this earthquake. Strong motion record obtained from the iterative modeling of the rupture plane has been compared with available strong motion records from near as well as far-field stations in terms of root mean square error between observed and simulated records. Several possibilities of nucleation point, rupture velocity, and dip of rupture plane have been considered in the present work and records have been simulated at near-field stations. Final selection of model parameters is based on root mean square error of waveform comparison. Final model confirms southward propagating rupture. Simulations at three near-field and twelve far-field stations have been made using final model. Comparison of simulated and observed record has been made in terms of peak ground acceleration and response spectra at 5?% damping. Comparison of simulated and observed record suggests that the method is capable of simulating record which bears realistic appearance in terms of shape and strong motion parameters. Present work shows that this technique gives records which matches in a wide frequency range for Sikkim earthquake and that too from simple and easily accessible parameters of the rupture plane.     

Tikhonov's regularization for deconvolution in the empirical Green function method and vertical directivity effect

Application of Tikhonov's technique, using input errors for the parameter of regularization estimation, enhances the accuracy and stability of the reconstruction of a source time function (STF) by the empirical Green function (EGF) method that gives us an opportunity to use simultaneously for analysis body and surface waves data, and to estimate the horizontal and vertical directivity effects. Knowledge of the last is particularly useful for the choice of an active nodal plane of earthquakes with the dip slip fault orientation that allows us to classify these earthquakes to the interplate or intraplate types and thereby to reach the better understanding of tectonic processes in the region of interest.By way of illustration, an attempt to estimate average parameters of faulting in a first approximation is made herein for two Russian Far East large events with opposite types of focal mechanism orientation, strike slip and dip slip. The former is not a matter of interest in the context of vertical directivity effect but enables us to test the method.The directivity analysis of pulse durations and inverse amplitudes of the relative source time functions (RSTFs) restored at eight globally distributed stations IRIS indicates that the destruction in the source of the Neftegorsk earthquake (05/27/1995 M_W=7.1) propagated roughly horizontally in the direction 8±11° during 19.2±0.4 s along the rupture extending 35.5±4.9 km. The calculated slip distribution along the rupture coincides within the error with the results of field geological measurements on the causal surface fault that proves that the Neftegorsk earthquake source is well described by the model of the linear unilateral fault and gives a good assessment of the method applied.The average parameters of faulting in the Kamchatka earthquake (03/08/1999 M_W=6.9) have been determined from data of 13 station IRIS. It was shown that the destruction in its source propagated downward at an angle of about 60° with horizon, in the direction about S156° E, during 13.4±0.2 s, along the rupture totaling 25.5±2.3 km in length. Therefore, the nodal plane, steeply dipped to the SE, was active and this event can be regarded as an intraplate type. Two asperities can be selected; the first with the maximum slip 3.3 m located at a distance of about 7 km from the onset of rupture, and the second with the maximum slip about 0.9 m centered at approximately 19 km from that.     

Asperity distribution of the 1964 Great Alaska earthquake and its relation to subsequent seismicity in the region

The 1964 Alaska earthquake was the second largest seismic events in the 20th century. The aim of this work is the use of surface deformation data to determine asperity and slip distributions on the fault plane of the Alaska earthquake: these distributions are calculated by a Monte Carlo method. To this aim, we decompose the fault plane in a large number of small square asperity units with a side of 25 km; this allows us to obtain plane surfaces with an irregular shape. In the first stage, each asperity unit is allowed to slip a constant amount or not to slip at all, providing the geometry of the dislocation surface that best reproduces the observed displacements. To this purpose, a large number of slip distributions have been tried by the use of the Monte Carlo method. The slip amplitude is the same for all the asperities and is equal to the average fault slip inferred from the seismic moment. In the second stage, we evaluate the slip distribution in the dislocation area determined by the Monte Carlo inversion: in this case, we allow unit cells to undergo different values of slip in order to refine the initial dislocation model. The results confirm the previous finding that the slip distribution of the great Alaska earthquake was essentially made of two dislocation areas with a higher slip, the Prince William Sound and the Kodiak asperities. Analysis of the post-1964 seismicity in the rupture region shows a strong correlation between the larger earthquakes (M_w≥6) and the distribution of locked asperities following the 1964 event, which can be considered as an independent test of the validity of the model. We do not find slip values higher than 25 m for any of the patches, and we determine two separate high-slip zones: one correspondent to the Prince William Sound asperity, and one (18 m slip) to the Kodiak asperity. The slip distribution connected with the 1964 shock appears to be consistent with the following seismicity in the region.     

Correlation between ground motion parameters and lining damage indices for mountain tunnels

In order to study the interdependencies between the overall lining damage indices for mountain tunnels and ground motion parameters, a total of 89 ground motion records were selected for analysis. These records were divided into three groups: (1) near-field ground motion with velocity pulses, (2) near-field ground motion without velocity pulses, and (3) far-field ground motion. Calculations of Pearson’s linear correlation coefficient have been used to clarify the grade of interrelationships between ground motion parameters. Further, nonlinear dynamic analyses were carried out on mountain tunnels to determine overall lining damage indices under seismic actions. Based on numerical results, linear correlation coefficients between ground motion parameters and overall lining damage indices for mountain tunnels were then calculated and analyzed. Overall lining damage indices were found to be highly correlated to the velocity-related seismic parameters but not so well correlated to spectra-related parameters.     

Scaling laws between seismo-electric/magnetic fields and earthquake magnitude

The mechanisms generating the charges/potentials and abnormal electromagnetic radiation that accompany earthquakes have been discussed based on observational, theoretical and laboratory studies, in particular over the last few decades. We previously conducted stick‐slip tests using precut granite cores with a contact electrode on the sample side, and proposed that surface charges released from charge trapping centres on sheared asperities are a possible source of the above‐mentioned charges or potentials. Here we discuss seismo‐electric and magnetic fields formed by the surface charges, and derive their scaling laws with the earthquake magnitude at a distance. According to our model, detectability of these fields at a point on the ground surface depends strongly on wet/dry condition at a fault zone (sheared asperities).     

Estimation of causative fault parameters of the Rudbar earthquake of June 20, 1990 from near field SH-wave data

We analyze the strong motion accelerograms recorded for the large (M_S=7.7, M_W=7.3, m_b=6.4) Rudbar earthquake of June 20, 1990. The earthquake had a complex source process. We have identified the imprints of rupture of three localized asperities on the major causative fault on the accelerograms. These asperities are interpreted to correspond to (i) the main shock that initiated the rupture process and was located in the domino block between the Kabateh and Zard Goli faults, (ii) a foreshock that occurred about 10 s earlier in the Kabateh fault and (iii) a later shock, on the western end of the Baklor fault, which terminated the bilateral rupture process at the western end. We estimate the strike, dip and slip of these causative sub-event rupture planes using the SH spectral amplitudes, based on a point source representation of sub-events and a non-linear least square formulation for inversion of the amplitude data. The results of our inversion of the near field data are comparable to other studies based on teleseismic data.     

Fault rupture directivity of large earthquakes in Himalaya

The paper examines the predominant fault rupture directivity during large earthquakes in different sectors of the Himalaya which influences strong ground motion and damage scenario. The nature of the faulting of earthquakes vis-à-vis their rupture directivity has been discussed. It is found that the rupture directivity near the Indo-Eurasian plate boundary varies from place to place i.e. either along the strike direction of the faults or at right angles to it. The secondary meizoseismal areas as observed for 1505 Dharchula, 1803 Uttarakhand, 1905 Kangra earthquakes in the Himalaya and 2001 Bhuj earthquake in stable continental region suggest that they are a fairly good indicator of predominant rupture directivity since the latter accentuates the site response up to a longer distance. The resulting larger ground motions, therefore, need to be incorporated in the design of engineering structures by suitable modifications in the BIS code.     

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.     

Seismic demands and analysis of site effects in the Marmara region during the 1999 Kocaeli earthquake

The characteristics of the strong ground motion accelerograms from the 1999 Kocaeli earthquake are investigated in detail in this study. The emphasis is on the comparison of the response spectra for the fault normal (FN) and fault parallel (FP) components of the ground motions. The results show that the near-fault records with directivity effects characterize themselves with increased base shear demands rather than increased displacement demands for both the FN and FP components and a narrower velocity sensitive region for the FN component. This study also shows that the effectiveness of base isolation may vary from site to site and for a given site, from component to component. The site effects in the Marmara region during the 1999 Kocaeli earthquake are examined. Site amplifications are predicted by the classical spectral ratio (CSR) and the receiver function (RF) methods. The CSR method gives higher estimates for the site amplifications compared to the RF method and is in better conformity with the observed damage during the Kocaeli earthquake. The districts of Istanbul that are especially susceptible to site amplification hazard are determined. It is apparent from the results that the site amplification hazard risk is the highest for Avcilar and Bakirkoy districts. This study also shows that for sites which have the risk of soil amplification for long-period structures, liquefaction may not be beneficial as a natural base isolator, and may result in shifting the eigenperiod of the low- and mid-rise structures to the critical periods with high site amplifications. This may be especially the case for Avcilar and Bakirkoy districts. In Fatih, Bakirkoy, and Cekmece districts, the predominant period of the ground motion is calculated to be very close to the eigenperiods of the typical residential buildings. Therefore, these three districts are expected to experience heavier damages in future earthquakes due to resonance effects.