The Zihuatanejo, Mexico, earthquake of 1994 December 10 (M= 6.6): source characteristics and tectonic implications

An analysis of the Zihuatanejo, Mexico, earthquake of 1994 December 10 ( M = 6.6), based on teleseismic and near-source data, shows that it was a normal-faulting, intermediate-depth ( H = 50 ± 5 km) event. It was located about 30 km inland, within the subducted Cocos plate. The preferred fault plane has an azimuth of 130°, a dip of 79° and a rake of −86°. The rupture consisted of two subevents which were separated in time by about 2 s, with the second subevent occurring downdip of the first. The measured stress drop was relatively high, requiring a Δσ of about a kilobar to explain the high-frequency level of the near-source spectra. A rough estimate of the thickness of the seismogenic part of the oceanic lithosphere below Zihuatanejo, based on the depth and the rupture extent of this event, is 40 km.
This event and the Oaxaca earthquake of 1931 January 15 ( M = 7.8) are the two significant normal-faulting, intermediate-depth shocks whose epicentres are closest to the coast. Both of these earthquakes were preceded by several large to great shallow, low-angle thrust earthquakes, occurring updip. The observations in other subduction zones show just the opposite: normal-faulting events precede, not succeed, updip, thrust shocks. Indeed, the thrust events, soon after their occurrence, are expected to cause compression in the slab, thus inhibiting the occurrence of normal-faulting events. To explain the occurrence of the Zihuatanejo earthquake, we note that the Cocos plate, after an initial shallow-angle subduction, unbends and becomes subhorizontal. In the region of the unbending, the bottom of the slab is in horizontal extension. We speculate that the large updip seismic slip during shallow, low-angle thrust events increases the buckling of the slab, resulting in an incremental tensional stress at the bottom of the slab and causing normal-faulting earthquakes. This explanation may also hold for the 1931 Oaxaca event.     

Upper Mississippi embayment shallow seismic velocities measured in situ

Vertical seismic compressional- and shear-wave (P-and S-wave) profiles were collected from three shallow boreholes in sediment of the upper Mississippi embayment. The site of the 60-m hole at Shelby Forest, Tennessee, is on bluffs forming the eastern edge of the Mississippi alluvial plain. The bluffs are composed of Pleistocene loess, Pliocene-Pleistocene alluvial clay and sand deposits, and Tertiary deltaic-marine sediment. The 36-m hole at Marked Tree, Arkansas, and the 27-m hole at Risco, Missouri, are in Holocene Mississippi river floodplain sand, silt, and gravel deposits. At each site, impulsive P- and S-waves were generated by man-made sources at the surface while a three-component geophone was locked downhole at 0.91-m intervals.

Consistent with their very similar geology, the two floodplain locations have nearly identical S-wave velocity (V_S) profiles. The lowest V_S values are about 130 m s⁻¹, and the highest values are about 300 m s⁻¹ at these sites. The shear-wave velocity profile at Shelby Forest is very similar within the Pleistocene loess (12 m thick); in deeper, older material, V_S exceeds 400 m s⁻¹.

At Marked Tree, and at Risco, the compressional-wave velocity (V_P) values above the water table are as low as about 230 m s⁻¹, and rise to about 1.9 km s⁻¹ below the water table. At Shelby Forest, V_P values in the unsaturated loess are as low as 302 m s⁻¹. V_P values below the water table are about 1.8 km s⁻¹. For the two floodplain sites, the V_P/V_S ratio increases rapidly across the water table depth. For the Shelby Forest site, the largest increase in the V_P/V_S ratio occurs at 20-m depth, the boundary between the Pliocene-Pleistocene clay and sand deposits and the Eocene shallow-marine clay and silt deposits.

Until recently, seismic velocity data for the embayment basin came from eartquake studies, crustal-scale seismic refraction and reflection profiles, sonic logs, and from analysis of dispersed earthquake surface waves. Since 1991, seismic data for shallow sediment obtained from reflection, refraction, crosshole and downhole techniques have been obtained for sites at the northern end of the embayment basin. The present borehole data, however, are measured from sites representative of large areas in the Mississippi embayment. Therefore, they fill a gap in information needed for modeling the response of the embayment to destructive seismic shaking.     

Recent Advances in Prediction and Processing of Strong Ground Motions

We present an overview of our recent results on utilizing small earthquakes in the earthquake engineering practice. Site-specific ground motion time-histories of large earthquakes can be successfully simulated using recordings of small earthquakes which are often referred to as 'empirical Green's functions' in seismology. Another important practical problem is whether and how these observations can be used in seismic risk studies which are based on empirical attenuation relations for ground motion parameters. We study a possibility of extrapolating attenuation relations for small earthquakes, to larger magnitudes using the data from the Garner Valley downhole array in Southern California. Finally we introduce efficient ground motion processing techniques in frequency- and time-domains and apply them to site response estimation.     

Seismic Ground Motion Expected for the Eastern District of Naples

The seismic ground motion of a test area in the eastern district of Naples is computed with a hybrid technique based on the mode summation and the finite difference methods. This technique allows us the realistic modelling of source and propagation effects, including local soil conditions. In the modelling, we consider the 1980 Irpinia earthquake, a good example of strong shaking for the area of Naples, which is located about 90 km from the epicenter.The detailed geological setting is reconstructed from a large number of drillings. The sub-soil is mainly formed by alluvial (ash, stratified sand and peat) and pyroclastic materials overlying a pyroclastic rock (yellow neapolitan tuff), representing the neapolitan bedrock. The detailed information available on mechanical properties of the sub-soil and its geometry warrants the application of the sophisticated hybrid technique.As expected, the sedimentary cover causes an increase of the signal's amplitudes and duration. If thin peat layers are present, the amplification effects are reduced, and the peak ground accelerations are similar to those observed for the bedrock model. This can be explained by the backscattering of wave energy at such layers, that tend to seismically decouple the upper from the lower part of the structure.For SH-waves, the influence of the variations of the S-wave velocities on the spectral amplification is studied, by considering locally measured velocities and values determined from near-by down-hole measurements. The comparison between the computed spectral amplifications confirms the key role of an accurate determination of the seismic velocities of the different layers.The comparison performed between a realistic 2-D seismic response and a standard 1-D response, based on the vertical propagation of waves in a plane layered structure, shows considerable difference, from which it is evident that serious caution must be taken in the modelling of expected ground motion at a specific site.     

水库诱发地震的地震动特征

利用新丰江、美国奥洛维尔和蒙蒂赛洛水库地震极近场地面运动资料，研究了水库诱发地震的地震动峰值加速度衰减。结果表明，水库地震的峰值加速度高，衰减快。利用国内水库地震资料研究了水库地震烈度。水库地震烈度与同级天然地震烈度相比，一般是烈度高，衰减慢。水库地震的地震动反应谱频率高，地震动持续时间短     

地震荷载下莫高窟围岩动态损伤特性研究

本以敦煌莫高窟为研究对象，阐述了地震荷载下洞窟围岩动态损伤的影响因素；采用动力有限元法。从地震动特性入手。分析研究了地震作用对洞窟围岩及其附属构筑物可能造成的损坏，为石窟物地震安全评估及防灾对策研究提供了依据。     

The GIS and analysis of earthquake damage distribution of the 1303 Hongtong M=8 earthquake

The geography information system of the 1303 Hongton M=8 earthquake has been established. Using the spatial analysis function of GIS, the spatial distribution characteristics of damage and isoseismal of the earthquake are studies. By comparing with the standard earthquake intensity attenuation relationship, the abnormal damage distribution of the earthquake is found, so the relationship of the abnormal distribution with tectonics, site condition and basin are analyzed. In this paper, the influence on the ground motion generated by earthquake source and the underground structures near source also are studied. The influence on seismic zonation, anti-earthquake design, earthquake prediction and earthquake emergency responding produced by the abnormal density distribution are discussed. Foundation item: National important fundamental research “The Basic Research of Important Project in Damage Environment” and The important project “The Seismic Hazard Assessment Research and Anti-earthquake Structure Research” from China Earthquake Administration during the 10th Five-year Plan. Contribution No. 04FE1008, Institute of Geophysics, China Earthquake Administration.     

Design Earthquakes in the U.K.

Three sites in the UK are taken, representative of low, medium and high hazard levels (by UK standards). For each site, the hazard value at 10⁻⁴ annual probability is computed using a generic seismic source model, and a variety of ground motion parameters: peak ground acceleration (PGA), spectral acceleration at 10 Hz and 1 Hz, and intensity. Disaggregation is used to determine the nature of the earthquakes most likely to generate these hazard values. It is found (as might be expected) that the populations are quite different according to which ground motion parameter is used. When PGA is used, the result is a rather flat magnitude distribution with a tendency to low magnitude events (\le 4.5 ML) which are probably not really hazardous. Hazard-consistent scenario earthquakes computed using intensity are found to be in the range 5.8–5.9 ML, which is more in accord with the type of earthquake that one expects to be a worst-case event in the UK. This revised version was published online in July 2006 with corrections to the Cover Date.     

On the Incorporation of the Effect of Crustal Structure into Empirical Strong Ground Motion Estimation

This article has two purposes. Firstly, a validation exercise of the modal summation technique for the computation of synthetic strong-motion records is performed for two regions of Europe (Umbria-Marche and south Iceland), using a variety of region specific crustal structure models, by comparing the predicted ground motion amplitudes with observed motions. It is found that the rate of decay of ground motions is well predicted by the theoretical decay curves but that the absolute size of the ground motions is underpredicted by the synthetic time-histories. This is thought to be due to the presence of low-velocity surface layers that amplify the ground motions but are not included in the crustal structure models used to compute the synthetic time-histories. Secondly, a new distance metric based on the computed theoretical decay curves is introduced which should have the ability to model the complex decay of strong ground motions. The ability of this new distance metric to reduce the associated scatter in empirically derived equations for the estimation of strong ground motions is tested. It is found that it does not lead to a reduction in the scatter but this is thought to be due to the use of crustal structure models that are not accurate or detailed enough for the regions studied. This revised version was published online in July 2006 with corrections to the Cover Date.