首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 953 毫秒
1.
In this paper we describe a stable automatic method to estimate in real time the seismic moment, moment magnitude and corner frequency of events recorded by a network comprising broad-band and accelerometer sensors. The procedure produces reliable results even for small-magnitude events $\hbox {M}_{\mathrm{W}}\approx 3$ . The real-time data arise from both the Transfrontier network at the Alps-Dinarides junction and from the Italian National Accelerometric Network (RAN). The data is pre-processed and the S-wave train identified through the application of an automatic method, which estimates the arrival times based on the hypocenter location, recording site and regional velocity model. The transverse component of motion is used to minimize conversion effects. The source spectrum is obtained by correcting the signals for geometrical spreading and intrinsic attenuation. Source spectra for both velocity and displacement are computed and, following Andrews (1986), the seismic moment and the first estimate of the corner frequency, $f_{0}$ , derived. The procedure is validated using the recordings of some recent moderate earthquakes (Carnia 2002; Bovec 2004; Parma 2008; Aquila 2009; Macerata 2009; Emilia 2012) and the recordings of some minor events in the SE Alps area for which independent seismic moment and moment magnitude estimates are available. The results obtained with a dataset of 843 events recorded by the Transfrontier and RAN networks show that the procedure is reliable and robust for events with $\hbox {M}_{\mathrm{W}}\ge 3$ . The estimates of $f_{0}$ are less reliable. The results show a scatter, principally for small events with $\hbox {M}_{\mathrm{W}}\le 3$ , probably due to site effects and inaccurate locations.  相似文献   

2.
3.
4.
Local site effect microzonation of Lorca town (SE Spain)   总被引:1,自引:0,他引:1  
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.  相似文献   

5.
Ground-motion prediction equations (GMPEs) are essential tools in seismic hazard studies to estimate ground motions generated by potential seismic sources. Global GMPEs which are based on well-compiled global strong-motion databanks, have certain advantages over local GMPEs, including more sophisticated parameters in terms of distance, faulting style, and site classification but cannot guarantee the local/region-specific propagation characteristics of shear wave (e.g., geometric spreading behavior, quality factor) for different seismic regions at larger distances (beyond about 80 km). Here, strong-motion records of northern Iran have been used to estimate the propagation characteristics of shear wave and determine the region-specific adjustment parameters for three of the NGA-West2 GMPEs to be applicable in northern Iran. The dataset consists of 260 three-component records from 28 earthquakes, recorded at 139 stations, with moment magnitudes between 4.9 and 7.4, horizontal distance to the surface projection of the rupture (R JB) less than 200 km, and average shear-wave velocity over the top 30 m of the subsurface (V S30) between 155 and 1500 m/s. The paper also presents the ranking results for three of the NGA-West2 GMPEs against strong motions recorded in northern Iran, before and after adjustment for region-dependent attenuation characteristics. The ranking is based on the likelihood and log-likelihood methods (LH and LLH) proposed by Scherbaum et al. (Bull Seismol Soc Am 94: 2164–2185, 2004, Bull Seismol Soc Am 99, 3234–3247, 2009, respectively), the Nash–Sutcliffe model efficiency coefficient (Nash and Sutcliffe, J Hydrol 10:282–290, 1970), and the EDR method of Kale and Akkar (Bull Seismol Soc Am 103:1069–1084, 2012). The best-fitting models over the whole frequency range are the ASK14 and BSSA14 models. Taking into account that the models’ performances were boosted after applying the adjustment factors, at least moderate regional variation of ground motions is highlighted. The regional adjustment based on the Iranian database reveals an upward trend (indicated as high Q factor) for the selected database. Further investigation to determine adjustment factors based on a much richer database of the Iranian strong-motion records is of utmost important for seismic hazard and risk analysis studies in northern Iran, containing major cities including the capital city of Tehran.  相似文献   

6.
We estimate the corner frequencies of 20 crustal seismic events from mainshock–aftershock sequences in different tectonic environments (mainshocks 5.7 < M W < 7.6) using the well-established seismic coda ratio technique (Mayeda et al. in Geophys Res Lett 34:L11303, 2007; Mayeda and Malagnini in Geophys Res Lett, 2010), which provides optimal stability and does not require path or site corrections. For each sequence, we assumed the Brune source model and estimated all the events’ corner frequencies and associated apparent stresses following the MDAC spectral formulation of Walter and Taylor (A revised magnitude and distance amplitude correction (MDAC2) procedure for regional seismic discriminants, 2001), which allows for the possibility of non-self-similar source scaling. Within each sequence, we observe a systematic deviation from the self-similar \( M_{0} \propto \mathop f\nolimits_{\text{c}}^{ - 3} \) line, all data being rather compatible with \( M_{0} \propto \mathop f\nolimits_{\text{c}}^{ - (3 + \varepsilon )} \) , where ε > 0 (Kanamori and Rivera in Bull Seismol Soc Am 94:314–319, 2004). The deviation from a strict self-similar behavior within each earthquake sequence of our collection is indicated by a systematic increase in the estimated average static stress drop and apparent stress with increasing seismic moment (moment magnitude). Our favored physical interpretation for the increased apparent stress with earthquake size is a progressive frictional weakening for increasing seismic slip, in agreement with recent results obtained in laboratory experiments performed on state-of-the-art apparatuses at slip rates of the order of 1 m/s or larger. At smaller magnitudes (M W < 5.5), the overall data set is characterized by a variability in apparent stress of almost three orders of magnitude, mostly from the scatter observed in strike-slip sequences. Larger events (M W > 5.5) show much less variability: about one order of magnitude. It appears that the apparent stress (and static stress drop) does not grow indefinitely at larger magnitudes: for example, in the case of the Chi–Chi sequence (the best sampled sequence between M W 5 and 6.5), some roughly constant stress parameters characterize earthquakes larger than M W ~ 5.5. A representative fault slip for M W 5.5 is a few tens of centimeters (e.g., Ide and Takeo in J Geophys Res 102:27379–27391, 1997), which corresponds to the slip amount at which effective lubrication is observed, according to recent laboratory friction experiments performed at seismic slip velocities (V ~ 1 m/s) and normal stresses representative of crustal depths (Di Toro et al. in Nature in press, 2011, and references therein). If the observed deviation from self-similar scaling is explained in terms of an asymptotic increase in apparent stress (Malagnini et al. in Pure Appl Geophys, 2014, this volume), which is directly related to dynamic stress drop on the fault, one interpretation is that for a seismic slip of a few tens of centimeters (M W ~ 5.5) or larger, a fully lubricated frictional state may be asymptotically approached.  相似文献   

7.
We study the geometrical and material conditions which lead to focusing of seismic waves traveling across a concave velocity interface representing the boundary of a sedimentary basin within a denser rock. We approximate, using geometrical analysis for plane-waves, the combination of interface eccentricities and velocity ratios for which the seismic rays converge to a near surface region of the basin. 2-D finite difference modeling is used to compute Peak Ground Velocity (PGV) and spectral amplification across the basin. We show that effective geometrical focusing occurs for a narrow set of eccentricities and velocity ratios, where seismic energy is converged to a region of $\pm $ 0.5 km from surface. This mechanism leads to significant amplification of PGV at the center of the basin, up to a factor of 3; frequencies of the modeled spectrum are amplified up to the corner frequency of the source. Finally, we suggest a practical method for evaluating the potential for effective geometrical focusing in sedimentary basins.  相似文献   

8.
To investigate the nature of Pn propagation, we have implemented the spectral-element method (SEM) for vertically and laterally varying media with and without attenuation. As a practical measure, essential features of the Pn waves are distilled into seismic attributes including arrival times, amplitudes and pulse frequencies. To validate the SEM simulations, we first compare the SEM results with reflectivity calculations of Braile and Smith (Geophys. J.R. Astr. Soc. 40, 145–176, 1975) and then to the asymptotic results of ?erveny and Ravindra (Theory of Seismic Headwaves, University of Toronto Press, pp. 235–250, 1971). Models with random, laterally varying Moho structures are then simulated, where the amplitude and pulse frequency characteristics are found to be stable to small Moho interface perturbations. SEM calculations for models with different upper-mantle velocity gradients are next performed where it is found that interference effects can strongly influence the Pn amplitudes and pulses frequencies. For larger-scale, laterally varying structures, SEM models similar to that found along the Hi-CLIMB array in Tibet are then performed. It is observed that large-scale structures, along with small-scale structures, upper-mantle velocity gradients and attenuation, can all significantly affect the Pn attributes. Ambiguities between upper-mantle velocity gradients and attenuation are also found when using Pn amplitudes and pulse frequency attributes. These ambiguities may be resolved, to some degree, by using the curvature of the travel times at longer regional distance, however, this would also be complicated by lateral variability.  相似文献   

9.
We estimate seismological fracture energies from two subsets of events selected from the seismic sequences of L’Aquila (2009), and Northridge (1994): 57 and 16 selected events, respectively, including the main shocks. Following Abercrombie and Rice (Geophys J Int 162: 406–424, 2005), we postulate that fracture energy (G) represents the post-failure integral of the dynamic weakening curve, which is described by the evolution of shear traction as a function of slip. Following a direct-wave approach, we compute mainshock-/aftershock-source spectral ratios, and analyze them using the approach proposed by Malagnini et al. (Pure Appl. Geophys., this issue, 2014) to infer corner frequencies and seismic moment. Our estimates of source parameters (including fracture energies) are based on best-fit grid-searches performed over empirical source spectral ratios. We quantify the source scaling of spectra from small and large earthquakes by using the MDAC formulation of Walter and Taylor (A revised Magnitude and Distance Amplitude Correction (MDAC2) procedure for regional seismic discriminants, 2001). The source parameters presented in this paper must be considered as point-source estimates representing averages calculated over specific ruptured portions of the fault area. In order to constrain the scaling of fracture energy with coseismic slip, we investigate two different slip-weakening functions to model the shear traction as a function of slip: (i) a power law, as suggested by Abercrombie and Rice (Geophys J Int 162: 406–424, 2005), and (ii) an exponential decay. Our results show that the exponential decay of stress on the fault allows a good fit between measured and predicted fracture energies, both for the main events and for their aftershocks, regardless of the significant differences in the energy budgets between the large (main) and small earthquakes (aftershocks). Using the power-law slip-weakening function would lead us to a very different situation: in our two investigated sequences, if the aftershock scaling is extrapolated to events with large slips, a power law (a la Abercrombie and Rice) would predict unrealistically large stress drops for large, main earthquakes. We conclude that the exponential stress evolution law has the advantage of avoiding unrealistic stress drops and unbounded fracture energies at large slip values, while still describing the abrupt shear-stress degradation observed in high-velocity laboratory experiments (e.g., Di Toro et al., Fault lubrication during earthquakes, Nature 2011).  相似文献   

10.
We determined depth variation of the 410- and 660-km discontinuities beneath southeastern China by common-converted-point stacking of \(\rm P\) -wave receiver functions of 121 permanent Chinese seismic stations. We then combined the results with seismic velocity variation to estimate temperature and water content variations in the mantle transition zone of the region. Previous tomographic studies have shown a stagnant slab in the mantle transition zone in eastern Asia that is connected to subduction of the western Pacific. Temperature variations obtained clearly outline the shape of the stagnant slab, with its western edge at 113.5 \(^\circ\) E and the southern edge at 28.5 \(^\circ\) N. The correlation between the location of the stagnant slab and surface tectonics suggests that the Cenozoic extension in eastern China is closely associated with the subduction of the western Pacific and its eastward migration. The water content of the stagnant slab is lower than in surrounding slabs, suggesting that the water has already been released from the subducting slab into the upper mantle.  相似文献   

11.
The main focus of this paper is on uncertainties associated with the selection of the incoherency coefficient and apparent propagation velocity in the simulation of asynchronous seismic excitations based on the provisions of CEN [Comité Européen de Normalisation] (Eurocode 8: Design of structures for earthquake resistance—part 2: bridges, BS EN 1998-2:2005, Brussels, Belgium 2005). To this effect, the importance of utilizing appropriate values of these two parameters in the simulation of spatially varying seismic excitations is highlighted, and practice-oriented recommendations are provided for their selection. In addition, preliminary probability distributions are proposed for the consideration of the effect of uncertainty in the selection of the apparent propagation velocity for reliability (fragility) analyses. The effect of multi-support input motions on the seismic response of bridges is also discussed.  相似文献   

12.
地震台阵常用频率-波数分析法来估计地震信号的后方位角和慢度。尽管有关的算法是众所周知的,但不同的实现在某些情况下可能导致不同的结果。如应用在宽频带的澳大利亚WRA台阵的记录时,标准的f-k分析方法往往给出不正确的结果。本文发现错误是由FFT的频谱泄漏效应引起的,如在进行FFT之前对原始数据进行高通滤波则可以有效地避免这样的错误。进一步对在时域中直接估算信号方位角和慢度的一种算法进行了分析,比较发现时域方法的计算速度和精度都不低于频域方法,且在某些低信噪比的情况下,前者可以给出更可靠的结果。  相似文献   

13.
—The specific barrier model is used for the first time to simulate ground motion accelerations for the purpose of probabilistic seismic hazard analyses at sites near a dominant fault system. It incorporates the simulation of fault geometry and the relationship between the stress drop and seismic moment to estimate the number of cracks on the fault for the specific barrier model. Radiated direct shear waves are established following Boore’s (1983) procedure. The simulated peak ground accelerations (PGA) are then calibrated by strong-motion data. Basically, the model is of uniform source, and the directivity of the source is not taken into consideration. The results show that the calibrated PGA values are not sensitive to the relationship between the stress drop and seismic moment. However, the calibrated PGA values may increase about 20 percent for sites near the fault when the cut-off frequency,?f max?, is raised from 5 Hz to 10 Hz. The variability of the simulated ground motion is, in general, smaller than that of the empirical strong-motion data shown in the literature. This may be improved by adding randomness into the parameter of ?f max and uncertainties into the empirical relationships adopted in the model. The simulated attenuation curves may be used to judge which types of conventional attenuation equations are better at representing the attenuation of PGA for sites near the fault, especially for large earthquake events.  相似文献   

14.
Using seismic data from regional earthquakes in Tibet recorded by the Hi-CLIMB experiment, Pn attributes are used to constrain the velocity gradient and attenuation structure of the Tibetan lithosphere under the Hi-CLIMB array. Numerical modeling is performed using the spectral-element method (SEM) for laterally varying upper-mantle velocity and attenuation, and the seismic attributes considered include the Pn travel-time, envelope amplitude, and pulse frequency. The results from the SEM modeling provide two alternative models for the upper-mantle beneath the Hi-CLIMB array in Tibet. The first model is derived from the 3D velocity model of Griffin et al. (Bull Seism Soc Am 101:1938–1947, 2011) with a constant upper-mantle velocity gradient, and laterally varying upper mantle attenuation. The second model has a laterally varying upper-mantle velocity gradient, and constant upper-mantle attenuation. In both cases, the Qiangtang terrane is distinguished from the Lhasa terrane by a change in Moho depth and upper-mantle velocities. The lower upper-mantle velocities, as well as higher Pn attenuation, suggest hotter temperatures beneath the Qiangtang terrane as compared to the Lhasa terrane. Although the fits to the Pn amplitude and pulse frequency data are comparable between the two models, the first model with the constant upper-mantle velocity gradient fits the travel times somewhat better in relation to the data errors.  相似文献   

15.
The seismic behaviour of caisson foundations supporting typical bridge piers is analysed with 3D finite elements, with due consideration to soil and interface nonlinearities. Single-degree-of freedom oscillators of varying mass and height, simulating heavily and lightly loaded bridge piers, founded on similar caissons are studied. Four different combinations of the static ( $\text{ FS }_\mathrm{V}$ FS V ) and seismic ( $\text{ FS }_\mathrm{E}$ FS E ) factors of safety are examined: (1) a lightly loaded ( $\text{ FS }_\mathrm{V}= 5$ FS V = 5 ) seismically under-designed ( $\text{ FS }_\mathrm{E} < 1$ FS E < 1 ) caisson, (2) a lightly loaded seismically over-designed ( $\text{ FS }_\mathrm{E} >1$ FS E > 1 ) caisson, (3) a heavily loaded ( $\text{ FS }_\mathrm{V} = 2.5$ FS V = 2.5 ) seismically under-designed ( $\text{ FS }_\mathrm{E} < 1$ FS E < 1 ) caisson and (4) a heavily loaded seismically over-designed caisson. The analysis is performed with use of seismic records appropriately modified so that the effective response periods (due to soil-structure-interaction effects) of the studied systems correspond to the same spectral acceleration, thus allowing their inelastic seismic performance to be compared on a fair basis. Key performance measures of the systems are then contrasted, such as: accelerations, displacements, rotations and settlements. It is shown that the performance of the lightly loaded seismically under-designed caisson is advantageous: not only does it reduce significantly the seismic load to the superstructure, but it also produces minimal residual displacements of the foundation. For heavily loaded foundations, however ( $\text{ FS }_{V} = 2.5$ FS V = 2.5 ), the performance of the two systems (over and under designed) is similar.  相似文献   

16.
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.  相似文献   

17.
— Surface-wave amplitudes from explosion sources show less variation for a given event han body wave amplitudes, so it is natural to expect that yield estimates derived from surface waves will be more accurate than yield estimates derived from body waves. However, yield estimation from surface waves is complicated by the presence of tectonic strain release, which acts like one or more earthquake sources superimposed on top of the explosion. Moment-tensor inversion can be used to remove the tectonic component of the surface waves, however moment-tensor inversion for shallow sources is inherently non-unique so the explosion isotropic moment cannot be determined with the necessary accuracy by this means. Explosions on an island or near a mountain slope can exhibit anomalous surface waves similar to those caused by tectonic strain release. These complications cause yield estimates derived from surface waves to be less accurate than yield estimates from body waves recorded on a well-calibrated network with good coverage. Surface-wave amplitudes can be expressed as a surface-wave magnitude M s , which is defined as the logarithm of the amplitude plus a distance correction, or as a path corrected spectral magnitude, log $M^{\prime}_0$ , which is derived from the surface-wave spectrum. We derive relations for M s vs. yield and log $M^{\prime}_0$ vs. yield for a large data set and estimate the accuracy of these estimates.  相似文献   

18.
In the last two decades, south-central Europe and the Eastern Alps have been widely explored by many seismic refraction experiments (e.g., CELEBRATION 2000, ALP 2002, SUDETES 2003). Although quite detailed images are available along linear profiles, a comprehensive, three-dimensional crustal model of the region is still missing. This limitation makes this region a weak spot in continental-wide comprehensive representations of crustal structure. To improve on this situation, we select and collect 37 published active-source seismic lines in this region. After geo-referencing each line, we sample them along vertical profiles—every 50?km or less along the line—and derive P-wave velocities in a stack of homogeneous layers (separated by discontinuities: depth of crystalline basement, top of lower crust, and Moho). We finally merge the information using geostatistical methods, and infer S-wave velocity and density using empirical scaling relations. We present here the resulting crustal model for a region encompassing the Eastern Alps, Dinarides, Pannonian basin, Western Carpathians and Bohemian Massif, covering the region within $45^{\circ}\text{--}51^{\circ}\hbox{N}$ and $11^{\circ} \text{--} 22^{\circ}\hbox{E}$ with a resolution of $0.2^{\circ} \times 0.2^{\circ}.$ We are also able to extend and update the map of Moho depth in a wider region within $35^{\circ}\text{--}51^{\circ}\hbox{N}$ and $12^{\circ}\text{--}45^{\circ}\hbox{E},$ gathering Moho values from the collected seismic lines, other published dataset and using the European plate reference EPcrust as a background. All the digitized profiles and the resulting model are available online.  相似文献   

19.
One of the major challenges related with the current practice in seismic hazard studies is the adjustment of empirical ground motion prediction equations (GMPEs) to different seismological environments. We believe that the key to accommodating differences in regional seismological attributes of a ground motion model lies in the Fourier spectrum. In the present study, we attempt to explore a new approach for the development of response spectral GMPEs, which is fully consistent with linear system theory when it comes to adjustment issues. This approach consists of developing empirical prediction equations for Fourier spectra and for a particular duration estimate of ground motion which is tuned to optimize the fit between response spectra obtained through the random vibration theory framework and the classical way. The presented analysis for the development of GMPEs is performed on the recently compiled reference database for seismic ground motion in Europe (RESORCE-2012). Although, the main motivation for the presented approach is the adjustability and the use of the corresponding model to generate data driven host-to-target conversions, even as a standalone response spectral model it compares reasonably well with the GMPEs of Ambraseys et al. (Bull Earthq Eng 3:1–53, 2005), Akkar and Bommer (Seismol Res Lett 81(2):195–206, 2010) and Akkar and Cagnan (Bull Seismol Soc Am 100(6):2978–2995, 2010).  相似文献   

20.
A damaging seismic sequence hit a wide area mainly located in the Emilia-Romagna region (Northern Italy) during 2012 with several events of local magnitude \(\hbox {M}_\mathrm{l} \ge 5\) , among which the \(\hbox {M}_\mathrm{l}\) 5.9 May 20 and the \(\hbox {M}_\mathrm{l}\) 5.8 May 29 were the main events. Thanks to the presence of a permanent accelerometric station very close to the epicentre and to the temporary installations performed in the aftermath of the first shock, a large number of strong motion recordings are available, on the basis of which, we compared the recorded signals with the values provided by the current Italian seismic regulations, and we observed several differences with respect to horizontal components when the simplified approach for site conditions (based on Vs30 classes) is used. On the contrary, when using the more accurate approach based on the local seismic response, we generally obtain a much better agreement, at least in the frequency range corresponding to a quarter wavelength comparable with the depth of the available subsoil data. Some unresolved questions still remain, such as the low frequency behaviour ( \(<\) 1 Hz) that could be due either to complex propagation at depth larger than the one presently investigated or to near source effects, and the behaviour of vertical spectra whose recorded/code difference is too large to be explained with the information currently available.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号