The 2011 Lorca seismic series: Temporal evolution,faulting parameters and hypocentral relocation

The Lorca 2011 seismic series was recorded by an unprecedented set of high quality on scale broadband seismograms and strong motion accelerograms. The waveforms from permanent and temporary broadband seismic networks deployed in the region by different institutions allowed to invert regional moment tensor for the fore, main and largest aftershock of the complete seismic sequence. Using double-difference algorithm we have performed a precise relocation of the seismic series, where body wave travel times from strong ground motion accelerograms were included in the data set. Regional moment tensor inversion for the three main events show similar oblique-reverse faulting regime with a northeast-southwest fault orientation. The scalar seismic moment, moment magnitude and focal depth retrieved from the inversion yield the following values for each event: \(\hbox {Mo}=6.5\times 10^{16}\) Nm (Mw = 5.2) for the mainshock, \(\hbox {Mo}= 9.6 \times 10^{15}\) Nm (Mw = 4.6) for the foreshock and \(\hbox {Mo}=7.3\times 10^{14}\) Nm (Mw = 3.9) for the large aftershock. The centroid depths range between 4 and 6 km. The double-difference relocation of the seismic series shows significant epicentral differences with the preliminary routine location. The epicentral solutions given by this relocation show a seismic sequence distributed following a NE–SW strike, subparallel to the Alhama de Murcia fault and compatible with the faulting parameters inverted from the moment tensor analysis. The hypocenters of the series generate a subvertical trend in depth distribution, being concentrated between 2 and 6 km. The depth distribution of the main events, which range from 4.6 to 5.5 km, is in good relationship with the faulting and depth parameters deduced from the moment tensor inversion technique. The regional moment tensor solutions for the three largest earthquakes, the epicentral distribution and the focal depths show good relationship with the surface geometry and tectonic regime of the Alhama de Murcia fault. The stress drop deduced for the mainshock gives a value ranging between 58 and 85 bars, which does not support the idea of a high stress drop release as a main factor contributing to the high ground acceleration recorded at Lorca. The PGA values observed at Lorca, which contributed to the high damage independently of structural deficiencies, could be generated mainly by shallowness and proximity to the seismic source together with a directivity effect in the seismic radiation.     

Local site effect microzonation of Lorca town (SE Spain)

Local site effect assessment based on subsurface ground conditions is often the key to evaluate urban seismic hazard. The site effect evaluation in Lorca town (south-eastern Spain) started with a classification of urban geology through the geological mapping at scale 1:10,000 and the use of geotechnical data and geophysical surveys. The 17 geological formations identified were classified into 5 geological/seismic formations according to their seismic amplification capacity obtained from ambient vibration measurements as well as from simultaneous strong motion records. The shear-wave velocity structure of each geological/seismic formation was evaluated by means of inversion of Rayleigh wave dispersion data obtained from vertical-component array records of ambient noise. Nakamura’s method was applied to determine a predominant period distribution map. The spectral amplification factors were fourfold the values recorded in a reference hard-rock site. Finally, the capability of this study for explaining the damage distribution caused by the May 11th, 2011 Lorca destructive earthquake (Mw \(=\) 5.2) was examined. The methods used in this work are of assistance to evaluate ground amplification phenomena in urban areas of complex geology as Lorca town due to future earthquakes with applicability on urban seismic risk management.     

Slope instabilities triggered by the 11th May 2011 Lorca earthquake (Murcia,Spain): comparison to previous hazard assessments and proposition of a new hazard map and probability of failure equation

The Lorca Basin has been the object of recent research aimed at studying the phenomena of earthquake-induced landslides and its assessment in the frame of different seismic scenarios. However, it has not been until the 11th May 2011 Lorca earthquakes when it has been possible to conduct a systematic approach to the problem. In this paper we present an inventory of slope instabilities triggered by the Lorca earthquakes which comprises more than 100 cases, mainly rock and soil falls of small size (1–100  \(\hbox {m}^{3}\) ). The distribution of these instabilities is here compared to two different earthquake-triggered landslide hazard maps: one considering the occurrence of the most probable earthquake for a 475-years return period in the Lorca Basin \((\hbox {M}_{\mathrm{w}}=5.0)\) based on both low- and high-resolution digital elevation model (DEM); and a second one matching the occurrence of the \(\hbox {M}_{\mathrm{w}}=5.2\) 2011 Lorca earthquake, which was performed using the higher resolution DEM. The most frequent Newmark displacements related to the slope failures triggered by the 2011 Lorca earthquakes are lower than 2 cm in both the hazard scenarios considered. Additionally, the predicted Newmark displacements were correlated to the inventory of slope instabilities to develop a probability of failure equation. The fit seems to be very good since most of the mapped slope failures are located on the higher probability areas. The probability of slope failure in the Lorca Basin for a seismic event similar to the \(\hbox {M}_{\mathrm{w}}\) 5.2 2011 Lorca earthquake can be considered as very low (0–4 %).     

Assessment of expected damage on buildings subjected to Lorca earthquake through an energy-based seismic index method and nonlinear dynamic response analyses

The city of Lorca (Spain) was hit on May 11th, 2011, by two consecutive earthquakes of magnitudes 4.6 and 5.2 M \(_\mathrm{w}\) , causing casualties and important damage in buildings. Many of the damaged structures were reinforced concrete frames with wide beams. This study quantifies the expected level of damage on this structural type in the case of the Lorca earthquake by means of a seismic index \(I_{v}\) that compares the energy input by the earthquake with the energy absorption/dissipation capacity of the structure. The prototype frames investigated represent structures designed in two time periods (1994–2002 and 2003–2008), in which the applicable codes were different. The influence of the masonry infill walls and the proneness of the frames to concentrate damage in a given story were further investigated through nonlinear dynamic response analyses. It is found that (1) the seismic index method predicts levels of damage that range from moderate/severe to complete collapse; this prediction is consistent with the observed damage; (2) the presence of masonry infill walls makes the structure very prone to damage concentration and reduces the overall seismic capacity of the building; and (3) a proper hierarchy of strength between beams and columns that guarantees the formation of a strong column-weak beam mechanism (as prescribed by seismic codes), as well as the adoption of counter-measures to avoid the negative interaction between non-structural infill walls and the main frame, would have reduced the level of damage from \(I_{v}\) \(\,=\,\) 1 (collapse) to about \(I_{v}\) \(\,=\,\) 0.5 (moderate/severe damage).     

Changes in dynamic characteristics of Lorca RC buildings from pre- and post-earthquake ambient vibration data

This paper estimates fundamental translational period and damping ratio parameters and examines the changes in dynamic characteristics of a set of low-to-medium rise buildings in Lorca town (SE of Spain) affected by the May 11th, 2011 earthquake. These building parameters have been calculated analysing structural dynamic response from ambient vibration measurements recorded at top RC buildings pre- and post earthquake, using the Fast Fourier Transform and the Randomdec technique. The empirical expression relating period \((T)\) and number of floor \((N)\) here obtained analysing ambient noise recorded on 59 healthy buildings before the earthquake is \(T= (0.054\pm 0.002)\, N\) , very similar to others empirical period–height relationships obtained for RC structures in the European built environment but quite different from code provisions. Measurements performed in 34 damaged buildings show a period elongation after the quake according to \(T^{*} =(0.075\pm 0.002)\,N\) expression. Moreover, we found a rise of the fundamental period with the EMS’s grade of damage of buildings. In contrast to natural frequency, damping ratio \((\xi )\) do not shows a significant variation with earthquake damage degree and the product \(\xi \, T\) remains near constant.     

Characterization of earthquake shaking severity in the town of Lorca during the May 11, 2011 event

Exceptional ground motions have been observed in the near-source region in relatively small and moderate earthquakes such as the Lorca earthquake of 2011. In the Mw 5.2 mainshock, horizontal peak ground acceleration reached 0.37 g at the LOR accelerograph station in the north-central part of the town. Maximum ground acceleration and velocity in this station was strongly polarised in the fault’s strike-normal axis, a feature interpreted as a near-field effect of rupture directivity. This quake hit the Lorca area causing serious structural damage, unevenly distributed in the urban area as result of a heterogeneous spread of building vulnerability and seismic shaking severity. This study investigates the characteristics of the earthquake shaking at the LOR station (on rock) and other 11 zones of Lorca with different site conditions. We have computed the 1-D transfer functions of the shallow structure at these 11 points with available local Vs models obtained with SPAC technique. These transfer functions were used to convolve the LOR accelerogram and simulate the ground motion which occurred at each site. A set of key engineering ground-motion parameters obtained here, indicate different shake levels at each site and all these values are clearly above the LOR station ones. These parameter values are in the range of typical expected values for European earthquakes with local intensity VIII. Their large values, especially those of acceleration response spectra and energy input spectra, more clearly for periods between 0.4 and 0.6 s, point out that the seismic demand has been quite high and help to explain the heavy damage spread on the town’s buildings.     

Surface dynamic deformation estimates from local seismicity: the Itoiz reservoir,Spain

We analyze the ground motion time histories due to the local seismicity near the Itoiz reservoir to estimate the near-source, surface 3D displacement gradients and dynamic deformations. The seismic data were obtained by a semipermanent broadband and accelerometric network located on surface and at underground sites. The dynamic deformation field was calculated by two different methodologies: first, by the seismo-geodetic method using the data from a three-station microarray located close to the dam, and second, by single station estimates of the displacement gradients. The dynamic deformations obtained from both methods were compared and analyzed in the context of the local free-field effects. The shallow 1D velocity structure was estimated from the seismic data by modeling the body wave travel times. Time histories obtained from both methods result quite similar in the time window of body wave arrivals. The strain misfits between methods vary from 1.4 to 35.0 % and rotational misfits vary from 2.5 to 36.0 %. Amplitudes of displacement gradients vary in the range of 10^?8 to 10^?7 strains. From these results, a new scaling analysis by numerical modeling is proposed in order to estimate the peak dynamic deformations for different magnitudes, up to the expected maximum M _w in the region (M5.5). Peak dynamic deformations due to local M _w5.5 earthquakes would reach amplitudes of 10^?5 strain and 10^?3 radians at the Itoiz dam. The single station method shows to be an adequate option for the analysis of local seismicity, where few three-component stations are available. The results obtained here could help to extend the applicability of these methodologies to other sites of engineering interest.     

Rupture process of the Lorca (southeast Spain) 11 May 2011 (M w = 5.1) earthquake

On 11 May 2011, a M _w ?=?5.1 earthquake shook the town of Lorca (SE Spain) causing a disproportionately large damage for its magnitude. In order to contribute to knowledge of the behavior of the active faults present in the region and define the parameters which control their motion, we made a detailed study of the rupture process of this earthquake from inversion of body waves at regional and teleseismic distances. Ground motion displacements obtained in this way are in agreement with near-field strong motion data and GPS observations recorded in Lorca. We have obtained a partly bilateral rupture propagating to WSW (238°, 54°, 59°) with 27 cm of maximum slip and shallow focus (4 km). The fault plane orientation corresponds to that of the Cejo de los Enamorados Fault located NE of the Lorca town and parallel to the Alhama de Murcia Fault. The distribution of slip on the fault plane can explain the lack of any observed surface rupture as we found that the rupture started at 4-km depth along a plane dipping at 54°, with motion propagating upward to stop at 1.5 km below the surface. The strong motion and GPS data recorded near the epicenter are in agreement with the maximum slip on the fault. Directivity effects and the extreme shallowness of the rupture could explain the considerable damage that the earthquake caused in the town of Lorca.     

Quasi-Static Strain and Stress Fields due to a Moment Tensor in Elastic–Viscoelastic Layered Half-Space

We derived explicit expressions in the time domain for 3-D quasi-static strain and stress fields, due to a point moment tensor source in an elastic surface layer overlying viscoelastic half-space under gravity. The expressions of strain in the elastic surface layer were directly obtained from the expressions of displacement in our previous paper. The conversion of strain into stress is easy, because the stress–strain relation of elastic material is linear. In the viscoelastic substratum, the expressions of strain were obtained by applying the correspondence principle of linear viscoelasticity to the associated elastic solution. The strain–stress conversion is not straightforward, as the stress–strain relation of viscoelastic material is usually given in a differential form. To convert strain into stress, we used an integral form of the stress–strain relation instead of the usual differential form. The expressions give the responses of elastic half-space at \( t = 0 \) , and the responses of an elastic plate floating on non-viscous liquid at \( t = \infty \) . The moment tensor is rationally decomposed into the three independent force systems, corresponding to isotropic expansion, shear faulting and crack opening, and so the expressions include the strain and stress fields for these force systems as special cases. As the first numerical example, we computed the temporal changes in strain and stress fields after the sudden opening of an infinitely long vertical crack cutting the elastic surface layer. Here, we observe that the stress changes caused by the sudden crack opening gradually decay with time and vanish at \( t = \infty \) everywhere. After the completion of stress relaxation, a characteristic pattern of shear strain remains in the viscoelastic substratum. Since the strain and stress fields at \( t = \infty \) can be read as the strain- and stress-rate fields caused by steady crack opening, respectively, this numerical example demonstrates the realization of a steady stress state supported by steady viscous flow in the asthenosphere, associated with steady seafloor spreading at mid-ocean ridges. For the second numerical example, we computed the temporal changes in strain and stress fields after the 2011 Tohoku-oki mega-thrust earthquake, which occurred at the North American-Pacific plate interface. In this numerical example, the stress changes caused by coseismic fault slip vanish at \( t = \infty \) in the viscoelastic substratum, but remain in the elastic surface layer. The coseismic stress changes (and also strain changes) in the elastic surface layer diffuse away from the source region with time, due to gradual stress relaxation in the viscoelastic substratum.     

Estimating ground motion incoherence through finite source simulation: a case study of the 1980 El-Asnam Earthquake

Spatial variability of ground motions has significant influence on dynamic response of extended structures such as bridges and tunnels. In this study, the widely used finite-source ground motion simulation approach, the so-called Empirical Green’s Function (EGF) method, is extended to synthesize seismic motions across an array of stations located at bedrock in the epicentral region of the 1980 El-Asnam region (North-West Algeria). The target event being simulated is the October 10 1980 \( M_{s} = 7.2 \) Earthquake, and the EGF is obtained from the ground motion recorded at Sogedia Factory station during the 8 November 1980 \( M_{L} = 5.6 \) aftershock. Coherency functions are then estimated from the simulated ground accelerations. A parametric study investigating the influence of shear wave velocity, earthquake magnitude, and epicentral distance is conducted by simulating ground acceleration for different scenarios using the Hybrid Green’s Function method. The main finding of the study is that finite source effects can cause significant loss in coherency at bedrock in the near-field. In the far-field, the source effect alone does not seem to produce incoherent motion, which implies that scattering and local site effects could be dominating there. Furthermore, coherency functions are found to be more sensitive to inter-station separation in the near-field than in the far-field. Increasing shear wave velocity seems to increase coherency functions, and larger earthquakes seem to produce more incoherent motion than smaller ones. The simulation method presented here produces incoherent motion mainly due to the finite source effect, while path effects are partially accounted for through the EGF, and local site effects are not considered. In this sense, the estimated coherency functions represent that of plane waves. A parametric model of plane wave coherency is calibrated and presented based on the simulation results. The results indicate that the parametric model can be used as a first approximation, and at least an upper bound of lagged coherency in the near-field region of the El-Asnam Earthquake scenario. This model could be useful in random vibration analysis or generation of spatially variable ground motion for time history analysis of lifeline structures in the study area.     

Consideration of aging and SSI effects on seismic vulnerability assessment of RC buildings

At present, the seismic vulnerability assessment of reinforced concrete (RC) buildings is made considering fixed base conditions; moreover, the mechanical properties of the building remain intact in time. In this study we investigate whether these two fundamental hypotheses are sound as aging and soil-structure interaction (SSI) effects might play a crucial role in the seismic fragility analysis of RC structures. Among the various aging processes, we consider the chloride-induced corrosion based on probabilistic modeling of corrosion initiation time and corrosion rate. Different corrosion aspects are considered in the analysis including the loss of reinforcement cross-sectional area, the degradation of concrete cover and the reduction of steel ultimate deformation. SSI is modeled by applying the direct one-step approach, which accounts simultaneously for inertial and kinematic interactions. Two-dimensional incremental dynamic analysis is performed to assess the seismic performance of the initial uncorroded ( \(\hbox {t}=0\) years) and corroded ( \(\hbox {t}=50\) years) RC moment resisting frame structures, having been designed with different seismic code levels. The time-dependent fragility functions are derived in terms of the spectral acceleration at the fundamental mode of the structure \(\hbox {S}_{\mathrm{a}}(\hbox {T}_{1}\) , 5 %) and the outcropping peak ground acceleration for the immediate occupancy and collapse prevention limit states. Results show an overall increase in seismic vulnerability over time due to corrosion highlighting the important influence of deterioration due to aging effects on the structural behavior. Moreover, the consideration of SSI and site effects may significantly alter the expected structural performance leading to higher vulnerability values.     

Simulation of Strong Ground Motion During the 1950 Great Assam Earthquake

In this paper, ground motion during the Independence Day earthquake of August 15, 1950 (M_w 8.6, Ben-Menahem et al., 1974) in the northeastern part of India is estimated by seismological approaches. A hybrid simulation technique which combines the low frequency ground motion simulated from an analytical source mechanism model with the stochastically simulated high-frequency components is used for obtaining the acceleration time histories. A series of ground motion simulations are carried out to estimate the peak ground acceleration (PGA) and spectral accelerations at important cities and towns in the epicentral region. One sample PGA distribution in the epicentral region encompassing the epicenter is also obtained. It is found that PGA in the epicentral region has exceeded 1 g during this earthquake. The estimated PGA’s are validated to the extent possible using the MMI values. The simulated acceleration time histories can be used for the assessment of important engineering structures in northeastern India.     

Full Moment Tensor Variations and Isotropic Characteristics of Earthquakes in the Gulf of California Transform Fault System

The full moment tensor is a mathematical expression of six independent variables; however, on a routine basis, it is a common practice to reduce them to five assuming that the isotropic component is zero. This constraint is valid in most tectonic regimes where slip occurs entirely at the fault surface (e.g. subduction zones); however, we found that full moment tensors are best represented in transform fault systems. Here we present a method to analyze source complexity of earthquakes of different sizes using a simple formulation that relates the elastic constants obtained from independent studies with the angle between the slip and the fault normal vector, referred to as angle \( \theta \) ; this angle is obtained from the full moment tensors. The angle \( \theta \) , the proportion of volume change \( \left( k \right) \) and the constant volume (shear) component \( \left( T \right) \) are numerical indicators of complexity of the source; earthquakes are more complex as \( \theta \) deviates from \( \pi /2 \) or as T and k deviate from zero as well. These parameters are obtained from the eigensolution of the full moment tensor. We analyzed earthquakes in the Gulf of California that exhibit a clear isotropic component and we observed that the constant volume parameter T is independent of scalar moments, suggesting that big and small earthquakes are equally complex. In addition, simple models of one single fault are not sufficient to describe physically all the combinations of \( \theta \) in a source type plot. We also found that the principal direction of the strike of the Transform Fault System in the Gulf of California is following the first order approximation of the normal surface of the full moment tensor solution, whereas for deviatoric moment tensors the principal direction does not coincide with the strike of the Transform Fault System. Our observations that small and large earthquakes are equally complex are in agreement with recent studies of strike-slip earthquakes.     

A numerical study on the 2D behavior of the single and layered clayey basins

In this numerical study the effects of the stratification and change in material type on the dynamic behavior of the clayey basins are investigated. For this purpose two kinds of single and double layer basin models with the total depth of 100 m and width of 2 km have been selected. The single layer basins consist of soft, medium and stiff clayey materials, while double layer basins consist of a stiff sandy sub-layer overlaid by one of the soft, medium and stiff clayey soil layer. A range of bedrock inclinations at the valley sides from slighter \(10^{\circ }\) and \(20^{\circ }\) to steeper \(30^{\circ }\) and \(40^{\circ }\) are selected. A numerical study using finite difference based nonlinear code which utilizes appropriate static and dynamic boundary conditions, and includes hysteresis damping formulation based on user defined degradation curves is conducted using real earthquake motions of different strength and frequency content. The results of the 2D analyses show the differences in dynamic behavior of single and double layer basins from the aspect of amount and position of the maximum amplification and resonance period. It was seen that the shallower lateral parts of the basins are sensitive to lower periods while inner parts are sensitive to higher periods. Among other differences, the most important difference between the behavior of single and double layer basins was seen at the resonance period of the inner parts of the basins.     

New Seismic Noise Models Obtained Using Very Broadband Stations

It has been two decades since the last comprehensive standard model of ambient earth noise was published Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The PETERSON model was updated by analyzing the absolute quietest conditions for stations within the GSN (Berger et al. in J Geophys Res 109, 2005; Mcnamara and Buland in Bull Seism Soc Am 94:1517–1527, 2004; Ringler et al. in Seismol Res Lett 81(4) doi:10.1785/gssrl.81.4.605, 2010). Unfortunately, both the original model and the updated models did not include any deployed station in North Africa and Middle East, which reflects the noise levels within the desert environment of those regions. In this study, a survey was conducted to create a new seismic noise model from very broadband stations which recently deployed in North Africa. For this purpose, 1 year of continuous recording of seismic noise data of the Egyptian National Seismic Network (ENSN) was analyzed in order to create a new noise model. Seasonal and diurnal variations in noise spectra were recorded at each station. Moreover, we constructed a new noise model for each individual station. Finally, we obtained a new cumulative noise model for all the stations. We compared the new high-noise model (EHNM) and new low-noise model (ELNM) with both the high-noise model (NHNM) and low-noise model (NLNM) of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The obtained noise levels are considerably lower than low-noise model of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993) at ultra long period band (ULP band), but they are still below the high-noise model of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The results of this study could be considered as a first step to create permanent seismic noise models for North Africa and Middle East regions.     

Correlations of three-component ground motions from the 23 June 2010 Val-des-Bois, Québec, earthquake

In this paper, we investigate the characteristics of three-component ground motions recorded during the $M_\text{w}\,5.2$ Val-des-Bois (Québec) earthquake, which occurred on the 23 June 2010. The earthquake is the largest recorded event in eastern Canada within the last decade. The records analyzed were provided by a strong motion monitoring network, comprising accelerometers located at sites with different soil conditions. The two orthogonal horizontal components and one vertical component at each recording station are uncorrelated to determine their principal directions, and the results obtained are used to characterize intensity ratios between the three uncorrelated components. A new hodograph representation is proposed to highlight the correlation between accelerations and displacement trajectories along various time increments at each recording station. The principal components are discussed in light of seismographic data, local site conditions, and trajectories. Time–frequency analyses of the uncorrelated records are also conducted to compare the distribution of spectral amplitudes and frequency content along the three principal components during the shaking. The results of this work shed more light on the characteristics of three-component ground motions from an important Eastern North America earthquake, and could be used to calibrate simulated multicomponent ground motions in this region.     

Flood basalts from dharwars of Mysore,India

The term flood basalt is redefined emphasizing the importance of the subaerial environment. Using the well established physical criteria of aerial extent, internal structures, time of extrusion and associations, flood basalt activity is distinguished in the Archeans (Dharwars) of Mysore from the geosynclinal volcanics. Study of chemical composition of the Dharwar and other Archean volcanics in the light ofSugimura’s (1968) SWS index, and plotting of the chemical analyses on theMacdonald andKatsura’s (1964) alkali-silica diagram,Kuno’s (1968) alkali-alumina-silica diagram and Scheynamann’s silica-Niggli qz diagram shows both geosynclinal and subaerial volcanics are mainly tholeiitic. Therefore in deciphering the environment of volcanism, it is suggested that the physical criteria take precedence over chemical composition.     

Ground motion simulations for the 23 October 2011 Van,Eastern Turkey earthquake using stochastic finite fault approach

The 23 October 2011 Van (Mw 7.1) earthquake that occurred in Eastern Turkey resulted in heavy damage particularly in the city of Van and town of Ercis. This paper presents ground motion simulations of Van earthquake by using stochastic finite fault method (EXSIM, Motazedian and Atkinson in Bull Seismol Soc Am 95:995–1010, 2005; Boore in Bull Seismol Soc Am 99:3202–3216, 2009) that provides a simple and effective tool to generate high frequency strong motion. The input parameters related to source, path, and site effects are calibrated on the basis of minimizing the error functions between simulations and observations both in time and frequency domain. Validated model parameters are used to produce synthetics in regional extent with the aim of understanding the level and distribution of the ground shaking particularly in the near fault region where no recordings are available within the 40 km of the epicenter. This paper evaluates the effect of two different slip models on ground motion intensity measures over the area of interest and addresses the variability in the near fault region associated with the source effect. The synthetics are compared with the corresponding estimations of ground motion prediction equations by Boore and Atkinson (Earthq Spectra 24:99–138, 2008), Akkar and Bommer (Seismol Res Lett 81:195–206, 2010) and Akkar and Cagnan (Bull Seismol Soc Am 100:2978–2995, 2010). Our results indicate that despite the limitation of the method for incorporating the directivity effect and inadequate representation of the soil conditions at the individual stations, a satisfactory match between synthetics and observations are obtained both in time and frequency domain. Spatial distributions of the synthetics in regional level also show reasonable correlation with ground motion prediction equations and damage observations.     

云南禄劝地震部分余震的矩张量反演

利用数字盒式磁带记录加速度仪组成的流动地震台网所记录的三分向近场加速度图,通过矩张量反演确定了1985年4月18日云南禄劝 Ms=6.1地震的部分余震的震源机制.以均匀半无限弹性介质的格林函数解释路径效应,并通过正演计算识别由两次积分得到的位移地震图中的直达 P 波、直达 S 波和 SP 转换波震相,然后用这些震相进行矩张量反演.反演结果表明,在解超定线性方程组时,采用适当的加权系数,可使上述直达波和转换波的理论计算值与观测值拟合得较好.尽管用以反演的三个余震大小不同(震级 ML48,3.2,3.5),震源位置也有差别,但它们的震源机制却非常接近,且与主震的震源机制相当一致.这一特征显示了余震的发生与主震发震构造的内在联系.这些实例说明,由震源球球面上分布适宜的数字地震台网取得的近场加速度资料,借助于即使是简单的介质模型,通过地震矩张量反演,不但可以得出这些地震震源的主要成份————剪切位错源,同时还可得出震源所含的其它成分,如膨胀源和补偿线性向量偶极.